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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ :==UDVP  
Lx|w~+k}  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 &Z#Vw.7U  
r%OrH-T  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 r>n8`W  
hg)!m\g  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： `K1PGibV  
_Eet2;9  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ e!O &~#'h}  
9 ayH:;  
以下是我的计算命令流，希望得到高手的帮助啊！ O:5ldI  
pZNlcB[Qn-  
)AZ`R8-A  
Fini oZ|{J  
/clear &q kl*#]  
/title,TUNNEL ANALYSIS dA3`b*nC  
/prep7 iX
&Z  
ET,1,PLANE2 Br?++\  
keyopt,1,3,2 ZVC
v(J  
B)LXxdkOn  
*SET,R,0.7 *GY,h$Ul  
*SET,H,10.0 y"{UNM|R  
*SET,X1,-2*R                                                    dW] Ej"W  
*SET,Y1,-2*R                                                    9 u6 g  
*SET,W1,4*R                                                   C 6:pY-  
*SET,H1,4*R                                                   ^+YGSg7  
                                                              >xk:pL*o`  
*SET,X2,-10.0                                                   `qQQQ.K7)z  
*SET,Y2,H-20.0                                                    2g`uC}  
*SET,W2,20.0                                                  Fp* &os  
*SET,H2,20.0 la6e`  
CYL4,0,0,0,45,R,405,0 aG8;,H=%,  
BLC4,X1,Y1,W1,H1                                              2n\i0?RD  
BLC4,X2,Y2,W2,H2 O>{t}6o  
AOVL,1,2,3                                              D0"yZp}  
NUMCMP,ALL                                                    JN9HT0  
NUMMRG,ALL,,,,LOW B@:XC&R^  
LANG,11,8,90,, wZ#~+ }T  
LANG,11,7,90,, `^@g2c+d  
LANG,10,7,90,, QyBK*uNdV  
LANG,10,6,90,,
$(!D/bvJ  
LANG, 9,6,90,, pNk,jeo
  
LANG, 9,5,90,, _16&K}<  
LANG,12,5,90,, 9fk\Ay1P  
LANG,12,8,90,, .,(uoK{  
L,1,7 kgib$t_7  
L,2,8 `XRb:d^  
L,3,5
7cQHRM+1  
L,4,6 _a:!U^4  
NUMCMP,ALL                                                    :D)&>{?  
NUMMRG,ALL,,,,LOW ocuNrkZ  
FLST,3,4,4,ORDE,2   >H]|A<9u(  
FITEM,3,29  gEe W1:AB  
FITEM,3,-32 X_Of k  
ASBL,2,P51X  {e!uvz,e  
FLST,3,8,4,ORDE,8   D=Yag!1  
FITEM,3,14  ~N+/ZVo&y  
FITEM,3,16  2&G1Q'!  
FITEM,3,18  uR[PKLh  
FITEM,3,20  #h|,GvmF<b  
FITEM,3,22  71,0v`Z<  
FITEM,3,24  mU0r"\**c3  
FITEM,3,26  ` TVcI\W  
FITEM,3,28  kkuQ"^<J  
ASBL,       3,P51X &B>uPZ]
 
NUMCMP,ALL                                                    [n@!=T  
NUMMRG,ALL,,,,LOW =Z$=-\<x0.  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 Eo3Aak o  
MP,EX,1,3.0E6                                    wh[:wE]eX  
MP,PRXY,1,0.33                                                F<'l'AsC-  
MP,DENS,1,1900.0 HZ`G)1&)  
TB,DP,1                                                       @R Yb-d  
TBDATA,1,2.4E4,13.6,0 c
bHn\m)J,  
;),O*Z|"v  
MP,EX,2,3.0E6                                    0jx~_zq-j  
MP,PRXY,2,0.33                                                OrqJo!FEg{  
MP,DENS,2,1900.1 8f`b=r(a>  
TB,DP,2                                                       %l$&_xV
-  
TBDATA,1,2.4E4,13.6,0 , $F0D  
ucm.~1G(  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 ff+9(P>*  
type,1 R S_lQ{'  
MAT,1  $5p'+bE  
lsel,S,line,,1,8,1 u:>3j,Cs  
lsel,A,line,,15,27,4     Ydd>A\v\;  
Lesize,ALL,,,4,,,,,0 -W"0,.Dvg  
lsel,S,line,,29,32,1     V<R+A*gY:  
Lesize,ALL,,,3,,,,,0 XQoT},C  
  UK9MWC5g9  
lsel,S,line,,9,13,2 d+| !6  
lsel,A,line,,18,20,2 f0`rJ?us  
lsel,A,line,,21 }(7QJk5 j  
lsel,A,line,,26,28,2     R~?;KJ  
Lesize,ALL,,,10,,,,,0 o_^d>Klb8  
lsel,S,line,,10 'f<_SKd  
lsel,A,line,,22 k=/|?%  
lsel,A,line,,24,25,1 %'g-%2C?  
Lesize,ALL,,,20-H-2*R,,,,,0 wjGjVTtHs  
lsel,S,line,,12,14,2 <$Kv^Y*  
lsel,A,line,,16,17,1 h[)aRo  
Lesize,ALL,,,H-2*R,,,,,0 &!EYT0=>p  
asel,s,area,,2,13,1  Nk~dfY<s  
amesh,all K@u."e
aD  
   |ssIUJ  
type,1 QZq9$;>dW  
MAT,2 v\tbf  
asel,s,area,,1  uNGxz*e  
amesh,all WHk/mAI-s  
nummrg,all,,,,low ^%/5-0?xE  
numcmp,all FwzA_ nn  
&1C9K>  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 ?cxK~Y\  
/solu N5_.m(:  
antype,static l,2z5p  
deltim,0.1,0.05,0.2 L/+J|_J)  
autots,on !使用自动时间步 g#2X'%&+  
pred,on !打开时间步长预测器 4E(5Ccb  
lnsrch,on !打开线性搜索 -"tgEC\tD  
nlgeom,on !打开大位移效果 NB#*`|qt  
nropt,full !设定牛顿-拉普森选项 hd BC ^n  
cnvtol,f,,0.02,0 aw~EK0yU   
Nsel,s,loc,x,-10 ~F6gF7]z  
Nsel,a,loc,x,10 ?
B!ZqJ#  
D,all,,,,,,UX,,,,, ^W05Z!}  
Nsel,s,loc,y,H-20 UN4)>\Y  
D,all,,,,,,Uy,,,,, `*!>79_2C  
time,1 YGmdiY:;1  
acel,,10  qiryC7.E  
nsel,all |6Z MxY  
esel,all 
w1G.^  
allsel,all YfU#kvE'  
solve ?YykCJJ ~@  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 0qUap*fvC  
time,2 ABDUp:  
esel,s,type,,2 6EWB3.x19  
ekill,all A>2p/iMc  
nsel,all OYy !4Fp  
esel,all Czxrn2p/  
allsel,all JB_`lefW,'  
solve E\N=p&g$  
fini EZfa0jJD  
4^^rOi0  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 [LDzR7vnf  
!/POST1   lM%fgyX  
!LCDEF,1,1,, 291|KG  
!LCDEF,2,2,, v('d H"Y  
!LCASE,2, Y/(-mcR  
!LCOPER,SUB,1, ,2, :vT%5CQ  
!Nsel,s,loc,Y,H qE:DJy<  
!PLVECT,U, , , ,VECT,ELEM,ON,0    Vgk,+l!4  
!PLNSOL, U,Y, 0,1.0 I,3!uogn  
%.U{):lNx  
m3-J0D<  
]<LU NxBR  
B'/ >Ax&  

看看谢谢楼主

观望中。。。。。

学习学习