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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ E^#|1Kpq  
GAtK1%nPD  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 :#
I8Cf  
mmr>"`5.  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 cd*y{Wt  
,LWM}L  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： m'oVqA&  
:h*20iP  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ Joq9.%7Q  
-5kq9Dy\,  
以下是我的计算命令流，希望得到高手的帮助啊！ 8+v6%,K2  
sVaWg?=qs'  
{Kd9}CDAZ  
Fini H>;km$b +  
/clear fx
%'7/+  
/title,TUNNEL ANALYSIS mkrvWZjZX  
/prep7 CG$S?  
ET,1,PLANE2 BAg*zYV7  
keyopt,1,3,2 M1Od%nz3  
<w.V!"!  
*SET,R,0.7 )Qb1$%r.  
*SET,H,10.0 |PDuvv!.f  
*SET,X1,-2*R                                                    SU2(XP]5  
*SET,Y1,-2*R                                                    hFj
.d]S  
*SET,W1,4*R                                                   (al7/EhY  
*SET,H1,4*R                                                   R 5bt~U  
                                                              ^=E4~
22q  
*SET,X2,-10.0                                                   DV*8Mkzg  
*SET,Y2,H-20.0                                                    j\!zz  
*SET,W2,20.0                                                  RAXqRP,iw  
*SET,H2,20.0 !0*=z~  
CYL4,0,0,0,45,R,405,0 6bo,x  
BLC4,X1,Y1,W1,H1                                              =EsKFt"  
BLC4,X2,Y2,W2,H2 T?^AllUZQR  
AOVL,1,2,3                                              B;hc|v{(  
NUMCMP,ALL                                                    SX/yY  
NUMMRG,ALL,,,,LOW 0%`\8  
LANG,11,8,90,, =
?vk n  
LANG,11,7,90,, q;tsA"l  
LANG,10,7,90,, f1hi\p0q  
LANG,10,6,90,, 
(fm\kV  
LANG, 9,6,90,, ldanM>5  
LANG, 9,5,90,, = J).(E89  
LANG,12,5,90,, >sPu*8D40a  
LANG,12,8,90,, 1\kOjF)l  
L,1,7 tN";o\!}  
L,2,8 J A4'e@  
L,3,5 fcD$km  
L,4,6 5|S|HZ8G  
NUMCMP,ALL                                                    u%VO'}Gz  
NUMMRG,ALL,,,,LOW qV9`  
FLST,3,4,4,ORDE,2   f![x7D$  
FITEM,3,29  `S{< $:D  
FITEM,3,-32 k[y{&f,  
ASBL,2,P51X  burEo
.=  
FLST,3,8,4,ORDE,8   6~;fj+S  
FITEM,3,14  q
,$UKg#i  
FITEM,3,16  a5L#c
=  
FITEM,3,18  wC'KI8-  
FITEM,3,20  :]=Y1*L\)  
FITEM,3,22  ZUAWSJ,s  
FITEM,3,24  _6^vxlF  
FITEM,3,26  sB-c'`,w`  
FITEM,3,28  7b:oz3?PI  
ASBL,       3,P51X (3{'GX2c  
NUMCMP,ALL                                                    |C7GI[P  
NUMMRG,ALL,,,,LOW
=u${2=  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 S)@R4{=e"V  
MP,EX,1,3.0E6                                    ?ti7iBz?  
MP,PRXY,1,0.33                                                >:X
zv  
MP,DENS,1,1900.0 }9<aX Y,  
TB,DP,1                                                       HBo^8wN
 
TBDATA,1,2.4E4,13.6,0 eV1O#FLbi  
!+9H=u  
MP,EX,2,3.0E6                                    H:d{Sru  
MP,PRXY,2,0.33                                                NVeb,Pf  
MP,DENS,2,1900.1 4#;rv$ {  
TB,DP,2                                                       i+Ob1B@w
 
TBDATA,1,2.4E4,13.6,0 T!(I\wz;Bo  
@Eqc&v!O  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 vlp]!7v  
type,1 g%1!YvS3v  
MAT,1  PIB|&I|p  
lsel,S,line,,1,8,1 91mXvQ:u  
lsel,A,line,,15,27,4     V|#B=W  
Lesize,ALL,,,4,,,,,0 Xdq2.:\  
lsel,S,line,,29,32,1     Qaq{UW
 
Lesize,ALL,,,3,,,,,0 T1\Xz
-1  
  7$,["cJX  
lsel,S,line,,9,13,2 }_@cqx:n^  
lsel,A,line,,18,20,2 L>xcgV7  
lsel,A,line,,21 i*W8_C:S  
lsel,A,line,,26,28,2     [UR+G8X21m  
Lesize,ALL,,,10,,,,,0 w v9s{I{P  
lsel,S,line,,10 5}e-\:J>B  
lsel,A,line,,22 e%(zjCA  
lsel,A,line,,24,25,1 ~;wSe[  
Lesize,ALL,,,20-H-2*R,,,,,0 ~9h6"0K!  
lsel,S,line,,12,14,2 1K09iB  
lsel,A,line,,16,17,1 XrFyN(p  
Lesize,ALL,,,H-2*R,,,,,0 8T$:^HW  
asel,s,area,,2,13,1  zP rT0  
amesh,all gC<\1AIu  
   JWlH(-U4|  
type,1 [)Xu60?Q  
MAT,2 Ud`V"X  
asel,s,area,,1  pWbzBgM?nU  
amesh,all p^5B_r:  
nummrg,all,,,,low iDp]lu  
numcmp,all xm/v:hl=  
E Jq=MP  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 }@SZ!-t%rD  
/solu H6bomp"  
antype,static ruB&&C6)v  
deltim,0.1,0.05,0.2 V1x
pJ  
autots,on !使用自动时间步 sZ]O&Za~  
pred,on !打开时间步长预测器 6\OSIxJZF  
lnsrch,on !打开线性搜索 mZ ONxR6q$  
nlgeom,on !打开大位移效果 &"Ua"H)  
nropt,full !设定牛顿-拉普森选项 N!9DZEcm  
cnvtol,f,,0.02,0 (@uQ>dR:  
Nsel,s,loc,x,-10 ^dYFFKQ  
Nsel,a,loc,x,10 g0cCw2S  
D,all,,,,,,UX,,,,,
ZJ=-cE
2n  
Nsel,s,loc,y,H-20 UyD=x(li  
D,all,,,,,,Uy,,,,, H Y.,f_m  
time,1 H,:Cg:E/^  
acel,,10  <4C`^p  
nsel,all b;9v.MZ4>g  
esel,all `$G7Ia_ $]  
allsel,all ogip#$A}3  
solve XRJ<1w:  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 o=q N+-N  
time,2 ;T-i+_  
esel,s,type,,2 {
~
b]6}O  
ekill,all kOFEH!9&  
nsel,all vb{&T<  
esel,all scJ`oc:<J  
allsel,all i,4  
solve \dbpCZ  
fini (<itE3P  
Vu^J'>X  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 ]/JE#  
!/POST1   +qD4`aI 
 
!LCDEF,1,1,, g]X4)e]  
!LCDEF,2,2,, o PR^Z pt  
!LCASE,2, oel3H5Nz  
!LCOPER,SUB,1, ,2, K ,f1c}  
!Nsel,s,loc,Y,H _o' jy^  
!PLVECT,U, , , ,VECT,ELEM,ON,0    #
s(B,`?N  
!PLNSOL, U,Y, 0,1.0 W{1=O)w  
<W|{zAyv  
Fl(+c0|kT  
]rZ"5y  
W\N-~9UA  

看看谢谢楼主

观望中。。。。。

学习学习