搜索
（1）创建物理环境/COM,Structural!指定结构分析/TITLE,TunnelConstructModelingAnalysis!定义工作标题/FILNAM,support,1!定义工作文件名（2）建立模型!进入前处理器/PREP7!定义分析参数*set,jd_s,0*set,jd_e,180-jd_s*set,num,30*set,jd,(jd_e-jd_s)/num*set,distance,10!溶洞与隧道距离*set,depth,80!隧道埋深*set,fricangle,37*set,dens,2200*set,d,12.8*set,cohesion,0.6E6*set,posionratio,0.32*set,elasticmoduli,3.6E9*set,t1,0.25*set,r_karst,3.6!溶洞半径*set,r1,2.5E-2*afun,deg!定义单元类型ET,1,BEAM3!定义衬砌支护单元KEYOPT,1,6,1ET,2,PLANE42!定义围岩单元KEYOPT,2,3,2ET,3,LINK1!定义LINK1单元!定义材料属性!衬砌支护MP,EX,1,2.95E10MP,PRXY,1,0.2MP,DENS,1,2500TB,DP,1!采用D-P模型TBDATA,1,2.42E6,54!围岩材料MP,EX,2,3.69E9MP,PRXY,2,0.32MP,DENS,2,2200TB,DP,2TBDATA,1,0.6E6,37!挖去土体材料MP,EX,3,3.69E9MP,PRXY,3,0.32MP,DENS,3,2200TB,DP,3TBDATA,1,0.6E6,37!锚杆MP,EX,4,17E10MP,PRXY,4,0.3MP,DENS,4,7960!定义实常数R,1,0.3,0.3*0.3*0.3/12,0.3,!衬砌支护实常数R,2,3.14*0.025*0.025/4,,!锚杆实常数!建立几何模型!创建隧道衬砌支护线K,1,,,,!创建关键点K,2,,11.37,,K,3,3.34,0,,K,4,-3.34,0,,K,201,0,,10,K,202,,11.37,10,K,203,3.34,0,10,K,204,-3.34,0,10,!创建圆circle,1,5.8,201circle,2,14.4,202circle,3,2.67,203circle,4,2.67,204lcsl,all!把四个圆在交点处打断lsel,s,line,,17,22,1!选择线lsel,a,line,,27,28,1lsel,a,line,,43lsel,a,line,,46,48,1cm,zh,line!定义一个线部件cmsel,s,zh,line!选择部件lsel,inve!反选当前的线ldel,all,,,1!删除所有选择的线!创建初期支护加固范围circle,1,8.8,201!画圆circle,2,17.4,202circle,3,5.67,203circle,4,5.67,204cmsel,s,zh,line!选择线部件lsel,invelcsl,alllsel,s,line,,5,6,1lsel,a,line,,31,34,1lsel,a,line,,37,42,1lsel,a,line,,44,45,1lsel,a,line,,49,58,1lsel,a,line,,60,61,1lsel,a,line,,64,65,1ldel,all,,,1allsel!绘制溶洞CYL4,0,-12.8,3.6!绘制圆实体!细分隧道分析线模型k,60,-52,-65!创建关键点k,61,-12,-65k,62,12,-65k,63,52,-65k,64,52,-17.6k,65,12,-17.6k,66,-12,-17.6k,67,-52,-17.6k,68,-52,-8k,69,-12,-8k,70,12,-8k,71,52,-8k,72,52,10k,73,12,10k,74,-12,10k,75,-52,10k,76,-52,80k,77,-12,80k,78,12,80k,79,52,80l,60,61!连接关键点60、61生成直线l,61,62l,62,63l,63,64l,64,65l,65,66l,66,67l,67,68l,68,69l,69,70l,70,71l,71,72l,72,73l,73,74l,74,75l,75,76l,76,77l,77,78l,78,79l,79,72l,78,73l,77,74l,75,68l,74,69l,73,70l,71,64l,70,65l,69,66l,67,60l,66,61l,65,62a,60,63,79,76!通过4个关键点生成一个面积区域asbl,all,all!通过线分割面生成新面adele,1!删除面NUMCMP,AREA!压缩面编号!生成锚杆asel,s,,,14!选择面14csys,1!激活柱坐标系wprota,,-90!工作平面绕X轴旋转-90度*do,i,1,30,1!循环控制wprota,,,-6!工作平面绕Y轴旋转-6度asbw,all!用工作平面切割所选择的所有面*enddoasel,s,,,14,15,1!选择面asel,a,,,19,20,1asel,a,,,22,25,3asel,a,,,26,28,2asel,a,,,31,32,1asel,a,,,34,50,2asel,a,,,53,54,1asel,a,,,56,59,3asel,a,,,60,62,2asel,a,,,65,67,1asel,a,,,69,71,2asel,a,,,75aadd,all!删除所选择的所有面allsel!保存几何模型save,Tunel-geom.db!划分网格生成有限元模型!划分梁单元mat,1!指定梁单元材料特性type,1asel,s,,,11lsla,s,esize,1!指定划分梁单元长度lmesh,all!划分所有线allsel!划分锚杆lsel,s,,,26,48,22!选择锚杆线lsel,a,,,61,71,10lsel,a,,,77,89,12lsel,a,,,92,95,3lsel,a,,,107,113,3lsel,a,,,125,131,6lsel,a,,,137,139,2lsel,a,,,143,151,8lsel,a,,,157,169,6lsel,a,,,172,175,3lsel,a,,,187,193,3lsel,a,,,205,211,3lsel,a,,,212,213,1mat,4!指定锚杆单元材料特性type,3lmesh,allallsel!划分开挖掉土体单元网格mat,3type,2mshkey,0!设定自由网格划分mshape,0!设定四边形网格划分esize,amesh,11!划分11号面!划分围岩网格!设置网格份数lsel,s,,,5,11,1!选择线lsel,a,,,13,15,1lsel,a,,,31,33,1lsel,a,,,37,39,1lsel,a,,,53,55,1lesize,all,,,6!把所选择线分为6份lsel,s,,,1,4,1lsel,a,,,50,52,1lsel,a,,,12lesize,all,,,4!把所选择线分为4份lsel,s,,,16,34,18lsel,a,,,40,42,1lsel,a,,,44,45,1lsel,a,,,49lesize,all,,,8!把所选择线分为8份mat,3!赋予围岩单元属性type,2mshkey,0!设定自由网格划分mshape,0!设定四边形网格划分allselasel,s,,,11!选择面11asel,inve!反选当前面amesh,all!划分所有面allsel!保存网格模型save,Tunel-grid.db（3）施加约束和荷载!施加约束csys,0nsel,s,loc,x,-52!选择X=-52线上所有节点nsel,a,loc,x,52!选择X=52线上所有节点d,all,ux!对所选择节点约束X方向位移allselnsel,s,loc,y,-65!选择Y=--65线上所有节点d,all,uy!对所选择节点约束Y方向位移!施加重力加速度acel,,9.8（4）求解/solu!求解设置antype,static!设定为静力求解nsubst,100!设定最大子步数为100pred,on!打开时间步长预测器nropt,full!设定牛顿-拉普森选项nlgeom,on!打开大位移效果lnsrch,on!打开线性搜索outres,all,all!输出所有项cnvtol,f,,0.02,2,0.5!力收敛准则设定cnvtol,m,,0.01,2,1!力矩收敛准则设定!初始应力模拟time,1!设定载荷步结束时间allselesel,s,mat,,1,4,3!选择材料号为1、4的单元ekill,all!杀死所选择的单元esel,allesel,s,live!选择活单元nsle,s!选择节点nsel,inve!反选择当前节点d,all,all!约束所选择节点、自由度allselsolve!求解save,tunnel-step1,db!把初始应力模拟求解结果保存!隧道开挖与支护模拟!把开挖模拟求解结果保存!列车荷载模拟time,3!设定载荷步结束时间f,6,fy,-525000!施加集中力f,41,fy,-525000allselsolve!求解save,Tunel-step3,db!把列车模拟求解结果保存（5）后处理/post1!初始应力模拟结果分析Resume,'tunnel-step1','db'!读入初始应力模拟数据set,1,last!读入本荷载步最后一个子步esel,s,live!选择活单元!显示位移云图plnsol,u,sum,!绘制总位移矢量云图plnsol,u,x!绘制X方向位移云图plnsol,u,y!绘制Y方向位移云图!显示应力云图plnsol,s,x!绘制X方向应力云图plnsol,s,y!绘制Y方向应力云图plnsol,s,z!绘制Z方向应力云图plnsol,s,1!绘制第1主应力云图plnsol,s,2!绘制第2主应力云图plnsol,s,3!绘制第3主应力云图plnsol,s,eqv!绘制等效应力云图!开挖模拟结果分析Resume,'tunnel-step2','db'!读入开挖模拟数据set,1,last!读入本荷载步最后一个子步esel,s,live!显示位移云图plnsol,u,sum,!绘制总位移矢量云图plnsol,u,x!绘制X方向位移云图plnsol,u,y!显示应力云图plnsol,s,x!绘制X方向应力云图plnsol,s,y!绘制Y方向应力云图plnsol,s,zplnsol,s,1plnsol,s,2!绘制第2主应力云图plnsol,s,3plnsol,s,eqv!绘制等效应力云图!显示梁支护内力esel,s,mat,,1!选择材料号为1的单元etable,IMOMEMT,SMISC,6!定义弯矩表etable,JMOMEMT,SMISC,12etable,ISHEAR,SMISC,2!定义剪力表etable,JSHEAR,SMISC,8etable,ZHOULI-I,SMISC,1!定义轴力表etable,ZHOULI-I,SMISC,7/TITLE,BENDINGMOMENTdistribution!定义弯矩分布标题PLLS,IMOMEMT,JMOMEMT,-0.5!绘制弯矩分布图/TITLE,SHEARforcedistribution!定义剪力分布标题PLLS,ISHEAR,JSHEART,0.5!绘制剪力分布图/TITLE,ZHOULIforcedistribution!定义轴力分布标题PLLS,ZHOULI-I,ZHOULI-J,0.1!绘制轴力分布图!显示锚杆内力esel,s,mat,,4!选择材料号为4的单元etab