曲墙式衬砌计算

1、拱形曲墙式衬砌结构计算一、基本资料某一级公路隧道，结构断面如下图，围岩级别为级，围岩容重=20KN/m3，围岩的弹性抗力系数K=0.2106 kN/m，衬砌材料C20混凝土，弹形模量Eh=2.6107kPa，重度h=23 KN/m3。 衬砌结构断面（尺寸单位：cm）二、荷载确定1、根据公路隧道设计规范的有关计算公式，围岩竖向均布压力：q=0.45 2S-1 式中：S围岩级别，此处S=5；围岩容重，此处=20 kN/m3；跨度影响系数，=1+i (B-5)，毛洞跨度lm=11.81+20.1=12.01m，式中0.1为一侧平均超挖量;lm=515m时，i=0.1,此处=1+0.1（12.01-5

2、）=1.701所以，有： q=0.4525-1201.701=244.944(kPa)此处超挖回填层重忽略不计。2、围岩水平均布压力：e=0.25q=0.25244.944=61.236（kPa）三、衬砌几何要素1、衬砌几何尺寸内轮廓线半径： r=5.4039m内径r所画圆曲线的终点截面与竖直轴的夹角：=113 拱顶截面厚度d0 =0.5 m ；拱底截面厚度dn=0.5m。外轮廓线半径： R=r+d0=5.9039m拱轴线半径： r=r+0.5d0=5.6539m拱轴线各段圆弧中心角：=1132、半拱轴线长度S及分段轴长S S =pr/180=1133.145.6539/180=11.1451

3、（m）将半拱轴长度等分为8段，每段轴长为：S=S/8=11.1451/8=1.3931（m）3、各分块接缝（截面）中心几何要素（1）与竖直轴夹角i1=/8=113/8=14.1252=1+=14.125+14.125=28.253=2+=28.25+14.125=42.3754=3+=42.375+14.125=56.55=4+=56.5+14.125=70.6256=5+=70.625+14.125=84.757=6+=84.75+14.125=98.8758=7+=98.875+14.125=113（2） 接缝中心点坐标计算X1=rsin1=5.6539sin14.125=1.3798(m

4、)X2=rsin2=5.6539sin28.25=2.6761(m)X3=rsin3=5.6539sin42.375=3.8106(m)X4=rsin4=5.6539sin56.5=4.7147(m)X5=rsin5=5.6539sin70.625=5.3337(m)X6=rsin6=5.6539sin84.75=5.6302(m)X7=rsin7=5.6539sin98.875=5.5862(m)X8=rsin8=5.6539sin113=5.2044(m)y1=r(1-cos1)=5.6539(1-cos14.125)=0.1709(m)y2=r(1-cos2)=5.6539(1-cos2

5、8.25)=0.6734(m)y3=r(1-cos3)=5.6539(1-cos42.375)=1.4771(m)y4=r(1-cos4)=5.6539(1-cos56.5)=2.5333(m)y5=r(1-cos5)=5.6539(1-cos70.625)=3.7782(m)y6=r(1-cos6)=5.6539(1-cos84.75)=5.1366(m)y7=r(1-cos7)=5.6539(1-cos98.875)=6.5262(m)y8=r(1-cos8)=5.6539(1-cos113)=7.8631(m)当然也可以在下图中直接量出xi、yi衬砌结构计算图示四、计算位移1、单位位移用

6、辛普生法近似计算，按计算列表进行，单位位移的计算见表1。单位位移值计算如下：11=S/Eh1/I=1.3931768/2.6107=41.150010-612=21=S/ Ehy/I=1.39312315.4976/2.6107=124.066110-622=S/ Ehy2/I=1.393111668.7445/2.6107=625.220310-6计算精度校核：11+212+22=（41.1500+2*124.0661+625.2203）10-6=914.502510-6SS=S/ Eh(1+y)2/I=1.393117067.7397/2.6107=914.502610-6闭合差0表1、单

7、位位移计算表注：I截面惯性矩，I=bd3/12,b取单位长度不考虑轴力的影响。2、载位移主动荷载在基本结构中引起的位移（1）每一块上的作用力(竖向力Q、水平力E、自重力G),分别由下面各式求得，竖向力： Qi =qbi水平压力： Ei=ehi自重力： Gi=(di-1+di)Sh/2式中：bi衬砌外缘相邻两截面间的水平投影长度hi衬砌外缘相邻两截面间的竖直投影长度di接缝i的衬砌截面厚度均由图2直接量得，其值见表。各集中力均通过相应图形的形心。 （2）外荷载在基本结构中产生的内力锲块上各集中力对下一接缝的力臂由图直接量得，分别记以aq 、ae、ag。内力按下式计算之：表2、载位移Mop计算表弯

8、矩： 轴力： 式中 xi、yi相邻两接缝中心点的坐标增值，按下式计算： xi=xi- xi-1 yi=yi- yi-1Moip和Noip的计算见表表3、载位移Noip计算表基本结构中，主动荷载产生的弯矩的校核为：M08q=-qB（x8-B/4）/2=-244.94411.8078(5.2044-11.8078/4)/2=-3257.3240M08e=-eH2/2=-61.2368.17682/2=-2047.1213M08g=-Gi(x8-xi+agi)=-16.0207(75.2044-1.3798-2.6761-3.8106-4.7147-5.3337-5.6302-5.5862+0.68

9、99+0.6482+0.5673+0.4521+0.3095+0.1483-0.022-0.1909)=-158.6354M08p=M08q+M08e+M08g=-3257.3240-2047.1213-158.6354=-5463.0807另一方面，从表中得到M08p=-5432.7146闭合差 =|5463.0807-5432.7146|/5432.7146=0.56（3）主动荷载位移计算过程见表表4、主动荷载位移计算表1p=S/EhMp0/I=-2257967.13381.3931/2.6107=-120983.61610-62p=S/ EhMp0y/I=-10457532.51341.

10、3931/2.6107=-560322.63610-6计算精度校核1p+2p=-（120983.616+560322.636）10-6=-681306.25210-6Sp=S/ EhMp0(1+y)/I=-12715499.64711.3931/2.610-6=-681306.25210-6闭合差=0。3、载位移单位弹性抗力及相应的摩擦力引起的位移（1）各接缝处的抗力强度抗力上零点假定在接缝3，3=42.375=b;最大抗力值假定在接缝5, 5=70.625=h。最大抗力值以上各截面抗力强度按下式计算：i=h（coS2b-coS2i）/（coS2b-coS2h）计算得，3=0，4=0.5534

11、h ，5=h 。最大抗力值以下各截面抗力强度按下式计算：=h1-(yi/ yh)2式中：yi所考察截面外缘点到h点的垂直距离；yh墙角外缘点到h点的垂直距离。由图量得：y6=1.4184m; y7=2.8695m; y8=4.2316m;则有： 6=h 1-（1.4184/4.2316）2= 0.8876h 7=h 1-（2.8695/4.2316）2= 0.5402h 8=0;按比例将所求得的抗力绘在图上。（2）各楔块上抗力集中力按下式近似计算：式中：楔块i外缘长度，可通过量取夹角，用弧长公式求得，的方向垂直于衬砌外缘，并通过楔块上抗力图形的形心。（3）抗力集中力与摩擦力之合力按下式计算:

12、式中 ：围岩与衬砌间的摩擦系数，此处取=0.2,则 =1.0198 其作用方向与抗力集中力的夹角为=arctg=11.3099。由于摩擦阻力的方向与衬砌位移方向相反，其方向朝上。画图时，也可取切向：径向=1：5的比例求出合力Ri的方向。的作用点即为与衬砌外缘的交点。 将的方向线延长，使之交于竖直轴。量取夹角k。将分解为水平与竖向两个分力：RH= Ri sink RV= Ri cosk以上计算列入表5。表5、弹性抗力及摩擦力计算表（4）计算单位抗力图及其相应的摩擦力在基本结构中产生的内力 弯矩 ： 轴力 ： 式中：rji -力Ri至接缝中心点ki的力臂，由图量得。计算见表6和表7。表6、M0计算

13、表表7、N0计算表（5）单位抗力及相应摩擦力产生的载位移计算过程见表8。表8、单位抗力及相应摩擦力产生的载位移计算表1=S/EhM01/I=1.3931(-2261.0895)/2.6107= -121.150910-62=S/ EhM0y/I=1.3931(-14274.0927)/2.6107= - 764.816910-6校核为： 1+2= -(121.1509+ 764.8169) 10-6=-885.967810-6 S=S/ EhM0(1+y)/I=1.3931(-16535.1822)/2.6107=-885.967810-6闭合差=0。4、墙底（弹性地基上的刚性梁）位移单位弯矩

14、作用下的转角：a=1/(KI8)= 96 /0.2106=48010-6主动荷载作用下的转角：ap=aM8p0=-5432.714648010-6 = -2607703.0110-6单位抗力及相应摩擦力作用下的转角：a=aM80=48010-6(-14.8845)=-7144.5610-6五、解力法方程衬砌矢高: f=y8=7.8631m计算力法方程的系数为：a11=11+a=（41.15+480）10-6=521.1510-6a12=12+fa=（124.0661+7.8631480）10-6=3898.354110-6a22=22+f2a=（625.2203+7.86317.8631480

15、）10-6=30302.824310-6a10=1p+ap+(1+a)h =-(120983.616+2607703.01+121.1509h+7144.56h) 10-6= -(2728686.63+7265.7109h)10-6a20=2p+fap+(2+fa)h = -(560322.636+7.86312607703.01+764.8169h+7.86317144.56h) = - (21064952.1+56943.2066h) 10-6以上将单位抗力图及相应摩擦力产生的位移乘以h，即为被动荷载的载位移。求解方程：X1 =(a12a20 - a22a10)/( a11a22 - a1

16、22)=(3898.3541(-21064952.1-56943.2066h)+30302.8243(2728686.63+7265.7109h)/(521.1530302.8243-3898.35413898.3541)=（954.8299-3.0467h）其中： X1p=954.8299, X1= -3.0467X2 = (a12a10 - a11a20)/( a11a22 - a122)=(3898.3541(-2728686.63-7265.7109h)+521.15(21064952.1+56943.2066h)/(521.1530302.8243-3898.35413898.354

17、1)=（572.3126+2.2711h）其中： X2p=572.3126, X2= 2.2711六、计算主动荷载和被动荷载（h =1）分别产生的衬砌内力计算公式为： 计算过程列入表9和表10中。表9、主、被动荷载作用下衬砌弯矩计算表表10、主、被动荷载作用下衬砌轴力计算表 七、最大抗力值的求解首先求出最大抗力方向内的位移。由式：计算过程列入表11。表11、最大抗力位移修正计算表位移值为： hp=1.3931466608.6752/2.6107=25001.251710-6 h=1.3931(-1615.2336)/2.6107= -86.545510-6最大抗力值为：h =hp /(1/K-

18、h)= 25001.251710-6/1/(0.2106)+86.545510-6=273.1019八、计算衬砌总内力按下式计算衬砌总内力：M=Mp+h MN=Np+h N计算过程列入表12表12、衬砌总内力计算表计算精度校核：根据拱顶切开点之相对转角和相对水平位移应为零的条件来检查。式中：=-1.393119830.6721/2.6107=-1062.542710-6a=M8a=2.23148010-6=1070.8810-6闭合差：=(1070.88-1062.5427)/1070.88100= 0.78%式中：=1.3931155928.2026/2.6107=8354.753010-6fa =7.86311070.8810-6=