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

后处理中，发现自重应力场很混乱，明显有问题，见图。 Bv/v4(G5g  
请问高人如何解决？万分感谢！ @n})oAC,  
命令如下： BDcl1f T  
finish \d v9:X$  
/CLEAR,start b%pLjvU  
/fILEname,slope analysis H#Vs3*VK  
/units,si  BeP0lZ  
/pnum,kp,1 ?K}/b[[0v  
/pnum,line,1 *UEo&B2+  
/number,0 m#8mU,7  
/PREP7 i\,#Z
!  
*afun,deg Xf$,ra"  
!定义单元类型 q}+9$v  
et,1,solid45 K _y;<a]  
!定义材料属性 a*??!  
!边坡围岩材料1的属性 bC"#.e  
mp,ex,1,6e7 o'Po<I  
mp,nuxy,1,0.29 ="Dmfy7  
mp,dens,1,23 =1LrU$\  
tb,dp,1 H%sbf& gi  
tbdata,1,3e4 !_zmm$bR  
tbdata,2,30 {AQ3y,sh  
tbdata,3,0 _L&C4 <e'  
!边坡围岩材料2的属性 , gz:2UY#  
mp,ex,2,5e7 MbjH\XRB  
mp,nuxy,2,0.21 =X>?Y,  
mp,dens,2,23.5 @'~v~3 $S  
tb,dp,2 @XB/9!  
tbdata,1,5e4 6i>xCb  
tbdata,2 3<?(1kSo>>  
tbdata,3,0 _%"/I96'  
!边坡围岩材料3的属性 M-0BQs`N  
mp,ex,3,7e8 K{"(|~=U  
mp,nuxy,3,0.34 wYNh0QlBH  
mp,dens,3,25.1 ].`i`.T  
tb,dp,3 .}]5y4UQ.  
tbdata,3,5e5 G9^!= v@  
tbdata,2,40 Q- j+#NGc  
tbdata,3,0 P>jlFm  
!边坡围岩材料4的属性 lZE x0  
mp,ex,4,1.5e9 6Pd;I,k  
mp,nuxy,4,0.3 Pm V:J9  
mp,dens,4,24.8 [%YCupr#  
tb,dp,4 >_tn7Z0L  
tbdata,4,8e4 QBDi;Xzb+  
tbdata,2,40 r@k&1*&  
tbdata,3,0 9<,\+}^{  
!创建边坡模型线 aq[kKS`  
k,1,101.95,69.87, 8CwgV  
k,2,103.17,55.92 <C]s\"o-`  
k,3,103.77,49.01 J`V7FlM  
k,4,104.57,39.96 ={8ClUV#  
k,5,105,35 )*$'e<?`  
k,6,140,35 r1$ O<3\  
k,7,140,0 /R@eOl}D  
k,8,105,0 9?@M Zh  
k,9,0,0 * UBU?  
k,10,0,35 E3sl"d;~  
k,11,0,49.11 f~8Xue,l"  
k,12,0,63.68 o{nBtxZ"  
k,13,0,66.77  :\gdQG  
k,14,0,69.87 JtvAi\52$  
l,1,2 } n_9d.  
l,2,3 OG/R6k.  
l,3,4 jD'\\jAUdm  
l,4,5 jSpmE  
l,5,6 I}$Y[Jve  
l,6,7 F0]NtKaH  
l,7,8 AfFFu\  
l,8,9 +-K-CXt  
l,9,10 2NqlE  
l,10,11 FL"7u2rh,  
l,11,12 "J3@Z,qW  
l,12,13 'I\bz;VT  
l,13,14 c9r, <TR9  
l,14,1 qeHb0G  
l,13,2 bp;b;f>  
l,12,3 0ir]  
l,11,4 @6(4}&sEdm  
l,10,5 6# ,2  
l,5,8 {3@/@jO?  
!创建边坡面模型 9;sebqC?  
al,14,1,15,13 h@@2vs2  
al,15,2,16,12 i|y8n7c  
al,16,3,17,11 /vBOf;L  
al,17,4,18,10 wHAoO#`wn5  
al,18,19,8,9 +{4ziqYj  
al,5,6,7,19 |^8l8u  
vext,1,6,1,,,35,,, DXsp 2
  
vsel,s,volu,,1 N'=8Dj  
vatt,1,1,1 M&:[3u-  
vsel,s,volu,,2 Rs8^
27  
vatt,2,2,1 z3[ J>  
vsel,s,volu,,3 hL+)XJu^J  
vatt,3,3,1 oZQ%
P  
vsel,s,volu,,4,6,1 oGu-:X=`9  
vatt,4,4,1 iQ~;to;Y  
allsel N/4`afiV.  
!划分网格 'n6D3Vse  
mshkey,1 3#!}W#xv  
mshape,1,2D +.! F]0ju  
vmesh,all 4w<U%57  
!固定面11，15，19，23，26，27，30，31 ncuqo'r  
da,11,all [daUtKz  
da,15,all U\*}}   
da,19,all ,Tvfn`;(  
da,23,all ^M0e0  
da,26,all j(AN]g:  
da,27,all g+bc4eU  
da,30,all {r yv7G  
da,31,all 9o@5:.b<j  
!施加重力 6F(;=iY8  
allsel,all A9`& Wnw?  
acel,,9.8,0 \_bk+}WJ]s  
!进入求解器 @A GM=v  
/solu z'z_6]5  
antype,static +fBbW::R^  
nusubst,100 u!S^lV@  
pred,on ABuK`(f.  
nropt,full ^_0zO$z,  
nlgeom,on nYhI0q  
insrch,on s{#rCc)  
outres,all,last [K1RP.  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 =]X_wA;%  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 Zi<Y?Vm/,O  
allsel P-[6'mw`  
solve h mC.5mY  
save yU<T_&
M  
finish ]q~_  
/post1  !进入后处理 !YYI{BJ7:N  
set,last TmZsC5  
plesol,s,y 7:awUoV8f  
!plesol,u,y,0,1 I?<5 %  
!plesol,s,xy,0,1 c5Fl:=h  
!plesol,epto,xy,0,1 { }Q!./5  
!pldisp,1                                      !绘制边坡模型变形图 !6{; z/Hy  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 SM;UNIRVE  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图 o/??w:'  

自己顶个，等待高人！

是不是密度的单位拟弄错了 'Z;8-1M?O  
怎么那么小

你这都收敛不了呀

finish                                     p<[gzmU9\b  
/CLEAR,start                           r%=}e++^%  
/fILEname,try                   [}ayaXXQ5  
/units,si                                 (D5sJ$&E@\  
/pnum,kp,1                                 Xg4iH5!E  
/pnum,line,1                             >v?&&FhHK<  
/number,0                               uPhL?s{  
/PREP7                                     qp##>c31X  
*afun,deg                                 :)*+aS"  
!定义单元类型                             >(.GIR  
et,1,solid45                               v^0*{7N'  
!定义材料属性                               E&z`BPd  
!边坡围岩材料1的属性                       C 9{8!fYp  
mp,ex,1,6e7                               k.5(d.*(  
mp,nuxy,1,0.29                             `W=3_  
mp,dens,1,23                               ahoXQ8c:\}  
tb,dp,1                                   XK+" x!   
tbdata,1,3e4,30,0                                                             NB^+Hcb$  
!边坡围岩材料2的属性                     r `;_ #&b  
mp,ex,2,5e7                               _/c1b>kcso  
mp,nuxy,2,0.21                             7$zeRYD+  
mp,dens,2,23.5                             #lDW?  
tb,dp,2                                   <6rc8jYz  
tbdata,1,5e4,0,0                                                               s;!_'1pi@  
!边坡围岩材料3的属性                     'X%5i2  
mp,ex,3,7e8                               j0wpaIp  
mp,nuxy,3,0.34                               )hy(0 D  
mp,dens,3,25.1                             ;n\$'"K&;  
tb,dp,3                           3bN]2\   
tbdata,1,5e5,40,0                                                               Q\moR
^>  
!边坡围岩材料4的属性                       uyqu n@
q  
mp,ex,4,1.5e9                             2j=i\B  
mp,nuxy,4,0.3                             |O';$a1S  
mp,dens,4,24.8                             iVSN>APe  
tb,dp,4                                   gVscdg5  
tbdata,1,8e4,40,0                                                           `m")v0n3  
!创建边坡模型线 5                           GZt L-  
k,1,101.95,69.87                           OaH1xZNOC`  
k,2,103.17,55.92            
l^k/Y ]  
k,3,103.77,49.01                         Q1[3C(  
k,4,104.57,39.96                         ABaK60.O[O  
k,5,105,35                               fw3P?_4;*  
k,6,140,35                               }N0$DqP  
k,7,140,0                               Hy~+|hLvh  
k,8,105,0                                   YQ-!>3/)-  
k,9,0,0                                   }nx5  
k,10,0,35                                   (~IoRhp^  
k,11,0,49.11                               '>:c:Tewy  
k,12,0,63.68                               ~W..P:wG5  
k,13,0,66.77                               zB68%  
k,14,0,69.87                               Da3Z>/S  
l,1,2                                   ascY E  
l,2,3                                   "xduh3/~=  
l,3,4                                   4 \Ig<C9  
l,4,5                                   p}e1!q;N  
l,5,6                                   +5i~}Q!  
l,6,7                                   [X"pOz  
l,7,8                                   3YO%$  
l,8,9                                   ZBuh(be  
l,9,10                                     M 8NWQ^Y  
l,10,11                                 xf4CM,Z7(  
l,11,12                                 UbDpSfub  
l,12,13                                 j;yKL-ycB  
l,13,14                                 nh'TyUd!  
l,14,1                                   dMw}4c3E  
l,13,2                                   ]tx/t^&/\u  
l,12,3                                   ,_M  
l,11,4                                   =9!|%j  
l,10,5                       i9+qU  
l,5,8                                     OpD%lRl  
!创建边坡面模型                         H3>49;`  
al,14,1,15,13                     >[0t@Tu,D  
al,15,2,16,12                     IgyoBfj\d  
al,16,3,17,11                     ,^1B"#0{C<  
al,17,4,18,10                     OB4nE}NO  
al,18,19,8,9                   /~1Ew  
al,5,6,7,19                               aoHAB<.C  
vext,1,6,1,,,35,,,                           ty\F~]Oo  
vsel,s,volu,,1                             1smKU9B2)  
vatt,1,1,1                                 wx*03(|j;  
vsel,s,volu,,2                           N8KH.P+  
vatt,2,2,1                                 Ib(G!oO:E-  
vsel,s,volu,,3                             !nDiAjj  
vatt,3,3,1                                 f\vMdY  
vsel,s,volu,,4,6,1                         J3QL%#  
vatt,4,4,1                           Am@:<J  
allsel,all                                     I~F]e|Ehqr  
!划分网格                                   S;^'Ek"Z.  
mshkey,1                                 
br,xwc  
mshape,1,3D                               ('z:XW96  
vmesh,all                                     QAzwNXE+  
da,11,all                                 z'01V8e  
da,15,all                                 eq(1'?7]`G  
da,19,all                                 z7_h$v  
da,23,all                                 zQ#2BOx1  
da,26,all                                 7 -S?U~s  
da,27,all                                 A~
X| vW  
da,30,all                                 0.kC|  
da,31,all                                 [NoOA  
!施加重力                                 <nU8.?\?~  
allsel,all                               {,B.OM)J  
acel,,9.8,0                               -&EU#Wqh  
!进入求解器                               *cP(3n3]R  
/solu                                     : "^/?Sd  
antype,static :/R>0n,  
nlgeom,on                             h?Lp9VF  
nsubst,500,500,200                                                                                                                                                         |OiM(E(  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定     VHm.uL_UW  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 2. q\!V}yQ  
allsel  

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

给你改了半天也没回来看 sI6I5  
更没有感谢 S$=caZ?  
FENT