/NOPR
KEYW,PR_STRUC,1 /PREP7 ET,1,BEAM3 ET,2,PLANE42 KEYOPT,2,3,3 R,1,0.15,0.000281,0.15 !菜单过滤设置 !保留结构分析部分菜单 !进入前处理器
!设置了梁单元类型，模拟喷射混凝土 !设置平面实体单元类型，模拟围岩 !设置平面应变模式 !设置梁单元几何常数 MP,EX,1,23E9 MP,PRXY,1,0.2 MP,DENS,1,2500 MP,EX,2,8E9 MP,PRXY,2,0.28 MP,DENS,2,2400 TB,DP,2
TBMODIF,1,1,1.2E6 TBMODIF,1,2,42 MP,EX,3,16E9 MP,PRXY,3,0.25 MP,DENS,3,2500 TB,DP,3
TBMODIF,1,1,2.4E6 TBMODIF,1,2,50 SAVE K,1,21,0,0
K,2,15,-1,0 K,3,27,-1,0 K,4,0,0,0 K,5,0,21,0 K,6,21,21,0 K,7,42,21,0
!输入弹性模量 !输入泊松比 !输入密度 !输入弹性模量 !输入泊松比 !输入密度
!输入 DP 模型的相应参数 !输入粘聚力 !输入内摩擦角
!以下输入加固圈（锚杆加固范围）参数
K,8,42,0,0
K,9,55,0,0 K,10,55,21,0 K,11,55,56,0 K,12,0,56,0 K,13,21,56,0 K,14,42,56,0 K,15,55,-14,0 K,16,42,-14,0 K,17,27,-14,0
K,18,0,-14,0 K,19,14.9896,6.0104 K,20,27.0104,6.0104 KSYMM,X,ALL,,,,0,0 NUMMRG,ALL,,,,LOW NUMCMP,ALL SAVE
K,,15,-14
K,,-15,-14 CIRCLE,1,6 CIRCLE,21,6 CIRCLE,21,8.5 CIRCLE,1,8.5 LDELE,3,4,1,1 LDELE,7,8,1,1 LDELE,11,12,1,1 LDELE,15,16,1,1 NUMMRG,ALL,,,,LOW NUMCMP,ALL SAVE
!创建将要划分六段的线对称复制所有关键点 !合并所有关键点 !压缩所有关键点编号!
L,31,30 L,30,12 L,12,13 L,13,14 L,29,28 L,31,25 L,30,24 L,12,5
L,13,6 L,14,7 L,11,10 L,25,24 L,24,5
L,5,6 L,6,7 L,28,27 L,25,26 L,25,36 L,24,46 L,43,21 L,35,5 L,5,4
L,5,19 L,6,49 L,40,1 L,20,7 L,7,8 L,9,10
L,8,48 L,4,50 L,4,45
L,26,47 L,27,32 L,26,33 L,4,18
L,8,16 L,9,15
!创建将要划分三段的线 L,29,31 L,14,11 L,28,25 L,7,10 L,27,26 L,8,9 L,33,32 L,15,16
!绘制将要划分五段的线 L,23,34 L,22,38 L,2,37 L,3,17
!绘制将要划分二段的线 L,44,47 L,43,46 L,42,45 L,50,41 L,40,49 L,39,48
!绘制将要划分一段的线
L,23,44 L,22,42 L,41,2
L,39,3
!绘制将要划分十二段的线 L,44,21 L,21,42 L,41,1 L,1,39
!绘制将要划分二十四段的线 L,22,23 L,34,38 L,2,3 L,37,17
!绘制将要划分八段的线 L,33,34 L,38,18 L,18,37 L,17,16 SAVE
!通过线创建面 AL,22,23,35,37,7,8,38,30 LSEL,S,LINE,,31,38,7 LSEL,A,LINE,,8 ASBL,1,ALL LSEL,S,LINE,,31,34,3 LSEL,A,LINE,,7,81,74 ASBL,3,ALL ALLSEL,ALL AL,20,21,30,39,5,6,40,25 LSEL,S,LINE,,26,40,14 LSEL,A,LINE,,6 ASBL,3,ALL LSEL,S,LINE,,26,29,3
LSEL,A,LINE,,5,83,78 ASBL,6,ALL NUMMRG,ALL,,,,LOW NUMCMP,ALL SAVE ALLSEL,ALL AL,2,32,69 AL,32,70,1 AL,79,80,61,2,60 AL,61,81,7,62,1 AL,37,60,65,55,77,42 AL,69,70,66,73,65 AL,73,56,74,55 AL,62,36,43,78,56,66 AL,50,44,52,43 AL,48,35,50,34 AL,46,18,48,17 AL,11,17,22,16 AL,10,16,21,15 AL,9,15,20,14 AL,8,14,19,13 AL,45,13,47,12 AL,47,24,49,23 AL,4,27,67 AL,27,68,3 AL,6,82,58,4,57 AL,83,5,59,3,58 AL,49,41,51,40 AL,39,57,63,53,75,41 AL,67,68,64,71,63 AL,71,54,72,53
AL,64,59,38,42,76,54 SAVE
!单元网格划分 ALLSEL,ALL LSEL,S,LINE,,5,44,1 LSEL,A,LINE,,79,83,1 LSEL,U,LINE,,27,32,5 LESIZE,ALL,,,6 LSEL,S,LINE,,45,52,1 LESIZE,ALL,,,3 LSEL,S,LINE,,53,56,1 LESIZE,ALL,,,5 LSEL,S,LINE,,57,62,1 LESIZE,ALL,,,2 LSEL,S,LINE,,63,66,1 LESIZE,ALL,,,1 LSEL,S,LINE,,67,70,1 LSEL,A,LINE,,1,4,1 LSEL,A,LINE,,27,32,5 LESIZE,ALL,,,12 LSEL,S,LINE,,71,74,1 LESIZE,ALL,,,24 LSEL,S,LINE,,75,78,1 LESIZE,ALL,,,8 SAVE
!加固圈单元划分 ALLSEL,ALL TYPE,2 MAT,3 REAL,1 LCCAT,79,80
LCCAT,81,7 LCCAT,6,82 LCCAT,5,83 AMAP,26,47,46,43,44 AMAP,27,46,45,42,43 AMAP,9,50,49,40,41 AMAP,10,49,48,39,40 LSEL,S,LCCAT LDELE,ALL
!以下划分开挖土体网格 ALLSEL,ALL TYPE,2 MAT,2 REAL,1
AMAP,24,44,43,21 AMAP,25,43,42,21 AMAP,7,41,40,1 AMAP,8,40,39,1 LCCAT,67,68 LCCAT,69,70 AMAP,30,44,42,22,23 AMAP,12,41,39,3,2 LSEL,S,LCCAT LDELE,ALL
!划分未开挖土体网格 ALLSEL,ALL LCCAT,19,20 LCCAT,82,83 LCCAT,21,22 LCCAT,80,81 AMAP,6,25,5,35,36
AMAP,4,5,7,20,19 LSEL,S,LCCAT LSELE,ALL AMAP,2,5,19,50,4 AMAP,1,20,7,8,48 AMAP,13,2,3,17,37 AMAP,15,8,9,15,16 AMAP,16,7,10,9,8 AMAP,17,14,11,10,7 AMAP,18,13,14,7,6 AMAP,19,12,13,6,5 AMAP,20,30,12,5,24 AMAP,21,31,30,24,25 AMAP,22,29,31,25,28 AMAP,23,28,25,26,27 AMAP,5,25,36,47,26 AMAP,3,35,5,4,45 AMAP,28,27,26,33,32 AMAP,31,23,22,38,34 LCCAT,39,57 LCCAT,63,53 LCCAT,54,64 LCCAT,59,38 LCCAT,37,60 LCCAT,55,65 LCCAT,56,66 LCCAT,62,36 AMAP,29,26,44,34,33 AMAP,32,42,4,18,38 AMAP,11,4,41,37,18 AMAP,14,39,8,16,17
LSEL,S,LCCAT LDELE,ALL SAVE
!生成初期支护喷射混凝土梁单元 TYPE,1 MAT,1 REAL,1
!右隧道 E,608,127 E,127,129 E,129,130 E,130,131 E,131,132 E,132,133 E,133,134 E,134,135 E,135,136 E,136,137 E,137,138 E,138,139 E,139,91 E,91,92 E,92,93 E,93,94 E,94,95 E,95,96 E,96,97 E,97,98 E,98,99 E,99,100 E,100,101
E,101,102 E,102,89 E,89,609
!左隧道 E,583,52 E,52,53 E,53,54 E,54,55 E,55,56 E,56,57 E,57,58 E,58,59 E,59,60 E,60,61 E,61,62 E,62,63 E,63,16 E,16,17 E,17,18 E,18,19 E,19,20 E,20,21 E,21,22 E,22,23 E,23,24 E,24,25 E,25,26 E,26,27 E,27,14 E,14,584 SAVE
!加载与初始地应力的求解
!对左右两侧边界施加 UX 方向的约束，在下边界施加 UY 方向的约束
/SOL ALLSEL,ALL
NSEL,S,LOC,X,-55.1,-54.9
NSEL,A,LOC,X,54.9,55.1 D,ALL,UX,0
NSEL,S,LOC,Y,-14.1,-13.9 D,ALL,UY,0 SAVE
!施加重力加速度 ACEL,0,10,0 SAVE !杀死梁单元 ESEL,S,TYPE,,1 EKILL,ALL SAVE !求解 !求解前设置 NROPT,FULL DELTIM,100,20,1000 ALLSEL,ALL OUTRES,ALL,ALL !初始应力求解 SOLVE FINISH SAVE
!自重应力场处理
/POST1 PLNSOL,U,Y,0,1
采用全牛顿—拉普森法进行求解 设置求解次数
!竖向位移
!!PLNSOL,S,2,0,1 PLNSOL,S,3,0,1 !左隧道土体开挖施工 !开挖左侧隧道土体 FINISH /SOL ANTYPE,,REST,1,1,0 ASEL,S,AREA,,24,25,1 ASEL,A,AREA,,30 ESLA,S EKILL,ALL
!第一主应力 !第三主应力
!进行初期支护，即激活喷射混凝土衬砌梁单元 ALLSEL,ALL
ESEL,S,ELEM,,1549,1574,1 EKILL,ALL
!杀死右隧道喷射混凝土梁单元 ALLSEL,ALL ESEL,S,ELEM,,1575,1600,1 EALIVE,ALL SAVE FINISH !施加节点荷载 /POST1 ALLSEL,ALL
!进入后处理器 !选择所有内容
ESEL,ALL ESEL,S,LIVE ESEL,S,P NSEL,S,P NFORCE,ALL FINISH /SOL
!人工手动选择开挖土体周围的一圈单元 !对选择的节点求解节点力
!施加节点力
F, 14 ,FX, F, 16 ,FX, F, 17 ,FX, F, 18 ,FX, F, 19 ,FX, F, 20 ,FX, F, 21 ,FX, F, 22 ,FX, F, 23 ,FX, F, 24 ,FX, F, 25 ,FX, F, 26 ,FX, F, 27 ,FX, F, 52 ,FX, F, 53 ,FX, F, 54 ,FX, F, 55 ,FX, F, 56 ,FX, F, 57 ,FX, F, 58 ,FX, F, 59 ,FX, F, 60 ,FX, F, 61 ,FX, F, 62 ,FX, F, 63 ,FX, F, 583 ,FX, F, 584 ,FX, F, 585 ,FX, F, 586 ,FX, F, 587 ,FX,
2.83E+04 7.48E+01 3.27E+03 6.42E+03 9.48E+03 1.24E+04 1.51E+04 1.75E+04 1.99E+04 2.19E+04 2.38E+04 2.51E+04 2.47E+04 -2.82E+04 -2.47E+04 -2.50E+04 -2.36E+04 -2.18E+04 -1.97E+04 -1.74E+04 -1.49E+04 -1.22E+04 -9.30E+03 -6.30E+03 -3.12E+03 -1.39E+04 1.39E+04 4.18E+03 2.19E+03 1.54E+03
F, 588 ,FX, F, 589 ,FX, F, 590 ,FX, F, 591 ,FX, F, 592 ,FX, F, 593 ,FX, F, 594 ,FX, F, 595 ,FX, F, 596 ,FX, F, 597 ,FX, F, 598 ,FX, F, 599 ,FX, F, 600 ,FX, F, 601 ,FX, F, 602 ,FX, F, 603 ,FX, F, 604 ,FX, F, 605 ,FX, F, 606 ,FX, F, 607 ,FX,
1.07E+03 7.81E+02 5.70E+02 4.13E+02 2.88E+02 1.86E+02 9.81E+01 1.85E+01 -5.62E+01 -1.31E+02 -2.11E+02 -2.99E+02 -4.02E+02 -5.28E+02 -6.86E+02 -8.98E+02 -1.19E+03 -1.66E+03 -2.31E+03 -4.31E+03
F, 14 ,FY, -2.06E+04 F, 16 ,FY, -8.10E+04
F, 17 ,FY, -8.04E+04 F, 18 ,FY, -7.86E+04 F, 19 ,FY, -7.56E+04 F, 20 ,FY, -7.08E+04 F, 21 ,FY, -6.54E+04 F, 22 ,FY, -5.88E+04 F, 23 ,FY, -5.12E+04 F, 24 ,FY, -4.28E+04 F, 25 ,FY, -3.40E+04
F, 26 ,FY, -2.55E+04
F, 27 ,FY, -1.85E+04 F, 52 ,FY, -2.04E+04 F, 53 ,FY, -1.83E+04 F, 54 ,FY, -2.53E+04 F, 55 ,FY, -3.38E+04 F, 56 ,FY, -4.27E+04 F, 57 ,FY, -5.10E+04 F, 58 ,FY, -5.87E+04 F, 59 ,FY, -6.54E+04 F, 60 ,FY, -7.08E+04 F, 61 ,FY, -7.56E+04 F, 62 ,FY, -7.86E+04 F, 63 ,FY, -8.04E+04 F, 583 ,FY, 2.86E+04 F, 584 ,FY, 2.85E+04 F, 585 ,FY, 6.60E+04 F, 586 ,FY, 6.24E+04 F, 587 ,FY, 6.18E+04 F, 588 ,FY, 6.12E+04 F, 589 ,FY, 6.12E+04 F, 590 ,FY, 6.06E+04 F, 591 ,FY, 6.06E+04 F, 592 ,FY, 6.06E+04 F, 593 ,FY, 6.06E+04 F, 594 ,FY, 6.06E+04 F, 595 ,FY, 6.06E+04 F, 596 ,FY, 6.06E+04 F, 597 ,FY, 6.06E+04 F, 598 ,FY, 6.06E+04 F, 599 ,FY, 6.06E+04
F, 600 ,FY, 6.06E+04
F, 601 ,FY, 6.06E+04 F, 602 ,FY, 6.12E+04 F, 603 ,FY, 6.12E+04 F, 604 ,FY, 6.12E+04 F, 605 ,FY, 6.18E+04 F, 606 ,FY, 6.24E+04
F, 607 ,FY, 6.60E+04 !左隧道开挖求解 ALLSEL,ALL SOLVE FINISH SAVE
!左隧道开挖模拟分析后处理
/POST1 PLNSOL,U,Y,0,1 PLNSOL,S,Y,0,1 PLNSOL,S,1,0,1 PLNSOL,S,3,0,1 PLNSOL,EPPL,EQV,0,1 ESEL,S,LIVE ESEL,R,TYPE,,1 ETABLE,,SMISC,6
ETABLE,,SMISC,12
ETABLE,,SMISC,1 ETABLE,,SMISC,7 ETABLE,,SMISC,2 ETABLE,,SMISC,8 PLLS,SMIS6,SMIS12,1,0 PLLS,SMIS1,SMIS7,1,0 PLLS,SMIS2,SMIS8,1,0
方向的位移 !Y 方向的应力 !第一主应力 !第三主应力 !塑性区
!6 和 12 表示弯矩!1 和 7 表示轴力
!2 和 8 表示剪力 !绘制弯矩图 !绘制轴力图 !绘制剪力图
!Y !右隧道土体开挖施工 !开挖右隧道土体 SAVE /SOL
ANTYPE,,REST,2,1,0 ASEL,S,AREA,, 7,8,1 ASEL,A,AREA,,12 ESLA,S EKILL,ALL
ASEL,S,AREA,,24,25,1
!重新启动求解器 !选择右隧道面
!选择右隧道土体单元 !对选择的单元给予“死单元”
!选择左隧道面，每次计算都要重新杀死单元 !选择左隧道土体单元 对选择的单元给予“死单元”
ASEL,A,AREA,,30 ESLA,S EKILL,ALL
ESEL,S,ELEM,,1549,1574,1 EALIVE,ALL SAVE FINISH !施加节点荷载 !节点力计算
!进行初期支护，即激活喷射混凝土衬砌梁单元
/POST1 ALLSEL,ALL ESEL,ALL ESEL,S,LIVE ESEL,S,P
!人工手动选择开挖土体周围的一圈单元
NFORCE,ALL
FINISH /SOL !节点力施加 F, 89 ,FX, F, 91 ,FX,
!进入求解器
2.82E+04 -7.48E+01
F, 92 ,FX, F, 93 ,FX, F, 94 ,FX, F, 95 ,FX, F, 96 ,FX, F, 97 ,FX, F, 98 ,FX, F, 99 ,FX, F, 100 ,FX, F, 101 ,FX, F, 102 ,FX, F, 127 ,FX, F, 129 ,FX, F, 130 ,FX, F, 131 ,FX, F, 132 ,FX, F, 133 ,FX, F, 134 ,FX, F, 135 ,FX, F, 136 ,FX, F, 137 ,FX, F, 138 ,FX, F, 139 ,FX, F, 608 ,FX, F, 609 ,FX, F, 610 ,FX, F, 611 ,FX, F, 612 ,FX, F, 613 ,FX, F, 614 ,FX, F, 615 ,FX,
3.12E+03 6.27E+03 9.31E+03 1.22E+04 1.49E+04 1.74E+04 1.97E+04 2.18E+04 2.37E+04 2.50E+04 2.47E+04 -2.83E+04 -2.47E+04 -2.51E+04 -2.37E+04 -2.19E+04 -1.98E+04 -1.75E+04 -1.51E+04 -1.24E+04 -9.46E+03 -6.42E+03 -3.27E+03 -1.39E+04 1.39E+04 -4.18E+03 -2.19E+03 -1.54E+03 -1.07E+03 -7.81E+02 -5.70E+02
F, 616 ,FX, F, 617 ,FX, F, 618 ,FX, F, 619 ,FX, F, 620 ,FX, F, 621 ,FX, F, 622 ,FX, F, 623 ,FX, F, 624 ,FX, F, 625 ,FX, F, 626 ,FX, F, 627 ,FX, F, 628 ,FX, F, 629 ,FX, F, 630 ,FX, F, 631 ,FX, F, 632 ,FX,
-4.13E+02 -2.88E+02 -1.86E+02 -9.81E+01 -1.85E+01 5.62E+01 1.31E+02 2.11E+02 2.99E+02 4.02E+02 5.28E+02 6.86E+02 8.98E+02 1.19E+03 1.66E+03 2.31E+03 4.31E+03
F, 89 ,FY, -2.04E+04 F, 91 ,FY, -8.12E+04
F, 92 ,FY, -8.05E+04 F, 93 ,FY, -7.86E+04 F, 94 ,FY, -7.54E+04 F, 95 ,FY, -7.09E+04 F, 96 ,FY, -6.53E+04 F, 97 ,FY, -5.87E+04 F, 98 ,FY, -5.10E+04 F, 99 ,FY, -4.26E+04 F, 100 ,FY, -3.38E+04 F, 101 ,FY, -2.53E+04 F, 102 ,FY, -1.83E+04 F, 127 ,FY, -2.06E+04
F, 129 ,FY, -1.85E+04
F, 130 ,FY, -2.55E+04 F, 131 ,FY, -3.40E+04 F, 132 ,FY, -4.28E+04 F, 133 ,FY, -5.12E+04 F, 134 ,FY, -5.88E+04 F, 135 ,FY, -6.55E+04 F, 136 ,FY, -7.10E+04 F, 137 ,FY, -7.54E+04 F, 138 ,FY, -7.87E+04 F, 139 ,FY, -8.06E+04 F, 608 ,FY, 2.85E+04 F, 609 ,FY, 2.86E+04 F, 610 ,FY, 6.60E+04 F, 611 ,FY, 6.25E+04 F, 612 ,FY, 6.16E+04 F, 613 ,FY, 6.11E+04 F, 614 ,FY, 6.10E+04 F, 615 ,FY, 6.08E+04 F, 616 ,FY, 6.08E+04 F, 617 ,FY, 6.08E+04 F, 618 ,FY, 6.08E+04 F, 619 ,FY, 6.08E+04 F, 620 ,FY, 6.08E+04 F, 621 ,FY, 6.08E+04 F, 622 ,FY, 6.08E+04 F, 623 ,FY, 6.08E+04 F, 624 ,FY, 6.08E+04 F, 625 ,FY, 6.08E+04 F, 626 ,FY, 6.08E+04 F, 627 ,FY, 6.09E+04
F, 628 ,FY, 6.10E+04
F, 629 ,FY, 6.12E+04 F, 630 ,FY, 6.16E+04 F, 631 ,FY, 6.26E+04
F, 632 ,FY, 6.61E+04 !右隧道开挖求解 ALLSEL,ALL SOLVE FINISH SAVE
!右隧道开挖模拟分析后处理
/POST1 PLNSOL,U,Y,0,1 PLNSOL,U,X,0,1 PLNSOL,S,Y,0,1 PLNSOL,S,X,0,1 PLNSOL,S,1,0,1 PLNSOL,S,2,0,1 PLNSOL,S,3,0,1 PLNSOL,EPPL,EQV,0,1 PLNSOL,S,EQV,0,1 ESEL,S,LIVE ESEL,R,TYPE,,1 ETABLE,,SMISC,6
ETABLE,,SMISC,12
ETABLE,,SMISC,1 ETABLE,,SMISC,7 ETABLE,,SMISC,2 ETABLE,,SMISC,8 PLLS,SMIS6,SMIS12,1,0 PLLS,SMIS1,SMIS7,1,0
方向的位移 !X 方向的位移 !Y 方向的应力 !X 方向的应力 !第一主应力 !第二主应力 !第三主应力 !塑性区
!6 和 12 表示弯矩!1 和 7 表示轴力
!2 和 8 表示剪力 !绘制弯矩图 !绘制轴力图
!Y PLLS,SMIS2,SMIS8,1,0
注：编程中出现的问题主要如下：
!绘制剪力图
1）在隧道开挖过程中，选取隧道开挖掌子面周边单元，并且求用 NFORCE 求解 单元上节点的节点力，方法及在隧道开挖轮廓周边节点上施加节点力的意义， 不施加对模拟开挖结果的影响。
2）整体模拟结束，为什么在显示衬砌的弯矩、剪力、轴力图示。只显示右边隧 道衬砌上的弯矩图、剪力图、轴力图。而左边隧道的弯矩图、剪力图、轴力 图均不显示。
3）NFORCE 命令的用法，如何求解所需节点的节点力?
4）
因篇幅问题不能全部显示，请点此查看更多更全内容