[求助]为什么自重应力场这么混乱？ [复制链接]

后处理中，发现自重应力场很混乱，明显有问题，见图。 -/f$s1  
请问高人如何解决？万分感谢！ c(Q@5@1y:  
命令如下： ZW4
f "  
finish N>mW64_H)  
/CLEAR,start [WwoGg*)mn  
/fILEname,slope analysis z^bS+0S5x!  
/units,si \&"C  
/pnum,kp,1 EX,>V,.UV  
/pnum,line,1 <jh4P!\&j  
/number,0 UU=]lWib  
/PREP7 '$ei3  
*afun,deg gOF^?M11x  
!定义单元类型 rN0<y4)!  
et,1,solid45 05d0p|},
 
!定义材料属性 QiK>]xJ'  
!边坡围岩材料1的属性 <PLAAh8  
mp,ex,1,6e7 IJ^~,+  
mp,nuxy,1,0.29 `v*HH}aDO  
mp,dens,1,23 X[BP0:`t  
tb,dp,1 `Ba]i)!  
tbdata,1,3e4 R1];P*>%gZ  
tbdata,2,30 BT7{]2?&V  
tbdata,3,0 "\4W])30  
!边坡围岩材料2的属性 ,!orD1,'  
mp,ex,2,5e7 K%q5:9m  
mp,nuxy,2,0.21 80lei  
mp,dens,2,23.5 Z?
)g'n  
tb,dp,2 :i.@d?  
tbdata,1,5e4  -WC0W  
tbdata,2 !XPjRdq  
tbdata,3,0 +,e#uuj$p  
!边坡围岩材料3的属性 +B*]RL[th  
mp,ex,3,7e8 ;aQ``B  
mp,nuxy,3,0.34 _ *f>UW*,  
mp,dens,3,25.1 2`o @L  
tb,dp,3 IO)B3,g  
tbdata,3,5e5 :ld~9  
tbdata,2,40 VLuHuih  
tbdata,3,0 5(>=};r+  
!边坡围岩材料4的属性 vJWBr:`L  
mp,ex,4,1.5e9 gAAC>{Wh  
mp,nuxy,4,0.3 D<WnPLA$g  
mp,dens,4,24.8 o/pw=R/):  
tb,dp,4 z,,"yVk`,  
tbdata,4,8e4 OBi(]l}^O  
tbdata,2,40 &8$v~  
tbdata,3,0 .}==p
&(  
!创建边坡模型线 Q'OtXs 80  
k,1,101.95,69.87, sfF~
k-  
k,2,103.17,55.92 ibkB>n{(  
k,3,103.77,49.01 ;Y7'U rn  
k,4,104.57,39.96 #Y7jNrxE
  
k,5,105,35 Oeh A3$|#  
k,6,140,35 $|K: 9  
k,7,140,0 Qs1p  
k,8,105,0 /VYT](  
k,9,0,0 Oe5rRQ$O  
k,10,0,35 H/o_?qK  
k,11,0,49.11 u*C*O4f>OC  
k,12,0,63.68 ~$J;yo
~  
k,13,0,66.77 JduO^Fit  
k,14,0,69.87 (*S<2HN5  
l,1,2 ng}C$d . I  
l,2,3 ~?\U];l  
l,3,4 q?!HzZ  
l,4,5 63l3WvoK  
l,5,6 ~F"S]  
l,6,7 2!#
g\"  
l,7,8 KqI:g*H'x7  
l,8,9 SBxpJsW>  
l,9,10 /4x\}qvU  
l,10,11 #/oH #/?  
l,11,12 ^4`aONydl  
l,12,13 0qS/>u*  
l,13,14 y3~`qq  
l,14,1 f@i#Znkf*?  
l,13,2 j_<qnBeQ  
l,12,3 Fm`*j/rq
 
l,11,4 \F|)w|v  
l,10,5 '+9<[]  
l,5,8 78l);/E{v  
!创建边坡面模型 yCQvo(V[F  
al,14,1,15,13 ]0UYxv%]  
al,15,2,16,12 -06G.;W\^  
al,16,3,17,11 I=!kPuw  
al,17,4,18,10 t3)6R(JC  
al,18,19,8,9 "xlR>M6e  
al,5,6,7,19
|[`YGA4  
vext,1,6,1,,,35,,, w)7y{ya$  
vsel,s,volu,,1 oJa}NH   
vatt,1,1,1 [* <x)  
vsel,s,volu,,2 d6n_Hpxw^  
vatt,2,2,1 ,U""m7   
vsel,s,volu,,3 Lm[,^k  
vatt,3,3,1 _r:Fmn_%-  
vsel,s,volu,,4,6,1 _=HaE&  
vatt,4,4,1 B_[efM<R$  
allsel $gr>Y2i  
!划分网格 '8 .JnCg  
mshkey,1 cKTjQJ#  
mshape,1,2D uWJJ\  
vmesh,all 2;gvo*k  
!固定面11，15，19，23，26，27，30，31 v]EMJm6d|  
da,11,all "^&H9.z,v  
da,15,all 8H7#[?F  
da,19,all k'8tcXs  
da,23,all )0/9 L  
da,26,all <tvLKx  
da,27,all KKC%!Xy  
da,30,all LK}
g<!o(  
da,31,all %`i*SF(gV  
!施加重力 J%?5d:iN+  
allsel,all F0!Z1S0g  
acel,,9.8,0 &g;!n&d zP  
!进入求解器 ,7d/KJ^7  
/solu -riX=K>$  
antype,static ,zgz7  
nusubst,100 ,sitOy}ks  
pred,on SCfp5W7~  
nropt,full '~Y@HRVL@|  
nlgeom,on LG8h@HY&L  
insrch,on N**)8(  
outres,all,last DQGrXMpV0  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 V/#Ra  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 5Npxs&Ea  
allsel CHGa_  
solve 7<su8*?  
save XP)^81i|  
finish PA'&]piPl:  
/post1  !进入后处理 0si1:+t-[+  
set,last -J3~j kf  
plesol,s,y 5n?P}kca)  
!plesol,u,y,0,1 4x6n,:;  
!plesol,s,xy,0,1 *Z#OfB4}  
!plesol,epto,xy,0,1 t*hy"e{*a  
!pldisp,1                                      !绘制边坡模型变形图 O14QlIk  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 [gT}<W  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图 .e7tq\k  

自己顶个，等待高人！

是不是密度的单位拟弄错了 Tc(=J7*r&  
怎么那么小

你这都收敛不了呀

finish                                     fz^j3'!\  
/CLEAR,start                           ;(?tlFc  
/fILEname,try                   t}eyfflZ  
/units,si                                 M/6Z,oOU  
/pnum,kp,1                                 gl+d0<Rzw  
/pnum,line,1                             U1r]e%df)  
/number,0                               UiG/Rn  
/PREP7                                     D#LV&4e>.E  
*afun,deg                                 +7b8ye  
!定义单元类型                             \<pr28  
et,1,solid45                               Ip'tB4Mq  
!定义材料属性                               CIIjZ)T  
!边坡围岩材料1的属性                       qf(mJlU  
mp,ex,1,6e7                               i#V(oSx  
mp,nuxy,1,0.29                             j"nOxs  
mp,dens,1,23                               zzZEX  
tb,dp,1                                   tfU*U>j  
tbdata,1,3e4,30,0                                                             r+bGZ  
!边坡围岩材料2的属性                     \Sd8PGl*'  
mp,ex,2,5e7                               WrGA7&!+  
mp,nuxy,2,0.21                             8NHm#Z3Ol  
mp,dens,2,23.5                             LJ[zF~4#  
tb,dp,2                                   !vnC-&G  
tbdata,1,5e4,0,0                                                               ^0.8-RT  
!边坡围岩材料3的属性                     LH4A!a]  
mp,ex,3,7e8                               |@AXW   
mp,nuxy,3,0.34                               lfj
5?y  
mp,dens,3,25.1                             To*+Z3Wd  
tb,dp,3                           lndz  
tbdata,1,5e5,40,0                                                               6&x\!+]F8  
!边坡围岩材料4的属性                       -!XG>Z  
mp,ex,4,1.5e9                             P{9:XSa%  
mp,nuxy,4,0.3                             ~W21%T+  
mp,dens,4,24.8                             [1vm~w'  
tb,dp,4                                   5e+j51  
tbdata,1,8e4,40,0                                                           |H t5a.  
!创建边坡模型线 5                           n>}Y@{<]/  
k,1,101.95,69.87                           :dt[ #  
k,2,103.17,55.92             Q~]oN  
k,3,103.77,49.01                         +fk*c[FG  
k,4,104.57,39.96                         vb=]00c  
k,5,105,35                               eS!]..%y  
k,6,140,35                               D<J'\
mo  
k,7,140,0                               <gwRE{6U  
k,8,105,0                                   E /H%q|q  
k,9,0,0                                   !
}!KT(%%  
k,10,0,35                                   6X@z(EEL  
k,11,0,49.11                               4

"2%mx:  
k,12,0,63.68                               y5sH7`2+5  
k,13,0,66.77                               .u)KP*_  
k,14,0,69.87                               t,9+G<)>H  
l,1,2                                   hTDGgSG^  
l,2,3                                   \< .BN;t{  
l,3,4                                   c6[m'cy  
l,4,5                                   mEJ7e#  
l,5,6                                   ;JkIZ8!  
l,6,7                                   DKS1Sm6d0  
l,7,8                                   "0!h-bQN  
l,8,9                                   ATkd#k%S  
l,9,10                                     iW6O9~  
l,10,11                                 t+KW=eW  
l,11,12                                 tKt}]KHV  
l,12,13                                 H7'42J@  
l,13,14                                 ZhY{,sy?QO  
l,14,1                                   X]y8-}Qf  
l,13,2                                   []A9j?_w  
l,12,3                                   @m5J%8>k  
l,11,4                                   rT`sY  
l,10,5                       fDsT@W,K  
l,5,8                                     aA4RC0'  
!创建边坡面模型                         - jZA
vb  
al,14,1,15,13                     B[ZQn]y  
al,15,2,16,12                     %:~Ah6R1  
al,16,3,17,11                     edMCj  
al,17,4,18,10                     pl>b 6 |  
al,18,19,8,9                   / <(|4e  
al,5,6,7,19                               w
<B S  
vext,1,6,1,,,35,,,                           >3/mV<g f  
vsel,s,volu,,1                             f$>_>E  
vatt,1,1,1                                 0hoMf=bb$  
vsel,s,volu,,2                           W_wC"?A%  
vatt,2,2,1                                 1#x5 o2n  
vsel,s,volu,,3                             C[,h!  
vatt,3,3,1                                 ,&,%B|gT]  
vsel,s,volu,,4,6,1                         0cG[<\qT  
vatt,4,4,1                           G|jHic!  
allsel,all                                     fc~fjtqwvz  
!划分网格                                   D]E=0+  
mshkey,1                                 ;:=j{,&dl[  
mshape,1,3D                               cgG*7E  
vmesh,all                                     gqv+|:#  
da,11,all                                 sN1H{W  
da,15,all                                 D N!V".m`J  
da,19,all                                 B5 /8LEWw  
da,23,all                                 1#N`elm  
da,26,all                                 Lz1KDXr`)+  
da,27,all                                 rCa]T@=  
da,30,all                                 l0G{{R0Y  
da,31,all                                 h1q3}-  
!施加重力                                 O h e^{:  
allsel,all                               8(n>99VVK  
acel,,9.8,0                               |vd|;" `  
!进入求解器                               K-]) RIM  
/solu                                     M;XU"8  
antype,static 0i}.l\  
nlgeom,on                             Dl.<(/  
nsubst,500,500,200                                                                                                                                                         P\dfxR;8%  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定     |\Gkhi>;  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 OD}Uc+;K  
allsel  

材料3和4定义错误，应该是TBDATA,1,.....不是TBDATA,3 ;LgMi5dN  
定义子步出错，命令流就是错的

给你改了半天也没回来看 \=RV?mI3?  
更没有感谢 elM<S3  
FENT