| 查看: 4499 | 回复: 11 | |||
| 【悬赏金币】回答本帖问题，作者zhangfeiniu将赠送您 5 个金币 | |||
[求助]
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  应力分布图.jpg |
回复此楼» 收录本帖的淘贴专辑推荐
 焊接有限元分析 |
» 猜你喜欢
想问一下有没有需要文章但是没时间带学生的研究生导师。
已经有22人回复
-
院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度
已经有11人回复
-
兄弟们帮我看看这两个结构怎么合成
已经有6人回复
-
前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函
已经有19人回复
-
東京大学 应用化学专业柳田研究室 招收修士/博士生
已经有3人回复
-
无under review，直接DIP，咋么办。
已经有10人回复
-
估计今年青基又没戏
已经有10人回复
-
特别资助审核状态
已经有8人回复
-
小孩7岁，上一年级，不自信，学习较差
已经有37人回复
-
刚博士毕业，入职普通一本高校
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于残余热应力中剪切应力和切向应力的区别
已经有6人回复
- 用ANSYS进行热结构分析时，温度场温度超过材料参数，对模拟结构有什么影响？
已经有4人回复
- ansys热应力间接耦合分析讨论
已经有12人回复
- 焊接残余应力与变形之间的联系问题
已经有3人回复
- 围护结构应力值的变化
已经有3人回复
- ansys热-结构耦合场分析，温度场是合理的，但是应力场出错了。
已经有11人回复
- ansys apdl 热-应力耦合算残余应力，材料参数是哪些
已经有19人回复
- 铜块真空降温过程中残余应力ANSYS模拟
已经有11人回复
- 如何利用abaqus模拟循环温度场热应力问题
已经有25人回复
- ANSYS用间接法进行热应力分析时，结构分析的问题
已经有12人回复
- ansys瞬态热分析
已经有9人回复
- ansys热残余应力
已经有14人回复
- ansys残余应力模拟
已经有9人回复
- 请教ANSYS热应力方面的问题：所模拟出的热应力出现振荡的现象。谢谢！
已经有13人回复
- 关于热应力计算的求助
已经有9人回复
- 关于热应力分析的耦合问题
已经有20人回复
- ANSYS算完如何提取残余应力
已经有8人回复
- ANSYS活塞热应力分析
已经有5人回复
- 用ANSYS模拟焊接温度应力场有没有什么具体的实例啊？
已经有5人回复
- 【求助】ansys热应力仿真问题
已经有5人回复
- 【求助】ANSYS求解水果在压力下的应力和应变分析
已经有5人回复
- 【求助】ANSYS分析应力过程中存在的问题
已经有3人回复
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 63
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
2楼2013-11-26 22:13:27
温度场每一步的结果我都看了的，没有问题啊LDREAD,TEMP,,,100+I,0,'g34','rth',' '还是这句程序有问题？？3楼2013-11-27 10:36:24
★
1592203609: 金币+1, 谢谢回帖 2013-11-28 12:52:26
1592203609: 金币+1, 谢谢回帖 2013-11-28 12:52:26
|
抱歉，最近太忙，没有时间试算你的问题。只能简单看一下你的命令流，提一点问题： 好像你的材料数据有问题： 1. 你这段命令流中似乎没有定义热膨胀系数 - APLX，不知其它部分是否定义过。你的应力图没有标尺，所以无法判断是否确实缺少了热膨胀系数的定义。 2. 高温时的弹性模量和切线模量似乎太大了，我看到的一些焊接热应力分析问题中，1300 度左右时，弹性模量和切线模量的量级比常温时低 4～5 个数量级，而你的弹性模量和切线模量，在 1300 度时只比常温降低了不到两个量级。这样，高温时变形偏小，降温后应力应变就不会大。 主要原因可能是第一条 (不知道是不是我弄错了)。 |
4楼2013-11-27 18:07:03
6楼2013-11-28 15:30:36
5楼2013-11-28 12:38:07
7楼2013-11-28 18:59:09
8楼2014-04-14 23:40:45
9楼2015-11-28 22:59:48
10楼2016-03-26 01:50:15
90