workbench焊接实例_[转载]Workbench的焊接模拟过程（高斯移动热源）

问题阐述：三维平板堆焊焊接残余应力分析，先分析焊接温度场，然后转化为结构分析焊接残余应力，板材模型为0.2×0.2×0.006m，焊接速度10mm/s。考虑其对称性，建立了板材的一半模型进行分析。

模型如图所示。

一半模型

靠近焊缝的部位网格要细化。

网格划分

对加热的表面named

selection为A1、A2，因为后面的高斯热源加载采用的是APDL，为后面加高斯热源作准备。

加载热源面

设置好瞬态热分析的步长和对流条件，插入Command定义高斯热源。

*DEL,_FNCNAME

*DEL,_FNCMTID

*DEL,_FNC_C1

*DEL,_FNC_C2

*DEL,_FNC_C3

*DEL,_FNCCSYS

*SET,_FNCNAME,'GAOSI'

*DIM,_FNC_C1,,1

*DIM,_FNC_C2,,1

*DIM,_FNC_C3,,1

*SET,_FNC_C1(1),2000

*SET,_FNC_C2(1),0.01

*SET,_FNC_C3(1),0.007

*SET,_FNCCSYS,0

! /INPUT,HANJIE.func,,,1

*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%

! Begin of equation: Qm*exp(-3*({X}^2

({Y}-V*{TIME})^2)/R^2)

*SET,%_FNCNAME%(0,0,1), 0.0,

-999

*SET,%_FNCNAME%(2,0,1), 0.0

*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%

*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%

*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%

*SET,%_FNCNAME%(6,0,1), 0.0

*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0,

0

*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0,

-1

*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2

*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3, 0, 0,

-3

*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1

*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0,

2

*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 17, -1

*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 1, 18, 3,

1

*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, 3, 2,

-1

*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 2, 0, 0, -4

*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -4, 17,

-1

*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -3, 1,

-5

*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, -2, 3,

-1

*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, 19

*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, 19, 17,

-1

*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -3, 4,

-2

*SET,%_FNCNAME%(0,17,1), 0.0, -1, 7, 1, -1, 0, 0

*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 1, 17, 3,

-1

*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -2, 0, 0

! End of equation: Qm*exp(-3*({X}^2

({Y}-V*{TIME})^2)/R^2)

D,A1,TEMP,%GAOSI%

D,A2,TEMP,%GAOSI%瞬态热分析设置完毕，然后就是求解，大概需要10分钟。

不同时刻时的温度场结果：

温度分布

焊接温度场的动画：

温度场动画

然后把瞬态热分析的结果转化为瞬态结构分析，把温度场作为焊接残余应力分析的热载荷。

在瞬态结构分析模块里，设置好焊接约束条件，设置好时间和步长，读取每一步温度载荷作为热载荷进行结构分析，要注意设置热载荷读取的时间和结构分析时间一致(默认条件下结构分析读取热分析的最后一步作为热载荷)。

求解时间大概半个小时，不同时刻的应力场分布：

应力场

应力场的动态过程：

应力场动画

不同时刻的焊接变形：

焊接变形

焊接变形的动态过程：

焊接变形动画