济南西客站市政配套工程站前广场基坑工程设计优化研究报告

济南西客站市政配套工程站前广场基坑工程设计优化研究研究报告委托单位济南市勘察测绘研究院完成单位山东大学完成地点济南完成时间2010年8月课题负责人徐帮树主要完成人徐帮树刘日成王育奎蔚立元公铭报告编写人刘日成、王育奎报告审核人徐帮树目录第一章绪论411研究依据412工程概况413场地工程环境514周边施工影响515工程地质与水文地质条件5151工程地质条件5152地形地貌及水文地质条件816支护体系（详见施工图）817降水设计（详见施工图）1018研究内容12第二章计算原理1421FLAC3D计算原理1422FLAC3D优点1523FLAC3D的本构模型1624FLAC3D计算中的问题1725渗流原理1826锚杆单元原理21第三章支护方案变形规律研究2431计算目的2432计算说明2433方案一2834方案二3635方案一和方案二对比分析4236方案三4337方案四5038方案三和方案四对比分析5539方案五57310本章小节64第四章降水设计方案优化研究6641降水方案设计概况（初步设计）6642研究内容6643计算方案6744计算参数及范围67441计算范围67442网格划分68443计算参数69444边界条件6945计算结果及分析70451每种工况下的降水漏斗曲线70452不同工况下的相同观测点地下水位值对比曲线74453不同工况下的渗流矢量与孔隙水压力云图7546本章小节82第五章总结83附表各种工况下的关键点地下水位值85附表1不设止水帷幕时各观测点地下水位值（M）85附表2止水帷幕深18M时各观测点地下水位值（M）87附表3止水帷幕深21M时各观测点地下水位值（M）89附表4止水帷幕深24M时各观测点地下水位值（M）91附表5止水帷幕深27M时各观测点地下水位值（M）93附表6止水帷幕深30M时各观测点地下水位值（M）95第一章绪论11研究依据（1）济南西客站站前广场岩土工程勘察报告（初勘），济南市勘察测绘研究院，2009年5月。（2）济南市西客站站前广场深基坑工程初步设计文件，济南市勘察测绘研究院，2009年8月。（3）济南市西客站站前广场深基坑降水方案设计优化研究合作意向书，济南，2010年4月。（4）济南西客站站前广场基坑工程方案设计优化研究建议书，济南市勘察测绘研究院和山东大学岩土与结构工程研究中心，2010年3月。（5）关于站前广场基坑工程方案设计优化研究建议书的意见，济南西客站市政配套工程设计总体组山东省建筑设计研究院，2010年4月。（6）济南西客站站前广场抽水试验报告，济南市勘察测绘研究院，2009年9月。（7）西客站片区场站一体化（站区配套设施）工程东广场基坑边坡支护，济南市勘察测绘研究院，2010年7月。（8）西客站站前广场基坑工程地下水控制，济南市勘察测绘研究院，2010年6月。12工程概况济南西客站是京沪高铁沿线五个始发站之一，济南西客站市政配套工程是济南市为配合高铁站正常使用而配套建设的工程。该项目包括东广场地下工程、过街地道工程、北综合体、南综合体、人防工程、轨道交通1、6号线土建预留工程等单体工程，共同形成站前广场，南北长约470M，东西长约346M，位于站西路以东，济西东路以南，站东路以西，站前路以北的围合区域。站前广场基坑工程挖深大部分13M，北区酒店综合楼双塔和南区高层办公楼区域挖深159M，预埋地铁1号线挖深1661M2141M，基坑支护范围约470M350M。根据建设单位开工安排，基坑开挖前南、北广场清除自然地面下约15M填土，地铁1线外伸段部分清除20M填土，填土清除范围为基坑开挖上口线外扩一倍基坑开挖深度，因此基坑边坡支护结构按挖深125M150M设计。本设计为东广场基坑工程的支护结构设计。13场地工程环境站前广场西侧紧邻京沪高铁主站房，主站房西侧为高铁路基。站前广场西侧预留济南轨道交通6号线地下结构外墙距高铁路基约为85M，站前广场西南角距拟建公交枢纽约为50M。其余各侧均为拆迁场地，无地下管线和既有建（构）筑物。14周边施工影响2011年3月31日济南轨道交通1、6号线预埋工程与站前广场同时达到±0000标高，意味济南轨道交通1、6号线预埋工程与站前广场同时建设。据了解，2010年5月公交枢纽开工，站前广场建设将与公交枢纽基坑同时施工。2010年8月高铁路基铺轨，自2010年8月站前广场建设期间需保证高铁路基“零沉降“。15工程地质与水文地质条件151工程地质条件勘区地层上部为第四系全新统河流冲积成因的软可塑粘性土、局部夹中、粗砂、卵石土，下部为第四系更新统山前冲洪积成因的粘性土，夹砂土及卵石土，钻探深度范围内地层可分为9层，自上而下分述如下填土分布于场地表层，分为素填土（层）和杂填土（1层）。素填土黄褐色，可塑，湿，含少量砖屑、灰渣。主要分布于该层中下部。1杂填土杂色，成分主要为建筑垃圾，砖块、灰渣等，结构松散，位于该层上部。层厚050620米，层底标高24263006米。粉质粘土褐黄色，局部黄褐色，可塑，湿，含铁锰氧化物，无摇振反应，刀切面较光滑，干强度和韧性中等。层厚040500米，层底深度120550米，层底标高21612511米。粉土褐黄色，中密，湿，含氧化铁斑点，摇振反应迅速，刀切面粗糙，干强度和韧性低。层厚030430米，层底深度250780米，层底标高22692849米。粉质粘土褐黄色、灰黄色，可塑，湿，含铁锰氧化物，无摇振反应，刀切面较光滑，干强度和韧性中等。层厚340990米，层底深度7301390米，层底标高17442156米。粉质粘土浅棕黄色，可塑，湿，含铁锰氧化物，少量姜石，无摇振反应，刀切面较光滑，干强度和韧性中等。该层内主要有4个亚层，分述如下1粉质粘土灰黄色，可塑，含铁锰氧化物、少量姜石，无摇振反应，刀切面较光滑，干强度和韧性中等；局部地段姜石含量可达1525,为粉质粘土混姜石。2砂层主要为中砂，浅棕黄色，中密，矿物成分为石英、长石，局部为中密粗砂、卵石土。砂层中局部地段为砂粒钙质胶结砂岩，钻探取样呈碎块或短柱状。3粘土浅棕黄色，可硬塑，湿，含铁锰氧化物，少量姜石，无摇振反应，刀切面光滑，干强度和韧性高。4泥钙质胶结砾岩浅棕黄色，呈泥钙质胶结，钻探岩芯呈碎块状或短柱状，零星分布。层厚11502080米，层底深度22003100米，层底标高083823米。粘土浅棕黄色，硬塑，湿，含铁锰氧化物及其结核，少量姜石，无摇振反应，刀切面光滑，干强度和韧性高。该层内主要有2个亚层，分述如下1粉质粘土浅棕黄色，硬可塑，湿，含铁锰氧化物，无摇振反应，刀切面较光滑，干强度和韧性较高。2砾岩浅棕黄色棕黄色，泥、钙质胶结，钻探岩芯呈硬块状。层厚1601000米，层底深度29003450米，层底标高448189米。粉质粘土浅棕黄色，硬可塑，湿，含铁锰氧化物，无摇振反应，刀切面较光滑，干强度和韧性中等。该层内主要有3个亚层，分述如下1粘土浅棕黄色，硬塑，湿，含铁锰氧化物，无摇振反应，刀切面光滑，干强度和韧性高。2细砂浅棕黄色，中密，矿物成分为石英、长石，局部含少量卵石。主要呈透镜体状零星分布。3砾岩棕黄色，泥、钙质胶结，钻探揭示呈碎块或短柱状。层厚7002200米，层底深度40005370米，层底标高23701000米。卵石、中砂夹粉质粘土该层以卵石为主，顶部局部分布中砂，中部夹棕黄色硬可塑粉质粘土、胶结砾岩。卵石土中密，其母岩成分为石灰岩，亚圆形，粒径26厘米，含量6080，混棕黄色粘性土，局部胶结为砾岩。1中砂棕黄色，中密，矿物成分为石英、长石，混少量卵石，局部为粗砂。2粉质粘土浅棕黄色，硬可塑，湿，含铁锰氧化物，无摇振反应，刀切面较光滑，干强度和韧性中等。主要呈薄层或透镜体分布于该层中部。层厚25803800米，层底深度76008500米，层底标高55364605米。辉长岩残积土黄绿色，密实，湿，原岩剧烈风化，结构全部破坏，具塑性，手捏呈土状。层厚200100米，层底深度80008600米，层底高程56055000米。强风化辉长岩黄绿色，密实，湿，矿物成分以普通辉石和斜长石为主。中粒结构。原岩强烈风化呈粗砾砂状和碎块状。该层未揭穿，最大揭露厚度1200米，最大揭露深度9200米。勘察报告提供相关参数如下表根据周边类似基坑施工经验，结合基坑规范规程及相关文献，考虑全封闭止水帷幕内土层的降水作用和勘察报告土层相关指标的实验方法，经与勘察单位协调沟通，将土层抗剪强度指标标准值（粘聚力C、内摩擦角）上调10，作为基坑支护设计力学参数152地形地貌及水文地质条件拟建场地位于黄河、小清河冲积平原的边缘相，场地以北有小清河，以东有腊山河。整个场地地形开阔，地势较平坦，勘探期间测得场地自然地面标高为29353104米。场区水文地质单元位于西郊玉符河隐伏冲积扇前缘砂、砂砾石富水区。地下水类型为第四系孔隙潜水。主要由大气降水和地下水渗流补给。勘探期间属平水期，在钻孔中测得地下水静止水位埋深340476米，相应标高25672650米，水位随季节性变化较大，变化幅度100200米。16支护体系（详见施工图）（1）北、南广场基坑支护11剖面挖深约125M，108放坡。采用分阶放坡的土钉墙支护型式，设置五道土钉，挂网喷砼；加强段采用预应力锚杆和土钉墙结合的复合土钉墙支护，设置二道土钉、三道预应力锚索，挂网喷砼；（2）北区酒店综合楼双塔和南区高层办公楼区域、预埋地铁1号线外伸端南北侧22剖面挖深约150M，108放坡。采用分阶放坡的土钉墙支护型式，设置六道土钉，挂网喷砼；加强段采用预应力锚杆和土钉墙结合的复合土钉墙支护，设置二道土钉、四道预应力锚索，挂网喷砼；（3）预埋地铁1号线外伸端东侧（横穿站东路，邻近临时道路）33剖面挖深约150M，105放坡。采用预应力锚杆及土钉墙结合的复合土钉墙支护，设置七道土钉、三道预应力锚索，挂网喷砼；（4）预埋地铁1号线站前广场内的基坑，在站前广场开挖深度125M的基础上，再下挖261M741M，采用10411坡率，挂网喷砼，见剖面设计44、55、66；（5）预埋地铁6号线站前广场内的结构，与站前广场同步建设，地铁6号线地下工程施工形成的一阶基坑东侧边坡将挖除，南、北侧将与站前南、北广场基坑边坡连接。保留先期施工一阶基坑南、北侧部分边坡，按115坡率挖除南、北两侧一阶基坑形成的边坡，挂网喷砼。（6）主要结构材料水泥P0425普通硅酸盐水泥钢筋HPB235强度标准值FYK235N/MM²HRB335强度标准值FYK335N/MM²HRB400强度标准值FYK400N/MM²钢绞线S152强度标准值FPTK1860N/MM²锚孔注浆材料水泥浆，强度不低于M20喷射混凝土强度等级C20（7）主要参数1）土钉土钉长度60150M，钻孔直径130MM，杆体材料20、28。取六根土钉做基本实验，获取土钉承载力设计值。2）预应力锚杆（索）锚杆长度120150M，钻孔直径150MM，杆体材料S152（钢绞线）。取三根预应力锚杆做基本实验和蠕变实验，获取锚杆承载力设计值。3）腰梁腰梁采用钢肩梁施工工艺，材料选用槽钢218A。4）面层土钉墙支护、放坡部分面层钢筋网按65250250布置，喷砼C20，厚度80MM。加强段和地铁1线外伸部分端部面层钢筋网按8200200布置，喷砼C20，厚度100MM。17降水设计（详见施工图）综合考虑场地工程地质、水文地质条件，基坑周边环境及基坑深度，依据基坑工程有关技术规范、规程，结合周边类似工程经验，本工程采用全封闭止水帷幕的截水方案，大口径管井坑内降水，坑外回灌实现基坑开挖和主体建设必须的地下水控制目标。171基坑设计分区及有关参数（1）本基坑分四个区域，（一）北广场基坑；（二）南广场基坑；（三）预埋地铁6号线与站前广场内的基坑；（四）站前广场内预埋地铁1号线基坑。（2）基坑侧壁安全等级为一级。（3）本工程为临时性工程，设计时限12个月。（4）基坑底边线距地下室外墙线按12M20M考虑。（5）基坑开挖边缘2M范围内无堆载，2M以外至一倍基坑开挖深度范围内按15KPA设计；若需超过以上荷载请通知设计修正。（6）地层综合渗透系数K12M/D。172地下水控制体系1721止水帷幕为避免水位下降给周边道路及环境造成不利影响，保证高铁路基“零“沉降，采用高压摆喷/旋喷止水帷幕，形成封闭式止水帷幕防渗板墙，构造有效的防护体系。（1）止水帷幕设计北区酒店综合体双塔区域、南区高层办公楼区域、地铁1线外伸部分的帷幕采用高压旋喷，桩径1100MM，间距800MM。桩顶位于自然地坪下40M，帷幕有效长度25M，自自然地平起算帷幕深度不小于25M，以封闭砂卵层，控制管涌、流砂；其他区域采用高压摆喷，摆角30度，间距1000MM，喷射帷幕墙厚度和搭接长度均不小于200MM。桩顶位于自然地坪下40M，帷幕有效长度20M，自自然地平起算帷幕深度不小于24M。（2）主要结构材料水泥PO425普硅水泥（3）主要参数水泥浆灌浆压力0506MPA、进浆量6080L/MIN、浆液比重150170G/CM3；或者控制水灰比0810。根据需要添加早强剂。高压水压力3640MPA、流量75L/MIN；压缩空气压力0708MPA、流量0812M3/MIN；三管提升速度812CM/MIN；摆角300；摆动速度810R/MIN；钻孔施工时钻孔偏斜率应控制在1以内；为保证止水帷幕施工效果，建议进行试喷试验，以有效验证、明确帷幕设计施工参数。1722降水方案场地地下水埋深较浅，地下水较丰富，在基坑开挖支护和基础施工过程中采用降水措施抽排地下水。（1）采用大口径管井降水井进行降水作业基坑周边布置降水井，井间距为20M；坑内设置疏干井，井间距为36M36M（4倍柱网距）。北广场周边布置降水井，共21眼，设计井深为20M。坑内设置疏干井，共31眼，北区酒店综合体塔楼区域（5、6、7疏干井）井深23M，其他均20M。南广场周边布置降水井，共23眼，设计井深均为20M。坑内设置疏干井，共32眼。南区高层办公楼区域（24、25、31、32疏干井）井深23M，其他均20M。预埋地铁6号线站前广场范围内区域，保留地铁6号线施工用井，按需降水，保证站前广场范围内地铁六号线建设需要。预埋地铁1号线共布置降水井23眼，其中18号井井深26M，923号井井深21M。降水井疏干井井径700MM，井管直径500MM，井管采用混凝土无砂滤水管，管壁外侧回填滤料。（2）基坑坑底、坑顶设置排水沟、集水盲沟、集水井，配合管井降水，控制地下水位。基坑底部排水沟按盲沟设置，宽深300300MM，集水井尺寸500500500MM，排水沟集水井均距坡脚300MM。（3）回灌措施保证帷幕外地下水位稳定为减少基坑降水引起帷幕外地下水绕流，导致地面不均匀沉降，止水帷幕外侧设置33眼回灌观测井；井深设计为130M，水平间距约30M，井径700MM，井管径500MM。（4）控制降水速度和强度，减少长期降水和绕渗对周边环境的不利影响。（5）做好雨季期间基坑坑顶地面硬化以及坑内、外积水的有组织排放；（6）疏干井避开主体建筑基础桩基位置，若位置与桩位冲突，降水井可在2M范围内移动。18研究内容（1）基坑复合土钉墙支护变形规律研究设计单位依据建筑基坑支护技术规程（JGJ12099），已对基坑支护设计内部稳定、外部稳定、局部抗拉和混凝土喷射面层进行了强度验算。西客站基坑工程是重要的临时工程，为合理评估基坑开挖引起的变形及地表沉降，确保基坑开挖对周围环境和建筑物的影响控制在允许范围之内。我们采用FLAC3D对放坡土钉墙、复合土钉墙设计方案进行验算。绘制地表沉降曲线，分析各设计剖面基坑开挖引起的地表沉降规律，沉降最大值及出现的位置。基坑开挖引起的地表沉降影响范围。绘制基坑边坡水平变形曲线，分析基坑侧向变形规律，找出各剖面基坑变形最大值。对一般支护设计剖面、支护设计剖面加强段进行对比分析，研究预应力锚索在加强控制基坑水平变形中的作用。研究坑底隆起变形规律。预测各设计剖面基坑开挖施工基坑隆起最大值及出现的位置。（2）降水方案设计优化研究根据济南西客站站前广场的基坑工程地质、水文地质、周边环境、地下管线、基坑设计深度等基础资料，确保高铁路基“零”沉降，开展降水方案设计优化研究根据设计规范和经验，初步提出几种止水帷幕、降水井设计方案及相应参数；针对不同设计方案，对不同止水帷幕深度、止水帷幕漏水情形，开展三维基坑开挖过程的三维地下水渗流数值模拟研究。对不同工况下计算得到地下水位变化趋势进行分析，评估降水效果，给出较优的基坑降水设计方案，确保高铁路基“零”沉降。第二章计算原理根据地质勘察报告、设计图纸和研究目标，我们选用FLAC3D开展数值模拟研究。三维连续体快速拉格朗日分析（FLAC3D）是由美国ITASCACONSULTINGGROUPINC开发的三维显式有限差分法程序，它可以模拟岩土或其他材料的三维力学行为。FLAC3D将计算区域划分为若干六面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。FLAC3D采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、流固耦合过程、大变形分析以及模拟施工过程等领域有其独到的优点。21FLAC3D计算原理1空间导数的有限差分近似在FLAC3D中采用了混合离散方法，区域被划分为常应变六面体单元的集合体，而在计算过程中，程序内部又将每个六面体分为以六面体角点为角点的常应变四面体的集合体，变量均在四面体上进行计算，六面体单元的应力、应变取值为其内四面体的体积加权平均。如图52所示四面体，节点编号为1到4，第N面表示与节点N相对的面，设其内任一点的速率分量为，则可由高斯公式可得IV21SJIVJIDN,其中V为四面体的体积，S为四面体的外表面，N为外表面的单位法向向量分量。2运动方程FLAC3D以节点为计算对象，将力和质量均集中在节点上，然后通过运动方程在时域内进行求解。节点运动方程可表示为如下形式22LILIMFTV其中为在T时刻L节点的在I方向的不平衡力分量，可由虚功原理导TFLI出。为L节点的集中质量，在分析静态问题时，采用虚拟质量以保证数值稳M定，而在分析动态问题时则采用实际的集中质量。图21四面体3阻尼力对于静态问题，FLAC3D在式3的不平衡力中加入了非粘性阻尼，阻尼力表示为23LILILIVSGNTFTF其中为阻尼系数，其默认值为08。4计算循环FLAC3D的计算循环如图53所示。22FLAC3D优点FLAC3D与有限元法相比有如下优点（1）FLAC3D采用混合离散方法来模拟材料的屈服或塑性流动行为，这种运动方程对每个节点由应力及外力利用虚功原理求节点不平衡力由不平衡力求节点速率本构方程对每个单元由节点速率求应变增量由应变增量求应力增量及总应力图22计算循环方法比有限元方法中通常采用的降阶积分更为合理。（2）FLAC3D利用动态的运动方程进行求解，即使问题在本质上是静力问题。这使得FLAC3D很容易模拟动态问题，如振动、失稳、大变形等。（3）FLAC3D采用显式方法进行求解，对显式法来说非线性本构关系与线性本构关系并无算法上的差别，对于已知的应变增量，可很方便地求出应力增量，并得到不平衡力，就同实际中的物理过程一样，可以跟踪系统的演化过程。而有限元采用的隐式解法对于非线性问题需进行迭代计算，比线性问题要耗费许多的计算时间，而且无法跟踪系统的演化过程。此外，显式法不形成刚度矩阵，每一步计算所需计算机内存很小，使用较少的计算机内存就可以模拟大量的单元，特别适于在微机上操作。而且在求解大变形过程中，因每一时步变形很小，可采用小变形本构关系，只需将各时步的变形叠加，即得到了大变形。这就避免了通常大变形问题中推导大变形本构关系及其应用中所遇到的麻烦，也使它的求解过程与小变形问题一样。23FLAC3D的本构模型1FLAC3D包含了10种弹塑性材料本构模型，即空单元模型。三种弹性模型各向同性、横向各向同性、正交各向异性。六种塑性模型DRUCKERPRAGER准则、摩尔库仑准则、应变硬化/软化模型、多节理模型、双线性应变硬化/软化多节理模型、修正的剑桥模型。每个单元可以有不同的材料模型或参数，材料参数可以为线性分布或随机分布。2FLAC3D有五种计算模式。静力模式。这是FLAC3D的默认模式，通过动态松弛方法得到静态解。动力模式。用户可以直接输入加速度、速度或应力波作为系统的边界条件或初始条件，边界可以是吸收边界和自由边界。动力计算可以与渗流问题相耦合。蠕变模式。有五种蠕变本构模型可供选择以模拟材料的应力应变时间关系MAXWELL模型、双指数模型、参考蠕变模型、粘塑性模型、脆盐模型。渗流模式。可以模拟地下水流、孔隙压力耗散以及可变形孔隙介质与其间的粘性流体的耦合。渗流服从各向同性达西定律，流体和孔隙介质均被看作可变形体。考虑非稳定流，将稳定流看作是非稳定流的特例。边界条件可以是固定孔隙压力或恒定流，可以模拟水源或深井。渗流计算可以与静力、动力或温度计算耦合，也可以单独计算。温度模式。可以模拟材料中的瞬态热传导以及温度应力。温度计算可以与静力、动力或渗流计算耦合，也可单独计算。3FLAC3D可以模拟多种结构形式。对于通常的岩体、土体或其他材料实体，用八节点六面体单元模拟。FLAC3D包含有四种结构单元梁单元、锚单元、桩单元、壳单元。可用来模拟岩土工程中的人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩等。FLAC3D的网格中可以有界面，这种界面将计算网格分割为若干部分，界面两边的网格可以分离，也可以发生滑动，因此，界面可以模拟节理、断层或虚拟的物理边界。4FLAC3D可以有多种边界条件。边界方位可以任意变化，边界条件可以是速度边界、应力边界，单元内部可以给定初始应力，节点可以给定初始位移、速度等，还可以给定地下水位以计算有效应力，所有给定量都可以具有空间梯度分布。5FLAC3D具有强大的内嵌程序语言FISH，使得用户可以定义新的变量或函数，以适应用户的特殊需要。24FLAC3D计算中的问题（1）地下工程中的模型建立地下工程的数值模拟中一个非常突出的问题就是单元的离散化比较困难。由于地下工程往往体形比较复杂，并且形状大小不一，使得网格的划分成为整个分析过程中的瓶颈。FLAC3D具有强大的自动三维网格生成器，内部定义了多种基本单元形态，用户还可以用FISH自定义单元形态，通过组合基本单元，可以生成非常复杂的三维网格，比如交叉隧洞段。大大降低了地下工程中的网格划分的难度，加快了模型建立的速度。（2）边界条件的处理应力边界条件FLAC3D中默认情况下，边界为自由表面（自由表面是一种特殊的应力边界条件）使用其内部命令，可以给出任一边界或部分边界的外力和应力情况。在地球坐标系（X、Y、Z）中，应力为SXX、SYY、SZZ、SXY、SXZ、SYZ，向量坐标中（D、S、N）中通过法向力、切向力给出应力。位移边界条件FLAC3D中不能直接控制位移。实际上，位移并不参与FLAC的计算。FLAC通过速度来控制位移。这里所说的速度与物理中的速度概念有一些差别。它实际上是单位计算时步（STEP）的位移平均值。例如位移为D，计算时步为N，则速度VD/N。控制位移为零，只需将速度取得很小，计算时步取得很大即可。在FLAC3D中，可以使用内部命令设定边界的位移条件。（3）初始地应力的模拟在进行开挖或建设前，岩体已经处在一种应力状态下。这种初始状态，就是岩体中的初始应力场。初始应力场主要由岩体自重和地质构造力产生。初始应力场问题很复杂，因为构造应力常常分布极不均匀，而费用昂贵的现场地应力测量只能给出计算范围中少数几个点的地应力值，很难给出准确的应力场。多数情况下，这种初始条件是没有具体的资料。这时，使用FLAC3D可以通过现有的资料来模拟初始应力场。初始的自重应力场，可以先在自重条件下计算至平衡，然后再将所有的节点位移置零。这时在单元中就已经加了初始的自重应力场。对于初始的构造应力场，需要利用FLAC3D中内嵌的FISH语言编写程序来添加。25渗流原理FLAC3D可以模拟多孔介质中的流体流动，比如地下水在土体中的渗流问题。FLAC3D既可以单独进行流体计算，只考虑渗流的作用，也可以将流体计算与力学计算进行耦合，也就是常说的流固耦合计算。比如土体的固结，就是一种典型的流固耦合现象，在土体固结过程中超孔隙水压力的逐渐消散导致了土体发生沉降，在这个过程包含两种力学效应（1）孔隙水压力的改变导致了有效应力的改变，从而影响土体的力学性能，如有效应力的减小可能使土体达到塑性屈服。（2）土体中的流体会对土体体积的改变产生反作用，表现为流体孔压的变化。FLAC3D具有强大的渗流计算功能，可以解决完全饱和及有地下水变化的渗流问题。对于地下水问题，FLAC3D认为地下水位以上的孔压为零，且不考虑气相的作用，这种近似方法对于可忽略毛细作用的材料是适用的。FLAC3D模拟多孔介质如土体中的流体流动时，既可以单独进行流体计算，只考虑渗流的作用，也可以将流体计算与力学计算进行耦合，即常说的流固耦合计算。流固耦合计算考虑孔隙水压力消散引起的岩土体位移变化，这一过程包含两种力学效果第一，孔隙水压力的变化引起体积应变的变化，进而影响有效应力；第二，应变的改变也会影响孔隙水压力。该软件具有强大的渗流计算功能，可以解决完全饱和及有地下水位变化的渗流问题。FLAC3D中的渗流计算具有以下特点1针对不同材料的渗流特点，提供了三种渗流模型各向同性、各向异性及不透水模型；2不同的单元可以赋予不同的渗流模型和渗流参数；3提供了丰富而又实用的流体边界条件，包括流体压力、涌入量、渗漏量、不可渗透边界、抽水井、点源或体积源等。4计算完全饱和岩土体中的渗流问题，可以采用显式差分法或者隐式差分法；5渗流问题可以和固体力学问题、热问题进行耦合；6流体和固体的耦合程度依赖于固体颗粒骨架的压缩程度，用BIOT系数表示颗粒的可压缩程度。渗流计算的控制微分方程主要有输运方程、质量守恒方程方程和本构方程。流体的运动用DARCY定律来描述，对于均质、各向同性固体和流体密度是常数的情况，这个方程具有如下形式24,I1,23IILFJQKSPXGL式中QI是渗透流量,P是孔隙压力，是介质关于饱和度S的相对机动系KS数，F是流体密度，GI是重力的三个分量。对于小变形，流体质量守恒方程为25TQVI,式中QV是流体单位体积源强度1/S，是单位体积孔隙介质的流体体积变化量。流体含量的变化与孔隙压力P、饱和度S、体积应变E和温度T的改变有关，孔隙流体的本构方程为26TTSTNM11式中M是BIOT模量N/M2，N是孔隙率，是BIOT系数，是不排水导热系数1/，用此来考虑流体和颗粒的热膨胀。FLAC3D在分析含有孔隙水压力的问题时，根据是否设置流体计算，有渗流模式和无渗流模式两种计算模式。在无渗流模式下，也可以在节点上设置孔隙水压力，但其值保持不变，固体单元的屈服判断由有效应力决定。在渗流模式中，可以进行瞬态渗流计算，孔隙水压力随着浸润线的改变而改变；还可以进行有效应力计算、不排水计算以及完全流固耦合计算。在完全流固耦合情况下，孔隙水压力的改变会产生力学变形，同时体积应变又会导致孔隙水压力的改变。流体模型设置好以后就可以进行参数的赋值，FLAC3D渗流计算中涉及的参数包括单元参数如渗透系数、孔隙率、流体密度和比奥系数，节点参数如流体体积模量、饱和度、抗拉强度和比奥模量。考虑颗粒压缩的情况下，选择比奥模量和比奥系数；而对于颗粒不可压缩的岩土体，则选用流体体积模量和孔隙率。对于理想多孔介质，比奥模量M与流体体积模量KF的关系为271/FFKN其中，K为固体体积模量，而对于室温下的纯水而言KF2109PA。这样，对于不可压缩的情况1，有下式28/FMKNFLAC3D默认情况下模型边界都是不透水边界，即边界上节点与外界没有流体交换，其孔压值可以自由变化。设置孔隙压力固定时表示透水边界条件，流体可以流入流出模型边界。对于特定的流体问题，在计算之前首先要进行分析，考察该问题的一些指标，如时标、施加扰动的属性及流固刚度比等，以便于采用合理的分析方法。下面分别进行介绍。时标的概念类似于时间，在流固耦合分析过程中需要考虑流体计算和力学计算两个时标。如果流体进程的时标和问题需要分析的时标的差别很大，则有可能使用简化的不耦合分析方法。力学过程和流体扩撒过程的特征时间分别定义为294/3MCCUTLKG2102FCT其中，LC为特征长度模型体积与表面积之比，KU是不排水体积模量，G为剪切模量，是密度，C是扩撒率。扩撒率定义为渗透系数与储水系数的比值，见下式2112/1/4/3WPKMCNGLKG饱和渗流模式非饱和渗流模式流固耦合模式扰动是指流固耦合问题中引起系统平衡状态改变的外界条件，包括流体边界条件如孔隙水压力的改变和力学边界条件如荷载的变化。如果问题中的扰动仅仅是由于孔隙水压力的改变引起的，那么流体进程和力学进程可以不耦合；若是由于力学扰动引起的，则二者耦合程度要考虑流固刚度比的影响。流固刚度比是指流体模量和固体模量之间的比值，定义为21224/3KMRKG如果固体的刚度很大，或流体具有高压缩性，则得出RK较小，称之为相对刚性骨架，可以不进行耦合计算；如果固体模量较小，流体不可压缩，则RK较大，称之为相对柔性骨架，这时需要进行耦合分析。26锚杆单元原理锚单元的轴向行为特性采用一维模型描述，轴向刚度K为LAEK213式中A为锚单元的横截面积，E为弹性模量，L为锚单元的长度。由轴向位移增量UT，计算轴向力增量FTFTKUT214其中，UT（U1BU1A）T1（U2BU2A）T2（U3BU3A）T3。式中UIB，UIA为节点位移，I1，2，3；A,B分别为节点A，B；T1，T2，T3分别为锚单元轴线方向的方向余弦。锚单元可以指定其拉伸屈服强度FT和压缩屈服强度FC，单元轴力不能超过强度极限，如图23所示。图23锚单元轴向力变形特性曲线采用弹簧滑块系统来描述锚杆（索）、砂浆及岩体直接的相互作用关系，如图24所示。该模型反映了锚杆（索）与岩体界面（浆体）直接的剪切特性。界面直接可能产生的最大剪切力取决于浆体的剪切刚度，单位厚度浆体的剪切刚度KG为/21LNDTGKG215式中G为浆体的剪切模量；T为浆体环的厚度；D为钢筋索的直径。浆体界面的剪切力G为/21LNU/TTGG216式中U为浆体与岩体间的相对位移。浆体单位长度所能承受的最大剪应力为GMGSPCLFTANAX217式中CG为浆体的粘结强度，G'为浆体摩擦角，PG为浆体与锚单元或岩体的实际周长，M为有效围压。界面材料特性假定为理想弹塑性，采用摩尔库仑准则作为屈服准则。每个锚单元均允许沿轴向产生变形并发生屈服。若锚单元两节点与网格节点重合，锚单元与实体单元联结成整体，不产生相对位移；锚单元节点与网格节点不重合时，锚单元与实体元之间产生相对位移，其位移大小与界面模型参数有关。锚杆与岩体间的界面可以发生剪切屈服、产生滑动直至拉拔破坏。图24锚杆单元中的界面力学模型第三章支护方案变形规律研究31计算目的综合考虑场地工程地质、水文地质条件，基坑周边环境及基坑深度，依据基坑工程有关技术规范、规程，结合周边类似工程经验，拟选择放坡和复合土钉墙支护相结合的形式，既充分利用现有场地尽量放坡，加快先期施工进度，又通过土钉墙、复合土钉墙支护增加边坡刚度。采用分阶放坡土钉墙结构完成东广场基坑边坡支护；采用复合土钉墙承担预埋地铁1号线站前广场东延部分端部基坑边坡支护；同时，考虑基坑空间效应，每隔60M左右设置基坑边坡加强段，增强边坡刚度，通过预应力锚索控制边坡变形。设计单位依据建筑基坑支护技术规程（JGJ12099），已对基坑支护设计内部稳定、外部稳定、局部抗拉和混凝土喷射面层进行了强度验算。西客站基坑工程是重要的临时工程，为合理评估基坑开挖引起的水平变形及地表沉降，确保基坑开挖对周围环境和建筑物的影响控制在允许范围之内。我们采用FLAC3D对放坡土钉墙、复合土钉墙设计方案变形规律进行研究。32计算说明1计算方案编号为了在报告中叙述方便，对各计算工况进行编号。各计算工况编号如表31。表31计算方案编号方案编号对应的设计剖面主要设计参数方案一11支护单元挖深约125M，108放坡。采用分阶放坡的土钉墙支护型式，设置五道土钉，挂网喷砼。方案二11支护单元加强段加强段采用预应力锚杆和土钉墙结合的复合土钉墙支护，设置二道土钉、三道预应力锚索，挂网喷砼。方案三22支护单元挖深约150M，108放坡。采用分阶放坡的土钉墙支护型式，设置六道土钉，挂网喷砼。方案四22支护单元加强段加强段采用预应力锚杆和土钉墙结合的复合土钉墙支护，设置二道土钉、四道预应力锚索，挂网喷砼。方案五33支护单元挖深约150M，105放坡。采用预应力锚杆及土钉墙结合的复合土钉墙支护，设置七道土钉、三道预应力锚索，挂网喷砼。（2）建模说明工程地质条件复杂，地层分布厚度不均匀。为方便建模工作，做适当的简化处理是允许的。假定地质分层均匀，每层厚度不变。地层土体采用MOHRCOULOMB弹塑性本构模型，用六结点实体单元模拟。喷层面板采用SHELL壳单元模拟，线弹性本构模型。土钉和锚杆采用CABLE杆单元模拟，线弹性本构模型。分布开挖，分步施作土钉、锚杆支护。每层开挖至锚杆位置向下05M，再施工土钉或锚杆。开挖顺序示意如图31。图31开挖与支护顺序示意图（3）初始地应力基坑工程埋深较浅，一般不考虑构造地应力的影响。初始地应力只考虑土体自重应力场。采取止水帷幕，基坑周边范围的土体地下水位随排水、开挖等施工过程，水位下降较多。基坑变形分析假定不受地下水位的影响，计算土体取饱和重度。（4）物理力学参数根据2009年5月济南西客站站前广场岩土工程勘察报告、2010年7西客站片区场站一体化站区配套设施工程东广场基坑边坡支护设计选取计算用土体物理力学参数，如表32。勘察报告土体力学参数没有提供变形模量、泊松比值。变形模量由压缩模量换算得到，如下式（31）SE0（32）1/21式中为变形模量（MPA）；为压缩模量（MPA）；为泊松比。0ES上述理论关系的与实测值不相符，实测的变形模量往往会大于相应的压缩模量，甚至达到数倍。究其原因，除了由于地基土并不是理想弹性体，其成层性、孔隙率和含水等因素对弹性性质有显著影响外，在室内试验中不可避免地存在土样扰动的问题，破坏了土的天然结构，使得测得的压缩模量偏小（梁发云，2004）。根据统计资料，值可能是值的几倍，一般说来，土愈坚硬则倍数愈大，0ESE而软土的值与比较接近（土力学，2005）。通常取，取2535SSE0倍。本文取3倍，得到计算用变形模量，如表32。表32土体物理力学参数表序号土类型土层厚M压缩模量ES12MPA变形模量MPA泊松比饱和容重KN/M3粘聚力KPA内摩擦角º1素填土1543913170271910102粉质粘土269020700301932631943粉土291127330281912572344粉质粘土557717310321962502195粉质粘土116832049030194248216粘土1088026400301963422497粉质粘土1210113033030194457229土钉采用HRB335钢筋，强度标准值FYK335N/MM²。锚杆钻孔注浆材料为水泥浆，强度不低于M20。计算用的土钉力学参数如表33。HRB335屈服强度335MPA、抗拉强度445MPA、断后伸长率17。弹性模量200GPA，泊松比取03。浆体单位长度粘聚力、浆体单位长度上水泥浆刚度参考FLAC3D使用手册（EXAMPLEAPPLICATIONEXCAVATIONANDSUPPORTFORASHALLOWTUNNEL，FLAC3D）。表33土钉物理力学参数序号土钉型号密度KG/M3弹性模量EGPA横截面积AMM2浆体粘聚力CGN/M浆体摩擦角G浆体刚度GKN/M/M浆体周长GPMM压缩强度CFKN拉伸强度TKN120780020031422105301751074082105261398222878002106158021053017510740822062927403锚杆采用预应力锚索，锚杆长度120150M，钻孔直径150MM，杆体材料S152（钢绞线）。计算用的锚杆力学参数如表34。弹性模量取195GPA，泊松比取03。钢绞线S152，强度标准值FPTK1860N/MM²。表34预就历程锚索物理力学参数序号土钉型号密度KG/M3弹性模量EGPA横截面积AMM2浆体粘聚力CGN/M浆体摩擦角G浆体刚度GKN/M/M浆体周长GPMM压缩强度CFKN拉伸强度KNT1自由段78001951813711147133735337352锚固段7800195181372105301751074713373533735土钉墙支护、放坡部分面层钢筋网按65250250布置，喷砼C20，厚度80MM。加强段和地铁1线外伸部分端部面层钢筋网按8200200布置，喷砼C20，厚度100MM。面层采用壳单元模拟，计算用的喷射混凝土面板力学参数如表35。表35喷射混凝土面板物理力学参数序号面板类型密度KG/M3弹性模量GPA泊松比厚度MM1一般25002603802加强25002603100（5）边界条件假定沿基坑周边走向无限长，按平面应变计算。模型左右边界、前后边界施工法向约束，底边界施加固定约束，如图32。图32模型边界约束条件33方案一331建模11支护单元剖面图，如图33。土钉设计参数，如表36。图35图37为开挖计算有限差分网格，最大开挖深度H125M，考虑基坑土质情况的开挖影响范围，计算深度和宽度为最大开挖深度的（23）H，即模型整体高度取30M，整体宽度为130M。土钉设计按梅花形布置，水平间距为15M。土钉间排距不变，按梅花形布置、或矩形布置，对基坑边坡整体变形影响不大。为简化计算，土钉按矩形网格布置，计算模型沿基坑周边走向方向取15M，如图34。图3311支护单元剖面图表3611支护剖面段土钉设计参数土钉编号杆体类型长度MM水平间距MM钻孔直径MMTD112060001500130TD212090001500130TD312090001500130TD412090001500130TD512090001500130图3411支护单元锚杆立面图图35计算网格图36计算网格和支护锚杆图37计算网格局部透视图232计算结果及分析提取地表沉降位移、边坡表面侧向水平位移和坑底隆起变形值，监测点位置如图38。图38变形监测点示意图（地表、边坡、坑底）方案一合位移云图及矢量图，如图39图314。由图可见，随基坑开挖基坑坑底和边坡有向上、向内反弹变形的趋势。最大变形位置出现在坑底，坑顶变形较小。地表沉降范围较小。图39方案一合位移云图及矢量图第一步开挖图310方案一合位移云图及矢量图第二步开挖图311方案一合位移云图及矢量图第三步开挖图312方案一合位移云图及矢量图第四步开挖图313方案一合位移云图及矢量图第五步开挖图314方案一合位移云图及矢量图第六步开挖方案一地面沉降曲线如图315。地表沉降曲线形如降水漏斗形。坑顶边缘出现变形向上反弹，向上位移值为3078MM。最大沉降值为400MM，出现在离基坑边缘约20M的位置。2MM影响范围在基坑外40M，1MM影响范围在基坑外52M。方案一坡面水平变形曲线，如图316。按108比例放坡，坡角较小，基坑开挖坡面出形整体向上反弹趋势。整个坡面变形是上小下大，最大变形出现中下部。坑顶向基坑外变形389MM。最大变形出在基坑深度12M处，接近坑脚位置，最大水平水移值3250MM。方案一坑底隆起曲线，如图317。隆起曲线形如抛物线。因为受到基坑外土体的约束作用，最小值出现在坑脚位置，隆起值为8228MM。最大值出现在距坡脚约12M处，最大值为11225MM。方案一地表沉降曲线1000500000500100015002000250030003500000200040006000800010000到坑顶水平距离（M）地表沉降（MM）图315方案一地表沉降曲线方案一基坑坡面变形曲线14001200100080060040020000010000001000200030004000水平位移（MM）深度（M）图316方案一坡面水平变形曲线方案一坑底向上隆起曲线60708090100110120000500100015002000250030003500到坑脚距离（M）坑底隆起（MM）图317方案一坑底向上降起曲线233土钉轴力结果及分析基坑开挖完成后，土钉轴力如图318。土钉轴力形如枣核形，中间大两端小。最大轴力出现下排土钉的中间部位，达到2464KN，远小于土钉拉伸屈服强度。图318方案一土钉轴力图34方案二341建模11支护单元加强段剖面图，如图319。土钉设计参数，如表37。预应力锚索设计参数如表38。土钉、锚索按梅花形布置，计算模型沿基坑周边走向方向取30M，如图320。计算网格如图321、图322。图31911支护单元加强段剖面图表3711支护剖面加强段土钉设计参数土钉编号杆体类型长度MM水平间距MM钻孔直径MMTD112060001500130TD212090001500130表3811支护剖面加强段预应力锚索设计参数土钉编号杆体类型水平间距MM自由段长度MM锚固段长度MM预应力锁定值KN钻孔直径MMMG12S152300060009000160150MG22S152300060009000160150MG32S152300060009000160150图32011支护单元加强段锚杆立面图图321计算网格图32211支护单元加强段计算网格和支护锚杆图32311支护单元加强段计算网格局部透视图232计算结果及分析提取地表沉降位移、边坡表面侧向水平位移和坑底隆起变形值，监测点位置如图38。方案二开挖完成后最终合位移云图及矢量图，如图324。由图可见，随基坑开挖基坑坑底和边坡有向上、向内反弹变形的趋势。图324方案二开挖完成后基坑合位移云图及矢量图方案二地面沉降曲线如图325。地表沉降曲线形如降水漏斗形。坑顶边缘出现变形向上反弹，向上位移值为3182MM。最大沉降值为410MM，出现在离基坑边缘约20M的位置。方案二坡面水平变形曲线，如图326。按108比例放坡，坡角较小，基坑开挖坡面出形整体向上反弹趋势。整个坡面变形是上小下大，最大变形出现中下部。坑顶向基坑外变形548MM。最大变形出在基坑深度12M处，接近坑脚的位置，最大水平水移值3089MM。方案二坑底隆起曲线，如图327。隆起曲线形如抛物线，最小值出现在坑脚位置，是因为受到基坑外土体的约束作用，隆起值为8307MM。最大值出现在距坡脚约12M处，最大值为11203MM。方案二地表沉降曲线1000500000500100015002000250030003500000200040006000800010000到坑顶边缘距离（M）地表沉降（MM）图325方案二地表沉降曲线方案二基坑坡面变形曲线1400120010008006004002000001000500000500100015002000250030003500水平位移（MM）深度（M）图326方案二坡面水平变形曲线方案二坑底向上反弹曲线60708090100110120000500100015002000250030003500到坑脚水平距离（M）坑底隆起（MM）图327方案二坑底向上降起曲线233土钉轴力结果及分析基坑开挖完成后，土钉、锚索轴力如图328。土钉轴力形如枣核形，中间大两端小。锚索预应力设计值为160KN，基坑开挖完成后有应力损失，最大值为1377KN。图328方案二土钉轴力图35方案一和方案二对比分析为分析方案一和方案二控制基坑开挖变形方面的效果，把方案一和方案二变形数据绘制在同一个曲线上，地表沉降、坡面变形、坑底隆起曲线，如图329图331。由图329可见，11剖面支护加强段预应力锚索减少基坑水平变形效果较明显，深度10M位置减小变形量802，深度8M位置减小变形量2000，深度6M位置减小变形量3478。由图330、图331可见，两条曲线几乎重合，两种方案对控制地表沉降及坑底隆起效果差别不大。基坑坡面水平变形曲线14001200100080060040020000010000001000200030004000水平位移（MM）深度（M）11剖面段11剖面加强段图329基坑坡面水平变形曲线（方案一、方案二）坑底向上反弹曲线60708090100110120000500100015002000250030003500到坑脚距离（M）坑底隆起（MM）11剖面段11剖面加强段图330坑底向上反弹曲线（方案一、方案二）地表沉降曲线1000500000500100015002000250030003500000200040006000800010000到坑顶边缘距离（M）地表沉降（MM）11剖面段11剖面加强段图331地表沉降曲线（方案一、方案二）36方案三361建模22支护单元剖面图，如图332。土钉设计参数，如表39。图333图335为开挖计算有限元差分网格，最大开挖深度H15M，考虑基坑土质情况的开挖影响范围，计算深度和宽度为最大开挖深度的（23）H，即模型整体高度取30M，整体宽度为130M。土钉设计按梅花形布置，水平间距为15M。土钉间排距不变，按梅花形布置、或矩形布置，对基坑边坡整体变形影响不大。为简化计算，土钉按矩形网格布置，计算模型沿基坑周边走向方向取15M，如图333。图33222支护单元剖面图表3922支护剖面段土钉设计参数土钉编号杆体类型长度MM水平间距MM钻孔直径MMTD112090001500130TD2120120001500130TD3120120001500130TD4120150001500130TD5120120001500130TD612090001500130图33322支护单元锚杆立面图图334计算网格图335计算网格和支护锚杆图336计算网格局部透视图362计算结果及分析提取地表沉降位移、边坡表面侧向水平位移和坑底隆起变形值，监测点位置如图337。图337变形观测点示意图（地表、边坡、坑底）方案三合位移云图及矢量图，如图338。由图可见，随基坑开挖基坑坑底和边坡有向上、向内反弹变形的趋势。图338方案三基坑开挖完成后合位移云图及矢量图方案三地面沉降曲线如图339。地表沉降曲线形如降水漏斗形。坑顶边缘出现变形向上反弹，向上位移值为2886MM。最大沉降值为480MM，出现在离基坑边缘约178M的位置。方案三坡面水平变形曲线，如图340。按108比例放坡，坡角较小，基坑开挖坡面出形整体向上反弹趋势。整个坡面变形是上小下大，最大变形出现中下部。坑顶向基坑外变形274MM。最大变形出在基坑深度125M处，接近坑脚的位置，最大水平水移值4059MM。方案三坑底隆起曲线，如图341。隆起曲线形如抛物线，最小值出现在坑脚位置，是因为受到基坑外土体的约束作用，隆起值为8543MM。最大值出现在距坡脚约12M处，最大值为10792MM。方案三地表沉降曲线1000500000500100015002000250030003500000200040006000800010000到坑顶边缘距离（M）地表沉降（MM）图339方案三地表沉降曲线方案三基坑坡面变形曲线160014001200100080060040020000010000001000200030004000500