在做隧道ANSYS简单开挖模拟时，出现了一些问题，望各位前辈给多多指教！ NWq [22X
|
O>{t}6o
隧道洞泾假定为1.0m（实际上不可能这么小，仅仅是数值实验），隧洞轴线左右及下部各取10m为边界，施加法向约束。 \n( 'KVbf
-bG#h)yj
通过APDL命令，改变隧道的轴线埋深，采取全断面一次性开挖，整个荷载步分为两步：第一步计算初始应力场，第二步是杀死 YjMbd ?v
a-o
hS= W
开挖体单元，计算开挖后引起的地表沉降（位移场，第二步减去第一步），计算结果发现： \p_8YC
/-cX(z
7
隧道的轴线埋深越浅，计算得到的地表沉降越小，这与工程实际不符，请问原因何在？ yM\tbT/l
lqwJ F &
以下是我的计算命令流，希望得到高手的帮助啊！ ce-m )o/
(,Zz&3
AV
+=lcN~U2
Fini YQw/[
/clear n$Oky-P"
/title,TUNNEL ANALYSIS HEW9YC"
/prep7 s`7
_J9
ET,1,PLANE2 o cuNrkZ
keyopt,1,3,2 7Ohu$5\
)l81R
*SET,R,0.7 88X*:Kf?:
*SET,H,10.0 ps4W wk(
*SET,X1,-2*R hwb(W?*
*SET,Y1,-2*R /m|&nl8"qe
*SET,W1,4*R _s#/f5<:B
*SET,H1,4*R -p]`(S%
smQpIB;
*SET,X2,-10.0 `
TVcI\W
*SET,Y2,H-20.0 \$0
x8B
*SET,W2,20.0 E6d8z=X(
*SET,H2,20.0 DqC }f#
CYL4,0,0,0,45,R,405,0 APOU&Wd
BLC4,X1,Y1,W1,H1 =M
km: '1r
BLC4,X2,Y2,W2,H2 5 <>agK]
AOVL,1,2,3 cbHn\m)J,
NUMCMP,ALL }sN9QgE
NUMMRG,ALL,,,,LOW Bo\dt@0;
LANG,11,8,90,, . I9] `Q
LANG,11,7,90,, f ;;
S
LANG,10,7,90,, , $F0D
LANG,10,6,90,, bT6)(lm
LANG, 9,6,90,, I3y4O^?
LANG, 9,5,90,, f 4 ;8?
LANG,12,5,90,, i(qYyO'
LANG,12,8,90,, JV*,!5
L,1,7 uL{~(?U $
L,2,8 B*~Bm.
L,3,5 1VM5W!}
L,4,6 d+| ! 6
NUMCMP,ALL 7bW!u*v-c
NUMMRG,ALL,,,,LOW >"My\o
FLST,3,4,4,ORDE,2 g(F*Y> hk
FITEM,3,29 C36.UZoc
FITEM,3,-32 k=/|?%
ASBL,2,P51X )jZ=/ xG
FLST,3,8,4,ORDE,8 2n r
UE
FITEM,3,14 #?d#s19s
FITEM,3,16 >Ti2E+}[M
FITEM,3,18 s 7cyo
]
FITEM,3,20 Y8xnvK*
FITEM,3,22 5x'y{S<
FITEM,3,24 X!+ a;wr
FITEM,3,26 P!&CH4+
FITEM,3,28 vrldRn'*9
ASBL, 3,P51X j24
NUMCMP,ALL ' ah0IYe
NUMMRG,ALL,,,,LOW .*ovIU8
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 }4ju2K
MP,EX,1,3.0E6 l , 2z5p
MP,PRXY,1,0.33 sN) .Jo
MP,DENS,1,1900.0 P
T.jR *
TB,DP,1 \ @t5S
TBDATA,1,2.4E4,13.6,0 ; S M^
a~TZ9yg+HL
MP,EX,2,3.0E6 ~"YNG? Rre
MP,PRXY,2,0.33 /'_ RI
MP,DENS,2,1900.1 swgBPJ"?
TB,DP,2 UN 4) >\Y
TBDATA,1,2.4E4,13.6,0 D}U<7=\3H
nY0UnlB`
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 D;n%sRq(Z
type,1 e[dRHl
MAT,1 ymCIk /\
lsel,S,line,,1,8,1 twS3J)UH
lsel,A,line,,15,27,4 w[G_ w:$a
Lesize,ALL,,,4,,,,,0 9 (.9l\h
lsel,S,line,,29,32,1
>U/g*[>
Lesize,ALL,,,3,,,,,0 cb,sb^-
$NRb'
lsel,S,line,,9,13,2 q,&T$Tw
lsel,A,line,,18,20,2 ( t['
lsel,A,line,,21 4^^rOi0
lsel,A,line,,26,28,2 0l@+xS;
Lesize,ALL,,,10,,,,,0 .4(f0RG
lsel,S,line,,10 >D #}B1(!
lsel,A,line,,22 W>nb9Isp
lsel,A,line,,24,25,1 :vT%5CQ
Lesize,ALL,,,20-H-2*R,,,,,0 K7VG\Ec
lsel,S,line,,12,14,2 04d$_1:}a
lsel,A,line,,16,17,1 % "^XxVJ*
Lesize,ALL,,,H-2*R,,,,,0 9r#{s Y
asel,s,area,,2,13,1 ^e RT 8I
amesh,all n "Vd"}sU.
>MrU^t
type,1 c/U6K
yiK
MAT,2 (W!$6+GT
asel,s,area,,1 Ny LnE
amesh,all !17Z\Ltqyj
nummrg,all,,,,low c `; LF'!
numcmp,all mK4|=Q
mh]$g<*m
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 =7,U qMl_
/solu uJ\ Nga<?
antype,static J,Ap9HJt
deltim,0.1,0.05,0.2 `ayc YoD
autots,on !使用自动时间步 .Tv(1HAc2l
pred,on !打开时间步长预测器 dhkpkt<G8
lnsrch,on !打开线性搜索 N
Z`hy>LF^
nlgeom,on !打开大位移效果 ,+9r/}K]/
nropt,full !设定牛顿-拉普森选项 < z)G& h@
cnvtol,f,,0.02,0 RkwY3 s"
Nsel,s,loc,x,-10 E%B Gf}h
Nsel,a,loc,x,10 V+z)B+
D,all,,,,,,UX,,,,,
+Bfi/ >
Nsel,s,loc,y,H-20 v2sU$M
D,all,,,,,,Uy,,,,, [+z:^a1?V
time,1 W
F<V2o{k
acel,,10 Y_sV e
nsel,all T%/w^27E
esel,all zmRK%a(
allsel,all c"1d#8 J
solve (=rv `1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 ;i>E @
time,2 koZp~W-
esel,s,type,,2 edx'p`%d5
ekill,all ^mQfXfuL
nsel,all qw1J{xoHW
esel,all q2vz #\A?
allsel,all xS1|Z|&
solve 4&fnu/,Z
fini [hbp# I~*[
2zu~#qU[)M
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 t@M] ec
!/POST1 J7o?h9
!LCDEF,1,1,, cU|jT8Q4H
!LCDEF,2,2,, 12dW:#[
!LCASE,2, x$DJ
!LCOPER,SUB,1, ,2, eCD,[At/
!Nsel,s,loc,Y,H +to9].O7y
!PLVECT,U, , , ,VECT,ELEM,ON,0 b%j :-^0V
!PLNSOL, U,Y, 0,1.0 vy2aNUmt
sb.SpF>
)W|jt/
;(6lN<i U
%$| k3[4V
