一、盾构隧道结构计算模型1、惯用法(自由圆环变形法)惯用法的想法早在1960年就提出了,在日本国内得到了广泛的应用。惯用法假设管片环是弯曲刚度均匀的环,不考虑管片接头部分的柔性特征和弯曲刚度下降,管片截面具有同样刚度,并且弯曲刚度均匀的方法。这种方法计算出的管片环变形量偏小,导致在软弱地基中计算出的管片截面内力过小,而在良好地基条件下计算出的内力又过大。地层反力假设仅在水平方向上下45°范围内按三角形规律分布,这种模型可以计算出解析解。P0土压qe2水压qw24545RcR0水压Pw1土压Pe1g垂直抗力Pe2+Pw2自重反力Pg=πg水压qw2土压qe2(水土分离算法)HwHkδ地基抗力kδ地基抗力qw1qe1qe1qw12、修正惯用法在采用惯用法的60年代,怎样评价错缝拼装效应是一个问题。如果错缝拼装管片,可弥补管片接头存在造成的刚度下降。于是,在对带有螺栓接头的管片环进行多次核对研究时,首次引入了η-ξ对错缝拼装的衬砌进行内力计算,即为修正惯用法。该法将衬砌视为具有刚度ηEI的均质圆环,将计算出的弯矩增大即(1+ξ)M,得到管片处的弯矩;将求出的弯矩减少即(1-ξ)M,得到接头处的弯矩。其中η称为弯曲刚度有效率,ξ称为弯矩增加率,它为传递给邻环的弯矩与计算弯矩之比。管片接头由于存在一些铰的作用,所以可以认为弯矩并不是全部经由管片接头传递,其一部分是利用环接头的剪切阻力传递给错缝拼装起来的邻接管片。M隧道纵向M2/2M1M2/2M2/2M2/2MM1+M2BB/2接头传递弯矩示意图二、管片计算荷载的确定1、荷载的分类衬砌设计所考虑的各种荷载,应根据不同的地质条件和设计方法进行假定并根据隧道的用途加以考虑。衬砌设计所考虑的各种荷载见表所示。衬砌设计荷载分类表荷载分类荷载名称永久荷载结构自重地表超载土压水压地层抗力设备重量可变荷载基本可变荷载地面车辆荷载地面车辆荷载引起的侧向土压力隧道内部车辆行人等引起的荷载其它可变荷载施工荷载(设备运输、施工机具及人员、盾构推进、压注浆等)偶然荷载地震荷载2、计算断面选择埋深最大断面埋深最小断面埋深一般断面水位3、水土压力计算对于粘性土层,如西安地铁黄土地层、成都地铁二号线膨胀土地层等,应采用水土压力合算的方式进行荷载计算。此时,地下水位以上地层荷载用湿容重计算,地下水位以下用饱和容重计算。对于透水性较好的砂性地层,如西安地铁粗砂、中砂地层,成都地铁卵石土地层等,应采用水土压力分算的方式进行荷载计算。此时地下水位以上地层荷载用湿容重计算,地下水位以下用浮容重计算。水土压力合算与分算,主要影响管片结构侧向荷载。一般水土分算时侧向压力更大。4、松弛土压力将垂直土压力作为作用于衬砌顶部的均布荷载来考虑。其大小宜根据隧道的覆土厚度、隧道的断面形式、外径和围岩条件等来决定。考虑长期作用于隧道上的土压力时,如果覆土厚度小于隧道外径,一般不考虑地基的拱效应而采用总覆土压力。但当覆土厚度大于隧道外径时,地基中产生拱效应的可能性比较大,可以考虑在计算时采用松弛土压力,一般采用泰沙基公式计算。HP0h0B1实际松动范围计算松动范围π/4+φ/2π/4-φ/2R图2-5松弛土压力示意图1010/.tan0/.tan0110)1(tan)/1(BHKBHKePeKBcBh1010/.tan0/.tan011)1(tan)/(BHKBHKvePeKBcB)22/4/cot(01RB式中:—土的松动高度—0h力—太沙基竖向松弛土压—v—内摩擦角——上覆荷载—0P—土的容重——土的粘聚力—c三、管片内力有限元计算单元选择1、beam3单元输入参数:(1)实常数(RealConstants)AREA–横截面积IZZ–横截面惯性矩HEIGHT–梁高(2)材料参数(MaterialProperties)EX——弹模DENS——密度PRXY——泊松比其它选项采用默认即可2、Link10单元(1)RealConstantsAREA–横截面积,取单位面积1。(2)MaterialPropertiesEX——弹模,需满足E*A=地层抗力刚度其它选项采用默认即可3、link1单元(1)RealConstantsAREA–横截面积,取单位面积1。(2)MaterialPropertiesEX——弹模,为了不影响计算结果取足够小,如1。四、结点力计算1、均布荷载时中间结点:2/)(11nnxxpFn两端结点:2/)(121xxpF2/)(1NNNxxpF2、梯形荷载时中间结点:))(2())(2(611111nnnnnnnnyyppyyppFn两端结点:))(2(6112211yyppF))(2(6111NNNNNyyppF五、圆形结构内力计算1、计算参数R=4.5m;B=1.5m;H=0.5m(厚);Pw1=0.1MPa;Pe1=0.14MPa;Pw2=0.19MPa;Pe2=0.05MPa;Q1=0.14MPaQ2=0.23MPa2、在自重情况下的结构内力分布计算命令流为:/prep7et,1,beam3et,2,link10keyopt,2,3,1!*****只受压******R,1,0.75,0.015625,0.5,,,,R,2,1,,!定义混凝土衬砌材料mp,ex,1,2.76e10mp,prxy,1,0.2mp,dens,1,2600mp,ex,2,5e7mp,prxy,2,0.01!建立几何模型k,1,0,0,0circle,1,4.5!设置划分单元大小为0.5m/个Lsel,allLESIZE,all,0.5,,,,,,,1Type,1Mat,1Real,1Lmesh,all!复制建立地层弹簧CSYS,1wpcsys,-1nsel,s,loc,x,4.4,4.6NGEN,2,500,all,,,1,,,1,allsel,all*do,i,1,60type,2real,2mat,2e,i,i+500*enddonsel,s,loc,x,5.4,5.6d,all,allnsel,s,,,32D,all,,,,,,UX,,,,,allsel,allcsys,0wpcsys,-1/solACEL,0,10,0,!********施加重力命令*********solve/post1ETABLE,,SMISC,1ETABLE,,SMISC,7ETABLE,,SMISC,2ETABLE,,SMISC,8ETABLE,,SMISC,6ETABLE,,SMISC,12工况一:πg+地层反力作为主动荷载施加工况二:πg+地层反力作为被动荷载施加结果比较:工况一:作为主动荷载施加工况二:作为被动荷载施加相差弯矩/kN*m正弯矩(内弯)197.8238.817.2%负弯矩139.2195.328.7%轴力/kN195019301%最大变形量/mm4.265.6624.7%适用地层软土岩层