地铁盾构隧道掘进过程ANSYS数值模拟分析

地铁盾构隧道掘进过程数值模拟分析
具体做法如下：
⑴采用在掘削面施加顶进压力的方法来模拟开挖面土体的移动；
⑵采用施加注浆压力的方法来模拟盾尾注浆过程；
⑶采用更换注浆层材料参数的方法来模拟盾尾注浆对盾尾空隙的填充效果。

由于目前计算软件的限制，难以模拟盾构机推进过程中对土体的扰动，这里简化处理。

即假定盾构隧道开挖后，随机进行注浆。

计算时，只需将开挖不断地向前推进，同时在后面进行注浆、换材料参数等操作，即可实现盾构隧道的动态开挖过程，详细的计算操作见后面的求解过程。

工程问题的描述
地铁盾构隧道管片衬砌内径为5.4m，外径为D=6m，埋深为12m。

自上至下，根据土层的物理性质参数不同将其分为3层，各层的材料参数和厚度如下：
第一层：厚8m，E=3.94Mpa，v=0.35，ρ=18.28KN/m3
第二层：厚18m，E=20.6Mpa，v=0.3，ρ=20.62KN/m3
第三层：厚15m，E=500Mpa，v=0.33，ρ=21.6KN/m3
施工中掘削面顶进压力为0.3Mpa，盾尾注浆压力为0.15Mpa
模型的建立
!进入前处理器
FINISH
/CLE
/PREP7 !进入前处理器
ET,1,SOLID45 !定义实体单元
ET,2,MESH200,6 !定义非求解单元，辅助面网格的划分!定义模型中的材料参数。

模型中共有5种材料，其中土体有3种，即地表浅层覆土、盾构隧道所在土层和基岩及管片衬砌和注浆层。

其中，管片衬砌为管片式的拼装结构，为了计算方便，将其等效为一均质体，等效时对原有刚度进行折减。

定义材料参数的命令流如下：!土体材料参数
MP,EX,1,3.94E6 !第一层土层材料参数
MP,PRXY,1,0.35
MP,DENS,1,1828
MP,EX,2,20.6E6 !第二层土层材料参数
MP,PRXY,2,0.30
MP,DENS,2,2160
MP,EX,3,500E6 !第三层土层材料参数
MP,PRXY,3,0.33
MP,DENS,3,2160
!管片材料参数，管片衬砌按各向同性计算
MP,EX,4,27.6E9
MP,PRXY,4,0.2
MP,DENS,4,2500
!注浆层，参数按水泥土取值
MP,EX,5,1E9
MP,PRXY,5,0.2
MP,DENS,5,2100
!建立平面内模型并划分单元
!在隧道中心线定义局部坐标，便于后来的实体选取LOCAL,11,0,0,0,0 !设置局部直角坐标系原点坐标(0,0,0) LOCAL,12,1,0,0,0 !设置局部柱坐标系原点坐标(0,0,0) CSYS,11 !将当前坐标系转化为局部直角坐标系WPCSYS,-1 !同时将工作平面恢复指定到局部直角坐标系下
CYL4,,,,,2.7,90 !创建开挖土体所在的面
CYL4,0,0,2.7,0,3,90 !创建管片
CYL4,0,0,3,0,3.2,90 !创建注浆加固圈
RECTNG,0,4.5,0,4.5
AOVLAP,ALL
NUMMRG,ALL !压缩编号
NUMCMP,ALL !合并图素
RECTNG,4.5,31.5,0,4.5
NUMMRG,ALL
NUMCMP,ALL
!划分单元
MSHAPE,0,2D !采用四边形单元划分网格MESHKEY,1 !映射网格划分
TYPE,2
LESIZE,1,,,6
LESIZE,2,,,6
LESIZE,3,,,6
AMESH,1
LESIZE,4,,,6
LESIZE,8,,,2
LESIZE,9,,,2
AMESH,2
LESIZE,5,,,6
LESIZE,10,,,1
LESIZE,11,,,1
AMESH,3
LESIZE,12,,,3
LESIZE,13,,,3
LESIZE,6,,,3
LESIZE,7,,,3
LESIZE,14,,,8,2
LESIZE,16,,,8,0.5
AMAP,4,7,6,8,10
AMAP,5,9,8,11,12
!利用对称性，得到下半部分模型
ARSYM,Y,ALL !通过坐标轴对称建立面NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
!建立隧道下方土层模型
RECTNG,0,4.5,-4.5,-26
RECTNG,4.5,31.5,-4.5,-26
NUMMRG,ALL
NUMCMP,ALL
LESIZE,28,,,3
LESIZE,29,,,5,0.5
LESIZE,30,,,5,2
LESIZE,32,,,5,0.5
LESIZE,31,,,8,2
ASEL,S,AREA,,11,12,1
AMESH,ALL
!建立隧道土方土层模型
RECTNG,0,4.5,4.5,15
RECTNG,4.5,31.5,4.5,15
NUMMRG,ALL
NUMCMP,ALL
LESIZE,34,,,3
LSEL,S,LINE,,35,36,1
LSEL,A,LINE,,33
LESIZE,ALL,,,4
LSEL,S,LINE,,37
LESIZE,37,,,8,0.5
AMESH,13
AMESH,14
NUMMRG,ALL
NUMCMP,ALL
!利用对称性得到平面内的全部模型
ALLSEL,ALL
ARSYM,X,ALL
NUMMRG,ALL
NUMCMP,ALL
!建立体模型。

将划分好的平面模型通过沿隧道轴线进行拉伸，得到隧道及其所在地层的三维实体模型
!沿隧道轴线定义一系列关键点，点间距为3m
*DO,I,1,20
K,,0,0,-3*I
*ENDDO
!沿着这些点建立线，便于后来体的拉伸
L,3,38
*DO,I,1,19
L,37+I,38+I
*ENDDO
*DO,I,1,20
LESIZE,64+I,,,1
*ENDDO
!拉伸，拉伸后建立完整的模型
TYPE,1
MAT,1
ESYS,11 !定义单元坐标系EXTOPT,ACLEAR,1
CSYS,12 !激活局部柱坐标系
*DO,I,1,20
ASEL,S,LOC,Z,-3*(I-1)
VDRAG,ALL,,,,,,64+I
*ENDDO
ALLSEL,ALL
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
!分层赋予土体材料参数。

拉伸成体模型时为了方便，统一将所有的材料参数号都设置为1，而该实例中考虑了土层特性的分层，如前面建立模型是将土层设置为3种类型的土体，因而需将不同地层赋予不同性质类型的材料参数
!赋予隧道所在地层的材料参数
CSYS,0 !激活整体坐标系
NSEL,S,LOC,Y,-10.5,5 !选择二号地层内所在节点
ESLN,S
MPCHG,2,ALL
!赋予基岩材料参数
NSEL,S,LOC,Y,-11,-26
ESLN,S
MPCHG,3,ALL
!加载与求解
!施加边界条件，并进行求解设置
FINISH
/SOLU
CSYS,0 !激活直角坐标系
/VIEW,1,-1,0.5,1
/REPLOT
EPLOT
NLGEOM,ON
NROPT,FULL
PRED,ON,,ON
OUTRES,ALL,LAST !控制写入数据库和结果文件的结果数据NSEL,S,LOC,X,-31.5
D,ALL,UX
NSEL,S,LOC,X,31.5
D,ALL,UX
NSEL,S,LOC,Y,-26
D,ALL,UY
NSEL,S,LOC,Z,0
NSEL,A,LOC,Z,-60
D,ALL,UZ
ALLSEL,ALL
EPLOT
/REPLOT
ACEL,0,10,0
STEP=0
!求解自重应力场
TIME,1 !时间步（荷载步结束的时间）
NSUBST,10 !子步数和时间步长
SOLVE
FINISH
!为了养成良好的习惯,自重应力场求解结束后,需对结果进行初步检查,判断准确无误后在进行下面的计算。

进入后处理器，分别绘制Y方向应力云图和Y方向位移云图。

/POST1
PLNSOL,S,Y,0,1
PLNSOL,U,Y,0,1
图0.1 建立的模型图0.2 附好材料属性的单元
图0.3初应力求解完成后Y方向的位移图0.4初应力求解完成后Y方向的应力!第一步开挖，进尺3m。

在进行开挖计算时首先区分两个概念，一个是开挖步，另一个是计算步。

一个开挖步也可以认为是一个循环进尺。

在计算实例中，第1-12步开挖，每步掘进3m，第13步余下24m。

计算中一个开挖步包含多个计算步：第一个计算步模拟开挖、施加盾尾注浆压力；第二个计算步模拟管片衬砌和注浆硬化。

!第一个计算步的命令流如下：
!*****第一步开挖，第一个计算步******
FINISH
/SOLU
RESCONTROL,FILE_SUMMARY !列表显示当前文件中的荷载步和荷载子步的信息
ANTYPE,,REST!重新启动介绍
TIME,2
NSUBST,10,1000,5
!DELTIM,0.1,0.001,0.2
CSYS,12 !激活局部柱坐标系
WPCSYS,-1 !将工作平面转移到当前局部坐标系ALLSEL,ALL
EPLOT
/REPLOT
VSEL,S,LOC,X,0,3.2
VSEL,R,LOC,Z,0,-3
VPOLT
/REPLOT
ESLV,S
EKILL,ALL
ASEL,S,LOC,X,0,3.2
ASEL, R, LOC, Z,-3
SFA, ALL, 1, PRES, 0.3E6 !施加掌子面顶进压力
ASEL,S,LOC,X,3.2
ASEL, R, LOC, Z, 0,-2.9
SFA, ALL, 1, PRES, 0.15E6 !施加注浆压力
ALLSEL, ALL
SOLVE
FINISH
!求解完成后，继续进行第二计算步的计算，对于隧道开挖模拟的重启动，需要注意之前第一步计算时杀死的单元，以及施加的顶进荷载和注浆压力在第二步计算时需要重新施加；单元的生死命令需要重新定义。

命令流如下：
/SOLU
RESCONTROL,FILE_SUMMARY
ANTYPE,,REST
TIME,3
NSUBST,10,1000,5
!DELTIM,0.1,0.001,0.2,1
CSYS,12
WPCSYS,-1
VSEL,S,LOC,X,0,2 !选择开挖的核心土
VSEL,R,LOC,Z,0,-3
ESLV,S
EKILL,ALL !杀死核心土
VSEL,S,LOC,X,2,2.3
VSEL,R,LOC,Z,0,-3
ESLV,S
MPCHG,4,ALL !施加管片
VSEL,S,LOC,X,2.3,2.7
VSEL,R,LOC,Z,0,-3
ESLV,S
MPCHG,5,ALL !施做注浆层
ASEL,S,LOC,X,0,3.19
ASEL,R,LOC,Z,-3
SFA,ALL,1,PRES,0.3E6
ASEL,S,LOC,X,3.2
ASEL,R,LOC,Z,0,-2.9
SFA,ALL,1,PRES,0
ALLSEL,ALL
SOLVE
FINISH
!进行循环开挖。

由于隧道开挖计算采用的是重启动分析，因此所有的参数和变量将重新回到初始值，因而不能采用*DO循环的方式来做，而由前面第一步开挖的求解可以看出，许多代码可以重用，将这些代码编成MAC文件，可以提高代码的简洁性。

!将第一步开挖，第一计算步中的命令流改写成以STEP（开挖步）为参数的命令流形式，并将其文件命名为EXCA1.MAC,代码如下：!******************EXCA1.MAC***************
!*CREATE,’EXCA1’,’MAC’,’’
!TIME,2*STEP !指定时间步
!RESCONTROL,FILE_SUMMARY
!NSUBST,10,1000,5 !设置荷载步,荷载子步,最大荷载步,
最小荷载步
!CSYS,12 !激活局部柱坐标系
!WPCSYS,-1 !将工作平面转换到局部坐标系!VSEL,S,LOC,X,0,2 !选择开挖核心土体
!VSEL,R,LOC,Z,0,-3*(STEP -1) !选择在0到-3范围内的核心土!ESLV,S !选择与体相关的单元
!EKILL,ALL !杀死核心土体单元
!VSEL,S,LOC,X,2,2.3 !选择盾构管片
!VSEL,R,LOC,Z,0,-3*(STEP-1) !选择在0到-3范围内的管片!ESLV,S !选择与体相关的单元
!MPCHG,4,ALL !将材料属性改为管片单元的材料属性!VSEL,S,LOC,X,2.3,2.7 !选择注浆加固区的体
!VSEL,R,LOC,Z,0,-3*(STEP-1) !选择0到-3范围注浆加固区体!ESLV,S !选择与体相关的单元
!MPCHG,5,ALL !将材料属性改为注浆加固区材料属性
!VSEL,S,LOC,X,0,3.2 !选择体
!VSEL,R,LOC,Z,-3*(STEP-1),-3*STEP !选择位置在-3到-6范围内的体
!ESLV,S !选择与体相关的单元
!EKILL,ALL !杀死单元，向前开挖土体
!ASEL,S,LOC,X,0,3.19 !选择面
!ASEL,S,LOC,Z,-3*(STEP-1)
!SFA,ALL,1,PRES,0 !将前一步计算时施加的注浆压力取消
!ASEL,S,LOC,X,0,3.19
!ASEL,R,LOC,Z,-3*STEP
!SFA,ALL,1,PRES,0.3E6 !施加盾构机顶进压力
!ASEL,S,LOC,X,3.2
!ASEL,R,LOC,Z,-3*(STEP-1)-0.1,-3*STEP+0.1
!SFA,ALL,1,PRES,0.15E6 !施加注浆压力
!ALLSEL,ALL
!SOLVE
!FINISH
!*END
!同样将第一步开挖，第二计算步中的命令流改写为以STEP（开挖步）为参数的命令流形式，并将其文件命名为EXCA2.MAC,代码如下
!*******************EXCA2.MAC************
!*CREATE,’EXCA2’,’MAC’,’’
!TIME,2*STEP+1
!RESCONTROL,FILE_SUMMARY
!NSUBST,10,1000,5 !设置荷载步,荷载子步,最大荷载步,最小荷载步
!CSYS,12
!WPCSYS,-1
!VSEL,S,LOC,X,0,2
!VSEL,R,LOC,Z,0,-3*STEP
!ESLV,S
!EKILL,ALL
!VSEL,S,LOC,X,2,2.3
!VSEL,R,LOC,Z,0,-3*STEP
!ESLV,S
!MPCHG,4,ALL
!VSEL,S,LOC,X,2.3,2.7
!VSEL,R,LOC,Z,0,-3*STEP
!ESLV,S
!MPCHG,5,ALL
!ASEL,S,LOC,X,0,3.19
!ASEL,R,LOC,Z,-3*STEP
!SFA,ALL,1,PRES,0.3E6
!ASEL,S,LOC,X,3.2
!ASEL,R,LOC,Z,-3*(STEP-1)-0.1,-3*STEP+1 !SFA,ALL,1,PRES,0
!ALLSEL,ALL
!SOLVE
!FINISH
!*END
!************第2步开挖*************
/SOLU
ANTYPE,,REST
STEP=2
EXCA1.MAC !调用名为EXCA1的宏文件
/SOLU
ANTYPE,,REST
STEP=2
EXCA2.MAC !调用名为EXCA2的宏文件
!******************第3步开挖******************** /SOLU
ANTYPE,,REST
STEP=3
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=3
EXCA2.MAC
!*******************第4步开挖***************
/SOLU
ANTYPE,,REST
STEP=4
/SOLU
ANTYPE,,REST
STEP=4
EXCA2.MAC
!******************第5步开挖******************** /SOLU
ANTYPE,,REST
STEP=5
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=5
EXCA2.MAC
!***************第6步开挖********************* /SOLU
ANTYPE,,REST
STEP=6
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=6
!*********************第7步开挖******************** /SOLU
ANTYPE,,REST
STEP=7
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=7
EXCA2.MAC
!******************第8步开挖********************
/SOLU
ANTYPE,,REST
STEP=8
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=8
EXCA2.MAC
!******************第9步开挖********************
/SOLU
ANTYPE,,REST
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=9
EXCA2.MAC
!******************第10步开挖******************** /SOLU
ANTYPE,,REST
STEP=10
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=10
EXCA2.MAC
!******************第11步开挖******************** /SOLU
ANTYPE,,REST
STEP=11
EXCA1.MAC
/SOLU
ANTYPE,,REST
EXCA2.MAC
!******************第12步开挖********************
/SOLU
ANTYPE,,REST
STEP=12
EXCA1.MAC
/SOLU
ANTYPE,,REST
STEP=12
EXCA2.MAC
!**********第13步开挖，开挖余下来的24m**************** /SOLU
RESCONTROL,FILE_SUMMARY
ANTYPE,,REST
STEP=13
TIME,2*STEP
NSUBST,10,1000,5
CSYS,12
WPCSYS,-1
VSEL,S,LOC,X,0,2
VSEL,R,LOC,Z,0,-3*(STEP-1)
EKILL,ALL
VSEL,S,LOC,X,2,2.3
VSEL,R,LOC,Z,0,-3*(STEP-1) ESLV,S
MPCHG,4,ALL
VSEL,S,LOC,X,2.3,2.7
VSEL,R,LOC,Z,0,-3*(STEP-1) ESLV,S
MPCHG,5,ALL
VSEL,S,LOC,X,0,3.2
VSEL,R,LOC,Z,-3*(STEP-1),-60 ESLV,S
EKILL,ALL
ASEL,S,LOC,X,0,3.19
ASEL,R,LOC,Z,-3*(STEP-1) SFA,ALL,1,PRES,0
ASEL,S,LOC,X,3.2
ASEL,R,LOC,Z,-3*(STEP-1)-0.1,-60 SFA,ALL,1,PRES,0.15E6 ALLSEL,ALL
SOLVE
!************开挖的最后一个计算步************** /SOLU
RESCONTROL,FILE_SUMMARY
ANTYPE,,REST
STEP=13
TIME,2*STEP+1
NSUBST,10,1000,5
CSYS,12
WPCSYS,-1
VSEL,S,LOC,X,0,2
VSEL,R,LOC,Z,0,-60
ESLV,S
EKILL,ALL
VSEL,S,LOC,X,2,2.3
VSEL,R,LOC,Z,0,-60
ESLV,S
MPCHG,4,ALL
VSEL,S,LOC,X,2.3,2.7
VSEL,R,LOC,Z,0,-60
ESLV,S
MPCHG,5,ALL
ASEL,S,LOC,X,3.2
ASEL,R,LOC,Z,-3*(STEP-1)-0.1,-60
SFA,ALL,1,PRES,0
ALLSEL,ALL
SOLVE
FINISH
!计算结果分析
!1.地层位移分析
本实例中分析地层位移仍然是相对于自重固结下的位移，因而分析地层位移时需采用命令流编写宏文件，对所关心的计算步进行工况组合运算，才能获得相对位移
!***************LOADCASE.MAC****************
!LCDEF,1,1 !将第一计算步结果定义为荷载工况1
!LCDEF,2,3 !将第N计算步的结果定义为荷载工况N
!LCASE,N !读入荷载工况2
!LCOPER,SUB,1 !将荷载工况2减去荷载工况1
!**********第1步开挖引起的地层位移*************
/POST1
N=3
LOADCASE.MAC
FINISH
!**********第2步开挖引起的地层位移*************
/POST1
N=5
LOADCASE.MAC
FINISH
!**********第3步开挖引起的地层位移************* /POST1
N=7
LOADCASE.MAC
FINISH
!**********第4步开挖引起的地层位移************* /POST1
N=9
LOADCASE.MAC
FINISH
!**********第5步开挖引起的地层位移************* /POST1
N=11
LOADCASE.MAC
FINISH
!**********第6步开挖引起的地层位移************* /POST1
N=13
LOADCASE.MAC
FINISH
!**********第7步开挖引起的地层位移************* /POST1
N=15
LOADCASE.MAC
FINISH
!**********第8步开挖引起的地层位移************* /POST1
N=17
LOADCASE.MAC
FINISH
!**********第9步开挖引起的地层位移************* /POST1
N=19
LOADCASE.MAC
FINISH
!**********第10步开挖引起的地层位移************* /POST1
N=21
LOADCASE.MAC
FINISH
!**********第11步开挖引起的地层位移************* /POST1
N=23
LOADCASE.MAC
FINISH
!**********第12步开挖引起的地层位移************* /POST1
N=25
LOADCASE.MAC
FINISH
!**********第13步开挖引起的地层位移************* /POST1
N=27
LOADCASE.MAC
FINISH
图1.1 第1步开挖引起的地层位移（m ） 图1.2 第2步开挖引起的地层位移（m ）
图1.3 第3步开挖引起的地层位移（m ） 图1.4 第4步开挖引起的地层位移（m ）
图1.5 第5步开挖引起的地层位移（m ） 图1.6 第6步开挖引起的地层位移（m ）
图1.7 第7步开挖引起的地层位移（m ） 图1.8 第8步开挖引起的地层位移（m ）
图1.9 第9步开挖引起的地层位移（m ） 图1.10 第10步开挖引起的地层位移（m ）
图1.11 第11步开挖引起的地层位移（m ） 图1.12 第12步开挖引起的地层位移（m ）
图1.13 第13步开挖引起的地层位移（m）综合分析地层位移云图，可以发现地层位移的最大值出现在隧道的拱顶和仰拱处，其中拱顶产生最大下沉约4.5cm，仰拱处产生最大隆起约4.9cm。

!2.地表沉降分析，利用宏文件LOADCASE,MAC文件，求解地表沉降
!*********第1次开挖引起的地表沉降*****************
FINISH
/POST1
N=3
LOADCASE.MAC
CSYS,0
WPCSYS,-1
NSEL,S,LOC,Y,15 !选择地表节点
ESLN,S !选择附在节点上的单元
EPLOT
PLNSOL,U,Y,0,1 !显示地表的Y方向的位移
!*********第2次开挖引起的地表沉降***************** FINISH
/POST1
N=5
LOADCASE.MAC
CSYS,0
WPCSYS,-1
NSEL,S,LOC,Y,15
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第3次开挖引起的地表沉降***************** FINISH
/POST1
N=7
LOADCASE.MAC
CSYS,0
WPCSYS,-1
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第4次开挖引起的地表沉降***************** FINISH
/POST1
N=9
LOADCASE.MAC
CSYS,0
WPCSYS,-1
NSEL,S,LOC,Y,15
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第5次开挖引起的地表沉降***************** FINISH
/POST1
N=11
LOADCASE.MAC
CSYS,0
WPCSYS,-1
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第6次开挖引起的地表沉降***************** FINISH
/POST1
N=13
LOADCASE.MAC
CSYS,0
WPCSYS,-1
NSEL,S,LOC,Y,15
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第7次开挖引起的地表沉降***************** FINISH
/POST1
N=15
LOADCASE.MAC
CSYS,0
WPCSYS,-1
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第8次开挖引起的地表沉降***************** FINISH
/POST1
N=17
LOADCASE.MAC
CSYS,0
WPCSYS,-1
NSEL,S,LOC,Y,15
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第9次开挖引起的地表沉降***************** FINISH
/POST1
N=19
LOADCASE.MAC
CSYS,0
WPCSYS,-1
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第10次开挖引起的地表沉降***************** FINISH
/POST1
N=21
LOADCASE.MAC
CSYS,0
WPCSYS,-1
NSEL,S,LOC,Y,15
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第11次开挖引起的地表沉降***************** FINISH
/POST1
N=23
LOADCASE.MAC
CSYS,0
WPCSYS,-1
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第12次开挖引起的地表沉降***************** FINISH
/POST1
N=25
LOADCASE.MAC
CSYS,0
WPCSYS,-1
NSEL,S,LOC,Y,15
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
!*********第13次开挖引起的地表沉降***************** FINISH
/POST1
N=27
LOADCASE.MAC
CSYS,0
WPCSYS,-1
ESLN,S
EPLOT
PLNSOL,U,Y,0,1
图2.1第1次开挖引起的地表沉降（m）图2.2 第2次开挖引起的地表沉降（m）
图2.3 第3次开挖引起的地表沉降（m ） 图2.4 第4次开挖引起的地表沉降（m ）
图2.5 第5次开挖引起的地表沉降（m ） 图2.6 第6次开挖引起的地表沉降（m ）
图2.7 第7次开挖引起的地表沉降（m ） 图2.8 第8次开挖引起的地表沉降（m ）
图2.9 第9次开挖引起的地表沉降（m ） 图2.10 第10次开挖引起的地表沉降（m ）
图2.11 第11次开挖引起的地表沉降（m ） 图2.12 第12次开挖引起的地表沉降（m ）
图2.13 第13次开挖引起的地表沉降 （m ）
在盾构机的顶进作用下，地表发生了前隆后沉的现象。

随着盾构
机的推进，开挖面后方的沉降槽越来越宽，最大沉降量也越来越大，隆起区域也逐步向前推进。

地表最大沉降量约为12mm，最大隆起量约为3mm。

!3.管片结构受力分析
进行管片受力分析，不需要进行工况组合，直接将每步计算结果读入即可。

!*******建立管片求解的宏文件GUANPIAN.MAC********** !CSYS,12
!WPCSYS,-1
!VSEL,S,LOC,X,2,2.3
!VSEL,R,LOC,Z,0,-3*STEP !选择管片
!ESLV,S !选择管片单元
!SET,2*STEP+1,LAST,1 !从结果文件中读出指定的数据组
!PLNSOL,S,EQV,0,1 !显示管片上的等效应力
!*********第1次开挖后管片上的等效应力************
FINISH
/POST1
STEP=1
GUANPIAN.MAC
!*********第2次开挖后管片上的等效应力************
FINISH
/POST1
STEP=2
GUANPIAN.MAC
!*********第3次开挖后管片上的等效应力************ FINISH
/POST1
STEP=3
GUANPIAN.MAC
!*********第4次开挖后管片上的等效应力************ FINISH
/POST1
STEP=4
GUANPIAN.MAC
!*********第5次开挖后管片上的等效应力************ FINISH
/POST1
STEP=5
GUANPIAN.MAC
!*********第6次开挖后管片上的等效应力************ FINISH
/POST1
STEP=6
GUANPIAN.MAC
!*********第7次开挖后管片上的等效应力************ FINISH
/POST1
STEP=7
GUANPIAN.MAC
!*********第8次开挖后管片上的等效应力************ FINISH
/POST1
STEP=8
GUANPIAN.MAC
!*********第9次开挖后管片上的等效应力************ FINISH
/POST1
STEP=9
GUANPIAN.MAC
!*********第10次开挖后管片上的等效应力************ FINISH
/POST1
STEP=10
GUANPIAN.MAC
!*********第11次开挖后管片上的等效应力************
FINISH
/POST1
STEP=11
GUANPIAN.MAC
!*********第12次开挖后管片上的等效应力************ FINISH
/POST1
STEP=12
GUANPIAN.MAC
!*********第13次开挖后管片上的等效应力************ FINISH
/POST1
STEP=13
CSYS,12
WPCSYS,-1
VSEL,S,LOC,X,2,2.3
VSEL,R,LOC,Z,0,-60 !选择管片
ESLV,S !选择管片单元
SET,2*STEP+1,LAST,1 !从结果文件中读出指定的数据组PLNSOL,S,EQV,0,1 !显示管片上的等效应力
图3.1 第1次开挖后管片等效应力（Pa ） 图3.2 第2开挖次后管片等效应力（Pa ）
图3.3 第3次开挖后后管片等效应力（Pa ） 图3.4 第4次开挖后后管片等效应力（Pa ）
图3.5 第5次开挖后后管片等效应力（Pa ） 图3.6 第6次开挖后后管片等效应力（Pa ）
图3.7 第7次开挖后后管片等效应力（Pa ） 图3.8 第8次开挖后后管片等效应力（Pa ）
图3.9 第9次开挖后后管片等效应力（Pa ） 图3.10 第10次开挖后后管片等效应力（Pa ）
图3.11第11次开挖后后管片等效应力（Pa ） 图3.12第12次开挖后后管片等效应力（Pa ）。