F,190,FY,-189900 ESEL,S,LIVE !选择具“生/SOLU !进入求解器 属性”的单元 /STATUS,SOLU PLDISP,1 !显示变形ALLS !选择所有内容 情况
SOLVE !求解计算 PLNSOL,U,X,0,1 !绘制水平方向位移云图 !求解 /SOLU !进入求解器 ANTYPE,,REST,2 !重新启动求解器 ESEL,S,,,172,350 !选择上台阶开
挖土体单元
EKILL,ALL !对所选单元赋予“死属性”
ESEL,S,,,5,8 !选择单元（上部开
挖区支护）
ESEL,A,,,1 !补选单元
ESEL,A,,,11,12 ESEL,A,,,16 ESEL,A,,,20,34,2 ESEL,A,,,35,49,2 ESEL,A,,,50 ESEL,A,,,52,53 ESEL,A,,,55,56 ESEL,A,,,59 ESEL,A,,,64,69 ESEL,A,,,73,75 ESEL,A,,,79,81 ESEL,A,,,85,150 ESEL,A,,,154,162 EALIVE,ALL !对所选单元赋予“生属性”
FINISH /SOLU !进入求解器
/STATUS,SOLU ALLS !选择所有
内容
SOLVE !进行求解
计算 save,tunnel-step2,db !把初
始应力模拟求解结果保存 /POST1 !进入后处理器 PLNSOL,U,Y,0,1 !绘制竖直方向位移云图 PLNSOL,S,X,0,1 !绘制水平方向应力云图 PLNSOL,S,Y,0,1 !绘制竖直方向应力云图 PLNSOL,S,1,0,1 !绘制第一主应力云图 PLNSOL,S,3,0,1 !绘制第三主应力云图
FINISH /POST1 !进入后处理器
ETABLE,IMFORX,SMISC, 1 ETABLE,JMFORX,SMISC, 7 ETABLE,IMFORY,SMISC,2 ETABLE,JMFORY,SMISC, 8 ETABLE,IMMOMZ,SMISC,6 ETABLE,JMMOMZ,SMISC,12 PLLS,IMFORX,JMFORX,0.2 !绘制上台阶开挖轴力图 PLLS,IMFORY,JMFORY,0.5 !绘制上台阶开挖剪力图 PLLS,IMMOMZ,JMMOMZ,-0.1 !绘制上台阶开挖弯矩图 ESEL,S,,,331,350 !选择单元 ESEL,A,,,456 !补选单元 ESEL,A,,,620,626,6 ESEL,A,,,459
ESEL,A,,,465,477,6 ESEL,A,,,480,588,6 ESEL,A,,,597
ESEL,A,,,600 ESEL,A,,,609
ESEL,A,,,615 ESEL,A,,,627 ESEL,A,,,621 NFORCE,ALL,ALL!提取所选单元的节点力
FINISH !进行隧道下台阶开挖 /SOLU !进入求解器 ANTYPE,,REST,3 !重新启动求解器 ESEL,S,,,172,350 !选择上台阶开
挖土体单元
EKILL,ALL !对所选单元赋予“死属性” ESEL,S,,,5,8 !选择单元（上部开挖区支护）
ESEL,A,,,1 !补选单元 ESEL,A,,,11,12 ESEL,A,,,16 ESEL,A,,,20,34,2 ESEL,A,,,35,49,2 ESEL,A,,,50 ESEL,A,,,52,53 ESEL,A,,,55,56 ESEL,A,,,59 ESEL,A,,,64,69 ESEL,A,,,73,75 ESEL,A,,,79,81 ESEL,A,,,85,150 ESEL,A,,,154,162
EALIVE,ALL !对所选单元赋予“生属性”
ESEL,S,,,351,449 !选择下台阶开挖土体单元
EKILL,ALL !对所选单元赋予“死属性”
ESEL,S,LIVE !选择具“生属性”的单元 FDELE,3,ALL !删除上下台阶分解处节点上所加的节点力 FDELE,21,ALL
FDELE,172,ALL,190 F,3,FX,-163860 F,4,FX,-139020 F,5,FX,-677.4 F,6,FX,-19.356 F,7,FX,22020 F,8,FX,66900 F,15,FX,-5695.2 F,16,FX,-4648.2 F,17,FX,3754.2 F,18,FX,1912.2 F,22,FX,-67080 F,23,FX,-104520 F,24,FX,-29340 F,25,FX,-9918 F,26,FX,-8268 F,27,FX,-7290 F,30,FX,-3839.4 F,32,FX,-2950.8 F,34,FX,-2502 F,36,FX,-1700.4 F,38,FX,-972.6 F,40,FX,-825.6 F,43,FX,558.12 F,45,FX,1156.2 F,47,FX,1553.4 F,49,FX,2032.8 F,51,FX,2641.2 F,55,FX,977.4 F,59,FX,102540 F,60,FX,153480 F,3,FY,-136740 F,4,FY,464.1 F,5,FY,-80760 F,6,FY,-82680 F,7,FY,-184980 F,8,FY,-97320 F,15,FY,-118320 F,16,FY,-105960 F,17,FY,-137280 F,18,FY,-161220 F,22,FY,-98880 F,23,FY,3067.2 F,24,FY,-172740 F,25,FY,-184500 F,26,FY,-155580 F,27,FY,-133860 F,30,FY,-96660 F,32,FY,-89880 F,34,FY,-85200 F,36,FY,-82020 F,38,FY,-80220 F,40,FY,-79800 F,43,FY,-86040 F,45,FY,-91080 F,47,FY,-98280 F,49,FY,-107940 F,51,FY,-120660 F,55,FY,-194280
F,59,FY,-15474 F,60,FY,40758
/SOLU !进入求解器 /STATUS,SOLU
ALLS !选择所有内容 SOLVE !进行求解计算
/SOLU !进入求解器 ANTYPE,,REST,4 !重新启动求解器 ESEL,S,,,172,350 !选择上台阶开挖土体单元
EKILL,ALL !对所选单元赋予“死属性”
ESEL,S,,,5,8 !选择单元（上部开挖区支护）
ESEL,A,,,1 !补选单元 ESEL,A,,,11,12 ESEL,A,,,16 ESEL,A,,,20,34,2 ESEL,A,,,35,49,2 ESEL,A,,,50 ESEL,A,,,52,53 ESEL,A,,,55,56 ESEL,A,,,59 ESEL,A,,,64,69 ESEL,A,,,73,75 ESEL,A,,,79,81 ESEL,A,,,85,150 ESEL,A,,,154,162
EALIVE,ALL !对所选单元赋予“生属性” ESEL,S,,,351,449 !选择下台阶开挖土体单元
EKILL,ALL !对所选单元赋予“死属性”
ESEL,S,,,2,4 !选择单元（下部开挖区支护） ESEL,A,,,9,10 !补选单元 ESEL,A,,,13,15 ESEL,A,,,17,19 ESEL,A,,,21,33,2
ESEL,A,,,36,48,2 ESEL,A,,,51 ESEL,A,,,54 ESEL,A,,,57,58 ESEL,A,,,60,63 ESEL,A,,,70,72 ESEL,A,,,76,78 ESEL,A,,,82,84 ESEL,A,,,151,153 ESEL,A,,,163,171
EALIVE,ALL !对所选单元赋予“生属性” FINISH
/SOLU !进入求解器 /STATUS,SOLU
ALLS !选择所有内容 SOLVE !进行求解计算
/POST1 !进入后处理器 ESEL,S,LIVE !选择具有“生属性”的单元 PLDISP,1 !显示变形情况 PLNSOL,U,X,0,1 !绘制水平方向位移云图 PLNSOL,U,Y,0,1 !绘制竖直方向位移云图 PLNSOL,S,X,0,1 !绘制水平方向应力云图 PLNSOL,S,Y,0,1 !绘制竖直方向应力云图 PLNSOL,S,1,0,1 !绘制第一主应力云图 PLNSOL,S,3,0,1 !绘制第三主应力云图 ETABLE,IMFORX,SMISC, 1 ETABLE,JMFORX,SMISC, 7 ETABLE,IMFORY,SMISC,2 ETABLE,JMFORY,SMISC, 8 ETABLE,IMMOMZ,SMISC,6 ETABLE,JMMOMZ,SMISC,12
PLLS,IMFORX,JMFORX,0.2 !绘制下台阶开挖后轴力图
PLLS,IMFORY,JMFORY,0.5 !绘制下台阶开挖后剪力图
PLLS,IMMOMZ,JMMOMZ,-0.1 !绘制下台阶开挖后弯矩图
