论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录 或 注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 2779阅读
  • 6回复

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

上一主题 下一主题
离线51200605
 

发帖
334
土币
46
威望
1025
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2009-10-20
后处理中,发现自重应力场很混乱,明显有问题,见图。 &E?TR A# E  
请问高人如何解决?万分感谢! sk%Xf,  
命令如下: X7~^D[ X  
finish /Kd9UQU  
/CLEAR,start cb3Q{.-.#  
/fILEname,slope analysis L*8U.{NY  
/units,si /&i6vWMhP  
/pnum,kp,1 mG*ER^Y@D  
/pnum,line,1 De-hHY{>  
/number,0 ez-jVi-Fi  
/PREP7 gX%"Ki7.  
*afun,deg w-j^jU><3  
!定义单元类型 @zC p/fo3  
et,1,solid45 6"L,#aKm^  
!定义材料属性 ^\f1zg9I  
!边坡围岩材料1的属性 : MEB] }  
mp,ex,1,6e7 G`e!WvC  
mp,nuxy,1,0.29 QM) ob  
mp,dens,1,23 R<<U(.E  
tb,dp,1 u*S-Pji,x  
tbdata,1,3e4 zk;'`@7  
tbdata,2,30 U[R[VY7  
tbdata,3,0 5Ic'6AIz  
!边坡围岩材料2的属性 f=EWr8mno  
mp,ex,2,5e7 @* <`*W  
mp,nuxy,2,0.21 Ql1J?9W  
mp,dens,2,23.5 /prR;'ks  
tb,dp,2 kf:Nub+h t  
tbdata,1,5e4 w7%.EA{N  
tbdata,2 A"}Ib'  
tbdata,3,0 1RgERj  
!边坡围岩材料3的属性 &}rmDx  
mp,ex,3,7e8 jhJ'fI  
mp,nuxy,3,0.34 Z}AhDIw!G  
mp,dens,3,25.1 KY'x;\0 g  
tb,dp,3 <r1/& RW,  
tbdata,3,5e5 &v/>P1Z G  
tbdata,2,40 W%wc@.P  
tbdata,3,0 e~ZxDAd  
!边坡围岩材料4的属性  ?8/T#ox  
mp,ex,4,1.5e9 t?(fDWd|-  
mp,nuxy,4,0.3 B'v~0Kau  
mp,dens,4,24.8 W; zzc1v  
tb,dp,4 3 ,f3^A  
tbdata,4,8e4 ?u4t;  
tbdata,2,40 xxQgX~'x  
tbdata,3,0 uqsVq0H  
!创建边坡模型线 V<i_YLYmJe  
k,1,101.95,69.87, b[2 #t  
k,2,103.17,55.92 W]oILL"d  
k,3,103.77,49.01 3Fg{?C_l  
k,4,104.57,39.96  8+,I(+  
k,5,105,35 'Ul^V  
k,6,140,35 47=YP0r?>T  
k,7,140,0 lD#S:HX  
k,8,105,0 S]Qf p,  
k,9,0,0 g7;OZ#\  
k,10,0,35 UrmnHc>}c  
k,11,0,49.11 iBt<EM]U/  
k,12,0,63.68 1wg#4h43l  
k,13,0,66.77 4M )oA|1w  
k,14,0,69.87 u- }@^Y$M  
l,1,2 $vLGX>H  
l,2,3 ;L@p|]fu  
l,3,4 98rO]rg  
l,4,5 O>LqpZ  
l,5,6 RI3GAd  
l,6,7 KIGMWS^^  
l,7,8 JKF/z@Vbe\  
l,8,9 0F%/R^mw  
l,9,10 N[zR%(YS  
l,10,11 !"LFeqI$lr  
l,11,12 o}=c (u  
l,12,13 0O!A8FA0  
l,13,14 D=jtXQF  
l,14,1 |4j'KM;U  
l,13,2 ^4~?]5Y\  
l,12,3 bIXD(5y  
l,11,4 ]^0mh["  
l,10,5 7NMy1'-q  
l,5,8 ANRZQpnXQ  
!创建边坡面模型 }3/|;0j$  
al,14,1,15,13 LL_@nvu}M  
al,15,2,16,12 6n:oEXM>  
al,16,3,17,11 ;w<r/dK   
al,17,4,18,10 ILIv43QKM(  
al,18,19,8,9 O9P4r*prA  
al,5,6,7,19 ELZ@0,  
vext,1,6,1,,,35,,, 0<)Ep~!  
vsel,s,volu,,1 @x@wo9<Fc  
vatt,1,1,1 J(Fk@{!F.*  
vsel,s,volu,,2 XE$;Z'Qhjm  
vatt,2,2,1 F|&%Z(@a  
vsel,s,volu,,3 %%T?LRv  
vatt,3,3,1 4d8}g25C  
vsel,s,volu,,4,6,1 v:gdG|n"  
vatt,4,4,1 +&4@HHU{G  
allsel (XNd]G  
!划分网格 Sw.Kl 0M  
mshkey,1 (5l'?7  
mshape,1,2D iLO,XW?d v  
vmesh,all Rr0]~2R  
!固定面11,15,19,23,26,27,30,31 o&)v{q  
da,11,all O& 1z-  
da,15,all 717OzrF}A?  
da,19,all w&>*4=^a  
da,23,all }1mkX\wWP  
da,26,all 8xt8kf*k  
da,27,all .^wBv 'Y  
da,30,all 4jw q$G  
da,31,all = G>Y9Sc  
!施加重力 _/NPXDL  
allsel,all ^8=e8O  
acel,,9.8,0 c{3P|O&.  
!进入求解器 *pYawT  
/solu 7?ILmYBw  
antype,static 0O?\0k;o  
nusubst,100 0C4Os p  
pred,on #('GGzL6c  
nropt,full AbL(F#{  
nlgeom,on \HL66%b[  
insrch,on }p>l,HD  
outres,all,last RN2z/F Uf  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定 2<E@f0BVAy  
cnvtol,u,,0.05,2,1     !位移收敛准则设定 Fu>;hx]s  
allsel wWVB'MRXB,  
solve {^Rr:+  
save tkP& =$  
finish %x8vvcO^t  
/post1  !进入后处理 re fAgS!=q  
set,last |,T"_R_K  
plesol,s,y juA}7   
!plesol,u,y,0,1 |)OC1=As  
!plesol,s,xy,0,1 ]$!7;P  
!plesol,epto,xy,0,1 #!C|~=  
!pldisp,1                                      !绘制边坡模型变形图 6~O;t'd  
!plesol,u,x                                    !绘制边坡模型水平方向位移云图 5^N y6t  
!plesol,eppl,eqv                               !绘制边坡模型塑性应变云图 f{-,"6Y1  
离线51200605

发帖
334
土币
46
威望
1025
原创币
0
只看该作者 1楼 发表于: 2009-10-20
自己顶个,等待高人!
离线lijunli1995

发帖
273
土币
311
威望
1970
原创币
0
只看该作者 2楼 发表于: 2009-10-20
是不是密度的单位拟弄错了 %3yrX>Js  
怎么那么小
robert li
离线lijunli1995

发帖
273
土币
311
威望
1970
原创币
0
只看该作者 3楼 发表于: 2009-10-20
你这都收敛不了呀
robert li
离线lijunli1995

发帖
273
土币
311
威望
1970
原创币
0
只看该作者 4楼 发表于: 2009-10-20
finish                                     p\}!uS4 (  
/CLEAR,start                           W~5gTiBZ]  
/fILEname,try                   Of#"nu  
/units,si                                 lNMJcl3  
/pnum,kp,1                                 tm.&k6%  
/pnum,line,1                             2RdpVNx\y  
/number,0                               cR/e Zfl  
/PREP7                                     wj8\eK)]L  
*afun,deg                                 Gh}* <X;N  
!定义单元类型                             1ZXRH;J40  
et,1,solid45                               g&\;62lV%  
!定义材料属性                               ~FVbL-2  
!边坡围岩材料1的属性                       SduUXHk  
mp,ex,1,6e7                               L+G i  
mp,nuxy,1,0.29                             f\;f&GI  
mp,dens,1,23                               r-Oz k$  
tb,dp,1                                   OI:T#uk5  
tbdata,1,3e4,30,0                                                             w+{{4<+cd  
!边坡围岩材料2的属性                     On}b|ev  
mp,ex,2,5e7                               9hgIQl  
mp,nuxy,2,0.21                             .uB[zJc  
mp,dens,2,23.5                             dKmPKeJM  
tb,dp,2                                   C't%e  
tbdata,1,5e4,0,0                                                               Lr Kx  
!边坡围岩材料3的属性                     ?lm<)y?I7+  
mp,ex,3,7e8                               y!x-R !3  
mp,nuxy,3,0.34                               orFB*{/Z  
mp,dens,3,25.1                             ]d*O>Pm  
tb,dp,3                           Z ZT2c0AK  
tbdata,1,5e5,40,0                                                               ;\]DZV4?)r  
!边坡围岩材料4的属性                       MfraTUxIo/  
mp,ex,4,1.5e9                             }Qh%Z)  
mp,nuxy,4,0.3                             O O-Obg^  
mp,dens,4,24.8                             *-'`Ea  
tb,dp,4                                   G+WCE*  
tbdata,1,8e4,40,0                                                           iZ0.rcQj'o  
!创建边坡模型线 5                           N6v*X+4JH  
k,1,101.95,69.87                           lUiO|  
k,2,103.17,55.92             y2PxC. -  
k,3,103.77,49.01                         `FK qVd  
k,4,104.57,39.96                         xR;z!Tg)  
k,5,105,35                               eGUe#(I /  
k,6,140,35                               )>]SJQ!k  
k,7,140,0                               .UU)   
k,8,105,0                                   @h5Q?I  
k,9,0,0                                   '.e 5Ku  
k,10,0,35                                   w$`u_P|@E:  
k,11,0,49.11                               r<;Y4<,BZ  
k,12,0,63.68                               I.o3Old  
k,13,0,66.77                               F#o{/u?T  
k,14,0,69.87                               Ok phbAX  
l,1,2                                   _k5$.f:Yj<  
l,2,3                                   7'5/T]Z  
l,3,4                                   {"0n^!  
l,4,5                                   JEfhr  
l,5,6                                   !v*#E{r"g=  
l,6,7                                   _+gpdQq\p  
l,7,8                                   ZJQkZ_9@2  
l,8,9                                   jRp @-S#V  
l,9,10                                     crJNTEz  
l,10,11                                 v%QC p  
l,11,12                                 PAiVUGp5[  
l,12,13                                 <#~n+,  
l,13,14                                  LNvkC4  
l,14,1                                   xzRC %  
l,13,2                                   :~p_(rE  
l,12,3                                   ]@qD4:  
l,11,4                                   \\/ !I   
l,10,5                       [n +(  
l,5,8                                     =|d5V%mK  
!创建边坡面模型                         ?[ D6|gp  
al,14,1,15,13                     nb@<UbabW}  
al,15,2,16,12                     R=W$3Ue~,  
al,16,3,17,11                     ZRUAw,T*  
al,17,4,18,10                     `Mo%)I<`=  
al,18,19,8,9                   .yb8<qs  
al,5,6,7,19                               z u1gP/  
vext,1,6,1,,,35,,,                           &A^2hPe}  
vsel,s,volu,,1                             !9^GkFR6n  
vatt,1,1,1                                 7>gW2 m  
vsel,s,volu,,2                           /sVmQqVY  
vatt,2,2,1                                 Koz0Xy  
vsel,s,volu,,3                             K,*IfHi6[  
vatt,3,3,1                                 ktv{-WG2_  
vsel,s,volu,,4,6,1                         SNV;s,  
vatt,4,4,1                           O>ZJOKe  
allsel,all                                     M<@9di7c  
!划分网格                                    XV !UeBq  
mshkey,1                                 r?x~`C  
mshape,1,3D                               HPK}Z|Vl  
vmesh,all                                     Q8DKU  
da,11,all                                 +PO& z!F  
da,15,all                                 aX~' gq>  
da,19,all                                 eh*F/Gu  
da,23,all                                 efh1-3f  
da,26,all                                 ^fM=|.?  
da,27,all                                 2Mu(GUe;  
da,30,all                                 5 d|+c<  
da,31,all                                 eoPoG C  
!施加重力                                 V.[b${  
allsel,all                               ?#__#  
acel,,9.8,0                               _K~?{".  
!进入求解器                               ~'):1}KN]  
/solu                                     b{ W ,wn  
antype,static 7RgnL<t~:8  
nlgeom,on                             7.C]ZcU  
nsubst,500,500,200                                                                                                                                                         P2)g%$ME  
cnvtol,f,,0.005,2,0.5  !力收敛准则设定     &0mhO+g   
cnvtol,u,,0.05,2,1     !位移收敛准则设定 &jDRRT3  
allsel  
robert li
离线lijunli1995

发帖
273
土币
311
威望
1970
原创币
0
只看该作者 5楼 发表于: 2009-10-20
材料3和4定义错误,应该是TBDATA,1,.....不是TBDATA,3 Q?%v b  
定义子步出错,命令流就是错的
robert li
离线lijunli1995

发帖
273
土币
311
威望
1970
原创币
0
只看该作者 6楼 发表于: 2009-10-21
给你改了半天也没回来看 H/ub=,Ej*  
更没有感谢 6:\z8fYD  
FENT
robert li
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

      工业和信息化部备案管理系统网站