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

后处理中，发现自重应力场很混乱，明显有问题，见图。 /6{P ?)]pE  
请问高人如何解决？万分感谢！ HT%'dZ1  
命令如下： w.o>G2u  
finish 3j7Na#<tL3  
/CLEAR,start  NIh?2w"\  
/fILEname,slope analysis IgyoBfj\d  
/units,si b!4Z~d0=  
/pnum,kp,1 s+<`iH9Hm  
/pnum,line,1 ?#~km0~F)  
/number,0 ){I!orQ  
/PREP7 7U1^=Y@t}  
*afun,deg 8G )O,F7z  
!定义单元类型 1:;S6
{oQ  
et,1,solid45 [pxC3{|d$  
!定义材料属性 oJ.5! Kg  
!边坡围岩材料1的属性 ?1.WF}X'  
mp,ex,1,6e7
7CwQmVe+  
mp,nuxy,1,0.29  SH6+'7  
mp,dens,1,23 s H(io  
tb,dp,1 JKTn  
tbdata,1,3e4 ,<s/K  
tbdata,2,30 J3QL%#  
tbdata,3,0 )}$]~ f4R  
!边坡围岩材料2的属性 U ATF}x  
mp,ex,2,5e7 )]tf|Mbu  
mp,nuxy,2,0.21 [x{Ai( /T^  
mp,dens,2,23.5 qYP;`L}o#  
tb,dp,2 Ws/\lD  
tbdata,1,5e4 ('z:XW96  
tbdata,2 {\vcwMUzZ  
tbdata,3,0 "t)$4gERK  
!边坡围岩材料3的属性 VOSq%hB  
mp,ex,3,7e8 Y !%2vOt  
mp,nuxy,3,0.34 R~"&E#C  
mp,dens,3,25.1 m*'hHt n  
tb,dp,3 zQ#2BOx1  
tbdata,3,5e5 Q2Dh(  
tbdata,2,40 }eAV8LU  
tbdata,3,0 'Px}#f0IR  
!边坡围岩材料4的属性 /hS
Em.<  
mp,ex,4,1.5e9 HChlkj'7w0  
mp,nuxy,4,0.3 Oq*;GR(Q  
mp,dens,4,24.8 |)*9BN  
tb,dp,4
(B0QBDj!  
tbdata,4,8e4 |,Kk#`lW<f  
tbdata,2,40 n
{I1ZlEeh  
tbdata,3,0 *cP(3n3]R  
!创建边坡模型线 
F4]=(T  
k,1,101.95,69.87, q.kDx_  
k,2,103.17,55.92 7g>|e  
k,3,103.77,49.01 t{-*@8Ke  
k,4,104.57,39.96 5t,X;  
k,5,105,35 u<JkP <"S  
k,6,140,35 VHm.uL_UW  
k,7,140,0 BG8`B'i  
k,8,105,0 @TJ2 |_s6]  
k,9,0,0 &~:+2  
k,10,0,35 #.
{ddY{  
k,11,0,49.11 $ftxid8  
k,12,0,63.68 *]:J@KGf  
k,13,0,66.77 WH_ W:  
k,14,0,69.87 \0n<6^y  
l,1,2 pIHpjx  
l,2,3 O>"T*   
l,3,4 88KQ) NU  
l,4,5 '{w[).c.  
l,5,6 3b?8<*  
l,6,7 ~gdnD4[G  
l,7,8 b!<)x}-t>  
l,8,9 WD@v<Wx)  
l,9,10 y!c<P,Lt3f  
l,10,11 xW|8-q  
l,11,12 M%Zh{  
l,12,13 9aBz%* xo  
l,13,14 Van=dzG  
l,14,1 0Ibe~!EiQJ  
l,13,2 to
;^'#B  
l,12,3 @E4ya$A)F  
l,11,4 {>Hn:jW<.  
l,10,5 !u;gGgQF  
l,5,8 ++ZP X'|  
!创建边坡面模型 vNHvuwK  
al,14,1,15,13 EXsVZg"#  
al,15,2,16,12 wb>"'%  
al,16,3,17,11 @<3kj R?j  
al,17,4,18,10 }wZsM[
NDB  
al,18,19,8,9 y~'%PUN  
al,5,6,7,19 hkOFPt&  
vext,1,6,1,,,35,,, sI6I5  
vsel,s,volu,,1 R|92T*h  
vatt,1,1,1 S$=caZ?  
vsel,s,volu,,2 -/:!AxIH  
vatt,2,2,1 O^n\lik  
vsel,s,volu,,3
,3 [FD9  
vatt,3,3,1 E|A~T7G=  
vsel,s,volu,,4,6,1 Ept=&mJPu  
vatt,4,4,1 d;nk>6<|  
allsel 7 lSR  
!划分网格 ?pWda<&  
mshkey,1 0oe2X1.%  
mshape,1,2D %~:@
}C%A  
vmesh,all WRrg5&._q  
!固定面11，15，19，23，26，27，30，31 V4&a+MJ@  
da,11,all `!xI!Y\  
da,15,all P%;lHC #i  
da,19,all mU]pK5  
da,23,all oB(9{6@N  
da,26,all LrT? ]o  
da,27,all L2c\i  
da,30,all ;ojJXH~$}  
da,31,all FX!Qd&kl1  
!施加重力 {v"Y!/ [z  
allsel,all A2!pbeG  
acel,,9.8,0 Yx&d\/9  
!进入求解器 m%nRHT0KAf  
/solu
"xNP"S  
antype,static < lUpvr  
nusubst,100 -p"}K~lt:  
pred,on l.(|&U~  
nropt,full ,/qY 9eh  
nlgeom,on 'nS>'yYH#  
insrch,on &dPI<HlM  
outres,all,last >2v<;.  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 m/q`k  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 {Iz"]Wh<f  
allsel "
ZL_  
solve }~p%e2<  
save qDcoccEf  
finish k;K-6<^h  
/post1  !进入后处理 1\)C;c
,  
set,last $rJgBN  
plesol,s,y >S8 n8U  
!plesol,u,y,0,1 n5#QQk2  
!plesol,s,xy,0,1 Fkg%_v$  
!plesol,epto,xy,0,1 |h\A5_0_  
!pldisp,1                                      !绘制边坡模型变形图 a)L|kux;l  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 FOteNQTj  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图 >2F9Tz,3  

自己顶个，等待高人！

是不是密度的单位拟弄错了 J*_^~t  
怎么那么小

你这都收敛不了呀

finish                                     2b
G92  
/CLEAR,start                           FS!9 j8  
/fILEname,try                   C,%Dp0  
/units,si                                 0:-z+`RHE  
/pnum,kp,1                                 4"\x#  
/pnum,line,1                             m2/S(f  
/number,0                               |onLJY7)  
/PREP7                                    
H=w6  
*afun,deg                                 Vk2%yw>  
!定义单元类型                             Spu;  
et,1,solid45                               .wUnN8crQ  
!定义材料属性                               zo("v*d*q  
!边坡围岩材料1的属性                       4.Fh4Y:$'  
mp,ex,1,6e7                               neE Zw#(Z  
mp,nuxy,1,0.29                             `gF`Sgz  
mp,dens,1,23                               $5\!ws<cZ  
tb,dp,1                                   bSn={O"M  
tbdata,1,3e4,30,0                                                             lHPnAaue@  
!边坡围岩材料2的属性                     nE?:nJ|%E  
mp,ex,2,5e7                               Y+kuj],h  
mp,nuxy,2,0.21                             [P0c,97_ H  
mp,dens,2,23.5                             Qmb+%z  
tb,dp,2                                   8k)*f+1o  
tbdata,1,5e4,0,0                                                               2 E?]!9T~|  
!边坡围岩材料3的属性                     SL`; `//  
mp,ex,3,7e8                               t,'J%)j  
mp,nuxy,3,0.34                               >DbG )0|  
mp,dens,3,25.1                             a8v\H8@X  
tb,dp,3                           wkx#WC  
tbdata,1,5e5,40,0                                                               6`%|-o :  
!边坡围岩材料4的属性                       :fnJp9c  
mp,ex,4,1.5e9                             *M$mAy<  
mp,nuxy,4,0.3                             ,izp^,`  
mp,dens,4,24.8                             \}%_FnP0ZU  
tb,dp,4                                   `Y+R9bd  
tbdata,1,8e4,40,0                                                           >;F}>_i  
!创建边坡模型线 5                           X$G:3uoN  
k,1,101.95,69.87                           1r*yYm'  
k,2,103.17,55.92             !I\eIV>0b  
k,3,103.77,49.01                         2pvby`P4  
k,4,104.57,39.96                         S4c-i2Rq  
k,5,105,35                               #xq3)B  
k,6,140,35                               G$)tp^%]  
k,7,140,0                               X\w["!B  
k,8,105,0                                   .r&CIL>  
k,9,0,0                                   '4D7:  
k,10,0,35                                   MmU`i ,z  
k,11,0,49.11                               v~W;&{  
k,12,0,63.68                               q2OF-.rE  
k,13,0,66.77                               q&O9W?E8dG  
k,14,0,69.87                               Q[%G`;e#  
l,1,2                                   Fe+(+ S  
l,2,3                                   MiKq|  
l,3,4                                   M+^ NF\  
l,4,5                                   O_GHvLO=  
l,5,6                                   $FDGHFM  
l,6,7                                   Wb#<ctM>  
l,7,8                                   B*E"yB\NV  
l,8,9                                   ,_/\pX0  
l,9,10                                     t-n'I/^5  
l,10,11                                 4E}Q<?UYSt  
l,11,12                                 H|,d`@U  
l,12,13                                 b6H7>x
 
l,13,14                                 DZLEx{cm  
l,14,1                                   LOu9#w"  
l,13,2                                   N})vrB;1  
l,12,3                                   +fKV/tSWi  
l,11,4                                   ttHRc!  
l,10,5                       TKbfZw  
l,5,8                                     KY0<N9{  
!创建边坡面模型                         pHye8v4fvi  
al,14,1,15,13                     T+t7/PwC;  
al,15,2,16,12                     U%Hcck'  
al,16,3,17,11                     S6C D
K:  
al,17,4,18,10                     UUM:*X  
al,18,19,8,9                   Fp"c {  
al,5,6,7,19                               @MoCEtt  
vext,1,6,1,,,35,,,                           pIug$Ke_%  
vsel,s,volu,,1                             ux*G*QZ  
vatt,1,1,1                                 <0m^b#hdG  
vsel,s,volu,,2                           xRO9o3  
vatt,2,2,1                                 k_Sm ep  
vsel,s,volu,,3                             =3rf}bl2  
vatt,3,3,1                                 K6)IBV;
 
vsel,s,volu,,4,6,1                         I2NMn5>  
vatt,4,4,1                           %Wm)  
allsel,all                                     a+CJJ3T-  
!划分网格                                   ??|,wIRz  
mshkey,1                                 rf 60'  
mshape,1,3D                               )7-mALyW  
vmesh,all                                     QNv5CQ&  
da,11,all                                 53.jx38xS  
da,15,all                                 S$P=;#r  
da,19,all                                 Tc>g+eS  
da,23,all                                 \}<J>R@  
da,26,all                                 qlg.\H:W~  
da,27,all                                 0r[a$p>`  
da,30,all                                 nB[B FVkU  
da,31,all                                 X+ybgB4(  
!施加重力                                 beT[7uVj_  
allsel,all                               7L6^IK  
acel,,9.8,0                               m;IKV,  
!进入求解器                               M0e&GR8<z>  
/solu                                     #,FXc~V  
antype,static YguW2R=6]  
nlgeom,on                             !}c\u  
nsubst,500,500,200                                                                                                                                                         a*_&[  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定      BPKrRex  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 gxeu2HG  
allsel  

材料3和4定义错误，应该是TBDATA,1,.....不是TBDATA,3 q[Ey!h)xq  
定义子步出错，命令流就是错的

给你改了半天也没回来看 V' i@N  
更没有感谢 <h<_''+  
FENT