24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1295  |  回复: 4

346507037

铁虫 (正式写手)

[交流] ansys温度场转应力场计算已有2人参与

麻烦大神看一下小弟的计算应力场命令流有什么问题吗?
计算一个载荷步的应力场,没有警告,不知道为什么没有计算出来应力值,麻烦大家看看命令流有什么问题吗?是不是求解设置的原因?
/com,structural
/prep7
alls
eplot
lsclear,all    !删除载荷,从数据库中清除实体模型载荷
etchg,tts       !将热单元转换为结构单元

MPTEMP,1,25,100,200,300,400,500,
MPTEMP,7,600,700,800,900,1000,1100, !定义材料属性温度范围
MPDATA,EX,1,1,1.96E11,1.92E11,1.84E11,1.76E11,1.68E11,1.60E11, !弹性模量(Pa)
MPDATA,EX,1,7,1.60E11,1.52E11,1.44E11,1.35E11,1.27E11,1.19E11,

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,15.4E-5,15.8E-5,16.3E-5,17E-5,17.9E-5,19.1E-5,   !热膨胀系数
MPDATA,ALPX,1,7,19.1E-5,15.4E-5,15.8E-5,16.3E-5,17E-5,17.9E-5,


TB,BISO,1,6,2,  
TBTEMP,25   
TBDATA,,3.2E8,1.96E9,,,,
TBTEMP,200  
TBDATA,,1.67E8,1.84E9,,,,
TBTEMP,600  
TBDATA,,1.23E8,1.52E9,,,,
TBTEMP,900
TBDATA,,1.05E8,1.27E9,,,,
TBTEMP,1200
TBDATA,,0.24E8,0.2E9,,,,
TBTEMP,1300
TBDATA,,0.105E8,0.07E9,,,,  


MPTEMP,1,25,100,200,300,400,500,
MPTEMP,7,600,700,800,900,1000,1100, !定义材料属性温度范围
MPDATA,EX,2,1,2.23E11,2.19E11,2.12E11,2.05E11,1.98E11,1.92E11,  !弹性模量(Pa)
MPDATA,EX,2,7,1.84E11,1.76E11,1.67E11,1.53E11,1.43E11,1.43E11,

MPDATA,NUXY,2,1,0.30,0.31,0.32,0.33,0.34,0.35,  !泊松比
MPDATA,NUXY,2,7,0.36,0.37,0.38,0.39,0.40,0.40,  
MPDATA,ALPX,2,1,1.02E-5,1.12E-5,1.21E-5,1.26E-5,1.28E-5,1.36E-5, !热膨胀系数
MPDATA,ALPX,2,7,1.38E-5,1.46E-5,1.54E-5,1.70E-5,1.81E-5,1.81E-5  

TB,BISO,2,6,2,  
TBTEMP,25   
TBDATA,,6.08E8,2.23E9,,,,
TBTEMP,200  
TBDATA,,5.8E8,2.12E9,,,,
TBTEMP,400
TBDATA,,5.12E8,1.98E9,,,,
TBTEMP,600  
TBDATA,,4.7E8,1.84E9,,,,
TBTEMP,800
TBDATA,,2.68E8,1.67E9,,,,   
TBTEMP,1000
TBDATA,,0.78E8,1.43E9,,,,







!*********************施加位移约束***************

csys,11
VSEL,S,MAT,,1
ESLV,S
NSLE,s,1
d,all,ux,0,,,,uy,uz                    
alls
eplot




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

!*****************杀死单元**************
VSEL,S,MAT,,2                           !*选择增材制造部分体积
ESLV,S                                  !*选择增材制造部分单元
CM,weld_elem,ELEM
EKIll,ALL
ESEL,s,elem,,weld_elem
NSLE,s                                  !*选择增材制造部分节点
CM,weld_node,node
ALLSEL,ALL
ESEL,S,LIVE
EPLOT


!************建立局部坐标系**********
LOCAL,11,0,X0,Y0,Z0,,,,1,1,
csys,11
NSEL,s,node,,weld_node
NSEL,r,loc,x,-a2,a1
NSEL,r,loc,y,-b,b
NSEL,r,loc,z,-0.04E-3,0        !*选择一个载荷步内的节点
CM,cur_node,node
ESLN,s
Ealive,all                   !*激活一个载荷步内的单元
NSEL,s,node,,weld_node
NSEL,u,node,,cur_node
CM,weld_node,node
NSEL,s,node,,cur_node
NSEL,r,loc,x,-a2,a1   !*选择前半热源内的节点
LDREAD,TEMP,,,time,,,RTH

/pbf,temp,,1                 !显示材料的温度
/replot  
TM=tm+tinc
time,tm
KBC,1
AUTOTS,ON
nsubst,10
ALLSEL,ALL

solve
FINISH  
SAVE
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
1  在你的奇特帖子里,已经告诉你: 如下的 NSEL 命令有问题:
      NSEL,s,node,,weld_node       !  其中两个连续的逗号需要去掉一个,否则可能选择不到任何节点。
    不知为什么你不远修改?
  2  你读取温度结果的命令:
      LDREAD,TEMP,,,time,,,RTH
    其中参数 time 没有赋过值,默认为 0,因此可能没有读取任何温度结果,或者即使读取了温度结果,也是错误的。你需要检查一下,是否正确读取了温度结果。
2楼2019-09-20 08:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

346507037

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by shgao20 at 2019-09-20 08:38:45
1  在你的奇特帖子里,已经告诉你: 如下的 NSEL 命令有问题:
      NSEL,s,node,,weld_node       !  其中两个连续的逗号需要去掉一个,否则可能选择不到任何节点。
    不知为什么你不远修改?
  2  你读取温 ...

首先感谢您的回答,这里可能忘记改了。疏忽了。time在计算温度场的时候已经赋值,在计算应力场的时候还的重新赋值吗?
3楼2019-09-20 12:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

激光增材

新虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
力的收敛值前边是感叹号

发自小木虫Android客户端
4楼2019-09-20 13:41:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
热分析时定义力的收敛准则是无效的,在应力分析时需要重新定义。你用惊叹号将它注销了,就不会再起作用了。
  你的 LDREAD 命令中的参数 time 与热分析时定义的计算终止时间 time 无关,在这里只是作为一个参数而已。由于你没有给它赋值过,默认为零,造成 LDREAD 命令读取不到温度结果,或者即使读取了温度结果也不是你需要的。
  在应力分析部分,如果你只计算一次 (某个指定时间的应力),是否定义 TIME 无关紧要;如果你要计算多个时间的应力,则需要使用 TIME 命令来定义每次的计算终止时间。不过需要将 LDREAD 命令放到每次计算之前,目前这个位置也是不对的 - 只读取一次指定时间的温度,不会随着 TIME 的修改而改变。
5楼2019-09-21 09:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 346507037 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见