焊接温度场有限元分析.docxVIP

焊接温度场有限元分析

第二类边界条件，已知边界上的热流密度分布，即，

屯＝qs（;

x,y,z,t)

(Z,4)

。绝热边界为第二类边界条件的特殊情况，在绝热表面，任一点的比热流量沿表面法线方向的温度梯度均为零，即；

。

-叮=

8n

(2.5)

第三类边界条件，已知边界上物体与周围介质间的热交换，即：

屯＝a(T－飞） （2.6)

当表面传热系数a非常大而导热系数入很小时，即a/J..-oo,表面温度接近千周围介质的温度，此时相当于等温边界条件：当表面传热系数a很小而导热系数入非常大时，即吐人-0,通过边界表面的热流趋千零，此时相当千绝热边界条件．

式中：

．

n为边界表面外法线方向(m) ;

q为单位面积上的外部热输入热流(W·m·2) ;

a为表面散热系数(W·m2．炉）；

飞为周围 介质温度(K)．

对焊接温度场的求解问题就是在初始条件和边界条件已知的倍况下，对非线

性瞬态热传导问题的控制方程进行求衅的过程。

焊接温度场的计算方法

在焊接热影响区晶粒长大模拟中，需要对焊接温度场进行计算i 获得焊接

.HAZ各点的执循环．解析法和有限元法依靠其各自的优势经常在焊接H辽晶粒长大模拟中被使用，这里主要介绍这两种方法c

2.2.］ 焊接温度场的解析计算

上世纪30年代，由于电子计算机尚未出现，对非线性问题的求解非常困难，

D.Ros邸thal和前苏联的H.H岱卡林院+等人在计算焊接温度场时，假设焊件尺寸无限大、材料的热物理参数恒定、忽略表面敕热、熔化冶热和相变潜热，阵低了求解难度，推导出一整套计算公式，并形成系统理论，莫定了解析法计笋焊按温度场的基础。尽管觯析法计算出的温度场与实践温度场有较大差异，但在远

离热源的热影响区具有很好的精戌，而且这种方法简单、计莽侯少，便于叙机应用，在模拟焊接HAZ晶粒长大时常被使用．

在固定坐标系中，无内部热源的无限大体传热方程用下式来表示：

护T 护T 护T 1T

(2.7)

式中：

T 为温度 (.C) :

t为时间(s)：

a为导品系数Cm2/s) ,

—ax?+－句，2.＋一8z2

＝－aa－t

如图 2.1所示，设点热源以速度v沿x方向移动，在固定坐标系中设A点的

坐标为（心，Oy，石），在r时间内热源由0 移至o（千

vt,0,0)，则A 点的传热公式为

T(x。,y。z。,,t)=f。.——

(x。一vV)2＋沽＋4

t I2qdtrexp(

t I2qdt

(2.8)

式中：

cp为体积热容(J/(m·K))

q为电弧功率(W)：

＇攷t,,ISh为板厚 (m)

＇

攷

t

,

,

IS

Cpl411Ct(t一t勹I2 4a(t-t)

｀II图2.1+元限大休廿动点热涛仿从役型

｀

I

I

fig.:2.1Modelofmovingheatsourceininfinitebulkpi.ccc

如果以热源当前所在的位囡为动坐标系的原点，经坐杯变换后，式 0 .8)

可改丐成：

T(iR.t)=To十盂exp气）0（忐 节）

(2.9)

,I，

,

I

，

式中：R=x斗户王．x＝xo-V1; y=yo．产邵J:(u. K}为超越函数

l(u,K)=i。..el(p(-2－枷

在到达 稳态时 ．节 ＝o.令1-+Xl, 可得 到下面的 实用公 式：

T(R)=T0＋志e.x［p一亡(x+R)]

由式(2.11)可推导出峰值温度的计算公式：

Tm=

Tm=飞＋

可evcpr

式中：1=/ +-2-,e为自然对数 常数．

222 有限元法模拟温度场

(2.10)

(2.11)

(2.12}

实际焊接过程是一个局部快速加热至高温，随后又快速冷却的过程。随若热源的移动．整个焊件的温度随时间和空间急剧变化，材料的热物理学参数也随温度剧烈变化，同时还存在熔化和相变时的潜热现象．因此，焊接温度场计觅屈千典型的非线性瞬态热传导问题。在应用解析法计算焊接温度场时所做的各种假设与焊接的实际情况不符合，因此所得的解与实际测定有一定的偏差，尤其是在焊接熔池附近的区域面且在计算中再板多层多逆焊接温度场时，应用解析法很难求解．有限元法以它在传热计算方面特有的优势可以更好地桢拟再现焊接传热过程，儿其是在校拟多层多道焊接传热过程时，有着不可替代地位，得到了国内外

广泛的应用（如图2.