应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ansys焊接热应力分析温度场是对的，残余应力分布是错的
北京石油化工学院2026年研究生招生接收调剂公告
121/2返回列表12下一页
查看: 5341  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhangfeiniu

新虫 (初入文坛)

[求助] ansys焊接热应力分析温度场是对的，残余应力分布是错的 已有2人参与

ansys模拟温度场的分析是正确的，可以模拟出来冷却到常温，但是用间接法做应力分析时，得到的应力分布图是错的，就是焊接过程结束后，冷却阶段应力分布没有明显变化，不知道是什么原因，下面是我的部分代码
/SOL
TUNIF,20
TREF,20
ANTYPE,TRANS           !瞬态分析类型  
TRNOPT,FULL            !指定瞬态分析方法为完全法
TIMINT,ON                  !打开瞬态效应，即打开自动时间加载
AUTOTS,ON                  !使用自动时间步长跟踪  
TINTP,,,,1,,,,,,           !设置一阶瞬态积分参数为1，以便瞬态时间积分采用欧拉向后差分算法
LNSRCH,ON   
NROPT, FULL                        !使用完全的Newton-Raphson（牛顿-拉普森选项）   
*SET,LSIZE,0.001                   !节点间距
*SET,U1,21                         !电压  
*SET,I1,120                      !电流
*SET,n1,0.70                       !热效率   
*SET,Q1,(U1*I1*n1)  
*SET,pi,3.1415926   
*SET,rb1,0.003                     !高斯函数参数  
*SET,QMAX1,( Q1/(2*pi*rb1**2))  
*SET,R1,0.0045                      !加热斑点有效半径   
*SET,TINC,0.2                   !时间步长   
*SET,TI,10                         !焊接时间
*SET,v1,0.01                      !MIG焊枪焊接速度
*SET,L,0.1                         !焊接距离
*SET,NUM_TIME,TI/TINC+1                 !热流密度表各维大小
*SET,NUM_X,L/LSIZE+1
*SET,NUM_z,61   
*DIM, FLUX1, TABLE, NUM_X,NUM_z,NUM_TIME,X,z,TIME      !定义数组变量  
*DO,K,1,NUM_TIME,1                         !开始焊接循环  
*DO,I,1,NUM_X,1                            !初始化热源移动方向上的X坐标
*SET,FLUX1(I,0,K),(I-1)*LSIZE   
*ENDDO  
*DO,J,1,NUM_z,1                            !初始化热源移动方向上的z坐标
*SET,FLUX1(0,J,K),(J-1)*LSIZE   
*ENDDO  
*SET,X1,v1*(K-1)*TINC                      !设置热源中心点的坐标，使热源移动起来  
*SET,z1,0.03
*DO,I,1,NUM_X,1
*DO,J,1,NUM_z,1
!计算焊接方向上各点到热源中心点的距离
*SET,DISTANCE1,SQRT(((I-1)*LSIZE-X1)**2+((J-1)*LSIZE-z1)**2)
*IF,DISTANCE1 ,LE,R1,THEN                  !在半径为R的圆形域内加载热流密度  
*SET,FLUX1(I,J,K),QMAX1/EXP(1/2*DISTANCE1**2/(rb1**2))  
*ELSE            !其他区域的高斯热流密度为0
*SET,FLUX1(I,J,K),0
*ENDIF  
*ENDDO  
*ENDDO  
*SET,FLUX1(0,0,K),(K-1)*TINC
*ENDDO               !瞬态计算整个表格的热流密度值
TUNIF,20,   
TREF,20
SFA,4,1,CONV,-1,20  
SFA,15,1,CONV,-1,20
SFA,8,1,CONV,-1,20  
SFA,14,1,CONV,-1,20
SFA,10,1,CONV,-1,20
SFA,6,1,CONV,-1,20  
SFA,5,1,CONV,-1,20  
SFA,12,1,CONV,-1,20
SFA,16,1,CONV,-1,20
SFA,1,1,CONV,-1,20  
SFA,13,1,CONV,-1,20
SFA,9,1,CONV,-1,20  
SFA,3,1,CONV,-1,20  
/solu   
TIME,10 !设置求解时间
AUTOTS,1
KBC,1   
DELTIM,0.2,0.2,0.2,1         !使用前一载荷步的终止时间为下一载荷步的其实时间
TSRES,ERASE                  !消除上一载荷步所加的热流载荷密度  
SFA,11,1,HFLUX, %FLUX1%     !给选定的面9施加热流密度 TSRES,ERASE
OUTRES,ALL,ALL,               !将所有的求解信息写入数据库  
solve   
SFADELE,11,,HFLUX   
TIME,20
AUTOTS,1
NSUBST,20,20,20
KBC,1   
TSRES,ERASE
solve   
TIME,50
AUTOTS,1
NSUBST,30,30,30
KBC,1   
TSRES,ERASE
solve   
TIME,350
AUTOTS,1
NSUBST,30,30,30
KBC,1   
TSRES,ERASE
solve   
save
FINISH  
/PREP7  
/COM,  Structural   
ETCHG,TTS   
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,300
MPTEMP,3,430
MPTEMP,4,530
MPTEMP,5,700
MPTEMP,6,1300   
MPDATA,EX,1,,4.5E10
MPDATA,EX,1,,2.5e10
MPDATA,EX,1,,1.93e10
MPDATA,EX,1,,1.37e10
MPDATA,EX,1,,4e9
MPDATA,EX,1,,4e9
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
TB,BKIN,1,6
TBTEMP,20,1
TBDATA,1,145E6,16E9
TBTEMP,300,2
TBDATA,1,25E6,9E9   
TBTEMP,430,3
TBDATA,1,16E6,0.8E9
TBTEMP,530,4
TBDATA,1,16E6,0.8E9
TBTEMP,700,5
TBDATA,1,15E6,0.2E9
TBTEMP,1300,6   
TBDATA,1,15E6,0.2E9
TUNIF,20
TREF,20
ANTYPE,TRANS           !瞬态分析类型  
TRNOPT,FULL            !指定瞬态分析方法为完全法
KBC,1                     !阶跃性加载  
TIMINT,ON                  !打开瞬态效应，即打开自动时间加载
AUTOTS,ON   
FINISH  
/SOLU   
FLST,2,3,4,ORDE,3   
FITEM,2,9   
FITEM,2,21  
FITEM,2,25  
!*  
/GO
DL,P51X, ,UX,   
*DO,I,1,50  
LDREAD,TEMP,,,0.2*I,0,'g34','rth',' ' !读入热分析的计算结果   
OUTRES,ALL,ALL,
TIME,0.2*I  
DELTIM,0.2,0.075,0.2,1  
SOLVE   
*ENDDO  
*DO,I,1,20  
LDREAD,TEMP,,,10+I*0.5,0,'g34','rth',' '
OUTRES,ALL,ALL,
TIME,10+I*0.5   
DELTIM,0.5,0.5,1,1  
SOLVE   
*ENDDO  
*DO,I,1,80  
LDREAD,TEMP,,,20+I,0,'g34','rth',' '
OUTRES,ALL,ALL,
TIME,20+I   
DELTIM,1,1,1,1  
SOLVE   
*ENDDO  
*DO,I,1,80  
LDREAD,TEMP,,,100+I,0,'g34','rth',' '   
OUTRES,ALL,ALL,
TIME,100+I  
DELTIM,1,1,1,1  
SOLVE   
*ENDDO  
FINISH
ansys焊接热应力分析温度场是对的，残余应力分布是错的
应力分布图.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

焊接有限元分析

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-26 16:07:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有2个 )

fangdian

木虫 (正式写手)
zhangfeiniu: 回帖置顶 2013-11-27 14:23:14
检查一下你每一步读入的温度场数据有问题吗？
赞一下 回复此楼
2楼2013-11-26 22:13:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfeiniu

新虫 (初入文坛)
zhangfeiniu: 回帖置顶 2013-11-27 14:06:07
引用回帖:
2楼: Originally posted by fangdian at 2013-11-26 22:13:27
检查一下你每一步读入的温度场数据有问题吗？

温度场每一步的结果我都看了的，没有问题啊LDREAD,TEMP,,,100+I,0,'g34','rth',' '还是这句程序有问题？？
赞一下 回复此楼
3楼2013-11-27 10:36:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

htbbzzg

铁杆木虫 (著名写手)
★
1592203609: 金币+1, 谢谢回帖 2013-11-28 12:52:26
抱歉，最近太忙，没有时间试算你的问题。只能简单看一下你的命令流，提一点问题：
好像你的材料数据有问题：
  1. 你这段命令流中似乎没有定义热膨胀系数 - APLX，不知其它部分是否定义过。你的应力图没有标尺，所以无法判断是否确实缺少了热膨胀系数的定义。
  2.  高温时的弹性模量和切线模量似乎太大了，我看到的一些焊接热应力分析问题中，1300 度左右时，弹性模量和切线模量的量级比常温时低 4～5 个数量级，而你的弹性模量和切线模量，在 1300 度时只比常温降低了不到两个量级。这样，高温时变形偏小，降温后应力应变就不会大。
    主要原因可能是第一条 (不知道是不是我弄错了)。
赞一下(4人) 回复此楼
4楼2013-11-27 18:07:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)
如果你确实已经定义过热膨胀系数，那就不是这个问题；不过你说的那本书的说法肯定是错误的，如果按照那种说法，高温下材料就不可能软化甚至熔化了。高温下的材料性能最好用实测值，如果取常温情况的 10％，也需要弄清楚所说的高温是多少度，不是任何温度都能这样做的。。
赞一下(2人) 回复此楼
6楼2013-11-28 15:30:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zhangfeiniu

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by htbbzzg at 2013-11-27 18:07:03
抱歉，最近太忙，没有时间试算你的问题。只能简单看一下你的命令流，提一点问题：
好像你的材料数据有问题：
  1. 你这段命令流中似乎没有定义热膨胀系数 - APLX，不知其它部分是否定义过。你的应力图没有标尺，所 ...

您好，热膨胀系数在温度场的模拟参数部分定义过了，第二点是，我在一本书上看到说，高温时的弹性模量和切线模量最好是常温的10%左右，取太小的话结果会不准确。
赞一下 回复此楼
5楼2013-11-28 12:38:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfeiniu

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by htbbzzg at 2013-11-28 15:30:36
如果你确实已经定义过热膨胀系数，那就不是这个问题；不过你说的那本书的说法肯定是错误的，如果按照那种说法，高温下材料就不可能软化甚至熔化了。高温下的材料性能最好用实测值，如果取常温情况的 10％，也需要弄 ...

哦哦，我自己记得乱七八糟，找到书再仔细看了一下，说的是：对于力学的熔点温度以上的屈服应力和加工硬化系数，为了计算的稳定性和精度，通常将它们设定为常温值的1%~10%左右，如果设置值过小，则计算时间和计算误差就会增大。
赞一下 回复此楼
7楼2013-11-28 18:59:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvss2007

新虫 (初入文坛)

【答案】应助回帖

你温度场计算冷却温度的时间步和你后续计算应力时候读取的温度结果时间步不一致，理论是可以对温度场结果通过线性差值得到，但是这样做就会有误差了，我觉得如果你想在应力场用那样的时间步计算的话，在温度场最后也用这样的时间步计算
赞一下 回复此楼
8楼2014-04-14 23:40:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suisuinian

新虫 (初入文坛)
你好，打扰了。我最近也准备做类似问题，想找些书籍先了解一下。想问问你都参考过什么书？谢谢。
赞一下 回复此楼
9楼2015-11-28 22:59:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李朋飞619

铜虫 (小有名气)
不知你解决了没有？我也是这个问题
赞一下 回复此楼
10楼2016-03-26 01:50:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangfeiniu 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料334求调剂 +18 Eecho# 2026-04-03 18/900 2026-04-06 00:56 by fmesaito
[考研] 化学调剂 +16 艾志恒 2026-04-03 17/850 2026-04-06 00:27 by 永字号
[考研] 285求调剂 +4 恶法大二的气味� 2026-04-05 5/250 2026-04-05 20:32 by 286640313
[考研] 一志愿生物与医药，296分，求调剂 +9 66鹿 2026-04-03 10/500 2026-04-05 20:11 by lys0704
[考研] 070300化学求调剂 +17 小黄鸭宝 2026-03-30 17/850 2026-04-05 12:03 by 宁馨哈哈
[考研] 调剂求助 +10 想换手机不想解� 2026-04-02 13/650 2026-04-05 09:41 by sam3303
[考研] 324求调剂 +9 想上学求调 2026-04-03 9/450 2026-04-04 23:57 by 果冻大王
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 材料科学与工程339求调剂 +12 hyz0119 2026-03-31 13/650 2026-04-03 18:33 by ls刘帅
[考研] 生物学硕341求调剂 +4 你笑起来像云朵 2026-04-03 4/200 2026-04-03 10:32 by macy2011
[考博] 材料工程专业硕士申博 +3 麟正宇 2026-03-30 3/150 2026-04-02 15:04 by greychen00
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考研] 安全工程 285 求调剂 +3 Xinyu56 2026-04-01 4/200 2026-04-01 21:50 by 静静静静静静静�
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 349求调剂 +6 吃的不少 2026-04-01 6/300 2026-04-01 17:55 by JYD2011
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 考研调剂求助 +7 13287130938 2026-03-31 7/350 2026-03-31 16:39 by 690616278
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
信息提示
关闭
请填处理意见
关闭