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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ t,0}}9%?  
q@VIFmqY!  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 hPGDN\#LD  
Y$)y:.2#  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 !Z{7X ^  
3A~<|<}t  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： ]-a
/)8  
G-]<+-Q$4  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ []yIz1P=j  
I3sf
O
U  
以下是我的计算命令流，希望得到高手的帮助啊！ +u5xK  
RF -c`C  
d|~'#:y@  
Fini 0@EI@X;q  
/clear 1~2+w]-kU  
/title,TUNNEL ANALYSIS x5;D'Y t"|  
/prep7 yzL9Ic  
ET,1,PLANE2 >){}nlQf  
keyopt,1,3,2 $/(``8li_  
S+mBVk"-~S  
*SET,R,0.7 r9b(d]  
*SET,H,10.0 42,dH
Ydt  
*SET,X1,-2*R                                                    K9VP@[zbJ  
*SET,Y1,-2*R                                                    CQ Ei(ty  
*SET,W1,4*R                                                   i=rA;2>  
*SET,H1,4*R                                                   @(N} {om  
                                                              86?~N  
*SET,X2,-10.0                                                   J_|%8N{[x  
*SET,Y2,H-20.0                                                    M\JAB ;A  
*SET,W2,20.0                                                  R `ob;>[Q  
*SET,H2,20.0 xH,e$t#@@~  
CYL4,0,0,0,45,R,405,0 OH]45bd &7  
BLC4,X1,Y1,W1,H1                                              ~-%z:Re'_  
BLC4,X2,Y2,W2,H2 ]+w 27!  
AOVL,1,2,3                                              BV/ ^S.~  
NUMCMP,ALL                                                    oz|+{b}%  
NUMMRG,ALL,,,,LOW ]z=dRq  
LANG,11,8,90,, ^)|tf\4  
LANG,11,7,90,, dPV<:uO  
LANG,10,7,90,,  zNn  
LANG,10,6,90,, QMX  
LANG, 9,6,90,, qM3(OvCt  
LANG, 9,5,90,, MA%g-}  
LANG,12,5,90,, v9f%IE4fX  
LANG,12,8,90,, d5{RIM|  
L,1,7 wk^$DM/KJ)  
L,2,8 "dYT>w  
L,3,5 dg+"G|nr  
L,4,6 Y$|KY/)H)  
NUMCMP,ALL                                                    fti0Tz'  
NUMMRG,ALL,,,,LOW mOyNl -f  
FLST,3,4,4,ORDE,2   w=ufJRj  
FITEM,3,29  Ia{t/IX\[  
FITEM,3,-32 ,#0#1k<Dm  
ASBL,2,P51X  o62GEl25  
FLST,3,8,4,ORDE,8   3fYfj  
FITEM,3,14  W~l.feW$i  
FITEM,3,16  *kj+6`:CPs  
FITEM,3,18  I.SMn,N  
FITEM,3,20  W55kR.X6M  
FITEM,3,22  N>sHT =_  
FITEM,3,24  {8%KO1xB  
FITEM,3,26  p"*xyex  
FITEM,3,28  ,Vz-w;oDn  
ASBL,       3,P51X 3YUF\L]yyw  
NUMCMP,ALL                                                    ^0
I"  
NUMMRG,ALL,,,,LOW Gsm.a  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 ,d^HAg^j  
MP,EX,1,3.0E6                                    lx~C{tl2  
MP,PRXY,1,0.33                                                o?#-Tkb  
MP,DENS,1,1900.0 Q } 0_}W  
TB,DP,1                                                       QV#HN"F/K  
TBDATA,1,2.4E4,13.6,0 ]\ZmK0q<:  
(~"#=fs.L  
MP,EX,2,3.0E6                                    /}eb1o  
MP,PRXY,2,0.33                                                01]W@\(  
MP,DENS,2,1900.1 D7jbo[GgS  
TB,DP,2                                                       qR%as0;  
TBDATA,1,2.4E4,13.6,0 "412w^5[T  
j:^gmZ;J  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 \t=#MzjR  
type,1 $g,v]MW  
MAT,1  N6*FlG-  
lsel,S,line,,1,8,1 7$R^u7DZ  
lsel,A,line,,15,27,4     _vAc/_N  
Lesize,ALL,,,4,,,,,0 {817Svp@  
lsel,S,line,,29,32,1     DW)81*~g  
Lesize,ALL,,,3,,,,,0 PT4iy<  
  '@^mesMG  
lsel,S,line,,9,13,2 TeJ=QpGW2  
lsel,A,line,,18,20,2 _ZK*p+u%  
lsel,A,line,,21 /f=31<+MtF  
lsel,A,line,,26,28,2     "837b/>/  
Lesize,ALL,,,10,,,,,0 <?Z]h]C^o  
lsel,S,line,,10 -V'Y
^Df  
lsel,A,line,,22 LGYg@DR  
lsel,A,line,,24,25,1 %9L+ Q1o  
Lesize,ALL,,,20-H-2*R,,,,,0 S.U#lAn(  
lsel,S,line,,12,14,2 %w/vKB"nO  
lsel,A,line,,16,17,1 _]0<G8|Rv  
Lesize,ALL,,,H-2*R,,,,,0 &IG*;$c!  
asel,s,area,,2,13,1  O`9vEovjs  
amesh,all 1V,DcolRY  
   A":cS }Ui  
type,1 ~svea>Fmr  
MAT,2 +Je%8jH  
asel,s,area,,1  }YU\}T-P  
amesh,all W*u$e8i7  
nummrg,all,,,,low "_+8z_  
numcmp,all Dh68=F0  
r#J_;P{U  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 iBY16_q  
/solu j:HIcCp  
antype,static z)u\(W*\iA  
deltim,0.1,0.05,0.2 X'/'r
.b6  
autots,on !使用自动时间步 `"#hhKG  
pred,on !打开时间步长预测器 60xL.Z   
lnsrch,on !打开线性搜索 @"aqnj>+  
nlgeom,on !打开大位移效果 wOEc~WOd  
nropt,full !设定牛顿-拉普森选项 wF*9%K'E  
cnvtol,f,,0.02,0 m mH xPd  
Nsel,s,loc,x,-10 wEU=R>j
.  
Nsel,a,loc,x,10 b4(,ls  
D,all,,,,,,UX,,,,, >h$Q%w{V  
Nsel,s,loc,y,H-20 o"p['m*g  
D,all,,,,,,Uy,,,,, =9:gW5F69  
time,1 @J'tPW<$  
acel,,10  2r@9|}La  
nsel,all 6|4ID"  
esel,all K[Vj+qdyl  
allsel,all E!=Iz5
 
solve Ns\};j?TU*  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 (543`dqAmC  
time,2 wZ_"@j<  
esel,s,type,,2 }V 4u`=  
ekill,all l,UOP[j  
nsel,all u Vo"_c w  
esel,all d@ ]N  
allsel,all [<wpH0lNoy  
solve 2g5jGe*0  
fini $gCN[%+j  
oE,TA2  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 x+[ATZ([  
!/POST1   VU+=b+B~m  
!LCDEF,1,1,, O;0VKNn['  
!LCDEF,2,2,, C;_00EQ=  
!LCASE,2, @yF>=5z:  
!LCOPER,SUB,1, ,2, blkPsp)m"  
!Nsel,s,loc,Y,H )OK"H^}f  
!PLVECT,U, , , ,VECT,ELEM,ON,0    TQQh:y  
!PLNSOL, U,Y, 0,1.0 7zowvE?#  
bR?xz-g%<3  
SErh"~[  
^Iqu^n?2.  
equi26jhr  

看看谢谢楼主

观望中。。。。。

学习学习