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

后处理中，发现自重应力场很混乱，明显有问题，见图。 lH.2H  
请问高人如何解决？万分感谢！ oDa{HP\O]W  
命令如下： 6Y^o8R  
finish ~l~g0J
 
/CLEAR,start |{T2|iJI  
/fILEname,slope analysis 8vK&d>  
/units,si PQ>JoRs  
/pnum,kp,1 -yeT$P&|  
/pnum,line,1 T!bu}KO  
/number,0 jLSZ#H  
/PREP7 x#D=?/~/Kv  
*afun,deg "RLb wm~  
!定义单元类型 7
}k8-:a%  
et,1,solid45 g:U ul4  
!定义材料属性 ^>fs  
!边坡围岩材料1的属性 s9iM hCu|  
mp,ex,1,6e7 iq$/6!t  
mp,nuxy,1,0.29 z{\.3G  
mp,dens,1,23 4KCxhJq  
tb,dp,1 q]x@q  
tbdata,1,3e4 zu#o<6E{  
tbdata,2,30 *rLs!/[Z_  
tbdata,3,0 pC6_ jIZ  
!边坡围岩材料2的属性 /7^~*  
mp,ex,2,5e7 }8HLyK,4  
mp,nuxy,2,0.21 YbKW;L&Ff  
mp,dens,2,23.5 .FUEF)  
tb,dp,2 ioggD  
tbdata,1,5e4 &,zeBFmc  
tbdata,2 :9}*p@  
tbdata,3,0 9\F^\h{  
!边坡围岩材料3的属性 -g0>>{M'  
mp,ex,3,7e8 jJ|;Nwm<[  
mp,nuxy,3,0.34 ]NCOi?Odx  
mp,dens,3,25.1 iwbjjQPr  
tb,dp,3 c,@6MeKHq  
tbdata,3,5e5 gA6C(##0  
tbdata,2,40 H< j+-u4b  
tbdata,3,0 s>5
 Z  
!边坡围岩材料4的属性 *q RQN+%  
mp,ex,4,1.5e9 |NI0zd  
mp,nuxy,4,0.3 @UKd0kxPN{  
mp,dens,4,24.8 Y*Pr

 
tb,dp,4 >KE(%9y~  
tbdata,4,8e4 5 F-Q&  
tbdata,2,40  X_lNnk  
tbdata,3,0 t`Sh!e  
!创建边坡模型线 nV,a|V5Xm  
k,1,101.95,69.87, SZ_hGD0  
k,2,103.17,55.92 <$ 5\^y,V  
k,3,103.77,49.01 DcOLK\  
k,4,104.57,39.96 o! 8X< o  
k,5,105,35 X\;y;pmRH  
k,6,140,35 !>Xx</iD1  
k,7,140,0 yC[}gHv  
k,8,105,0 gnQd#`  
k,9,0,0 $Ix^Rm9c  
k,10,0,35 gisZmu0  
k,11,0,49.11 <*[(t;i  
k,12,0,63.68 J8jbtL O'  
k,13,0,66.77 ps2C8;zT  
k,14,0,69.87 n3(HA  
l,1,2 PF.HYtZqK  
l,2,3 +mJAIjH  
l,3,4 KnuqU2< {  
l,4,5 (f)QEho7  
l,5,6 B-RaAiE@  
l,6,7 iY="M_kQ_  
l,7,8 8:f(PN  
l,8,9 u%FA.  
l,9,10 hYLu   
l,10,11 fA8 ,wy|>  
l,11,12 s SDBl~g  
l,12,13 ?IK[]=!  
l,13,14 %n^]1R#  
l,14,1 OA_ %%A;o  
l,13,2 k@?<Aw8_X  
l,12,3 B#=dz,}  
l,11,4 F J)la9  
l,10,5 ="V6z$N  
l,5,8 >">-4L17m  
!创建边坡面模型 .L}ar7  
al,14,1,15,13 C`fQ` RL\  
al,15,2,16,12 ujaaO6oZ7  
al,16,3,17,11 E11"uWk`  
al,17,4,18,10 oZQu&O'  
al,18,19,8,9 Lr`yl$6  
al,5,6,7,19 Yv>% 5`  
vext,1,6,1,,,35,,, 1'ZBtX~A  
vsel,s,volu,,1 um/iK}O  
vatt,1,1,1 zJPzI{-w|  
vsel,s,volu,,2 !(/dbHB  
vatt,2,2,1 sPut@4[S  
vsel,s,volu,,3 vV%w#ULxE~  
vatt,3,3,1 [L:,A{rve  
vsel,s,volu,,4,6,1 O%OeYO69  
vatt,4,4,1 E;yP.<PW  
allsel )gZ yW  
!划分网格 uKK+V6}!kj  
mshkey,1 env]*gx+=  
mshape,1,2D sq_ f[!  
vmesh,all $aEL>,X  
!固定面11，15，19，23，26，27，30，31 {a(TT)d  
da,11,all BV|LRB}G  
da,15,all GujmBb  
da,19,all PAJt M  
da,23,all |(]XZ!{  
da,26,all WmeV[iI  
da,27,all +5voAx!  
da,30,all "5=Gu1  
da,31,all d4~!d>{n|c  
!施加重力 "~XAD(T6  
allsel,all Vf0m7BJc3  
acel,,9.8,0 t\ ym4`"  
!进入求解器 `j9 ;9^  
/solu *4]}_ .rG#  
antype,static >o.4sN@  
nusubst,100 Pa+%H]vB  
pred,on V3 2F  
nropt,full b:2#3
;)  
nlgeom,on v#TU7v?~  
insrch,on  t
mKHT  
outres,all,last 13kl\<6  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 )m|)cLT&  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 RGsgT^  
allsel bZLY#g7L"  
solve nH_M#  
save =A{F&:+a]  
finish @kd`9Yw  
/post1  !进入后处理 ,?S1e#  
set,last XkDIP4v%  
plesol,s,y /V0[Urc@  
!plesol,u,y,0,1 8u6*;*o  
!plesol,s,xy,0,1 MaN6bM  
!plesol,epto,xy,0,1 !cW6dc^  
!pldisp,1                                      !绘制边坡模型变形图 Qhy!:\&1  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 %py3fzg  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图 W^k95%zBM  

自己顶个，等待高人！

是不是密度的单位拟弄错了 9#TD1B/  
怎么那么小

你这都收敛不了呀

finish                                     )
w0x{_  
/CLEAR,start                           Gsv<Rjj:  
/fILEname,try                   arDl2T,igF  
/units,si                                 PQDLbSe)\  
/pnum,kp,1                                 `mzlOB  
/pnum,line,1                             o>_})WM1[  
/number,0                               R|n  
/PREP7                                     =n
HKTB>  
*afun,deg                                 i+h*<){X  
!定义单元类型                             k?Z:=.YW  
et,1,solid45                               V_:`K$  
!定义材料属性                               l3sF/zkH  
!边坡围岩材料1的属性                       7Y9#y{v1  
mp,ex,1,6e7                               ao2^3e  
mp,nuxy,1,0.29                             O;~e^ <*  
mp,dens,1,23                               is4}s,]$6  
tb,dp,1                                   hE|P|0U,n  
tbdata,1,3e4,30,0                                                             sqrLys_S  
!边坡围岩材料2的属性                     IplOXD  
mp,ex,2,5e7                               4p,:}h  
mp,nuxy,2,0.21                             E +_n@t"  
mp,dens,2,23.5                             S3\NB3@qC&  
tb,dp,2                                   \IE![=p\w  
tbdata,1,5e4,0,0                                                               5Jbwl$mZ  
!边坡围岩材料3的属性                     3]l)uoNt/  
mp,ex,3,7e8                               G IN|cv=  
mp,nuxy,3,0.34                               x="Wqcnj{  
mp,dens,3,25.1                             oT-gZedW(  
tb,dp,3                           f`n4'dG  
tbdata,1,5e5,40,0                                                               TI'v /=;)  
!边坡围岩材料4的属性                       s0/O/G?  
mp,ex,4,1.5e9                             QWk3y"5n<  
mp,nuxy,4,0.3                             UJL2IF-x  
mp,dens,4,24.8                             Ipk;Nq  
tb,dp,4                                   }:0_%=)N<  
tbdata,1,8e4,40,0                                                           f' ?/P~[  
!创建边坡模型线 5                           hZx&j{  
k,1,101.95,69.87                           p?i.<Z  
k,2,103.17,55.92             _~z oMdT!  
k,3,103.77,49.01                         ( zWBrCX  
k,4,104.57,39.96                         IMkE~0x4</  
k,5,105,35                               |NuMDVd+s  
k,6,140,35                               +/ZIs|B4,z  
k,7,140,0                               XOvJlaY)'.  
k,8,105,0                                   qYC&0`:H  
k,9,0,0                                   7%y$^B7{  
k,10,0,35                                   zmo2uUEd  
k,11,0,49.11                               +>"s)R43  
k,12,0,63.68                               gQ_<;'m)2  
k,13,0,66.77                               >v]^nJl  
k,14,0,69.87                               7*l$i/!  
l,1,2                                   ubwM*P  
l,2,3                                   Q;]JVT1  
l,3,4                                   3`bQ0-D;  
l,4,5                                   *?FVLE  
l,5,6                                   :W.H#@'(  
l,6,7                                   l7|z]v-  
l,7,8                                   ^

%r6+ey  
l,8,9                                   Y4rxnXGw  
l,9,10                                     "`>6M&`U  
l,10,11                                 /eV)5`V  
l,11,12                                 32wtN8kx  
l,12,13                                
@G4
Z  
l,13,14                                 /j~~S'sw  
l,14,1                                   NiZfaC6V  
l,13,2                                   "w:h  
l,12,3                                   @7Ec(]yp  
l,11,4                                   ^Hx}.?1  
l,10,5                       @IL04' \  
l,5,8                                     (VeK7cU  
!创建边坡面模型                         -8 uS#  
al,14,1,15,13                     B!wN%>U  
al,15,2,16,12                     $94lF~  
al,16,3,17,11                     B4y_{V  
al,17,4,18,10                     2FMmANH0ev  
al,18,19,8,9                   0t7N yKU  
al,5,6,7,19                               i!a!qE.1  
vext,1,6,1,,,35,,,                           L^JU{\C  
vsel,s,volu,,1                             dla_uXtM6  
vatt,1,1,1                                 C~&E7w  
vsel,s,volu,,2                           ?w:\0j5~  
vatt,2,2,1                                 ?MH=8Cl1w  
vsel,s,volu,,3                             2u]G]:ml  
vatt,3,3,1                                 Y/sav;  
vsel,s,volu,,4,6,1                         j9fBl:Fr  
vatt,4,4,1                           f Fi=/}  
allsel,all                                     bJ:5pBJ3  
!划分网格                                   x{pj`'J)  
mshkey,1                                 .#&)%}GC  
mshape,1,3D                               sl$6Zv-l%0  
vmesh,all                                     T\gs  
da,11,all                                 F@<cp ?dR  
da,15,all                                 sVh)Ofn  
da,19,all                                 (dg,w*t'  
da,23,all                                 #N][-i  
da,26,all                                 g:!U,<C^a  
da,27,all                                 K@p9_K8  
da,30,all                                 02,t  
da,31,all                                 -53c0g@X  
!施加重力                                 &b (*  
allsel,all                               }1 O"?6  
acel,,9.8,0                               Qb {[xmc  
!进入求解器                               fq>{5ODO  
/solu                                     ;MQl.?vj  
antype,static Xm,fyk>  
nlgeom,on                             H'i\N?VL  
nsubst,500,500,200                                                                                                                                                         r5gqRh}+  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定     G ]h  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 N8nt2r<h  
allsel  

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

给你改了半天也没回来看 >-V632(/{o  
更没有感谢 5`3x(=b  
FENT