后处理中,发现自重应力场很混乱,明显有问题,见图。 o8Vtxnkg
请问高人如何解决?万分感谢! ,WSK
'
命令如下: g S xK9P
finish 1henQiIO
/CLEAR,start 8s(?zK\
/fILEname,slope analysis golr,+LSo
/units,si QFt7L
/pnum,kp,1 /J!C2
/pnum,line,1 :ZG^`H/X1d
/number,0 \i_y(;
/PREP7 )NAC9:8!
*afun,deg =2!AK[KxX
!定义单元类型 xEdCGwgp#
et,1,solid45 /)fx(u#
!定义材料属性 9 f$S4O5
!边坡围岩材料1的属性 9ZFvN*Zf'
mp,ex,1,6e7 &Uq++f6
mp,nuxy,1,0.29 FdwT
mp,dens,1,23 k.=67L
tb,dp,1 ]AkHNgW
tbdata,1,3e4 :&:IZkO
tbdata,2,30 =Lx*TbsFYt
tbdata,3,0 5dOA^P@`,M
!边坡围岩材料2的属性 JA0$Fz
mp,ex,2,5e7 !9|)v7}
mp,nuxy,2,0.21 4bXAA9"
mp,dens,2,23.5 nn><
k"
tb,dp,2 ~3-2Iu^F
tbdata,1,5e4 [GKSQt{)
tbdata,2 .
+_IpygQ
tbdata,3,0 6&V4W"k
!边坡围岩材料3的属性 zWvG];fsN
mp,ex,3,7e8 kw)@[1U
mp,nuxy,3,0.34 `s74g0h
mp,dens,3,25.1 !#cZ!
tb,dp,3 s!S,;H
tbdata,3,5e5 k<W n
tbdata,2,40 9Br2}!Ny
tbdata,3,0 ^,S\-Uy9
!边坡围岩材料4的属性 St/<\Y,wr
mp,ex,4,1.5e9 r"#h6lYK&
mp,nuxy,4,0.3 Ij;=
mp,dens,4,24.8 Cr0
\7
tb,dp,4 q@ !p
tbdata,4,8e4 lsmzy_gV7
tbdata,2,40 2y;
|6`
tbdata,3,0 pl62mp!
!创建边坡模型线
/;LteBoY
k,1,101.95,69.87, mw}obblR
k,2,103.17,55.92 |rJ_
k,3,103.77,49.01 &We1i&w
k,4,104.57,39.96 h.0Y!'?
k,5,105,35 48,Aq*JFw
k,6,140,35 9D74/3b*
k,7,140,0 ~:7AHK2
k,8,105,0 L8-
k,9,0,0 il^SGH
k,10,0,35 ,Yhy7w
k,11,0,49.11 R<W#.mpo6
k,12,0,63.68 0[6llcuj
k,13,0,66.77 7O5`&Z'-
k,14,0,69.87 ( 8+ _~_
l,1,2 FI@2KM
l,2,3 g*8LdH6mq
l,3,4 'Rk~bAX
l,4,5 #sit8k`GR8
l,5,6 $t5
0<1
l,6,7 v8g3]MVj3
l,7,8 QrDrdA
l,8,9 Sd'Meebu
l,9,10 }D.?O,ue
l,10,11 A(AyLxB47*
l,11,12 1xK'T_[
l,12,13 $eQ_!7Gom$
l,13,14 g3{)AX[Uy
l,14,1 j}YZl@dYV
l,13,2 20 gPx;
l,12,3 .K^'Q|?
l,11,4 I%xrDiK97
l,10,5 QI>yi&t
l,5,8 L;.6j*E*
!创建边坡面模型 } {/4sll
al,14,1,15,13 AEkgm^t.{
al,15,2,16,12 K8*QS_*
al,16,3,17,11 &|<~J(L;
al,17,4,18,10 g-xbb&]
al,18,19,8,9 +T[3wL~
al,5,6,7,19 /:dVW"A|
vext,1,6,1,,,35,,, :K{`0U&l5
vsel,s,volu,,1 6Nx T W
vatt,1,1,1 bTC2Ya
vsel,s,volu,,2 QQ@, v@j5
vatt,2,2,1 s<0yQ-=.?N
vsel,s,volu,,3 %kD WUJZ
vatt,3,3,1 E*Vx^k$
vsel,s,volu,,4,6,1 :7Mo0,Bw,
vatt,4,4,1 {B|U8j[
allsel wbI(o4rXE
!划分网格 mm%w0dOb"
mshkey,1 ZTd_EY0 q
mshape,1,2D OB3AZH$
vmesh,all ,G1|]
~
!固定面11,15,19,23,26,27,30,31 Y">tfLIL_
da,11,all U0u @[9!
da,15,all y5?T`ts,#
da,19,all Faa:h#
da,23,all eG# (9
da,26,all Rz`<E97-
da,27,all P&| =
da,30,all 9<<$uf.B
da,31,all @FkNT~OZ
!施加重力 O89<IXk
allsel,all lsW.j#yE!
acel,,9.8,0 do*`-SDy
!进入求解器 u]E% R&
/solu ;~z>GJox
antype,static Ndi'b_Sh\
nusubst,100 3f's>+,#%
pred,on En6fmEn&;o
nropt,full H*|Bukgt/M
nlgeom,on Pb} &c
insrch,on ?>q=Nf^ Q.
outres,all,last MS^,h>KI
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 DS>s_3V
cnvtol,u,,0.05,2,1 !位移收敛准则设定 8q7KqYu
allsel -j3 -H&
solve ?-<t-3%hyV
save K1O/>dN_\O
finish 9b]*R.x:$&
/post1 !进入后处理 EN.yU!N.4
set,last <!X'- >i%q
plesol,s,y G+[>or}
!plesol,u,y,0,1 2 )RW*Qu;+
!plesol,s,xy,0,1 iBJ*6orz
!plesol,epto,xy,0,1 oZ6xHdPc4
!pldisp,1 !绘制边坡模型变形图 \Qk:\aLR
!plesol,u,x !绘制边坡模型水平方向位移云图 WbH/K]/1)h
!plesol,eppl,eqv !绘制边坡模型塑性应变云图 ;~~Oc