竖井命令流

finish finish /clear,start /filname,process of shujing /title,process of shujing /prep7 r1=5.1 !竖井半径（井壁内径） r2=r1+0.35 !二支半径（井壁外径）,0.35为壁厚 r3=r2+0.12 !初支半径,0.12为初支厚度 r4=r2+2.5 !围岩中锚杆区的半径,2.5为锚杆长度（包括初支和锚杆区） r5=12*r1 !边界范围半径 h=50 !基岩上部竖井深度 et,1,solid45 et,2,surf154 !定义三维 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 面效应单元 （井壁摩擦） KEYOPT,2,2,0 !指定表面效应单元的K2=0,所加荷载与单元坐标系方向相同 KEYOPT,2,4,1 !指定表面效应单元的K4=0,去掉边中点,成为四结点表面单元 !围岩参数 ex1=1500e6 !弹性模量（Pa） nuxy1=0.4 !泊松比 dens1=1900 !密度（kg*m3） c1=100e3 !粘聚力（Pa） jiaodu1=25 !内摩擦家（°度） mp,ex,1,ex1 mp,nuxy,1,nuxy1 mp,dens,1,dens1 tb,dp,1,,,0 !德鲁克——普拉格准则 tbdata,1,c1,jiaodu1,jiaodu1 !锚杆区 ex2=1600e6 nuxy2=0.38 dens2=1900 c2=120e3 jiaodu2=26 mp,ex,2,ex2 mp,nuxy,2,nuxy2 mp,dens,2,dens2 tb,dp,2,,,0 tbdata,1,c2,jiaodu2,jiaodu2 !初支 ex3=2500e6 nuxy3=0.167 dens3=2245 mp,ex,3,ex3 mp,nuxy,3,nuxy3 mp,dens,3,dens3 !井壁 ex4=2800e6 nuxy4=0.1667 dens4=2500 mp,ex,4,ex4 mp,nuxy,4,nuxy4 mp,dens,4,dens4 /view,1,1,1,1 !建立模型 cyl4,0,0,0,0,r1,90,10 !井壁内径 cyl4,0,0,0,0,r2,90,10 !井壁外径 cyl4,0,0,0,0,r3,90,10 !初支区 cyl4,0,0,0,0,r4,90,10 !锚杆区 vptn,all vgen,5,all,,,,,10,,, cyl4,0,0,0,0,r5,90,50 vptn,all nummrg,all numcmp,all allsel !划分网格 local,11,1,,, !摩擦力 csys,11 asel,s,loc,x,r2 aatt,4,,2,0 !面材料4 类型2即surf154 实常数没有 lesize,70,,,10 !圆弧 lesize,96,,,4 !为线指定网格尺寸 将选择出来的面的两条面划分 lesize,22,,,4 lesize,99,,,4 lesize,111,,,4 lesize,117,,,4 mshape,0 !指定单元形状 0是四边形（2d） mshkey,1 !指定划分方式 1是映射 0是自由 amesh,all csys,0 vsel,all vatt,1,1,1 !井筒 lesize,67,,,4 lesize,68,,,4 lesize,66,,,10 !圆弧 lsel,s,,,84,93,3 !体扫掠的预划分 lsel,a,,,9 lesize,all,,,4 allsel vsweep,6 !体扫掠 vsweep,5 vsweep,4 vsweep,3 vsweep,1 !井壁 lesize,69,,,1 lesize,71,,,1 lesize,70,,,10 !圆弧 vsweep,18 vsweep,15 vsweep,9 vsweep,8 vsweep,7 !初支 lesize,72,,,1 lesize,74,,,1 lesize,73,,,10 !圆弧 vsweep,19 vsweep,16 vsweep,13 vsweep,2 vsweep,11 !锚杆区 lesize,75,,,1 lesize,77,,,1 lesize,76,,,10 !圆弧 vsweep,20 vsweep,17 vsweep,14 vsweep,12 vsweep,10 !围岩区 lesize,124,,,10 lesize,125,,,10 lesize,122,,,10 lesize,123,,,10 lesize,79,,,10 !圆弧 vsweep,21 allsel !出现WARNING 网格划分大小有问题 如何设置LIMITs !施加边界条件 asel,s,loc,z,0 !底面 da,all,uz,0 asel,s,loc,x,0 !侧面 da,all,symm !symm表示对称 节点处的法向位移（UX）为0 绕两个切线方向的转角（ROTYZ）也为0 asel,s,loc,y,0 da,all,symm local,11,1,,, !外边界圆弧 csys,11 asel,s,loc,x,r5 nsla,s,1 !选择与面相关的节点 nrotat,all !使节点坐标系方向与当前坐标系方向一致 d,all,ux,0 !只约束UX waring在这个模型上已施加实体模型和有限元模型的边界条件由于实体载荷已转换到节点allsel !或单元上面，它们（指这些载荷）可以直接替换作用力 !!!!! /solu antype,static,new time,1 deltim,0.1,0.05,0.2 !时间步设置,总长0.1,最小0.05,最大0.2 autots,on !使用自动时间步 pred,on !打开时间步长预测器 lnsrch,on !打开线性搜索 nlgeom,on !打开大位移效果 nropt,full !设定牛顿-拉普森选项 sstif,on cnvtol,f,,0.02,2,0.5 ! cnvtol,lab,value,toler,norm,minref toler是个收敛系数value是 !参考值，收敛判断为：如果不平衡力小于等于value*toler，则收敛。 !设定力收敛条件 esel,s,type,,2 !选中实体上表面的表面效应单元 ekill,all esel,all !选择所有元素 esel,s,live !选择活单元,即所有围岩实体单元 nsle,s !选择生单元上的节点 nsel,invert !反向选择,即选择了死单元上的节点 !为什么不直接esel,s,type,,2 nsle,s d,all,all! d,all,all !将死单元上的节点约束所有位移,使其不参与 矩阵运算 nsel,all !选择所有节点 esel,all !选择所有单元 allsel csys,0 acel,0,0,10!自重应力场 time1 solve !开挖第一步 time2 vsel,s,,,6,18,12 eslv !选择已选体（v）上的单元 ekill,all esel,all esel,s,live nsle nsel,invert d,all,all esel,s,type,,2 ealive,all !选中实体上表面的表面效应单元以方便加荷载 sfe,all,1,pres,,-100000!在面内加Z向荷载,大小为100KPa,z为1,x为2 allsel time,2 solve !开挖第二步 time3 vsel,s,,,5,15,10 eslv ekill,all vsel,s,,,18 eslv ealive,all vsel,s,,,18 eslv mpchg,4,all vsel,s,,,19 eslv mpchg,3,all vsel,s,,,20 eslv mpchg,2,all allsel esel,all !选择所有元素 esel,s,live !选择活单元,即所有围岩实体单元 nsle,s !选择生单元上的节点 nsel,invert !反向选择,即选择了死单元上的节点 d,all,all !将死单元上的节点约束所有位移,使其不参与矩阵运算 矩阵运算 nsel,all !选择所有节点 esel,all !选择所有单元 time,3 solve