24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1156  |  回复: 2
【悬赏金币】回答本帖问题,作者Destiny625将赠送您 50 个金币

Destiny625

铁虫 (正式写手)

[求助] 热应力分析的结果比正常值大了一个数量级已有1人参与

本人的课题方向是激光熔覆方向的,最近在做激光熔覆的温度场和应力场的仿真,温度场做出来的结果正常,然后用间接热结构耦合的方法,先转换模型的单元类型,再将温度场的计算结果作为载荷进行应力场的仿真,仿真时施加的位移约束是自由位移约束,但是应力场仿真的结果算出来的应力大小比正常值大了一个数量级,我看其他论文里人家算出来的正常值应该是10的9次方,也就是几百兆帕,但是我算出来的是10的10次方,就是几千兆帕了,这很明显太大了,而且实际情况中基板熔覆后两边应该是往上翘的,可是我的模拟结果却是往下翘的,请问各位大佬这是什么情况,怎么解决,下面附上我应力场部分的命令流。仿真结果见下图。
FINISH
/COM,Structural
/PREP7
toffst,273
ETCHG,TTS
  
MPTEMP,,,,,,,,
MPTEMP,1,0,100,200,300,400,600
MPDATA,EX,1,1,1.985E11,1.93E11,1.85E11,1.76E11,1.67E11,1.59E11
MPDATA,ALPX,1,1,17E-6,17.4E-6,18E-6,18.6E-6,19.1E-6,19.6E-6
MPDATA,PRXY,1,1,0.294,0.295,0.301,0.310,0.318,0.326
MPTEMP,7,800,1200,1300,1500,2000
MPDATA,EX,1,7,1.51E11,0.6E11,0.2E11,0.1E11,0.1E11
MPDATA,ALPX,1,7,20.2E-6,20.7E-6,21.1E-6,21.6E-6,21.6E-6
MPDATA,PRXY,1,7,0.333,0.339,0.342,0.388,0.388

MPTEMP,1,20,100,200,300,400,500
MPDATA,EX,2,1,1.951E11,1.912E11,1.857E11,1.796E11,1.726E11,1.645E11
MPDATA,ALPX,2,1,15.24E-6,15.8E-6,16.43E-6,16.97E-6,17.44E-6,17.85E-6
MPDATA,PRXY,2,1,0.267,0.273,0.29,0.31,0.322,0.313
MPTEMP,7,600,700,800,900,1000,1100
MPDATA,EX,2,7,1.55E11,1.441E11,1.314E11,1.168E11,1.001E11,0.811E11
MPDATA,ALPX,2,7,18.21E-6,18.54E-6,18.83E-6,19.11E-6,19.38E-6,19.66E-6
MPDATA,PRXY,2,7,0.296,0.282,0.262,0.24,0.229,0.223
MPTEMP,13,1200,1420,1460,2000
MPDATA,EX,2,13,0.595E11,0.02E11,0.02E11,0.02E11
MPDATA,ALPX,2,13,19.95E-6,20.7E-6,20.7E-6,20.7E-6
MPDATA,PRXY,2,13,0.223,0.223,0.223,0.223

FINISH
/SOLU    !求解
ANTYPE,4         !分析类型:瞬态
TRNOPT,FULL                                         !求解方法:完全的N-R方法
TIMINT,ON    !打开瞬态效应
KBC,0
LUMPM, 0                                             !不使用集中质量矩阵
NLGEOM,ON                        !使用大变形效应
pstres,1               !打开预应力选项  
neqit,50
TREF,20                !定义热应力计算参考温度

NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,BaseHeight
D,ALL,UX,0
D,ALL,UY,0
D,ALL,UZ,0
ALLSEL,ALL

NSEL,S,LOC,X,baselength
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,BaseHeight
D,ALL,UY,0
D,ALL,UZ,0
ALLSEL,ALL

NSEL,S,LOC,X,0
NSEL,R,LOC,Y,BaseWidth
NSEL,R,LOC,Z,BaseHeight
D,ALL,UZ,0
ALLSEL,ALL

*DO,K,1,101,1
LDREAD,TEMP,,,0.1*K,,'3','rth'
OUTRES,ALL,ALL,
TIME,0.1*K                                            
DELTIM,,,,1         !指定时间步长大小
SOLVE
*ENDDO

热应力分析的结果比正常值大了一个数量级
微信图片_20200623184051.png
回复此楼

» 猜你喜欢

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

Destiny625

铁虫 (正式写手)

2楼2020-06-25 10:00:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

应力偏大的可能原因很多,比如;
     使用的单元有问题,比如使用了没有中间节点的单元,在局部区域生成了 4 节点的四面体单元,其计算结果是不准确的
     几何模型有问题 (如结构突变,内部尖角等),造成了局部应力集中;
     材料属性有问题;
     温度结果偏大,特别是温度梯度偏大;
     点约束 (造成应力集中);
     其它问题,如求解设置过于简单,计算子步数太少,等。
  总之,必须了解你的详细情况,有时甚至需要进行试算,才能分析和解决问题。
3楼2020-06-27 09:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Destiny625 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见