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

在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ ^VC7C~NZ!M  
@?j@yRe  
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 "sSjVu  
lxCX-a`@p  
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 M@5KoMsB9  
,Os7T 1>  
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： Yo@m50s$  
qU}lGf!dVn  
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ /H.(d 4C  
 }Olr  
以下是我的计算命令流，希望得到高手的帮助啊！ &s2#1  
dqL-'  
w64.R
4e  
Fini a$laRtId7  
/clear 3a/[."W u  
/title,TUNNEL ANALYSIS *u<rU,C8  
/prep7 #?A]v>I;C  
ET,1,PLANE2 }=xI3;7  
keyopt,1,3,2 nP5d?  
p9k4w% ~:  
*SET,R,0.7 $:PF9pY(  
*SET,H,10.0 ./';P<)  
*SET,X1,-2*R                                                    2z[r@}3  
*SET,Y1,-2*R                                                    ^#}dPGm  
*SET,W1,4*R                                                   u~ VswXc4  
*SET,H1,4*R                                                   .Kg|f~InO  
                                                              @'@s*9Nr  
*SET,X2,-10.0                                                   XTIu(f|d_;  
*SET,Y2,H-20.0                                                    J&n ^y  
*SET,W2,20.0                                                  Y>(ZsHu  
*SET,H2,20.0 mL8A2>Gig  
CYL4,0,0,0,45,R,405,0 xcAF  
BLC4,X1,Y1,W1,H1                                              r94j+$7  
BLC4,X2,Y2,W2,H2 .jLMl*6%:  
AOVL,1,2,3                                              %_]=i@Y~  
NUMCMP,ALL                                                    [vZfH!vLP  
NUMMRG,ALL,,,,LOW u'#`yTB6b  
LANG,11,8,90,, JTSq{NN  
LANG,11,7,90,, rsd2v9  
LANG,10,7,90,, M-}j9,oR`  
LANG,10,6,90,, 782[yLyv  
LANG, 9,6,90,, C+X)">/+L  
LANG, 9,5,90,, 2Px$0&VN  
LANG,12,5,90,, XhQw+j~1.  
LANG,12,8,90,, 9k714bnMLX  
L,1,7 L7i}Ga!
8  
L,2,8 J
slk  
L,3,5 200L  
L,4,6 HVtr,jg  
NUMCMP,ALL                                                    aW52.X z%8  
NUMMRG,ALL,,,,LOW bbfDt^  
FLST,3,4,4,ORDE,2   N |OMj%Uk  
FITEM,3,29  5HWwl.
D  
FITEM,3,-32 ?Ry%c6(}  
ASBL,2,P51X  :;" aUHU'  
FLST,3,8,4,ORDE,8   SJ1w1^#Pz  
FITEM,3,14  e6tH/`Uln  
FITEM,3,16  TBoM{s=.  
FITEM,3,18  n25irCD`  
FITEM,3,20  ORV}j,Ym  
FITEM,3,22  5KL9$J9k  
FITEM,3,24  Vy?R/ Uu  
FITEM,3,26  r7
8u=r  
FITEM,3,28  /}h71V!  
ASBL,       3,P51X d,h~u{  
NUMCMP,ALL                                                    {bAWc.  
NUMMRG,ALL,,,,LOW B2ek&<I7N  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 g:G%Ei~sF  
MP,EX,1,3.0E6                                    &4w\6IR  
MP,PRXY,1,0.33                                                cq'}2pob  
MP,DENS,1,1900.0 _A@fP[C  
TB,DP,1                                                       { r<(t#  
TBDATA,1,2.4E4,13.6,0 Q0 uP8I}n  
~zXG<}n  
MP,EX,2,3.0E6                                    [4K9|/J  
MP,PRXY,2,0.33                                                T% /xti5$!  
MP,DENS,2,1900.1 h|(ZXCH  
TB,DP,2                                                       dv1x78xG>  
TBDATA,1,2.4E4,13.6,0  VJ~X#Q  
t 1Ir4  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 C8bGae(  
type,1 aTL8l.c2  
MAT,1  b0~H>cnA
 
lsel,S,line,,1,8,1 Dz"u
8 f  
lsel,A,line,,15,27,4     E<X{72fb>  
Lesize,ALL,,,4,,,,,0
X/- W8  
lsel,S,line,,29,32,1     +d6Aw}*
 
Lesize,ALL,,,3,,,,,0 s2(w#n)  
  STFQ";z$  
lsel,S,line,,9,13,2 cJ7{4YK_#/  
lsel,A,line,,18,20,2 Mp^OL7p^^  
lsel,A,line,,21 . Ce&9l  
lsel,A,line,,26,28,2     Yj^avO=;  
Lesize,ALL,,,10,,,,,0 AHp830\  
lsel,S,line,,10 :{TmR3.  
lsel,A,line,,22 /ZIJ<#o[
 
lsel,A,line,,24,25,1 '[M^f+H|  
Lesize,ALL,,,20-H-2*R,,,,,0 xN{"%>Mx  
lsel,S,line,,12,14,2 Dr8WV\4@  
lsel,A,line,,16,17,1 2yEO=SN,(  
Lesize,ALL,,,H-2*R,,,,,0 RpmBP[  
asel,s,area,,2,13,1  A-Q{*{^#  
amesh,all 8xD<A|  
   4."o.:8x  
type,1 !i~x"1  
MAT,2 CN8@c!mB  
asel,s,area,,1  z^a6%N  
amesh,all
D^N#E>,  
nummrg,all,,,,low i&}zcGC  
numcmp,all g "K#&  
lgv-)5|O+H  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 Th&Wq  
/solu niBjq#bJi  
antype,static |5<&r]xN  
deltim,0.1,0.05,0.2 /A[oj2un  
autots,on !使用自动时间步 *D09P%  
pred,on !打开时间步长预测器 7T|J[WO  
lnsrch,on !打开线性搜索 q5'yD;[hE  
nlgeom,on !打开大位移效果 Y?&DEKFbD  
nropt,full !设定牛顿-拉普森选项 U'Ja\Ek/f  
cnvtol,f,,0.02,0 XL.CJ5y>  
Nsel,s,loc,x,-10 l~Ie#vak  
Nsel,a,loc,x,10 ;mo\ yW1  
D,all,,,,,,UX,,,,, Wd^F%)(  
Nsel,s,loc,y,H-20 PcN
fTB{  
D,all,,,,,,Uy,,,,, f'zU^/$rf  
time,1 fzkCI  
acel,,10  @sgT[P*ut  
nsel,all #2lvfR|  
esel,all T$.-{I  
allsel,all O%fp;Y{`  
solve 9znx1AsN  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 w.T=Lzp  
time,2 .j:.WnW  
esel,s,type,,2 q
j3bt_F!x  
ekill,all Kl46CZs#8  
nsel,all ~f]r>jQM  
esel,all uJ)=+Exii
 
allsel,all 8[`^(O#\E  
solve +/~\b/  
fini ${I@YSU  
HgSmAziv  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 f,`}hFD  
!/POST1   Z$#ZYD  
!LCDEF,1,1,, eMm~7\ R  
!LCDEF,2,2,, zI2KIXcc  
!LCASE,2, @7HOL-i  
!LCOPER,SUB,1, ,2, VM\Z<}C  
!Nsel,s,loc,Y,H f5=t*9_-[  
!PLVECT,U, , , ,VECT,ELEM,ON,0     wc+N  
!PLNSOL, U,Y, 0,1.0 T956L'.+G  
Dt~}9HrU  
x@ s`;qz  
PkCeV]`w  
Hdq/E>
u  

看看谢谢楼主

观望中。。。。。

学习学习