基于数值模拟的高温合金立式连铸流动传热与应力特性研究

基于数值模拟的高温合金立式连铸流动传热与应力特性研究一、引言1.1研究背景与意义高温合金，作为一类以铁、镍、钴为基，能在600℃以上高温及一定应力作用下长期工作的金属材料，凭借其优异的高温强度、抗氧化和抗热腐蚀性能、良好的疲劳性能以及断裂韧性等综合性能，在现代工业中占据着举足轻重的地位。因其合金化程度较高，又被称为“超合金”，是广泛应用于航空、航天、石油、化工、舰船等领域的关键材料。在航空航天领域，高温合金是制造航空发动机热端部件的核心材料，用量占发动机总重量的40%-60%。例如，在现代先进航空发动机中，燃烧室、导向叶片、涡轮叶片、涡轮盘等关键热段零部件多由高温合金制成，其性能直接影响发动机的工作效率、可靠性和使用寿命，进而决定飞行器的性能。在能源领域，无论是发电用燃气轮机还是舰船用燃气轮机，高温合金都是燃烧室、过渡导管、导向叶片、涡轮工作叶片以及涡轮盘等关键部件的不可或缺的材料。在核电领域，高温合金是核岛关键热交换器件核电蒸发器的重要基础原料，用于制造燃料元件包壳材料、结构材料和燃料棒定位格架等。在石油化工领域，高温合金用于制造高温高压管道、阀门、热交换器等设备，以承受高温高压和腐蚀性介质的侵蚀。连铸工艺作为一种先进的金属成型工艺，在高温合金的生产中具有重要作用。它将液态金属直接浇铸成预定尺寸的铸坯，省去了传统铸造、锻造等多道工序，显著提高了生产效率。同时，连铸工艺能够使金属被迅速冷却，结晶致密，组织均匀，从而提高产品的机械性能。此外，连续铸造时，铸件上没有浇注系统的冒口，在轧制时不用切头去尾，节约了金属，提高了收得率。连铸工艺还简化了工序，减轻了劳动强度，所需生产面积也大为减少，且易于实现机械化和自动化，铸锭时还能实现连铸连轧，进一步提高生产效率。然而，连铸过程是一个极为复杂的物理过程，涉及到高温钢液的流动、传热、凝固以及应力应变等多个方面。在这一系列冷却、凝固过程中，铸坯的传热尤为关键。冷却不合理会造成铸坯裂纹、脱方、夹杂物等缺陷，严重影响铸坯的质量。例如，若铸坯在冷却过程中各部分温度不均匀，产生较大的热应力，就容易导致裂纹的产生；拉坯速度与冷却速度不匹配，可能使铸坯出现脱方缺陷；钢液中的夹杂物若不能有效上浮去除，会留在铸坯内部，降低铸坯的性能。而且连铸过程伴随着高温，运用人眼法、物理试验法等传统研究方法难以全面、准确地观测和分析其中的物理现象。人眼法依赖操作人员的经验，难以获取标准数据；物理试验法在常温条件下利用介质代替钢液进行模拟，无法真实反映连铸过程中的高温、复杂工况，导致试验数据不准确。数值模拟技术的出现为解决这些问题提供了有效的手段。数值模拟可以对连铸过程中的危险、超越正常条件的、待开发的过程等进行模拟，具有成本低、速度快、能模拟真实条件或理想条件的能力，结果直观、易于理解。通过建立数学模型，对连铸过程中的流动、传热、应力等现象进行数值求解，可以深入了解连铸过程的物理本质，预测铸坯的质量缺陷，为优化连铸工艺参数提供理论依据。例如，通过数值模拟可以研究不同浇铸温度、拉坯速度、冷却强度等工艺参数对铸坯温度分布、凝固进程、应力分布的影响，从而找到最优的工艺参数组合，提高铸坯质量，降低生产成本。因此，对高温合金立式连铸流动传热和应力进行模拟研究具有重要的理论意义和实际应用价值。1.2高温合金概述高温合金，作为现代工业中至关重要的金属材料，是指以铁、镍、钴为基，能在600℃以上的高温及一定应力作用下长期工作的一类金属材料。因其合金化程度较高，也被称为“超合金”。高温合金具有一系列优异的性能，使其在众多领域中发挥着不可或缺的作用。高温合金具有较高的高温强度，在高温环境下能保持稳定的力学性能，有效承受各种应力。其具备良好的抗氧化和抗热腐蚀性能，可在高温氧化环境和燃气腐蚀条件下长期工作。再者，高温合金为单一奥氏体组织，在各种温度下都拥有良好的组织稳定性和使用可靠性，还具有良好的疲劳性能、断裂韧性等综合性能。按照基体元素来划分，高温合金主要分为铁基、镍基、钴基等高温合金。铁基高温合金的使用温度一般只能达到750-780℃，对于在更高温度下使用的耐热部件，则需采用镍基和难熔金属为基的合金。镍基高温合金含镍量在一半以上，适用于1000℃以上的工作条件，在整个高温合金领域占据特殊重要的地位，广泛应用于制造航空喷气发动机、各种工业燃气轮机最热端部件。钴基高温合金以钴为基体，钴含量大约占60%，同时需要加入Cr、Ni等元素来提升耐热性能。虽然钴基高温合金耐热性能较好，但由于各个国家钴资源产量比较少，加工比较困难，因此用量不多，通常用于高温条件（600-1000℃）和较长时间受极限复杂应力的高温零部件，例如航空发动机的工作叶片、涡轮盘、燃烧室热端部件和航天发动机等。若依据合金强化类型来分类，高温合金可分为固溶强化型高温合金和时效沉淀强化合金。根据材料成型方式，又可将高温合金分为铸造高温合金（包括普通铸造合金、单晶合金、定向合金等）、变形高温合金、粉末冶金高温合金（包含普通粉末冶金和氧化物弥散强化高温合金）。铸造高温合金的发展历程丰富而曲折。自20世纪30年代后期起，英、德、美等国就开启了高温合金的研究之旅。第二次世界大战期间，为满足新型航空发动机的需求，高温合金的研究和使用迎来蓬勃发展时期。40年代初，英国率先在80Ni-20Cr合金中加入少量铝和钛，形成γ相以进行强化，成功研制成第一种具有较高高温强度的镍基合金。同一时期，美国为适应活塞式航空发动机用涡轮增压器发展的需要，开始用Vitallium钴基合金制作叶片，还研制出Inconel镍基合金，用于制作喷气发动机的燃烧室。此后，冶金学家为进一步提高合金的高温强度，在镍基合金中加入钨、钼、钴等元素，增加铝、钛含量，研制出一系列牌号的合金，如英国的“Nimonic”，美国的“Mar-M”和“IN”等；在钴基合金中，加入镍、钨等元素，发展出多种高温合金，如X-45、HA-188、FSX-414等。由于钴资源缺乏，钴基高温合金的发展受到一定限制。40年代，铁基高温合金也取得了发展，50年代出现A-286和Incoloy901等牌号，但因高温稳定性较差，从60年代以来发展较为缓慢。苏联于1950年前后开始生产“ЭИ”牌号的镍基高温合金，后来又生产“ЭП”系列变形高温合金和ЖС系列铸造高温合金。70年代美国采用新的生产工艺制造出定向结晶叶片和粉末冶金涡轮盘，研制出单晶叶片等高温合金部件，以适应航空发动机涡轮进口温度不断提高的需求。发展至今，国际市场每年高温金属合金消费量在30万吨左右，广泛应用于航空航天、能源、石油、化工、舰船等各个领域。在国内，自1956年第一炉高温合金GH3030试炼成功以来，我国高温合金的研究、生产和应用已历经多年的发展历程。这一历程可大致分为三个阶段。第一个阶段从1956年至20世纪70年代初，是我国高温合金的创业和起始阶段，主要以仿制前苏联高温合金为主体的合金系列，如GH4033、GH4049、GH2036、GH3030、K401和K403等。第二个阶段从20世纪70年代中至90年代中期，是我国高温合金的提高阶段，主要任务是试制欧美型号的发动机，同时提高高温合金生产工艺技术和产品质量控制水平。第三阶段从20世纪90年代中至今，是我国高温合金的全新发展阶段，主要成果是应用和开发了一批新工艺，研制和生产了一系列高性能、高档次的新合金。目前，我国的高温合金研究主要研究单位有钢铁研究总院、北京航空材料研究院、中国科学院金属研究所、北京科技大学、东北大学、西北工业大学等，主要生产企业包括中航工业、钢研高纳、炼石有色、抚顺特钢、高钢特钢和第二重型机械集团万航模锻厂（二重）等。在此基础上，我国已具备了高温合金新材料、新工艺自主研发和研究的能力。然而，尽管我国高温金属合金材料已发展多年，但行业发展仍处于成长期。由于高温金属合金材料领域具有较高技术含量，该行业企业拥有较深护城河。我国高温金属合金每年需求量在2万吨以上，国内年生产量在1万吨左右，市场容量超过80亿元，其中进口占比较大，在技术水平、生产规模和产品质量等方面与国际先进水平仍存在一定差距，在高端产品领域，如航空航天用高性能单晶高温合金等，对进口的依赖程度较高，一些关键生产设备和技术仍需从国外引进，在生产过程中的质量控制和稳定性方面也有待进一步提高。1.3连铸工艺及研究现状连铸工艺，即连续铸造工艺，是一种先进的金属成型工艺，广泛应用于钢铁、有色金属等材料的生产。其基本原理是将液态金属连续地浇入一个特定的结晶器中，通过结晶器的强制冷却作用，使液态金属在结晶器内初步凝固形成一定形状的坯壳，然后在拉坯力的作用下，将带有液芯的铸坯从结晶器中连续拉出，并在二次冷却区进一步冷却，直至完全凝固，最后切割成所需长度的铸坯。以钢铁连铸为例，具体流程如下：首先，装有精炼好钢水的钢包被运至连铸机上方的大包回转台，回转台转动到浇注位置后，钢包底部的钢水缓慢流入中间包。中间包起到稳定钢流、均匀温度和去除夹杂物的作用，通过塞棒等装置精确控制钢水流量，将钢水按照设定的流量要求注入到结晶器中。结晶器是连铸工艺的核心部件，通常为水冷铜质结构，它使钢水迅速凝固成坯壳。为了避免漏钢，结晶器在振动装置的作用下上下振动，帮助铸坯脱模。带有液芯的坯壳离开结晶器后，进入二次冷却区，通过喷水等方式进行冷却，使其进一步凝固。在拉矫机的作用下，铸坯被连续拉出并矫直，最后由火焰枪或其他切割设备将铸坯切割成一定长度的定尺铸坯。根据后续加工机组的不同，铸轧流程也有所差异。薄板坯连铸连轧工艺，铸坯直接经过隧道保温炉后送到热连轧机进行轧制；常规热连轧机则需要将连铸坯再经过加热炉加热后，送至粗轧机进行粗轧，然后再进行精轧；铸轧全连续工艺不需要将板坯切割成定尺，直接将连铸坯送到粗轧和精轧机进行全连续轧制。在连铸过程中，流动传热和应力对铸坯质量有着至关重要的影响。钢液在结晶器和二次冷却区内的流动状态，直接影响着钢液的温度分布、溶质分布以及夹杂物的运动和上浮。合理的钢液流动可以使温度均匀分布，促进夹杂物上浮去除，减少铸坯内部缺陷。如果钢液流动不合理，可能导致铸坯局部过热或过冷，产生缩孔、疏松、偏析等缺陷。传热过程决定了铸坯的凝固速度和凝固方式，进而影响铸坯的组织结构和性能。若铸坯冷却不均匀，会产生较大的热应力，导致铸坯出现裂纹等缺陷。在连铸过程中，铸坯由于温度变化、凝固收缩以及受到拉坯力、摩擦力等外力作用，会产生复杂的应力应变。过大的应力可能使铸坯产生裂纹、变形等质量问题，严重影响铸坯的质量和后续加工性能。连铸凝固模拟研究旨在通过建立数学模型和数值计算方法，对连铸过程中的流动传热和凝固现象进行模拟分析，从而深入了解连铸过程的物理本质，预测铸坯质量缺陷，为连铸工艺的优化提供理论依据。自20世纪60年代以来，随着计算机技术的飞速发展，连铸凝固模拟研究取得了显著的进展。早期的研究主要集中在对传热过程的模拟，采用简单的数学模型和有限差分法等数值方法，对二维或简单三维的连铸过程进行模拟分析。随着计算技术的不断进步和对连铸过程认识的深入，研究逐渐从单纯的传热模拟扩展到包括流动、传热、凝固和应力应变等多物理场耦合的模拟。在流动模拟方面，采用计算流体力学（CFD）方法，对钢液在结晶器和二次冷却区内的流动进行模拟，研究流动特性及其对铸坯质量的影响。在传热模拟中，考虑了相变潜热、辐射传热、对流换热等多种因素，提高了模拟的准确性。在凝固模拟中，引入了凝固动力学模型，对铸坯的凝固过程进行更精确的描述。在应力应变模拟方面，建立了热弹塑性力学模型，分析铸坯在连铸过程中的应力应变分布，预测裂纹等缺陷的产生。近年来，连铸凝固模拟研究呈现出以下几个发展趋势：一是多物理场耦合模拟的深入发展，将流动、传热、凝固、应力应变以及微观组织演变等多个物理过程进行更加全面和精确的耦合模拟，以更真实地反映连铸过程的复杂性。二是模型精度的不断提高，通过改进数学模型、采用更精细的数值算法和提高计算资源，提高模拟结果的准确性和可靠性。三是与实验研究的紧密结合，通过实验数据验证和修正模拟模型，使模拟结果更具实际应用价值。四是模拟软件的商业化和工程化应用，开发出一系列功能强大、易于操作的连铸凝固模拟软件，广泛应用于钢铁企业的生产实践中，为连铸工艺的优化和铸坯质量的提高提供有力支持。例如，ANSYS、FLUENT等软件在连铸过程模拟中得到了广泛应用，能够对连铸过程中的各种物理现象进行详细的模拟分析。1.4研究内容与方法本研究聚焦于高温合金立式连铸过程，综合运用数值模拟与实验验证的方法，深入探究其流动传热与应力分布规律，旨在为工艺优化提供坚实的理论与实践依据。具体研究内容涵盖以下几个关键方面：建立数学模型：基于高温合金立式连铸的实际物理过程，综合考虑质量守恒、动量守恒、能量守恒等基本定律，建立精确的流动传热和应力分析数学模型。在流动模型中，充分考虑钢液的粘性、湍流特性以及自由表面的复杂变化；传热模型全面涵盖结晶器、二次冷却区的传导、对流和辐射传热，精确处理相变潜热；应力模型则充分考虑铸坯的热膨胀、收缩以及凝固过程中的力学行为，准确描述材料的热弹塑性本构关系。通过这些模型，能够全面、准确地描述连铸过程中各物理量的变化规律。数值模拟分析：运用专业的数值模拟软件ANSYS，对所建立的数学模型进行高效求解。深入研究不同工艺参数，如浇铸温度、拉坯速度、冷却强度等，对高温合金在连铸过程中的流动形态、温度分布、凝固进程以及应力应变分布的显著影响。通过模拟结果的详细分析，揭示各工艺参数与铸坯质量之间的内在联系，为后续的工艺优化提供关键的理论指导。例如，通过模拟不同浇铸温度下铸坯的温度分布和凝固进程，找出最有利于铸坯质量的浇铸温度范围；分析不同拉坯速度对应力应变分布的影响，确定合理的拉坯速度，以减少铸坯裂纹等缺陷的产生。实验验证：设计并开展高温合金立式连铸实验，通过精心布置的热电偶、应变片等测量设备，准确获取连铸过程中铸坯的温度、应力等关键数据。将实验数据与数值模拟结果进行细致对比，深入分析两者之间的差异，进一步验证数值模拟模型的准确性和可靠性。同时，通过实验研究，还可以发现一些在数值模拟中未考虑到的因素，为模型的进一步完善提供实际依据。例如，在实验中观察铸坯表面的质量缺陷，分析其产生原因，与数值模拟结果相结合，找出影响铸坯质量的关键因素，从而对模拟模型进行修正和优化。工艺参数优化：基于数值模拟和实验研究的结果，运用先进的优化算法，对高温合金立式连铸的工艺参数进行全面优化。以铸坯质量最优为目标，以工艺参数为变量，建立优化模型，通过求解该模型，确定最佳的浇铸温度、拉坯速度、冷却强度等工艺参数组合。将优化后的工艺参数应用于实际生产中，显著提高铸坯的质量和生产效率，降低生产成本，实现高温合金立式连铸工艺的高效、优质生产。在研究方法上，本研究采用数值模拟与实验研究相结合的方式。数值模拟方面，选用ANSYS软件，该软件具备强大的多物理场耦合分析能力，在连铸过程模拟中得到广泛应用。通过对模型的合理简化与参数设置，确保模拟结果的准确性与可靠性。实验研究则依据实际连铸生产条件，设计搭建实验平台，严格控制实验过程，保证实验数据的真实性与有效性。对实验数据进行科学的统计分析和误差处理，提高数据的可信度。二、高温合金立式连铸流动传热和应力模拟的理论基础2.1流动传热理论基础在高温合金立式连铸过程中，流动传热现象极为复杂，涉及到多个物理过程的相互作用。为了深入理解和准确模拟这一过程，需要基于流体力学和传热学的基本理论，建立相应的数学模型。流体力学基本方程是描述流体运动的基础，主要包括连续性方程、动量守恒方程和能量守恒方程。连续性方程基于质量守恒定律，反映了流体在流动过程中质量的连续性。对于不可压缩流体，其连续性方程的微分形式为：\frac{\partialu_i}{\partialx_i}=0其中，u_i表示流体在x_i方向上的速度分量（i=1,2,3分别对应x,y,z方向）。该方程表明，在单位时间内，流入和流出某一微元体的流体质量相等，即流体在流动过程中不会出现质量的积累或亏损。动量守恒方程，也称为Navier-Stokes方程，描述了流体在流动过程中的动量变化与所受外力之间的关系。其一般形式为：\rho\frac{\partialu_i}{\partialt}+\rhou_j\frac{\partialu_i}{\partialx_j}=-\frac{\partialp}{\partialx_i}+\mu\frac{\partial^2u_i}{\partialx_j\partialx_j}+\rhog_i其中，\rho为流体密度，t为时间，p为压力，\mu为动力粘度，g_i为重力加速度在x_i方向上的分量。方程左边第一项表示流体的非定常惯性力，第二项表示对流惯性力；右边第一项为压力梯度力，第二项为粘性力，第三项为重力。Navier-Stokes方程全面考虑了流体在流动过程中所受到的各种力的作用，是研究流体流动的核心方程之一。能量守恒方程描述了流体在流动过程中的能量变化，包括内能、动能和势能等。在连铸过程中，主要考虑的是热传导、对流和辐射等传热方式引起的能量变化。其一般形式为：\rhoc_p\frac{\partialT}{\partialt}+\rhoc_pu_j\frac{\partialT}{\partialx_j}=\frac{\partial}{\partialx_j}(\lambda\frac{\partialT}{\partialx_j})+Q其中，c_p为定压比热容，T为温度，\lambda为导热系数，Q为内热源项，如凝固潜热等。方程左边第一项表示流体的非定常热容量变化，第二项表示对流热传递；右边第一项为热传导，第二项为内热源产生的热量。能量守恒方程反映了在连铸过程中，高温合金液在流动过程中的能量转换和传递关系，对于理解铸坯的凝固过程和温度分布具有重要意义。在高温合金立式连铸过程中，钢液的流动通常处于湍流状态。湍流是一种高度复杂的流动状态，其速度、压力等物理量在时间和空间上呈现出不规则的脉动。为了准确模拟湍流流动，需要选择合适的湍流模型。常用的湍流模型包括雷诺平均Navier-Stokes（RANS）模型、大涡模拟（LES）模型和直接数值模拟（DNS）模型等。RANS模型是目前工程中应用最为广泛的湍流模型之一。其基本思想是将湍流瞬时物理量分解为平均量和脉动量，通过求解平均运动方程和湍流模型方程来获得湍流的统计特性。在RANS模型中，常用的有k-\varepsilon模型、k-\omega模型和剪切应力传输（SST）模型等。k-\varepsilon模型基于湍动能k和湍动能耗散率\varepsilon的输运方程来封闭方程组，适用于完全发展的湍流流动，计算效率较高，但在处理近壁区流动和复杂流动时存在一定的局限性。k-\omega模型则以湍动能k和比耗散率\omega为变量，对近壁区流动有较好的模拟效果，但对自由流边界条件较为敏感。SST模型结合了k-\varepsilon模型和k-\omega模型的优点，在近壁区采用k-\omega模型，在远场采用k-\varepsilon模型，能够更好地处理复杂流动问题，如分离流、逆压梯度流等。LES模型通过直接求解大尺度涡的运动方程，而对小尺度涡采用亚网格尺度模型进行模拟。这种方法能够捕捉到湍流的瞬态特性和大尺度结构，模拟结果更为准确，但计算量较大，对计算资源要求较高。LES模型适用于对湍流结构和瞬态特性要求较高的研究，如航空航天领域中飞行器周围的复杂湍流流动模拟等。DNS模型则是直接对Navier-Stokes方程进行数值求解，不引入任何湍流模型，能够精确地模拟湍流的所有尺度结构和瞬态特性。然而，由于DNS模型需要极高的计算分辨率来捕捉湍流的小尺度脉动，计算量随雷诺数的增加呈指数增长，目前仅适用于低雷诺数、简单几何形状的流动模拟，主要用于湍流基础研究，为其他湍流模型的发展提供理论依据。在高温合金立式连铸过程中，凝固是一个关键环节。凝固过程伴随着潜热的释放，这对铸坯的温度分布和凝固进程有着重要影响。因此，在数值模拟中需要合理处理潜热。常用的潜热处理方法包括等价比热容法、温度回升法和热焓法等。等价比热容法是将凝固潜热等效为比热容的增加，通过定义等价比热容（或有效比热容、当量比热容）来处理潜热。假设单位体积、单位时间内固相率的增量为\frac{\partialf_s}{\partialt}，潜热为L，则等价比热容c_{eq}可表示为：c_{eq}=c_p+\frac{L}{\DeltaT}\frac{\partialf_s}{\partialt}其中，c_p为实际比热容，\DeltaT为温度变化范围。在计算中，用等价比热容代替实际比热容代入能量守恒方程进行求解。该方法适用于结晶温度范围较宽的合金，计算相对简单，目前应用较为广泛。温度回升法的基本思路是在凝固开始的一段时间内，将固相增多所释放的潜热折算成所能补偿的温度下降，加到温度计算中去。先不考虑潜热的释放进行温度计算，求出微小时间内的温度降低量\DeltaT_1，如果\DeltaT_1大于某一阈值，表明产生凝固，由于释放潜热，使温度回升到某一值\DeltaT_2（假定没有过冷）。该方法思路明确，实现相对简单，但在程序设计中需要分配大量的存储空间，适应范围相对较窄。热焓法是通过定义热焓来处理潜热。凝固过程中金属的焓H可表示为：H=\int_{T_0}^{T}c_p(T')dT'+f_sL其中，T_0为参考温度。对热焓关于温度求导，得到热焓与温度的关系，再代入能量守恒方程进行求解。这种方法与等价比热容法类似，适用于有一定结晶温度范围的合金，但引入了热焓这一物理量，计算相对复杂。综上所述，流动传热理论基础是高温合金立式连铸流动传热和应力模拟的重要理论依据。通过合理选择流体力学基本方程、湍流模型和潜热处理方法，能够建立准确的数学模型，为深入研究连铸过程中的物理现象和优化连铸工艺提供有力的支持。2.2应力分析理论基础在高温合金立式连铸过程中，铸坯会受到多种力的作用，产生复杂的应力应变。为了准确分析铸坯的应力分布和变形情况，需要基于弹塑性力学的基本理论，建立相应的应力分析模型。弹塑性力学是研究固体材料在弹性和塑性状态下的力学行为的学科，其基本假设包括连续介质假设、小变形假设和均匀性假设等。连续介质假设认为物体内部连续分布着物质，不存在空隙，这使得我们可以用连续函数来描述物体内各点的力学量。小变形假设假定物体在受力后产生的变形远小于物体的原始尺寸，在分析过程中可以忽略因变形引起的几何形状变化，从而简化计算。均匀性假设则表示物体内各点的材料性质相同，不随位置变化而改变，便于建立统一的力学模型。屈服准则是判断材料是否进入塑性状态的重要依据。当材料所受应力满足屈服准则时，材料开始发生塑性变形。常用的屈服准则有Tresca屈服准则和Mises屈服准则。Tresca屈服准则，又称最大剪应力屈服准则，由法国工程师H.Tresca于1864年提出。该准则认为，当材料中的最大剪应力达到某一极限值时，材料开始屈服。其数学表达式为：\tau_{max}=\frac{\sigma_1-\sigma_3}{2}=k其中，\sigma_1、\sigma_3分别为最大主应力和最小主应力，k为材料的屈服常数，可通过实验测定。在平面应力状态下，若\sigma_1\geq\sigma_2\geq\sigma_3，则屈服条件可表示为\sigma_1-\sigma_3=2k；若\sigma_1\geq\sigma_3\geq\sigma_2，则屈服条件为\sigma_1-\sigma_2=2k；若\sigma_2\geq\sigma_1\geq\sigma_3，则屈服条件为\sigma_2-\sigma_3=2k。Tresca屈服准则形式简单，物理意义明确，在分析一些简单受力情况时较为方便，但它没有考虑中间主应力的影响，在某些情况下与实际情况存在一定偏差。Mises屈服准则，由德国力学家R.vonMises于1913年提出。该准则认为，当材料的弹性形变能达到某一极限值时，材料开始屈服。其数学表达式为：\sqrt{\frac{1}{2}[(\sigma_1-\sigma_2)^2+(\sigma_2-\sigma_3)^2+(\sigma_3-\sigma_1)^2]}=\sigma_s其中，\sigma_s为材料的屈服强度。Mises屈服准则考虑了所有主应力的影响，更符合大多数金属材料的屈服行为，在工程中得到了广泛应用。在平面应力状态下，Mises屈服准则的表达式可简化为(\sigma_1-\sigma_2)^2+(\sigma_2-\sigma_3)^2+(\sigma_3-\sigma_1)^2=2\sigma_s^2。本构关系描述了材料在受力时的应力-应变关系，是弹塑性力学中的核心内容。对于弹性材料，其应力-应变关系遵循胡克定律，即应力与应变成线性关系。在各向同性弹性材料中，胡克定律的一般形式为：\sigma_{ij}=\lambda\varepsilon_{kk}\delta_{ij}+2\mu\varepsilon_{ij}其中，\sigma_{ij}为应力张量分量，\varepsilon_{ij}为应变张量分量，\lambda和\mu为拉梅常数，\delta_{ij}为克罗内克符号（当i=j时，\delta_{ij}=1；当i\neqj时，\delta_{ij}=0），\varepsilon_{kk}=\varepsilon_{11}+\varepsilon_{22}+\varepsilon_{33}为体积应变。对于弹塑性材料，其本构关系较为复杂，需要考虑材料的塑性变形特性。常用的弹塑性本构关系模型有增量理论和全量理论。增量理论，如Prandtl-Reuss理论，认为塑性应变增量与应力偏量增量之间存在一定的关系，通过逐步加载的方式来描述材料的弹塑性变形过程。其表达式为：d\varepsilon_{ij}^p=\frac{3}{2}\frac{d\lambda}{s_{ij}}其中，d\varepsilon_{ij}^p为塑性应变增量张量分量，s_{ij}为应力偏量张量分量，d\lambda为塑性乘子，可根据屈服准则和加载条件确定。全量理论，如Hencky-Nadai理论，则试图建立应力与应变全量之间的关系，但由于塑性变形的不可逆性，全量理论存在一定的局限性，一般适用于简单加载情况。在连铸过程中，坯壳主要受到热应力、机械应力和凝固收缩应力等的作用。热应力是由于铸坯在冷却过程中温度分布不均匀，各部分热胀冷缩不一致而产生的应力。例如，铸坯表面冷却速度快，收缩较大，而内部冷却速度慢，收缩较小，从而在表面产生拉应力，在内部产生压应力。机械应力主要包括拉坯力、摩擦力和支撑力等。拉坯力是使铸坯从结晶器中拉出的力，若拉坯力过大，可能导致铸坯产生裂纹；摩擦力是铸坯与结晶器壁、二次冷却区支撑辊等接触时产生的力，会影响铸坯的表面质量和内部应力分布；支撑力是支撑辊对铸坯的作用力，若支撑不均匀，会使铸坯局部受力过大，产生变形和应力集中。凝固收缩应力是由于高温合金在凝固过程中体积收缩受到阻碍而产生的应力。在凝固初期，液态金属的收缩可以通过液体的流动来补充，但随着凝固的进行，固态坯壳逐渐形成，收缩受到限制，就会产生收缩应力。热应力的计算通常基于热弹性理论，通过求解温度场和应力场的耦合方程来得到。首先根据传热学原理计算铸坯的温度分布，然后根据热膨胀系数和弹性模量等材料参数，利用热弹性力学的基本方程计算热应力。例如，对于一维热传导问题，温度分布可由热传导方程\rhoc_p\frac{\partialT}{\partialt}=\frac{\partial}{\partialx}(\lambda\frac{\partialT}{\partialx})求解得到，热应力则可通过公式\sigma=E\alpha\DeltaT计算，其中E为弹性模量，\alpha为热膨胀系数，\DeltaT为温度变化量。机械应力的计算较为复杂，需要考虑铸坯的受力状态和边界条件。对于拉坯力，可以根据连铸机的工艺参数和铸坯的质量等因素进行估算；摩擦力可通过摩擦系数和接触压力等参数计算；支撑力则需要根据支撑辊的布置和铸坯的变形情况来确定。在实际计算中，通常采用有限元等数值方法，将铸坯离散为多个单元，对每个单元进行受力分析，然后通过组装得到整个铸坯的应力分布。综上所述，应力分析理论基础是高温合金立式连铸应力模拟的重要理论依据。通过合理选择屈服准则、本构关系和应力计算方法，能够建立准确的应力分析模型，为深入研究连铸过程中的应力分布和变形情况，预测铸坯的质量缺陷提供有力的支持。2.3数值模拟方法数值模拟是本研究中深入探究高温合金立式连铸流动传热和应力分布规律的关键手段。在数值模拟过程中，有限元方法（FEM）作为一种高效且广泛应用的数值计算方法，发挥着核心作用。有限元方法的基本原理是将连续的求解域离散为一组有限个、且按一定方式相互连接在一起的单元的组合体。通过在每一个单元内假设近似函数，来分片地表示求解域上待求的未知场函数。以二维传热问题为例，对于一个具有复杂几何形状的高温合金铸坯，首先将其划分成众多三角形或四边形单元。在每个单元内，假设温度场可以用一个简单的线性函数来近似表示，如对于三角形单元，温度T可表示为T=a+bx+cy，其中a、b、c为待定系数，可通过单元节点的温度值来确定。这样，将整个求解域上复杂的温度分布问题转化为对每个单元内简单函数的求解。有限元方法的基本流程主要包括网格划分、单元选择、数值求解等关键步骤。在网格划分环节，根据高温合金立式连铸模型的几何形状和物理过程特点，选择合适的网格类型和大小进行划分。对于形状规则的区域，如结晶器的主体部分，可采用结构化网格，这种网格具有排列整齐、计算效率高的优点；对于形状复杂的区域，如铸坯的拐角处或与冷却介质接触的边界区域，采用非结构化网格，它能更好地适应复杂的几何形状，提高网格划分的质量和计算精度。网格的大小对计算结果的精度和计算效率有着显著影响。较小的网格尺寸能够更精确地描述物理量的变化，但会增加计算量和计算时间；较大的网格尺寸虽然计算速度快，但可能会导致计算精度下降。因此，需要通过网格无关性验证来确定合适的网格尺寸。例如，进行不同网格密度下的模拟计算，对比计算结果，当网格密度增加到一定程度后，计算结果的变化小于某个设定的阈值，此时的网格尺寸即可认为是合适的，既能保证计算精度，又能控制计算成本。在单元选择方面，根据具体的物理问题和计算需求，选用合适的单元类型。对于流动传热问题，常用的单元类型有四边形单元和三角形单元。四边形单元具有较高的计算精度和计算效率，适用于大部分规则区域的模拟；三角形单元则更灵活，能够适应复杂的几何形状，在处理边界和局部复杂区域时具有优势。对于应力分析问题，可选用具有更高阶位移模式的单元，如六面体单元或四面体单元，以更准确地描述应力和应变的分布。不同单元类型在计算精度和计算效率上存在差异。例如，高阶单元虽然能够提供更高的计算精度，但计算过程更为复杂，计算量也更大；低阶单元计算简单、效率高，但在描述复杂物理场时精度可能相对较低。因此，在实际应用中，需要综合考虑问题的复杂程度、计算精度要求和计算资源等因素，选择合适的单元类型。数值求解是有限元方法的核心步骤。在完成网格划分和单元选择后，将控制方程离散化，得到一组代数方程组，然后通过迭代求解的方式得到数值解。在迭代过程中，需要设置合理的收敛准则，以确保计算结果的准确性和可靠性。收敛准则通常基于残差的大小来确定，当残差小于某个预设的阈值时，认为计算结果收敛，迭代过程结束。常见的求解器有直接求解器和迭代求解器。直接求解器通过直接对系数矩阵进行分解和求解，能够得到精确的解，但对于大规模问题，其计算量和存储量较大；迭代求解器则通过不断迭代逼近精确解，计算量相对较小，适用于大规模问题的求解。例如，在求解大型高温合金连铸模型的流动传热问题时，由于涉及的网格数量众多，采用迭代求解器如共轭梯度法或广义极小残差法（GMRES），可以在合理的计算时间内得到满足精度要求的解。在数值模拟中，边界条件和初始条件的处理方法至关重要，它们直接影响着模拟结果的准确性。边界条件主要包括速度边界条件、温度边界条件和应力边界条件等。对于速度边界条件，在浇口处，根据浇铸工艺设定钢液的入口速度；在铸坯与结晶器壁接触的边界，考虑到铸坯与结晶器壁之间的摩擦力，可采用无滑移边界条件，即铸坯表面速度与结晶器壁速度相同。对于温度边界条件，在结晶器壁，通过冷却介质的流量和温度，利用对流换热公式计算结晶器壁与铸坯之间的换热系数，从而确定边界的热流密度；在二次冷却区，根据冷却水流速、温度和喷淋方式等参数，确定冷却区边界的对流换热系数和热流密度。在应力边界条件方面，在铸坯的底部，考虑拉坯力的作用，设定相应的拉力边界条件；在铸坯与支撑辊接触的位置，根据支撑辊的布置和压力分布，设定接触压力边界条件。初始条件是模拟开始时各物理量的分布状态。对于温度场，将浇铸温度作为初始温度，即整个铸坯在初始时刻的温度均为浇铸温度。对于速度场，初始时刻钢液在结晶器内的速度分布可根据浇铸工艺和初始状态假设进行设定，如假设钢液在结晶器内的初始速度为均匀分布或根据经验公式计算得到的分布。对于应力场，由于在初始时刻铸坯尚未受到明显的外力作用，可假设初始应力为零。准确设定边界条件和初始条件是保证数值模拟结果准确性的关键。如果边界条件设定不合理，如在温度边界条件中，换热系数的取值不准确，会导致铸坯温度场的模拟结果与实际情况偏差较大；初始条件设定不当，如初始速度场与实际浇铸情况不符，会影响钢液的流动形态和传热过程的模拟结果。综上所述，数值模拟方法中的有限元方法通过合理的网格划分、单元选择和数值求解，以及准确的边界条件和初始条件处理，能够有效地对高温合金立式连铸流动传热和应力进行模拟分析，为深入研究连铸过程的物理本质和优化连铸工艺提供有力的支持。三、高温合金立式连铸流动传热模拟分析3.1模型建立与参数设置为了深入研究高温合金立式连铸过程中的流动传热现象，首先需要建立精确的几何模型。考虑到实际连铸过程的复杂性，对模型进行了合理的简化。假设结晶器为二维矩形结构，铸坯在结晶器内沿竖直方向向下运动，忽略结晶器的圆角和其他次要结构，以简化计算过程并突出主要物理现象。结晶器的尺寸为长L=100\text{mm}，宽W=50\text{mm}，高度H=300\text{mm}。铸坯初始尺寸与结晶器内部尺寸相同，在连铸过程中，铸坯逐渐凝固并向下拉出结晶器。这样的尺寸设定是基于实际连铸生产中的常见规格，同时也考虑到计算资源和模拟精度的平衡。通过合理选择尺寸，可以在保证模拟结果准确性的前提下，减少计算量，提高模拟效率。在材料参数方面，选用常见的镍基高温合金作为研究对象。其密度\rho=8.2\text{g/cm}^3，导热系数\lambda=11\text{W/(m·K)}（在800℃时），定压比热容c_p=0.42\text{J/(g·K)}（在1000℃时）。这些参数是通过查阅相关文献资料以及实验测量得到的，它们反映了镍基高温合金在高温环境下的热物理性质。密度决定了合金的质量分布，导热系数影响热量的传导速度，定压比热容则体现了合金吸收和释放热量的能力，这些参数对于准确模拟连铸过程中的流动传热现象至关重要。冷却介质采用水，其密度\rho_w=1000\text{kg/m}^3，导热系数\lambda_w=0.6\text{W/(m·K)}，定压比热容c_{p,w}=4.2\text{kJ/(kg·K)}。水作为一种常见且高效的冷却介质，在连铸过程中起着关键的冷却作用。其热物理参数决定了冷却效果和热量传递的效率，准确设定这些参数对于模拟结晶器内的冷却过程和铸坯的温度分布具有重要意义。边界条件和初始条件的设置直接影响模拟结果的准确性。在温度边界条件方面，结晶器壁与冷却介质之间通过对流换热进行热量传递，对流换热系数h=5000\text{W/(m}^2\text{·K)}，这一数值是根据实际连铸生产中的经验数据和相关研究确定的。在结晶器入口处，高温合金液的温度设定为浇铸温度T_0=1550\text{℃}，这是高温合金在液态下进入结晶器的初始温度，对后续的凝固和传热过程有着重要影响。对于速度边界条件，在结晶器入口处，高温合金液以一定的速度流入，入口速度v_0=0.05\text{m/s}，该速度是根据连铸工艺的要求和实际生产经验设定的，它决定了高温合金液在结晶器内的流动形态和填充速度。在铸坯与结晶器壁接触的边界，假设铸坯与结晶器壁之间无相对滑动，即铸坯表面速度与结晶器壁速度相同，这一假设符合实际连铸过程中的物理现象，有助于准确模拟铸坯在结晶器内的运动和传热过程。在热流密度边界条件方面，考虑到结晶器壁与高温合金液之间的热交换，以及铸坯在凝固过程中的潜热释放，通过能量守恒原理确定热流密度。在结晶器壁，根据对流换热公式q=h(T-T_w)计算热流密度，其中T为铸坯表面温度，T_w为冷却介质温度。在铸坯内部，考虑凝固潜热的释放，采用等价比热容法处理潜热，将凝固潜热等效为比热容的增加，通过定义等价比热容来处理潜热，从而准确模拟铸坯在凝固过程中的热量变化和温度分布。初始条件设定为：在模拟开始时刻，整个计算区域内高温合金液的温度均为浇铸温度T_0=1550\text{℃}，速度为零。这一初始条件反映了连铸过程开始时的状态，为后续的模拟计算提供了基础。通过合理设置边界条件和初始条件，可以使模拟结果更接近实际连铸过程，为深入研究连铸过程中的流动传热现象提供可靠的依据。3.2流动特性分析在高温合金立式连铸过程中，钢液的流动特性对铸坯质量有着至关重要的影响。通过数值模拟，得到了常规工况下铸坯的速度矢量图和流线图，以此分析其流场分布规律。从速度矢量图（图1）可以清晰地看到，在结晶器入口处，高温合金液以设定的入口速度v_0=0.05\text{m/s}流入，速度方向垂直向下。随着钢液在结晶器内的流动，由于受到结晶器壁的约束和冷却作用，靠近结晶器壁的钢液速度逐渐减小，形成了明显的速度梯度。在结晶器的中心区域，钢液速度相对较大，且分布较为均匀，呈现出较为稳定的流动状态。这种速度分布差异导致钢液在结晶器内的流动呈现出不均匀的特性，对铸坯的凝固过程和质量产生影响。[此处插入速度矢量图，图1：常规工况下铸坯的速度矢量图]流线图（图2）进一步展示了钢液在结晶器内的流动路径。可以观察到，钢液从结晶器入口流入后，沿着特定的流线向下流动。在结晶器的上部区域，流线较为密集，表明钢液流动速度较快；随着钢液向下流动，流线逐渐稀疏，速度也相应减小。在结晶器的底部，由于受到拉坯力和铸坯凝固的影响，钢液的流动方向发生了一定的改变，流线呈现出弯曲的形态。这种流动路径的变化反映了钢液在连铸过程中的复杂流动特性，对铸坯的温度分布和凝固进程有着重要的影响。[此处插入流线图，图2：常规工况下铸坯的流线图]为了深入探究不同工艺参数对流场的影响，分别研究了拉速和浇注温度的变化对流场的作用。当拉速从0.05\text{m/s}增加到0.1\text{m/s}时，钢液在结晶器内的整体流速明显增大。从速度矢量图（图3）可以看出，入口处钢液的速度增加，且在结晶器内的速度分布更加不均匀，靠近结晶器壁的速度梯度进一步增大。这是因为拉速的提高使得钢液在结晶器内的停留时间缩短，需要更快地向下流动以满足拉坯的需求。这种流速的变化会导致钢液与结晶器壁之间的摩擦力增大，可能影响铸坯的表面质量，同时也会改变钢液的温度分布和凝固进程，增加铸坯内部出现缺陷的风险。[此处插入拉速为0.1m/s时的速度矢量图，图3：拉速为0.1m/s时铸坯的速度矢量图]当浇注温度从1550\text{℃}升高到1600\text{℃}时，钢液的粘度降低，流动性增强。在速度矢量图（图4）中可以观察到，结晶器内钢液的速度略有增加，且速度分布更加均匀。这是因为温度升高使得钢液的分子热运动加剧，内摩擦力减小，从而更容易流动。然而，浇注温度过高也可能带来一些问题，如铸坯凝固时间延长，容易导致铸坯内部出现缩孔、疏松等缺陷；同时，过高的温度还会增加结晶器的热负荷，影响结晶器的使用寿命。[此处插入浇注温度为1600℃时的速度矢量图，图4：浇注温度为1600℃时铸坯的速度矢量图]综上所述，常规工况下铸坯的流场分布呈现出明显的不均匀性，不同工艺参数如拉速和浇注温度的变化会显著影响流场特性。在实际连铸生产中，需要合理控制这些工艺参数，以优化钢液的流动状态，提高铸坯质量。3.3传热特性分析在高温合金立式连铸过程中，传热特性对铸坯的凝固进程和质量有着关键影响。通过数值模拟，深入分析了常规工况下铸坯的温度场分布规律，得到了铸坯在不同时刻的等温线图（图5）。从图中可以看出，在连铸初期，高温合金液刚进入结晶器时，铸坯整体温度较高，等温线较为稀疏且分布相对均匀，表明此时温度梯度较小。随着连铸过程的进行，铸坯与结晶器壁接触，热量迅速通过结晶器壁传递给冷却介质，铸坯表面温度急剧下降，等温线在靠近结晶器壁的区域变得密集，形成了较大的温度梯度。在铸坯的中心区域，由于热量传递相对较慢，温度仍然较高，等温线分布较为稀疏。这一温度分布特征导致铸坯在凝固过程中，表面先凝固形成坯壳，而中心区域仍为液态，随着时间的推移，坯壳逐渐增厚，直至中心区域完全凝固。[此处插入不同时刻的等温线图，图5：不同时刻铸坯的等温线图（a）t=10s（b）t=30s（c）t=50s]为了更直观地了解铸坯温度随时间的变化情况，绘制了铸坯不同位置处的温度随时间变化曲线（图6）。选取铸坯表面、距表面10mm处和中心位置三个代表性位置进行分析。从曲线可以明显看出，铸坯表面温度下降最快，在短时间内就从浇铸温度1550℃迅速降至较低温度，这是因为表面直接与结晶器壁接触，散热条件良好。随着时间的增加，表面温度逐渐趋于稳定，但仍低于铸坯内部温度。距表面10mm处的温度下降速度次之，其温度变化曲线介于表面和中心位置之间。铸坯中心位置的温度下降最慢，在较长时间内仍保持较高温度，这是由于中心区域热量传递到表面需要经过较长的路径，且受到周围已凝固坯壳的阻碍。随着连铸过程的持续，中心位置的温度也逐渐降低，最终达到凝固温度。[此处插入温度随时间变化曲线，图6：铸坯不同位置处的温度随时间变化曲线]进一步研究不同工艺参数对温度场的影响，发现冷却水量和冷却水温的变化对铸坯温度场有着显著作用。当冷却水量从10L/min增加到20L/min时，从铸坯的等温线图（图7）可以观察到，靠近结晶器壁的区域等温线更加密集，温度梯度进一步增大。这是因为冷却水量的增加，使得结晶器壁与冷却介质之间的换热能力增强，铸坯表面的热量能够更快地被带走，从而导致表面温度下降更快，温度梯度增大。这种温度梯度的变化会影响铸坯的凝固速度和凝固方式，可能导致铸坯内部产生更大的热应力，增加裂纹产生的风险。[此处插入冷却水量为20L/min时的等温线图，图7：冷却水量为20L/min时铸坯的等温线图]当冷却水温从20℃降低到10℃时，铸坯的整体温度明显下降，等温线向低温区域移动（图8）。这是因为较低的冷却水温提供了更大的温度差，增强了铸坯与冷却介质之间的换热驱动力，使得铸坯的散热速度加快，整体温度降低。然而，冷却水温过低也可能带来一些问题，如铸坯凝固速度过快，导致内部气体和夹杂物来不及排出，从而影响铸坯质量；同时，过快的凝固速度还可能使铸坯产生较大的收缩应力，增加裂纹产生的可能性。[此处插入冷却水温为10℃时的等温线图，图8：冷却水温为10℃时铸坯的等温线图]综上所述，常规工况下铸坯的温度场分布呈现出明显的不均匀性，不同工艺参数如冷却水量和冷却水温的变化会显著影响温度场特性。在实际连铸生产中，需要合理控制这些工艺参数，以优化铸坯的温度分布，确保铸坯质量。3.4坯壳生长分析在高温合金立式连铸过程中，坯壳生长是一个至关重要的环节，它直接影响铸坯的质量和性能。通过数值模拟，深入分析坯壳生长的过程和规律，探讨不同工艺参数对坯壳生长的影响，对于优化连铸工艺具有重要意义。坯壳生长的过程是一个动态的凝固过程。在连铸初期，高温合金液刚进入结晶器，由于结晶器壁的强制冷却作用，靠近结晶器壁的高温合金液迅速散热，温度降低，开始凝固形成坯壳。随着时间的推移，坯壳逐渐增厚，液态金属不断向固态转变。在这个过程中，热量通过坯壳向结晶器壁传递，再由结晶器壁传递给冷却介质。同时，铸坯在拉坯力的作用下不断向下运动，新的高温合金液不断补充到结晶器中，继续凝固形成坯壳。从坯壳厚度随时间变化曲线（图9）可以清晰地看出，坯壳厚度随时间呈近似线性增加的趋势。在开始阶段，坯壳生长速度较快，这是因为此时高温合金液与结晶器壁之间的温差较大，热量传递速度快，凝固速度也相应较快。随着时间的推移，坯壳厚度逐渐增加，坯壳内部的温度梯度逐渐减小，热量传递速度变慢，坯壳生长速度也逐渐减缓。[此处插入坯壳厚度随时间变化曲线，图9：坯壳厚度随时间变化曲线]坯壳厚度在铸坯横截面上的分布也呈现出一定的规律。从坯壳厚度沿铸坯宽度方向的分布曲线（图10）可以看出，在铸坯的边缘处，坯壳厚度较大，这是因为边缘处与结晶器壁的接触面积大，散热条件好，凝固速度快。而在铸坯的中心区域，坯壳厚度较小，这是由于中心区域热量传递到结晶器壁需要经过较长的路径，且受到周围已凝固坯壳的阻碍，散热速度慢，凝固速度也相应较慢。这种坯壳厚度的不均匀分布会导致铸坯在凝固过程中产生应力，若应力过大，可能会导致铸坯出现裂纹等缺陷。[此处插入坯壳厚度沿铸坯宽度方向的分布曲线，图10：坯壳厚度沿铸坯宽度方向的分布曲线]不同工艺参数对坯壳生长有着显著的影响。结晶器锥度是影响坯壳生长的重要参数之一。结晶器锥度是指结晶器从上口到下口的尺寸变化率，它对铸坯与结晶器壁之间的间隙和传热效果有着重要影响。当结晶器锥度较小时，铸坯与结晶器壁之间的间隙较大，传热效果变差，坯壳生长速度减慢，可能导致铸坯表面质量下降，甚至出现漏钢等问题。相反，当结晶器锥度较大时，铸坯与结晶器壁之间的间隙过小，摩擦力增大，可能会引起铸坯表面划伤，同时也会影响坯壳的均匀生长。通过模拟不同结晶器锥度下的坯壳生长情况（图11），发现当结晶器锥度为0.8%时，坯壳生长较为均匀，铸坯表面质量较好。这是因为在这个锥度下，铸坯与结晶器壁之间的间隙适中，既能保证良好的传热效果，又能避免过大的摩擦力。[此处插入不同结晶器锥度下坯壳生长情况对比图，图11：不同结晶器锥度下坯壳生长情况对比图（a）锥度0.5%（b）锥度0.8%（c）锥度1.0%]保护渣性能对坯壳生长也有着重要作用。保护渣在连铸过程中起着润滑、保温和吸收夹杂物等作用。保护渣的熔点、粘度和传热性能等参数会影响其在结晶器内的行为，进而影响坯壳生长。当保护渣熔点过高时，在结晶器内不易熔化，无法形成良好的液渣层，导致铸坯与结晶器壁之间的润滑效果变差，摩擦力增大，影响坯壳生长，甚至可能导致铸坯表面出现裂纹。若保护渣粘度过大，液渣层的流动性变差，难以均匀地分布在铸坯与结晶器壁之间，同样会影响润滑效果和传热效果，使坯壳生长不均匀。通过模拟不同保护渣性能下的坯壳生长情况（图12），发现当保护渣熔点为1100℃，粘度为0.2Pa・s时，保护渣能够在结晶器内良好地熔化和铺展，形成均匀的液渣层，有效地降低铸坯与结晶器壁之间的摩擦力，促进坯壳均匀生长，提高铸坯表面质量。[此处插入不同保护渣性能下坯壳生长情况对比图，图12：不同保护渣性能下坯壳生长情况对比图（a）熔点1150℃，粘度0.3Pa・s（b）熔点1100℃，粘度0.2Pa・s（c）熔点1050℃，粘度0.1Pa・s]综上所述，坯壳生长在高温合金立式连铸过程中呈现出特定的过程和规律，坯壳厚度随时间和位置的变化具有一定的特征。不同工艺参数如结晶器锥度和保护渣性能对坯壳生长有着显著影响，在实际连铸生产中，需要合理控制这些工艺参数，以促进坯壳均匀生长，提高铸坯质量。四、高温合金立式连铸应力模拟分析4.1应力模型建立与参数设置在完成高温合金立式连铸流动传热模拟的基础上，进一步添加应力分析模块，构建完整的应力分析模型，以深入研究连铸过程中铸坯的应力分布和变形情况。为实现这一目标，将应力分析模块与已有的流动传热模型进行紧密耦合。在耦合过程中，充分考虑流动传热对铸坯应力应变的影响，以及应力应变反过来对流动传热的作用。例如，铸坯在凝固过程中的收缩变形会改变其与结晶器壁之间的接触状态，进而影响传热系数和热流密度；而温度分布的不均匀会导致热应力的产生，影响铸坯的应力应变分布。通过这种双向耦合的方式，能够更真实地反映连铸过程中各物理量之间的相互作用关系。为了准确模拟铸坯的应力应变，需要确定高温合金在不同温度和应变率下的力学性能参数。查阅相关文献资料以及实验测量数据，获取了高温合金在不同温度下的弹性模量、屈服强度等关键力学性能参数。在1000℃时，弹性模量E=150\text{GPa}，屈服强度\sigma_y=350\text{MPa}；在1200℃时，弹性模量E=120\text{GPa}，屈服强度\sigma_y=280\text{MPa}。这些参数会随着温度和应变率的变化而发生改变，在模拟过程中，采用合适的本构模型来描述这种变化关系。选用热弹塑性本构模型，该模型能够考虑材料在高温下的弹性变形、塑性变形以及温度对材料性能的影响。其数学表达式为：\sigma_{ij}=D_{ijkl}(\varepsilon_{kl}-\varepsilon_{kl}^p-\alpha\DeltaT\delta_{ij})其中，\sigma_{ij}为应力张量分量，D_{ijkl}为弹性刚度张量，\varepsilon_{kl}为总应变张量分量，\varepsilon_{kl}^p为塑性应变张量分量，\alpha为热膨胀系数，\DeltaT为温度变化量，\delta_{ij}为克罗内克符号。通过该本构模型，能够准确地描述高温合金在连铸过程中的应力应变行为。在连铸过程中，铸坯受到多种力的作用，包括热应力、机械应力和凝固收缩应力等。为了准确模拟这些应力，需要确定合理的应力边界条件。在铸坯的底部，考虑拉坯力的作用，设定拉力边界条件，拉坯力F=5000\text{N}，这一数值是根据连铸机的实际拉坯能力和铸坯的质量等因素确定的。在铸坯与结晶器壁接触的边界，考虑摩擦力的作用，采用库仑摩擦定律来描述摩擦力，摩擦系数\mu=0.1，该数值是根据铸坯与结晶器壁的材料特性和表面状态等因素确定的。在铸坯与支撑辊接触的位置，根据支撑辊的布置和压力分布，设定接触压力边界条件，通过实验测量和理论分析确定接触压力的大小和分布情况。综上所述，通过在流动传热模型基础上添加应力分析模块，并合理设置模型耦合方式、力学性能参数和应力边界条件，建立了准确的高温合金立式连铸应力分析模型，为后续的应力模拟分析奠定了坚实的基础。4.2热应力分析在高温合金立式连铸过程中，热应力是影响铸坯质量的关键因素之一。热应力的产生主要源于铸坯在冷却过程中温度分布的不均匀性，导致各部分热胀冷缩程度不一致。通过数值模拟，对铸坯凝固过程中的热应力分布规律进行深入分析，结果如图13所示。[此处插入铸坯不同时刻的热应力分布云图，图13：铸坯不同时刻的热应力分布云图（a）t=10s（b）t=30s（c）t=50s]从热应力分布云图可以清晰地看出，在连铸初期（t=10s），铸坯整体温度较高，温度梯度相对较小，热应力分布较为均匀，数值也相对较低。随着连铸过程的推进（t=30s），铸坯表面与结晶器壁接触，热量迅速散失，温度急剧下降，而内部温度仍然较高，这使得表面与内部之间产生较大的温度梯度，从而在铸坯表面产生较大的拉应力，在内部产生压应力。在铸坯的拐角处，由于散热条件的差异，温度梯度更为明显，热应力集中现象较为突出，热应力数值明显高于其他部位。到了连铸后期（t=50s），铸坯大部分已经凝固，温度梯度逐渐减小，但在铸坯的中心区域和表面附近，仍然存在一定的温度差，导致热应力仍然存在。此时，热应力分布呈现出从表面向中心逐渐减小的趋势，表面的拉应力仍然较大，若超过铸坯材料的抗拉强度，就可能导致表面裂纹的产生。为了更直观地了解热应力随时间和位置的变化情况，绘制了铸坯表面和中心位置的热应力随时间变化曲线（图14）。从曲线可以看出，铸坯表面的热应力在开始阶段迅速增大，这是因为表面温度下降速度快，与内部形成较大的温度梯度。随着时间的推移，表面温度逐渐趋于稳定，温度梯度减小，热应力也逐渐减小，但始终保持在一定的数值范围内。铸坯中心位置的热应力变化相对较为平缓，在连铸初期，由于温度较高，热应力较小。随着凝固的进行，中心温度逐渐降低，与表面的温度差增大，热应力逐渐增大，在凝固后期，随着温度梯度的减小，热应力也逐渐减小。[此处插入铸坯表面和中心位置的热应力随时间变化曲线，图14：铸坯表面和中心位置的热应力随时间变化曲线]热应力对铸坯质量有着显著的影响。当热应力超过铸坯材料的屈服强度时，铸坯会发生塑性变形；若热应力继续增大，超过材料的抗拉强度，就会导致铸坯产生裂纹。裂纹一旦产生，不仅会降低铸坯的力学性能，还可能在后续加工过程中进一步扩展，影响最终产品的质量和使用寿命。热应力还可能导致铸坯出现变形、翘曲等缺陷，影响铸坯的尺寸精度和表面质量。在实际连铸生产中，需要采取有效的措施来控制热应力，如优化冷却制度，使铸坯各部分冷却均匀；合理调整工艺参数，如拉坯速度、浇铸温度等，以减小温度梯度；对铸坯进行适当的热处理，消除残余热应力。综上所述，铸坯凝固过程中的热应力分布呈现出明显的规律性，随时间和位置的变化而变化。热应力的产生对铸坯质量有着重要影响，在高温合金立式连铸过程中，必须充分考虑热应力的因素，采取相应的措施来降低热应力，提高铸坯质量。4.3机械应力分析在高温合金立式连铸过程中，机械应力是影响铸坯质量的重要因素之一。铸坯在拉坯、矫直等过程中会受到多种机械力的作用，这些力会导致铸坯内部产生复杂的机械应力分布。通过数值模拟，深入分析铸坯在这些过程中的机械应力分布规律，对于优化连铸工艺、提高铸坯质量具有重要意义。在拉坯过程中，铸坯受到拉坯力的作用，拉坯力通过拉矫机传递到铸坯上，使铸坯克服结晶器壁的摩擦力以及自身的凝固收缩阻力，从结晶器中连续拉出。从拉坯过程中铸坯的机械应力分布云图（图15）可以看出，在铸坯的底部，由于直接受到拉坯力的作用，机械应力较大，呈现出明显的拉应力状态。随着离底部距离的增加，拉坯力逐渐传递和分散，机械应力逐渐减小。在铸坯的中心区域，机械应力相对较小，分布较为均匀。在铸坯与结晶器壁接触的部位，由于摩擦力的作用，会产生一定的剪切应力，该剪切应力在接触面上呈现出不均匀分布，靠近拉坯方向的一侧剪切应力相对较大。[此处插入拉坯过程中铸坯的机械应力分布云图，图15：拉坯过程中铸坯的机械应力分布云图]为了更直观地了解机械应力随时间和位置的变化情况，绘制了铸坯底部和中心位置的机械应力随时间变化曲线（图16）。从曲线可以看出，铸坯底部的机械应力在拉坯开始时迅速增大，达到一定值后保持相对稳定。这是因为在拉坯初期，铸坯需要克服较大的初始阻力，随着拉坯过程的进行，阻力逐渐稳定，拉坯力也相应保持稳定。铸坯中心位置的机械应力在拉坯过程中变化相对较小，始终保持在较低的水平。这是由于中心位置受到拉坯力的影响相对较小，主要受到周围坯壳的约束作用。[此处插入铸坯底部和中心位置的机械应力随时间变化曲线，图16：铸坯底部和中心位置的机械应力随时间变化曲线]在矫直过程中，铸坯需要从弯曲状态转变为直线状态，这会导致铸坯内部产生弯曲应力和剪切应力。从矫直过程中铸坯的机械应力分布云图（图17）可以观察到，在铸坯的弯曲部位，弯曲应力较大，呈现出明显的拉应力和压应力分布。在铸坯的外侧，受到拉伸作用，产生拉应力；在铸坯的内侧，受到压缩作用，产生压应力。随着矫直过程的进行，弯曲应力逐渐减小，铸坯逐渐趋于直线状态。在矫直过程中，铸坯还会受到剪切应力的作用，剪切应力主要分布在铸坯的中性层附近，其大小和分布与矫直方式和矫直参数密切相关。[此处插入矫直过程中铸坯的机械应力分布云图，图17：矫直过程中铸坯的机械应力分布云图]机械应力产生的原因主要包括拉坯力、摩擦力、矫直力等外力的作用，以及铸坯自身的凝固收缩和热膨胀等因素。拉坯力是为了克服铸坯与结晶器壁之间的摩擦力、铸坯的凝固收缩阻力以及重力等，将铸坯从结晶器中拉出。摩擦力则是由于铸坯与结晶器壁、支撑辊等接触表面之间的相对运动而产生的。矫直力是在铸坯矫直过程中，为了使铸坯从弯曲状态转变为直线状态而施加的外力。铸坯的凝固收缩和热膨胀会导致其内部产生应力，这些应力与外力相互作用，进一步加剧了机械应力的复杂性。机械应力对铸坯质量有着显著的影响。过大的机械应力可能导致铸坯产生裂纹、变形等缺陷。在拉坯过程中，如果拉坯力过大，超过铸坯的抗拉强度，会在铸坯底部或薄弱部位产生裂纹。在矫直过程中，若弯曲应力和剪切应力过大，会使铸坯在弯曲部位或中性层附近产生裂纹。机械应力还可能导致铸坯的尺寸精度下降，影响后续加工和使用性能。在实际连铸生产中，需要合理控制拉坯力、矫直力等工艺参数，优化铸坯与结晶器壁、支撑辊等的接触条件，以减小机械应力，提高铸坯质量。4.4综合应力分析在高温合金立式连铸过程中，铸坯所受的应力是热应力、机械应力等多种应力相互叠加的结果，综合应力的分布情况对铸坯质量有着至关重要的影响。通过数值模拟，深入分析铸坯在连铸过程中的综合应力分布规律，结果如图18所示。[此处插入铸坯不同时刻的综合应力分布云图，图18：铸坯不同时刻的综合应力分布云图（a）t=10s（b）t=30s（c）t=50s]从综合应力分布云图可以看出，在连铸初期（t=10s），铸坯整体温度较高，热应力和机械应力相对较小，综合应力分布较为均匀，数值也较低。随着连铸过程的推进（t=30s），铸坯表面与结晶器壁接触，温度迅速下降，热应力逐渐增大，同时在拉坯力和摩擦力等机械力的作用下，机械应力也有所增加，综合应力分布变得不均匀。在铸坯的拐角处和表面附近，由于热应力和机械应力的叠加作用，综合应力明显增大，呈现出应力集中的现象。到了连铸后期（t=50s），铸坯大部分已经凝固，温度梯度逐渐减小，热应力有所降低，但机械应力仍然存在，综合应力分布呈现出从表面向中心逐渐减小的趋势。在铸坯的底部，由于受到拉坯力的直接作用，综合应力较大，若超过铸坯材料的强度极限，就可能导致铸坯产生裂纹等缺陷。为了更直观地了解综合应力随时间和位置的变化情况，绘制了铸坯表面和中心位置的综合应力随时间变化曲线（图19）。从曲线可以看出，铸坯表面的综合应力在开始阶段迅速增大，这是由于表面温度下降快，热应力迅速增加，同时受到机械力的作用，综合应力急剧上升。随着时间的推移，表面温度逐渐趋于稳定，热应力减小，但机械应力仍然存在，综合应力在一定数值范围内波动。铸坯中心位置的综合应力变化相对较为平缓，在连铸初期，由于温度较高，热应力和机械应力较小，综合应力也较小。随着凝固的进行，中心温度逐渐降低，热应力和机械应力逐渐增大，综合应力也相应增大，在凝固后期，随着温度梯度的减小，热应力减小，综合应力也逐渐减小。[此处插入铸坯表面和中心位置的综合应力随时间变化曲线，图19：铸坯表面和中心位置的综合应力随时间变化曲线]综合应力的大小和分布直接影响铸坯的质量和性能。当综合应力超过铸坯材料的屈服强度时，铸坯会发生塑性变形；若综合应力继续增大，超过材料的抗拉强度，就会导致铸坯产生裂纹。裂纹的产生会严重降低铸坯的力学性能，影响后续加工和使用。在实际连铸生产中，需要采取有效的措施来控制综合应力，如优化冷却制度，使铸坯各部分冷却均匀，减小热应力；合理调整拉坯力、矫直力等工艺参数，减小机械应力；对铸坯进行适当的热处理，消除残余应力。通过这些措施，可以降低综合应力，提高铸坯质量，减少缺陷的产生，提高生产效率和产品质量。综上所述，铸坯在连铸过程中的综合应力分布呈现出明显的规律性，随时间和位置的变化而变化。综合应力对铸坯质量有着重要影响，在高温合金立式连铸过程中，必须充分考虑综合应力的因素，采取相应的措施来降低综合应力，确保铸坯质量。五、结果验证与讨论5.1模拟结果验证为了确保数值模拟结果的准确性和可靠性，精心设计并开展了高温合金立式连铸实验。实验装置主要包括感应熔炼炉、立式连铸机、结晶器、冷却系统、温度测量系统和应力测量系统等。感应熔炼炉用于熔炼高温合金，将原料加热至预定的浇铸温度；立式连铸机则负责实现高温合金的连续铸造过程，通过拉坯装置将铸坯从结晶器中拉出；结晶器采用水冷铜质结构，为高温合金液的凝固提供强制冷却条件；冷却系统配备了循环水泵、水箱和冷却管道，确保冷却介质能够稳定地循环流动，带走铸坯凝固过程中释放的热量；温度测量系统采用高精度热电偶，在铸坯的不同位置进行布置，实时测量铸坯在连铸过程中的温度变化；应力测量系统则使用应变片，粘贴在铸坯的关键部位，测量铸坯在连铸过程中所受的应力。在实验过程中，严格控制工艺参数，使其与数值模拟中的参数保持一致。浇铸温度设定为1550℃，拉坯速度为0.05m/s，冷却水量为10L/min。通过热电偶和应变片，准确测量了铸坯在不同位置和不同时刻的温度和应力数据。在铸坯的表面、距表面10mm处和中心位置分别布置热电偶，测量这些位置的温度随时间的变化；在铸坯的底部和表面等关键部位粘贴应变片，测量拉坯和矫直过程中的机械应力。将实验测量得到的温度和应力数据与数值模拟结果进行详细对比，结果如图20和图21所示。从温度对比结果（图20）可以看出，实验测量值与模拟值在整体趋势上高度吻合。在连铸初期，铸坯温度迅速下降，模拟值和实验值都呈现出相似的下降趋势，且数值差异较小。随着连铸过程的进行，铸坯温度逐渐趋于稳定，模拟值和实验值的偏差也保持在较小的范围内。在铸坯的表面，实验测量的温度略低于模拟值，这可能是由于实验过程中存在一定的散热损失，以及热电偶测量存在一定的误差。[此处插入温度对比图，图20：实验测量与模拟的铸坯温度对比图]从应力对比结果（图21）可以观察到，实验测量的应力值与模拟值在变化趋势上也基本一致。在拉坯过程中，铸坯底部的机械应力迅速增大，模拟值和实验值都反映了这一变化特征，且应力大小的差异在可接受范围内。在矫直过程中，铸坯弯曲部位的应力变化也得到了较好的模拟，模拟值和实验值能够较好地匹配。然而，在某些局部位置，实验测量值与模拟值存在一定的偏差，这可能是由于实验过程中铸坯的实际受力情况较为复杂，存在一些难以准确测量和模拟的因素，如铸坯与结晶器壁之间的摩擦力分布不均匀、支撑辊的接触压力存在波动等。[此处插入应力对比图，图21：实验测量与模拟的铸坯应力对比图]通过对实验测量值与模拟值的偏差进行深入分析，发现主要的偏差来源包括测量误差、模型简化和参数不确定性等。测量误差主要来自热电偶和应变片的精度限制，以及测量过程中的干扰因素，如电磁干扰、机械振动等。模型简化方面，在建立数值模拟模型时，对一些复杂的物理现象进行了简化处理，如忽略了结晶器的圆角和其他次要结构，以及对钢液流动和传热过程中的一些细节进行了简化，这些简化可能导致模拟结果与实际情况存在一定的偏差。参数不确定性则是由于高温合金的材料参数在不同