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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ [!-gb+L  
IhYR4?e  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 ZcQu9XDIt  
_Uq' N0U  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 FutS  
7K:V<vX5  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： mpMAhm
:  
G7Sw\wW  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ Wi'}d6c  
]MosiMJF  
以下是我的计算命令流，希望得到高手的帮助啊！ .  yg#  
#c>GjUJ.w  
5$D"uAp<V  
Fini |f#hGk6  
/clear E<[ s+iX  
/title,TUNNEL ANALYSIS J5I@*f)l  
/prep7 qul#)HI  
ET,1,PLANE2 voRry6Q;  
keyopt,1,3,2 EN
2SI
+  
AXW.`~ 4  
*SET,R,0.7 pB
8D  
*SET,H,10.0 =I@t%Y  
*SET,X1,-2*R                                                    ]
f0'YLG  
*SET,Y1,-2*R                                                    @we1#Vz.  
*SET,W1,4*R                                                   R=lw}jH[Z  
*SET,H1,4*R                                                   q!e
E~O;A  
                                                              y]m: {  
*SET,X2,-10.0                                                   L ,dh$F  
*SET,Y2,H-20.0                                                    /[.V(K D  
*SET,W2,20.0                                                  !S#3mT-  
*SET,H2,20.0 7|DG1p9C  
CYL4,0,0,0,45,R,405,0 7GYf#} N  
BLC4,X1,Y1,W1,H1                                              <f>w"r  
BLC4,X2,Y2,W2,H2 V0>X2&.A  
AOVL,1,2,3                                              +*]$PVAFA  
NUMCMP,ALL                                                    EIg~^xK  
NUMMRG,ALL,,,,LOW T8x)i\<  
LANG,11,8,90,, 7"!b5(4=  
LANG,11,7,90,, b-sN#'TDg  
LANG,10,7,90,, f +{=##'0  
LANG,10,6,90,, gwRB6m$  
LANG, 9,6,90,, alaL/p{O  
LANG, 9,5,90,, (WyNO QO'  
LANG,12,5,90,, l<f9$l^U
 
LANG,12,8,90,, 10Ik_L='  
L,1,7 ERpAV-Zf  
L,2,8 oLT#'42+H  
L,3,5 j|k/&q[St  
L,4,6 ahg:mlaob  
NUMCMP,ALL                                                    nn_O"fZi  
NUMMRG,ALL,,,,LOW %8Ff
P5#  
FLST,3,4,4,ORDE,2   /\w4k  
FITEM,3,29  %^IQ<  
FITEM,3,-32 )8@-  
ASBL,2,P51X  I'
4(Ibl+  
FLST,3,8,4,ORDE,8   N3n]  
FITEM,3,14  &Q883A J  
FITEM,3,16  Lt't   
FITEM,3,18  i6'=]f'{  
FITEM,3,20  lFnls6dp  
FITEM,3,22  0Z1ksfLU  
FITEM,3,24  &zdS9e-fF  
FITEM,3,26  b>;?{  
FITEM,3,28  ]#)(D-i  
ASBL,       3,P51X t0m;tb bg  
NUMCMP,ALL                                                     3%kUj  
NUMMRG,ALL,,,,LOW #'_#t/u  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 e= IdqkJ%  
MP,EX,1,3.0E6                                    79v+ze  
MP,PRXY,1,0.33                                                .Vbd-jr'M  
MP,DENS,1,1900.0 0K%okq|n  
TB,DP,1                                                       *[Z`0AgP  
TBDATA,1,2.4E4,13.6,0 Z-PBCU  
WNGX`V,d  
MP,EX,2,3.0E6                                    MXaFqK<Y  
MP,PRXY,2,0.33                                                pl`4&y%Me  
MP,DENS,2,1900.1 ,8@<sFB'  
TB,DP,2                                                       D$hQ-K  
TBDATA,1,2.4E4,13.6,0 9zwD%3Ufn  
^wIB;!W  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 MzIDeZ  
type,1 d@R7b^#g  
MAT,1  P0l fK}  
lsel,S,line,,1,8,1
1<r!9x9G  
lsel,A,line,,15,27,4     oy^-?+  
Lesize,ALL,,,4,,,,,0 FS1\`#Bm)  
lsel,S,line,,29,32,1     4YfM.~ 6  
Lesize,ALL,,,3,,,,,0 q`*.F#/4c  
  rmX*s}B  
lsel,S,line,,9,13,2 5Z>a}s_i  
lsel,A,line,,18,20,2 Q,3kaR@O  
lsel,A,line,,21 *D?=Ts  
lsel,A,line,,26,28,2     6!\V|  
Lesize,ALL,,,10,,,,,0 w2 Y%yjCV  
lsel,S,line,,10 uy{O  
lsel,A,line,,22 U@Od
QAX  
lsel,A,line,,24,25,1 @+b$43^  
Lesize,ALL,,,20-H-2*R,,,,,0 N%dY.Fk  
lsel,S,line,,12,14,2 C+NN.5No  
lsel,A,line,,16,17,1 FRQ.ix2  
Lesize,ALL,,,H-2*R,,,,,0 jY=y<R_oK  
asel,s,area,,2,13,1  Bb/if:XS  
amesh,all {`!6w>w0  
   \3JCFor/  
type,1 d{S'6*`D  
MAT,2 IuXgxR%  
asel,s,area,,1  qp})4XTv  
amesh,all Wk%|%/
:  
nummrg,all,,,,low ?Lg<)B9   
numcmp,all X _ZO)|  
:^)?AO#J  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 <:9ts@B  
/solu t2dsYU/  
antype,static pJHdY)Cz  
deltim,0.1,0.05,0.2 ]?y~;-^  
autots,on !使用自动时间步 5iA>Z!sP[  
pred,on !打开时间步长预测器 cW%)C.M  
lnsrch,on !打开线性搜索 jFN0xGZ  
nlgeom,on !打开大位移效果 +#Pb@^6"m  
nropt,full !设定牛顿-拉普森选项 *Jsb~wta  
cnvtol,f,,0.02,0 h#YO;m2wd  
Nsel,s,loc,x,-10 (<Cq_Kw  
Nsel,a,loc,x,10 ;J W]b]  
D,all,,,,,,UX,,,,, |vs5N2_  
Nsel,s,loc,y,H-20 DTezG':  
D,all,,,,,,Uy,,,,, w/#7G\U  
time,1 UR'[?  
acel,,10  s<YN*~  
nsel,all @[r[l#4yUi  
esel,all kLgkUck8]  
allsel,all >?^oxB"<Gc  
solve \IL)~5d  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 DT]4C!dh  
time,2 RL`E}:V  
esel,s,type,,2
v2;E Wp  
ekill,all h@ lz  
nsel,all al[^pPK
Z  
esel,all nah?V" ?Y  
allsel,all -}6ew@GE  
solve IWo'{pk  
fini u9"kF  
kkl'D!z2g  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 l[E^nh>  
!/POST1   .YF1H<gwa  
!LCDEF,1,1,, }b2YX+/e$f  
!LCDEF,2,2,, +R*DE5dz  
!LCASE,2, ai<qK3!O  
!LCOPER,SUB,1, ,2, s{^B98d+W  
!Nsel,s,loc,Y,H %6Gg&Y$j!  
!PLVECT,U, , , ,VECT,ELEM,ON,0    _#{*I(l  
!PLNSOL, U,Y, 0,1.0 T_*inPf  
yY#h1  

&PC6C<<f  
>w.;A%|N  
m"!SyN}&9?  

看看谢谢楼主

观望中。。。。。

学习学习