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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ Ce~ a(J|"  
.R!
/?eN  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 L32[IL|  
2z+-vT%  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 miv)R  
{CH *?|t  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： ~IIlCmMl,  
p*l]I*x'<  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ bxFDB^  
KvtX>3#qM  
以下是我的计算命令流，希望得到高手的帮助啊！ CgxGvM4  
_1EWmHZ?  
cE SSSH!m  
Fini SGK=WLGM
8  
/clear ){S/h<4m  
/title,TUNNEL ANALYSIS W{js9$oJ  
/prep7 ^`< %P
k  
ET,1,PLANE2 6O9?":3;  
keyopt,1,3,2
+GDT@,/  
x}(p\Efx  
*SET,R,0.7 ,LpGE>s  
*SET,H,10.0 KH_~DZU*5  
*SET,X1,-2*R                                                    Z<M?_<3  
*SET,Y1,-2*R                                                    3<vw#]yL  
*SET,W1,4*R                                                   9EIOa/*  
*SET,H1,4*R                                                   klON6<w  
                                                              =bDy :yY}  
*SET,X2,-10.0                                                   R3\oLT4  
*SET,Y2,H-20.0                                                    2C@s-`b   
*SET,W2,20.0                                                  10`]&v]T  
*SET,H2,20.0 {L9WeosQ  
CYL4,0,0,0,45,R,405,0 J3P)oM[  
BLC4,X1,Y1,W1,H1                                              v$_YZm{!<  
BLC4,X2,Y2,W2,H2 %V2A}78  
AOVL,1,2,3                                              fBZAO  
NUMCMP,ALL                                                    NL7CeHs5  
NUMMRG,ALL,,,,LOW 6'YsSde".  
LANG,11,8,90,, yWkg4  
LANG,11,7,90,,  p;k7\7  
LANG,10,7,90,, co-dq\P  
LANG,10,6,90,, @GrQ/F7  
LANG, 9,6,90,, }k7t#O  
LANG, 9,5,90,, nJ{vO{N  
LANG,12,5,90,, !caY  
LANG,12,8,90,, y^xEZD1X6-  
L,1,7 {=)g?!zC  
L,2,8 5q Y+^jO]o  
L,3,5 }63Qh}_Y  
L,4,6 0S}ogU[k  
NUMCMP,ALL                                                    ^w1&A3=6  
NUMMRG,ALL,,,,LOW l*(Ml= O{  
FLST,3,4,4,ORDE,2    WPnw  
FITEM,3,29  $~s|%>@  
FITEM,3,-32 {~u#.
(  
ASBL,2,P51X  xLp<G(;  
FLST,3,8,4,ORDE,8  
dvLO#o{  
FITEM,3,14  5l]G1+  
FITEM,3,16  gZ b+m  
FITEM,3,18  Y(D&JKx  
FITEM,3,20  ljJi|+^$  
FITEM,3,22  "jFRGgd79  
FITEM,3,24  riuG,$
EX  
FITEM,3,26  %40uw3  
FITEM,3,28  %n7Y5|Uh  
ASBL,       3,P51X S]#xG+$<  
NUMCMP,ALL                                                    T)Ohk(jK1  
NUMMRG,ALL,,,,LOW |)'6U3  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 VP\'p1a  
MP,EX,1,3.0E6                                    w,/&oe5M+  
MP,PRXY,1,0.33                                                m5w ZS>@  
MP,DENS,1,1900.0 dH8^\s .F  
TB,DP,1                                                       ,k=1'7d  
TBDATA,1,2.4E4,13.6,0 Yc]  
dYP-QUM$7  
MP,EX,2,3.0E6                                    J#OiY  
MP,PRXY,2,0.33                                                miCW(mbO8  
MP,DENS,2,1900.1 g~#HiBgWq[  
TB,DP,2                                                       iJH;OV;P  
TBDATA,1,2.4E4,13.6,0 PRo;NE  
h0v4!`PQ-  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 U!xOJ  
type,1 AR}q<k6E  
MAT,1  E?@batIrf  
lsel,S,line,,1,8,1 {TV6eV  
lsel,A,line,,15,27,4     \8 ~`NF  
Lesize,ALL,,,4,,,,,0 } eL*gy  
lsel,S,line,,29,32,1     `EjPy>kM  
Lesize,ALL,,,3,,,,,0 }%m:^*@$9  
  iR`c/  
lsel,S,line,,9,13,2 !~Am1\02  
lsel,A,line,,18,20,2 $: %U`46%s  
lsel,A,line,,21 UCv9G/$  
lsel,A,line,,26,28,2     luYa+E0  
Lesize,ALL,,,10,,,,,0 )X dpzWod  
lsel,S,line,,10 rElbzL"&<  
lsel,A,line,,22 Q8^fgI|  
lsel,A,line,,24,25,1 xbeVqP  
Lesize,ALL,,,20-H-2*R,,,,,0 x<-n}VK\  
lsel,S,line,,12,14,2 kq m$a  
lsel,A,line,,16,17,1 ?cF-w!>o8  
Lesize,ALL,,,H-2*R,,,,,0 6,D)o/_  
asel,s,area,,2,13,1  4jyDM68i  
amesh,all N8s2v W  
   j!oD9&W4~  
type,1 ls<7
Qe"a  
MAT,2 ^9g+\W  
asel,s,area,,1  1E]|>)$  
amesh,all j|N8"8"  
nummrg,all,,,,low tk"+PTGJT  
numcmp,all ~|7jz;$V  
[^U#ic>cT  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 3*DwXH+  
/solu QnWM<6xK"  
antype,static b)5z
'zQu  
deltim,0.1,0.05,0.2 kCq]#e~wq  
autots,on !使用自动时间步 mM&*_#( 6  
pred,on !打开时间步长预测器 "HuV'  
lnsrch,on !打开线性搜索 HhmC
+3w.7  
nlgeom,on !打开大位移效果 | Q Y_ci  
nropt,full !设定牛顿-拉普森选项 0NN{2"M$p  
cnvtol,f,,0.02,0 tPT\uD#t  
Nsel,s,loc,x,-10 XYKWOrkQqa  
Nsel,a,loc,x,10 m]1=o7  
D,all,,,,,,UX,,,,, :uAW  
Nsel,s,loc,y,H-20 0]F'k8yLN  
D,all,,,,,,Uy,,,,, p8>.Q/4  
time,1
,C#Mf@b  
acel,,10  MZ(TST"  
nsel,all N}Ol`@@#h  
esel,all v}=3  
allsel,all :J3ZTy
jb  
solve $GzTDq Y9@  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 "%c\i-&t  
time,2 0L34)W  
esel,s,type,,2 6AZJ,Q\E@  
ekill,all {[~ !6&2(k  
nsel,all idm!6]  
esel,all <p48?+K9  
allsel,all 0hKF)b  
solve t|gEMDGa3  
fini d9zI A6y  
wAr (5nEbx  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 dQ=mg#(  
!/POST1   Xd_86q8o  
!LCDEF,1,1,, C043h?x  
!LCDEF,2,2,, >IW0YIQy,  
!LCASE,2, /t-m/&>  
!LCOPER,SUB,1, ,2, SR9Cl  
!Nsel,s,loc,Y,H 8oSndfV  
!PLVECT,U, , , ,VECT,ELEM,ON,0    G>?hojvi  
!PLNSOL, U,Y, 0,1.0 tHS
e>*eC  
>EFjyhVE  
Do7&OBI~  
&zsaVm8  
bk E4{P"  

学习学习

观望中。。。。。

看看谢谢楼主