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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ d7X&3L%Oq  
:-k|jt  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 EbQLMLD%  
;Y`Y1  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 '3(^Zv  
fo@^=-4A-  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： G-Tmk7m  
pD732L@q  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ }#O!GG{  
9RaO[j`  
以下是我的计算命令流，希望得到高手的帮助啊！ oY18a*_>M1  
(G>[A}-  
}p7iv:P=3  
Fini Z nc(Q  
/clear }6c>BU}DF  
/title,TUNNEL ANALYSIS eyJ07  
/prep7 `um,S  
ET,1,PLANE2 ^[tE^(|T  
keyopt,1,3,2 pFd8p@m_2  
Bvn3:+(47  
*SET,R,0.7 "n!yK
  
*SET,H,10.0 YZBh}l6t  
*SET,X1,-2*R                                                    d]l8ei@>h  
*SET,Y1,-2*R                                                    kW g.-$pp  
*SET,W1,4*R                                                   e{P v:jl  
*SET,H1,4*R                                                   c0@8KW[,  
                                                              -Ks>s  
*SET,X2,-10.0                                                   lS.Adl^k  
*SET,Y2,H-20.0                                                    w6% Q"%rp  
*SET,W2,20.0                                                  aKbmj  
*SET,H2,20.0 #dA$k+3  
CYL4,0,0,0,45,R,405,0 6gg8h>b  
BLC4,X1,Y1,W1,H1                                              }9'`3vsJ  
BLC4,X2,Y2,W2,H2 $E\|\g  
AOVL,1,2,3                                              :jLL IqhB  
NUMCMP,ALL                                                     Z.6dL  
NUMMRG,ALL,,,,LOW -1g:3'% P  
LANG,11,8,90,, hi0HEm\  
LANG,11,7,90,, 8-#%l~dr  
LANG,10,7,90,, G21o@38e  
LANG,10,6,90,, g\JJkXjD#  
LANG, 9,6,90,, yp.K-  
LANG, 9,5,90,, V0\[|E;F  
LANG,12,5,90,, H$(bSw$  
LANG,12,8,90,, Wp=3heCa6  
L,1,7 &N|$G8\CY  
L,2,8 ATNOb  
L,3,5 Iry$z^  
L,4,6 1PkCWRpR  
NUMCMP,ALL                                                    R2~Tr$:  
NUMMRG,ALL,,,,LOW 5yy:JTAH5  
FLST,3,4,4,ORDE,2   iEr,ly  
FITEM,3,29  `C+<!)2  
FITEM,3,-32 []>'Dw_r  
ASBL,2,P51X  @!#e\tx  
FLST,3,8,4,ORDE,8   FRBu8WW0L  
FITEM,3,14  T pkSY`T  
FITEM,3,16  n{;j  
FITEM,3,18  qos7u91z  
FITEM,3,20  )u)=@@k21  
FITEM,3,22  u*l|MIi6J  
FITEM,3,24  &7aWVKon  
FITEM,3,26  <$(B[T  
FITEM,3,28  e`D}[G#  
ASBL,       3,P51X ^/2I)y]W0  
NUMCMP,ALL                                                    L[G O6l  
NUMMRG,ALL,,,,LOW @p+;iS1}  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 9oWU]A\k>  
MP,EX,1,3.0E6                                    %iN>4;T8  
MP,PRXY,1,0.33                                                !+T1kMP+l  
MP,DENS,1,1900.0 &v{Ehkr*  
TB,DP,1                                                       |W">&Rb<t#  
TBDATA,1,2.4E4,13.6,0 zH8E,)  
@c3xUK  
MP,EX,2,3.0E6                                    NI >%v  
MP,PRXY,2,0.33                                                &_ekA44E  
MP,DENS,2,1900.1 4>hHUz[_  
TB,DP,2                                                       QU|_ r2LM  
TBDATA,1,2.4E4,13.6,0 iq(PC3e`V  
a:h<M^n049  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 'pdTV:]zA  
type,1 cj/`m$  
MAT,1  XIHN6aQ{X  
lsel,S,line,,1,8,1 I{`70  
lsel,A,line,,15,27,4     _!\d?]Ya  
Lesize,ALL,,,4,,,,,0 wHc my  
lsel,S,line,,29,32,1     Lx^ eaP5  
Lesize,ALL,,,3,,,,,0 HGDrH   
  /U~|B
.z@6  
lsel,S,line,,9,13,2 l
90mM'[  
lsel,A,line,,18,20,2 \*xB<mq  
lsel,A,line,,21 :]^e-p!z  
lsel,A,line,,26,28,2     >Jz9wo`  
Lesize,ALL,,,10,,,,,0 ~&?bU]F  
lsel,S,line,,10 y>^^.  
lsel,A,line,,22 &XP(D5lf`B  
lsel,A,line,,24,25,1 IHl q27O  
Lesize,ALL,,,20-H-2*R,,,,,0 Bh>L"'.2  
lsel,S,line,,12,14,2 ^OR0Vp>L  
lsel,A,line,,16,17,1 d8j1L/e  
Lesize,ALL,,,H-2*R,,,,,0 N@q}eGe  
asel,s,area,,2,13,1  P#,u9EIJ  
amesh,all g`7XE  
   QHEtG2  
type,1 "F<CGSo  
MAT,2 kmI0V
[Y  
asel,s,area,,1  BX,)G HE  
amesh,all >ZgzE  
nummrg,all,,,,low Aw o)a8e  
numcmp,all z$32rt8{`v  
RK>Pe3<  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 k_al*iM>H  
/solu K7+yU3  
antype,static 1j<(?MT-  
deltim,0.1,0.05,0.2 WSkGVQu  
autots,on !使用自动时间步 z^gJy,T  
pred,on !打开时间步长预测器 6o3 bq|  
lnsrch,on !打开线性搜索 K}VCFV  
nlgeom,on !打开大位移效果 mPV<a&U  
nropt,full !设定牛顿-拉普森选项 O !L`0 =%c  
cnvtol,f,,0.02,0 kSQ8kU_w+  
Nsel,s,loc,x,-10 VM"cpC_8  
Nsel,a,loc,x,10 Ccf/hA#mb  
D,all,,,,,,UX,,,,, AZtS4]4G)  
Nsel,s,loc,y,H-20 +eM${JyXH  
D,all,,,,,,Uy,,,,, a|aVc'j  
time,1 /0h *(nL  
acel,,10  rH+OXGoB  
nsel,all <j'V}|3  
esel,all 3FEJ 9ZyG  
allsel,all l]nt@0+  
solve b'H'QY   
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 _FLEz|%~  
time,2 RpHlq  
esel,s,type,,2 ^.SY
AwL  
ekill,all ;?-AFd\i
 
nsel,all RgE`Hr  
esel,all o`?rj!\  
allsel,all "/#JC}]  
solve woYD &Oml  
fini tT$OnZu&  
?9b9{c'an  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 &1xCPKIr  
!/POST1   +]db-  
!LCDEF,1,1,, xvr5$x|h  
!LCDEF,2,2,, *2'8d8>R%]  
!LCASE,2, %`vzQt`>  
!LCOPER,SUB,1, ,2, K"}fD;3  
!Nsel,s,loc,Y,H w2)Ro:G  
!PLVECT,U, , , ,VECT,ELEM,ON,0    _]Hna<Ly  
!PLNSOL, U,Y, 0,1.0 ou|emAV  
BT$p~XB  
DX>a0-Xj  
n/H OP  
L[` l80  

看看谢谢楼主

观望中。。。。。

学习学习