24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2531  |  回复: 3

346507037

铁虫 (正式写手)

[交流] ansys应力计算不收敛已有2人参与

大神ANSYS只读一个载荷步温度的计算应力场就出现了不收敛的现象,通过这个警告能和计算应力收敛图判断出应力不收敛的 原因吗?
工艺是激光增材制造
以下是命令流
/com,structural
/prep7
alls
eplot
lsclear,all    !删除载荷,从数据库中清除实体模型载荷
etchg,tts       !将热单元转换为结构单元




MPTEMP,1,20,100,250,400,550,700,
MPTEMP,7,850,1000,1200,1400,1600,1750, !不锈钢定义材料属性温度范围
MPDATA,EX,1,1,2.08E11,2.05E11,2.02E11,1.94E11,1.86E11,1.79E11, !弹性模量(Pa)
MPDATA,EX,1,7,1.70E11,1.62E11,1.27E11,0.17E11,0.15E11,0.10E11,

MPDATA,NUXY,1,1,0.30,0.30,0.31,0.31,0.31,0.31,  !泊松比
MPDATA,NUXY,1,7,0.31,0.32,0.32,0.32,0.33,0.33,

MPDATA,ALPX,1,1,12.8E-6,12.8E-6,13.5E-6,13.9E-6,14.2E-6,14.4E-6,   !热膨胀系数
MPDATA,ALPX,1,7,15.1E-6,16.1E-6,16.1E-6,16.2E-6,16.5E-6,17.0E-6,


TB,BISO,1,6,2,  
TBTEMP,25   
TBDATA,,3.2E8,3.2E5,,,,
TBTEMP,300  
TBDATA,,3.08E8,3.08E5,,,,
TBTEMP,600  
TBDATA,,2.46E8,2.46E5,,,,
TBTEMP,800
TBDATA,,2.26E8,2.26E5,,,,
TBTEMP,1000
TBDATA,,2.07E8,2.07E5,,,,
TBTEMP,1750
TBDATA,,1.14E8,1.14E5,,,,  


NSEL,S,LOC,Z,0,jg
NSEL,R,LOC,Z,0,
d,all,ux,0,,,,uy                  
alls
eplot



/solu
antype,4
trnopt,full            !完全瞬态积分法
lumpm,0                !使用与单元相关的质量矩阵公式
nlgeom,off             !打开大变形选项
pred,on                !打开预测校正
sstif,1                !打开应力刚化功能
solcontrol,on          !自动控制
nropt,full,,on         !完全牛顿-拉普森方法
cnvtol,u,,0.5,l2,-1    !设置收敛值,位移
cnvtol,f,,0.05,l2,-1   !设置收敛值,力
cutcontrol,plslimit,10,!设置在一个非线性求解中时间步的减少
eqslv,sparse,,2        !求解器类型
lnsrch,on              !线性搜索
pstres,1               !打开预应力选项   
autots,on              !自动步长        
tintp,,,,1,,           !向后积分        
timint,on              !瞬态效果  
tref,25                !参考温度
kbc,0                  !连续性载荷
neqit,50               !规定每个子步中最大迭代次数,默认为25
timint,on


LOCAL,11,0,X0,Y0,Z0,,,,1,1,
csys,11
NSEL,s,node,,weld_node
NSEL,r,loc,z,-CH,0
ESLN,s,1
Ealive,all
ESEL,S,LIVE
EPLOT        !*激活整层的单元

TM=tm+tinc
time,tm
KBC,1
AUTOTS,ON
nsubst,10
LDREAD,TEMP,,,time,,,RTH
/pbf,temp,,1                 !显示材料的温度
/replot  
ALLSEL,ALL
solve
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
有几个问题可以考虑一下:
  1  打开了大变形选项,就不要打开应力钢化和预应力选项;
  2  你的材料性能中,切线模量有点小,几乎是理想塑性材料了,难以收敛,可以将 E5 改为 E8;
  3  检查一下结构的几何模型,看是否存在内尖角,如果有内尖角,很容易产生应力集中,导致不收敛;
  4  你在应力分析时读取了时刻 time 的温度,但是 参数 time 没有赋过值。你可以检查一下读入的温度值是什么情况;
  5  检查一下温度结果是否有过大的温度和温度梯度,也会造成应力分析是不收敛。

» 本帖已获得的红花(最新10朵)

2楼2019-10-13 10:03:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

346507037

铁虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by shgao20 at 2019-10-13 10:03:32
有几个问题可以考虑一下:
  1  打开了大变形选项,就不要打开应力钢化和预应力选项;
  2  你的材料性能中,切线模量有点小,几乎是理想塑性材料了,难以收敛,可以将 E5 改为 E8;
  3  检查一下结构的几何模型 ...

好的,大神。我再试试
3楼2019-10-13 10:07:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyang617730

木虫 (著名写手)

微信公众号:CAE初行者,ID:


小木虫: 金币+0.5, 给个红包,谢谢回帖
我一般在wb里面做,不收敛通常有以下做法:
1,有塑性的,大变形打开
2,不收敛时可以加弱弹簧(weak spring),会有所改善
3,缩短载荷步,还可以设置KBC,0,1两个状态,选ramping up的试试
4,网格,单元长宽比很重要,通常2.5以下都还好
发自小木虫Android客户端
微信公众号:CAE初行者,ID:cae_cxz
4楼2019-10-13 14:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 346507037 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见