电渣重熔空心铸锭温度场与热应力的数值模拟及影响因素探究

电渣重熔空心铸锭温度场与热应力的数值模拟及影响因素探究一、引言1.1研究背景与意义在现代工业生产中，电渣重熔空心铸锭凭借其独特的性能优势，在众多领域得到了广泛应用。电渣重熔技术利用电流通过熔渣产生的电阻热作为热源，将自耗电极熔化，金属液滴在穿过渣池的过程中实现提纯和精炼，最终在水冷结晶器中凝固形成铸锭。这种工艺能够有效去除金属中的有害杂质元素和非金属夹杂物，降低气体含量，显著改善铸锭的结晶组织，使其成分更加均匀，组织更加致密。空心铸锭由于其内部空心的结构特点，在减轻重量的同时，还能提高材料的比强度和比刚度，因此在航空航天、能源、机械制造等领域具有重要的应用价值。例如，在航空航天领域，空心铸锭被用于制造飞机发动机的轴类零件、火箭的壳体等关键部件，能够有效减轻飞行器的重量，提高其性能和燃油效率；在能源领域，空心铸锭可用于制造核电站的压力管道、风力发电机的主轴等，满足其在高温、高压、高应力等恶劣工况下的使用要求；在机械制造领域，空心铸锭可用于制造大型机械的曲轴、液压缸等零件，提高零件的综合性能和使用寿命。然而，电渣重熔空心铸锭的生产过程涉及到复杂的多物理场耦合现象，包括电磁场、温度场、流场以及凝固过程等，这些物理场之间相互作用、相互影响，使得铸锭的质量受到多种因素的制约。例如，温度场的分布不均匀会导致铸锭内部出现缩孔、疏松等缺陷；冷却过程中的热应力过大则可能引发铸锭的开裂，严重影响铸锭的质量和性能。传统的实验研究方法虽然能够直观地获取一些实验数据，但受到实验条件和测试手段的限制，难以全面深入地揭示电渣重熔空心铸锭过程中各物理量的分布规律和变化趋势。此外，实验研究成本高、周期长，且难以对一些关键参数进行精确控制和独立研究，这在一定程度上制约了电渣重熔技术的进一步发展和优化。随着计算机技术和数值计算方法的飞速发展，数值模拟作为一种强大的研究工具，为电渣重熔空心铸锭的研究提供了新的途径。通过建立合理的数学模型，利用数值计算方法对电渣重熔空心铸锭过程进行模拟，可以在计算机上重现整个重熔过程，获得丰富的物理信息，如温度分布、电流密度分布、流速分布以及凝固界面的移动等。数值模拟不仅能够弥补实验研究的不足，还可以为实验方案的设计提供理论指导，降低实验成本和周期。通过数值模拟，能够深入分析各工艺参数（如电流强度、电压、电极插入深度、渣系成分等）对电渣重熔空心铸锭过程和产品质量的影响规律，从而为优化电渣重熔工艺参数、提高铸锭质量提供科学依据，具有重要的工程应用价值。综上所述，开展电渣重熔空心铸锭的温度场与冷却过程热应力的数值模拟研究，对于深入理解电渣重熔空心铸锭过程的物理本质，揭示各物理场之间的耦合机制，优化电渣重熔工艺参数，提高空心铸锭的质量和性能，推动相关工业领域的发展具有重要的理论意义和实际应用价值。1.2国内外研究现状电渣重熔技术自20世纪40年代提出以来，在国内外得到了广泛的研究与应用。早期的研究主要集中在工艺探索和设备开发上，随着计算机技术和数值计算方法的兴起，数值模拟逐渐成为研究电渣重熔过程的重要手段。在国外，诸多科研团队和学者在电渣重熔数值模拟领域取得了丰硕成果。美国、俄罗斯、德国、日本等国家的研究起步较早，他们对电渣重熔过程中的电磁场、温度场、流场等进行了深入研究。例如，美国的一些研究机构通过建立三维数学模型，对电渣重熔过程中的电流分布、焦耳热产生以及温度场演变进行了详细模拟，分析了不同工艺参数对温度场的影响规律，发现电流强度和电极插入深度的变化会显著改变渣池和金属熔池的温度分布，进而影响铸锭的凝固过程和质量。俄罗斯的学者则侧重于研究电渣重熔过程中的冶金反应和凝固组织形成，通过数值模拟揭示了熔渣成分、温度场与凝固组织之间的内在联系，为优化电渣重熔工艺提供了理论依据。德国的研究团队利用先进的计算流体力学软件，对电渣重熔过程中的流场进行了精确模拟，分析了金属液滴在渣池中的运动轨迹和传热传质过程，以及流场对温度场均匀性的影响。在国内，电渣重熔技术的研究也取得了长足进展。东北大学、北京科技大学、上海大学等高校在电渣重熔数值模拟方面开展了大量研究工作。东北大学的研究团队针对不同钢种的电渣重熔空心铸锭过程，建立了考虑多物理场耦合的数学模型，研究了工艺参数对温度场和凝固组织的影响，发现合适的渣系成分和电流参数能够改善温度场分布，减少铸锭内部的缺陷。北京科技大学的学者通过数值模拟与实验相结合的方法，研究了电渣重熔过程中的热应力分布规律，分析了冷却速度、铸锭尺寸等因素对热应力的影响，提出了通过优化工艺参数来降低热应力、防止铸锭开裂的措施。上海大学的研究人员则对电渣重熔空心铸锭过程中的电磁搅拌作用进行了研究，通过数值模拟分析了电磁搅拌对温度场和流场的调控效果，以及对铸锭质量的改善作用。然而，目前国内外关于电渣重熔空心铸锭的温度场与冷却过程热应力的数值模拟研究仍存在一些不足之处。一方面，虽然已经建立了多种数学模型来描述电渣重熔过程，但这些模型往往对一些复杂的物理现象进行了简化，如渣金界面的传热传质、熔渣的物化性质变化等，导致模拟结果与实际情况存在一定偏差。另一方面，对于多物理场之间的耦合机制研究还不够深入，尤其是电磁场、温度场、流场以及热应力场之间的相互作用关系尚未完全明确。此外，实验研究与数值模拟的结合还不够紧密，缺乏足够的实验数据来验证和完善数值模型，限制了数值模拟结果的准确性和可靠性。综上所述，尽管国内外在电渣重熔空心铸锭的数值模拟研究方面取得了一定成果，但仍有许多问题有待进一步深入研究和解决。开展更加精确、全面的数值模拟研究，深入揭示电渣重熔空心铸锭过程中的物理本质和多物理场耦合机制，对于提高电渣重熔工艺水平和空心铸锭质量具有重要意义。1.3研究内容与方法1.3.1研究内容本研究旨在深入探究电渣重熔空心铸锭过程中的温度场分布规律以及冷却过程热应力的演变情况，具体研究内容如下：建立电渣重熔空心铸锭的数学模型：综合考虑电渣重熔过程中的电磁场、温度场、流场以及凝固过程等多物理场耦合现象，基于电磁学、传热学、流体力学以及金属凝固理论等相关知识，建立三维瞬态数学模型。该模型需准确描述各物理场的控制方程以及它们之间的耦合关系，为后续的数值模拟提供理论基础。例如，根据麦克斯韦方程组确定电磁场的分布，依据能量守恒定律建立温度场的控制方程，利用Navier-Stokes方程描述流场的运动，同时考虑金属凝固过程中的潜热释放对温度场的影响。模拟电渣重熔空心铸锭的温度场：运用建立的数学模型，对电渣重熔空心铸锭过程中的温度场进行数值模拟。分析在不同工艺参数（如电流强度、电压、电极插入深度、渣系成分等）下，渣池、金属熔池以及铸锭的温度分布随时间的变化规律。研究温度场的不均匀性对铸锭凝固过程和质量的影响，如温度梯度对凝固界面形态和晶粒生长方向的影响，以及局部高温或低温区域可能导致的缩孔、疏松等缺陷。分析冷却过程热应力：在温度场模拟的基础上，考虑材料的热物理性能随温度的变化以及铸锭的几何形状和尺寸，建立热应力计算模型。模拟电渣重熔空心铸锭冷却过程中的热应力分布，分析热应力的产生机制和影响因素。研究热应力过大导致铸锭开裂的风险，以及如何通过优化工艺参数来降低热应力，提高铸锭的质量和可靠性。例如，分析冷却速度、铸锭的约束条件等因素对热应力大小和分布的影响。实验验证与模型优化：设计并开展电渣重熔空心铸锭实验，通过在实验过程中布置温度传感器和应力测量装置，获取实际的温度和应力数据。将实验数据与数值模拟结果进行对比分析，验证数学模型的准确性和可靠性。根据实验结果对模型进行优化和修正，进一步提高模型的精度，使其能够更准确地预测电渣重熔空心铸锭过程中的温度场和热应力分布。1.3.2研究方法本研究采用数值模拟与实验研究相结合的方法，具体如下：数值模拟方法：采用有限元法对电渣重熔空心铸锭过程进行数值模拟。有限元法是一种将连续的求解域离散为有限个单元的数值计算方法，它能够有效地处理复杂的几何形状和边界条件。在本研究中，利用有限元软件对电渣重熔系统进行网格划分，将其离散为大量的小单元。通过对每个单元的控制方程进行求解，得到整个求解域内的物理量分布。在数值模拟过程中，对不同的物理场采用相应的控制方程和边界条件。例如，对于电磁场，根据麦克斯韦方程组建立控制方程，并考虑电流的输入和边界的电磁条件；对于温度场，依据传热学原理建立能量守恒方程，并考虑渣池与金属熔池之间的传热、结晶器的水冷边界条件以及铸锭与周围环境的热交换等。通过迭代计算，求解多物理场耦合问题，得到电渣重熔空心铸锭过程中的温度场和热应力分布。实验研究方法：设计并搭建电渣重熔实验装置，进行电渣重熔空心铸锭实验。实验装置包括电渣炉、自耗电极、水冷结晶器、渣系供给系统以及电气控制系统等。选择合适的实验材料和渣系，按照设定的工艺参数进行实验。在实验过程中，利用热电偶测量渣池、金属熔池和铸锭不同位置的温度，采用应变片或应力传感器测量铸锭冷却过程中的应力。同时，对铸锭进行宏观和微观组织分析，观察其凝固组织和缺陷情况。通过实验，获取实际的温度、应力和组织数据，为数值模拟结果的验证和模型优化提供依据。软件工具：选用专业的多物理场耦合模拟软件，如ANSYS、COMSOLMultiphysics等，进行电渣重熔空心铸锭过程的数值模拟。这些软件具有强大的物理建模和数值计算功能，能够方便地实现电磁场、温度场、流场等多物理场的耦合模拟。在使用软件时，根据电渣重熔过程的特点和研究需求，合理设置模型参数和边界条件，确保模拟结果的准确性和可靠性。例如，在ANSYS软件中，利用其电磁分析模块建立电磁场模型，利用热分析模块建立温度场模型，通过耦合场分析实现多物理场的耦合计算；在COMSOLMultiphysics软件中，通过选择相应的物理场接口和模块，构建电渣重熔过程的多物理场耦合模型，并进行数值求解。二、电渣重熔空心铸锭基本原理与过程2.1电渣重熔技术原理电渣重熔技术是一种利用电流通过熔渣产生电阻热作为热源，对金属进行二次熔炼的先进特种熔炼方法。其基本工作原理如图1所示，在一个铜制水冷结晶器内，盛有处于熔融状态的炉渣。自耗电极通常为经过初炼的金属铸锭或锻坯，其一端插入到熔渣之中。自耗电极、渣池、金属熔池、钢锭以及底水箱通过短网导线和变压器形成一个闭合回路。当电流通过渣池时，由于熔渣具有一定的电阻，根据焦耳定律Q=I^2Rt（其中Q为焦耳热，I为电流，R为渣池电阻，t为时间），电流通过渣池产生焦耳热，使得渣池温度迅速升高，一般可达1600℃-1800℃，该温度足以将自耗电极端头逐渐熔化。熔化后的金属液滴在重力的作用下，穿过渣池，落入结晶器底部，形成金属熔池。由于结晶器采用水冷方式，金属熔池底部的液态金属在强制冷却作用下迅速凝固，逐渐形成钢锭。在整个重熔过程中，电极的熔化和金属的凝固同时进行，渣池始终覆盖在金属熔池上方，发挥着多重重要作用。首先，渣池作为发热体，将电能转化为电阻热，为金属的熔化和过热提供所需的能量。其次，渣池起到净化剂的作用，在重熔过程中，熔化的金属形成金属熔滴穿过渣池时，渣与金属液滴接触面积很大，利用渣的化学特性可以有效地去除钢中有害元素硫、磷等，去除有害气体氢、氧、氮等，吸收、溶解熔融金属与渣界面的非金属夹杂物，使金属得到提纯、净化。再者，渣池还充当保护剂，在铸锭和结晶器间形成一层渣皮，这层渣皮起到了绝缘、隔热、润滑作用，使金属不与结晶器直接接触，防止了电分流，使热流主要向底水箱方向传导，结晶趋向纵轴方向，并有利于铸锭和结晶器壁的相对运动。另外，在金属熔池上方的熔融渣池使金属液避免直接与大气接触，起到了防止金属氧化及贮热保温作用。在电渣重熔过程中，金属始终在渣的包覆下熔化、凝固。在电渣重熔过程中，熔渣的选择至关重要，它直接影响到重熔过程的稳定性和铸锭的质量。熔渣应具备多种特性，以满足重熔过程的不同需求。首先，熔渣需具有适当的电导率，以产生足够的热量，保证金属熔化、过热及精炼的进行，并提高电渣重熔的电效率，降低电耗。一般要求在2000℃时，熔渣的电导k≤3Ω^{-1}・cm^{-1}。其次，熔渣的熔点应低于重熔金属熔点，熔渣成分力求选在低熔共晶点附近，这样可以减少渣皮凝固时的偏析，防止渣成分的变化造成工艺参数难以稳定的现象。通常渣的熔点应比重熔金属的熔点低100-200℃。此外，熔渣应有适当的黏度，并且在高温下黏度值随温度变化越小越好。具有良好流动性的渣才能保证高温下渣池热对流，使铸锭径向温度均匀，利于去气脱硫等物化反应的进行，利于形成薄渣皮，提高铸锭表面质量。一般熔渣在1800℃时，黏度η≤0.5Pa・s。为了保证重熔过程脱硫良好，熔渣还应具有较高的碱度，一般碱度B＞1，B=CaO/(SiO_{2}+1/2Al_{2}O_{3})。同时，熔渣与熔融金属之间应有足够大的表面张力，以有助于渣壳分离，防止熔渣卷入金属内；而熔渣与非金属夹杂物之间的表面张力越小越好，以利于熔渣对非金属夹杂物的浸润、吸附和溶解。常用渣的基本形态是以CaF_{2}为基础，配入适当的CaO、Al_{2}O_{3}、MgO、SiO_{2}、TiO_{2}等氧化物组成。例如，CaF_{2}能降低渣的熔点、黏度和表面张力，但和其他组元相比，CaF_{2}的电导率较高，纯CaF_{2}在1650℃时电导率达4.54Ω^{-1}・cm^{-1}，且渣中CaF_{2}含量高，熔炼中易放出有害气体和烟尘，造成环境污染；CaO可增大渣的碱度，提高脱硫效率，在CaO加入量为40\%情况下，脱硫率最高可达到85\%，而且CaO的加入能够降低渣的电导率，但CaO吸水性强，易带入氢和氧，造成钢增氢增氧；Al_{2}O_{3}能明显降低渣的电导率，减少电耗，提高生产率，但渣中Al_{2}O_{3}增加，将使渣的熔化温度和黏度升高，并将降低渣的脱硫效果，另外会使重熔过程难以建立和稳定，一般Al_{2}O_{3}的含量不大于50\%。电渣重熔技术具有诸多显著优点。在精炼提纯方面效果显著，能够有效去除钢中的有害杂质元素和非金属夹杂物，降低钢中的气体含量，提高钢的纯净度。例如，经过电渣重熔处理后，钢中的硫含量可降低至0.005\%以下，非金属夹杂物的数量和尺寸也显著减少。同时，该技术能使钢锭成分均匀性高，减少宏观偏析的产生。此外，电渣重熔过程中形成的定向凝固组织优良，水冷结晶器的强制冷却作用使得金属熔池底部首先凝固，然后凝固界面逐渐向上移动，形成自下而上的定向凝固方式，这种定向凝固组织使得钢锭内部的晶粒沿着热流方向生长，组织致密，柱状晶发达，显著提高了钢锭的力学性能，尤其是韧性和抗疲劳性能。但该技术也存在一些缺点，如渣成分不在低熔共晶点上，熔点偏高；渣的电导率高，重熔电耗高；渣中CaF_{2}含量高，熔炼中易污染环境；钢中夹杂物以脆性铝酸盐及刚玉为主，影响钢的塑性和韧性等。为了克服这些缺点，近年来开发了多种多元渣系，如CaF_{2}-CaO-Al_{2}O_{3}三元渣系，CaF_{2}-CaO-Al_{2}O_{3}-MgO四元渣系及CaF_{2}-CaO-Al_{2}O_{3}-MgO-SiO_{2}五元渣系等。这些渣系的成分基本选择在三元共晶线上，渣的熔点较低，渣皮凝固时偏析现象较轻微，渣成分稳定，铸锭表面光洁，成型良好。2.2空心铸锭的制备过程空心铸锭的制备是在特定的电渣重熔设备中进行，该设备主要由水冷结晶器、自耗电极、渣系供给系统、电气控制系统等关键部分组成。水冷结晶器是整个制备过程的核心部件之一，通常采用铜材质制成，其具有良好的导热性能，能够实现高效的冷却作用。结晶器的内部形状与所需制备的空心铸锭的外形相适配，一般为环形结构。在结晶器的内部，还设置有用于支撑和引导自耗电极的装置，确保自耗电极能够准确地插入到渣池中。自耗电极通常是由经过初炼的金属材料制成，其化学成分和性能需满足空心铸锭的质量要求。渣系供给系统负责将预先配置好的渣料输送到结晶器内，形成具有特定成分和性能的渣池。电气控制系统则用于控制电流、电压等电参数，确保电渣重熔过程的稳定进行。在空心铸锭的制备过程中，首先要进行渣池的准备工作。将按特定比例配置好的渣料加入到水冷结晶器内，通过电极通电产生的焦耳热将渣料熔化，形成高温的渣池。在熔化过程中，需要密切监控渣池的温度、电导率等参数，确保渣池达到合适的工作状态。一般来说，渣池的温度需达到1600℃-1800℃，以保证后续自耗电极的熔化和金属液滴的精炼。例如，在熔化含CaF_{2}、CaO、Al_{2}O_{3}等成分的渣料时，CaF_{2}能降低渣的熔点和黏度，使渣料更容易熔化，但渣中CaF_{2}含量过高会导致电导率升高，增加电耗；CaO可增大渣的碱度，提高脱硫效率，但吸水性强，易带入氢和氧；Al_{2}O_{3}能降低渣的电导率，减少电耗，但会使渣的熔化温度和黏度升高。因此，需要根据实际情况精确控制渣料的成分比例。渣池准备就绪后，将自耗电极缓慢插入到渣池中，使其一端与渣池充分接触。接通电源后，电流通过自耗电极、渣池、金属熔池、钢锭以及底水箱形成闭合回路。根据焦耳定律Q=I^2Rt，电流通过渣池时产生大量的焦耳热，使渣池温度急剧升高。在高温作用下，自耗电极端头开始逐渐熔化，形成金属液滴。这些金属液滴在重力的作用下，穿过渣池，落入结晶器底部，汇聚形成金属熔池。在金属液滴穿过渣池的过程中，渣池与金属液滴之间发生一系列复杂的物理化学反应。渣池中的熔渣能够吸附和溶解金属液滴中的有害杂质元素，如硫、磷等，同时也能去除金属液滴中的有害气体，如氢、氧、氮等，从而实现对金属液滴的精炼和净化。例如，在渣钢界面上，硫元素会与熔渣中的碱性氧化物发生反应，生成硫化物进入渣相，降低钢中的硫含量。随着自耗电极的不断熔化和金属液滴的持续落入，金属熔池的体积逐渐增大。由于结晶器采用水冷方式，金属熔池底部的液态金属在强制冷却作用下迅速凝固，形成铸锭的初始凝固层。随着凝固过程的进行，凝固层逐渐向上生长，铸锭的高度不断增加。在整个凝固过程中，需要严格控制冷却速度和温度分布，以确保铸锭的质量。如果冷却速度过快，可能导致铸锭内部产生应力集中，引发裂纹等缺陷；如果冷却速度过慢，则会影响生产效率，且可能导致铸锭组织粗大。同时，温度分布不均匀也会导致铸锭内部出现缩孔、疏松等缺陷。为了控制冷却速度和温度分布，通常会在结晶器的不同部位设置多个冷却通道，通过调节冷却介质的流量和温度来实现对铸锭冷却过程的精确控制。在空心铸锭的制备过程中，还需要注意一些关键环节。一是要确保自耗电极的稳定熔化和金属液滴的均匀滴落，这需要精确控制电流、电压等电参数，以及自耗电极的插入深度和下降速度。二是要保证渣池的稳定性和均匀性，避免渣池出现波动或局部过热等情况。三是要及时清理渣池表面的浮渣和杂质，防止其进入金属熔池，影响铸锭质量。四是要对铸锭进行实时监测，通过布置温度传感器、应力传感器等设备，获取铸锭在制备过程中的温度、应力等数据，以便及时调整工艺参数，确保铸锭质量。2.3工艺参数对重熔过程的影响在电渣重熔空心铸锭过程中，电流强度、电压、电极插入深度、渣系成分等工艺参数对重熔过程有着至关重要的影响，它们直接关系到铸锭的质量和性能。电流强度是影响电渣重熔过程的关键参数之一。根据焦耳定律Q=I^2Rt，电流强度的变化会显著影响渣池产生的焦耳热，进而改变渣池和金属熔池的温度分布。当电流强度增大时，渣池产生的焦耳热增加，渣池和金属熔池的温度升高。例如，在对G20轴承钢的电渣重熔研究中发现，随着电流强度的增加，整个体系电流密度增大，温度提高。然而，过高的电流强度可能导致渣池和金属熔池温度过高，使金属液滴的过热程度增加，这不仅会增加金属中元素的挥发损失，还可能导致铸锭内部组织粗大，出现缩孔、疏松等缺陷。相反，若电流强度过小，渣池产生的热量不足，会使金属熔化速度减慢，生产效率降低，同时可能导致金属液滴在渣池中停留时间过长，增加夹杂物的吸收几率。此外，电流强度还会影响金属熔滴的形成和滴落行为。研究表明，熔滴的尺寸和数量随着输入电流的增大而减小，一些尺寸较大的熔滴在滴落过程中受到电磁力、重力以及浮力等力的作用会被分裂成几个尺寸较小的熔滴。电压作为电渣重熔过程中的另一个重要电参数，同样对重熔过程产生显著影响。电压的变化会改变渣池的电阻和电流分布，从而影响焦耳热的产生和温度场的分布。在恒电流操作下，提高电压会使渣池电阻增大，焦耳热增加，渣池和金属熔池的温度升高；反之，降低电压则会使温度降低。电压还会影响电渣重熔过程的稳定性。如果电压过高，可能会导致渣池内的电流分布不均匀，出现局部过热现象，引发渣池的波动和喷溅，影响重熔过程的正常进行；而电压过低则可能使重熔过程难以稳定维持，金属熔化速度不稳定，影响铸锭的质量。例如，在实际生产中，当电压波动较大时，铸锭的表面质量会受到明显影响，可能出现表面不平整、渣皮厚度不均匀等问题。电极插入深度对电渣重熔过程也有着重要作用。电极插入深度的改变会影响电流在渣池中的分布路径和渣池的有效电阻。随着电极插入深度的增加，电流在渣池中的分布区域发生变化，渣池的有效电阻减小。相关研究表明，当电极插入深度从10mm增加到30mm时，体系的温度呈下降趋势。这是因为电极插入深度增加，电流通过渣池的路径缩短，渣池产生的焦耳热减少，导致体系温度降低。电极插入深度还会影响金属熔滴的形成和下落过程。如果电极插入过深，金属熔滴在渣池中下落的距离增加，与渣池的接触时间变长，可能会导致金属熔滴过度精炼，一些有益元素被过多去除；而电极插入过浅，则可能使金属熔滴不能充分与渣池接触，无法有效去除夹杂物和有害元素。合适的电极插入深度对于保证电渣重熔过程的稳定性和铸锭质量至关重要，在实际生产中需要根据具体情况进行优化调整。渣系成分是影响电渣重熔过程的关键因素之一，不同的渣系成分会导致渣池具有不同的物理化学性质，从而对重熔过程产生多方面的影响。渣系的电导率直接关系到渣池产生焦耳热的能力。例如，CaF_{2}的电导率较高，纯CaF_{2}在1650℃时电导率达4.54Ω^{-1}・cm^{-1}，在渣系中增加CaF_{2}的含量会提高渣池的电导率，使渣池产生更多的焦耳热，提高金属的熔化速度，但同时也会增加电耗；而Al_{2}O_{3}能明显降低渣的电导率，减少电耗，例如在CaF_{2}-Al_{2}O_{3}渣系中，当Al_{2}O_{3}含量从10%增加到30%时，电导率将从3.44Ω^{-1}・cm^{-1}降为1.75Ω^{-1}・cm^{-1}，但渣中Al_{2}O_{3}增加会使渣的熔化温度和黏度升高，影响渣池的流动性和重熔过程的稳定性。渣系的熔点和黏度对重熔过程也有着重要影响。渣的熔点应低于重熔金属熔点，通常渣的熔点应比重熔金属的熔点低100-200℃，这样可以保证渣池在合适的温度下保持液态，有效发挥其作用。熔渣应有适当的黏度，一般熔渣在1800℃时，黏度η≤0.5Pa・s，良好的流动性有助于渣池内的热对流，使铸锭径向温度均匀，利于去气脱硫等物化反应的进行，并有利于形成薄渣皮，提高铸锭表面质量。例如，在渣系中加入适量的CaF_{2}可以降低渣的熔点和黏度，但含量过高会导致其他问题，如污染环境等。渣系的碱度和表面张力也会影响重熔过程。为保证重熔过程脱硫良好，熔渣应具有较高的碱度，一般碱度B＞1，熔渣与熔融金属之间应有足够大的表面张力，以有助于渣壳分离，防止熔渣卷入金属内；而熔渣与非金属夹杂物之间的表面张力越小越好，以利于熔渣对非金属夹杂物的浸润、吸附和溶解。三、温度场数值模拟3.1物理模型建立为了准确模拟电渣重熔空心铸锭过程中的温度场，需要构建一个合理的物理模型。该模型涵盖了电渣重熔系统中的关键组成部分，包括自耗电极、渣池、金属熔池、空心铸锭和结晶器，这些部分相互关联，共同决定了电渣重熔过程的物理现象。自耗电极通常由经过初炼的金属材料制成，其化学成分和性能对最终铸锭的质量有着重要影响。在电渣重熔过程中，自耗电极作为金属的来源，在渣池产生的高温作用下逐渐熔化，为金属熔池提供液态金属。自耗电极的形状一般为圆柱形，其直径和长度根据实际生产需求和设备条件进行选择。例如，在某些大型电渣重熔设备中，自耗电极的直径可达数百毫米，长度数米，以满足大规模生产的要求。渣池是电渣重熔过程的核心区域之一，它不仅是产生焦耳热的热源，还在金属液滴的精炼和铸锭的凝固过程中发挥着重要作用。渣池由特定成分的熔渣组成，这些熔渣具有合适的电导率、熔点、黏度和碱度等物理化学性质，以确保电渣重熔过程的顺利进行。如前文所述，常用的渣系以CaF_{2}为基础，配入适量的CaO、Al_{2}O_{3}、MgO等氧化物，通过调整各成分的比例来满足不同的工艺要求。渣池的形状近似为环形，环绕在金属熔池周围，其高度和宽度也会根据具体的工艺参数和设备结构而有所不同。金属熔池是液态金属汇聚的区域，在电渣重熔过程中，自耗电极熔化形成的金属液滴穿过渣池后落入金属熔池，随着重熔过程的进行，金属熔池的体积逐渐增大。金属熔池的形状和尺寸受到多种因素的影响，如电流强度、电极插入深度、渣池的热传递等。在重力和电磁力的作用下，金属熔池内的液态金属会发生流动，这种流动对温度场的分布和铸锭的凝固过程有着重要影响。例如，金属熔池内的对流运动会使温度分布更加均匀，促进热量的传递，影响凝固界面的形态和晶粒的生长方向。空心铸锭是电渣重熔的最终产品，其内部为空心结构，外部为凝固的金属层。空心铸锭的形状和尺寸根据实际应用需求而定，通常具有一定的壁厚和内径。在电渣重熔过程中，空心铸锭从底部开始逐渐凝固，随着金属熔池内液态金属的不断补充，铸锭的高度逐渐增加。空心铸锭的凝固过程是一个复杂的传热和相变过程，涉及到热量的传递、潜热的释放以及晶体的生长等多个物理现象。结晶器是实现空心铸锭凝固的关键部件，通常采用水冷方式，以提供快速的冷却条件，促使金属熔池底部的液态金属迅速凝固。结晶器一般由铜材质制成，具有良好的导热性能。其内部形状与空心铸锭的外形相匹配，为空心铸锭的凝固提供了一个特定的空间。在结晶器的内壁与空心铸锭之间，会形成一层渣皮，这层渣皮不仅起到了隔热和润滑的作用，还对铸锭的表面质量和凝固过程产生影响。例如，渣皮的厚度和均匀性会影响铸锭的冷却速度和温度分布，进而影响铸锭的表面平整度和内部组织。在构建物理模型时，需要考虑这些部件之间的相互作用和边界条件。例如，自耗电极与渣池之间存在着热量传递和质量传递，渣池与金属熔池之间发生着复杂的物理化学反应和热交换，金属熔池与空心铸锭之间存在着凝固界面的移动和热量的传导，空心铸锭与结晶器之间通过渣皮进行热量传递和力的作用。同时，还需要考虑外界环境对模型的影响，如结晶器的水冷边界条件、铸锭与周围环境的热辐射和对流换热等。通过合理地处理这些相互作用和边界条件，可以使构建的物理模型更加接近实际的电渣重熔过程，为后续的数值模拟提供可靠的基础。3.2数学模型及控制方程在电渣重熔空心铸锭过程中，涉及到电磁场、温度场等多个物理场，这些物理场相互作用、相互影响，共同决定了铸锭的质量和性能。为了准确模拟这一复杂过程，需要建立相应的数学模型，并给出各物理场的控制方程。3.2.1电磁场数学模型及控制方程电渣重熔过程中的电磁场是产生焦耳热的根源，其分布决定了渣池和金属熔池的加热方式和热量分布。在电渣重熔体系中，电磁场的行为遵循麦克斯韦方程组，其在低频情况下的形式如下：\nabla\times\vec{H}=\vec{J}+\frac{\partial\vec{D}}{\partialt}\tag{1}\nabla\times\vec{E}=-\frac{\partial\vec{B}}{\partialt}\tag{2}\nabla\cdot\vec{B}=0\tag{3}\nabla\cdot\vec{D}=\rho\tag{4}其中，\vec{H}是磁场强度（A/m），\vec{J}是电流密度（A/m²），\vec{D}是电位移矢量（C/m²），\vec{E}是电场强度（V/m），\vec{B}是磁感应强度（T），\rho是电荷密度（C/m³），t是时间（s）。在电渣重熔过程中，由于电流频率较低，位移电流\frac{\partial\vec{D}}{\partialt}通常可以忽略不计。同时，假设材料为各向同性且电导率\sigma为常数，根据欧姆定律\vec{J}=\sigma\vec{E}，以及\vec{B}=\mu\vec{H}（\mu为磁导率，H/m），\vec{D}=\epsilon\vec{E}（\epsilon为介电常数，F/m），可以对麦克斯韦方程组进行简化。将\vec{J}=\sigma\vec{E}代入式（1），可得\nabla\times\vec{H}=\sigma\vec{E}。再结合式（2），消去\vec{E}，得到关于磁场强度\vec{H}的方程：\nabla\times(\frac{1}{\sigma}\nabla\times\vec{H})+\mu\frac{\partial\vec{H}}{\partialt}=0\tag{5}对于电渣重熔空心铸锭的轴对称模型，采用柱坐标系(r,\theta,z)更为方便。在柱坐标系下，磁场强度\vec{H}可以表示为\vec{H}=H_r\vec{e}_r+H_{\theta}\vec{e}_{\theta}+H_z\vec{e}_z，其中\vec{e}_r、\vec{e}_{\theta}、\vec{e}_z分别是径向、周向和轴向的单位矢量。由于轴对称性，H_{\theta}=0，且所有物理量与\theta无关。此时，式（5）在柱坐标系下的分量形式为：\frac{\partial}{\partialr}(\frac{1}{\sigma}\frac{\partialH_z}{\partialr})+\frac{1}{r}\frac{\partial}{\partialr}(\frac{r}{\sigma}\frac{\partialH_r}{\partialz})+\frac{\partial}{\partialz}(\frac{1}{\sigma}\frac{\partialH_r}{\partialz})+\mu\frac{\partialH_r}{\partialt}=0\tag{6}\frac{\partial}{\partialr}(\frac{1}{\sigma}\frac{\partialH_r}{\partialr})+\frac{1}{r}\frac{\partial}{\partialr}(\frac{r}{\sigma}\frac{\partialH_r}{\partialz})+\frac{\partial}{\partialz}(\frac{1}{\sigma}\frac{\partialH_z}{\partialz})+\mu\frac{\partialH_z}{\partialt}=0\tag{7}通过求解上述方程，可以得到磁场强度\vec{H}的分布，进而根据\vec{J}=\sigma\vec{E}=\sigma(\nabla\times\vec{H})计算出电流密度\vec{J}的分布。电流密度的分布对于确定焦耳热的产生位置和大小至关重要，焦耳热密度q_{J}可由下式计算：q_{J}=\vec{J}\cdot\vec{E}=\sigma|\vec{J}|^2\tag{8}3.2.2温度场数学模型及控制方程温度场是电渣重熔空心铸锭过程中另一个关键的物理场，它直接影响着金属的熔化、凝固以及铸锭的质量。在电渣重熔过程中，温度场的变化遵循热传导方程，考虑到焦耳热的产生以及金属凝固过程中的潜热释放，热传导方程可以表示为：\rhoC_p\frac{\partialT}{\partialt}=\nabla\cdot(k\nablaT)+q_{J}+q_{s}\tag{9}其中，\rho是材料的密度（kg/m³），C_p是材料的定压比热容（J/(kg・K)），T是温度（K），k是材料的热导率（W/(m・K)），q_{J}是焦耳热密度（W/m³），q_{s}是潜热释放率（W/m³）。在金属凝固过程中，潜热的释放会对温度场产生重要影响。通常采用焓法来处理潜热问题，定义材料的焓H为：H=H_0+\int_{T_0}^{T}C_pdT+\DeltaH_ff_s\tag{10}其中，H_0是参考温度T_0下的焓（J/kg），\DeltaH_f是材料的熔化潜热（J/kg），f_s是固相率。固相率f_s与温度T的关系可以通过材料的凝固特性曲线来确定，一般采用Scheil方程或其他合适的模型来描述。将式（10）代入式（9），得到以焓为变量的热传导方程：\rho\frac{\partialH}{\partialt}=\nabla\cdot(k\nablaT)+q_{J}\tag{11}对于电渣重熔空心铸锭的轴对称模型，在柱坐标系下，式（11）的分量形式为：\rho\frac{\partialH}{\partialt}=\frac{1}{r}\frac{\partial}{\partialr}(rk\frac{\partialT}{\partialr})+\frac{\partial}{\partialz}(k\frac{\partialT}{\partialz})+q_{J}\tag{12}在求解温度场时，需要考虑边界条件。在电渣重熔系统中，主要的边界条件包括：水冷结晶器壁面的对流换热边界条件，可表示为-k\frac{\partialT}{\partialn}=h(T-T_{\infty})，其中h是对流换热系数（W/(m²・K)），T_{\infty}是冷却介质的温度（K），n是壁面的法向方向；铸锭与周围环境的热辐射边界条件，可表示为-k\frac{\partialT}{\partialn}=\epsilon\sigma_{SB}(T^4-T_{env}^4)，其中\epsilon是铸锭表面的发射率，\sigma_{SB}是斯蒂芬-玻尔兹曼常数（5.67\times10^{-8}W/(m²·K^4)），T_{env}是周围环境的温度（K）。通过求解上述热传导方程，并结合相应的边界条件，可以得到电渣重熔空心铸锭过程中温度场的分布情况，为进一步分析铸锭的凝固过程和热应力分布提供基础。3.3边界条件与初始条件设定准确设定边界条件与初始条件是确保数值模拟结果准确性和可靠性的关键环节，它们为数学模型的求解提供了必要的约束和起点。在电渣重熔空心铸锭的数值模拟中，边界条件与初始条件主要涉及电流、电压、热流密度、温度等物理量。在电渣重熔过程中，电流通过自耗电极、渣池、金属熔池和铸锭形成闭合回路。为了确定电流在各区域的分布，需要设定电流边界条件。通常假设自耗电极顶部为电流入口，电流密度均匀分布，其值根据实验或实际生产中的电流强度以及自耗电极的横截面积来确定。例如，若已知电流强度为I，自耗电极半径为r_{e}，则电流密度J_{e}=\frac{I}{\pir_{e}^{2}}。对于铸锭底部，可将其设定为电流出口，电流均匀流出。同时，考虑到电渣重熔系统的轴对称性，对称轴上的电流密度满足轴对称条件，即轴向电流密度J_{z}关于对称轴呈反对称分布，径向电流密度J_{r}=0。电压边界条件的设定与电流边界条件密切相关。在自耗电极顶部，施加已知的电压值V_{e}，该电压值根据实际的供电情况确定。由于电流通过渣池会产生电压降，根据欧姆定律，渣池内的电压分布可通过电流密度和渣池的电导率来计算。在铸锭底部，电压设定为零电位，作为参考电位。此外，在渣池与金属熔池的界面处，由于渣金之间存在一定的电位差，需要考虑界面的电位变化，这一电位差与渣金之间的物理化学反应以及电荷转移有关，可通过相关的电化学理论进行估算。热流密度边界条件主要涉及到电渣重熔系统与外界环境之间的热交换，以及系统内部各区域之间的热传递。水冷结晶器壁面与冷却介质之间存在对流换热，其热流密度可通过牛顿冷却定律来描述：q_{conv}=h(T-T_{\infty})，其中q_{conv}为对流换热热流密度（W/m^{2}），h为对流换热系数（W/(m^{2}·K)），T为结晶器壁面温度（K），T_{\infty}为冷却介质温度（K）。对流换热系数h的取值与冷却介质的流速、性质以及结晶器的结构等因素有关，可通过实验测定或经验公式估算。铸锭表面与周围环境之间存在热辐射换热，其热流密度可根据斯蒂芬-玻尔兹曼定律计算：q_{rad}=\epsilon\sigma_{SB}(T^{4}-T_{env}^{4})，其中q_{rad}为热辐射热流密度（W/m^{2}），\epsilon为铸锭表面的发射率，\sigma_{SB}是斯蒂芬-玻尔兹曼常数（5.67×10^{-8}W/(m^{2}·K^{4})），T为铸锭表面温度（K），T_{env}为周围环境温度（K）。发射率\epsilon的取值与铸锭的材料和表面状态有关，一般在0-1之间。在渣池与金属熔池的界面处，存在着热传导和热对流的综合作用，热流密度可通过界面两侧的温度差以及界面的热阻来计算。初始条件是数值模拟开始时各物理量的初始值。在电渣重熔空心铸锭的模拟中，通常假设初始时刻自耗电极、渣池、金属熔池和铸锭的温度均匀分布，其值可根据实际生产中的开炉温度确定。例如，在实际生产中，开炉前渣池通常需要预热到一定温度，假设渣池的初始温度为T_{01}，自耗电极、金属熔池和铸锭的初始温度为T_{02}，这些初始温度值将作为温度场控制方程求解的初始条件。对于电磁场和流场，初始时刻的磁场强度和流速通常设为零。但在实际重熔过程中，随着电流的接通和金属液滴的下落，电磁场和流场会逐渐建立起来。通过合理设定边界条件与初始条件，能够更准确地模拟电渣重熔空心铸锭过程中的物理现象，为深入研究温度场和热应力分布提供可靠的基础。在实际模拟过程中，还需要根据具体的研究对象和实验数据，对边界条件和初始条件进行适当的调整和优化，以提高模拟结果的精度和可靠性。3.4数值模拟结果与分析利用所建立的温度场数学模型和数值计算方法，对电渣重熔空心铸锭过程进行数值模拟，得到不同时刻的温度场分布云图，如图1所示。通过对这些云图的分析，可以深入了解温度分布规律以及熔池形状、深度的变化情况。【此处插入不同时刻温度场分布云图，如t1时刻、t2时刻、t3时刻等，t1<t2<t3，图中需清晰标注自耗电极、渣池、金属熔池、空心铸锭和结晶器等部分，并给出温度标尺】从温度场分布云图中可以看出，在电渣重熔初期，自耗电极端部与渣池接触区域温度迅速升高，这是由于电流通过渣池产生的焦耳热主要集中在该区域。随着时间的推移，热量逐渐向渣池和金属熔池传递，渣池整体温度升高，且温度分布呈现出一定的不均匀性，靠近电极端部和结晶器壁的区域温度相对较高，而渣池中心区域温度略低。这是因为电极端部是电流的主要入口，焦耳热产生较多，而结晶器壁由于水冷作用，热量散失较快，导致靠近壁面的渣池温度升高以维持热平衡。金属熔池的温度分布也具有明显的特征。金属熔池底部温度较低，这是由于结晶器的强制冷却作用，使得底部液态金属迅速散热凝固。随着离底部距离的增加，金属熔池温度逐渐升高，在金属熔池与渣池的界面附近，温度达到较高值。这是因为金属液滴从渣池落下时，携带了渣池的热量，且渣池与金属熔池之间存在强烈的热交换。同时，金属熔池内存在着一定的对流运动，这也对温度分布产生影响，使得温度分布更加均匀。例如，在熔池内部，由于电磁力和浮力的作用，液态金属会形成对流，热的液态金属从熔池中心向上运动，冷的液态金属从熔池边缘向下运动，从而促进了热量的传递和温度的均匀化。在整个电渣重熔过程中，熔池的形状和深度不断发生变化。在重熔初期，熔池深度较浅，形状较为扁平。随着自耗电极的不断熔化和金属液滴的持续落入，熔池深度逐渐增加，形状也逐渐变得更加饱满。这是因为随着重熔的进行，金属熔池内的液态金属不断积累，使得熔池的体积增大，深度增加。同时，由于结晶器壁的约束和热传递的影响，熔池的形状也逐渐发生改变。在熔池深度增加的过程中，其底部的形状也会发生变化，从初期的较平坦逐渐变为略带弧形。这是因为底部液态金属在凝固过程中，受到周围液态金属的挤压和热传递的影响，导致凝固界面呈现出一定的弯曲。为了更直观地分析熔池深度的变化，绘制了熔池深度随时间的变化曲线，如图2所示。从曲线中可以看出，在电渣重熔开始阶段，熔池深度迅速增加，这是由于自耗电极初始熔化速度较快，大量金属液滴进入熔池。随着重熔过程的持续进行，熔池深度的增加速度逐渐减缓，最终趋于稳定。这是因为随着熔池深度的增加，金属熔池的散热面积增大，热量散失加快，同时自耗电极的熔化速度也逐渐稳定，使得熔池深度的增加速度逐渐减小，当达到一定时间后，熔池深度基本保持不变，此时电渣重熔过程进入相对稳定的阶段。【此处插入熔池深度随时间变化曲线，横坐标为时间，纵坐标为熔池深度，曲线需清晰标注各阶段变化趋势，并可适当添加数据点进行说明】熔池形状和深度的变化对空心铸锭的凝固过程和质量有着重要影响。熔池深度过深可能导致铸锭内部出现缩孔、疏松等缺陷。这是因为熔池深度过大时，底部液态金属在凝固过程中，由于补缩困难，容易形成缩孔；同时，较大的温度梯度也会导致铸锭内部组织不均匀，产生疏松缺陷。熔池形状的不规则也可能影响铸锭的凝固质量，例如熔池边缘温度过高或过低，会导致铸锭边缘凝固速度不一致，从而产生应力集中，增加铸锭开裂的风险。因此，在电渣重熔空心铸锭过程中，需要合理控制工艺参数，以获得合适的熔池形状和深度，保证铸锭的质量。四、冷却过程热应力数值模拟4.1热应力产生机制在电渣重熔空心铸锭的冷却过程中，热应力的产生是一个复杂的物理现象，其主要源于温度梯度以及材料热物理性能的变化。温度梯度是导致热应力产生的关键因素之一。在电渣重熔空心铸锭的冷却过程中，由于铸锭各部分与冷却介质（如水冷结晶器）的接触情况不同，热量传递的速率存在差异，从而导致铸锭内部形成显著的温度梯度。例如，铸锭外层靠近结晶器壁面，受到水冷作用，热量迅速散失，温度较低；而铸锭内层热量传递相对较慢，温度较高。这种温度的不均匀分布使得铸锭各部分材料的热膨胀和收缩程度不一致。根据热膨胀原理，材料的热膨胀量\DeltaL与温度变化\DeltaT成正比，即\DeltaL=L_0\alpha\DeltaT（其中L_0为材料的初始长度，\alpha为材料的线膨胀系数）。由于铸锭内外层温度不同，它们的热膨胀量也不同，外层材料收缩较多，内层材料收缩较少。然而，铸锭是一个整体，各部分材料之间相互约束，无法自由地进行热膨胀和收缩。这种相互约束使得铸锭内部产生内应力，即热应力。当热应力超过材料的屈服强度时，铸锭可能会发生塑性变形；若热应力进一步超过材料的抗拉强度，则会导致铸锭开裂。材料热物理性能随温度的变化也是热应力产生的重要原因。在电渣重熔空心铸锭的冷却过程中，材料的热导率、比热容、线膨胀系数等热物理性能并非恒定不变，而是随着温度的降低而发生显著变化。以线膨胀系数为例，许多金属材料的线膨胀系数在不同温度区间呈现出不同的数值。在高温阶段，材料的线膨胀系数相对较大，随着温度的降低，线膨胀系数逐渐减小。这种线膨胀系数的变化会导致材料在冷却过程中的热膨胀和收缩行为变得更加复杂。当铸锭内部存在温度梯度时，不同温度区域的材料由于线膨胀系数的差异，其热膨胀和收缩程度的不一致性进一步加剧，从而产生更大的热应力。热导率和比热容的变化也会影响热量在铸锭内部的传递速度和分布情况，进而对热应力的产生和分布产生影响。例如，热导率的降低会使热量传递变慢，导致温度梯度增大，从而增加热应力的产生；比热容的变化则会影响材料吸收和释放热量的能力，进而影响温度变化的速率，间接影响热应力的大小。综上所述，在电渣重熔空心铸锭的冷却过程中，温度梯度和材料热物理性能的变化相互作用，共同导致了热应力的产生。深入理解热应力的产生机制，对于准确预测铸锭冷却过程中的热应力分布，采取有效的措施降低热应力，提高铸锭质量具有重要意义。4.2热应力计算模型热应力的计算基于弹性力学理论，考虑材料的热膨胀和力学性能随温度的变化。在弹性力学中，物体的应力-应变关系由本构方程描述，对于考虑热效应的情况，其本构方程需要进行相应的修正。在小变形假设下，物体的总应变\varepsilon_{ij}由弹性应变\varepsilon_{ij}^e和热应变\varepsilon_{ij}^T组成，即\varepsilon_{ij}=\varepsilon_{ij}^e+\varepsilon_{ij}^T。其中，热应变\varepsilon_{ij}^T与温度变化\DeltaT相关，对于各向同性材料，热应变可表示为\varepsilon_{ij}^T=\alpha\DeltaT\delta_{ij}，\alpha为材料的线膨胀系数，\delta_{ij}为克罗内克符号（当i=j时，\delta_{ij}=1；当i\neqj时，\delta_{ij}=0）。根据广义胡克定律，弹性应变与应力之间的关系为：\varepsilon_{ij}^e=\frac{1}{2G}(\sigma_{ij}-\frac{\nu}{1+\nu}\sigma_{kk}\delta_{ij})\tag{13}其中，\sigma_{ij}为应力张量分量，G为剪切模量，\nu为泊松比，\sigma_{kk}=\sigma_{11}+\sigma_{22}+\sigma_{33}。将式（13）和\varepsilon_{ij}^T=\alpha\DeltaT\delta_{ij}代入\varepsilon_{ij}=\varepsilon_{ij}^e+\varepsilon_{ij}^T，得到考虑热效应的本构方程：\sigma_{ij}=2G\varepsilon_{ij}+\frac{2G\nu}{1-2\nu}\varepsilon_{kk}\delta_{ij}-(3\lambda+2G)\alpha\DeltaT\delta_{ij}\tag{14}其中，\lambda=\frac{G\nu}{1-2\nu}为拉梅常数。在电渣重熔空心铸锭的冷却过程中，温度场是随时间变化的，因此热应力也是时间的函数。为了求解热应力，需要将上述本构方程与温度场的控制方程以及相关的边界条件联立求解。在实际计算中，通常采用有限元方法将连续的铸锭离散为有限个单元，对每个单元应用上述本构方程，并根据单元之间的连接关系和边界条件建立方程组，通过数值求解该方程组得到铸锭内各点的热应力分布。在求解热应力时，还需要考虑材料的力学性能随温度的变化。例如，材料的弹性模量E、剪切模量G、泊松比\nu等在不同温度下会发生变化。一般来说，随着温度的降低，材料的弹性模量和剪切模量会增大，泊松比会有所变化。这些力学性能的变化会对热应力的大小和分布产生影响。在数值模拟中，通常通过实验测量或查阅相关材料手册获取材料在不同温度下的力学性能数据，并将其作为输入参数应用于热应力计算模型中。例如，对于某特定的合金钢材料，在高温阶段，其弹性模量可能为E_1，泊松比为\nu_1；随着温度降低，在低温阶段，弹性模量变为E_2，泊松比变为\nu_2。在计算热应力时，根据铸锭各点的实时温度，选取相应温度下的力学性能参数进行计算，以更准确地反映材料在实际冷却过程中的力学行为。4.3模拟结果与热应力分析通过数值模拟，得到了电渣重熔空心铸锭冷却过程中不同时刻的热应力分布云图，如图3所示。从云图中可以清晰地观察到热应力在铸锭内部的分布情况。【此处插入不同时刻热应力分布云图，如t1时刻、t2时刻、t3时刻等，t1<t2<t3，图中需清晰标注自耗电极、渣池、金属熔池、空心铸锭和结晶器等部分，并给出热应力标尺】在冷却初期，热应力主要集中在铸锭的外层靠近结晶器壁的区域。这是因为铸锭外层首先与水冷结晶器接触，冷却速度快，温度迅速降低，而内层温度下降相对较慢，从而在铸锭内外层之间形成较大的温度梯度。根据热应力的产生机制，这种温度梯度导致铸锭外层收缩较快，内层收缩较慢，由于内外层之间的相互约束，使得外层承受较大的拉应力，内层承受较大的压应力。随着冷却时间的延长，热应力的分布范围逐渐向铸锭内部扩展，且热应力的大小也逐渐发生变化。在铸锭的中心区域，热应力相对较小。这是因为中心区域温度变化相对较为均匀，温度梯度较小，热应力的产生也相对较弱。为了更准确地分析热应力的大小和分布规律，提取了铸锭不同位置处的热应力随时间的变化曲线，如图4所示。从曲线中可以看出，在冷却初期，铸锭外层的热应力迅速增大，达到一个较高的值。例如，在铸锭外层某点处，热应力在冷却开始后的短时间内就达到了σ_1MPa，这是由于该点冷却速度快，温度梯度大所致。随着冷却的继续进行，热应力逐渐趋于稳定，但仍保持在一定的水平。而在铸锭内层，热应力的增长速度相对较慢，在冷却初期热应力值较低。随着时间的推移，内层热应力逐渐增大，这是因为内层热量逐渐向外传递，温度梯度逐渐增大。在铸锭中心，热应力始终保持在较低的水平，基本在σ_2MPa以下，这与热应力分布云图的分析结果一致。【此处插入铸锭不同位置热应力随时间变化曲线，横坐标为时间，纵坐标为热应力，曲线需清晰标注各位置，并可适当添加数据点进行说明】热应力对铸锭质量有着重要的影响。当热应力超过铸锭材料的屈服强度时，铸锭会发生塑性变形。在铸锭冷却过程中，如果塑性变形过大，可能会导致铸锭的尺寸精度下降，影响其后续的加工和使用。例如，在某些对尺寸精度要求较高的航空航天零部件制造中，铸锭的塑性变形可能会导致零部件无法满足设计要求，从而造成废品。若热应力进一步超过铸锭材料的抗拉强度，铸锭就会出现裂纹。裂纹的产生会严重降低铸锭的力学性能，使铸锭的强度、韧性等指标大幅下降。裂纹还可能在后续的加工和使用过程中扩展，导致铸锭的断裂失效，严重影响产品的质量和安全性。在大型机械制造中，如发动机曲轴的制造，如果铸锭内部存在裂纹，在发动机高速运转时，曲轴可能会发生断裂，引发严重的安全事故。因此，在电渣重熔空心铸锭过程中，必须充分考虑热应力的影响，采取有效的措施来降低热应力，如优化冷却工艺、调整工艺参数等，以保证铸锭的质量和性能。五、温度场与热应力的影响因素分析5.1工艺参数对温度场和热应力的影响在电渣重熔空心铸锭过程中，工艺参数的变化对温度场和热应力有着显著的影响，深入研究这些影响规律对于优化工艺、提高铸锭质量具有重要意义。电流强度作为关键工艺参数之一，对温度场的影响十分显著。根据焦耳定律Q=I^2Rt，电流强度的改变会直接影响渣池产生的焦耳热，进而改变整个系统的温度分布。当电流强度增大时，渣池产生的焦耳热大幅增加，使得渣池和金属熔池的温度显著升高。通过数值模拟研究发现，在其他条件不变的情况下，将电流强度从I_1增大到I_2（I_2>I_1），渣池的平均温度升高了\DeltaT_1，金属熔池的最高温度升高了\DeltaT_2。这是因为电流强度增大，电流密度增大，渣池电阻不变时，焦耳热与电流强度的平方成正比，导致渣池温度升高，进而通过热传导和对流等方式使金属熔池温度也升高。在实际生产中，若电流强度过大，可能导致渣池和金属熔池温度过高，使得金属液滴的过热程度增加。这不仅会增加金属中元素的挥发损失，例如某些易挥发元素如锌、铅等在高温下挥发加剧，影响铸锭的化学成分；还可能导致铸锭内部组织粗大，出现缩孔、疏松等缺陷。这是因为高温会使金属液的凝固速度减慢，在凝固过程中，由于补缩不充分，容易形成缩孔；同时，高温下晶粒生长速度加快，导致组织粗大。相反，若电流强度过小，渣池产生的热量不足，会使金属熔化速度减慢，生产效率降低，同时可能导致金属液滴在渣池中停留时间过长，增加夹杂物的吸收几率。例如，在生产中发现，当电流强度低于某一临界值时，金属液滴在渣池中长时间停留，会吸附更多的渣中夹杂物，降低铸锭的纯净度。电流强度对热应力也有重要影响。随着电流强度的增大，温度场的不均匀性加剧，进而导致热应力增大。这是因为温度升高会使材料的热膨胀量增加，而铸锭不同部位的温度变化不一致，导致热膨胀的差异增大，从而产生更大的热应力。例如，在铸锭的外层靠近结晶器壁的区域，由于冷却速度快，温度较低，而内部由于焦耳热的作用温度较高。当电流强度增大时，这种温度差异进一步增大，使得外层和内层之间的热应力增大。通过数值模拟计算可知，当电流强度从I_1增大到I_2时，铸锭外层的热应力从\sigma_1增大到\sigma_2，且\sigma_2>\sigma_1。过大的热应力可能导致铸锭发生塑性变形甚至开裂，影响铸锭的质量和性能。电极插入深度同样对温度场和热应力产生重要影响。随着电极插入深度的增加，电流在渣池中的分布路径发生变化，渣池的有效电阻减小。研究表明，当电极插入深度从h_1增加到h_2（h_2>h_1）时，体系的温度呈下降趋势。这是因为电极插入深度增加，电流通过渣池的路径缩短，渣池产生的焦耳热减少，导致体系温度降低。例如，通过数值模拟发现，电极插入深度每增加一定距离，渣池的平均温度会降低\DeltaT_3。电极插入深度还会影响金属熔滴的形成和下落过程。如果电极插入过深，金属熔滴在渣池中下落的距离增加，与渣池的接触时间变长，可能会导致金属熔滴过度精炼，一些有益元素被过多去除；而电极插入过浅，则可能使金属熔滴不能充分与渣池接触，无法有效去除夹杂物和有害元素。在实际生产中，需要根据具体情况选择合适的电极插入深度，以保证电渣重熔过程的稳定性和铸锭质量。在热应力方面，电极插入深度的变化会改变温度场的分布，从而影响热应力的大小和分布。当电极插入深度增加时，由于温度降低，铸锭不同部位的热膨胀差异减小，热应力也会相应减小。例如，在铸锭的某一位置，当电极插入深度从h_1增加到h_2时，热应力从\sigma_3减小到\sigma_4。但如果电极插入深度变化导致温度场分布不均匀性加剧，也可能会使热应力增大。因此，在优化电极插入深度时，需要综合考虑温度场和热应力的变化，以达到最佳的铸锭质量。渣池深度对温度场和热应力也有不可忽视的影响。渣池深度增加，渣池的热容量增大，能够储存更多的热量。在电渣重熔过程中，这会使渣池的温度更加稳定，波动较小。例如，当渣池深度从d_1增加到d_2（d_2>d_1）时，渣池温度的波动范围从\DeltaT_4减小到\DeltaT_5。这是因为较大的热容量使得渣池在吸收和释放热量时，温度变化相对缓慢。渣池深度还会影响热量向金属熔池和铸锭的传递。随着渣池深度的增加，热量传递到金属熔池和铸锭的路径变长，传递效率可能会降低。这可能导致金属熔池的温度分布更加不均匀，进而影响铸锭的凝固过程和质量。例如，在渣池深度较大时，金属熔池底部的温度可能会相对较低，导致凝固速度加快，而金属熔池上部的温度相对较高，凝固速度较慢，从而在铸锭内部形成较大的温度梯度。在热应力方面，渣池深度的变化会影响铸锭的冷却速度和温度梯度，进而影响热应力的大小和分布。当渣池深度增加时，铸锭的冷却速度可能会减慢，这是因为渣池的隔热作用增强，热量散失变慢。冷却速度的减慢会使铸锭内部的温度梯度减小，从而降低热应力。例如，在铸锭的某一部位，当渣池深度从d_1增加到d_2时，热应力从\sigma_5减小到\sigma_6。但如果渣池深度增加导致温度场分布不均匀性加剧，也可能会使热应力增大。因此，在确定渣池深度时，需要综合考虑温度场和热应力的变化，以获得最佳的铸锭质量。5.2材料特性对温度场和热应力的影响材料的热导率、比热容、热膨胀系数等特性在电渣重熔空心铸锭过程中对温度场和热应力有着显著的影响，这些特性的差异会导致铸锭在重熔和冷却过程中呈现出不同的物理行为。热导率是材料传导热量的能力，对温度场分布起着关键作用。热导率较高的材料能够更快速地传导热量，使得温度分布更加均匀。在电渣重熔空心铸锭过程中，若铸锭材料的热导率较大，渣池产生的热量能够迅速传递到铸锭内部，减少温度梯度。例如，对于热导率为k_1的材料，在相同的热传递条件下，其内部温度变化相对较小，温度分布较为平缓。相反，热导率较低的材料会阻碍热量的传递，导致温度梯度增大。如热导率为k_2（k_2<k_1）的材料，热量在其内部传递缓慢，靠近热源的区域温度较高，而远离热源的区域温度较低，从而在材料内部形成较大的温度梯度。通过数值模拟研究发现，当铸锭材料的热导率从k_1降低到k_2时，铸锭内部某位置的温度梯度从\DeltaT_1增大到\DeltaT_2。这种温度梯度的变化会对热应力产生重要影响。根据热应力的产生机制，温度梯度越大，热应力越大。因此，热导率较低的材料在冷却过程中更容易产生较大的热应力，增加铸锭开裂的风险。比热容是单位质量的材料温度升高1K所吸收的热量，它对温度场的变化速率有着重要影响。比热容较大的材料在吸收或释放相同热量时，温度变化相对较小。在电渣重熔空心铸锭的冷却过程中，比热容大的材料能够储存更多的热量，冷却速度相对较慢。例如，材料A的比热容为C_{p1}，材料B的比热容为C_{p2}（C_{p2}>C_{p1}），在相同的冷却条件下，材料B的温度下降速度比材料A慢。这是因为材料B吸收相同的热量时，温度升高的幅度较小，在冷却时释放相同的热量，温度下降的幅度也较小。通过数值模拟计算可知，在冷却时间为t时，材料A的温度从T_1下降到T_2，而材料B的温度从T_1下降到T_3（T_3>T_2）。冷却速度的差异会影响热应力的大小。冷却速度较慢时，铸锭内部各部分的温度变化较为均匀，热应力相对较小；而冷却速度过快，会导致铸锭内部各部分温度变化不一致，热应力增大。因此，比热容较大的材料在冷却过程中产生的热应力相对较小，有利于提高铸锭的质量。热膨胀系数是衡量材料在温度变化时体积或长度变化的参数，它与热应力的产生密切相关。热膨胀系数较大的材料在温度变化时，其体积或长度的变化更为显著。在电渣重熔空心铸锭的冷却过程中，由于铸锭各部分温度不同，热膨胀系数大的材料会产生更大的热膨胀差异。例如，材料C的热膨胀系数为\alpha_1，材料D的热膨胀系数为\alpha_2（\alpha_2>\alpha_1），当温度从T_4降低到T_5时，材料C的长度变化量为\DeltaL_1，材料D的长度变化量为\DeltaL_2（\DeltaL_2>\DeltaL_1）。这种热膨胀差异会在铸锭内部产生较大的热应力。因为铸锭是一个整体，各部分之间相互约束，热膨胀差异会导致内部产生内应力。当热膨胀系数过大时，热应力可能超过材料的屈服强度，导致铸锭发生塑性变形；若热应力进一步超过材料的抗拉强度，铸锭就会出现裂纹。因此，在选择铸锭材料时，需要考虑其热膨胀系数，尽量选择热膨胀系数较小的材料，以降低热应力，提高铸锭的质量。5.3冷却条件对温度场和热应力的影响冷却条件在电渣重熔空心铸锭过程中对温度场和热应力有着重要的影响，这些条件主要包括冷却介质流量、温度以及结晶器结构等方面，它们相互作用，共同决定了铸锭的冷却过程和最终质量。冷却介质流量的变化对温度场和热应力有着显著的影响。当冷却介质流量增加时，其带走热量的能力增强。在电渣重熔空心铸锭过程中，水冷结晶器通过冷却介质将铸锭的热量带走，促使铸锭凝固。例如，当冷却介质流量从Q_1增大到Q_2（Q_2>Q_1）时，铸锭表面与冷却介质之间的热交换更加剧烈，铸锭表面温度迅速降低。通过数值模拟发现，铸锭表面温度在冷却介质流量增大后，在相同时间内下降了\DeltaT_1。由于铸锭表面冷却速度加快，温度梯度增大，使得铸锭内部产生的热应力也相应增大。这是因为铸锭表面和内部的温度变化不一致，表面收缩较快，而内部收缩相对较慢，这种收缩差异导致了热应力的增加。相反，当冷却介质流量减少时，铸锭的冷却速度减慢，温度梯度减小，热应力也随之减小。但冷却介质流量过小，会使铸锭冷却时间延长，生产效率降低，同时可能导致铸锭内部组织粗大，影响铸锭质量。冷却介质温度同样对温度场和热应力产生重要影响。较低的冷却介质温度能够提供更大的温差驱动力，加速热量从铸锭传递到冷却介质中。在实际生产中，将冷却介质温度从T_1降低到T_2（T_2<T_1），铸锭的冷却速度明显加快。数值模拟结果表明，冷却介质温度降低后，铸锭的平均冷却速度提高了v_1。由于冷却速度加快，铸锭内部的温度梯度增大，热应力也会相应增大。过高的热应力可能导致铸锭出现裂纹等缺陷。而冷却介质温度过高，会使铸锭冷却速度过慢，不利于铸锭的凝固过程，可能导致铸锭内部出现缩孔、疏松等缺陷。例如，在某一电渣重熔空心铸锭生产中，当冷却介质温度过高时，铸锭内部出现了明显的缩孔缺陷，严重影响了铸锭的质量。结晶器结构的设计对温度场和热应力的分布也起着关键作用。结晶器的形状、尺寸以及冷却通道的布置等因素都会影响冷却效果和温度分布。不同形状的结晶器，如圆形、方形等，会导致铸锭在冷却过程中的温度分布不同。圆形结晶器在冷却时，热量传递相对均匀，铸锭周向温度差异较小；而方形结晶器在角部区域热量集中，温度梯度较大。通过数值模拟对比圆形和方形结晶器对铸锭温度场的影响发现，方形结晶器角部的温度梯度比圆形结晶器大\DeltaT_2，这会导致方形结晶器角部的热应力增大。结晶器的尺寸也会影响冷却效果。较大尺寸的结晶器，其冷却面积相对较大，能够更有效地带走热量，但同时也可能导致铸锭不同部位的冷却速度差异增大。例如，在大型电渣重熔空心铸锭生产中，采用大尺寸结晶器时，铸锭中心和边缘的冷却速度差异明显，容易产生较大的热应力。冷却通道的布置方式也会影响冷却介质的流动和热量传递。合理的冷却通道布置可以使冷却介质均匀地分布在结晶器周围，提高冷却效率，减小温度梯度和热应力。如采用螺旋式冷却通道，能够使冷却介质在结晶器内形成螺旋状流动，增强热交换效果，使铸锭温度分布更加均匀，热应力减小。六、实验验证与结果对比6.1实验方案设计为了验证数值模拟结果的准确性，设计并开展电渣重熔空心铸锭实验。本次实验的主要目的是通过实际测量电渣重熔空心铸锭过程中的温度和热应力数据，与数值模拟结果进行对比分析，从而评估数值模型的可靠性，并进一步深入了解电渣重熔空心铸锭过程中的物理现象。实验设备选用一台中型电渣重熔炉，其主要组成部分包括水冷结晶器、自耗电极、渣系供给系统、电气控制系统等。水冷结晶器采用铜材质制作，具有良好的导热性能，其内部尺寸为内径D_1、外径D_2、高度H，能够满足本次实验对空心铸锭尺寸的要求。自耗电极选用特定成分的合金钢材料，其化学成分经过精确检测，符合实验要求，电极直径为d，长度为L。渣系供给系统能够精确控制渣料的加入量和加入速度，确保渣池的成分和性能稳定。电气控制系统可实现对电流、电压等电参数的精确调节和监控，本次实验设定电流强度为I，电压为V。实验材料方面，渣料选用以CaF_