关于ansys模拟隧道开挖时，杀死的结构单元激活后形状的讨论 [复制链接]

看了好几个论坛，也没找到这个问题的直接讨论，我在做论文时发现：在模拟开挖过程中，激活初始自重时杀死的结构单元，开挖后发现拱顶位移反而减小了，我就猜想：激活的单元其实是有使围岩的变形退回到初始值的趋势，我弄了个简单的例子试了一下,下面是命令流： wA1Ey:q  
主要思路就是在求自重时杀死结构单元，自重求解后，激活结构单元，再进行计算，发现果真回弹了很多 sX!3_'-  
finish HoL~j({  
/clear X&.:H~xS+  
/PREP7 H7}g!n?  
*afun,deg Q-3r}jJe  
!定义单元类型 >~^`5a`$uI  
ET,1,BEAM3 !定义衬砌支护单元 ~f .y:Sbb  
KEYOPT,1,6,1 !局部坐标系输出力、弯矩 iel-<(~  
ET,2,PLANE42 !定义围岩单元 FvAbh]/4  
KEYOPT,2,3,2 ! Plane strain YQx?* gZS  
ET,3,LINK1 !定义LINK1单元 s!aO*\[<h  
1]Lhk?4t  
MP,EX,1,2.95E10 {dBB{.hX  
MP,PRXY,1,0.2 BPh".RJ  
MP,DENS,1,2500 ^8Z@^
M&O"  
TB,DP,1 !采用D-P模型 '9"%@AFxZ  
TBDATA,1,2.42E6,54 ]2PQ X4t0  
!围岩材料 {=qEBbM  
MP,EX,2,0.8E9 eX@v7i,}  
MP,PRXY,2,0.38 [bsXF#  
MP,DENS,2,2000 ot0U-G(  
TB,DP,2 wePI*."]  
TBDATA,1,0.15E6,23 ovbEmb  
R,1,0.28,0.28*0.28*0.28/12,0.28, @wMQC\Z  
K,1,,,, !创建关键点 +\srZ<67  
K,2,10,, @Jm.HST#S8  
K,3,10,10,, 3jXR"@Z-  
K,4,0,10,, {x9j_/R  
a,1,2,3,4 !. p  
Xout:dn  
esize,1 'H!V54 \j  
mat,2 !赋予围岩单元属性 _tA7=*@8  
type,2 TqXge{r  
mshkey,0 !设定自由网格划分 %6N)G!P  
mshape,0 !设定四边形网格划分 {wHvE4F2  
allsel G~S))p  
asel,s,,,1 !选择面11 a^(2q{*  
amesh,all brj[c>
ID  
n 3h^VQ*]G  
mat,1 !赋予围岩单元属性 aj?2jU~Pq  
type,1 6N&|2:U  
e,27,85 8<Xq=*J+  
e,85,84 ovB
=Zm  
e,84,83 ry
kj2/O  
e,83,82 5 >'6
6gZ  
e,82,81 8-A:k E  
e,81,80 ]I8]mUiUH  
e,80,79  w"BIv9N  
e,79,78 NtqFnxm
/  
e,78,77 t@6w$5:}  
e,77,7 &jt02+Hj'  
*.:!Ax  
/solu ygMd$0:MN  
!施加约束 1y 1_6TZ+  
csys,0 }\>+H  
nsel,s,loc,x,0 !选择X=-52线上所有节点 Q7L)f71i  
nsel,a,loc,x,10 !选择X=52线上所有节点 H<$pHyxU  
d,all,ux !对所选择节点约束X方向位移 ^]i" H|(x  
allsel !ZX&r{pJp  
nsel,s,loc,y,0 !选择Y=--65线上所有节点 )U]:
9)  
d,all,uy eVrnVPkM  
acel,,9.8 Etw~*  
antype,static !设定为静力求解 )=y.^@UT@  
nsubst,100 !设定最大子步数为100 5`{=`  
pred,on !打开时间步长预测器
k~{Fnkt  
nropt,full !设定牛顿-拉普森选项 r1+c/;TpZ  
nlgeom,on !打开大位移效果 >n1h^AW  
lnsrch,on !打开线性搜索 `t%|.=R  
outres,all,all ed=n``P~}  
allsel e
~3]/BL  
esel,s,mat,,1 !选择材料号为1、4的单元 IeH^Wm&^  
ekill,all @`5QG2  
allsel `|&\e_"DE  
solve KM5jl9Vv  
s:3aRQ%  
esel,s,mat,,1 !选择材料号为1的单元 'frWu6]< 4  
ealive,all g%ZdIKj!  
allsel q?(A!1(u  
solve 7:z>+AM[r  
}M^_Z#|,  
' 4,y  
自重计算后的竖向位移图（杀死梁单元） xUQdVrFU  
激活梁单元计算后竖向位移图（是在自重求解后） t?KUK>>w  
图见附件

那么在用单元生死去模拟隧道开挖过程时，先进行初始应力求解，再开挖，在开挖面上施加反向节点荷载，应力释放40%，计算，再激活结构单元后，应力全部释放，结果是全部释放的位移值比前者还要小，也就模拟不了隧道初期支护后，围岩的继续变形，而好多人的硕士论文就是这么做的，不知道该怎么着，希望大家对讨论，多给建议!

刚看了一篇硕士的论文，看到了答案，困扰好几天的问题解决了，解决的方案是，在激活单元后，打开大位移选项，一定要在激活后再打开大位移选项，下面是计算的结果.

这种情况，我在模拟隧道开挖也是这样。期待高手指点。

嗯，老大问我这个上浮是什么个状况，我都不知道怎么答

解决的方案是，在激活单元后，打开大位移选项，一定要在激活后再打开大位移选项,这样被激活的单元就会调整形状来适应当前的变形，这个在一篇硕士论文和ANSYS中文帮助中都有，我当时迷惑了好久，发了帖子那天就找到答案，呵呵

不错，值得研究一下

值得研究一下

生死单元一直是我比较难懂的，学习下

看起来好复杂哦，呵呵。

楼主，请问对于多施工步来说如何进行呢？我的qq1003672678，加我请教下您

又学到一招，谢谢