ansys竖井开挖

本人在做一个地铁竖井施工的模型 本人在做一个地铁竖井施工的模型。计算过程中发现开挖后，应力释放的结果不理想，请各位帮帮忙！！ 附：命令流： 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,154      !定义三维表面效应单元  （井壁摩擦） 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 lesize,70,,,10  !圆弧 lesize,96,,,4 lesize,22,,,4 lesize,99,,,4 lesize,111,,,4 lesize,117,,,4 mshape,0 mshkey,1 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 !施加边界条件 asel,s,loc,z,0 !底面 da,all,uz,0 asel,s,loc,x,0  !侧面 da,all,symm asel,s,loc,y,0 da,all,symm local,11,1,,,          !外边界圆弧 csys,11 asel,s,loc,x,r5 nsla,s,1 nrota,all d,all,ux,0 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     !设定力收敛条件 esel,s,type,,2                !选中实体上表面的表面效应单元 ekill,all esel,all         !选择所有元素 esel,s,live    !选择活单元，即所有围岩实体单元 nsle,s     !选择生单元上的节点 nsel,invert    !反向选择，即选择了死单元上的节点 d,all,all     !将死单元上的节点约束所有位移，使其不参与 矩阵运算 nsel,all            !选择所有节点 esel,all            !选择所有单元 allsel csys,0 acel,0,0,10 !自重应力场time1 solve !开挖第一步time2 vsel,s,,,6,18,12 eslv 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