隧道工程课程设计报告.doc

《隧道工程课程设计报告.doc》由会员分享，可在线阅读，更多相关《隧道工程课程设计报告.doc（13页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

2、的相互作用是通过弹性支承对支护构造施加约束来实现的。图1 计算模型计算模型中，二衬构造采用弹性平面梁单元模拟，弹性抗力以及隧底地基均采用弹簧单元模拟。组合荷载根据不同作用方向分别转换成等效节点力施加在相应的单元结点上。具体计算模型见图1。2、荷载计算围岩压力计算参照课本中有关我国铁路隧道推荐的方法进展确定双线隧道或参照?铁 路隧道设计标准?，深浅埋分别计算。按破坏阶段设计计算垂直压力公式：q=r x hq = 0.45 x 2(s-1) x r x w式中： hq等效荷载高度值 S围岩级别 r围岩的容重 w宽度影响系数，其值为w=1+iB-505=2.082816e6N/m水平均布松动压力系数

3、取0.3，那么e=0.3q=N/m3、ANSYS操作命令流 !荷载构造方法计算(马蹄形断面finish !退出当前处理程序 /clear !去除以前数据,重新开场一个新的分析/COM,Structural !定义分析类型,构造分析(热分析、流体分析等) /prep7 !进入前处理器*AFUN,deg !定义角度单位为度缺省为弧度，RAD! 定义建模及材料参数的一些变量值*set,Py,2.082816e5 !定义垂直围岩压力大小假设有地表荷载加地表荷载值*set,cylxsh,0.3 !定义侧压力系数*set,cyl,Py*cylxsh !水平侧压力*set,CQHD,0.47 !定义初支或二

4、衬厚度*set,CQDYCD,0.17 !定义梁单元长度参数及弹簧单元面积梁单元长度与弹簧单元面积相等*set,CQETXML,31e9 !定义衬砌初支或二衬的弹性模量*set,CQUBSB,0.2 !定义衬砌的泊松比*set,WYTXKL,100e6 !定义围岩的弹性抗力系数*set,WYMD,1600 !定义围岩的密度*set,CQMD,2500 !定义衬砌的密度!定义单元类型及材料属性及单元实常数et,1,beam3 !定义1号单元为梁单元mp,ex,1,CQETXML !定义1号材料的弹性模量 mp,prxy,1,CQUBSB !定义1号材料的泊松比mp,dens,1,CQMD !定义

5、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,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.

6、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,3,2 !三点画圆弧：起点-终点中间点 最好全按逆时针画，弯矩及轴力的方向一致larc,3,5,4larc,5,7,6larc,7,9,8larc,9,11,10larc,11,13,12larc,13,15,14larc,15,1,16!设置线单元材料属性，划分单元(初支单元)lsel,s,1,8 !L

7、SEL,Type,Item,Comp,VMIN, VMAX,VINClATT,1,1,1 !给线单元付材料号、实常数、单元类型号 LATT, MAT, REAL, TYPE, ESYSMSHKEY,1 !设置单元为映射单元单元比拟规那么lesize,1,CQDYCD !LESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACElesize,2,CQDYCDlesize,3,CQDYCDlesize,4,CQDYCDlesize,5,CQDYCDlesize,6,CQDYCDlesize,7,CQDYCDlesize,8,CQDYCDlsel,alllmesh,all !划分线

8、单元 nummrg,all,low !所有号数从1开场numcmp,node,eord !压缩节点号并排序!施加水平弹簧支撑单元在梁单元的每个节点上分别产生水平方向弹簧单元*do,i,2,118 !do循环，要对应下面一个*EDNDOa1=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循环，

9、要对应下面一个*EDNDOa1=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 !弹簧单元,节点号,平移,弹性系数,弹簧水平长度 产生水平弹簧*enddoa1=abs(ny(235)-ny(1)*WYTXKL !计算节点294上的弹性抗力kx=a1/2PSPRNG,236,TRAN,kx,-0.2,!施加竖直弹簧支撑单元在梁单元的每个节点上分别产生竖直方向弹簧单元b1=abs(nx(236)-nx

10、(2)*WYTXKL !计算节点1上的弹性抗力ky=b1/2PSPRNG,1,TRAN,ky,-0.2,*do,i,2,57 !do循环，要对应下面一个*EDNDOb1=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循环，要对应下面一个*EDNDOb1=abs(nx(i-1)-nx(i+1)

11、*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,181,235 !do循环，要对应下面一个*EDNDOb1=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

12、 !弹簧单元,节点号,平移,弹性系数,弹簧水平长度 产生水平弹簧*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 !在59181号节点施加竖向围岩压力a=-Py*0.5*abs(nx(i-1)

13、-nx(i+1) !计算作用在节点上的节点力f,i,fy,a*enddo !水平向荷载X方向压力，全部节点施加水平荷载，换算成节点荷载，在节点上施加*do,i,2,118 !在2118号节点施加水平围岩压力f,i,fX,a*enddo*do,i,119,235 !在1120235号节点施加水平围岩压力f,i,fX,a*enddof,236,fX,a!设定重力加速度ACEL,0,9.8,0, !重力加速度以向下位置，不用加负号ALLSEL,ALL finish!删除轴力为拉力的弹簧单元/prep7 *do,i,239,261 EDELE, iNDELE, i*enddo*do,i,297,415

14、 EDELE, iNDELE, i*enddo *do,i,451,471 EDELE, iNDELE, i*enddo *do,i,529,651 EDELE, 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,ALLso

15、lve !求解save,jg,db !将第一时间步的计算结果保存!进入后处理器，对结果进展后处理/POST1ESEL,s,TYPE,1NSLE,s,ALL etabel,ni,smisc,1 !定义梁i节点的轴力 etabel,nj,smisc,7 !定义梁j节点的轴力etabel,mi,smisc,6 !定义梁i节点的弯矩etabel,mj,smisc,12 !定义梁j节点的弯矩etable,Qi,smisc,2 !定义梁i节点的剪力etable,Qj,smisc,8 !定义梁j节点的剪力ESEL,s,TYPE,2etabel,tanhuangzhouli,smisc,1 plls,ni,n

16、j,1plls,mi,mj,-14、计算结果 衬砌内力弯矩图N.m 衬砌内力轴力N5.衬砌截面强度检算1隧道构造截面抗压验算： K N = a Ra b h 式中：K N轴力MN Ra混凝土抗压极限强度MPa 构件纵向弯曲系数，取=1 a轴向力偏心影响系数，可由e0/h值查?铁道隧道设 计标准?得到2偏心受压构件的抗拉强度验算 K N = j式中： R1混凝土极限抗拉强度注：当e00.2h，由抗拉强度控制承载能力，不必验算抗压。3验算最大弯矩截面eo = M / N = 2514E+06/2381E67 0.2h = 0.2*0.42 =0.084由抗拉强度控制承载能力, 不必验算抗压。查上表10.2.3 取 j = 1.0；K N = j420)/(6x-1)=e6那么K = e6/2381e6=3.6 所以衬砌不满足截面抗拉要求。第 13 页