为什么这个简单的模型ansys计算都不收敛？

xhemers

这是一个平面问题，一个竖直与水平力，材料用的是dp模型，怎么计算都不收敛，请问是怎么回事？

lscookie

是不是只施加了力边界没有给位移相关的边界条件

xhemers

lscookie 发表于 2013-4-17 16:27
是不是只施加了力边界没有给位移相关的边界条件

这就是实际巷道模型，我把一定范围内的岩体等效为竖直力与水平力了，这个还有什么位移边界，我没有加重力加速度！

lscookie

不是说你加的不对，是ansys求解器中，只含有力边界条件的情况下求解不容易收敛，有节点位移边界条件才好收敛

xhemers

哦，理解了，谢谢指导！

xuehunzupo

力是不是大了，导致隧道破坏，试试变小力

xhemers

谢谢，解答，也不是这个问题！我把力已经改小了几个数量级，同样地还是不收敛！

flypanzxl

是不是没有约束所有点的Z向位移为零？

xhemers

谢谢解答！也不是这个问题，我是平面问题，当然没有z方向的位移！我贴出命令流，请高手指点下

COM
fini
/cle
*set,d,4
*set,a,3
*set,b,4
*set,s1,0.4
*set,x,20
*set,y,20
*set,H,2.9e7!垂直应力
*set,V,8.7e7!水平应力
*set,k,2.9e3!密度
*set,e,50e9!弹性模量
*set,u,0.3!泊松比
*set,c,0.1e6!粘结力
*set,j,30!内摩擦角
/prep7
ET,1,PLANE182   
MPTEMP,1,0  
MPDATA,DENS,1,,k!密度
!dp模型
MPDATA,EX,1,,e
MPDATA,PRXY,1,,u  
!TB,DP,1,,,  
!TBMODIF,1,1,c
!TBMODIF,1,2,j  
!TBMODIF,1,3,

!建模
WPCSYS,-1,0
RECTNG,,x,,y,
WPCSYS,-1,0
CYL4,x/2,y/2,0,0,d/2,180  


wpoff,x/2,y/2
RECTNG,-b/2,b/2,0,-a,



ASEL,S, , ,       2
ASEL,A, , ,       3
CM,m,AREA      
ALLSEL,ALL  
ASBA,       1,m

ALLSEL,ALL

!切分
WPCSYS,-1,0
wpoff,x/2,y/2
wpoff,-d/2-1
wpro,,,90
ASBW,       all

WPCSYS,-1,0
wpoff,x/2,y/2
wpro,,,90
ASBW,       all

WPCSYS,-1,0
wpoff,x/2,y/2
wpoff,d/2+1
wpro,,,90
ASBW,       all

WPCSYS,-1,0
wpoff,x/2,y/2
wpro,,90
ASBW,       all

WPCSYS,-1,0
wpoff,x/2,y/2
wpoff,,d/2+1
wpro,,90
ASBW,       all
WPCSYS,-1,0
wpoff,x/2,y/2-a
wpro,,90
ASBW,       all
WPCSYS,-1,0
wpoff,x/2-b/2,y/2-a
wpro,,,90
ASBW,       20

WPCSYS,-1,0
wpoff,x/2+b/2,y/2-a
wpro,,,90
ASBW,       19
NUMCMP,KP   
NUMCMP,LINE
NUMCMP,AREA


!网格划分
TYPE,   1   
MAT,       1
ESIZE,s1,0,
AMESH,15!巷道拱顶部分
AMESH,16!

AMESH,4
AMESH,5
AMESH,1
AMESH,2
AMESH,13
AMESH,14


ESIZE,s1,0,
AMESH,9!巷道直墙部分
AMESH,10!

ESIZE,s1,0,
AMESH,7!巷道直墙部分地板
AMESH,12!
AMESH,3
AMESH,11
AMESH,6
AMESH,8
AMESH,17
AMESH,18
WPCSYS,-1,0
!约束边界条件

NSEL,S,LOC,X,0  
NSEL,A,LOC,X,x   
CM,V,NODE  
ALLSEL,ALL  

NSEL,S,LOC,Y,0
NSEL,A,LOC,Y,Y   
CM,H,NODE  
ALLSEL,ALL  



!施加荷载
SF,H,PRES,H
SF,V,PRES,V

NROPT,full,,
/SOL
ANTYPE,0
NLGEOM,1

flypanzxl

不知道为什么，复制你的命令流的时候总是有乱码，因此就直接看了下你设置的边界条件
我感觉是不是应该这样设置边界条件呢，如下：
首先，面沿纵中切开，切开线上的节点做对称约束(X向位移约束)
其次，下表面固定约束，上表面照样施加载荷(固定处支反力一样，但这样约束了Y向位移)
我感觉这样应该就可以了
希望对你有所帮助

xhemers

谢谢，你的解答，我重新发个txt命令流！根据你的建议，我增加了底部约束，力是收敛了，但是位移不收敛？

flypanzxl

从你的施加边界条件开始，替换成以下程序：
alls
nsel,s,loc,x,0.5*x
dsym,symm,x
nsel,s,loc,y,0
d,all,all
nsel,s,loc,y,y
SF,all,PRES,H
nsel,s,loc,x,0
nsel,a,loc,x,x
SF,all,PRES,V
alls
/solu
nlgeom,on
nsubst,20
outres,all,all
solve
我解下来是收敛的啊，而且我分别用直接解和非线性解两种方式，发现几乎无区别，毕竟这里没什么由于弯曲变形引起的迭加计算
希望对你有所帮助