后处理中,发现自重应力场很混乱,明显有问题,见图。 /ig:9R
请问高人如何解决?万分感谢! \'CDRr"uw
命令如下: OnK~3j
finish D;_ MPN[
/CLEAR,start 8'f4 Od ?
/fILEname,slope analysis lhw ,J]0*
/units,si VxXzAeM
/pnum,kp,1 DBT&DS
/pnum,line,1 ^9ePfF)5
/number,0 -*m+(7G\
/PREP7 }b0; 0j
*afun,deg <_XWWT%
!定义单元类型 t$A%*JBKm
et,1,solid45 %"af748!+D
!定义材料属性 -$WYj"
!边坡围岩材料1的属性 l?Fb ='#
mp,ex,1,6e7 @)-$kk*
mp,nuxy,1,0.29 &d5ia+#
mp,dens,1,23 tWoh''@#
tb,dp,1 mGGsB5#w>
tbdata,1,3e4 kjtjw1\o
tbdata,2,30 Hv\-_>}K
tbdata,3,0 "sl1vzRN
!边坡围岩材料2的属性 ]@0NO;bK>F
mp,ex,2,5e7 :P@rkT3Q t
mp,nuxy,2,0.21 ]- 4QNc=
mp,dens,2,23.5 cg*)0U-_(
tb,dp,2 m/qbRk68s
tbdata,1,5e4 /Ne<V2AX
tbdata,2 61jI
tbdata,3,0 ")!,ZD
!边坡围岩材料3的属性 #*g5u{k'P
mp,ex,3,7e8 I<8sI%,s
mp,nuxy,3,0.34 |7}CQU
mp,dens,3,25.1 ZG du|
tb,dp,3 6'RrQc=q
tbdata,3,5e5 H03jDM8Q
tbdata,2,40 &ZX{R#[L
tbdata,3,0 aN $}?
!边坡围岩材料4的属性 +C(/Lyo}
mp,ex,4,1.5e9 '8T=~R6
mp,nuxy,4,0.3 ty*@7g0k
mp,dens,4,24.8 pTyi!:g3W
tb,dp,4 L0tAgW!@
tbdata,4,8e4 A^2Uzmzl?
tbdata,2,40 mK [0L
tbdata,3,0 -atGlu2
!创建边坡模型线 ^+m+zd_
k,1,101.95,69.87, !Wy[).ZAf
k,2,103.17,55.92 zdEPDdB
k,3,103.77,49.01 p$x{yz3
k,4,104.57,39.96 E)9yH\$6
k,5,105,35 Sf,R^9#|
k,6,140,35 kr9gK~
k,7,140,0 !\,kZ|#>
k,8,105,0 ;XDz)`c
k,9,0,0 +5&wOgx
k,10,0,35 k!KDWb
k,11,0,49.11 -~QHqU.
k,12,0,63.68 Gsc\/4Wx
k,13,0,66.77 0sh/|`\
k,14,0,69.87 zWb4([P;
l,1,2 NSFs\a@1
l,2,3 ~~6^Sh60g
l,3,4 .^m>AKC0cX
l,4,5 q=DN
{a:
l,5,6 f~_th @K
l,6,7 /2HN>{F^Y
l,7,8 Cc, `}SP
l,8,9 7zv1wb
l,9,10 ]+m/;&0
l,10,11 jOyvDY9\
l,11,12 PGARXw+
l,12,13 ^_%kE%I
l,13,14 F1Hh7
F
l,14,1 'D%w|Pe?Q
l,13,2 M!tXN&V]
l,12,3 A?oXqb
l,11,4 @mBX~ ?=Z3
l,10,5 ?"KC-u|
l,5,8 a+9*@z2
!创建边坡面模型 AT\qiznvP
al,14,1,15,13 F|HJH"2*&q
al,15,2,16,12 5 XA=G
al,16,3,17,11 ]l(wg]
al,17,4,18,10 5&e<#"
al,18,19,8,9 &k1T08C*
al,5,6,7,19 >"@?ir
vext,1,6,1,,,35,,, Cua%1]"4w
vsel,s,volu,,1 e[Jem5C
vatt,1,1,1 E3*\
^Q_
vsel,s,volu,,2 {"
4e+y
vatt,2,2,1 ad_`x
vsel,s,volu,,3 2]c{P\
vatt,3,3,1 ee/&/Gt
vsel,s,volu,,4,6,1 #%FN>v3e
vatt,4,4,1 n?:=
allsel zA;@@)hwR
!划分网格 $DP&a1'g
mshkey,1 hHPs&EA.p
mshape,1,2D ?9X&tK)E-
vmesh,all P@]8pIB0d^
!固定面11,15,19,23,26,27,30,31 Hku!bJ
da,11,all 6y5A"-
da,15,all 5&f{1M6l>
da,19,all P/ oXDI8
da,23,all rO:u6."_
da,26,all : ~ A%#
da,27,all wO>L#"X^v
da,30,all !P$'#5mr
da,31,all \i[BP
!施加重力 Z^9/v
allsel,all er.CDKD%L
acel,,9.8,0 \)48904^
!进入求解器 ^o !O)D-q
/solu A~dQ\M
antype,static L}yyaM)
nusubst,100 /n4pXT
pred,on o|j*t7
nropt,full /S\cU`ZVe
nlgeom,on JNFIT;L
insrch,on fN
"tA
outres,all,last P &)1Rka
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 (LtkA|:
cnvtol,u,,0.05,2,1 !位移收敛准则设定 X{g%kf,D=
allsel gLSA!#[h
solve ((rv]f{
save A`u$A9[
finish '?Jxt:<
/post1 !进入后处理 f2,1<^{
set,last P=5NKg
plesol,s,y V
>,Z-&.%
!plesol,u,y,0,1 <q,+ON\'
!plesol,s,xy,0,1 Cj*-[EL<
!plesol,epto,xy,0,1 IAOcKQ3
!pldisp,1 !绘制边坡模型变形图
pAu72O?
!plesol,u,x !绘制边坡模型水平方向位移云图 Oc&),ru2l
!plesol,eppl,eqv !绘制边坡模型塑性应变云图 v[lnw} =m9