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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ @g|v
;B|{  
seba9y  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 CYt?,qk-r  
N'F77 .  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 tCw<Ip  
%3s1z<;R[S  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： Z3%}ajPu[  
.7O*pJ2(H  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ 3D6RLu
  
4M^=nae  
以下是我的计算命令流，希望得到高手的帮助啊！ <3L5"77G6  
Dxtp2wu%t  
S};#+ufgTt  
Fini B[qzUD*P_n  
/clear WEAT01  
/title,TUNNEL ANALYSIS f@6QvkIa  
/prep7 Jf@H/luW  
ET,1,PLANE2 ]z;P9B3@&  
keyopt,1,3,2 )tG\vk=@  
NxfOF  
*SET,R,0.7 O],T,Z?z  
*SET,H,10.0 "x%Htq@  
*SET,X1,-2*R                                                    )z>|4@,  
*SET,Y1,-2*R                                                    i)\L:qF5  
*SET,W1,4*R                                                   2L!u1  
*SET,H1,4*R                                                   OwuE~K7b{  
                                                              Fzm*Pz3  
*SET,X2,-10.0                                                   ;:iY)}  
*SET,Y2,H-20.0                                                    3N8t`N  
*SET,W2,20.0                                                  "WlZ)wyF%  
*SET,H2,20.0 ~cWAl,(B<F  
CYL4,0,0,0,45,R,405,0 _,
kj:R.  
BLC4,X1,Y1,W1,H1                                              CMe 06^U  
BLC4,X2,Y2,W2,H2 rr/B=O7  
AOVL,1,2,3                                              ^0Zf,40  
NUMCMP,ALL                                                    B)$| vK=  
NUMMRG,ALL,,,,LOW y@1+I~@  
LANG,11,8,90,, #HYr0Tw6`  
LANG,11,7,90,, Nv$R\'3  
LANG,10,7,90,, W'els)WJ|x  
LANG,10,6,90,, U1Q:= yD  
LANG, 9,6,90,, Yy5h"r  
LANG, 9,5,90,, ?xA:@:l/  
LANG,12,5,90,, '-QwssE  
LANG,12,8,90,, (XVw"m/ye  
L,1,7 6e(Qwt  
L,2,8 \Bz_p'[G  
L,3,5 Q[Z8ok  
L,4,6 `K*Q5n  
NUMCMP,ALL                                                    % BVs47g  
NUMMRG,ALL,,,,LOW [2>yYr s_=  
FLST,3,4,4,ORDE,2   %W'v}p  
FITEM,3,29  60P<4
 
FITEM,3,-32 j 7fL7:,T  
ASBL,2,P51X  \$R_YKGf1G  
FLST,3,8,4,ORDE,8   >/{@C  
FITEM,3,14  w2Us!<x  
FITEM,3,16  >f^r^P  
FITEM,3,18  UMv.{iEj  
FITEM,3,20  5
]~'_V  
FITEM,3,22  c>,KZ!  
FITEM,3,24  {SOr#{1z*  
FITEM,3,26  ~3^ 8>d/  
FITEM,3,28  8Pfb~&X^Ws  
ASBL,       3,P51X !]42^?GH  
NUMCMP,ALL                                                    9FoHD  
NUMMRG,ALL,,,,LOW A8x
vo/n$  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 LU \i0|i|  
MP,EX,1,3.0E6                                    s kvGU(G}  
MP,PRXY,1,0.33                                                E{`kaWmC&~  
MP,DENS,1,1900.0 D!@Ciw  
TB,DP,1                                                       @#<D ^"  
TBDATA,1,2.4E4,13.6,0 ]"r&]qx7  
B#exHf8  
MP,EX,2,3.0E6                                    q6@Lp^f  
MP,PRXY,2,0.33                                                tI]Q%S,  
MP,DENS,2,1900.1 \j+1V1t9  
TB,DP,2                                                       0\H\lKcK  
TBDATA,1,2.4E4,13.6,0 *RDn0d[  
d+<G1w&z  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 wYdb*"R  
type,1 tW.9yII  
MAT,1  bNm#tmSt  
lsel,S,line,,1,8,1 u&^b~#T  
lsel,A,line,,15,27,4     ICpAt~3[M  
Lesize,ALL,,,4,,,,,0 UG'Q]S#!  
lsel,S,line,,29,32,1     }bs+-K  
Lesize,ALL,,,3,,,,,0 QaS1
Dh  
  r~E=4oB7  
lsel,S,line,,9,13,2 x%s-
+&  
lsel,A,line,,18,20,2 XywE1}3  
lsel,A,line,,21 O $LfuL  
lsel,A,line,,26,28,2     \<.+rqa!  
Lesize,ALL,,,10,,,,,0 1!ii;s^e  
lsel,S,line,,10 rTM}})81  
lsel,A,line,,22 R"4Vtww  
lsel,A,line,,24,25,1 hmvfw:Nq4  
Lesize,ALL,,,20-H-2*R,,,,,0 z[LNf.)}  
lsel,S,line,,12,14,2 kC WEtbz1  
lsel,A,line,,16,17,1 5rwu!Y;7*  
Lesize,ALL,,,H-2*R,,,,,0 &[_g6OL  
asel,s,area,,2,13,1  -]L6=  
amesh,all Jk&3%^P{m  
   P !AEf#1  
type,1 neB\q[
k  
MAT,2 3("_Z%  
asel,s,area,,1  grS:j+_M2m  
amesh,all f6EZ( v  
nummrg,all,,,,low y.anl  
numcmp,all W~6
EEyD%  
@u.58H& }R  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 A]<y:^2])C  
/solu WeJl4wF  
antype,static !4]TXH0f  
deltim,0.1,0.05,0.2 ` w=>I  
autots,on !使用自动时间步 O80<Z#%j`  
pred,on !打开时间步长预测器 Z"Zmo>cV4  
lnsrch,on !打开线性搜索 @>u]4Jn  
nlgeom,on !打开大位移效果 3Ko/{f  
nropt,full !设定牛顿-拉普森选项 ;}7Rjl#  
cnvtol,f,,0.02,0 YO#M/%^j  
Nsel,s,loc,x,-10 -mAUo;O  
Nsel,a,loc,x,10 =w;F<M|Y  
D,all,,,,,,UX,,,,, Q8C_9r/:N>  
Nsel,s,loc,y,H-20 :Uz|3gq  
D,all,,,,,,Uy,,,,, WM Fb4SUR  
time,1 gRrL[z  
acel,,10  C`K?7v3$m  
nsel,all |^0XYBxQ  
esel,all nv GF2(;l  
allsel,all H]P. x!I  
solve 4<9=5
q]  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 J cPtwa;q@  
time,2 IW>T}@ |  
esel,s,type,,2 *,3SGcYdJj  
ekill,all ;t'5},(FP  
nsel,all +GG9^:<yr  
esel,all ,
qA(\[  
allsel,all ;>#wU'  
solve 
^.1)};i  
fini < nXL  
yGI;ye'U  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 4^:\0UF  
!/POST1   u0 P|0\  
!LCDEF,1,1,, 4Z1ST;  
!LCDEF,2,2,, bmJ5MF]_fG  
!LCASE,2, a<@1-j<  
!LCOPER,SUB,1, ,2, 3%l*N&gsg:  
!Nsel,s,loc,Y,H 2!0c4a^z  
!PLVECT,U, , , ,VECT,ELEM,ON,0    dpJ_r>NI  
!PLNSOL, U,Y, 0,1.0  ;ZH3{  
m/Oh\KlIl  
0f~C#/[t7  
4 kn|^  
:a^t3s  

看看谢谢楼主

观望中。。。。。

学习学习