在做隧道ANSYS简单开挖模拟时,出现了一些问题,望各位前辈给多多指教! (B\Kb4m
.x$V~t
隧道洞泾假定为1.0m(实际上不可能这么小,仅仅是数值实验),隧洞轴线左右及下部各取10m为边界,施加法向约束。 *`OXgkQ
v_5O*F7)
通过APDL命令,改变隧道的轴线埋深,采取全断面一次性开挖,整个荷载步分为两步:第一步计算初始应力场,第二步是杀死 0084`&Ki
86%k2~L
开挖体单元,计算开挖后引起的地表沉降(位移场,第二步减去第一步),计算结果发现: /;_$:`|/
<2*+Y|Lk2
隧道的轴线埋深越浅,计算得到的地表沉降越小,这与工程实际不符,请问原因何在? mTH[*Y,
>%"TrAt
以下是我的计算命令流,希望得到高手的帮助啊! 0uKm)t/
6 w!qZ4$
z .lb(xQ
Fini AYNz {9
/clear `LFT"qnp
/title,TUNNEL ANALYSIS Il/`#b@h
/prep7 ^/2O_C
ET,1,PLANE2 9on$0
keyopt,1,3,2 qw4wg9w5p
o
^w^dgJ
*SET,R,0.7 P3bRv^
*SET,H,10.0 (q"S0{
*SET,X1,-2*R U;&s=M0[
*SET,Y1,-2*R ~i-n_7 +
*SET,W1,4*R <mLU-'c@
*SET,H1,4*R F\+9u$=
\ bNDeA&l
*SET,X2,-10.0 jhG6,;1zMI
*SET,Y2,H-20.0 `iuQ.I
*SET,W2,20.0 (N}\Wft%
*SET,H2,20.0 El`G<esX
CYL4,0,0,0,45,R,405,0 +SV!QMIg
BLC4,X1,Y1,W1,H1 a&)$s;
BLC4,X2,Y2,W2,H2 6wiuNGZb
AOVL,1,2,3 2QHu8mFU
NUMCMP,ALL eJOo~HIWQ
NUMMRG,ALL,,,,LOW L #vk77
LANG,11,8,90,, OwC{ Ad{
LANG,11,7,90,, #SLiv
LANG,10,7,90,, ;XY#Jl>tg
LANG,10,6,90,, oz'jt} ?
LANG, 9,6,90,, %|}7YH41
LANG, 9,5,90,, MMRO@MdfV
LANG,12,5,90,, "=f,4Zbj
LANG,12,8,90,, ORo +]9)Yv
L,1,7 1\f8-:C
L,2,8 Sr10ot&ox
L,3,5 :v-,-3AG
L,4,6 G* mLb1
NUMCMP,ALL ^sZHy4-yK#
NUMMRG,ALL,,,,LOW .@(MNq{"6
FLST,3,4,4,ORDE,2 :fG9p`
FITEM,3,29 wu0JXB%&^
FITEM,3,-32 kD+B8TrW
ASBL,2,P51X NLWj5K)1P
FLST,3,8,4,ORDE,8 h7\EN
FITEM,3,14 Y%rC\Ij/i
FITEM,3,16 >*w(YB]/$V
FITEM,3,18 +DT)7koA
FITEM,3,20 Qa16x<Xlm
FITEM,3,22 vP<8,XG
FITEM,3,24 h1_KZ[X
FITEM,3,26 wCr+/"t
FITEM,3,28 \ U`rF
ASBL, 3,P51X a67NWH
NUMCMP,ALL 9T;l*
NUMMRG,ALL,,,,LOW 8-vNXvl
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!初始参数 M\9p-%"L
MP,EX,1,3.0E6 A^7Y%
MP,PRXY,1,0.33 bO{wQ1)Z_
MP,DENS,1,1900.0 Yv;18j*<
TB,DP,1 c;7ekj
TBDATA,1,2.4E4,13.6,0 W^nG\"T^
Y1FP |
MP,EX,2,3.0E6 /Big^^u
MP,PRXY,2,0.33 _EZrZB
MP,DENS,2,1900.1 0/%VejZ'
TB,DP,2 _q\w9gN
TBDATA,1,2.4E4,13.6,0 {wfe!f
*j=
whdw%J
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!剖分单元 blB00
type,1 2}R)0][W
MAT,1 &@D,|kHk
lsel,S,line,,1,8,1 n|iO)L\9aB
lsel,A,line,,15,27,4 [^-DFq5@
Lesize,ALL,,,4,,,,,0 ddjaM/.E
lsel,S,line,,29,32,1 q:v&wb%
Lesize,ALL,,,3,,,,,0 uod&'g{N
0~XZ
lsel,S,line,,9,13,2 *XlnEHv
lsel,A,line,,18,20,2 aTY\mKk
lsel,A,line,,21 ygp NMq#?X
lsel,A,line,,26,28,2 {*PB+WGe
Lesize,ALL,,,10,,,,,0 -z4pI=
lsel,S,line,,10 S(*SUH
lsel,A,line,,22 J4ltHk.|
lsel,A,line,,24,25,1 /e}NZo{)g
Lesize,ALL,,,20-H-2*R,,,,,0 O}#h^AU-BS
lsel,S,line,,12,14,2 :>g*!hpb
lsel,A,line,,16,17,1 RoWGQney
Lesize,ALL,,,H-2*R,,,,,0 "h}miVArS
asel,s,area,,2,13,1 ;HeUD5Nt6F
amesh,all )W/mt[;
e'z[JG=
type,1 _6 |lw&o07
MAT,2 <J QvuC
asel,s,area,,1 T/&4lJ^2l^
amesh,all y\&`A:^[ A
nummrg,all,,,,low u>.qhtm[
numcmp,all 5}4r'P$m:
^J]~&.l
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算初始应力场 xeX Pc7JG
/solu _0<EbJ8Z
antype,static (TV ye4Z
deltim,0.1,0.05,0.2 7x^P 74
autots,on !使用自动时间步
um[nz
pred,on !打开时间步长预测器 O?9&6x
lnsrch,on !打开线性搜索 :)\<
nlgeom,on !打开大位移效果 KZw"?%H[
nropt,full !设定牛顿-拉普森选项 m. G}#/
cnvtol,f,,0.02,0 (W|Eg
Nsel,s,loc,x,-10 l(@UpV-
Nsel,a,loc,x,10 Vx^+Z,y&QP
D,all,,,,,,UX,,,,, ^U.8grA
Nsel,s,loc,y,H-20 :YU_ \EV
D,all,,,,,,Uy,,,,, COa"zg
time,1 #xS8
acel,,10 V\"x#uB
nsel,all &!#a^d+` 0
esel,all +tNu8M@xFo
allsel,all %klC&
_g~_
solve uf (`I
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开挖 shOQ/
time,2 T,,,+gPx
esel,s,type,,2 "3A.x1uQ
ekill,all :Uf\r
`a9
nsel,all Ax4nx!W,
esel,all V&E)4KBOs
allsel,all Yqy7__vm
solve 3EN?{T<yf
fini ~1Q$FgLk
t%wC~1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!地表沉降 wvum7K{tI
!/POST1 PYYOC"$
!LCDEF,1,1,, _
a|zvH
!LCDEF,2,2,, CfA^Xp@vc
!LCASE,2, C{<qc,!4
!LCOPER,SUB,1, ,2, {i)k# `
!Nsel,s,loc,Y,H hTZaI *
!PLVECT,U, , , ,VECT,ELEM,ON,0 y_:i'Ri.
!PLNSOL, U,Y, 0,1.0 o })k@-oL
e.%`
tK3J
Re**)3#gn
eDR4c%
C d)j%