隧道工程课程设计

武汉理工大学课程设计1初始条件某高速公路隧道通过=3\*ROMANIII类围岩(即=4\*ROMANIV级围岩)，埋深H=30m，隧道围岩天然容重γ=23KN/m3，计算摩擦角ф=35o,变形模量E=6GPa,采用矿山法施工；衬砌材料采用C25喷射混凝土，材料容重，变形模量。2隧道洞身设计2.1隧道建筑界限及内轮廓图的确定该隧道横断面是根据两车道高速公路=4\*ROMANIV级围岩来设计的，根据《公路隧道设计规范》确定隧道的建筑限界如下：W—行车道宽度；取3.75×2mC—余宽；因设置检修道，故余宽取为0mJ—检修道宽度；双侧设置，取为1.0×2mH—建筑限界高度；取为5.0m2—左侧向宽度；取为1.0m—右侧向宽度；取为1.5m—建筑限界左顶角宽度；取1.0m—建筑限界右顶角宽度；取1.0mh—检修道高度；取为0.25m隧道净宽为1.0+1.0+7.50+1.50+1.0=12m设计行车速度为120km/h,建筑限界左右顶角高度均取1m；隧道轮廓线如下图：图1隧道内轮廓限界图根据规范要求，隧道衬砌结构厚度为50cm（一次衬砌为15cm和二次衬砌35cm）通过作图得到隧道的尺寸如下：图2隧道内轮廓图得到如下尺寸:3隧道衬砌结构设计3.1支护方法及衬砌材料根据《公路隧道设计规范》（JTG-2004），本设计为高速公路，采用复合式衬砌，复合式衬砌是由初期支护和二次衬砌及中间防水层组合而成的衬砌形式。复合式衬砌应符合下列规定：1初期支护宜采用锚喷支护，即由喷射混凝土，锚杆，钢筋网和钢筋支架等支护形式单独或组合使用，锚杆宜采用全长粘结锚杆。2二次衬砌宜采用模筑混凝土或模筑钢筋混凝土结构，衬砌截面宜采用连结圆顺的等厚衬砌断面，仰拱厚度宜与拱墙厚度相同。由规范8.4.2-1，对于两车道=4\*ROMANIV级围岩：初期支护：拱部边墙的喷射混凝土厚度为12-15cm，拱墙的锚杆长度为2.5-3m，锚杆间距为1.0-1.2m；二次衬砌厚度：拱墙混凝土厚度为35cm因此确定衬砌尺寸及规格如下：深埋隧道外层初期支护，根据规范规定，采用锚喷支护，锚杆采用普通水泥砂浆锚杆，规格HRB20×2.5m，采用梅花型局部布设，采用C25喷射混凝土。初次衬砌：采用C25喷射混凝土，厚度为14.8cm。防水板：采用塑料防水板及无纱布，且无纺布密度为300g/m2,防水板应采用铺满的EVA板防水层，厚度为2.0mm，搭接长度为150mm。二次衬砌：根据《公路隧道设计规范》，采用C25号模筑防水混凝土，厚度为35cm。整个衬砌厚度为14.8+0.2+35=50cm。3.2隧道深浅埋的确定及围岩压力的计算隧道的内轮廓尺寸为B=12.37m，H=13.69m因为=4\*ROMANIV级围岩需要预留变形量，查《公路隧道设计规范》8.4.1知=4\*ROMANIV级围岩需预留变形量为50-80mm，衬砌厚度为50cm，又每侧超挖量为10cm，故隧道的近似开挖尺寸为：由于是=4\*ROMANIV级围岩—深浅埋隧道的分界深度—等效荷载高度—跨度影响系数；i—围岩压力增减率，当取i=0.1埋深故为深埋隧道。又可用公式计算均布垂直压力：因为该隧道围岩级别为=4\*ROMANIV围岩水平均不压力为：4衬砌结构内力计算4.1基本资料公里等级高速公路围岩级别=4\*ROMANIV级围岩容重r=23KN/m3弹性抗力系数K=0.18×106KN/m变形模量材料容重材料变形模量衬砌厚度d=0.5m图3衬砌结构断面图4.2荷载的确定4.2.1围岩压力的确定隧道的内轮廓尺寸为B=12.37m，H=13.69m因为=4\*ROMANIV级围岩需要预留变形量，查《公路隧道设计规范》8.4.1知=4\*ROMANIV级围岩需预留变形量为50-80mm，衬砌厚度为50cm，又每侧超挖量为10cm，故隧道的近似开挖尺寸为：由于是=4\*ROMANIV级围岩—深浅埋隧道的分界深度—等效荷载高度—跨度影响系数；i—围岩压力增减率，当取i=0.1埋深故为深埋隧道。又可用公式计算均布垂直压力：考虑一衬后的围岩释放变形折减系数0.44.2.2衬砌自重（1）全部垂直荷载q=61.768+11=72.768KN/m2围岩水平均布压力e=0.4×72.768=29.1KN/m24.3衬砌几何要素4.3.1衬砌的几何尺寸内轮廓线半径：内径所画圆曲线端点截面与竖直轴线的夹角：拱顶截面厚度，拱底截面厚度4.3.2半拱轴线长度S及分段轴长将半拱轴长度等分为8段，则4.4计算位移4.4.1单位位移用辛普生法近似计算，按计算列表进行，单位位移的计算列表见表4-1表4-1单位位移计算表截

面αsinαcosαxyd1/Iy/Iy2/I(1+y)2/I积

分

系

数1/300.0000.0001.0000.0000.0000.50096.0000.0000.00096.0001115.2500.2630.9651.5750.2210.50096.00021.2164.689143.1214230.5000.5080.8622.9920.8420.50096.00080.83268.061325.7252345.7500.7160.6984.2021.7860.50096.000171.456306.220745.1324461.0000.8750.4855.2343.2170.50096.000308.832993.5131707.17772576.2500.9710.2386.1685.6100.50096.000538.5603021.32224194.44224691.5001.000-0.0266.4298.1760.50096.000784.8966417.31008083.102227106.750.958-0.2886.00710.7220.50096.0001029.312211036.288313190.900748122.000.848-0.5304.97213.0200.50096.0001249.920016273.955818869.79981∑768.003547.072229568.599337430.7337注：1I—-截面的惯性矩，I=bd3/12,b取单位长度2不考虑轴力影响单位位移值计算如下：计算精度校核：闭合差△=04.4.2载位移——主动荷载在基本结构中引起的位移1）每一块上的作用力(竖向力Q、水平力E、自重力G),分别由下面各式求得，Qi=q×biEi=e×hiGi=(di-1+di)/2×△S×rh其中：bi——衬砌外缘相邻两截面间的水平投影长度hi——衬砌外缘相邻两截面间的竖直投影长度di——接缝i的衬砌截面厚度均由图3直接量得，其值见表4-2。各集中力均通过相应图形的形心表4-2载位移Mop计算表截面投影长度集中力S-Qaq-Gag-GagbhQGEaqagae00.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.00011.5750.224114.61023.1007.4010.7880.7880.1120.112-18.191-0.82921.4240.612103.62223.10020.2200.7120.7120.306-73.779-16.447-6.18731.2100.95088.04923.10031.3880.6050.6050.475-53.270-13.976-14.90941.0200.97974.22323.10032.3460.5100.5100.490-37.854-11.781-15.83350.9332.84567.89323.10093.9990.4670.4671.423-31.672-10.776-133.713360.2672.56619.42923.10084.7810.1340.1341.283-2.594-3.084-108.774470.4172.54630.34423.10084.1200.2090.2091.273-6.327-4.816-107.085580.7852.29857.12323.10075.9260.3930.3931.149-22.421-9.067-87.239续表4-2∑i-1（Q+G）∑i-1ExyΔxΔy-Δx∑i-1（Q+G）-Δy∑i-1EMoip0.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.0001.5310.1941.5310.1940.0000.000-18.908137.7107.4012.9640.7651.4330.571-197.3388-4.226-316.8855264.43127.6214.2101.6761.2460.911-329.4811-25.163-753.6844375.58159.0095.1882.8690.9781.193-367.3188-70.398-1256.8669472.90491.3565.8324.2520.6441.383-304.5500-126.3455-1863.9225563.896185.3546.1525.7640.3201.512-180.4477-280.2566-2439.0779606.425270.1356.1437.309-0.0091.5455.458-417.3599-2969.2007659.870354.2555.8068.817-0.3371.508222.376-534.2166-3399.77742）外荷载在基本结构中产生的内力块上各集中力对下一接缝的力臂由图直接量得，分别记以aq、ae、ag。内力按下式计算之：弯矩：轴力：式中Δxi、Δyi——相邻两接缝中心点的坐标增值。Δxi=xi-xi-1Δyi=yi-yi-1Moip和Noip的计算见表4-3及表4-4。表4-3载位移Noip计算表截

面αsinαcosα∑i（Q+G）∑iEsinα*∑ii（Q+G）cosα*∑iEENop00.0000.0001.0000.0000.0000.0000.0000.000115.2500.2630.965137.7107.40136.2227.14029.082230.5000.5080.862264.43127.621134.20923.799110.410345.7500.7160.698375.58159.009269.02941.176227.853461.0000.8750.485472.90491.356413.61144.290369.321576.2500.9710.238563.896185.354547.73644.056503.680691.5001.000-0.026606.425270.135606.218-7.071613.2897106.7500.958-0.288659.870354.255631.872-102.0955733.9678122.0000.848-0.530740.093430.181627.634-227.9611855.595主动荷载位移计算过程见表4-4表4-4主动荷载位移计算表截面Mp01/Iy/I1+yMp0/IyMp0/IMp0(1+y))/I积分系数1/300.00096.0000.0001.0000.0000.0000.00011-18.90896.00018.6241.194-1815.1883-352.1466-2167.322942-316.885596.00073.4401.765-30420.9982-23272.0051-53693.003423-753.684496.000160.8962.676-72353.6692-121264..788-193618..48044-1256.866996.000275.4243.869-120659..391-346171..792-466831..18325-1863.922596.000408.1925.252-178936..786-760839..215-939776..00246-2439.077996.000553.3446.764-234151..560-13496499.594-15838011.15427-2969.200796.000701.6648.309-285043..895-20833855.832-23684299.72848-3399.777496.000846.4329.817-326378..291-28776777.394-32040555.6861∑-10831466.795-60597444860△1p=△S/Eh×∑Mp0/I=8.4×10-5×（-1083146.795）=-9098.426×10-5△2p=△S/Eh×∑Mp0y/I=8.4×10-8×(-6059744.064)=-50901.850×10-5计算精度校核△Sp=△1p+△2p△Sp=△S/Eh×∑Mp0(1+y)/I因此，△Sp=8.4×10-5×（-7142890.860）=-60000.276×10-5△1p+△2p=-（9098.426+50901.850）×10-6=-60000.276×10-5闭合差△=0.000。4.4.3载位移——单位弹性抗力及相应的摩擦力引起的位移1）各接缝处的抗力强度按假定拱部弹性抗力的上零点位于与垂直轴接近450的第3截面，α3=45.75°=αb;最大抗力位于第5截面,α5=76.25°=αh;拱部各截面抗力强度，按镰刀形分布，最大抗力值以上各截面抗力强度按下式计算：σi=σh（cos2αb-cos2αi）/（cos2αb-cos2αh）计算得，σ3=0，σ4=0.5847σh，σ5=σh。边墙截面弹性抗力计算公式为：σ=σh[1-(yiˊ/yhˊ)2]式中yiˊ——所求抗力截面与外轮廓线交点到最大截面抗力截面的垂直距离；yhˊ——墙底外边缘cˊ到最大抗力截面的垂直距离。(yiˊ和yhˊ在图3中可量得)y6ˊ=2.47m;y7ˊ=4.93m;y8ˊ=6.87m;则有：σ6=σh[1-（2.47/6.87）2]=0.8707σhσ7=σh[1-（4.93/6.87）2]=0.4850σhσ8=0;按比例将所求得的抗力绘在图4上。图4结构抗力图2）各楔块上抗力集中力按下式近似计算：式中，——楔块i外缘长度，由图3量得。的方向垂直于衬砌外缘，并通过楔块上抗力图形的形心。3）抗力集中力与摩擦力之合力按近似计算:式中μ——围岩与衬砌间的摩擦系数。取μ=0.2,则=1.0198其作用方向与抗力集中力的夹角为β=arctanμ=11.301°。由于摩擦阻力的方向与衬砌位移方向相反，其方向朝上。的作用点即为与衬砌外缘的交点。将的方向线延长，使之交于竖直轴。量取夹角ψk(自竖直轴反时针方向量度)。将分解为水平与竖向两个分力：RH=RisinψkRV=Ricosψk以上计算例入表4-5中，并参见图3。表4-5弹性抗力及摩擦力计算表截面σ(σh)(σi-1+σi)//2△S外(σh)R(σh)ψksinψkcosψk30.00000.0000.00000.00000.0000.0001.00040.58470.2922.16700.646164.6760.9040.42851.00000.7922.16701.751079.9260.9850.17560.87070.9352.16702.067095.1760.996-0.09070.48500.6782.16701.4980110.4260.937-0.34980.00000.2432.16700.5359125.6760.812-0.583续表4-5RH(σh)RV(σh)vhRi(σh)0.0000.0000.0000.5840.2760.2760.5840.6341.7240.3060.5832.3081.7172.059-0.1860.3964.3672.0271.404-0.523-0.1275.7701.4690.435-0.313-0.4396.2060.5254）计算单位抗力图及其相应的摩擦力在基本结构中产生的内力弯矩轴力式中rKi-----力Ri至接缝中心心点K的力臂，由由图3量得，计算算见表4-6和表4-7表4-6Mσ0计算算表截面号R4=0.6461σσhR5=1.7510σσhR6=2.0670σσhR7=1.4980σσhR8=0.5359σσhMoσ(σh)r4i-R4r4ir5i-R5r5ir6i-R6r6ir7i-R7r7ir8i-R8r8i40.8236-0.532-0.53252.5623-1.6560.7563-1.324-2.98064.0536-2.6192.4634-4.3130.8280-1.711-8.64475.3654-3.4673.6842-6.4512.3729-4.9050.7853-1.176-15.99987.2631-4.6935.3648-9.3943.8727-8.0052.4563-3.6800.9365-0.502-26.273表4-7Nσ0计算算表截

面

号αsinαcosαΣRV(σh)ΣRH(σh)sinαΣRV((σh)cosαΣRH((σh)Noσ(σh)461.00000.87440.48530.27630.58400.24160.2834-0.0418576.00000.97010.24260.58262.30800.56520.55990.0054691.50000.9997-0.02540.39624.36660.3960-0.11080.50687107.000000.9566-0.2915-0.12665.7704-0.1211-1.68191.56078122.000000.8486-0.5290-0.43926.2057-0.3727-3.28292.91025）单位抗力及相应摩摩擦力产生的的载位移计算过程见表4--8。表4-8单位抗力及及相应摩擦力力产生的载位位移计算表截面号Mσ0（σh）1/Iy/I(1+y)Mσ01/I（σh）Mσ0y/I（σh）Mσ0(1+y)/I（σh）积分系数

1/34-0.53296.000275.4243.869-51.084-146.5611-197.645525-2.98096.000408.1925.252-286.0599-1216.3224-1502.388346-8.64496.000553.3446.764-829.8088-4783.0116-5612.822427-15.99996.000701.6648.309-1535.8889-11225.8813-12761.770348-26.27396.000846.4329.817-2522.1888-22238.1134-24760.33221Σ-3857.2556-27288.6612-31145.8868△1σ=△S/Eh×∑Mσσ01/I=8.4×100-8×(-38557.2566)=-322.40100×10-5△2σ=△S/Eh×∑MMσ0y/I=8.4×100-5×(-272288.6112)=--229..2243××10-5校核为：△1σ+△2σ=-(32.44010+229.22243)×10-11=-261..6253××10-6△Sσ=△S/Eh×∑MMσ0(1+y)/I=8.4×100-11×(-311145.8668)=-2611.62533×10-6闭合差△=0。4.4.4墙底（弹弹性地基上的的刚性梁）位位移1）单位弯矩作用下下的转角：β1=1/(KI8))=-966.000/0.188×106=-533..33×10-62）主动荷载作用下下的转角：βp=β1M8p0=-33999.7744×533.333×10-6=-188132011.467××10-63）单位抗力及相应应摩擦力作用用下的转角：：βσ=β1M8σ0=533.33×110-6×(-26..273)=-114012..179×10-64.5解力法方程衬砌矢高f==y8=11.7776m计算力法方程的系系数：a11=δ11+ββ1=（64.5122+533.333）×10-6=597.8842×100-6a12=δ12+ffβ1=（297.9554+11.7776×533.333）×10-6=6578..448×110-6a22=δ22+ff2β1=（2483..762+88.817××11.77762×533.333）×10-6=764422.860××10-6a10=△1p+ββp+(△1σ+βσ)×σh=-(909844.26+18132201.4667+324.0010σh+140122.179σσh)×100-6=-(19044185.7727+143366.189σσh)×10-66a20=△2p+ffβp+(△2σ+fβσ)×σh=-(5090018.500+11.7776×18113201..467+2292..243σh+11.7776×140122.179σσh)=-(2188612788.98+1672999.1633σh)×100-6以上将单位抗力图图及相应摩擦擦力产生的位位移乘以σh倍，即被动动荷载的载位位移。求解方程：X1=(a12a220-a22a10)/(a11a22-a122)=（594.53500-1.7788σh）其中：X1pp=594.55350,XX1σ=-1.7778X2=(a12a110-a11a20)/(a11a22-a122)=（243.8255++2.3411σh）其中：X2pp=243.88255,X2σ=2.34414.6计算主动荷载和被被动荷载（σh=1）分别产生生的衬砌内力力计算公式为：和计算过程列入表44-9和表4-10中。表4-9主、被动荷荷载作用下衬衬砌弯矩计算算表截面MopX1pyX2p*y[Mp]Moσ(σh)X1σ(σh)X2σ*y(σh)[Moσ](σhh)00.000594.5350.0000.000594.5350.000-1.7780.000-1.7781-18.809594.5350.19447.302623.0280.000-1.7780.454-1.3242-315.5455594.5350.765186.527465.5170.000-1.7781.7910.0133-747.5655594.5351.676408.652255.6210.000-1.7783.9242.1464-1240.4666594.5352.869699.53553.604-0.532-1.7786.7164.4065-1816.5111594.5354.2521036.7466-185.2300-2.980-1.7789.9545.1966-2345.2773594.5355.7641405.4100-345.3288-8.644-1.77813.4943.0727-2812.8662594.5357.3091782.1211-436.2077-15.999-1.77817.110-0.6668-3169.3221594.5358.8172149.8099-424.9766-26.273-1.77820.641-7.410表4-10主、被动动荷载作用下下衬砌轴力计计算表截面NopcosaX2pcosφ［Np］Noσ(σh)X2σcosφ(σh))［Nσ］(σh)00.0001.000243.826243.8260.0002.3412.341138.5000.946230.644269.1440.0002.2142.2142128.4140.827201.721330.1350.0001.9371.9373251.4740.650158.591410.0650.0001.5231.5234393.0430.428104.293497.336-0.0421.0010.9605526.6470.17542.650569.2970.0050.4090.4156625.417-0.090-21.997603.4200.507-0.2110.2967727.270-0.349-85.094642.1761.561-0.8170.7448822.162-0.583-142.1999679.9632.910-1.3651.5454.7计算最大抗力值首先求出最大抗力力方向内的位位移。由式：并考虑接缝5的径径向位移与水水平方向有一一定的偏离，因因此将其修正正如下计算过程列入表33-11。表4-11最大抗力位位移修正计算算表截面号Mp/IMσ/Iyiy5-yiMp/I(Y5-Yii)Mσ/I

(Y5-Yii)积分系数

1/33078292.6449-234.14000.0005.616-439691.5519-2314.9294182044.7991-174.33440.2185.398442877.7783-1941.0532261302.61191.6940.8424.774292658.77048.0884333662.0779282.5371.7863.830128925.77641082.1177247058.9666580.2403.2042.41217026.22251399.540045-24392.3384684.2665.6160.0000.0000.0001Σ207860.2245-1782.3559位移值为：δhp=8.4×110-8×2078660.2455×0.97=1693..645×110-5δhσ=8.4×10-88×(-17882.3599)×0..97=-14.9772×10-5则可得最大抗力σh=δhpp/(1//K-δhσ)=16993.6455×10-5/[1//(0.188×106)+14.9972×10-5]=109.0743.8计算衬砌总总内力按下式进行计算：：M=Mp+σhMσN=Np+σhNσ计算过程列入表44-12表4-12衬砌总内力力计算表计算精度校核：截

面

号MpMσ[M]M/IMy/INpNσ[N]e积分系数

1/30329.456-300.236629.2202805.12000.000243.826202.340446.1650.065511286.546-274.455512.0911160.69771351.0511269.144191.401460.5440.026342146.320-158.7588-12.438-1194.0009-5480.5000330.135167.399497.534-0.025023-24.00025.8321.832175.8821768.6655410.065131.608541.6730.003444-243.2766247.602