筒体内壁堆焊焊后消氢感应加热温度场的数值模拟

来源：核心期刊咨询网时间：2019-11-13 10:2312

摘要：摘要：采用有限元方法对堆芯补水箱筒体内壁堆焊焊后消氢感应加热温度场进行了数值分析。针对筒体的结构特点，设计了一种加宽的C型感应线圈。研究了感应线圈的匝数、线圈与筒体的距离、线圈电流、电流频率对筒体温度分布的影响。结果表明，感应线圈的形状决定

　　摘要：采用有限元方法对堆芯补水箱筒体内壁堆焊焊后消氢感应加热温度场进行了数值分析。针对筒体的结构特点，设计了一种加宽的C型感应线圈。研究了感应线圈的匝数、线圈与筒体的距离、线圈电流、电流频率对筒体温度分布的影响。结果表明，感应线圈的形状决定了筒体温度场分布的形态，线圈中心对应筒体的温度比四周的温度高。线圈匝数、线圈电流和电流频率越大，筒体升温速度越快，温度分布越均匀。当线圈电流为215 A、电流频率为10 kHz、线圈匝数为15匝、线圈与筒体之间的距离为25 mm时，筒体加热时的溫度分布均匀性较好，能够满足焊后消氢的工艺要求。筒体感应加热试验的结果表明，模拟温度与试验温度变化趋势相同，温度值较为吻合。

　　关键词：感应加热;筒体;消氢;数值分析;温度场

　　0 前言

　　核电站一回路系统中的堆芯补水箱筒体需要在高温、高压、腐蚀性环境下长期工作，通常采取在筒体内壁堆焊镍基或不锈钢合金来提高筒体的耐蚀性和耐热性[1-2]。堆焊作业后氢气容易在焊接区域聚集析出，形成氢致裂纹，因此需要对堆焊后的筒体进行整体消氢处理[3]。常规的消氢工艺是焊后将筒体加热到250～400 ℃，保温2 h。

　　目前，工业中对大型工件焊后消氢的加热方式有火焰加热、电阻加热和感应加热。火焰加热操作简单，但是能量利用率较低，成本高且环境污染严重;电阻加热是通过电阻丝或电阻片来加热工件，工件升温速度慢，热效率低，维护成本高[4-6];电磁感应加热具有加热效率高、控温精确、绿色环保等优点[7]。因此，采用电磁感应加热的方式对堆焊后的筒体进行整体消氢热处理。

　　国内外众多学者对感应加热过程中电磁场和温度场进行了大量的数值模拟研究。章德斌[8]用Marc有限元分析软件模拟了轴类零件感应加热温度场的分布情况，结果表明，电流控制温度的加热深度随着加热时间的延长，零件内部温度逐渐升高。A. O. Glebov 等人[9]采用有限元分析的方法，研究了铁磁板在感应加热过程中的温度场分布，提出了一种用线性微分方程计算铁磁体中涡流三维场的方法，使计算时间减少了一个数量级以上。Zhu Hongda等人[10]建立了电磁场和温度场耦合的数学模型来计算钢坯热挤压感应加热过程。将模拟结果与物理建模过程中记录的温度测量值进行比较，验证了该数学模型的有效性。

　　文中以堆芯补水箱筒体内壁堆焊的电磁感应加热消氢工艺为研究对象，基于ANSYS有限元软件建立三维瞬态电磁-热耦合模型，对筒体内壁堆焊焊后感应加热温度场进行数值模拟，分析感应加热工艺参数对筒体温度场的影响规律。

　　1 有限元模型的建立

　　1.1 实体模型

　　堆芯补水箱筒体的内径为2 580 mm，壁厚120 mm，长度为3 640 mm，材质为SA5083钢。由于筒体长度较大，故分节进行感应加热，每节长1 500 mm。筒体圆周方向安装4个相同的横截面為10 mm×10 mm的正方形感应线圈。线圈采用TP2紫铜线制作，长度、宽度各为1 000 mm，呈C型且覆盖1/4筒体表面。在筒体和感应线圈之间敷设一层20 mm厚的硅酸铝纤维保温棉，用于保温和防护线圈受热。考虑到线圈和工件的对称性，为了简化模型提高计算效率，只取1/4筒体进行计算。线圈、保温棉和筒体的相对位置，如图1所示。

　　1.2 数学模型和边界条件

　　感应加热过程中，工件中涡流产生的单位体积发热功率为：

　　1.3 材料属性

　　电磁感应加热过程中各个材料的物理参数随着温度变化呈非线性变化，因此需要不断更新相关材料的物理属性。感应线圈的相对磁导率为0.99 H/m，电阻率为1.75×10-6 Ω·m。保温棉和空气的相对电阻率为恒定值1，热导率见表1，SA5083钢的热物理性能参数见表2。

　　1.4 单元类型及网格划分

　　物理模型建立后，需要进行单元类型选择和网格划分。在计算电磁场时，筒体、线圈和保温棉网格的单元类型均设为Solid236。由于只研究筒体的温度分布，在计算温度场时，将不考虑部分的网格设置为空单元Null，将筒体网格设为Solid70。感应加热时，感应涡流仅透入工件几毫米深，热量沿工件横向传导并加热工件，因此需细化工件沿横向划分的网格，而线圈和保温棉的网格可以疏松一些。模型整体的网格划分为108 925个节点和80 426个单元，如图2所示。

　　1.5 电磁-热耦合计算

　　划分完网格需要对模型加载计算条件，文中使用的是APDL命令流的形式进行模拟分析，主要的加载参数是电流大小、频率和加热时间。首先设置磁力平衡条件，如果没有error提示，模型进入电磁场分析环节，当电磁场收敛后进入温度场耦合求解，否则继续磁场计算，直到收敛。计算结束后，得到温度场(.rth)和电磁场(.rst)两个文件。

　　2 感应加热温度场模拟结果分析

　　2.1 感应磁场

　　感应磁场通过控制感应涡流的大小影响工件温度场，因此分析温度场前，要先分析电磁场。对线圈通入电流215 A、电流频率10 kHz、加热3 600 s后，分析线圈15匝时对磁场分布的影响规律，如图3所示。由图3a可知，磁场以辐射状向四周扩散，强度由中间往四周逐渐升高，在两端达到最高。由图3b可知，磁场主要集中在工件的表面。

　　2.2 感应涡流

　　感应加热过程中，筒体升温的内热源是感应涡流，涡流回路产生焦耳热对筒体进行加热。当线圈电流215 A、电流频率10 kHz、线圈匝数15匝时，加热3 600 s后，筒体感应涡流的分布，如图4所示。

　　通过比较图3和图4可知，涡流和磁场的分布形状和线圈形状相似，其中涡流在筒体外壁表面形成回路产生热量，且中心位置的涡流较小。因此，线圈结构决定磁场和涡流形状，从而决定筒体加热后的温度分布。

　　2.3 温度场分析

　　线圈结构决定了筒体加热后温度分布。因此本节研究线圈匝数、线圈电流、电流频率和线圈与筒体加热距离对温度分布的影响，通过参数优化，使筒体温度达到250～400 ℃的消氢技术要求。

　　2.3.1 线圈匝数对温度场的影响

　　线圈匝数是控制线圈加热范围和加热均匀性的重要参数，研究在其它参数不变的情况下，线圈匝数对温度场的影响。图5为不同匝数下筒体外壁温度分布。所选取线圈电流215 A，电流频率10 kHz，分别研究线圈匝数为5匝、10匝和15匝时，加热3 600 s后，筒体表面温度分布结果。由图5可知，线圈匝数为5匝时，筒体中心温度为140 ℃左右，四周温度高于169 ℃，高温区温度达到187 ℃，温度分布不均匀，说明感应涡流没有形成较好的回路。线圈的匝数越大，筒体的加热温度越高，温度分布越规则。当线圈匝数为15匝时，筒体的高温区集中在线圈正下方，最高温度为360 ℃，四周温度最低为280 ℃，筒体温度满足焊后消氢工艺的要求。

　　为了便于分析，将筒体沿轴向剖开，并沿轴向依次取10个均匀分布于筒体表面的测温点，观察筒体轴向的温度分布曲线，如图6所示。由图6可知，线圈匝数为5匝时，筒体温度偏低，在200 ℃以下;线圈匝数为10匝时，筒体温度有一部分在250 ℃以下;线圈匝数为15匝时，筒体温度在250～400 ℃范围内，满足焊后消氢的工艺要求。

　　综上所述，线圈的匝数越多，在筒体表

　　面产生的磁通密度越大，感应涡流越大，工件温度升高，筒体中心区域温度扩展越快，加热效率越高。

　　2.3.2 线圈电流对温度场的影响

　　线圈电流也是电磁感应加热的重要参数，通过改变电流的大小可以控制感应加热的时间和温度。上节分析表明，线圈匝数为15匝的加热效果最佳。故本节选取线圈匝数15匝、电流频率10 kHz、线圈电流分别为115 A，165 A，215 A，265 A和315 A，加热3 600 s后，分析不同电流时筒体外壁和内壁感应加热温度分布的规律，如图7～8所示。由图7和图8可知，线圈电流115 A时，筒体表面温度很低。随着线圈电流增大，内壁面温度越来越高，温度分布也逐渐均匀;当线圈电流为215 A时，筒体高温区温度为362 ℃，低温区温度为271 ℃，内、外壁温度分布均匀;而当线圈电流增加到315 A时，温度分布不均匀且已经超过了工艺要求。

　　推荐阅读：《电焊机》(月刊)创刊于1971年，由成都电焊机研究所主办。

转载请注明来自：http://www.qikan2017.com/lunwen/lig/14846.html