《焊接数值模拟技术课件.ppt》由会员分享，可在线阅读，更多相关《焊接数值模拟技术课件.ppt（34页珍藏版）》请在三一办公上搜索。
1、,焊接数值模拟技术,0.前言 焊接过程是复杂的物理化学冶金过程，焊接数值模拟技术也涉及到多方面知识。包括：电弧物理、传热、冶金和力学等。焊接数值模拟技术涉及到焊接传热、传质、焊接接头组织及性能预测、焊接力学等方面。,1焊接热过程数值分析1.1 焊接热过程基本特点 焊接热过程，既在进行焊接时，被焊金属在焊接热源作用下将发生被加热和局部熔化的过程。焊接热过程的主要特点：（1）不均匀的局部加热（2）加热的瞬时性（3）非稳定加热（4）多种传热方式共存,1.2 焊接熔池温度场数值模拟 熔池温度场数值模拟意义 焊接温度场的模拟数学模型就是傅立叶 热传导方程。平面（x,y）导热问题 密度；Cp比热容；k导热
2、系数；Q外部热源,凝固潜热的表达式 Q相变潜热；L单位质量液体凝固潜热（J/kg）；fs固相分数 第一类边界条件，已知边界上各点的温度值 Ts=Ts(x,y,z,t)第二类边界条件，已知边界上的热流密度分布 第三类边界条件，已知边界上物质与周围介质的热交换 n边界表面外法线方向；表面换热系数；Ta为周围环境温度,1.2 焊接熔池流体动力学模拟 建模的基本假设条件：（1）焊接热源的能量分布符合高斯分布，热源的传热系数恒定（2）焊接材料的热物理参数为常数，但在液相和固相时热物理参数不同（3）流体为层流不可压缩流体，材料密度为常量,对流换热问题数学描述（1）换热微分方程 对流换热系数（2）连续性方程
3、 单位时间流入、流出微元体质量相等。（3）动量微分方程 作用于微元体表面和内部所有外力的总和，等于微元体中流体动量的变化率。,（4）能量微分方程 由导热进入微元体的热量与由对流进入微元体的热量之和等于微元体中流体的热焓增量。Q1+Q2=H焊接熔池流体动力学分析边界条件（1）上表面 pv蒸汽压力；表面张力；kr自由表面曲率,（2）下表面 qconv=hc（TT）k热导率；u径向速度；v轴向速度；qconv热流密度；hc 界面换热系数（3）对称轴 r 表示径向,（4）侧面（5）液相分数 Tl熔化点；Ts凝固点,焊接熔池传质、传热数值模拟研究主要成果（1）熔池内液体金属流动影响焊接熔深 熔池的表面张
4、力影响液体金属流动 1）如随温度升高，表面张力增加，则焊接熔深大 2）如随温度升高，表面张力减小，则焊接熔深浅（2）焊接电流线发散，增加熔深（3）浮力对熔池内流体流动的作用较小（4）熔滴对熔池的冲击力对熔深影响较小,（5）焊接热源导致熔池表面金属蒸发对熔池表 面温度的影响 1）激光焊接熔池表面金属蒸发是影响熔池 表面温度分布的主要原因 2）一般电弧焊表面张力引起的对流是影响 温度的主要原因,2焊接冶金和焊接接头组织性能的预测（1）建立预测焊缝组织发展的模型 模拟过程涉及的问题包括：1）焊接熔池中的热化学反应 2）焊缝的凝固过程 3）焊缝中的固态相变,（2）焊接接头（力学）性能预测 1）铁素体钢
5、：评价相变行为、淬硬性、强度、韧性、脆化、裂纹敏感性 2）奥氏体不锈钢：评价显微组织、热裂纹倾 向、点状腐蚀、力学性能（3）焊接热影响区相变、组织性能预测 1）基于CCT图预测焊接热影响区组织及硬度 2）利用相变动力学模型，预测焊接热影响区组 织及硬度 3）温度、相变及热应力耦合模拟,I,3焊接应力与变形数值模拟 基于有限元技术和焊接热弹塑性理论 主要研究如下问题：（1）焊接应力的发生机制和残余应力分布形态（2）焊接裂纹及其力学性能指标（3）高精度焊接变形预测（4）焊接应力、变形对焊接接头强度影响,4其它焊接应方法数值模拟（1）电阻点焊数值模拟（2）陶瓷与金属扩散焊模拟（3）激光焊数值模拟,4
6、应用电阻点焊数值模拟（1）有限元模型 1）物理模型 电阻点焊的基本特点：具有轴对称性 稳定的导电通道 电位分布和电阻分布交互影响产热 液态熔核的温度趋向均匀,点焊物理模型描述：电位及温度场轴对称分布 表面辐射及对流散热用传热系数表示 假设熔核温度最高温度,2）基本电、热方程 电势方程 U（r，z）电位函数；E电阻率 r圆柱坐标系的径向坐标 z轴向坐标,热传导方程 T（r，z）温度函数；qV单位体积单位时间内热源生成热量 热导率，CP比热容，密度 内热计算,3）有限元数学模型 整体组集方程式 T(t)未知节点温度列向量；F(t)节点温度载荷列向量；KT 整体温度刚度矩阵；C 整体变温矩阵,（2）
7、边界条件和初始条件,（3）物性参数处理 电极:CrZrCu，工件:65Mn 1）表面传热系数 水强制冷却传热系数hw 表面传热系数ha Tw材料表面温度（K）；材料热辐射率,2）材料电阻率 铜电极材料电阻率:E=16.8+0.043(T20)20 T 1000 工件材料电阻率:热导率与电导率关系:魏德曼弗兰兹定律 L=2.45108 W/K2,3）材料热导率 电极材料:=395.20.04375(T20)20 T 100=391.70.093(T100)100 T 300=373.10.097(T300)300 T 600=344.00.07733(T600)600 T 900 工件材料:65
8、钢,4）材料密度 电极材料:8900 kg/m3 工件材料:=79070.3310(T+273)20 T 780=81180.5045(T+273)780 T 1410=84560.8437(T+273)1410 T,5）材料比热容及相变潜热 电极材料比热容：Cp=356.4+0.0988(T+273)20 T 1000 工件材料比热容：按铁的比热容计算 相变潜热处理:等效比热容法 工件材料相变潜热:H=1.607104 727 T 780 Hl=2.467105 1410 T 1485 H 相变潜热；Hl 固体液体转变潜热 相变潜热折算(为比热容)公式:,6）材料硬度 利用硬度H与抗拉强度b存在线性关系,（4）有关计算 1）接触电阻 接触电阻两种类型：收缩电阻；薄膜电阻,收缩电阻计算:E 接触材料电阻率 两种材料接触时E=（E1+E2）/2 多点接触（A=n2）：考虑硬度和强度关系,2）动态电阻 总电阻：电极、工件的内电阻 工件/工件接触电阻 工件/电极接触电阻 计算方法：首先计算出各个单元的电阻，然后利用串 并联法求出总电阻。计算单元电阻：确定单元内电阻率（单元内电阻率相 同），由单元形心温度Te确定单元电 阻率。,3）电极间电压 交流点焊，电流按正弦波处理,（5）数值模拟结果的验证 1）建立正确、合理、简化的数学模型 2）合理设定边界条件 3）准确的物性参数 4）实验对比,
