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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ NO6.qWl  
{Y\hr+A  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 OP
Km^}  
ku`'w;5jT  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 #7'ww*+  
%f, 9  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： :X@;XEol~
 
YB)1dzU  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ M@csB.'  
Lw-j#}&6E  
以下是我的计算命令流，希望得到高手的帮助啊！ oYOf<J  
l*%?C*  
zW+Y{^hf  
Fini o#w6]Fmc  
/clear
U:8^>_  
/title,TUNNEL ANALYSIS s31^9a  
/prep7 \@I.K+hj$  
ET,1,PLANE2 !y*V;J  
keyopt,1,3,2 e2Xx7*vS  
{|>~#a49h  
*SET,R,0.7 S;*,V|#QD  
*SET,H,10.0 }V9146  
*SET,X1,-2*R                                                    Wx:_F;  
*SET,Y1,-2*R                                                    WhK?>u  
*SET,W1,4*R                                                   62Jn8DwAT  
*SET,H1,4*R                                                   >%d]"]  
                                                              36lIV,YnU  
*SET,X2,-10.0                                                   G::6?+S  
*SET,Y2,H-20.0                                                    +n)(\k{
 
*SET,W2,20.0                                                  cL=P((<K?  
*SET,H2,20.0 8f29Hj+  
CYL4,0,0,0,45,R,405,0 m<OxO\Mpf  
BLC4,X1,Y1,W1,H1                                              l;?.YtMg  
BLC4,X2,Y2,W2,H2 >)5rOU  
AOVL,1,2,3                                              B&L{/.v_z\  
NUMCMP,ALL                                                    n>iPAD  
NUMMRG,ALL,,,,LOW L{(r@Vu  
LANG,11,8,90,, @W^g(I(w  
LANG,11,7,90,, r$0=b -  
LANG,10,7,90,, u1nv'\*  
LANG,10,6,90,, yv'mV=BMJ!  
LANG, 9,6,90,, vtr:{   
LANG, 9,5,90,, -3KB:K<  
LANG,12,5,90,, 21cIWvy  
LANG,12,8,90,, R HF;AX n  
L,1,7 R%}k52`  
L,2,8 ?g<*1N?:  
L,3,5 !3T x\a`?/  
L,4,6 D Zh6/n#q  
NUMCMP,ALL                                                    /=Xen mmS  
NUMMRG,ALL,,,,LOW 6`+DBr  
FLST,3,4,4,ORDE,2   XFBk:~}sI  
FITEM,3,29  -'YX2!IU,  
FITEM,3,-32 w7%N=hL1
 
ASBL,2,P51X  5[*MT%ms  
FLST,3,8,4,ORDE,8   vtvr{Uqo@  
FITEM,3,14  TF3Tha]  
FITEM,3,16  vTB*J,6.  
FITEM,3,18  {5_*f)$[H  
FITEM,3,20  _
&BnET  
FITEM,3,22  gXf_~zxS  
FITEM,3,24  sIELkF?.  
FITEM,3,26  ClfpA?vv  
FITEM,3,28  4.VEE~sH$  
ASBL,       3,P51X AM Rj N;  
NUMCMP,ALL                                                    9jvg[H  
NUMMRG,ALL,,,,LOW :0srFg?X  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 [N$@nA-d  
MP,EX,1,3.0E6                                    g%\e80~1(  
MP,PRXY,1,0.33                                                O|gb{  
MP,DENS,1,1900.0 Sb{S^w\m0  
TB,DP,1                                                       4N
Ek#n  
TBDATA,1,2.4E4,13.6,0 U&B~GJT+  
%X.Q\T  
MP,EX,2,3.0E6                                    o}Xp-P   
MP,PRXY,2,0.33                                                Ex*g>~e  
MP,DENS,2,1900.1 WG^D$L
:  
TB,DP,2                                                       W]y$6P  
TBDATA,1,2.4E4,13.6,0 gi"v${R  
[s7I.rdGzz  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 4wNxn lP  
type,1 `zep`j&8^  
MAT,1  _Juhl^LM;  
lsel,S,line,,1,8,1 *w%;$\^  
lsel,A,line,,15,27,4     &1Az`[zKGW  
Lesize,ALL,,,4,,,,,0 uAJ_`o[  
lsel,S,line,,29,32,1     ob*2V!"  
Lesize,ALL,,,3,,,,,0 Gn_v}31d%  
  bFflA  
lsel,S,line,,9,13,2 -7+Fb^"L  
lsel,A,line,,18,20,2 }  g  
lsel,A,line,,21 m\~[^H~g  
lsel,A,line,,26,28,2     \SMH",u  
Lesize,ALL,,,10,,,,,0 j5ui  
lsel,S,line,,10 WcUeWGC>  
lsel,A,line,,22 :A 1,3g  
lsel,A,line,,24,25,1 GZ4{<QG  
Lesize,ALL,,,20-H-2*R,,,,,0 bUWtlg  
lsel,S,line,,12,14,2 mKn[>M1  
lsel,A,line,,16,17,1 r@n%  
Lesize,ALL,,,H-2*R,,,,,0 :;|x'[JoE?  
asel,s,area,,2,13,1  HX <;=m  
amesh,all IGI2).$[  
   r+Z+x{  
type,1 ;eA~z"g  
MAT,2 ev+NKUi=  
asel,s,area,,1  f<zh-Gq  
amesh,all "Dy&`  
nummrg,all,,,,low Fs=nAn#  
numcmp,all 6Es-{u(,  
4e5Ka{# <  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 L'Wcb =;  
/solu 8T2$0  
antype,static |E}N8\Gr  
deltim,0.1,0.05,0.2 +-{HT+W  
autots,on !使用自动时间步 WIbU^WJ0  
pred,on !打开时间步长预测器 %iX+"  
lnsrch,on !打开线性搜索 H-S28%.  
nlgeom,on !打开大位移效果 zz$*upxK  
nropt,full !设定牛顿-拉普森选项 \"uR&D  
cnvtol,f,,0.02,0 LOOv8'%O8  
Nsel,s,loc,x,-10
=izB :  
Nsel,a,loc,x,10 '8W}|aF  
D,all,,,,,,UX,,,,, 5&VLq  
Nsel,s,loc,y,H-20 ?Jx8z`(  
D,all,,,,,,Uy,,,,, d:j$!@o  
time,1 $60+}B`m  
acel,,10  Ig=4Z*au!g  
nsel,all 6x.#K9@q4  
esel,all ~+|p.(I  
allsel,all <Siz5qQI4  
solve y6Rg@L&U  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 %f^TZ,q$  
time,2 b0 5h,  
esel,s,type,,2 ;,F}!R  
ekill,all ho0T$hB  
nsel,all D!y Cnq=8  
esel,all rv\<Q-uQ8  
allsel,all 2R)Y}*VX  
solve <7)@Jds\  
fini kgA')]  
W;o\}irep  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 'g$|:bw/  
!/POST1   `4$" mO>+  
!LCDEF,1,1,, 2'/ ip@  
!LCDEF,2,2,, ty
"k  
!LCASE,2, jaDZPX-yS  
!LCOPER,SUB,1, ,2, Ypp>7J/  
!Nsel,s,loc,Y,H J\Sewg9  
!PLVECT,U, , , ,VECT,ELEM,ON,0    6>KDK<5NQ  
!PLNSOL, U,Y, 0,1.0 gTs5xDvJ  
oS]XE!
^M  
\=nY&Ml  
J%%nv5y  
6@nEcr  

看看谢谢楼主

观望中。。。。。

学习学习