某跨江大桥地基基础稳定性分析中的若干问题探讨 旱 鑫兽 仲 ｒ．飞 （江苏省建苑岩土工程勘测有限公司，江苏南京２１００２９） 摘 要：以某跨江大桥为例，针对地基基础设计中存在的诸如分析方 目前大多数研究群桩一土一承台结构共同作用 分析方法的出发点都基于弹性半空间基础上建立的 布辛奈斯克（Ｂｏｕｓｓｉｎｅｓｑ）解和明德林（Ｍｉｎｄｌｉｎ）解。两者 法、本构模型和接触面单元、计算参数以及施工仿真等方面的有争议 的部分问题进行了分析，提出了解决途径并明确了在指导设计中的意 一义，为指导设计、施工提供理论指导。 关键词：地基基础；稳定性分析；本构模型；施工仿真 中图分类号：Ｕ４４５．５５ 文献标志码：Ｂ 文章编号：１００７—７３５９（２０１６）０６—０１２４—０３ ＤＯｌ：１０．１６３３ｏ４．ｃｎｋｉ．１００７－７３５９．２０１６．０６．０４７ １工程概况 跨江大桥建设规模大、投资高、安全ｌ生问题突出， 大桥基础的稳定性问题一直是交通道桥领域研究的 热点和难点问题。某跨江大桥地基第四纪覆盖层厚度 达３００ｍ，桥墩下采用大直径钻孑Ｌ灌注桩后注浆。由于 覆盖层厚度大，桩端实际位于相对较软弱的土基上。 对于此类超大直径的群桩基础是否稳定、桩身轴力分 布、桩承台内应力、桩基础沉降、水平位移与各桥墩的 位移等问题都应该明确。 大型垮江桥梁的基础对沉降和变形以及受力情 况有着极其严格的要求。某跨江大桥索塔基础采用直 径２．５ｍ，桩长１２０ｍ，平均桩间距６．２５ｍ的超大、超长 的钻孔灌注桩基础，对于大直径、超长度以及复杂水 文地质条件下的桩基础设计是否能够按照现有传统 的基础设计规范进行设计，是设计人员重点关注的问 题。另外，基础、结构以及地基土层作为共同作用体， 传统的规范方法无法模拟三者的共同作用，基于现行 规范的施工图设计的可靠性和可信性存在一定的问 题。基于上述考虑，为保证大桥在基础施工及运营期 的安全，有必要对大桥地基基础稳定性分析中的几个 关键性技术问题进行明确。 ２群桩基础与土体共同作用研究 ２．１现有分析方法及不足 级工程师 卅云飞 卯 ’男’江苏 海安人’毕业于。可’海大 硕士’高 的相同点是是把地基土视为各向同性弹性体，忽略了 土的非均质性等因素，区别在于布辛奈斯克解所针对 的应力应变关系都是在弹性半空间的表面，而明德林 解可以在弹性半空间体内任何一点。对此规范设计中 作了若干简化假定，但简化和假定主要是基于中小型 桩基的大量实践经验，实践表明，按此计算，与实际观 测报告相比，沉降明显偏大，设计相对偏于安全的，经 济型有待商榷。在规范推荐公式的基础上部分学者和 研究人员进行了修正和补充完善，取得了大量有益的 成果，比如Ｂｕｔｔｅｒｆｉｅｌｄ和Ｂａｎｅｒｊｅｅｔ 睬用明德林位移基 本解建立的弹性分析方法和Ｃｈｏｗｔ２］将理论ｔ—ｚ曲线［３１ 非线性分析法。对该跨江大桥基础等特大型群桩基 础，国内外的实例和经验均少，不能完全依赖于经验 的、简化的常规设计方法，有必要引入新的思路和分 析方法。 围绕大直径超长群桩基础的承载性能研究，国内 外学者和研究人员主要从模型实验和数值模拟分析 人手，在此基础上得到符合性验证结果并推荐相应的 计算方法。鉴于超大型项目如本文所提及的大桥，在 超厚软弱地层和复杂的水文地质条件下研究和工程 实例在国内外还很少见，故有必要对设计中常见的几 个问题进行明确。 ２．２基础设计中的主要问题 某跨江大桥两个主塔基础场地覆盖层厚度３００ｍ 左右，桩端位于地层分布相对较稳定及层厚较均匀的 沉积地层中，桩长１２０ｍ左右。基础的面积、体积及刚 度大，同时，群桩效应导致桩基的受力不均，桩基沉降 大，超大桩长的摩擦桩的受力机理非常复杂［４１。设计施 工中面临的主要的技术问题是超长大直径钻孔灌注 桩承载ｌ生能复杂，桩的受力特性、桩基沉降、轴力的分 配规律等；在不同时间效应、不同荷载工况下，复杂的 水文条件下群桩的受力、沉降更为复杂。 ３解决思路 １． 、・ 日 大型复杂群桩基础通常根据变形进行基础设计， 对此，有限元数值分析法提供了较好的解决思路和方 法，该采用数值模拟法，假设地基土水平成层，同一土 层均质，各项同性，土体为理想弹塑性材料，桩为线弹 性材料。将桩、土、承台一起作为脱离体分析其受力变 形，在求解中自然地体现它们之间的共同作用。从原 则上讲，能够模拟桩基的实际工作性能，如尺寸效应， 刚度影响荷载作用方式及大小。该法可以充分考虑土 层的应力应变非线性特性，土层的非均质性（成层 性），考虑复杂的荷载状态和复杂的桩土界面特性条 件，桩土响应的时间效应施工效应及群桩效应等。通 过计算，还可以给出每根桩、每个局部，在加荷的每个 阶段的受力变形状况。因此有必要通过三维非线性有 限元计算分析，可以更好地回答大桥桩基础设计中的 变形等稳定性分析问题，进一步论证设计的合理性， 更恰当地评价桩基础的安全性，或者提出实施中需要 注意的薄弱环节。目前有限元数值模拟方法的难点主 要在于桩、土体本构模型的合理选择和本构模型的参 数选取方面，文章对几个关键性问题进行相关探讨。 ４土体本构模型及接触面单元的选取 ４．１土体本构模型的选择 土体应力变形的非线性特性对能否准确模拟土 体及结构物的应力变形特性有重要影响，从而影响到 土与结构共同作用的性质，应当选择恰当的土体本构 模型对这种非线性特征进行描述。针对目前的设计方 案，采用非线性三维有限元方法考虑桥墩群桩、承台、 地基土的共同作用，模拟实际的荷载条件、边界条件、 土体的非线性特性，计算其受力变形。演算过程表明， 土的本构模型采用邓肯一张双曲线非线性弹性模型 吻合度较好。 在模型建立方面，应充分考虑结构刚度对基础变 形和受力的影响，将桩基础施工分为多个不同阶段进 器 行模拟，在模拟过程中应当考虑桩土接触问题、施工 的时间效应、结构的刚度和注浆加固等多种情形下基 础沉降和应力分布的规律。总体来讲，模型的建立以 省时、可行和可靠为目的。 ４．２接触面单元的选择 由于桩体材料为混凝土，混凝土模量远高于桩侧 土体的模量，在水平或竖直荷载作用下，桩土之间必 然存在位移，同时考虑到由于水下作业钻孔灌注桩采 用护壁泥浆必然对桩周土体产生扰动，桩周存在与桩 间土性质不同的区域，在此条件下，因此设置接触面 单元更符合实际边界条件。 接触面单元主要有两种型式：＠Ｇｏｏｄｍａｎ单元ｌ ５ｌ， 无厚度，形式简单，应用便捷，能够考虑到接触面变形 的非线性特质，缺陷是实际工程中两侧单元可能相互 嵌入，法向位移可能产生偏差，通用解决方案是受压 时弹性系数取大值，受拉时取小值，但是物理意义不 明确，实际使用时尚应结合水文条件、施工工况等条 件时应加以修正。②有厚度的薄单元，如Ｄｅｓａｉ单元ｆ６３ 等。Ｄｅｓａｉ单元概念与实际工况比较吻合，能够较好的 模拟结构面的错动、剪切等接触状态，更符合实际。但 是依据Ｄｅｓａｉ的建议薄单元的厚度，可取单元长度的 ０．０１～０．１取值，范围过大。实际应用中，取值过大，与 实际单元接近，失去模拟意义；取值过小，计算精度下 降。为了了解接触面单元厚度对计算的影响，应该对 几种不同厚度的接触面单元作对比计算，并对模型进 行相应修正。 ４．３计算参数的取值 数值模拟计算参数的合理确定关系到计算成果 的合理与否。模型主要涉及的桩土材料力学参数有： 桩体的弹性模量、泊松比、土体的弹性模量、泊松比、 粘聚力、内摩擦角等。其中桩体的参数较为明确，可以 有规范或试验值参考。但是土体弹性模量等取值难以 确定，对计算结果影响甚大。设计人员常常依据经验 公式，以勘察报告提供的压缩模量进行换算取值，人 为因素大，缺乏理论依据，导致不同取值下的沉降计 算结果相差过大。大比例模型试验精度高，试验后可 得到本构模型和相应参数，但是对于一般情况，鉴于 现场模拟难度大、费用高，故目前设计人员仍然多采 用室内试验数据如三轴仪固结排水试验、固结试验、 直剪试验等的结果近似确定，这就对勘察报告的质量 提出了很高的要求。由于室内试验与岩土的原始实际 受力条件通常不符，依据岩土工程勘察报告提供的三 轴固结排水曲线、固结排水剪试验数据和直剪试验数 据进行整理取值，模拟结果沉降将相当大。这是由于 室内试验由于土的扰动，无法保存原来的结构性，实 际土体在百米深处存在相当大的自重应力，取出后应 力释放。但如果试验时并未先加荷到原来的压力后再 开始逐级加荷，就会造成了压缩模量偏小。为此，室内 土工试验数据应尽可能还原土的原始状态，试验前可 先进行自重应力下的预固结等前期试验，可以需要利 用固结试验和直剪试验资料确定对应土层的参数。沉 降计算时仍应结合已有计算经验并参照其它类似工 一 程，实际使用加以数据修正。 有厚度接触面单元也需要强度指标，实际上弹性 模量和泊松比也是非线形的。考虑到桩基施工时对桩 周土的扰动，以及泥皮效应会使桩基及土的强度和模 旱 鲤兽 水流冲刷对桩基的影响，跨江大桥水文地质条件 复杂，随着桩基施工的进行导致桩周土体被水流反复 冲刷，影响桩基础的承载性能，鉴于深水环境下大型 桥桩基础从建设到施工、使用期河床冲刷等经验较 少，对此积累的冲刷条件下桩基础的承载性能研究亦 量不同程度的降低。对于接触面强度降低的问题。水 利大坝工程作过不少接触面变形试验，主要用大比例 模型结合大型直剪仪进行，模拟精度较高。故在大型 桥梁基础设计与三维模拟时，可参照类似水利工程中 少见。目前仍然多采用经验或半经验公式，但已有的 研究成果和资料表明其基本的规律是随着冲刷深度 的增加，基础存在回弹效应［７１，对应的存在临界深度， ＿ 的参数进行取值，实践证明效果较好。 对此，有待进一步深入研究。 ５施工过程的仿真 参考文献 ［１］Ｂｕｔｔｅｒｆｉｅｌｄ Ｒ，Ｂａｎｅ ｒｊｅｅ Ｐ Ｋ．Ｔｈｅ Ｐｒｏｂｌｅｍ ｏｆ Ｐｉｌｅ－Ｇｒｏｕｐ－Ｐｉｌｅ Ｃａｐ 某垮江大型桥梁基础的承台设计厚度大，一般为 Ｉｎｔｅｒａｃｔｉｏｎ．Ｇｅｏｔｅｃｈｎｉｑｕｅ，１９７１，２１（２）：１３５－１４２． ［２］２ Ｃｈｏｗ Ｙ Ｋ．Ａｎａｌｙｓｉｓ ｏｆ Ｖｅｒｔｉｃａｌｌｙ Ｌｏａｄｅｄ Ｐｉｌｅ Ｇｒｏｕｐｓ．Ｉｎｔｅｒｎａｔｉｏｎａｌ １２ｍ，远超一般中小规模项目承台，在计算模拟时应 Ｊｏｕｒｎａｌ ｆｏｒ Ｎｕｍｅｒｉｃａｌ ａｎｄ Ａｎａｌｙｔｉｃａｌ Ｍｅｔｈ－ｏｄｓ ｉｎ Ｇｅｏｍｅｃｈａｎｉｃｓ， 当将其分级，模拟时应随着荷载增加逐步增加单元。 １９８６，１０：５９—７２． 同时应模拟混凝土浇筑、硬化过程，浇筑初期，刚度强 ［３】Ｋｒａｆｔ Ｌ Ｍ，Ｒａｙ Ｒ Ｐ，Ｋａｇａｗａ Ｔ．Ｔｈｅｏｒｅｔｉｃａｌ ｔ－ｚ Ｃｕｒｖｅｓ．Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ 度均低，作为柔性材料，无法调整各桩之间的荷载分 Ｇｅｏｔｅｃｈｎｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｄｉｖｉｓｉｏｎ，ＡＳＣＥ，１９７８，１０４（ＧＴ１２）： 配，只能将荷载传到紧挨其下的桩体上。如果桩的布 １４６５—１４８８． 置是均匀的，荷载也是均匀的，那么各桩承担的荷载 【４］程哗．超长大直径钻孔灌注桩承载性能研究【Ｄ］．南京：东南大学， ２００５． 就大体一致。当承台混凝土强度和模量都达到设计值 【５］钱家欢，殷宗泽．土工原理与计算【Ｍ】．北京：水利电力出版社，１９９６． 后，承台将均匀下沉，承台所受的桩反力就变得不均 ［６】殷宗泽．土工数值计算［Ｍ］．北京：水利电力出版社，１９９４． 匀了，也就是承台将调整桩轴力的分布。因此，承台是 窗体底端 柔性还是刚性对桩群轴力的分布有显著影响。设计计 ［７］袁灯平．大直径超长桥梁桩基与群桩基础竖向承载特性及沉降控 算中应当模拟这种混凝土由软到硬的过程。 制研究【Ｄ］．上海：同济大学，２００７． （上接第１１２页） ５结论 ［２］王文正．公路隧道施工过程的数值模拟及ＡＮＳＹＳ实现［Ｊ］＿西部交 通科技，２００７ｆ４）． 【３】梁勇旗，刘丰铭．施工顺序对浅埋偏压黄土隧道稳定性的影响分析 通过对本隧道不同施工方法数值模拟，得出以下 ｆＪ】．地下空间与工程学报，２０１３（￥２）． 结论： 【４Ｊ郭建，王起才，等．浅埋暗挖黄土隧道施工过程数值模拟分析ｆＪ】．铁 ①两种不同的施工方法，其围岩沉降分布大致一 道建筑，２０１０（８）． 样，位移最大处均位于拱顶部分；其中三台阶七步法 ［５］刘赖．郑西客运专线大断面黄土隧道施工方法研究［Ｊ】．现代隧道技 开挖引起的洞身下沉量较ＣＤ法更大。 术，２００７（６）． ②隧道开挖顺序对支护结构内力的分布区域影 ［６】张金柱，郝文广．郑西大断面黄土隧道施工方法模拟分析［Ｊ】．隧道 建设，２００７（￥２）． 响较大，三台阶七步工法主要集中在拱顶处，ＣＤ工法 ［７］李宇翔．三车道公路隧道分部开挖及支护技术研究【Ｄ】．重庆：重庆 主要集中在边墙处。 交通大学，２０１３． 本文结合某隧道具体实际施工情况，对两种施工 ［８】Ｅｉｋｉ Ｍｉｚｕｎｏ，ｈｓｕｏ Ｈｉｒａｎｏ，Ｋｕｎｉａｋｉ Ｋａｔａｏｋａ．Ｅｘｃａｖａｔｉｎｇ ｔｕｎｎｅｌｓ ｏｆ ｌａｒｇｅ 方法有针对性的进行分析、计算和研究，提出合理的 ｃｒｏｓｓ－ｓｅｃｔｉｏｎ ｗｉｔｈ ａ ｔｈｉｎ ｏｖｅｒｂｕｒｄｅｎ ｂｅｎｅａｔｈ ｄｅｎｓｅｌｙ ｐｏｐｕｌａｔｅｄ ａｒｅａｓ（Ｋｏｎａｎ Ｔｕｎｎｅｌ ｏｎ Ｒｉｎｇ Ｒｏａｄ Ｎｏ．２ ｉｎ Ｙｏｋｏｈａｍａ Ｃｉｔｙ）． 施工方案，并且通过了实际工程的验证取得了良好的 Ｔｕｎｎｅｌｌｉｎｇ ａｎｄ Ｇｒｏｕｎｄ Ｃｏｎｄｉｔｉｏｎｓ，１９９４． 实践应用效果。 【９］李权ＡＮＳＹＳ在土木工程中的应用 ｂ京：人民邮电出版社，２００５． ［１０】黄志钦．浅埋暗挖法隧道综合加固施工技术［Ｊ】．福建交通科技， 参考文献 ２０１４ｆ４）． 【１】王晓州．大断面黄土隧道建设技术［Ｍ】．北京：中国铁道出版社，２００９