硅钢薄带双辊铸轧过程的数值模拟：工艺参数与凝固组织的深度解析

硅钢薄带双辊铸轧过程的数值模拟：工艺参数与凝固组织的深度解析一、引言1.1研究背景与意义在现代钢铁工业中，硅钢作为一种重要的软磁材料，被广泛应用于电力、电子等众多领域。其优异的电磁性能，如高磁导率、低铁损等，使其成为变压器、电机等电气设备不可或缺的核心材料，在国民生产中有着十分重要的地位，被誉为钢铁工业的艺术品。随着社会经济的快速发展以及各行业对电气设备性能要求的不断提升，对硅钢的质量和性能也提出了更为严苛的要求。传统的硅钢生产工艺往往存在流程冗长、能耗巨大以及成本高昂等诸多问题，这在一定程度上限制了硅钢产业的高效发展。双辊铸轧技术作为一种极具创新性的近终形加工技术，为硅钢生产带来了新的变革契机。该技术以液态金属为起始原料，通过两个反向旋转的铸辊作为结晶器，直接将液态金属转化为薄带钢，实现了铸轧工序的高度一体化。与传统生产工艺相比，双辊铸轧技术展现出诸多显著优势。一方面，它极大地缩短了生产流程，取消了传统工艺中复杂的热轧环节，避免了铸带加热工序，从而有效降低了能源消耗，减少了碳排放，充分契合了当前钢铁行业绿色低碳发展的理念；另一方面，其独特的亚快速凝固特性（凝固速度可达100-1000℃/s），能够对铸带的微观组织进行灵活调控，显著提升硅钢的组织性能，例如可以有效控制等轴晶、柱状晶形态，减少偏析现象，抑制第二相的析出，进而提高硅钢最终产品的有利组织和织构，增强磁感，降低铁损。这些优势使得双辊铸轧技术在提高生产效率、降低生产成本以及提升产品质量等方面具有巨大的潜力，正逐渐成为钢铁行业研究与发展的重点方向，有望为新一代高性能硅钢产品开发提供解决方案。然而，双辊铸轧过程涉及到复杂的多物理场耦合现象，包括电磁场、流场、温度场以及应力应变场等，各物理场之间相互作用、相互影响，使得铸轧过程的内在机理极其复杂。同时，工艺参数如浇注温度、铸轧速度、熔池液面高度，以及设备结构参数如水口的出口角度和浸入深度、铸辊的材质与结构等众多因素，都会对铸轧过程中金属熔体的流动、传热和凝固行为产生显著影响，进而决定着最终产品的质量和性能。在实际生产中，由于铸轧过程工序高度凝聚，工艺参数间联系紧密，一些重要的物理现象和过程难以通过现有的实验条件进行全面、准确的测量和观察，这给深入理解铸轧过程的本质规律以及优化工艺参数和设备结构带来了极大的困难。数值模拟技术作为一种强大的研究工具，能够在计算机虚拟环境中对双辊铸轧过程进行全面、细致的模拟分析。通过建立精确的数学模型，数值模拟可以深入研究不同工艺参数和设备结构条件下，硅钢薄带双辊铸轧过程中金属熔体的流动形态、热量传递规律、凝固进程以及应力应变分布等关键物理现象。这不仅有助于揭示双辊铸轧过程的内在机理，还能为实验研究提供重要的理论指导，帮助优化实验方案，减少实验次数和成本。更为重要的是，通过数值模拟可以预测不同工艺参数组合下的铸轧结果，提前评估产品质量，为实际生产提供科学依据，从而实现工艺参数和设备结构的优化，提高铸坯质量，降低生产成本，增强产品在市场中的竞争力。因此，开展硅钢薄带双辊铸轧过程的数值模拟研究具有至关重要的理论意义和实际应用价值。在理论层面，它能够深化对多物理场耦合作用下双辊铸轧过程复杂机理的认识，丰富和完善金属材料加工的理论体系；在实际应用方面，可为硅钢双辊铸轧技术的工业化应用和推广提供坚实的技术支撑，推动钢铁行业朝着高效、节能、绿色的方向发展。1.2国内外研究现状自19世纪中叶亨利・贝塞麦发明双辊铸轧薄带钢专利以来，双辊铸轧技术经历了漫长的发展历程。早期，由于受到技术水平和工业基础的限制，该技术在工业化应用方面进展缓慢。直到近20年来，随着材料科学、计算机技术、自动化控制技术等相关领域的飞速发展，双辊铸轧技术才取得了实质性的突破，多条试验线和半工业生产线相继建成，逐渐接近工业化生产水平，并在全球范围内引起了广泛关注。国外在双辊铸轧技术的研究和应用方面起步较早，投入了大量的人力、物力和财力进行技术开发和理论研究，取得了众多令人瞩目的成果。美国纽柯公司作为该领域的佼佼者，其双辊薄带铸轧生产线于2002年正式投产，截至目前已累计生产近50万t薄带，所生产的最薄钢带厚度可达0.7mm，产品在一定程度上能够替代冷轧钢板，在市场上占据了重要份额。此外，欧洲、日本等工业发达国家也纷纷开展双辊铸轧技术的研究工作，建立了多个先进的研究机构和实验平台，对铸轧过程中的工艺参数优化、设备结构改进、产品质量提升等方面进行了深入研究，为双辊铸轧技术的工业化应用奠定了坚实的基础。在国内，双辊铸轧技术的研究虽然起步相对较晚，但近年来发展迅速。东北大学轧制技术及连轧自动化国家重点实验室在双辊薄带钢铸轧领域开展了大量的基础研究工作，取得了一系列创新性成果。在工艺研究方面，对薄带铸轧无取向硅钢、取向硅钢和特殊用途硅钢三类钢种进行了探索性研究，系统研究了铸带特殊取向晶粒在冷轧过程中的晶体转动和剪切带变形行为，提出了控制铸轧参数获得理想退火织构和磁性能的工艺条件；在取向硅钢研究中，系统研究了薄带连铸条件下Goss织构的演变规律，创新性提出了取向硅钢超低碳的成分设计，取消了后续脱碳退火，进一步简化了生产流程；在特种硅钢研究中，通过薄带铸轧独有的{100}取向晶粒可控异常长大现象，结合张力往复轧制和简单热处理工艺实现了Cube取向晶粒的比例调控，从而获得优异磁性能。这些研究成果为薄带铸轧硅钢工业化技术研究和实践提供了重要的理论指导。数值模拟技术作为研究双辊铸轧过程的重要手段，在国内外都得到了广泛的应用。通过数值模拟，可以深入研究铸轧过程中的多物理场耦合现象，揭示金属熔体的流动、传热和凝固规律，为工艺参数优化和设备结构设计提供科学依据。在国外，许多研究机构和企业利用先进的数值模拟软件，如ANSYS、FLUENT等，对双辊铸轧过程中的电磁场、流场、温度场以及应力应变场等进行了深入研究。例如，[具体文献]利用数值模拟方法研究了浇注温度、铸轧速度等工艺参数对熔池内流场和温度场的影响规律，为工艺参数的优化提供了理论支持。在国内，众多科研团队也积极开展双辊铸轧过程的数值模拟研究。东北大学的研究人员采用三维流热耦合有限元分析理论对双辊铸轧薄带钢凝固过程中流场和温度场进行了数值模拟，给出了浇注温度、铸轧速度、熔池液面高度、水口的出口角度、水口的浸入深度及辊套材质等参数对熔池内流场和温度场的影响规律；大连理工大学的学者建立了描述双辊薄带连铸过程的三维流热耦合模型，研究了不同工艺参数和布流系统下的连铸过程，分析了各参数对流场、温度场、液面波动及侧封板处冲击情况的影响。然而，当前双辊铸轧技术的研究仍存在一些不足之处。一方面，虽然数值模拟技术在铸轧过程研究中得到了广泛应用，但由于双辊铸轧过程的复杂性，多物理场耦合模型的建立还不够完善，一些关键的物理现象和相互作用机制尚未得到充分揭示，导致模拟结果与实际生产情况存在一定的偏差。另一方面，在工艺参数优化和设备结构设计方面，虽然已经取得了一些成果，但仍缺乏系统性和综合性的研究，难以满足工业化生产对产品质量和生产效率的更高要求。此外，对于双辊铸轧过程中的一些关键技术问题，如侧封技术、铸辊的使用寿命和冷却方式等，还需要进一步深入研究和改进。1.3研究内容与方法本研究将围绕硅钢薄带双辊铸轧过程展开，主要研究内容如下：双辊铸轧过程数值模型的建立：综合考虑双辊铸轧过程中涉及的复杂物理现象，如金属熔体的流动、传热、凝固以及应力应变等，运用计算流体力学（CFD）和有限元方法（FEM），建立能够准确描述双辊铸轧过程的三维多物理场耦合数值模型。在模型中，详细定义各物理场的控制方程，包括连续性方程、动量守恒方程、能量守恒方程等，并合理确定模型的初始条件和边界条件，如浇注温度、铸轧速度、冷却条件等，确保模型能够真实反映实际铸轧过程。工艺参数对铸轧过程的影响分析：利用所建立的数值模型，采用控制变量法，系统研究浇注温度、铸轧速度、熔池液面高度等关键工艺参数对双辊铸轧过程中金属熔体的流场、温度场、凝固行为以及最终产品质量的影响规律。通过改变单一工艺参数，固定其他参数不变，对不同工况下的铸轧过程进行数值模拟，分析模拟结果，揭示各工艺参数与铸轧过程物理现象之间的内在联系。例如，研究浇注温度的变化如何影响金属熔体的流动性和凝固速度，进而影响铸带的微观组织和性能；分析铸轧速度的改变对熔池内的流场分布、温度梯度以及铸带的应力应变状态的影响等。模型的实验验证与优化：为了确保数值模型的可靠性和准确性，设计并开展双辊铸轧实验。在实验过程中，严格控制实验条件，使其与数值模拟中的参数设置保持一致，测量铸轧过程中的关键物理量，如铸带的温度分布、厚度变化、微观组织等。将实验测量结果与数值模拟结果进行对比分析，验证模型的正确性。若发现模拟结果与实验数据存在偏差，深入分析原因，对模型进行优化和改进，如调整模型的参数设置、完善控制方程等，直到模型能够准确预测铸轧过程的物理现象和产品质量。在研究方法上，本研究将主要采用数值模拟与实验研究相结合的方式。数值模拟方面，选用专业的数值模拟软件，如ANSYS、FLUENT等，充分利用其强大的计算功能和丰富的物理模型库，对双辊铸轧过程进行全面、深入的模拟分析。在实验研究中，搭建双辊铸轧实验平台，配备先进的测量设备，如热电偶、红外测温仪、电子万能试验机等，对铸轧过程中的物理量进行精确测量。同时，运用控制变量法，严格控制实验条件，确保实验结果的准确性和可靠性。通过数值模拟与实验研究的相互验证和补充，深入揭示硅钢薄带双辊铸轧过程的内在机理，为工艺参数优化和设备结构改进提供科学依据。二、双辊铸轧技术与数值模拟理论基础2.1双辊铸轧技术概述2.1.1双辊铸轧工艺原理双辊铸轧技术作为一种先进的近终形加工技术，其工艺原理独特而高效。该技术以液态金属为起始原料，通过特定的浇注系统将其精准地注入到一对反向旋转且内部通有冷却水的铸辊之间。这对铸辊犹如两个高效的结晶器，液态金属在进入辊缝后，迅速与低温的铸辊表面接触，开始快速散热并凝固。在凝固的同时，由于铸辊的持续旋转，对正在凝固的金属施加了压力，使其在结晶凝固的瞬间承受压力加工和塑性变形，从而在极短的时间内完成从液态金属到固态薄带的全部过程，实现了铸轧工序的一体化。双辊铸轧技术与传统的铸造成形工艺相比，具有显著的优势。在传统铸造中，液态金属先在特定的模具中凝固成铸坯，然后需要经过加热、热轧等多个工序才能最终加工成所需的板材。而双辊铸轧技术直接将液态金属转化为薄带，取消了传统工艺中复杂的热轧环节，避免了铸带加热工序，大大缩短了生产流程。这不仅减少了生产设备的投入和占地面积，还降低了能源消耗和生产成本。同时，由于双辊铸轧过程中的凝固速度极快，通常可达100-1000℃/s，这种亚快速凝固特性使得铸带的微观组织更加均匀、细小，偏析现象得到有效抑制，从而显著提升了产品的性能和质量。从生产流程的角度来看，双辊铸轧技术的高效性更加明显。传统的薄带生产工艺往往需要经过连铸、加热、热轧、冷轧等多个步骤，每个步骤都需要消耗大量的时间、能源和资源。而双辊铸轧技术将这些步骤高度浓缩，从液态金属到薄带成品的生产过程一气呵成，极大地提高了生产效率。以某钢铁企业为例，采用传统工艺生产薄带时，从钢水到成品需要经过多个车间的流转，生产周期长达数天；而采用双辊铸轧技术后，生产周期缩短至数小时，生产效率大幅提升。在能源消耗方面，双辊铸轧技术也具有突出的优势。传统工艺中，铸坯的加热和热轧过程需要消耗大量的能源，而双辊铸轧技术取消了这些环节，直接将液态金属铸轧成薄带，能源消耗仅为传统流程的1/8左右。这不仅降低了企业的生产成本，还减少了碳排放，符合当前全球倡导的绿色低碳发展理念。双辊铸轧技术以其独特的工艺原理，在缩短生产流程、降低能源消耗、提高产品质量等方面展现出巨大的优势，正逐渐成为钢铁行业发展的重要方向。2.1.2硅钢薄带双辊铸轧的特点及应用硅钢薄带双辊铸轧技术作为双辊铸轧技术在硅钢生产领域的应用，具有一系列独特的特点，使其在高性能硅钢的生产中展现出巨大的优势。该技术最显著的特点之一是能够生产高牌号硅钢，满足高端市场对硅钢性能的严苛要求。传统的硅钢生产工艺在生产高牌号硅钢时，往往面临诸多挑战，如晶粒粗大、织构不理想等问题，这些问题会导致硅钢的磁性能下降，无法满足变压器、电机等高端电气设备的需求。而双辊铸轧技术凭借其亚快速凝固特性，能够有效控制铸带的微观组织，细化晶粒，抑制第二相的析出，从而提高硅钢的有利组织和织构。例如，在生产取向硅钢时，双辊铸轧技术可以使Goss织构更加均匀、稳定，显著提高硅钢的磁感强度，降低铁损，满足变压器铁芯对高磁感、低铁损硅钢的需求。在微观组织控制方面，双辊铸轧硅钢薄带表现出明显的优势。与传统工艺相比，其铸带的晶粒更加细小、均匀，偏析现象得到有效抑制。研究表明，双辊铸轧硅钢薄带的晶粒尺寸可以控制在几十微米以下，而传统工艺生产的硅钢晶粒尺寸往往在几百微米以上。这种细小均匀的晶粒结构使得硅钢在受力时能够更加均匀地分布应力，提高了硅钢的强度和韧性，同时也有利于改善硅钢的磁性能。在实际应用中，硅钢薄带双辊铸轧产品在变压器铁芯、电机制造等领域发挥着重要作用。在变压器铁芯制造中，高磁感、低铁损的双辊铸轧硅钢薄带能够有效提高变压器的效率，降低能源损耗。据统计，使用双辊铸轧硅钢薄带制造的变压器，其铁损可比传统硅钢变压器降低10%-20%，这对于电力系统的节能降耗具有重要意义。在电机制造领域，双辊铸轧硅钢薄带能够提高电机的效率和功率密度，使电机更加小型化、轻量化，满足现代工业对高效、节能电机的需求。以新能源汽车电机为例，采用双辊铸轧硅钢薄带制造的电机，能够提高电机的效率，延长汽车的续航里程，推动新能源汽车产业的发展。硅钢薄带双辊铸轧技术以其独特的特点，在高性能硅钢的生产中具有不可替代的优势，其产品在变压器铁芯、电机制造等领域的广泛应用，为推动相关产业的发展做出了重要贡献。随着技术的不断进步和完善，硅钢薄带双辊铸轧技术有望在更多领域得到应用，为国民经济的发展提供更加强有力的支持。2.2数值模拟基本理论2.2.1控制方程在双辊铸轧过程中，金属熔体的流动、传热和凝固过程涉及到多个物理场的相互作用，需要通过一系列控制方程来进行描述。这些控制方程是基于质量守恒、动量守恒和能量守恒等基本物理定律建立的，它们能够准确地反映双辊铸轧过程中的物理现象，为数值模拟提供了重要的理论基础。连续性方程是描述质量守恒的基本方程，其表达式为：\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0其中，\rho表示流体的密度，t为时间，\vec{v}是速度矢量。该方程表明，在单位时间内，流入和流出控制体的质量之差等于控制体内质量的变化率。在双辊铸轧过程中，连续性方程确保了金属熔体在流动过程中的质量守恒，对于理解熔池内金属的流动形态和分布具有重要意义。动量守恒方程，也称为Navier-Stokes方程，描述了流体在流动过程中的动量变化。在直角坐标系下，其表达式为：\frac{\partial(\rhov_i)}{\partialt}+\nabla\cdot(\rhov_i\vec{v})=-\frac{\partialp}{\partialx_i}+\nabla\cdot(\mu\nablav_i)+\rhog_i+S_{M,i}其中，v_i是速度矢量\vec{v}在i方向（i=x,y,z）上的分量，p为压力，\mu是动力粘度，g_i是重力加速度在i方向上的分量，S_{M,i}是动量源项。该方程考虑了惯性力、压力梯度、粘性力、重力以及其他外力对流体动量的影响。在双辊铸轧过程中，动量守恒方程用于描述金属熔体在铸辊间的流动行为，分析熔池内的速度分布和压力变化，对于研究铸轧过程中的轧制力、铸带的变形等问题具有重要作用。能量守恒方程用于描述传热过程，其表达式为：\frac{\partial(\rhoh)}{\partialt}+\nabla\cdot(\rhoh\vec{v})=\nabla\cdot(k\nablaT)+S_{E}其中，h是比焓，k是热导率，T为温度，S_{E}是能量源项。该方程表明，单位时间内控制体内能量的变化等于通过传导、对流和辐射等方式传递的热量以及其他能量源的输入。在双辊铸轧过程中，能量守恒方程用于计算金属熔体的温度分布，分析铸轧过程中的热量传递和凝固现象，对于控制铸带的凝固质量和组织性能具有重要意义。在双辊铸轧过程中，金属熔体的流动通常呈现出湍流状态。为了准确描述湍流流动，需要引入湍流模型。常用的湍流模型有k-\epsilon模型，该模型通过求解湍流动能k和湍流动能耗散率\epsilon的输运方程来封闭控制方程组。湍流动能k的方程为：\frac{\partial(\rhok)}{\partialt}+\nabla\cdot(\rhok\vec{v})=\nabla\cdot(\frac{\mu_t}{\sigma_k}\nablak)+G_k-\rho\epsilon其中，\mu_t是湍流粘度，\sigma_k是湍流动能的普朗特数，G_k是由平均速度梯度产生的湍流动能。湍流动能耗散率\epsilon的方程为：\frac{\partial(\rho\epsilon)}{\partialt}+\nabla\cdot(\rho\epsilon\vec{v})=\nabla\cdot(\frac{\mu_t}{\sigma_{\epsilon}}\nabla\epsilon)+C_{1\epsilon}\frac{\epsilon}{k}G_k-C_{2\epsilon}\rho\frac{\epsilon^2}{k}其中，\sigma_{\epsilon}是湍流动能耗散率的普朗特数，C_{1\epsilon}和C_{2\epsilon}是经验常数。在双辊铸轧过程中，金属熔体的凝固过程是一个复杂的相变过程，涉及到潜热的释放和固相率的变化。为了描述凝固过程，通常采用焓-孔隙率法。该方法将液相和固相视为连续介质，通过引入固相率f_s来描述凝固过程。固相率f_s的定义为：f_s=\frac{\rho_s}{\rho}其中，\rho_s是固相的密度。在凝固过程中，能量守恒方程中的源项S_{E}包括凝固潜热的释放，其表达式为：S_{E}=-\rhoL\frac{\partialf_s}{\partialt}其中，L是凝固潜热。通过上述控制方程，可以全面、准确地描述双辊铸轧过程中金属熔体的流动、传热和凝固过程。这些方程相互耦合，共同决定了双辊铸轧过程的物理行为。在数值模拟中，需要对这些方程进行离散化处理，并采用合适的数值算法进行求解，以获得双辊铸轧过程中各物理量的分布和变化规律。2.2.2关键问题处理在双辊铸轧过程的数值模拟中，固相率、凝固潜热和湍流粘度等关键问题的处理方法对于准确模拟铸轧过程起着至关重要的作用。这些问题涉及到复杂的物理现象和相互作用，需要采用合理的模型和算法进行处理，以确保模拟结果的准确性和可靠性。固相率作为描述金属凝固程度的重要参数，其准确计算对于理解铸轧过程中的凝固行为和预测铸带质量至关重要。在实际计算中，通常采用Scheil方程来计算固相率。Scheil方程基于溶质再分配理论，考虑了凝固过程中溶质在固相和液相之间的分配行为，能够较为准确地描述非平衡凝固条件下的固相率变化。其表达式为：f_s=1-\left(1-\frac{T-T_s}{T_l-T_s}\right)^{\frac{1}{1-k}}其中，T是当前温度，T_s和T_l分别是固相和液相的温度，k是溶质分配系数。凝固潜热是金属在凝固过程中释放的热量，它对铸轧过程中的温度场和凝固进程有着显著影响。为了准确处理凝固潜热，在数值模拟中通常采用焓-孔隙率法。该方法将凝固潜热视为能量源项，通过求解能量守恒方程来考虑其对温度场的影响。具体来说，在能量守恒方程中，将凝固潜热项S_{E}=-\rhoL\frac{\partialf_s}{\partialt}加入到能量源项中，其中\rho是密度，L是凝固潜热，f_s是固相率。通过这种方式，可以准确地模拟凝固潜热的释放对温度场的影响，从而更真实地反映铸轧过程中的传热和凝固现象。湍流粘度是描述湍流流动特性的重要参数，它反映了湍流对动量传输的影响。在双辊铸轧过程中，金属熔体的流动通常呈现出湍流状态，因此准确计算湍流粘度对于描述熔体的流动行为至关重要。在k-\epsilon湍流模型中，湍流粘度\mu_t与湍流动能k和湍流动能耗散率\epsilon相关，其表达式为：\mu_t=C_{\mu}\rho\frac{k^2}{\epsilon}其中，C_{\mu}是经验常数。通过求解湍流动能k和湍流动能耗散率\epsilon的输运方程，可以得到它们在流场中的分布，进而计算出湍流粘度的分布。这样可以准确地描述湍流对熔体流动的影响，提高模拟结果的准确性。固相率、凝固潜热和湍流粘度等关键问题的准确处理是保证双辊铸轧过程数值模拟准确性的关键。通过采用合理的模型和算法，能够更真实地反映铸轧过程中的物理现象，为工艺参数优化和设备结构设计提供可靠的理论依据。2.2.3有限元分析理论双辊铸轧过程是一个涉及复杂多物理场耦合的过程，包括流动、传热和凝固等多个方面。为了准确模拟这一过程，需要运用有限元分析理论，将连续的物理场离散为有限个单元，通过对每个单元的求解来逼近整个物理场的解。在双辊铸轧过程的有限元分析中，首先需要将求解区域划分为有限个单元，常用的单元类型有四面体单元、六面体单元等。以三维问题为例，假设求解区域被划分为N个单元，每个单元内的物理量可以通过节点值进行插值得到。对于速度场\vec{v}、温度场T等物理量，通常采用线性插值或高阶插值函数进行逼近。在离散化之后，需要对每个单元建立相应的方程。以动量守恒方程为例，在每个单元内，根据变分原理，可以将动量守恒方程转化为弱形式：\int_{\Omega_e}\rho\frac{\partial\vec{v}}{\partialt}\cdot\vec{\deltav}d\Omega+\int_{\Omega_e}\rho(\vec{v}\cdot\nabla)\vec{v}\cdot\vec{\deltav}d\Omega=-\int_{\Omega_e}\nablap\cdot\vec{\deltav}d\Omega+\int_{\Omega_e}\mu\nabla\vec{v}\cdot\nabla\vec{\deltav}d\Omega+\int_{\Omega_e}\rho\vec{g}\cdot\vec{\deltav}d\Omega+\int_{\Omega_e}\vec{S}_{M}\cdot\vec{\deltav}d\Omega其中，\Omega_e表示第e个单元的体积，\vec{\deltav}是速度的变分。通过对每个单元的弱形式方程进行组装，可以得到整个求解区域的代数方程组。对于能量守恒方程和其他相关方程，也可以采用类似的方法建立单元方程并进行组装。在求解这些代数方程组时，通常采用迭代法，如高斯-赛德尔迭代法、共轭梯度法等。通过不断迭代，逐步逼近方程的解，直到满足收敛条件为止。在双辊铸轧过程的数值模拟中，还需要考虑多物理场之间的耦合关系。例如，温度场的变化会影响金属熔体的物性参数，如密度、粘度等，而这些物性参数的变化又会反过来影响流动场和凝固过程。因此，在数值模拟中需要采用合适的耦合算法，如顺序耦合算法、强耦合算法等，来处理多物理场之间的相互作用。有限元分析理论为双辊铸轧过程的数值模拟提供了有效的方法。通过合理的单元划分、方程建立和求解，以及多物理场耦合的处理，可以准确地模拟双辊铸轧过程中金属熔体的流动、传热和凝固行为，为工艺参数优化和设备结构设计提供有力的支持。三、硅钢薄带双辊铸轧过程数值模型建立3.1模型假设与简化3.1.1物理模型简化在建立硅钢薄带双辊铸轧过程的数值模型时，为了便于计算和分析，需要对实际的物理模型进行合理的简化。实际的双辊铸轧过程涉及到众多复杂的物理现象和因素，若全面考虑所有因素，模型的计算量将极为庞大，甚至可能超出当前计算机的计算能力。因此，通过忽略一些对铸轧过程影响较小的次要因素，可以在保证模型准确性的前提下，显著提高计算效率。在双辊铸轧过程中，金属熔体从浇铸系统进入铸辊间的熔池，然后在铸辊的作用下凝固成薄带。在这个过程中，熔池内的金属熔体流动、传热以及凝固行为是研究的重点。为了简化模型，首先假设铸辊为刚性体，不考虑铸辊在轧制力作用下的弹性变形。这是因为在实际铸轧过程中，铸辊通常采用高强度的材料制造，其弹性变形相对较小，对铸轧过程的影响可以忽略不计。通过这一假设，可以大大简化模型的计算过程，提高计算效率。同时，忽略熔池内的表面张力和浮力的影响。表面张力主要影响熔池表面的形状和稳定性，而浮力则主要影响金属熔体在熔池内的垂直分布。在双辊铸轧过程中，由于铸轧速度较快，熔池内的金属熔体流动较为剧烈，表面张力和浮力的作用相对较小，对铸轧过程的整体影响可以忽略。在边界条件方面，将铸轧区的入口设定为速度入口，即给定金属熔体进入铸轧区的速度。这一速度通常根据实际生产中的铸轧工艺参数确定，它直接影响着金属熔体在熔池内的流动状态和传热过程。在出口处，设定为压力出口，给定出口压力。出口压力的大小会影响铸轧过程中的轧制力和铸带的质量，因此需要根据实际情况合理设定。对于铸辊与金属熔体的接触界面，假设为无滑移边界条件，即金属熔体与铸辊表面的速度相同。这一假设符合实际铸轧过程中金属熔体在铸辊表面的附着和跟随运动的情况。同时，考虑铸辊与金属熔体之间的热传递，设定接触界面的换热系数。换热系数的大小反映了铸辊与金属熔体之间的传热能力，它受到铸辊材质、表面粗糙度以及金属熔体的物性等多种因素的影响。在初始条件设定上，给定金属熔体的初始温度、速度和成分分布。初始温度通常根据浇注温度确定，它是影响金属熔体凝固过程的重要因素之一。初始速度与入口速度相同，而初始成分分布则根据所研究的硅钢的化学成分确定。这些初始条件的准确设定对于数值模拟结果的准确性至关重要。通过对物理模型的简化，合理确定边界条件和初始条件，可以建立起既能准确反映双辊铸轧过程主要物理现象，又具有较高计算效率的数值模型，为后续的研究工作奠定坚实的基础。3.1.2数学模型假设在建立双辊铸轧过程的数学模型时，为了能够运用已有的数学理论和方法对复杂的铸轧过程进行求解，需要提出一些合理的假设。这些假设不仅能够简化数学模型的形式，使其更易于处理，同时也能够在一定程度上反映铸轧过程的本质特征，确保模型的合理性和有效性。将液态金属视为连续不可压缩牛顿流体是一个重要的假设。从物理本质来看，在双辊铸轧过程中，液态金属的流动速度相对较低，且在铸轧过程中其密度变化较小。根据流体力学理论，当流体的流速远低于声速时，可近似认为流体是不可压缩的。在双辊铸轧过程中，金属熔体的流动速度通常在每秒数米的量级，远低于金属中的声速，因此将其视为不可压缩流体是合理的。同时，牛顿流体的假设意味着液态金属的粘性应力与应变速率呈线性关系。在大多数情况下，液态金属在常温下的粘性特性符合牛顿流体的定义。尽管在双辊铸轧过程中，液态金属的温度较高，但其粘性特性并没有发生根本性的改变，仍然可以用牛顿流体的模型来描述。这一假设使得我们可以运用经典的Navier-Stokes方程来描述液态金属的流动行为，大大简化了数学模型的建立和求解过程。忽略液态金属中的电磁效应也是常见的假设之一。在实际的双辊铸轧过程中，虽然可能存在一定的电磁现象，如感应电流、磁场等，但这些电磁效应通常较弱，对液态金属的流动和凝固过程的影响相对较小。一方面，双辊铸轧过程中并没有外部强磁场的作用，内部产生的感应电流和磁场强度都非常有限；另一方面，与流动、传热和凝固等主要物理过程相比，电磁效应的作用在量级上相对较小。因此，忽略电磁效应可以简化数学模型，同时不会对模拟结果的准确性产生显著影响。假设铸轧过程是稳态过程，即不考虑铸轧过程中的瞬态变化。在实际生产中，双辊铸轧过程在达到稳定运行状态后，其主要物理参数，如温度、速度、压力等，在一定时间内基本保持不变。虽然在铸轧过程的启动和停止阶段会存在瞬态变化，但这些阶段通常持续时间较短，对整个铸轧过程的影响相对较小。因此，在研究铸轧过程的稳态特性时，假设铸轧过程是稳态过程是合理的。这一假设使得我们可以在模型中忽略时间变量的导数项，简化了数学模型的求解过程，提高了计算效率。这些数学模型假设在简化模型的同时，充分考虑了双辊铸轧过程的物理本质和实际情况，具有一定的合理性和有效性。通过这些假设，可以建立起简洁而准确的数学模型，为深入研究双辊铸轧过程提供有力的工具。3.2模型构建与参数设置3.2.1几何模型构建为了深入研究硅钢薄带双辊铸轧过程，利用专业建模软件ANSYSICEM构建了双辊铸轧的几何模型。该模型全面涵盖了轧辊、侧封板和熔池等关键部分，旨在准确模拟实际铸轧过程中的物理现象。在构建轧辊模型时，充分考虑了其实际尺寸和形状。轧辊直径设定为[X]mm，这一尺寸是根据实际生产设备和工艺要求确定的，对铸轧过程中的轧制力、铸带的变形和冷却效果等都有着重要影响。轧辊长度为[X]mm，足够的长度能够保证铸带在轧制过程中的稳定性，避免出现边缘翘曲等缺陷。同时，轧辊表面被设置为光滑状态，以减少金属熔体与轧辊之间的摩擦阻力，使金属熔体能够在轧辊表面顺利流动并凝固。侧封板模型的构建同样至关重要。侧封板位于熔池两侧，其主要作用是防止金属熔体在铸轧过程中从侧面流出，确保熔池的稳定性。侧封板的高度为[X]mm，这一高度能够有效阻挡金属熔体的溢出，同时又不会对熔池内的流动和传热产生过大的阻碍。侧封板与轧辊之间的间隙被精确控制在[X]mm，合适的间隙既能保证金属熔体与侧封板之间的良好接触，又能避免因间隙过小而导致的过大摩擦和磨损。熔池模型作为金属熔体流动和凝固的关键区域，其构建也经过了精心设计。熔池的长度根据铸轧工艺参数和实际生产经验确定为[X]mm，这一长度能够确保金属熔体在进入辊缝之前有足够的时间进行充分的流动和混合，从而保证铸带的质量均匀性。熔池的宽度与轧辊长度相同，为[X]mm，以保证金属熔体在整个宽度方向上的均匀分布。熔池的深度则根据金属熔体的浇注量和铸轧速度等因素进行调整，在本模型中设定为[X]mm，合适的深度能够保证金属熔体在熔池内的稳定流动，同时也有利于热量的传递和凝固过程的进行。通过合理构建轧辊、侧封板和熔池等部分的几何模型，并精确设置各部分的尺寸参数，本研究建立的双辊铸轧几何模型能够较为准确地模拟实际铸轧过程中的物理现象，为后续的数值模拟分析提供了可靠的基础。3.2.2材料参数设定在硅钢薄带双辊铸轧过程的数值模拟中，准确设定硅钢和铸辊等材料的热物理参数是确保模拟结果准确性的关键。这些参数直接影响着金属熔体在铸轧过程中的流动、传热和凝固行为。硅钢作为铸轧的主要材料，其热物理参数随温度的变化而呈现出一定的规律。在液态阶段，硅钢的密度约为[X]kg/m³，这一密度值决定了金属熔体在熔池内的流动特性和浮力作用。随着温度的降低，硅钢逐渐凝固，其密度会发生一定的变化，在固态时密度约为[X]kg/m³。硅钢的比热容在液态时为[X]J/(kg・K)，它反映了硅钢吸收或释放热量的能力，对铸轧过程中的温度分布和凝固速度有着重要影响。在凝固过程中，比热容会随着固相率的增加而逐渐减小，在固态时比热容约为[X]J/(kg・K)。硅钢的导热系数在液态时为[X]W/(m・K)，它决定了热量在硅钢中的传递速度，较高的导热系数有助于热量的快速散发，促进凝固过程的进行。在固态时，导热系数约为[X]W/(m・K)，相对液态时有所降低。铸辊作为双辊铸轧过程中的关键部件，其材料的热物理参数同样对铸轧过程有着重要影响。铸辊通常采用具有良好导热性能和强度的材料，如铜合金。铜合金铸辊的密度约为[X]kg/m³，这一密度保证了铸辊在高速旋转过程中的稳定性。铸辊的比热容为[X]J/(kg・K)，较低的比热容使得铸辊能够快速吸收金属熔体的热量，从而实现高效的冷却效果。铸辊的导热系数高达[X]W/(m・K)，这使得热量能够迅速从铸辊表面传递到内部，进而传递到冷却介质中，保证铸轧过程的顺利进行。除了上述主要参数外，硅钢和铸辊材料的其他热物理参数，如热膨胀系数、潜热等，也在数值模拟中进行了合理设定。硅钢的热膨胀系数在不同温度区间有所不同，在室温至熔点之间，热膨胀系数约为[X]×10⁻⁶/K，这一参数影响着硅钢在加热和冷却过程中的尺寸变化。铸辊的热膨胀系数相对较小，约为[X]×10⁻⁶/K，以保证铸辊在高温和高速旋转条件下的尺寸稳定性。在凝固过程中，硅钢会释放大量的潜热，潜热的数值约为[X]J/kg，这部分热量的释放对凝固过程的温度场和凝固速度有着显著影响。通过对硅钢和铸辊等材料热物理参数的准确设定，能够更真实地反映双辊铸轧过程中的物理现象，为数值模拟结果的可靠性提供有力保障。3.2.3边界条件与初始条件设定在硅钢薄带双辊铸轧过程的数值模拟中，合理设定边界条件和初始条件是确保模型能够准确反映实际铸轧过程的关键环节。这些条件的设定直接影响着金属熔体在铸轧过程中的流动、传热和凝固行为。对于铸辊的冷却条件，考虑到实际生产中铸辊内部通有冷却水进行冷却，在模型中设定铸辊表面与冷却水之间的换热系数为[X]W/(m²・K)。这一换热系数的大小直接影响着铸辊向冷却水传递热量的速度，进而影响金属熔体在铸辊表面的凝固速度和温度分布。通过查阅相关文献和实际生产数据，确定该换热系数能够较为准确地模拟铸辊的冷却过程。同时，假设冷却水的温度恒定为[X]℃，这一温度条件保证了铸辊在整个铸轧过程中能够持续稳定地吸收金属熔体的热量，实现高效冷却。在熔池的入口条件方面，给定金属熔体的入口速度为[X]m/s，这一速度根据实际铸轧工艺参数确定，它直接影响着金属熔体在熔池内的流动状态和填充情况。入口温度设定为[X]℃，该温度略高于硅钢的熔点，确保金属熔体在进入熔池时处于液态，能够顺利进行铸轧过程。入口处的金属熔体成分也根据所研究的硅钢具体牌号进行准确设定，以保证模拟结果的准确性。熔池的出口条件设定为压力出口，出口压力为[X]Pa。这一压力条件影响着金属熔体在铸轧过程中的流动阻力和轧制力的大小。通过调整出口压力，可以模拟不同的铸轧工况，研究其对铸轧过程和铸带质量的影响。在初始条件设定上，整个计算区域的初始温度设定为[X]℃，这一温度接近金属熔体的浇注温度，反映了铸轧过程开始时的温度状态。初始速度场设定为零，即假设在铸轧开始瞬间，金属熔体尚未开始流动。随着模拟的进行，在边界条件的作用下，金属熔体逐渐开始流动并凝固。通过合理设定铸辊的冷却条件、熔池的入口和出口条件以及初始温度场和速度场等边界条件和初始条件，能够使建立的数值模型更真实地反映硅钢薄带双辊铸轧的实际过程，为后续的模拟分析和结果讨论提供可靠的基础。3.3模型验证3.3.1实验方案设计为了验证所建立的硅钢薄带双辊铸轧过程数值模型的准确性和可靠性，精心设计了一系列实验。实验设备选用了专门的双辊铸轧实验机，该实验机具备精确控制铸轧工艺参数的能力，能够满足实验研究的需求。实验机的主要参数包括：铸辊直径为[X]mm，与数值模型中的轧辊直径一致，以保证实验条件与模拟条件的一致性；铸辊长度为[X]mm，能够满足硅钢薄带的铸轧宽度要求；最大铸轧速度可达[X]m/min，可在一定范围内调整铸轧速度，研究其对铸轧过程的影响。在材料选择上，采用了实际生产中常用的硅钢材料，其化学成分和性能参数与数值模拟中所使用的硅钢材料相同。具体的化学成分如下：硅含量为[X]%，碳含量为[X]%，锰含量为[X]%，磷含量为[X]%，硫含量为[X]%，其余为铁及微量杂质。这些化学成分的精确控制对于保证实验结果的准确性至关重要。实验过程中，对工艺参数进行了严格控制。浇注温度设定为[X]℃，这一温度略高于硅钢的熔点，确保金属熔体在进入铸轧区时处于液态，能够顺利进行铸轧过程。铸轧速度分别设置为[X]m/min、[X]m/min和[X]m/min，通过改变铸轧速度，研究其对铸轧过程中金属熔体的流动、传热和凝固行为的影响。熔池液面高度保持在[X]mm，以保证熔池内金属熔体的稳定性和均匀性。为了准确测量铸轧过程中的关键物理量，在实验中布置了多个热电偶，用于测量铸带的温度分布。热电偶采用了高精度的K型热电偶，其测量精度可达±1℃，能够满足实验测量的要求。在铸带的不同位置，如铸带的中心、边缘以及不同的厚度方向上，分别布置了热电偶，以获取铸带在不同位置的温度变化情况。同时，使用激光测厚仪实时监测铸带的厚度变化，激光测厚仪的测量精度可达±0.01mm，能够准确测量铸带在铸轧过程中的厚度波动。通过这些测量设备，能够获取铸轧过程中的关键物理量，为模型验证提供可靠的数据支持。3.3.2实验结果与模拟结果对比分析将实验测量得到的铸带温度分布和厚度变化数据与数值模拟结果进行了详细对比分析，以验证数值模型的准确性和可靠性。在铸带温度分布方面，实验结果与模拟结果表现出较好的一致性。图1展示了在浇注温度为[X]℃、铸轧速度为[X]m/min的条件下，铸带中心线上的温度分布对比情况。从图中可以看出，实验测量的温度曲线与模拟计算得到的温度曲线在趋势上基本吻合，在铸轧区的入口处，由于金属熔体刚进入铸轧区，温度较高，随着铸轧过程的进行，金属熔体与铸辊表面接触，热量迅速传递给铸辊，温度逐渐降低。在铸轧区的出口处，铸带的温度已经接近铸辊的温度。在整个铸轧过程中，实验测量的温度与模拟结果的偏差在可接受范围内，最大偏差不超过[X]℃。这表明数值模型能够较为准确地预测铸带在铸轧过程中的温度变化情况。[此处插入铸带中心线温度分布对比图]在铸带厚度变化方面，实验结果与模拟结果也具有较高的一致性。图2为在不同铸轧速度下，铸带厚度的实验测量值与模拟计算值的对比情况。从图中可以看出，随着铸轧速度的增加，铸带的厚度略有增加。这是因为铸轧速度的增加导致金属熔体在铸轧区的停留时间缩短，凝固速度减慢，从而使得铸带的厚度增加。实验测量的铸带厚度与模拟计算结果的偏差较小，在铸轧速度为[X]m/min时，偏差为[X]mm；在铸轧速度为[X]m/min时，偏差为[X]mm；在铸轧速度为[X]m/min时，偏差为[X]mm。这些偏差在实际生产中是可以接受的，说明数值模型能够准确地预测铸轧速度对铸带厚度的影响。[此处插入不同铸轧速度下铸带厚度对比图]尽管实验结果与模拟结果总体上吻合较好，但仍存在一些细微的差异。这些差异可能是由多种因素导致的。在实验过程中，由于测量设备的精度限制，可能会引入一定的测量误差。热电偶的测量精度虽然可达±1℃，但在实际测量过程中，由于热电偶与铸带的接触情况、测量环境等因素的影响，可能会导致测量结果存在一定的偏差。激光测厚仪的测量精度可达±0.01mm，但在测量过程中，也可能会受到铸带表面的平整度、振动等因素的影响，从而产生测量误差。实验条件与数值模拟条件的不完全一致也可能导致结果的差异。在实验中，虽然尽量控制工艺参数与数值模拟中的参数相同，但在实际操作中，很难做到完全精确的控制。浇注温度可能会存在±[X]℃的波动，铸轧速度可能会有±[X]m/min的偏差，这些微小的差异都可能对铸轧过程产生一定的影响，从而导致实验结果与模拟结果的不同。数值模型本身也存在一定的局限性。在建立模型时，对一些复杂的物理现象进行了简化和假设，这些简化和假设可能会导致模型与实际情况存在一定的偏差。在模型中忽略了熔池内的表面张力和浮力的影响，虽然在大多数情况下这些影响可以忽略不计，但在某些特殊情况下，它们可能会对铸轧过程产生一定的作用。模型中的一些参数，如换热系数、粘度等，可能会受到实际生产条件的影响而发生变化，从而影响模型的准确性。通过对实验结果与模拟结果的对比分析，验证了所建立的硅钢薄带双辊铸轧过程数值模型的准确性和可靠性。虽然存在一些细微的差异，但这些差异可以通过进一步优化实验条件和改进数值模型来减小。该数值模型为深入研究硅钢薄带双辊铸轧过程提供了有力的工具，能够为工艺参数优化和设备结构设计提供可靠的理论依据。四、工艺参数对双辊铸轧过程的影响4.1浇注温度的影响4.1.1对熔池温度场的影响浇注温度作为双辊铸轧过程中的关键工艺参数之一，对熔池温度场有着显著的影响。通过数值模拟不同浇注温度下的双辊铸轧过程，深入分析其对熔池温度分布的作用机制，有助于揭示铸轧过程中的传热规律，为优化工艺参数提供理论依据。图1展示了在其他工艺参数保持不变的情况下，不同浇注温度（T_1、T_2、T_3，且T_1\ltT_2\ltT_3）下熔池的温度分布云图。从图中可以清晰地看出，随着浇注温度的升高，熔池的整体温度显著上升。在较低的浇注温度T_1下，熔池内的温度相对较低，靠近铸辊表面的区域温度迅速降低，形成了较大的温度梯度。这是因为在较低的浇注温度下，金属熔体进入熔池后，与低温的铸辊表面接触，热量迅速传递给铸辊，导致温度急剧下降。随着浇注温度升高到T_2，熔池内的温度分布更加均匀，温度梯度有所减小。这是由于较高的浇注温度使得金属熔体具有更多的热量，在与铸辊接触时，热量传递相对缓慢，从而使熔池内的温度分布更加均匀。当浇注温度进一步升高到T_3时，熔池的整体温度进一步升高，且高温区域明显扩大。此时，金属熔体的热量充足，在熔池内的流动过程中，能够保持较高的温度，使得高温区域向熔池内部扩展。[此处插入不同浇注温度下熔池温度分布云图]为了更直观地分析浇注温度对熔池温度场的影响，绘制了不同浇注温度下熔池中心线上的温度变化曲线，如图2所示。从图中可以看出，在熔池入口处，由于金属熔体刚进入熔池，温度与浇注温度相近，因此不同浇注温度下的温度曲线在此处几乎重合。随着金属熔体在熔池内的流动，与铸辊表面的接触逐渐增多，热量不断传递给铸辊，温度逐渐降低。在较低的浇注温度下，温度下降的速率较快，曲线斜率较大；而在较高的浇注温度下，温度下降的速率相对较慢，曲线斜率较小。这表明浇注温度越高，金属熔体在熔池内的散热速度越慢，能够保持较高的温度。在熔池出口处，较低浇注温度下的铸带温度明显低于较高浇注温度下的铸带温度，这将对铸带的凝固和质量产生重要影响。[此处插入不同浇注温度下熔池中心线温度变化曲线]过高或过低的浇注温度都会对铸轧过程产生不利影响。当浇注温度过低时，金属熔体的流动性变差，难以在熔池内均匀分布，容易导致铸带厚度不均匀，甚至出现冷隔、夹渣等缺陷。过低的浇注温度还会使铸带在出口处的温度过低，凝固过快，增加铸轧力，对铸辊和设备造成较大的负荷，影响设备的使用寿命。相反，当浇注温度过高时，金属熔体在熔池内的散热困难，凝固时间延长，可能导致铸带在出口处尚未完全凝固，出现漏钢现象，严重影响铸轧过程的稳定性和铸带的质量。过高的浇注温度还会使铸带的晶粒粗大，组织性能下降，降低产品的质量。4.1.2对凝固过程和铸带质量的影响浇注温度不仅对熔池温度场有着重要影响，还直接关系到铸轧过程中的凝固行为和铸带质量。通过数值模拟和理论分析，深入研究浇注温度对凝固区大小、凝固速度和铸带质量的影响规律，对于优化双辊铸轧工艺、提高铸带质量具有重要意义。随着浇注温度的升高，凝固区大小和凝固速度会发生显著变化。图3展示了不同浇注温度下凝固区的分布情况。从图中可以看出，在较低的浇注温度下，金属熔体进入熔池后，由于与铸辊表面的温度差较大，热量迅速传递给铸辊，凝固速度较快，凝固区相对较小。随着浇注温度的升高，金属熔体的初始温度升高，在熔池内的散热速度变慢，凝固速度相应降低，凝固区逐渐扩大。这是因为较高的浇注温度使得金属熔体具有更多的热量，需要更长的时间才能将热量传递给铸辊，从而延缓了凝固过程。[此处插入不同浇注温度下凝固区分布情况图]为了定量分析浇注温度对凝固速度的影响，绘制了不同浇注温度下凝固速度随时间的变化曲线，如图4所示。从图中可以看出，在铸轧开始阶段，不同浇注温度下的凝固速度都较快，这是由于金属熔体与铸辊表面接触，热量迅速传递，导致凝固速度迅速增加。随着时间的推移，较低浇注温度下的凝固速度逐渐趋于稳定，而较高浇注温度下的凝固速度下降较为缓慢。这表明浇注温度越高，凝固速度在后期下降得越慢，凝固过程持续的时间越长。[此处插入不同浇注温度下凝固速度随时间变化曲线]凝固过程的变化会直接影响铸带的质量。在较低的浇注温度下，由于凝固速度较快，铸带内部的溶质来不及均匀扩散，容易导致成分偏析，影响铸带的性能均匀性。快速的凝固还可能使铸带产生较大的内应力，导致铸带出现裂纹等缺陷。而在较高的浇注温度下，虽然成分偏析现象可能得到一定程度的改善，但过长的凝固时间会使铸带的晶粒长大，降低铸带的强度和韧性。过高的浇注温度还可能导致铸带表面出现氧化、气孔等缺陷，影响铸带的表面质量。为了获得高质量的铸带，需要合理控制浇注温度。根据数值模拟结果和实际生产经验，对于本研究中的硅钢薄带双辊铸轧过程，合适的浇注温度范围为[T_{min},T_{max}]。在这个温度范围内，既能保证金属熔体具有良好的流动性，在熔池内均匀分布，又能使凝固过程适中，减少成分偏析和内应力，控制晶粒尺寸，从而获得质量优良的铸带。当浇注温度低于T_{min}时，应适当提高浇注温度，以改善金属熔体的流动性和凝固条件；当浇注温度高于T_{max}时，应降低浇注温度，以避免铸带出现过热缺陷。4.2铸轧速度的影响4.2.1对熔池流场的影响铸轧速度作为双辊铸轧过程中的关键工艺参数之一，对熔池内金属液的流动速度和流线分布有着显著的影响，进而对铸带质量产生重要作用。通过数值模拟不同铸轧速度下的双辊铸轧过程，深入分析其对熔池流场的影响机制，对于优化铸轧工艺、提高铸带质量具有重要意义。图5展示了在其他工艺参数保持不变的情况下，不同铸轧速度（v_1、v_2、v_3，且v_1\ltv_2\ltv_3）下熔池内金属液的速度矢量分布。从图中可以清晰地看出，随着铸轧速度的增加，熔池内金属液的流动速度明显增大。在较低的铸轧速度v_1下，金属液在熔池内的流动相对缓慢，靠近铸辊表面的区域速度较低，形成了较小的速度梯度。这是因为在较低的铸轧速度下，铸辊对金属液的拖拽作用相对较弱，金属液的惯性较小，难以快速填充熔池，导致流动速度较慢。随着铸轧速度升高到v_2，金属液的流动速度显著增加，速度梯度也有所增大。这是由于较高的铸轧速度使得铸辊对金属液的拖拽作用增强，金属液在惯性的作用下，能够更快速地填充熔池，从而使流动速度加快。当铸轧速度进一步升高到v_3时，金属液的流动速度进一步增大，且在熔池内形成了较为复杂的流动形态。此时，金属液在高速流动的过程中，受到铸辊表面摩擦力、熔池边界条件以及自身惯性等多种因素的影响，形成了多个漩涡和回流区域，使得熔池内的流动更加紊乱。[此处插入不同铸轧速度下熔池内金属液速度矢量分布图]为了更直观地分析铸轧速度对熔池内流线分布的影响，绘制了不同铸轧速度下熔池内的流线图，如图6所示。从图中可以看出，在较低的铸轧速度下，流线分布相对较为规则，金属液从熔池入口进入后，沿着铸辊表面逐渐向下流动，在熔池底部形成较为平缓的流动轨迹。随着铸轧速度的增加，流线分布变得更加复杂，出现了明显的漩涡和回流现象。在较高的铸轧速度下，漩涡和回流区域的范围增大，数量增多，这表明金属液在熔池内的流动更加不均匀，容易导致金属液在熔池内的停留时间不一致，从而影响铸带质量。[此处插入不同铸轧速度下熔池内流线分布图]熔池内金属液的流动状态对铸带质量有着至关重要的影响。当铸轧速度较低时，金属液在熔池内的流动速度较慢，容易在熔池底部形成堆积，导致铸带厚度不均匀，甚至出现冷隔、夹渣等缺陷。金属液在熔池内的停留时间较长，容易造成热量散失过多，使铸带在出口处的温度过低，凝固过快，增加铸轧力，对铸辊和设备造成较大的负荷，影响设备的使用寿命。当铸轧速度过高时，金属液在熔池内的流动速度过快，形成的漩涡和回流现象会使金属液在熔池内的分布不均匀，导致铸带内部的组织和成分不均匀，影响铸带的性能。过高的铸轧速度还可能导致金属液在熔池出口处的流速过大，使铸带表面出现波动，影响铸带的表面质量。4.2.2对凝固组织和力学性能的影响铸轧速度不仅对熔池流场有着重要影响，还直接关系到铸轧过程中的凝固组织和铸带的力学性能。通过数值模拟和实验研究，深入分析铸轧速度对铸带凝固组织和力学性能的影响规律，对于优化双辊铸轧工艺、提高铸带质量具有重要意义。随着铸轧速度的增加，铸带的凝固组织会发生显著变化。图7展示了不同铸轧速度下铸带的微观组织形貌。从图中可以看出，在较低的铸轧速度下，铸带的晶粒尺寸较大，柱状晶比例较高。这是因为在较低的铸轧速度下，金属液在熔池内的停留时间较长，有足够的时间进行形核和长大，导致晶粒尺寸增大。同时，由于铸辊表面的冷却作用，靠近铸辊表面的金属液优先凝固，形成了柱状晶。随着铸轧速度的增加，铸带的晶粒尺寸逐渐减小，柱状晶比例降低，等轴晶比例增加。这是由于较高的铸轧速度使得金属液在熔池内的停留时间缩短，形核和长大的时间不足，导致晶粒尺寸减小。高速流动的金属液会产生强烈的搅拌作用，使熔池内的温度更加均匀，抑制了柱状晶的生长，促进了等轴晶的形成。[此处插入不同铸轧速度下铸带微观组织形貌图]为了定量分析铸轧速度对晶粒尺寸和柱状晶比例的影响，绘制了晶粒尺寸和柱状晶比例随铸轧速度的变化曲线，如图8所示。从图中可以看出，随着铸轧速度的增加，晶粒尺寸逐渐减小，呈现出近似线性的下降趋势。柱状晶比例也随着铸轧速度的增加而逐渐降低，在铸轧速度较低时，柱状晶比例下降较为明显，随着铸轧速度的进一步增加，柱状晶比例的下降趋势逐渐变缓。[此处插入晶粒尺寸和柱状晶比例随铸轧速度变化曲线]凝固组织的变化会直接影响铸带的力学性能。在较低的铸轧速度下，由于晶粒尺寸较大，柱状晶比例较高，铸带的强度和硬度相对较低，塑性和韧性较差。较大的晶粒尺寸使得晶界面积减小，晶界对位错运动的阻碍作用减弱，导致铸带在受力时容易发生塑性变形，强度和硬度降低。柱状晶的存在使得铸带在不同方向上的性能存在差异，容易出现各向异性，降低了铸带的塑性和韧性。而在较高的铸轧速度下，由于晶粒尺寸减小，等轴晶比例增加，铸带的强度和硬度显著提高，塑性和韧性也得到改善。较小的晶粒尺寸增加了晶界面积，晶界对位错运动的阻碍作用增强，使得铸带在受力时能够承受更大的应力，从而提高了强度和硬度。等轴晶的均匀分布使得铸带在各个方向上的性能更加一致，减少了各向异性，提高了塑性和韧性。为了获得良好的凝固组织和力学性能，需要合理控制铸轧速度。根据数值模拟结果和实际生产经验，对于本研究中的硅钢薄带双辊铸轧过程，合适的铸轧速度范围为[v_{min},v_{max}]。在这个速度范围内，既能保证铸带具有细小均匀的晶粒组织，又能使铸带具有良好的力学性能。当铸轧速度低于v_{min}时，应适当提高铸轧速度，以细化晶粒，改善凝固组织；当铸轧速度高于v_{max}时，应降低铸轧速度，以避免铸带出现组织和性能缺陷。4.3水口参数的影响4.3.1水口浸入深度对流场和温度场的影响水口浸入深度作为双辊铸轧过程中的重要参数之一，对熔池内的流场和温度场有着显著的影响，进而对铸带质量产生重要作用。通过数值模拟不同水口浸入深度下的双辊铸轧过程，深入分析其对熔池内物理场的影响机制，对于优化铸轧工艺、提高铸带质量具有重要意义。图9展示了在其他工艺参数保持不变的情况下，不同水口浸入深度（h_1、h_2、h_3，且h_1\lth_2\lth_3）下熔池内金属液的速度矢量分布。从图中可以清晰地看出，随着水口浸入深度的增加，熔池内金属液的流动状态发生明显变化。在较小的水口浸入深度h_1下，金属液从水口流出后，在熔池表面形成较强的冲击，导致熔池表面附近的流速较大，而熔池底部的流速相对较小，形成了较大的速度梯度。这是因为在较小的浸入深度下，水口流出的金属液直接冲击熔池表面，使得表面金属液获得较大的动能，而底部金属液受到的冲击较小，流动速度较慢。随着水口浸入深度增加到h_2，金属液的冲击位置逐渐下移，熔池表面的流速有所减小，底部的流速有所增加，速度梯度减小。这是由于金属液在更深的位置进入熔池，冲击能量在熔池内部得到更均匀的分布，使得熔池内的流动更加均匀。当水口浸入深度进一步增加到h_3时，金属液在熔池内的流动更加平稳，速度分布更加均匀。此时，金属液在较深的位置进入熔池，避免了对熔池表面的强烈冲击，减少了表面波动，使得熔池内的流场更加稳定。[此处插入不同水口浸入深度下熔池内金属液速度矢量分布图]水口浸入深度的变化还会对熔池内的温度场产生重要影响。图10展示了不同水口浸入深度下熔池的温度分布云图。从图中可以看出，在较小的水口浸入深度下，由于金属液直接冲击熔池表面，使得表面温度分布不均匀，高温区域集中在水口出口附近，而熔池底部的温度相对较低，形成了较大的温度梯度。随着水口浸入深度的增加，金属液在熔池内的分布更加均匀，热量传递更加充分，熔池内的温度分布逐渐趋于均匀，温度梯度减小。在较大的水口浸入深度下，熔池内的温度场更加稳定，有利于铸带的均匀凝固。[此处插入不同水口浸入深度下熔池温度分布云图]熔池内流场和温度场的变化会直接影响铸带质量。当水口浸入深度过小时，熔池表面的流速过大，容易产生较大的液面波动，导致铸带表面出现缺陷，如波纹、夹杂等。表面温度分布不均匀还可能导致铸带表面的凝固速度不一致，产生应力集中，增加铸带出现裂纹的风险。而当水口浸入深度过大时，虽然熔池内的流场和温度场更加稳定，但金属液在熔池内的停留时间可能过长，导致热量散失过多，铸带在出口处的温度过低，凝固过快，增加铸轧力，对铸辊和设备造成较大的负荷，影响设备的使用寿命。4.3.2水口角度对液面波动和侧封板冲击的影响水口角度作为双辊铸轧过程中的关键参数之一，对熔池内的液面波动和侧封板冲击有着重要影响，进而对铸带质量和生产稳定性产生显著作用。通过数值模拟不同水口角度下的双辊铸轧过程，深入分析其对熔池内物理现象的影响机制，对于优化铸轧工艺、提高铸带质量具有重要意义。图11展示了在其他工艺参数保持不变的情况下，不同水口角度（\theta_1、\theta_2、\theta_3，且\theta_1\lt\theta_2\lt\theta_3）下熔池内的液面波动情况。从图中可以清晰地看出，随着水口角度的变化，熔池内的液面波动呈现出明显的差异。在较小的水口角度\theta_1下，金属液从水口流出后，以较小的角度冲击熔池表面，导致熔池表面的局部流速较大，形成较大的液面波动。这是因为较小的水口角度使得金属液的冲击能量集中在较小的区域，对熔池表面产生较强的扰动，从而引起较大的液面波动。随着水口角度增加到\theta_2，金属液的冲击角度增大，冲击能量在熔池表面的分布更加均匀，液面波动有所减小。此时，金属液以更合理的角度进入熔池，对表面的扰动相对较小，使得液面波动得到一定程度的抑制。当水口角度进一步增加到\theta_3时，虽然冲击能量分布更加均匀，但由于金属液的冲击方向与熔池表面的夹角过大，可能导致金属液在熔池表面产生较大的横向流动，反而使液面波动再次增大。[此处插入不同水口角度下熔池内液面波动情况图]水口角度的变化还会对侧封板的冲击产生重要影响。图12展示了不同水口角度下侧封板受到的冲击压力分布。从图中可以看出，在较小的水口角度下，侧封板靠近水口一侧受到的冲击压力较大，这是因为金属液以较小的角度冲击熔池表面后，部分能量传递到侧封板上，且冲击位置相对集中，导致局部压力较大。随着水口角度的增加，冲击压力在侧封板上的分布逐渐均匀，最大冲击压力有所减小。然而，当水口角度过大时，侧封板另一侧可能受到较大的冲击压力，这是由于金属液的横向流动增强，对侧封板的另一侧产生了较大的冲击力。[此处插入不同水口角度下侧封板受到的冲击压力分布图]过大的液面波动和侧封板冲击会对铸带质量和生产稳定性产生不利影响。当液面波动过大时，铸带表面容易出现波纹、夹杂等缺陷，影响铸带的表面质量和尺寸精度。液面波动还可能导致铸带在凝固过程中出现不均匀的情况，增加铸带内部的应力，降低铸带的力学性能。而侧封板受到过大的冲击压力，可能导致侧封板磨损加剧，甚至出现泄漏现象，影响铸轧过程的稳定性和铸带的质量。为了获得稳定的铸轧过程和高质量的铸带，需要合理控制水口角度。根据数值模拟结果和实际生产经验，对于本研究中的硅钢薄带双辊铸轧过程，合适的水口角度范围为[\theta_{min},\theta_{max}]。在这个角度范围内，既能保证熔池内的液面波动较小，又能使侧封板受到的冲击压力在可承受范围内，从而确保铸轧过程的稳定进行，提高铸带质量。当水口角度低于\theta_{min}时，应适当增大水口角度，以减小液面波动和侧封板冲击；当水口角度高于\theta_{max}时，应降低水口角度，以避免铸带出现质量问题和生产不稳定的情况。五、双辊铸轧过程中的凝固组织模拟5.1凝固组织模拟模型5.1.1元胞自动机-有限元模型原理在双辊铸轧过程的凝固组织模拟中，基于ProCast的元胞自动机-有限元（CAFE）模型是一种强大的工具，它将元胞自动机（CA）的微观模拟能力与有限元（FE）的宏观分析优势相结合，能够深入揭示凝固过程中微观组织的演变规律。元胞自动机是一种离散的动态系统，它由规则排列的元胞组成，每个元胞具有有限的离散状态。在凝固组织模拟中，这些元胞可以代表材料中的微观区域，其状态可以表示为液态、固态或正在凝固的状态。元胞的状态更新遵循一定的局部规则，这些规则通常基于元胞自身及其邻居元胞的状态来确定。在凝固过程中，晶核的形成和生长可以通过元胞状态的转变来模拟。当某个元胞的温度降低到凝固点以下，且满足一定的形核条件时，该元胞的状态可以从液态转变为固态，形成晶核。随着时间的推移，晶核周围的液态元胞会根据生长规则逐渐转变为固态，从而实现晶体的生长。元胞自动机的并行性和自组织性使其能够有效地模拟微观组织的复杂演变过程，捕捉到微观尺度上的细节信息。有限元方法则是一种用于求解连续介质力学问题的数值方法。在双辊铸轧凝固组织模拟中，有限元方法主要用于求解宏观的温度场、流场和应力场等。通过将求解区域离散为有限个单元，有限元方法可以将复杂的连续介质问题转化为一系列简单的单元问题进行求解。在温度场计算中，有限元方法根据能量守恒原理，求解热传导方程，得到整个铸轧区域的温度分布。在流场计算中，通过求解动量守恒方程和连续性方程，得到金属熔体在铸轧过程中的流动速度和压力分布。这些宏观物理场的计算结果为元胞自动机模拟提供了重要的边界条件和初始条件。CAFE模型将元胞自动机和有限元方法有机结合。在模拟过程中，有限元方法首先计算宏观的温度场、流场等物理场，然后将这些物理场的结果传递给元胞自动机。元胞自动机根据接收到的物理场信息，按照自身的规则进行微观组织的模拟，包括晶核的形成、生长和晶粒的相互作用等。在晶核形成阶段，元胞自动机根据有限元计算得到的温度场，判断哪些元胞满足形核条件，从而形成晶核。在晶体生长阶段，元胞自动机根据温度梯度、溶质浓度等因素，确定晶体的生长方向和速度，实现晶粒的生长。元胞自动机模拟得到的微观组织演变结果，又可以反馈给有限元方法，用于更新宏观物理场的计算，考虑微观组织对宏观物理性质的影响。通过这种双向耦合的方式，CAFE模型能够更加准确地模拟双辊铸轧过程中的凝固组织演变，为研究铸轧过程中的微观组织控制提供了有力的手段。5.1.2模型参数确定在基于ProCast的元胞自动机-有限元模型中，准确确定模型参数是保证模拟结果准确性的关键。这些参数包括形核率、生长速度、界面能等，它们对凝固组织的模拟结果有着重要影响。形核率是指单位时间、单位体积内形成的晶核数量，它是凝固组织模拟中的一个重要参数。形核率的大小直接影响着晶核的数量和分布，进而影响着最终的凝固组织。在本研究中，采用经典的均匀形核理论和异质形核理论来计算形核率。对于均匀形核，根据热力学原理，形核率与过冷度、原子扩散系数等因素有关，其计算公式为：I_{hom}=N_0\nu\exp\left(-\frac{\DeltaG^{*}_{hom}}{kT}\right)其中，I_{hom}是均匀形核率，N_0是单位体积内的原子数，\nu是原子振动频率，\DeltaG^{*}_{hom}是均匀形核的临界自由能，k是玻尔兹曼常数，T是温度。对于异质形核，考虑到实际凝固过程中通常存在外来质点作为形核核心，其形核率与均匀形核率以及外来质点的性质、数量等因素有关，计算公式为：I_{het}=I_{hom}f(\theta)其中，I_{het}是异质形核率，f(\theta)是与接触角\theta有关的函数，反映了外来质点对形核的促进作用。接触角\theta可以通过实验测量或理论计算得到，它取决于外来质点与液态金属之间的界面能以及液态金属的表面能等因素。生长速度是指晶体在凝固过程中生长的速率，它也是影响凝固组织的重要参数。晶体的生长速度与温度梯度、溶质浓度梯度等因素密切相关。在本研究中，采用KGT模型来计算晶体的生长速度。KGT模型考虑了溶质扩散和界面动力学对晶体生长的影响，其生长速度的计算公式为：v=\frac{D}{r}\left(\frac{C_0-C_s}{C_l-C_s}\right)\left(1-\exp\left(-\frac{vr}{D}\right)\right)其中，v是晶体生长速度，D是溶质在液相中的扩散系数，r是晶体的曲率半径，C_0是合金的原始成分，C_s和C_l分别是固相和液相中的溶质浓度。界面能是指固液界面单位面积上的能量，它对晶核的形成和生长起着重要的作用。界面能的大小影响着晶核的临界尺寸和形核功，进而影响形核率。在本研究中，界面能通过查阅相关文献和实验数据来确定。对于硅钢材料，其固液界面能与温度、成分等因素有关。通过参考已有的研究成果，确定在本研究的温度和成分范围内，硅钢的固液界面能为[具体数值]。同时，考虑到界面能在凝固过程中的变化，采用温度相关的界面能模型，以更准确地反映其对凝固组织的影响。除了上述主要参数外，模型中还涉及到其他一些参数，如溶质分配系数、扩散系数等。这些参数同样对凝固组织模拟结果有着重要影响。溶质分配系数反映了溶质在固相和液相中的分配比例，它可以通过实验测量或热力学计算得到。扩散系数则决定了溶质在液相和固相中扩散的快慢，通常通过实验测定或理论模型计算得到。在本研究中，根据硅钢的成分和凝固条件，通过查阅相关文献和实验数据，确定了溶质分配系数和扩散系数的具体数值。通过合理确定形核率、生长速度