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

后处理中，发现自重应力场很混乱，明显有问题，见图。 V:27)]q  
请问高人如何解决？万分感谢！ .l|$dE/E  
命令如下： ExM,g'7  
finish !+njS  
/CLEAR,start DJ%PWlK5  
/fILEname,slope analysis |'.  
/units,si h$=
2p5'-  
/pnum,kp,1 
8[>zG
2  
/pnum,line,1 W`&hp6Jq  
/number,0 W l16`9  
/PREP7 -DCbko  
*afun,deg yBRC*0+Vy  
!定义单元类型 U3kyraj  
et,1,solid45 7rPF$ \#  
!定义材料属性 8] ikygt"  
!边坡围岩材料1的属性 J=L5=G7(  
mp,ex,1,6e7 bxWa oWE0  
mp,nuxy,1,0.29 +O5hH8<&b  
mp,dens,1,23 V+~Nalm O  
tb,dp,1 +>9Q/E  
tbdata,1,3e4 L]Mo;kT<Q  
tbdata,2,30 *qMY22X  
tbdata,3,0 2[CdZ(k]5  
!边坡围岩材料2的属性 6r_)sHf  
mp,ex,2,5e7 mqJ_W[y7  
mp,nuxy,2,0.21 -![|}p
X  
mp,dens,2,23.5 /@Zrq#o zx  
tb,dp,2 v3qA":(w+(  
tbdata,1,5e4 (ik\|y% A  
tbdata,2 rGkyGz8>  
tbdata,3,0 c)tfAD(N8x  
!边坡围岩材料3的属性 uGt-l4  
mp,ex,3,7e8 FC*[*  
mp,nuxy,3,0.34 >3_Gw4S*H  
mp,dens,3,25.1 oE~Bq/p  
tb,dp,3 .~}1+\~5  
tbdata,3,5e5 X jX2]  
tbdata,2,40 s{" 2L{,$  
tbdata,3,0 d7i]FV  
!边坡围岩材料4的属性 X7wKy(g  
mp,ex,4,1.5e9 rbpSg7}Q  
mp,nuxy,4,0.3 ty`DJO=Omj  
mp,dens,4,24.8 ;6wA"  
tb,dp,4 qw8Rlws%  
tbdata,4,8e4 n(|^SH4$b  
tbdata,2,40 g*"P:n71  
tbdata,3,0 M[u
A@  
!创建边坡模型线 ]~nKK@Rw  
k,1,101.95,69.87, :aQt
;C6Z>  
k,2,103.17,55.92 :yjFQ9^?&  
k,3,103.77,49.01 $kKjgQS(  
k,4,104.57,39.96
eY\yE"3  
k,5,105,35 >*n0n!vF  
k,6,140,35 yWya&|D9  
k,7,140,0 gO^gxJ'0t  
k,8,105,0 E!#WnSpnK  
k,9,0,0 -gWZwW/lD  
k,10,0,35 PT9*)9<L  
k,11,0,49.11 h}EPnC}  
k,12,0,63.68 AK#1]i~  
k,13,0,66.77 '=6\v!  
k,14,0,69.87 aATA9V  
l,1,2 e$pV%5=  
l,2,3 <9%R\_@$H  
l,3,4 j)GtEP<n#  
l,4,5 +]50DxflA  
l,5,6 IMfqiH)  
l,6,7 )Z VD+X  
l,7,8 `KQvJjA6  
l,8,9 TU7'J  
l,9,10 CA#,THty  
l,10,11 u4_9)P`]0  
l,11,12 g4@ lM"|S  
l,12,13 ow#1="G,=  
l,13,14 h-D}'R  
l,14,1 YDFyX){  
l,13,2 h*Pc=/p  
l,12,3 f`/x"@~H5  
l,11,4 w\brVnt  
l,10,5 B+0hzkPY  
l,5,8 ]d%8k}U  
!创建边坡面模型 eN~=*Mn(za  
al,14,1,15,13 3{h_&Gbo'D  
al,15,2,16,12 ,{q;;b9  
al,16,3,17,11 .}`Ix'
.  
al,17,4,18,10 6(e>P)  
al,18,19,8,9 l0hlM#  
al,5,6,7,19 _7)n(1h[3b  
vext,1,6,1,,,35,,, N&V`K0FU  
vsel,s,volu,,1 O<e{  
vatt,1,1,1 e*n@j  
vsel,s,volu,,2 W,-g=6,  
vatt,2,2,1 xp9pl[l  
vsel,s,volu,,3 M|[oaanY'  
vatt,3,3,1 t.'!`5G  
vsel,s,volu,,4,6,1 }#E[vRf  
vatt,4,4,1 =kqt  
allsel :Lug7bUVD  
!划分网格 X~i<g?]  
mshkey,1 "x /OIf  
mshape,1,2D _Y[bMuUb=  
vmesh,all Ip]KPrwp  
!固定面11，15，19，23，26，27，30，31 1>&]R=  
da,11,all O,A{3DAe0  
da,15,all  ^Va1f'g  
da,19,all Lu0x (/  
da,23,all F*K_+ ?m  
da,26,all $DUZ!zaH!  
da,27,all 4YX3+oS  
da,30,all &l[$*<P5V  
da,31,all w8D"CwS1Rx  
!施加重力 A_#DJJMm  
allsel,all lUiL\~Gq  
acel,,9.8,0 f>Jr|#k  
!进入求解器 ;xs"j-r/  
/solu *r% c  
antype,static O<;3M'y\  
nusubst,100 0,8okAH  
pred,on vFK<J Sk!  
nropt,full "2T#MO/  
nlgeom,on /Qk4  
insrch,on kn"(A.R  
outres,all,last f0aKlhEC  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 uc"P3,M  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 XEZF{lP  
allsel E\2%
E@0#  
solve PIpi1v*qz  
save wuJ4kW$  
finish ;{o|9x|  
/post1  !进入后处理 7 ^mL_SMj  
set,last lo!+f"7ym\  
plesol,s,y dmN&+t  
!plesol,u,y,0,1 AjgF6[B  
!plesol,s,xy,0,1 -8rjgB~."/  
!plesol,epto,xy,0,1 aCLqk'  
!pldisp,1                                      !绘制边坡模型变形图 mju>>
\9  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 Nl(3Xqov  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图 fe#\TNeQJ[  

自己顶个，等待高人！

是不是密度的单位拟弄错了 ~B(4qK1G  
怎么那么小

你这都收敛不了呀

finish                                     q,eVjtF  
/CLEAR,start                           aC.~&MxFC  
/fILEname,try                   9dUravC7  
/units,si                                 O,h;hQZ  
/pnum,kp,1                                 :|8M`18lZ  
/pnum,line,1                             <r`2)[7N  
/number,0                               zY!j:FT1HY  
/PREP7                                     FfPar:PHj  
*afun,deg                                 vVe';|8v  
!定义单元类型                             Ab"@714@  
et,1,solid45                               ~-J]W-n  
!定义材料属性                               >R!jB]5  
!边坡围岩材料1的属性                       sAD}#Zw$  
mp,ex,1,6e7                               )i^<r;_z  
mp,nuxy,1,0.29                             vv+z'(l  
mp,dens,1,23                               H_X [t*2  
tb,dp,1                                   w{@o^rs  
tbdata,1,3e4,30,0                                                             xZwLlY  
!边坡围岩材料2的属性                     hUMf"=q+  
mp,ex,2,5e7                               |! E)GahM  
mp,nuxy,2,0.21                             :'l^kSP_*C  
mp,dens,2,23.5                             NI [ pp`  
tb,dp,2                                   C-MjJ6D<  
tbdata,1,5e4,0,0                                                               zvH8^1yzG  
!边坡围岩材料3的属性                     4'A!; ]:  
mp,ex,3,7e8                               2=`o
_<P'"  
mp,nuxy,3,0.34                               l6 H|PR{  
mp,dens,3,25.1                             \(Y\|zC'0$  
tb,dp,3                           {I#]@,  
tbdata,1,5e5,40,0                                                               \EtQ5T*u  
!边坡围岩材料4的属性                       EbE-}>7OO  
mp,ex,4,1.5e9                             Rt!FPoN,y  
mp,nuxy,4,0.3                             m6CI{Sa](l  
mp,dens,4,24.8                             iJ3e1w$  
tb,dp,4                                   s<eb;Z2D  
tbdata,1,8e4,40,0                                                           C$D-Pt"+  
!创建边坡模型线 5                           ?9\EN|O^  
k,1,101.95,69.87                           a
(b#  
k,2,103.17,55.92             lqZ5?BD1  
k,3,103.77,49.01                         #[~pD:qqM  
k,4,104.57,39.96                         Zk"eA'"\  
k,5,105,35                               /}  WDU  
k,6,140,35                               7Vo$(kj  
k,7,140,0                               h+&OQ%e=8  
k,8,105,0                                   `FTy+8mw  
k,9,0,0                                   DBD%6o>]K  
k,10,0,35                                   &NoS=(s,  
k,11,0,49.11                               8UyMVY  
k,12,0,63.68                               ?!cvf{a  
k,13,0,66.77                               +M$Q =6/  
k,14,0,69.87                               QpA/SmJ  
l,1,2                                   71gT.E  
l,2,3                                   E!l!OtFL  
l,3,4                                   $5<#n@  
l,4,5                                   $#S&QHyEe  
l,5,6                                   ,TN 2  
l,6,7                                   w6GyBo{2O_  
l,7,8                                   cm[&?  
l,8,9                                   Dq5j1m.  
l,9,10                                     $gy*D7  
l,10,11                                 X4E%2-m@'  
l,11,12                                 W!&'p
g  
l,12,13                                 f@DYN!Z_m  
l,13,14                                 48qV>Gwf  
l,14,1                                   \6<=$vD  
l,13,2                                   M .JoHH  
l,12,3                                   bc)~k:  
l,11,4                                   xt%7@/hiE  
l,10,5                       L3--r  
l,5,8                                     C=It* j55  
!创建边坡面模型                         7/f3Z1g  
al,14,1,15,13                     G) 7;;  
al,15,2,16,12                     S.m{eur!,E  
al,16,3,17,11                     *q".-u!D[  
al,17,4,18,10                     dEA6  
al,18,19,8,9                   @&:ar
 
al,5,6,7,19                               DV-;4AxxRq  
vext,1,6,1,,,35,,,                           "C SC  
vsel,s,volu,,1                             5b[jRj6  
vatt,1,1,1                                  4/1d&Sg  
vsel,s,volu,,2                           RyK\uv  
vatt,2,2,1                                 -~ Mb  
vsel,s,volu,,3                             5Z\#0":e  
vatt,3,3,1                                 ws|;`  
vsel,s,volu,,4,6,1                         .#Z%1U%P.  
vatt,4,4,1                           #9xd[A: N  
allsel,all                                     Rh{zH~oZ  
!划分网格                                   7-T{a<g  
mshkey,1                                 A1#%`^W9  
mshape,1,3D                               s
q[iY  
vmesh,all                                     x`mN U  
da,11,all                                 {{MRELipW  
da,15,all                                 tj4VWJK  
da,19,all                                 dhr3,&+T2  
da,23,all                                 {(wHPzq  
da,26,all                                 ac.Ms(D  
da,27,all                                 @$c\dvO  
da,30,all                                 W"'iIh)z `  
da,31,all                                 <$~mE9a6  
!施加重力                                 i Ae<&Ms  
allsel,all                               lM{ +!-G,  
acel,,9.8,0                               NchXt6$i9  
!进入求解器                               (B_\TdQ  
/solu                                     "xHgqgFyO  
antype,static ;)e2@'Agl  
nlgeom,on                             "WK{ >T  
nsubst,500,500,200                                                                                                                                                         o=?C&f{  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定     5HO9+i  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 QxOjOKAG  
allsel  

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

给你改了半天也没回来看 P2'DD 3   
更没有感谢 !0C^TCuG  
FENT