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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ zdYjF|  
%A`+WYeuX  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 vt8By@]:  
#4PN"o@  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 $]d^-{|  
rH-23S  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： L_T5nD^D  
}OUtsh]y  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ ]/Pn EU[  
+whDU2 "  
以下是我的计算命令流，希望得到高手的帮助啊！ ,prf;|e?  
>Ry01G]_/h  
w>gYx(8b  
Fini 2ESo2  
/clear 3S{/>1Y  
/title,TUNNEL ANALYSIS .]^?<bG  
/prep7 wT@og|M  
ET,1,PLANE2 Dtk=[;"k2a  
keyopt,1,3,2 _XBd3JN@
 
UN
Yqft4  
*SET,R,0.7 5+vaE 2v  
*SET,H,10.0 AH^/V}9H  
*SET,X1,-2*R                                                    A"]YM'.  
*SET,Y1,-2*R                                                    p{_" bB  
*SET,W1,4*R                                                   Y4-t7UlS;  
*SET,H1,4*R                                                   Ac@VGT:9  
                                                              3BI1fXT4=j  
*SET,X2,-10.0                                                   7! Nsm  
*SET,Y2,H-20.0                                                    OXA7w.^  
*SET,W2,20.0                                                  %EH)&k  
*SET,H2,20.0 & 21%zPm  
CYL4,0,0,0,45,R,405,0 h`q1  
BLC4,X1,Y1,W1,H1                                              ]M=&+c>H~  
BLC4,X2,Y2,W2,H2 *@5@,=d  
AOVL,1,2,3                                              f|5co>Hk  
NUMCMP,ALL                                                    4Z*/WsCv  
NUMMRG,ALL,,,,LOW :;}P*T*PU  
LANG,11,8,90,, i5Ggf"![  
LANG,11,7,90,, ye&;(30Oq  
LANG,10,7,90,, T)/eeZ$  
LANG,10,6,90,, fhiM U8(&  
LANG, 9,6,90,, M3\AY30L  
LANG, 9,5,90,, ?s01@f#  
LANG,12,5,90,, Cdn J&N{  
LANG,12,8,90,, [y(MCf19  
L,1,7 [0!(xp^  
L,2,8 /uflpV|  
L,3,5 9[4xFE?|  
L,4,6 e'~3oqSvR  
NUMCMP,ALL                                                    WWY6ha  
NUMMRG,ALL,,,,LOW 7Q 3k7  
FLST,3,4,4,ORDE,2   ?<!|  
FITEM,3,29  )lkjqFQ(  
FITEM,3,-32 kiEa<-]  
ASBL,2,P51X  7CysfBF0g  
FLST,3,8,4,ORDE,8   i!Ba]n   
FITEM,3,14  6nn*]|7  
FITEM,3,16  t@(HF-4~=  
FITEM,3,18  4#D,?eA7  
FITEM,3,20  }BEB1Q}L  
FITEM,3,22  6ujWNf  
FITEM,3,24  \fOEqe*5SM  
FITEM,3,26  Rq-ZL{LR7  
FITEM,3,28  97Vtn4N3  
ASBL,       3,P51X c<~H(k'+c  
NUMCMP,ALL                                                    ).O)p9  
NUMMRG,ALL,,,,LOW *. t^MP  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 ?ub35NLa  
MP,EX,1,3.0E6                                    @7}W=HB  
MP,PRXY,1,0.33                                                4g/d
P^  
MP,DENS,1,1900.0 3]hWfj1m2  
TB,DP,1                                                       kJU2C=m@e2  
TBDATA,1,2.4E4,13.6,0 X}]-*T|a  
 7GGUV  
MP,EX,2,3.0E6                                    +@UV?"d  
MP,PRXY,2,0.33                                                ?dTD\)%A  
MP,DENS,2,1900.1
rv;3~'V  
TB,DP,2                                                       BtZyn7a  
TBDATA,1,2.4E4,13.6,0 'u658Tj  
crCJrN=  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 z?zL97H  
type,1 4-w{BZuS  
MAT,1  lZ0 =;I  
lsel,S,line,,1,8,1 `cO:<^%  
lsel,A,line,,15,27,4     iU-j"&L5  
Lesize,ALL,,,4,,,,,0 kAx4fE[c  
lsel,S,line,,29,32,1     bt *k.=p  
Lesize,ALL,,,3,,,,,0 Kw}'W 8`c  
  2%1hdA<  
lsel,S,line,,9,13,2 }]TxlSp!;  
lsel,A,line,,18,20,2 *Ex|9FCt$  
lsel,A,line,,21 L];b<*d  
lsel,A,line,,26,28,2     =a!=2VN9y  
Lesize,ALL,,,10,,,,,0 g0H[*"hj  
lsel,S,line,,10 8L XHk l  
lsel,A,line,,22 $>gFf}#C  
lsel,A,line,,24,25,1 $'TM0Yu,  
Lesize,ALL,,,20-H-2*R,,,,,0 w!CNRtM:~  
lsel,S,line,,12,14,2 mOSv9w#,  
lsel,A,line,,16,17,1 N!32 wJ  
Lesize,ALL,,,H-2*R,,,,,0 !k%#R4*>  
asel,s,area,,2,13,1  [lAp62i5  
amesh,all K,]=6Rj  
   Ru XC(qcq  
type,1
0V]s:S  
MAT,2 -di o5a  
asel,s,area,,1  !wNO8;(  
amesh,all 67TwPvh  
nummrg,all,,,,low D$N/FJ8|G  
numcmp,all ,Q,^3*HX9}  
$(9U@N9E  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 v<;Md-<  
/solu >7r!~+B"9'  
antype,static #g=XUZ/"  
deltim,0.1,0.05,0.2 l~.-e^p?  
autots,on !使用自动时间步 %S
I'BJ  
pred,on !打开时间步长预测器 bcR_E5x$  
lnsrch,on !打开线性搜索 "3hMq1NQ`g  
nlgeom,on !打开大位移效果 ] - .aL  
nropt,full !设定牛顿-拉普森选项 '|4!5)/K  
cnvtol,f,,0.02,0 *H122njH+T  
Nsel,s,loc,x,-10 SaCh 7 ^  
Nsel,a,loc,x,10 1}37Q&2  
D,all,,,,,,UX,,,,, R3!t$5HG  
Nsel,s,loc,y,H-20 C?Ucu]cW  
D,all,,,,,,Uy,,,,, 2oU_2P  
time,1 G`zm@QL  
acel,,10  ccnK#fn v  
nsel,all {[(h[MW#  
esel,all xpI wrJO  
allsel,all ,47qw0=C  
solve Y~Ifj,\  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 dd["dBIZ '  
time,2 Wf<LR3  
esel,s,type,,2 jIF |P-  
ekill,all |'.  
nsel,all w:l"\Tm  
esel,all 6Iw\c  
allsel,all e*!kZAf  
solve x:7IIvP  
fini <1pEwI~  
aP`P)3O6)1  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 +O5hH8<&b  
!/POST1   Q &t<Y^B  
!LCDEF,1,1,, ^J d r>@  
!LCDEF,2,2,, 2[CdZ(k]5  
!LCASE,2, PzGWff!*n  
!LCOPER,SUB,1, ,2, Gc!x|V;T  
!Nsel,s,loc,Y,H }-fl$j?9E  
!PLVECT,U, , , ,VECT,ELEM,ON,0    2-b6gc7  
!PLNSOL, U,Y, 0,1.0 X?$_Sd"G+5  

Sc   
4OX^(  
Z6p
UZ[j,  
fT{Yg /j  

看看谢谢楼主

观望中。。。。。

学习学习