数值分析中的问题 [复制链接]

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ 'QV4=h`  
gvU6p[D  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 V+Tj[:ok  
LTY.i3  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 d)G-K+&B  
JV/,QWar  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： )T3wU~%  
_N]yI0k(  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ w}1)am&pD  
m:41zoV  
以下是我的计算命令流，希望得到高手的帮助啊！ =0:hrg+Zgx  
{'ZnxK'  

;I7Z*'5!  
Fini 7uDUZdJy  
/clear /.?\P#9)  
/title,TUNNEL ANALYSIS an7N<-?
 
/prep7 RTU:J67E  
ET,1,PLANE2 zV&3l9?U  
keyopt,1,3,2 67g/(4&  
fBF}-{VX(  
*SET,R,0.7 k@Tt,.];  
*SET,H,10.0 UB1/FM4
~  
*SET,X1,-2*R                                                    B 3h<K}  
*SET,Y1,-2*R                                                    + mcN6/  
*SET,W1,4*R                                                   UEeqk"t^  
*SET,H1,4*R                                                   hB.dqv]^  
                                                              s=nds"J  
*SET,X2,-10.0                                                   !^7:Rr_  
*SET,Y2,H-20.0                                                    Kq?7#,_  
*SET,W2,20.0                                                  'IP!)DS  
*SET,H2,20.0 f65Sr"qB3  
CYL4,0,0,0,45,R,405,0 y,x~S\>+  
BLC4,X1,Y1,W1,H1                                              MQ+ek4  
BLC4,X2,Y2,W2,H2 4mAtYm  
AOVL,1,2,3                                              `WB|h)Y  
NUMCMP,ALL                                                    QXz!1o+"  
NUMMRG,ALL,,,,LOW 4(&'V+o  
LANG,11,8,90,, ~4^e a  
LANG,11,7,90,, O92Yd$S  
LANG,10,7,90,, !+6l.`2WI  
LANG,10,6,90,, |l|]Tw  
LANG, 9,6,90,, xpxm9ySwu  
LANG, 9,5,90,, 4{?x(~  
LANG,12,5,90,, Lv4=-mWv&0  
LANG,12,8,90,, 'SieZI
m)  
L,1,7 OS 6 )`  
L,2,8 QD<4(@c5|  
L,3,5 ?*@h]4+k'  
L,4,6 _e3'f:  
NUMCMP,ALL                                                    ?MYD}`Cv  
NUMMRG,ALL,,,,LOW EyeLC6u  
FLST,3,4,4,ORDE,2   uM}O8N  
FITEM,3,29  yU.0'r5uR  
FITEM,3,-32 v(^rq  
ASBL,2,P51X  f+6l0@K2  
FLST,3,8,4,ORDE,8   x\DkS,O  
FITEM,3,14  Ae#6=]V+^  
FITEM,3,16  NIfc/%  
FITEM,3,18  54{"ni2a  
FITEM,3,20  $T4PC5.  
FITEM,3,22  $ZU(bEUOG  
FITEM,3,24  J]0#M:w&  
FITEM,3,26  hvyN8We  
FITEM,3,28  k +Oq$Pi  
ASBL,       3,P51X {dwV-qz  
NUMCMP,ALL                                                    !D3}5A1,  
NUMMRG,ALL,,,,LOW D'h2 DP!  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 3$8}%?i  
MP,EX,1,3.0E6                                    +7^Ul6BB#K  
MP,PRXY,1,0.33                                                07|NPS  
MP,DENS,1,1900.0 4.'EEuRw\}  
TB,DP,1                                                       o^?{j*)g  
TBDATA,1,2.4E4,13.6,0 zEW:Xe)  
;$QJnQ"R  
MP,EX,2,3.0E6                                    m}9V@@  
MP,PRXY,2,0.33                                                }'W^Ki$  
MP,DENS,2,1900.1 aK{\8L3]  
TB,DP,2                                                       >{S$0D  
TBDATA,1,2.4E4,13.6,0 yLx.*I^6  
FQFENq''B  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 pqSE|3*l  
type,1 $D5U
#  
MAT,1  I[UA' ~f  
lsel,S,line,,1,8,1 X$9QW3.M  
lsel,A,line,,15,27,4     cPcV[6)5K9  
Lesize,ALL,,,4,,,,,0 ?/(K7>`  
lsel,S,line,,29,32,1     Xoy1Gi?  
Lesize,ALL,,,3,,,,,0 F4x7;?W{*  
  6iG<"{/U5  
lsel,S,line,,9,13,2 lt4jnV2"a  
lsel,A,line,,18,20,2 ,[\(U!Z7:%  
lsel,A,line,,21 0UJ`<Bfd  
lsel,A,line,,26,28,2     @<G/H|f  
Lesize,ALL,,,10,,,,,0 hD1AK+y  
lsel,S,line,,10 _6!/}Fm  
lsel,A,line,,22 -y?Z}5-rs  
lsel,A,line,,24,25,1 `tG_O  
Lesize,ALL,,,20-H-2*R,,,,,0 >!1f`  
lsel,S,line,,12,14,2 SN(=e#ljE  
lsel,A,line,,16,17,1 p2vBj.*J  
Lesize,ALL,,,H-2*R,,,,,0 a*j <TR  
asel,s,area,,2,13,1  +{%(
_<  
amesh,all u;Eu<jU1  
   x\.i`ukx  
type,1 C|\^uR0  
MAT,2 &8_;:  
asel,s,area,,1  +^?
-}v  
amesh,all 7omGg~!k(  
nummrg,all,,,,low =h::VB}Lv  
numcmp,all &ZN'Ey?  
J&:W4\ m  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 *\KMkx  
/solu GB?#1|,  
antype,static @$d\5Q(G  
deltim,0.1,0.05,0.2 t`T\d\  
autots,on !使用自动时间步 i*R:WTw#  
pred,on !打开时间步长预测器 4~A#^5J  
lnsrch,on !打开线性搜索 -]\E}Ti  
nlgeom,on !打开大位移效果 S7ehk*`  
nropt,full !设定牛顿-拉普森选项 )8 :RiG2B  
cnvtol,f,,0.02,0 C;q}3c*L  
Nsel,s,loc,x,-10 W8$=a
  
Nsel,a,loc,x,10 P0ltN  
D,all,,,,,,UX,,,,, )O@^H  
Nsel,s,loc,y,H-20 ~f10ZB_k>'  
D,all,,,,,,Uy,,,,, BhW]Oq&  
time,1 gAA %x7  
acel,,10  V`,[=u?c  
nsel,all xdf
vme[  
esel,all -`dxx)x  
allsel,all v*~%x  
solve J(0=~Z[  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 {]\uR-a(o  
time,2 rtL
9cw5  
esel,s,type,,2 ./nYXREO|  
ekill,all on;sq8;  
nsel,all *{dD'9Bg  
esel,all ZqbM%(=z(`  
allsel,all 1mn$Rh&dO  
solve xrVZxK:!  
fini n5k^v$'  
s$>m0^  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 1gShV ]2  
!/POST1   i %hn  
!LCDEF,1,1,, B:- KZuO  
!LCDEF,2,2,, /-{O\7-D  
!LCASE,2, wGzXp5 dl  
!LCOPER,SUB,1, ,2, ONN{4&7@<  
!Nsel,s,loc,Y,H ,rFLpQl  
!PLVECT,U, , , ,VECT,ELEM,ON,0    |<3Q+EB^  
!PLNSOL, U,Y, 0,1.0 O&gwr  
"UoCT7X  
*KjVPs  
]u4>;sa  
udEb/7ZL  

看看谢谢楼主

观望中。。。。。

学习学习