24小时热门版块排行榜    

查看: 1448  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 纯娱乐,不喜欢勿喷 +4 Tide man 2026-08-04 7/350 2026-08-04 11:23 by Tide man
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
信息提示
请填处理意见