隧道工程课程设计报告(完整)

1、隧道工程课程设计一、 工程概况3某地区一暗挖双线马蹄形隧道， 埋深 h=125m，围岩等级为 v 级，地层平均容重 16.0 kN/m3。 宽度 B=13.08m，隧道采用复合式衬砌形式，衬砌厚度为0.42m ，配筋采用 22200m，m钢材采用 HRB335，钢筋保护层厚度 50mm。二、 计算1、 衬砌结构的计算模型隧道工程建筑物是埋置于地层中的结构物， 它的受力和变形与围岩密切相关， 支护结构 与围岩作为一个统一的受力体系相互约束，共同工作。这种共同作用正是地下结构与地面 结构的主要区别。根据本工程浅埋及松散地层的特点，使用阶段结构安全性检算采用“荷 载结构”模式，即将支护和围岩分开考虑

2、，支护结构是承载主体，围岩作为荷载的来源 和支护结构的弹性支承。支护结构与围岩的相互作用是通过弹性支承对支护结构施加约束 来实现的。计算模型中， 二衬结构采用弹性平面梁单元模拟， 弹性抗力以及隧底地基均采用弹簧单 元模拟。组合荷载根据不同作用方向分别转换成等效节点力施加在相应的单元结点上。具 体计算模型见图 1。图 1 计算模型2、荷载计算围岩压力计算参照课本中有关我国铁路隧道推荐的方法进行确定 路隧道设计规范 ，深浅埋分别计算。按破坏阶段设计计算垂直压力公式：q=r x h q = 0.45 x 2(s-1) x r x w双线隧道) 或参照 铁式中： h q 等效荷载高度值S 围岩级别r

3、围岩的容重w宽度影响系数，其值为 w=1+i （B-5）计算得， q=0.45x2(5-1)x16000x1.805=2.082816e6N/m水平均布松动压力系数取0.3 ，则 e=0.3q=0.0634e6N/m3、ANSYS操作命令流! 荷载结构方法计算 finish!/clear!/COM,Structural ! /prep7 ! *AFUN,deg !（ 马蹄形断面）退出当前处理程序清除以前数据 , 重新开始一个新的分析定义分析类型 , 结构分析 （ 热分析、流体分析等 ） 进入前处理器 定义角度单位为度（缺省为弧度， RAD）定义建模及材料参数的一些变量值*set,Py,2.08

4、2816e5*set,px1,0.0634e6*set,px2,0.0634e6定义垂直围岩压力大小（若有地表荷载加地表荷载值）*set,cylxsh,0.3 !*set,cyl,Py*cylxsh*set,CQHD,0.47 *set,CQDYCD,0.17 相等)*set,CQETXML,31e9*set,CQUBSB,0.2*set,WYTXKL,100e6*set,WYMD,1600*set,CQMD,2500定义侧压力系数水平侧压力定义初支或二衬厚度 定义梁单元长度参数及弹簧单元面积梁单元长度与弹簧单元面积定义衬砌（初支或二衬）的弹性模量 定义衬砌的泊松比 定义围岩的弹性抗力系数 定

5、义围岩的密度 定义衬砌的密度! 定义单元类型及材料属性及单元实常数et,1,beam3 !定义 1 号单元为梁单元mp,ex,1,CQETXML!定义 1 号材料的弹性模量mp,prxy,1,CQUBSB!定义 1 号材料的泊松比mp,dens,1,CQMD!定义 1 号材料的密度R,1,CQHD,1/12*CQHD*CQHD*CQHD,CQHD !1-实常数号；第一个参数为梁截面的面积；第二 个参数为梁单元的转动惯量；第三个参数为梁高。! 建立几何模型 ! 创建关键点k,1,0,-4.65025, !k- 点命令； 1- 点号；坐标（ x,y,z） k,2,2.36643,-4.45519,

6、k,3,4.66945,-3.87810,k,4,5.73443,-3.21404,k,5,6.39434,-2.14650,k,6,6.65675,-1.08747,k,7,6.745,0,k,8,4.76943,4.76943,k,9,0,6.745, k,10,-4.76943,4.76943, k,11,-6.745,0, k,12,-6.65675,-1.08747, k,13,-6.39434,-2.14650, k,14,-5.73443,-3.21404, k,15,-4.66945,-3.87810, k,16,-2.36643,-4.45519, ! 二衬轴线 larc,1

7、,3,2 ! 三点画圆弧： 起点 - 终点中间点 最好全按逆时针画， 弯矩及轴力的 方向一致 larc,3,5,4 larc,5,7,6 larc,7,9,8 larc,9,11,10 larc,11,13,12 larc,13,15,14 larc,15,1,16 ! 设置线单元材料属性，划分单元 （ 初支单元 ） lsel,s,1,8 !LSEL,Type,Item,Comp,VMIN, VMAX,VINC lATT,1,1,1 ! 给线单元付材料号、 实常数、 单元类型号 LATT, MAT, REAL, TYPE, ESYS MSHKEY,1 ! 设置单元为映射单元（单元比较规则） l

8、esize,1,CQDYCD !LESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE lesize,2,CQDYCD lesize,3,CQDYCD lesize,4,CQDYCD lesize,5,CQDYCD lesize,6,CQDYCD lesize,7,CQDYCD lesize,8,CQDYCD lsel,all lmesh,all ! 划分线单元nummrg,all,low ! 所有号数从 1 开始 numcmp,node,eord ! 压缩节点号并排序! 施加水平弹簧支撑单元（在梁单元的每个节点上分别产生水平方向弹簧单元） *do,i,2,118 !do

9、 循环，要对应下面一个 *EDNDO a1=abs(ny(i-1)-ny(i+1)*WYTXKL ! 计算节点 i-1 和节点 i+1 之间 y 的变化量 * 弹性抗力 kx=a1/2PSPRNG,i,TRAN,kx,0.2, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM! 弹簧单元 , 节点号 , 平移 , 弹性系数 , 弹簧水平长度 产生 水平弹簧*enddo*do,i,121,235 !do 循环，要对应下面一个 *EDNDO a1=abs(ny(i-1)-ny(i+1)*WYTXKL ! 计算节点 i-1 和节点 i+1 之间 y 的变化量 * 弹性

10、抗力 kx=a1/2PSPRNG,i,TRAN,kx,-0.2, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM! 弹簧单元 , 节点号 , 平移 , 弹性系数 , 弹簧水平长度 产生 水平弹簧*enddo a1=abs(ny(235)-ny(1)*WYTXKL ! 计算节点 294 上的弹性抗力 kx=a1/2PSPRNG,236,TRAN,kx,-0.2, ! 施加竖直弹簧支撑单元（在梁单元的每个节点上分别产生竖直方向弹簧单元）计算节点 1 上的弹性抗力b1=abs(nx(236)-nx(2)*WYTXKL ! ky=b1/2PSPRNG,1,TRAN,k

11、y,-0.2,*do,i,2,57 !do循环，要对应下面一个 *EDNDO.资料.b1=abs(nx(i-1)-nx(i+1)*WYTXKL !计算节点 i-1 和节点 i+1 之间 x 的变化量 * 弹性抗力 ky=b1/2PSPRNG,i,TRAN,ky,-0.2, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM! 弹簧单元 , 节点号 , 平移 , 弹性系数 , 弹簧水平长度产生水平弹簧*enddo *do,i,58,180 !do 循环，要对应下面一个 *EDNDO b1=abs(nx(i-1)-nx(i+1)*WYTXKL ! 计算节点 i-1 和

12、节点 i+1 之间 x 的变化量 * 弹性抗力 ky=b1/2PSPRNG,i,TRAN,ky,0.2, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM! 弹簧单元 , 节点号 , 平移 , 弹性系数 , 弹簧水平长度 产生 水平弹簧*enddo*do,i,181,235 !do 循环，要对应下面一个 *EDNDO b1=abs(nx(i-1)-nx(i+1)*WYTXKL ! 计算节点 i-1 和节点 i+1 之间 x 的变化量 * 弹性抗力 ky=b1/2PSPRNG,i,TRAN,ky,-0.2, !PSPRNG, NLOC, TYPE, K, DX,

13、DY, DZ, ELEM! 弹簧单元 , 节点号 , 平移 , 弹性系数 , 弹簧水平长度 产生 水平弹簧*enddob1=abs(nx(235)-nx(1)*WYTXKL ! 计算节点 1 上的弹性抗力 ky=b1/2PSPRNG,235,TRAN,ky,-0.2,! 施加边界条件与围岩压力 /solu ! 进入求解器计算!nsel,s,1 !选择约束的节点 , 仰拱中间节点!d,all,ux !施加水平方向的约束allselfcum,add, ! 一定要，使荷载能叠加! 竖向荷载（ y 方向压力，最大跨度以上节点施加竖向荷载，换算成节点荷载，在节点上施 加）*do,i,58,180 !在

14、59 181 号节点施加竖向围岩压力a=-Py*0.5*abs(nx(i-1)-nx(i+1) !计算作用在节点上的节点力f,i,fy,a*enddo! 水平向荷载( X 方向压力，全部节点施加水平荷载，换算成节点荷载，在节点上施加)*do,i,2,118 !在 2 118 号节点施加水平围岩压力a=-(px1+(px2-px1)/(ny(148)-ny(1)*(ny(148)-ny(i)*abs(nY(i-1)-nY(i+1)*0.5 f,i,fX,a*enddo*do,i,119,235 !在 1120 235 号节点施加水平围岩压力a=(px1+(px2-px1)/(ny(119)-ny

15、(1)*(ny(119)-ny(i)*abs(nY(i-1)-nY(i+1)*0.5f,i,fX,a*enddo a=(px1+(px2-px1)/(ny(119)-ny(1)*(ny(119)-ny(236)*abs(nY(235)-nY(1)*0.5 f,236,fX,a重力加速度以向下位置，不用加负号! 设定重力加速度 ACEL,0,9.8,0, ALLSEL,ALL finish! 删除轴力为拉力的弹簧单元/prep7*do,i,239,261EDELE, iNDELE, i*enddo*do,i,297,415EDELE, iNDELE, i*enddo*do,i,451,471ED

16、ELE, iNDELE, i*enddo*do,i,529,651EDELE, iNDELE, i*enddo! 进入求解器计算/solutime,1 !定义时间步OUTRES,ALL,LAST, !设置输出选项solcontrol,0,0 !设置是否为非线性求解nlgeom,1 !大变形开关PIVCHECK,1 !如果不收敛，计算退出neqit,100, !非线性求解，平衡力最大迭代次数autots,1 !自动步长nsubst,5 !子步数kbc,0 !指定荷载按线性增加ESEL,ALLNSEL,ALLsolve !求解save,jg,db !将第一时间步的计算结果保存! 进入后处理器，对结

17、果进行后处理/POST1ESEL,s,TYPE,1节点的轴力 节点的轴力 节点的弯矩 节点的弯矩 节点的剪力 节点的剪力NSLE,s,ALLetabel,ni,smisc,1 !定义梁etabel,nj,smisc,7 !定义梁etabel,mi,smisc,6 !定义梁etabel,mj,smisc,12 !定义梁etable,Qi,smisc,2 !定义梁etable,Qj,smisc,8 !定义梁ESEL,s,TYPE,2 etabel,tanhuangzhouli,smisc,1 plls,ni,nj,1plls,mi,mj,-14、计算结果衬砌内力弯矩图（ N.m）衬砌内力轴力（ N）5. 衬砌截面强度检算（ 1）隧道结构截面抗压验算：K N <= a Ra b h 式中： K安全系数，取 2.4N轴力（ MN）Ra 混凝土抗压极限强度（ MPa） 构件纵向弯曲系数，取 =1a 轴向力偏心影响系数， 可由 e0/h 值查铁道隧道设 计规范 得到 （ 2）偏心受压构件的抗拉强度验算错误！未找到引用源。 式中： R1 混凝土极限抗拉强度 注：当 e0<=0.2h ，由抗压强度控制承载能