[求助] 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

回复此楼

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

焊接有限元分析

» 猜你喜欢

• 不自信的我 已经有5人回复
• 磺酰氟产物，毕不了业了！ 已经有4人回复
• 论文终于录用啦！满足毕业条件了 已经有16人回复
• 求个博导看看 已经有19人回复
• 投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球 已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

• 关于残余热应力中剪切应力和切向应力的区别 已经有6人回复
• 用ANSYS进行热结构分析时，温度场温度超过材料参数，对模拟结构有什么影响？ 已经有4人回复
• 焊接残余应力与变形之间的联系问题 已经有3人回复
• 围护结构应力值的变化 已经有3人回复
• ansys apdl 热-应力耦合算残余应力，材料参数是哪些 已经有19人回复
• ansys瞬态热分析 已经有9人回复
• ansys热残余应力 已经有14人回复
• ansys残余应力模拟 已经有9人回复
• 请教ANSYS热应力方面的问题：所模拟出的热应力出现振荡的现象。谢谢！ 已经有13人回复
• 关于热应力计算的求助 已经有9人回复
• 关于热应力分析的耦合问题 已经有20人回复
• ANSYS活塞热应力分析 已经有5人回复
• 用ANSYS模拟焊接温度应力场有没有什么具体的实例啊？ 已经有5人回复
• 【求助】ansys热应力仿真问题 已经有5人回复
• 【求助】ANSYS求解水果在压力下的应力和应变分析 已经有5人回复
• 【求助】ANSYS分析应力过程中存在的问题 已经有3人回复