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

后处理中，发现自重应力场很混乱，明显有问题，见图。 (< +A w7  
请问高人如何解决？万分感谢！ (W}F\P  
命令如下： [
ZWAXl $  
finish B+W7zv  
/CLEAR,start :ld~9  
/fILEname,slope analysis nI:M!j5s`  
/units,si t2[/eM.G  
/pnum,kp,1 gAAC>{Wh  
/pnum,line,1 D<WnPLA$g  
/number,0 o/pw=R/):  
/PREP7 ]]j^  
*afun,deg s6H.Q$3L  
!定义单元类型 Fp\;j\pfw  
et,1,solid45 \<{a=@_k9  
!定义材料属性 __=53]jGE  
!边坡围岩材料1的属性 1!z{{H;W  
mp,ex,1,6e7 ]az(w&vqg2  
mp,nuxy,1,0.29 e=i X]%^  
mp,dens,1,23 ~[;r) g\  
tb,dp,1 O= S[n  
tbdata,1,3e4 ,Lig6Z`  
tbdata,2,30 JK$3qUDnI  
tbdata,3,0 "&6vFmr  
!边坡围岩材料2的属性 $d<NN2  
mp,ex,2,5e7 K43%9=sM  
mp,nuxy,2,0.21 ~$J;yo
~  
mp,dens,2,23.5 JduO^Fit  
tb,dp,2 (*S<2HN5  
tbdata,1,5e4 w;'XqpP$*|  
tbdata,2 <6UXk[y  
tbdata,3,0 m23"xnRB  
!边坡围岩材料3的属性 P< OH{l  
mp,ex,3,7e8 E'x
"EN  
mp,nuxy,3,0.34 g+#<;Gbpe  
mp,dens,3,25.1 H^d?(
Svh  
tb,dp,3 ZJDV'mC}  
tbdata,3,5e5  2%@tnk|@  
tbdata,2,40 P~a@{n*8  
tbdata,3,0 2uj .*  
!边坡围岩材料4的属性 Y 9$jJ1V  
mp,ex,4,1.5e9 UarLxPQ  
mp,nuxy,4,0.3 QAkK5,`vV.  
mp,dens,4,24.8 XvI~"}  
tb,dp,4 aC1z.?!U  
tbdata,4,8e4 ]IV{;{E)  
tbdata,2,40 $hm[x$$  
tbdata,3,0 o`YBz~2  
!创建边坡模型线 m.D8@[y  
k,1,101.95,69.87, $Cz2b/O  
k,2,103.17,55.92 mG[jR*JW  
k,3,103.77,49.01 iT'doF  
k,4,104.57,39.96 5}%R  
k,5,105,35 ?lC>E[  
k,6,140,35 0YiTv;mq;  
k,7,140,0 a
2P)@R  
k,8,105,0 "
rBB&l  
k,9,0,0 )Rn\6ka  
k,10,0,35 ,`f]mv l  
k,11,0,49.11 )+N{D=YM  
k,12,0,63.68 A4zI1QF  
k,13,0,66.77 'g%:/lwA  
k,14,0,69.87 z"f@iJX?2  
l,1,2 NK0'\~7&  
l,2,3 J _rrc;F  
l,3,4 Sr \y1nt  
l,4,5 4X^$"lM  
l,5,6 Y_CYx  
l,6,7 (\ab%M
  
l,7,8 >6S7#)0T  
l,8,9 8UU
L=  
l,9,10 (.UU40:t  
l,10,11 F!z ^0+H(  
l,11,12 6Z|h>H5a  
l,12,13 J%?5d:iN+  
l,13,14 F0!Z1S0g  
l,14,1 &g;!n&d zP  
l,13,2 ,7d/KJ^7  
l,12,3 -riX=K>$  
l,11,4 ,zgz7  
l,10,5 yksnsHs}d  
l,5,8 A:m+v{*`4  
!创建边坡面模型 '&cH,yc;b  
al,14,1,15,13 S%- k
N;  
al,15,2,16,12 rt^~ I\V  
al,16,3,17,11 /ZyMD(_J  
al,17,4,18,10 Jg$<2CR&  
al,18,19,8,9 %,WH*")  
al,5,6,7,19  26p[x'W  
vext,1,6,1,,,35,,, D('.17  
vsel,s,volu,,1 P2vG)u  
vatt,1,1,1 &|x7T<,)  
vsel,s,volu,,2 >ZJ]yhbhK  
vatt,2,2,1 Hs)Cf)8u  
vsel,s,volu,,3  x+cL(R  
vatt,3,3,1 lJ;Wi  
vsel,s,volu,,4,6,1 #@oB2%&X?  
vatt,4,4,1 [W3X$r~-  
allsel pg%(6dqK4  
!划分网格 Y#}qXXZ>]  
mshkey,1 i D9 */  
mshape,1,2D GQYR`;>  
vmesh,all 'Q=;I  
!固定面11，15，19，23，26，27，30，31 h/n(  
da,11,all =803rNe  
da,15,all a <C?- g|  
da,19,all +j{Cfv$do  
da,23,all 8?iI;(  
da,26,all P4@`C{F5m  
da,27,all {dZ]+2Z~+  
da,30,all 1 tOslP@  
da,31,all nF'YG+;|@  
!施加重力 pq`uB  
allsel,all x|m9?[ !_  
acel,,9.8,0 yx Om=V  
!进入求解器 4Kch=jt4#  
/solu :.5l  
antype,static 36Lkcda[  
nusubst,100 X&WP.n)  
pred,on %M`48TW)  
nropt,full u@dvFzc  
nlgeom,on 0Fb];:a  
insrch,on K T0t4XPM  
outres,all,last X.>=&~[  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 |R
L#BKC`  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 `,6|6.8#  
allsel 'Ou C[$Z  
solve &jqylX  
save y@!M<#SEzG  
finish 0Agse)  
/post1  !进入后处理 ]L~NYe9  
set,last RiQ]AsTtl  
plesol,s,y OK]
QDb  
!plesol,u,y,0,1 Aax;0qGbH  
!plesol,s,xy,0,1 d/R!x{$-f  
!plesol,epto,xy,0,1 Q3WI@4  
!pldisp,1                                      !绘制边坡模型变形图 opa}z-7>^  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 N" L&Z4Z  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图 y.lWyH9  

自己顶个，等待高人！

是不是密度的单位拟弄错了 ZHK>0>;  
怎么那么小

你这都收敛不了呀

finish                                     8
QTry%  
/CLEAR,start                           R0=f`;  
/fILEname,try                   i%9vZ  
/units,si                                 U"7o;q  
/pnum,kp,1                                  Gk~aTO  
/pnum,line,1                             2V@5:tf  
/number,0                               4fKvB@O@.  
/PREP7                                     y[XD=j  
*afun,deg                                 >B{qPrmI  
!定义单元类型                             ;JkIZ8!  
et,1,solid45                               DKS1Sm6d0  
!定义材料属性                               @4;'>yr(  
!边坡围岩材料1的属性                       I#MPJ@*WT  
mp,ex,1,6e7                               \>5sW8P]H`  
mp,nuxy,1,0.29                             sLA.bp.O  
mp,dens,1,23                               i?/Q7D<P  
tb,dp,1                                   >\1twd{u]  
tbdata,1,3e4,30,0                                                             6!]@S|vDX  
!边坡围岩材料2的属性                     BHY-fb@R]H  
mp,ex,2,5e7                               aE'nW_f  
mp,nuxy,2,0.21                             Yd3lL:M  
mp,dens,2,23.5                             ]S%_&ZMCM  
tb,dp,2                                   :=B.)]F.)  
tbdata,1,5e4,0,0                                                               t5E$u(&+'B  
!边坡围岩材料3的属性                     SPV+ O{  
mp,ex,3,7e8                               PaSwfjOnqr  
mp,nuxy,3,0.34                               <Z-Pc?F&(k  
mp,dens,3,25.1                             OH>.N"IG  
tb,dp,3                           7SHllZ  
tbdata,1,5e5,40,0                                                               'aEK{#en  
!边坡围岩材料4的属性                       *B3f ry  
mp,ex,4,1.5e9                             q(qm3OxYo  
mp,nuxy,4,0.3                             {LiJ=Ebt  
mp,dens,4,24.8                             \NNA"  
tb,dp,4                                   %O9Wm_%  
tbdata,1,8e4,40,0                                                           @S3L%lOH  
!创建边坡模型线 5                           1R}9k)JQ  
k,1,101.95,69.87                           +~V_^-JG&  
k,2,103.17,55.92             >l 0aME@-0  
k,3,103.77,49.01                         D]E=0+  
k,4,104.57,39.96                         ;:=j{,&dl[  
k,5,105,35                               c$1
u  
k,6,140,35                               .h <=C&Yg  
k,7,140,0                               sN1H{W  
k,8,105,0                                   D N!V".m`J  
k,9,0,0                                   B5 /8LEWw  
k,10,0,35                                   1#N`elm  
k,11,0,49.11                               gba1R  
k,12,0,63.68                               r#WqXh_uk  
k,13,0,66.77                               <9&GOaJ  
k,14,0,69.87                               C
t `)R  
l,1,2                                   f1:>H.m`  
l,2,3                                   "S#$:92  
l,3,4                                   Ei:m@}g  
l,4,5                                   B-LV/WJ_  
l,5,6                                   L&+k`b
  
l,6,7                                   N_ ODr]L  
l,7,8                                   n}Z%-w$K#  
l,8,9                                   ~EmK;[Z  
l,9,10                                     .*N,x(V  
l,10,11                                 B4un6-<i  
l,11,12                                 ED8
{  
l,12,13                                 F7<M{h5s  
l,13,14                                 7a_8007$l  
l,14,1                                   pImq<Z  
l,13,2                                   tfHr'Qy BC  
l,12,3                                   ]A[}:E 5}  
l,11,4                                   7J3A]>q
U  
l,10,5                       :Jyr^0`J  
l,5,8                                     >0SG]er@  
!创建边坡面模型                         z>+CMH5L)  
al,14,1,15,13                     @wD#+Oz  
al,15,2,16,12                     ;LgMi5dN  
al,16,3,17,11                     8cfsl lI  
al,17,4,18,10                     {KSLB8gtL  
al,18,19,8,9                   Pjq9BK9p  
al,5,6,7,19                               jA^yUd-  
vext,1,6,1,,,35,,,                           L__{U_p  
vsel,s,volu,,1                             v3[ZPc;;  
vatt,1,1,1                                 bo/U5p
 
vsel,s,volu,,2                           G-D}J2r=F  
vatt,2,2,1                                 d3a!s  
vsel,s,volu,,3                             R[j'<gd.  
vatt,3,3,1                                 A PSkW9H  
vsel,s,volu,,4,6,1                         GF@`~im  
vatt,4,4,1                           r!w4Br0  
allsel,all                                     Zah<e6L  
!划分网格                                   \gT({XU?  
mshkey,1                                 x%)oL:ue  
mshape,1,3D                               &iD&C>;pf  
vmesh,all                                     -^_^ByJe  
da,11,all                                 xX%{i0E  
da,15,all                                 
~rrl"a>  
da,19,all                                 .HTRvE`X  
da,23,all                                 k3LHLJZ#  
da,26,all                                 %a8e_  
da,27,all                                 KKk<wya&O  
da,30,all                                 F+5 5p8  
da,31,all                                 }x6)}sz7  
!施加重力                                 bfb9A+]3'  
allsel,all                               DIYR8l}x  
acel,,9.8,0                               A~{vja0?  
!进入求解器                               9Pvv6WyKy  
/solu                                     BSkmFd(*  
antype,static dM}c-=w`  
nlgeom,on                             B{`K?e0  
nsubst,500,500,200                                                                                                                                                         Q-iBK*-w  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定     j7Zv"Vq@  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 p1GP@m,^n0  
allsel  

材料3和4定义错误，应该是TBDATA,1,.....不是TBDATA,3 F1}d@^K 7d  
定义子步出错，命令流就是错的

给你改了半天也没回来看 .2|(!a9W  
更没有感谢 FBrh!vQ<  
FENT