小弟做的是三维隧道核心土开挖,在算地应力就遇到了问题,求高手解决下面是命令流 T'9I&h%\
finish$/clear$/filname,rtest2$/prep7 \\F^uM7,
et,1,mesh200,2 !3-D线单元2节点 :USN`"
et,2,mesh200,6 !3-D面单元4节点 oJC-?
et,3,SHELL181 !用于模拟支护结构的壳单元 _E eH
et,4,SOLID45 !用于模拟围岩的三维实体单元 m(QGP\Ya
r,1,0.4,0.4,0.4,0.4 !壳单元的厚度,单位 \u@4eBAV
r,2,0.6,0.6,0.6,0.6 :0,q>w
mp,ex,1,30e9 !支护结构材料属性,弹性模型,单位Pa ~_WsjD0O
mp,prxy,1,0.2 !泊松比 7Iu^l4=2
mp,dens,1,2500 !密度 WF0%zxg ]
mptemp,,,,,,,, !围岩(清除前一次的温度设定 ) %VG;vW\V
mptemp,1,0 !开始此次的设定 CZB!vh0
mp,ex,2,0.022e9 !围岩材料属性 d (Ufj|;
mp,prxy,2,0.38 !泊松比,无单位 BrYU*aPW;
mp,dens,2,2000 (9:MIP
tb,dp,2 !dp模型 ,4oYKJ$+h
tbmodif,1,1,0.035e6 !输入凝聚力 6@pPaq6
tbmodif,1,2,22 !输入摩擦角 ^)ouL25Z*2
mptemp,,,,,,,, L3:dANG
mptemp,1,0 Jlw<%}r
mp,ex,3,0.022e9 !围岩材料属性 1qWIku
mp,prxy,3,0.38 !泊松比,无单位 9{{QdN8
mp,dens,3,2001 K*;e>{p
tb,dp,3 !dp模型 DQ7+
tbmodif,1,1,0.035e6 !输入凝聚力 :.kc1_veYS
tbmodif,1,2,22 !输入摩擦角 Va^Y3/
mptemp,,,,,,,, dz
fR ^Gv
mptemp,1,0 T^LpoN/T
mp,ex,4,0.022e9 !围岩材料属性 f8`K8Y]4
mp,prxy,4,0.38 !泊松比,无单位 ^'EeJN
mp,dens,4,2001 XMb]&VvH
tb,dp,4 ,"?h_NbF
tbmodif,1,1,0.035e6 :uhU<H<,f
tbmodif,1,2,22 R;,&s!\<
mptemp,,,,,,,, 1e`/N+6u
mptemp,1,0 n4Fh*d ixg
mp,ex,5,0.022e9 !围岩材料属性 x`8rR;N!
mp,prxy,5,0.38 !泊松比,无单位 8A/;a{
mp,dens,5,2003 ( ONn{12Q
tb,dp,5 Wyu$J
tbmodif,1,1,0.035e6 L[\m{gN
tbmodif,1,2,22 0r?975@A
mptemp,,,,,,,, n1OxT"tD
mptemp,1,0 Oo'IeXQ9(
mp,ex,6,0.022e9 !围岩材料属性 4|Z;EAFx
mp,prxy,6,0.38 !泊松比,无单位 n%;wQ^
mp,dens,6,2004 %H=d_Nm{
tb,dp,6 c$?(zt;
tbmodif,1,1,0.035e6 C?@vBM}
tbmodif,1,2,22 utIR\e#:B
mptemp,,,,,,,, 9(?9yFbj5
mptemp,1,0 pvL)BD
mp,ex,7,0.022e9 !围岩材料属性 Cz=HxU80J
mp,prxy,7,0.38 !泊松比,无单位 )N[9r{3
mp,dens,6,2005 )cA#2mlS'1
tb,dp,7 _@I8B
tbmodif,1,1,0.035e6 :3`6P:^
tbmodif,1,2,22 +2%ih!
mptemp,,,,,,,, uq<kT [
mptemp,1,0 A4RA5N/}
mp,ex,8,0.022e9 !围岩材料属性 v"M5';ZS>
mp,prxy,8,0.38 !泊松比,无单位 XWH{+c"
mp,dens,8,2006 0TA{E-A
tb,dp,8 Il(p!l<Xz#
tbmodif,1,1,0.035e6 DBDHe-1[+
tbmodif,1,2,22 om%L>zfB
mptemp,,,,,,,, &YQ
mptemp,1,0 hE>ux"_2/
mp,ex,9,20e9 !支护结构材料属性,弹性模型,单位Pa 40TS=evG
mp,prxy,9,0.2 y<7C!E#b8
mp,dens,9,2000 tfzIem
NUMSTR,KP,1 !设置起始关键点的编号为1 Ay7I_"%
NUMSTR,LINE,1 xWk:7 ,/
NUMSTR,AREA,1 qfp,5@p
NUMSTR,VOLU,1 %:I\M)t}k
CSYS,0 !激活笛卡尔坐标系 b&:>v9U
k,1,0,0 <<#-IsT
k,2,0,11.37 4yhan/zA
k,3,3.34,0 _'9("m V
k,4,-3.34,0 ^LfN6{
k,101,0,0,10 ^):m^w.
k,102,0,11.37,10 (_S`9Z8=
k,103,3.34,0,10 $hexJzX
k,104,-3.34,0,10 x]
[/9e
circle,1,5.8,101 ~B!O
X
circle,2,14.4,102 u6o:~=WwM
circle,3,2.67,103 ; MU8@?yN
circle,4,2.67,104 RlH|G
lcsl,all C[f'1O7
lsel,s,,,17,22,1 *?|LE
C
lsel,a,,,46,48,1 Xup rl2+
lsel,a,,,27,28,1 V1Yab#
lsel,a,,,43 R=uzm=&nR
lsel,inve :1h1+b@,
ldele,all,,,1 -2NXQ+m ;
lsel,all @Qw~z0PE<l
lcomb,17,18 *^j'G^n
lcomb,19,21 ^(<Ecdz(
lcomb,27,43 R `}C/'Ty
lcomb,28,46 e~#;ux
cm,zh,line 7_Yxz$m
circle,1,8.8,101 -#0(Jm'
circle,2,17.4,102 Xv[5)4N
circle,3,5.67,103 @c&}\#;
circle,4,5.67,104 eF]`?AeWQ
cmsel,s,zh,line ovv<7`
lsel,inve yWI30hW
lcsl,all %LBa;M
lsel,s,,,23,25,2 !u@XEN>/
lsel,a,,,18,21,3 S/YT
V
lsel,a,,,26,62,36 3IXai)6U
lsel,a,,,24,59,35 j#^EZ/
lsel,a,,,33,55,22 k
I{)"
lsel,a,,,34,58,24 O$QtZE61
cmsel,a,zh,line l,cnMr^.W
lsel,inve kev|AU (WX
ldele,all,,,1 ks92-%;:
lsel,all 6H+'ezM
lcomb,23,25 z~[:@mGl
lcomb,18,21 Rf *we+
lcomb,33,55 4 .7YIM
lcomb,34,58 RTN?[`
allsel npsDy&