焊接数值模拟.ppt

2-5不稳定导热的有限元解法二维热传导 3、区域离散化 将一个矩形区域，划分成多个直角三角形。 设直角边长为h，(?x =?y=h) 节点x=rh,y=sh (r, s为正整数) 此节点记为(r,s)，（相当于（x,y）点） 2-5不稳定导热的有限元解法/二维热传导 4、温度差值函数的建立 对于三角形单元 T=f(Ti,Tj,Tm) T=NiTi+NjTj+NmTm 5、单元变分的计算 将求解区域分成有限个单元后，泛函I(T)变成各个单元内泛函的积分。 2-5不稳定导热的有限元解法/二维热传导 5、单元变分 2-5不稳定导热的有限元解法二维热传导 （5、单元变分） 2-5不稳定导热的有限元解法/二维热传导 5、单元变分 2-5不稳定导热的有限元解法/二维热传导 5、单元变分 2-5不稳定导热的有限元解法/二维热传导 5、单元变分 2-5不稳定导热的有限元解法/二维热传导 5、单元变分 2-5不稳定导热的有限元解法/二维热传导 在时间上采用向前差分： 6、总体合成 i= 1, 2, 3, ……n 上式包含若干线性方程组。对于每一个方程来说，都是对所有 单元求和而成。现以i(r,s)为例，进行求解。 2-5不稳定导热的有限元解法/二维热传导 i(r,s)点涉及六个单元Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ，所以 其它单元中不含有节点i(r,s)，它们的泛函对Ti变