三维热传导模型.docVIP

热传导模型的控制方程 在移动直角坐标系内，热传导方程为 (1) 为了处理问题的方便，改用移动柱坐标系（r,θ,z）。对方程（1）作坐标变换x=rcosθ，y=rsinθ，z=z可变为 (2) 整理后，可得 (3) 式中 α=k/(ρcp) 热传导模型的边界条件 (1) 离光斑无穷远处，工件的温度维持室温Ta。 该边界条件的数学表述如下 当r→∞（0≤θ≤2π）时，T=Ta。 (2) 设工件表面光斑大小为rb，则工件表面被激光直接辐照的区域内（即工件表面的光斑内）的点的温度可以通过下式得出： (4) 上式中，A为工件表面对激光的吸收系数，I为加工中用的激光光束的功率分布函数，设P为入射到工件表面的激光功率,对高斯光束而言, 距离光斑中心r的点的激光功率密度为 (5) 把(5)式代入(4),有 (6) 热传导模型的有限差分方程的建立 (1) 计算区域的确定及网格划分 由于有对称性（关于x轴），我们只计算x轴上部（即y≥0）区域的温度分布。 网格划分及计算区域如图1所示。 (2) 有限差分方程的建立 采用有限差分法对模型进行数值求解。为此，必须首先把控制方程化为有限差分方程。对图1所示的网格，假设r、θ和z方向的网格步长分别为Δr、Δθ和Δz，那么采用中心差分格式时，有 (7) (8) (9) (10) (1