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

后处理中，发现自重应力场很混乱，明显有问题，见图。 Cq(dj^/~m  
请问高人如何解决？万分感谢！ py6|uGN  
命令如下： #qtAFIm'  
finish q~48lxDU  
/CLEAR,start Lo5itW  
/fILEname,slope analysis gu7mGHn-  
/units,si 5IE2&V  
/pnum,kp,1 yF8 av=<{  
/pnum,line,1 1>r7s*  
/number,0 ^OnU;8IC
 
/PREP7 ]v{f!r=}  
*afun,deg %O#)=M~  
!定义单元类型 4p\<b8(9>  
et,1,solid45 0GK<l  
!定义材料属性 %QCh#v=ks  
!边坡围岩材料1的属性 Y!N*J  
mp,ex,1,6e7 3Un q 9  
mp,nuxy,1,0.29 01uMbtM  
mp,dens,1,23 33z)F  
tb,dp,1 CkKr@.dV  
tbdata,1,3e4 L}T:Y).  
tbdata,2,30 7O^ySy"l  
tbdata,3,0 WwG78b-OA  
!边坡围岩材料2的属性 71_{FL8  
mp,ex,2,5e7 rXPq'k'h#-  
mp,nuxy,2,0.21 o{f|==<t3#  
mp,dens,2,23.5 ;}"_hLX  
tb,dp,2 ^V#9{)B  
tbdata,1,5e4 <csz4tL}P  
tbdata,2 >RZ]t[)y  
tbdata,3,0 $6
e&sDJ  
!边坡围岩材料3的属性 eKPxSN Z  
mp,ex,3,7e8 `$vT
GkGpY  
mp,nuxy,3,0.34 fl#gWAM  
mp,dens,3,25.1 BQu_)@  
tb,dp,3 Cp` [0v~0  
tbdata,3,5e5 A[`c2v-hF  
tbdata,2,40 {/#^v?,  
tbdata,3,0 >qk[/\^O  
!边坡围岩材料4的属性 Y&!M#7/'J3  
mp,ex,4,1.5e9 xOnbYU  
mp,nuxy,4,0.3 ?f<JwF<  
mp,dens,4,24.8 1
BO$xq  
tb,dp,4 M0zJGIT~b  
tbdata,4,8e4 KJ M:-z@  
tbdata,2,40 A{3Aw|;  
tbdata,3,0 ()a(PvEO
 
!创建边坡模型线 5N<f\W,  
k,1,101.95,69.87, KE!aa&g  
k,2,103.17,55.92 3omFd#EP  
k,3,103.77,49.01 7N9~nEU  
k,4,104.57,39.96 bL soKe  
k,5,105,35 $fES06%  
k,6,140,35 qDHiyg^u  
k,7,140,0 {\3k(NdEX  
k,8,105,0 &D0suK#  
k,9,0,0 x(pq!+~K  
k,10,0,35 K(@QKRZ7[  
k,11,0,49.11 })RT2zw}  
k,12,0,63.68 Vz[E)(QX-`  
k,13,0,66.77 |W4 \  
k,14,0,69.87 golr,+LSo  
l,1,2 5mxHOtvtWM  
l,2,3 RP 'VEJ  
l,3,4 [-C-+jC  
l,4,5 k]?M^jrm  
l,5,6 f'P}]_3(  
l,6,7 M0VC-\W7f  
l,7,8 U'(}emh}  
l,8,9 E`s9SE  
l,9,10 9 f$S4O5  
l,10,11 ~_<I}!j/B  
l,11,12 &Uq++f6  
l,12,13 |:yQOq|  
l,13,14 <{!^  
l,14,1 ]AkHNgW  
l,13,2 2ED^uc: 0S  
l,12,3 Dlz||==  
l,11,4 uis;S)+  
l,10,5 mw0#Dhyy1=  
l,5,8 7LMad%  
!创建边坡面模型 YF/@]6j  
al,14,1,15,13 DwM4/m  
al,15,2,16,12 8&snLOU -Q  
al,16,3,17,11 \|^fG9M~  
al,17,4,18,10 N'hj  
al,18,19,8,9 bU'{U0lM  
al,5,6,7,19 w}n:_e  
vext,1,6,1,,,35,,,
CY2DxP%  
vsel,s,volu,,1 eBiP\  
vatt,1,1,1 KY'"Mg^!  
vsel,s,volu,,2 xHEkmL`)4  
vatt,2,2,1 $[9,1.?C  
vsel,s,volu,,3 obWBX'  
vatt,3,3,1 \%\b*OO  
vsel,s,volu,,4,6,1 $__e7  
vatt,4,4,1 <qZ"W6&&  
allsel ~@a7RiE@  
!划分网格 lOp/kGmn+  
mshkey,1 K^z-G=|N  
mshape,1,2D w6tb vhcmU  
vmesh,all fq-$u;~h  
!固定面11，15，19，23，26，27，30，31 :()(P9?  
da,11,all {MX_t/o=f  
da,15,all `< 8Fc`;[  
da,19,all N,v4SIC@  
da,23,all T8Sgu6:*R  
da,26,all !R,9Pg*Ey  
da,27,all x" 2
1 Jh  
da,30,all }G n2%  
da,31,all Ht^MY  
!施加重力 W}B4^l
 
allsel,all ]:TX> X!  
acel,,9.8,0 ]P/eg$u'I  
!进入求解器 d5<@WI:wz  
/solu INwc@XB  
antype,static `K[:<p}  
nusubst,100 {m5tgVi&  
pred,on  (2vR8  
nropt,full .N~qpynY  
nlgeom,on Yh/-6wg  
insrch,on 4f,D3e%T|  
outres,all,last oH~ZqX.3  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 y$@d%U*rW^  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 YLk; ^?  
allsel rN#ydw:9  
solve 2,|*KN*e`W  
save kfb+OE:7  
finish (?-5p;  
/post1  !进入后处理 AfvTStwr  
set,last rC]jz$sle  
plesol,s,y aM'0O![d  
!plesol,u,y,0,1 ^EU&6M2  
!plesol,s,xy,0,1 6suB!XF;  
!plesol,epto,xy,0,1 =DT
n9}u  
!pldisp,1                                      !绘制边坡模型变形图 ^5vFF@to  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 }8WpX2U  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图
F=   

自己顶个，等待高人！

是不是密度的单位拟弄错了 /6F 1=O(c>  
怎么那么小

你这都收敛不了呀

finish                                     4)2*|w  
/CLEAR,start                           qa^x4xZM  
/fILEname,try                   ygpC1nN  
/units,si                                 a}E8ADyC  
/pnum,kp,1                                 TOMvJ>bF  
/pnum,line,1                             3}g?d/^E3  
/number,0                               1:YDN.*  
/PREP7                                     AG}' W  
*afun,deg                                 .NjdkHYR  
!定义单元类型                             1b3Lan_2  
et,1,solid45                               2l5@gDk5  
!定义材料属性                               w0/W=!_  
!边坡围岩材料1的属性                       "G6d'xkP  
mp,ex,1,6e7                               w?M*n<) O  
mp,nuxy,1,0.29                             08k1 w,6W  
mp,dens,1,23                               O-T/H-J`  
tb,dp,1                                   qx0o,oZN!  
tbdata,1,3e4,30,0                                                             -<z'f){gb  
!边坡围岩材料2的属性                     Rz=]K
eZu  
mp,ex,2,5e7                               7C2/^x P  
mp,nuxy,2,0.21                             5@[%P=  
mp,dens,2,23.5                             MiRH i<g0  
tb,dp,2                                   73}k[e7e  
tbdata,1,5e4,0,0                                                               T>f-b3dk  
!边坡围岩材料3的属性                     [9(B;;R@  
mp,ex,3,7e8                               E=s`$ A  
mp,nuxy,3,0.34                               /u:Sn=SPd  
mp,dens,3,25.1                             ;o^m"I\y  
tb,dp,3                           L6 _Sc-sU  
tbdata,1,5e5,40,0                                                               6w1:3~a  
!边坡围岩材料4的属性                       SmR*b2U  
mp,ex,4,1.5e9                             vMRKs#&8  
mp,nuxy,4,0.3                             4zf#zJw  
mp,dens,4,24.8                             gP1$#KgU  
tb,dp,4                                   mz\m^g3  
tbdata,1,8e4,40,0                                                           `SW`d<+L  
!创建边坡模型线 5                           *"zE,Bp"  
k,1,101.95,69.87                           Wnf`Rf)1z  
k,2,103.17,55.92             *DNH_8m  
k,3,103.77,49.01                         ;4qalxzu  
k,4,104.57,39.96                         M$&>"%Oi  
k,5,105,35                               cQ+,F2  
k,6,140,35                               qA30z%#z_  
k,7,140,0                               GtGToI  
k,8,105,0                                   n<hwstk  
k,9,0,0                                   m9e$ZZG$  
k,10,0,35                                   wk
m SIN:  
k,11,0,49.11                               b/5~
VY*T  
k,12,0,63.68                               0e7v ?UT  
k,13,0,66.77                               "R9kF-  
k,14,0,69.87                               ~{NDtB)  
l,1,2                                   i-ogeR?  
l,2,3                                   &H+<uYV  
l,3,4                                   Q o=  
l,4,5                                   T2(+
HI2  
l,5,6                                   EvECA,!i  
l,6,7                                   R>0ta
Q  
l,7,8                                   O/-OW: 03  
l,8,9                                   3
QR-8  
l,9,10                                     dO8Z {
wfs  
l,10,11                                 iTK1I0  
l,11,12                                 HE,wEKp  
l,12,13                                 }Ql;%7  
l,13,14                                 A23Z)`  
l,14,1                                   pEaH^(I*  
l,13,2                                   JGdBpj:  
l,12,3                                   Sv>aZ  
l,11,4                                   dF?pEet?2  
l,10,5                       U)l>#gf8  
l,5,8                                     |ITg-t  
!创建边坡面模型                        
(f.A5~e  
al,14,1,15,13                     4%6Q+LS']Q  
al,15,2,16,12                     bFJn-g n  
al,16,3,17,11                     q;AQ6k(  
al,17,4,18,10                     v0xi(Wu  
al,18,19,8,9                  
g,W#3b6>j  
al,5,6,7,19                               "1nd~ BBOw  
vext,1,6,1,,,35,,,                           2'}2r ~6  
vsel,s,volu,,1                             ypifXO;m7  
vatt,1,1,1                                 eo,]b1C2n  
vsel,s,volu,,2                           3&*%>)  
vatt,2,2,1                                 _j2`#|oG  
vsel,s,volu,,3                             kN) pi "  
vatt,3,3,1                                 4C[gW  
vsel,s,volu,,4,6,1                         ;kFp)*i  
vatt,4,4,1                           p,0 \NUC  
allsel,all                                     #,9|Hr%  
!划分网格                                   @6tx5D?  
mshkey,1                                 +I~?8*  
mshape,1,3D                               xAO]u[J  
vmesh,all                                     'qD9kJ`  
da,11,all                                 tw(2V$J  
da,15,all                                 a3)#tt=rA  
da,19,all                                 qqAsh]Z  
da,23,all                                 tO0+~Wm  
da,26,all                                 ejQCMG7  
da,27,all                                 JjCf<ktE.  
da,30,all                                 D|BN_ai9  
da,31,all                                 b$eJH  
!施加重力                                 -Tt}M#W  
allsel,all                               [/kO>  
acel,,9.8,0                               JU6np4  
!进入求解器                               0CR;t`M@  
/solu                                     t2gjhn^p  
antype,static Qp54(`  
nlgeom,on                             GhT7:_r~  
nsubst,500,500,200                                                                                                                                                         .edZKmC6  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定     ^ Q}1&w%  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 RAxz+1JT  
allsel  

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

给你改了半天也没回来看 jJvNN -^  
更没有感谢 m0"\3@kB  
FENT