基于ANSYS的铸件充型凝固过程数值模拟：理论、方法与应用

基于ANSYS的铸件充型凝固过程数值模拟：理论、方法与应用一、绪论1.1研究背景与意义铸造作为一种重要的金属热加工工艺，在现代工业生产中占据着举足轻重的地位。它能够将液态金属转化为具有特定形状和性能的铸件，广泛应用于汽车、航空航天、机械制造、能源等众多领域，为各行业提供关键的零部件。例如，汽车发动机的缸体、缸盖，航空发动机的涡轮叶片等，这些高精度、复杂形状的铸件对于设备的性能和可靠性起着决定性作用。在铸造过程中，充型和凝固是两个最为关键的环节，直接关乎铸件的质量和性能。充型过程是指液态金属在重力、压力等外力作用下，填充铸型型腔的过程。若充型过程不平稳，液态金属可能会卷入气体，形成气孔缺陷；若充型时间过长或充填顺序不合理，还可能导致冷隔、浇不足等缺陷，使铸件无法达到设计要求，严重影响其使用性能。而凝固过程则是液态金属逐渐冷却转变为固态的过程，在此过程中，铸件的组织结构和性能得以形成。如果凝固过程控制不当，铸件内部可能会产生缩孔、缩松、热应力、裂纹等缺陷，降低铸件的强度、韧性和耐腐蚀性。传统的铸造工艺设计主要依赖于经验和反复的试验，这种方式不仅耗费大量的人力、物力和时间，而且难以准确预测铸件在充型和凝固过程中可能出现的缺陷，导致铸件质量不稳定，废品率较高。随着计算机技术和数值计算方法的飞速发展，数值模拟技术逐渐成为铸造工艺研究和优化的重要手段。通过数值模拟，可以在计算机上虚拟再现铸造过程中液态金属的流动、热量传递、应力应变等物理现象，预测铸件可能出现的缺陷，并通过优化工艺参数来减少或消除这些缺陷。这不仅可以大大缩短新产品的研发周期，降低生产成本，还能提高铸件的质量和性能，增强企业的市场竞争力。ANSYS作为一款功能强大的通用工程模拟软件，具有丰富的物理模型库、高效的数值求解器和友好的用户界面。它能够对铸造过程中的多物理场进行耦合分析，包括流场、温度场、应力场等，为铸造充型凝固过程的数值模拟提供了有力的工具。基于ANSYS进行铸件充型凝固过程的数值模拟，能够更加准确地揭示铸造过程中的物理本质，为铸造工艺的优化提供科学依据，对于推动铸造行业的技术进步和发展具有重要的现实意义。1.2国内外研究现状1.2.1充型过程数值模拟发展状况国外对于充型过程数值模拟的研究起步较早，在20世纪80年代初期，HWANGWS等学者通过将计算流体动力学的结果与浇铸和填充问题相结合，开启了浇注和充型过程数值模拟计算的研究先河。此后，越来越多的学者投身于这一领域，研究重点逐渐从二维过程转移到三维过程。例如，LINHJ等通过结合使用SOLA-VOF方法与传热学，实现了对填充缺陷的预测；PARADIESCJ等则通过N-S方程及非牛顿表观粘度模型，深入分析了半固态合金的流动特性。进入21世纪，相关研究变得更加深入、详细和准确。LEEDY模拟了半固态镁合金的填充过程，并深入分析了熔融金属粘度对充型过程的影响；MEBRIDED等开发了一种新的计算模型，成功解决了离心铸造方面的问题，能够在填充复杂铸件时精确捕获流体自由流动的表面细微之处；HUMPHREYSNJ通过对TiAl基合金离心铸造与重力铸造充型过程中自由表面的研究，为离心铸造充型过程的模拟与铸件缺陷预测提供了坚实的理论基础；CLEARYWP证明了SPH作为预测复杂薄壁模具填充行为计算工具的适用性；SARKARS等建立了三维单相和两相磁流体动力流动和湍流模型，用于表征叠加磁场影响下连铸模具中的单相和多相流体动力学。我国在充型过程数值模拟领域的研究起步虽晚，但发展迅速。1987年，沈阳铸造研究所的王君卿在研究中采用SIMPLE、SMAC、SOLA-VOF三种方法进行计算，并对几种求解带自由表面的数值计算方法进行了比较分析。此后，清华大学的荆涛、柳百成利用SOLA-VOF法对充型过程进行模拟，并研究了充型过程对浇注完成后铸型内初始温度场的影响。90年代中后期，华中理工大学的袁浩扬、陈立亮等人以SOLA-VOF法为基础，结合他们提出的三维自由表面边界速度确定方法，成功实现了铸造充型流动过程的三维数值模拟，并开发了适用于微机的充填过程模拟软件。近年来，国内学者不断在数值算法、模型改进等方面进行创新研究，如赵静雨等建立了基于SPH方法的铸造充型过程的三维数学模型，并通过水力学算例验证了模型的正确性；杨铭等针对熔模铸造蜡模充型过程数值模拟软件自主开发，采用Cross-WLF粘度模型描述蜡料的流变性能，基于ProjectionVOF方法进行蜡料流动过程速度场、压力场以及自由表面的求解。1.2.2凝固过程数值模拟发展概况凝固过程数值模拟的研究旨在预测铸件在凝固过程中的温度分布、相变过程以及由此产生的缩孔、缩松、热应力等缺陷。在国外，相关研究已经取得了丰硕的成果。早期主要集中在温度场的模拟，随着研究的深入，逐渐发展到对凝固过程中应力场、微观组织演变等多物理场的耦合模拟。例如，一些学者通过建立热弹塑性模型来计算凝固过程中的热应力，考虑了材料在高温下的力学性能变化以及相变对应力的影响。在微观组织模拟方面，相场法成为研究的热点，通过将相场方程与温度场、溶质场、流速场及其它外部场耦合，能够对金属液的凝固过程进行真实的模拟，如多个晶粒生长时多元相场的耦合、枝晶生长过程中相场与温度场或溶质场的耦合等。国内在凝固过程数值模拟领域也取得了显著的进展。大连理工大学较早地开展了铸件和连铸坯的应力数值模拟研究，采用自行开发有限元应力分析程序来处理研究连铸、半连铸及铸锭等能够简化为一维或二维的问题。近年来，国内学者在借鉴国外先进技术的基础上，不断进行创新和改进。例如，在缩松、缩孔预测方面，李嘉荣等提出了球墨铸铁缩松、缩孔形成预测的“收缩膨胀动态叠加法(DECAM)”，李文珍等从微观形核和生长的角度建立了球墨铸铁在凝固过程中的体积变化模型，并提出了基于微观模拟的球墨铸铁缩松、缩孔定量预测方法微观模拟法(MMM)。在应力场数值模拟方面，除了采用传统的有限元方法外，还尝试将FDM(有限差分法)和FEM(有限元法)方法相结合，利用FDM分析流动与传热，用FEM计算应力，以提高模拟的精度和效率。尽管凝固过程数值模拟取得了很大的进展，但目前仍然面临一些挑战。例如，在模拟复杂铸件的凝固过程时，如何准确考虑铸件的几何形状、材料特性、边界条件以及多物理场之间的相互作用，仍然是一个亟待解决的问题。此外，对于一些新型铸造材料和工艺，如铝合金半固态铸造、镁合金铸造等，现有的模拟方法和模型还需要进一步完善和优化。1.2.3目前铸造数值模拟存在问题及课题提出目前，铸造数值模拟虽然在充型和凝固过程的研究中取得了一定的成果，但仍然存在一些问题，限制了其在实际工程中的广泛应用。一方面，许多研究将充型过程和凝固过程分开进行模拟。只进行充型过程数值分析时，得到的温度计算结果往往不能直接用于预测缺陷，缺乏工程实用价值；而只进行凝固过程数值分析时，温度初始条件多建立在“瞬间充型、温度均布”的假设基础上，这与实际铸造过程存在较大差异，从而难以得到可靠的温度场、应力场计算结果。建立在这种不准确温度场结果基础上的缩孔、缩松、热裂等缺陷的判断也难以准确，尤其是对于尺寸厚大的零件，这种误差更为明显。另一方面，现有的铸造模拟软件大多是针对特定的铸造工艺或铸件类型开发的，通用性较差。不同软件之间的数据兼容性和互操作性也存在问题，导致在实际应用中难以实现多物理场的全面耦合分析和不同模拟环节之间的无缝衔接。此外，一些复杂的物理现象，如液态金属在充型过程中的紊流、凝固过程中的相变潜热释放、热应力与组织演变的耦合等，在数值模拟中还难以准确描述和模拟。针对以上问题，本课题提出基于ANSYS这一通用模拟软件平台进行铸造充型与凝固全过程数值模拟的研究。ANSYS具有强大的多物理场耦合分析能力和丰富的材料模型库，能够综合考虑铸造过程中各种复杂的物理现象和因素。通过深入研究ANSYS在铸造领域的应用，开发适用于铸造充型凝固过程的数值模拟方法和流程，可以提高铸造数值模拟的准确性和可靠性，为铸造工艺的优化和铸件质量的提升提供更加有效的技术支持，推动铸造数值模拟技术向工程实用化方向发展。1.3相关数值模拟方法1.3.1充型过程数值模拟方法充型过程数值模拟主要涉及流场模拟和自由表面处理两个关键方面。在流场模拟方面，常用的方法有有限差分法（FDM）、有限元法（FEM）和有限体积法（FVM）。有限差分法是将求解区域划分为差分网格，用差商代替微商，将控制方程转化为差分方程进行求解。该方法计算简单，程序实现容易，但对复杂几何形状的适应性较差。有限元法是将连续的求解区域离散为有限个单元的组合体，通过对单元内的未知函数进行插值逼近，将控制方程转化为代数方程组进行求解。有限元法对复杂几何形状具有良好的适应性，能够处理各种边界条件，但计算量较大，对计算机硬件要求较高。有限体积法是将计算区域划分为一系列不重复的控制体积，使每个网格点周围都有一个控制体积，将控制方程对控制体积进行积分，得到离散方程。有限体积法具有守恒性好、对复杂几何形状适应性强等优点，在充型过程流场模拟中得到了广泛应用。对于充型过程中的自由表面处理，常用的方法有VOF（VolumeofFluid）法、LevelSet法和SPH（SmoothedParticleHydrodynamics）法等。VOF法是通过引入体积函数来追踪自由表面，计算每个控制体积内流体的体积分数，根据体积分数来确定自由表面的位置。该方法计算效率较高，能够较好地处理大变形自由表面问题，但在自由表面曲率较大时可能会出现数值振荡。LevelSet法是通过定义一个符号距离函数来描述自由表面，该函数在自由表面上的值为零，在流体内部为正值，在流体外部为负值。通过求解符号距离函数的演化方程来追踪自由表面的运动。LevelSet法具有较高的精度和稳定性，能够准确处理自由表面的拓扑变化，但计算量较大，对初始条件的设置较为敏感。SPH法是一种无网格的拉格朗日方法，通过将流体离散为一系列相互作用的粒子，根据粒子的运动和相互作用来描述流体的流动。SPH法对复杂几何形状和自由表面的处理具有独特的优势，能够自然地处理自由表面的破碎和融合等现象，但计算精度相对较低，计算效率有待提高。不同的流场模拟方法和自由表面处理方法各有优缺点，在实际应用中需要根据具体问题的特点和要求选择合适的方法。例如，对于简单几何形状的铸件充型模拟，有限差分法或有限体积法结合VOF法可能是较为合适的选择；而对于复杂几何形状和自由表面拓扑变化较大的铸件充型模拟，有限元法结合LevelSet法或SPH法可能更能满足计算精度和可靠性的要求。1.3.2凝固过程数值模拟方法凝固过程数值模拟主要关注热应力计算和热力耦合分析。在热应力计算方面，常用的力学模型有纯弹性模型、弹塑性模型、弹塑性-蠕变模型、弹-粘塑性模型、统一变量模型和流变学模型等。纯弹性模型假设材料在整个凝固过程中始终处于弹性状态，不考虑材料的塑性变形和蠕变等行为，计算简单，但不能准确描述材料在高温下的力学性能变化。弹塑性模型考虑了材料的塑性变形，能够较好地模拟铸件在凝固过程中由于温度变化和收缩受阻而产生的塑性应变，但对于一些高温合金，蠕变行为对热应力的影响不可忽略，此时弹塑性模型的模拟精度会受到一定限制。弹塑性-蠕变模型综合考虑了材料的塑性变形和蠕变行为，能够更准确地预测铸件在高温下的热应力，但模型较为复杂，计算量较大。弹-粘塑性模型则将粘性和塑性相结合，用于描述材料在高温下的非线性力学行为。统一变量模型和流变学模型由于更接近于合金在铸造中的物理本质而逐渐成为研究的热点。流变学模型用简单的流体模型和力学模型来描述铸造合金在固液两相区的流动及变形规律，从而准确地反映流动变形随时间的变化。在热力耦合数值计算方面，常用的方法有顺序耦合和直接耦合。顺序耦合是先计算温度场，然后将温度场的结果作为载荷施加到应力场计算中，依次进行求解。这种方法计算过程相对简单，易于实现，但由于没有考虑温度场和应力场之间的实时相互作用，可能会导致一定的计算误差。直接耦合是将温度场和应力场的控制方程联立求解，同时考虑热传导、热对流、热辐射以及力学平衡等因素，能够更准确地模拟热力耦合过程，但计算量较大，对计算机硬件和求解器的要求较高。在实际应用中，需要根据铸件的材料特性、几何形状、凝固工艺以及对计算精度和效率的要求，选择合适的热应力计算模型和热力耦合数值计算方法。例如，对于一些对热应力要求不高的普通铸件，可以采用较为简单的弹塑性模型和顺序耦合方法；而对于航空航天等领域中对热应力和变形控制要求严格的高性能铸件，则需要采用更精确的流变学模型和直接耦合方法进行模拟分析。1.4ANSYS平台软件简介ANSYS软件是一款功能强大的通用工程模拟软件，广泛应用于机械、航空航天、汽车、电子、能源等众多领域。它提供了丰富的物理场分析模块，包括结构力学、流体力学、热分析、电磁学等，能够对各种复杂的工程问题进行多物理场耦合分析。在铸造数值模拟领域，ANSYS具有独特的优势。其流体分析模块ANSYSFluent可以对铸造充型过程中的液态金属流动进行精确模拟，通过求解Navier-Stokes方程，考虑流体的粘性、紊流等因素，能够准确预测液态金属的流速、压力分布以及自由表面的变化情况。同时，ANSYSFluent还支持多种湍流模型和多相流模型，可以模拟不同铸造工艺下的复杂流动现象。在热分析方面，ANSYS软件能够精确计算铸件在充型和凝固过程中的温度场分布，考虑热传导、热对流和热辐射等多种传热方式，以及凝固潜热的释放和吸收。通过与流体分析模块的耦合，可以实现铸造充型过程的热流耦合计算，更真实地反映实际铸造过程中的物理现象。ANSYS的结构力学分析模块ANSYSMechanical则可用于计算铸件在凝固过程中的热应力场。它提供了多种力学本构模型，如前面提到的弹塑性模型、弹塑性-蠕变模型等，可以根据铸件材料的特性和实际工况选择合适的模型进行热应力分析。通过将温度场计算结果作为载荷输入到结构力学分析模块中，能够准确预测铸件在凝固过程中由于温度变化和收缩受阻而产生的热应力和变形，为评估铸件的质量和可靠性提供重要依据。此外，ANSYS软件还具有友好的用户界面和强大的前后处理功能。用户可以方便地创建几何模型、划分网格、设置边界条件和材料参数等。在后处理方面，ANSYS提供了丰富的可视化工具，能够以直观的方式展示模拟结果，如温度场分布云图、应力场分布云图、流速矢量图等，帮助用户快速理解和分析模拟结果，从而进行铸造工艺的优化和改进。二、铸件充型过程流场与温度场耦合计算2.1控制方程2.1.1连续性方程连续性方程是描述流体流动过程中质量守恒的基本方程。在铸件充型过程中，液态金属作为流体，其质量在流动过程中不会凭空产生或消失，只能从一个区域转移到另一个区域，连续性方程正是对这一物理现象的数学表达。其一般形式为：\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0其中，\rho表示流体的密度，t为时间，\vec{v}是流体的速度矢量，\nabla\cdot(\rho\vec{v})是质量流量的散度，表示单位时间内通过单位体积表面的质量通量，\frac{\partial\rho}{\partialt}是密度随时间的变化率。在稳态流动情况下，流体的密度不随时间变化，即\frac{\partial\rho}{\partialt}=0，此时连续性方程简化为\nabla\cdot(\rho\vec{v})=0。对于不可压缩流体，其密度\rho为常数，连续性方程进一步简化为\nabla\cdot\vec{v}=0，这意味着速度场的散度为零，即流入某一控制体积的流体质量等于流出该控制体积的流体质量，反映了质量守恒这一基本物理原理。在铸件充型过程中，若液态金属可近似看作不可压缩流体，该简化形式的连续性方程将在数值模拟中发挥重要作用，用于确定液态金属在不同位置的流速关系，为后续分析充型过程的流动特性奠定基础。2.1.2动量守恒方程动量守恒方程在分析流体受力和运动状态方面具有重要应用，它是牛顿第二定律在流体力学中的具体体现，描述了流体动量的变化与作用在流体上的力之间的关系。对于牛顿流体，动量守恒方程的一般形式为Navier-Stokes方程：\rho\left(\frac{\partial\vec{v}}{\partialt}+\vec{v}\cdot\nabla\vec{v}\right)=-\nablap+\mu\nabla^2\vec{v}+\vec{f}其中，\rho为流体密度，t为时间，\vec{v}是流体速度矢量，p为压力，\mu为动态粘度，\vec{f}为作用在流体上的体积力（如重力、电磁力等）。方程左边\rho\left(\frac{\partial\vec{v}}{\partialt}+\vec{v}\cdot\nabla\vec{v}\right)表示单位体积流体的动量变化率，其中\frac{\partial\vec{v}}{\partialt}是速度对时间的偏导数，表示非定常项，反映了流体速度随时间的变化；\vec{v}\cdot\nabla\vec{v}是对流项，表示由于流体的流动而引起的动量变化。方程右边-\nablap表示压力梯度力，压力的变化会促使流体产生运动；\mu\nabla^2\vec{v}是粘性力项，体现了流体内部粘性对动量传递的影响，粘性使得相邻流体层之间产生剪切应力，阻碍流体的相对运动；\vec{f}则表示体积力，例如在重力场中，\vec{f}=\rho\vec{g}，\vec{g}为重力加速度。在铸件充型过程中，液态金属在铸型型腔中的流动受到重力、压力以及铸型壁面的摩擦力等多种力的作用。通过求解动量守恒方程，可以得到液态金属在不同时刻、不同位置的速度和压力分布，进而分析液态金属的充型路径、流速变化以及可能出现的流动不稳定现象，如紊流、漩涡等，这些信息对于预测铸件的质量缺陷（如气孔、冷隔等）具有重要意义。该方程的推导基于牛顿第二定律，即力等于质量乘以加速度。在流体动力学中，选择一个流体的微小控制体，对控制体应用牛顿第二定律，得到作用在控制体上的力等于控制体动量的变化率。同时，考虑流体的连续性来描述流体的质量守恒，在高速流动中，还需考虑流体的可压缩性，即流体密度的变化。将这些因素整合到流体动力学的基本方程中，最终得到动量守恒方程。通过这样的推导过程，从基本物理原理出发建立了描述流体运动的数学模型，为深入研究铸件充型过程中的流体动力学行为提供了理论依据。2.1.3能量方程能量方程在研究流体能量传递和转化中具有至关重要的作用，它基于热力学第一定律，即能量守恒定律，描述了流体在流动过程中各种能量形式的相互转化和传递关系。对于非绝热、不可压缩流动，能量方程的一般形式为：\rhoc_p\left(\frac{\partialT}{\partialt}+\vec{v}\cdot\nablaT\right)=\nabla\cdot(k\nablaT)+\Phi其中，\rho为流体密度，c_p是定压比热容，T是温度，t为时间，\vec{v}是流体速度矢量，k是热传导系数，\Phi是由于粘性耗散产生的热。方程左边\rhoc_p\left(\frac{\partialT}{\partialt}+\vec{v}\cdot\nablaT\right)表示单位体积流体的内能变化率，其中\frac{\partialT}{\partialt}是温度对时间的偏导数，表示非定常项，反映了温度随时间的变化；\vec{v}\cdot\nablaT是对流项，表示由于流体的流动而引起的温度变化，即热对流。方程右边\nabla\cdot(k\nablaT)表示热传导项，体现了由于温度梯度引起的热量传递，热量会从高温区域向低温区域传导；\Phi表示粘性耗散产生的热，当流体流动时，粘性力会做功，将机械能转化为热能，从而导致流体温度升高。在铸件充型过程中，液态金属在填充铸型型腔的同时，会与铸型壁面以及周围环境进行热量交换，并且由于液态金属内部的粘性作用，也会产生能量耗散。通过求解能量方程，可以得到液态金属在充型过程中的温度分布以及温度随时间的变化情况。这对于预测铸件的凝固过程、凝固时间以及可能出现的热应力、缩孔缩松等缺陷具有重要意义。例如，如果充型过程中液态金属温度分布不均匀，可能导致铸件不同部位的凝固速度不一致，从而产生热应力，当热应力超过材料的强度极限时，就会引发裂纹等缺陷。因此，准确掌握铸件充型过程中的温度场变化，是优化铸造工艺、提高铸件质量的关键。2.1.4紊流模型及方程在铸件充型过程中，紊流现象较为常见，它对充型过程有着重要影响。紊流中流体的运动呈现出不规则的脉动特性，各层流体之间相互混掺，流体质点作不规则运动。这种不规则运动使得流体的动量、热量和质量传递过程变得更加复杂，增加了充型过程的不确定性。与层流相比，紊流具有更强的扩散能力，能够加快液态金属在铸型型腔中的填充速度，但同时也容易导致液态金属卷入气体，形成气孔等缺陷。此外，紊流还会影响液态金属的温度分布和凝固过程，进而对铸件的质量产生影响。为了准确模拟紊流对充型过程的影响，需要引入合适的紊流模型。常用的紊流模型有雷诺平均纳维-斯托克斯方程（RANS）模型、大涡模拟（LES）模型和直接数值模拟（DNS）模型等。DNS模型直接求解Navier-Stokes方程，不进行任何平均化处理，能够精确地描述紊流的所有细节，从理论上讲，它是最准确的紊流模拟方法。然而，由于紊流运动包含的尺度范围非常广泛，从大尺度的涡旋到极小尺度的耗散涡，要精确解析所有尺度的运动，需要极其精细的网格和极小的时间步长，这对计算机的存储能力和运算速度提出了极高的要求，目前的计算机硬件水平还难以满足实际工程中对DNS模型的计算需求。LES模型则是一种折中的方法，它通过对大尺度涡旋进行直接模拟，而对小尺度涡旋采用亚网格尺度模型进行模拟。大尺度涡旋包含了紊流运动的主要能量和动量传递，对流动的整体特性起着关键作用，通过直接模拟可以较为准确地捕捉其运动特征。小尺度涡旋的运动相对较为均匀和各向同性，对大尺度运动的影响可以通过亚网格尺度模型进行近似描述。LES模型在一定程度上降低了对计算机资源的需求，同时又能够获得比RANS模型更详细的紊流信息，在一些对紊流细节要求较高的研究中得到了应用。RANS模型是目前工程中应用最为广泛的紊流模型。它将Navier-Stokes方程对时间进行平均，得到雷诺时均方程，通过引入一些假定和模型来封闭时均方程组，从而求解时均流速、时均压力等物理量。根据对雷诺应力的处理方式不同，RANS模型又分为基于涡粘性假设的模型和应力输运模型两类。基于涡粘性假设的模型引入一个涡粘度，将紊流中的雷诺应力与时均流速梯度建立联系，通过求解涡粘度来确定雷诺应力，这类模型包括零方程模型（如混合长模型）、单方程模型（如k方程模型）和双方程模型（如k-ε模型、k-ω模型）等。应力输运模型则是通过建立雷诺应力的输运方程来直接求解雷诺应力，理论上更加严谨，但计算复杂度也相对较高。在本研究中，考虑到计算效率和对铸件充型过程模拟精度的要求，选择k-ε双方程模型作为紊流模型。k-ε模型引入了紊动能k和紊动能耗散率ε两个变量，通过求解这两个变量的输运方程来确定涡粘度，进而封闭雷诺时均方程。其紊动能k的方程为：\frac{\partial(\rhok)}{\partialt}+\frac{\partial(\rhoku_j)}{\partialx_j}=\frac{\partial}{\partialx_j}\left(\frac{\mu_t}{\sigma_k}\frac{\partialk}{\partialx_j}\right)+G_k-\rho\varepsilon其中，u_j是速度分量，x_j是坐标分量，\mu_t是涡粘度，\sigma_k是紊动能k对应的普朗特数，G_k是由平均速度梯度产生的紊动能，\rho\varepsilon是紊动能的耗散项。紊动能耗散率ε的方程为：\frac{\partial(\rho\varepsilon)}{\partialt}+\frac{\partial(\rho\varepsilonu_j)}{\partialx_j}=\frac{\partial}{\partialx_j}\left(\frac{\mu_t}{\sigma_{\varepsilon}}\frac{\partial\varepsilon}{\partialx_j}\right)+C_{1\varepsilon}\frac{\varepsilon}{k}G_k-C_{2\varepsilon}\rho\frac{\varepsilon^2}{k}其中，\sigma_{\varepsilon}是紊动能耗散率ε对应的普朗特数，C_{1\varepsilon}和C_{2\varepsilon}是经验常数。选择k-ε模型的依据主要有以下几点：首先，k-ε模型在工程应用中经过了大量的验证，对于各种常见的紊流流动具有较好的模拟精度，能够满足对铸件充型过程模拟的基本需求。其次，该模型计算相对简单，对计算机资源的要求较低，在保证一定模拟精度的前提下，可以大大提高计算效率，适合处理复杂的铸件几何形状和大规模的计算区域。此外，k-ε模型的参数具有明确的物理意义，并且通过大量实验数据得到了较为准确的经验常数，使得模型的可靠性和可重复性较高。因此，综合考虑计算效率、模拟精度和模型的可靠性，k-ε模型是本研究中模拟铸件充型过程紊流现象的合适选择。2.2控制方程的有限体积法离散形式2.2.1流场模拟的各种数值方法比较在流场模拟中，常用的数值方法有有限差分法（FDM）、有限元法（FEM）和有限体积法（FVM），它们在原理、应用特点和适用场景等方面存在差异。有限差分法是将求解区域划分为差分网格，用差商代替微商，将控制方程转化为差分方程进行求解。其优点是计算原理简单直观，程序实现相对容易，对于规则的计算网格和简单的几何形状，能够快速得到计算结果。然而，有限差分法对复杂几何形状的适应性较差，当计算区域存在不规则边界时，网格划分较为困难，且边界条件的处理也相对复杂，容易引入数值误差。此外，有限差分法在处理守恒性问题时表现不佳，对于一些需要严格满足质量、动量和能量守恒的物理问题，其计算结果可能存在一定偏差。有限元法是将连续的求解区域离散为有限个单元的组合体，通过对单元内的未知函数进行插值逼近，将控制方程转化为代数方程组进行求解。有限元法对复杂几何形状具有良好的适应性，能够方便地处理各种不规则边界条件，可以精确地模拟具有复杂外形的物体周围的流场。它在处理流体与结构相互作用等问题时具有独特的优势，能够考虑结构变形对流体流动的影响。但是，有限元法的计算量通常较大，尤其是在处理大规模问题时，生成的代数方程组规模庞大，求解过程需要消耗大量的计算机内存和计算时间，对计算机硬件要求较高。此外，有限元法基于变分原理，在保守性方面可能不如有限体积法，对于一些对守恒性要求严格的流场模拟问题，其计算精度可能受到一定影响。有限体积法是将计算区域划分为一系列不重复的控制体积，使每个网格点周围都有一个控制体积，将控制方程对控制体积进行积分，得到离散方程。有限体积法的显著优点是具有严格的守恒性，它直接基于守恒定律，能够准确地保证质量、动量和能量在每个控制体积内的守恒，这对于流场模拟至关重要，能够确保模拟结果在物理上的合理性。同时，有限体积法对复杂几何形状也具有较强的适应性，可以灵活地处理各种不规则边界条件，通过合理选择控制体积的形状和积分格式，能够在保证计算精度的前提下，有效地提高计算效率。此外，有限体积法生成的离散方程形式简单，易于编程实现，并且在湍流模拟方面表现出色，能够较好地处理雷诺平均纳维-斯托克斯方程（RANS），因此在流场模拟中得到了广泛应用。在铸件充型过程的流场模拟中，由于铸型型腔的几何形状通常较为复杂，且需要精确地模拟液态金属的流动过程，保证质量和动量的守恒，因此有限体积法相较于有限差分法和有限元法具有明显的优势。它能够在处理复杂几何形状的同时，准确地捕捉液态金属的流动特性，为后续分析充型过程中的各种物理现象提供可靠的数值解。2.2.2计算区域离散将计算区域划分为离散的控制体积是有限体积法的关键步骤。在对铸件充型过程进行数值模拟时，首先需要根据铸型型腔的几何形状和尺寸，利用专业的网格划分软件（如ANSYSICEMCFD等）对计算区域进行离散。划分的原则主要包括以下几点：一是准确性原则。为了准确地模拟液态金属在铸型型腔中的流动，控制体积的划分应能够精确地描述铸型型腔的几何形状，尤其是在一些关键部位，如薄壁处、拐角处等，需要加密网格，以提高计算精度，确保能够捕捉到液态金属在这些区域的流动细节和物理现象。二是一致性原则。控制体积的大小和形状应尽量保持一致，避免出现过大或过小的控制体积，以减少数值计算中的误差。同时，在不同区域之间的网格过渡应尽量平滑，避免出现网格突变，保证计算结果的稳定性和可靠性。三是高效性原则。在保证计算精度的前提下，应尽量减少控制体积的数量，以提高计算效率。可以根据液态金属的流动特性和速度分布，在流速变化较大的区域适当加密网格，而在流速变化较小的区域适当稀疏网格，实现计算资源的合理分配。划分的方法通常有结构化网格划分和非结构化网格划分两种。结构化网格是指网格节点在空间上呈规则排列，网格线之间具有明确的拓扑关系，如四边形、六面体网格等。结构化网格划分简单，计算效率高，数据存储和处理方便，但其对复杂几何形状的适应性较差，对于一些不规则的铸型型腔，可能需要进行复杂的网格拼接和处理。非结构化网格则是指网格节点和单元的排列没有固定的规则，单元形状可以是三角形、四面体等。非结构化网格对复杂几何形状具有良好的适应性，能够根据铸型型腔的形状进行灵活划分，无需进行复杂的网格拼接，但非结构化网格的生成算法相对复杂，计算量较大，数据存储和处理也相对困难。在实际应用中，常常根据铸型型腔的具体情况，综合采用结构化网格和非结构化网格划分方法，如在铸型型腔的主体部分采用结构化网格，以提高计算效率；在一些复杂的局部区域，如薄壁、拐角等，采用非结构化网格，以准确描述几何形状。通过这样的混合网格划分策略，可以在保证计算精度的同时，提高计算效率，满足铸件充型过程数值模拟的需求。2.2.3动量方程的离散以笛卡尔坐标系下的二维动量方程为例，对其进行离散推导。二维动量方程三、凝固过程数值分析3.1凝固过程温度场数学模型及定解条件3.1.1数学模型在铸件凝固过程中，温度场的变化是一个复杂的传热过程，涉及热传导、对流和相变潜热等多种因素。热传导是热量在物体内部从高温区域向低温区域传递的过程，遵循傅里叶定律，其表达式为\vec{q}=-k\nablaT，其中\vec{q}是热流密度矢量，k是热导率，\nablaT是温度梯度。对流则是由于流体的宏观运动而引起的热量传递，在铸件凝固过程中，液态金属的流动会导致热量的对流传递。相变潜热是指金属在凝固过程中，从液态转变为固态时释放的热量，这部分热量对温度场的分布有着重要影响。考虑这些因素，建立凝固过程温度场的数学模型，其通用形式为：\rhoc\frac{\partialT}{\partialt}=\nabla\cdot(k\nablaT)+Q+\rhocv\cdot\nablaT其中，\rho是材料密度，c是比热容，T是温度，t是时间，k是热导率，Q是单位体积内的热源项（主要用于考虑相变潜热等因素，当考虑相变潜热时，Q=-L\frac{\partialf_s}{\partialt}，L是相变潜热，f_s是固相分数），v是流体速度矢量（在凝固过程中，若存在液态金属的流动，该项体现对流对温度场的影响）。此模型综合考虑了热传导、对流和相变潜热等因素，全面地描述了铸件凝固过程中温度场的变化。热传导项\nabla\cdot(k\nablaT)体现了热量在材料内部的传导，它反映了材料内部不同位置之间的温度差异导致的热量传递；对流项\rhocv\cdot\nablaT考虑了液态金属流动对热量传递的影响，当液态金属流动时，会携带热量，从而改变温度场的分布；热源项Q则主要用于考虑相变潜热的影响，在凝固过程中，相变潜热的释放会使温度场的变化变得更加复杂。通过这个数学模型，可以准确地模拟铸件凝固过程中温度场的动态变化，为后续分析凝固过程中的各种物理现象提供基础。3.1.2定解条件定解条件是求解数学模型的关键，它包括初始条件和边界条件，这些条件对于准确模拟铸件凝固过程中的温度场至关重要。初始条件定义了铸件在凝固开始瞬间的温度分布状态，即t=0时，T(x,y,z,0)=T_0(x,y,z)，其中T_0(x,y,z)表示铸件在初始时刻的温度分布。准确设定初始条件对于模拟结果的准确性有着重要影响，因为初始温度分布会直接影响后续凝固过程中温度场的变化趋势。例如，如果初始温度分布不均匀，那么在凝固过程中，温度较高的区域会先开始凝固，而温度较低的区域则会后凝固，这将导致铸件内部的凝固过程不一致，进而影响铸件的质量。边界条件则描述了铸件与周围环境之间的热交换情况，常见的边界条件有三类：第一类边界条件是给定物体表面温度随时间的变化关系，即T(x_s,y_s,z_s,t)=T_s(x_s,y_s,z_s,t)，其中(x_s,y_s,z_s)是物体表面的坐标，T_s(x_s,y_s,z_s,t)是给定的表面温度。在铸件凝固过程中，当铸件表面与恒温热源或冷源接触时，可采用第一类边界条件。例如，在一些特殊的铸造工艺中，会对铸型进行预热或冷却处理，此时铸件表面的温度就可以看作是一个已知的随时间变化的函数。第二类边界条件是给出通过物体表面的比热流随时间的变化关系，即\vec{q}\cdot\vec{n}=q_s(x_s,y_s,z_s,t)，其中\vec{n}是表面的外法线方向，q_s(x_s,y_s,z_s,t)是给定的比热流。当铸件表面存在特定的热流输入或输出时，如铸件表面受到辐射加热或冷却，就可以使用第二类边界条件。在实际铸造过程中，铸件表面与周围环境之间的辐射换热就可以通过这种方式来描述。第三类边界条件是给出物体周围介质温度以及物体表面与周围介质的换热系数，即-k\frac{\partialT}{\partialn}=h(T-T_f)，其中h是换热系数，T_f是周围介质温度。在大多数铸造情况下，铸件表面与周围空气或铸型之间的热交换都可以用第三类边界条件来描述。例如，铸件在空气中自然冷却时，空气的温度和铸件表面与空气之间的换热系数是已知的，通过这个边界条件可以准确地模拟铸件表面的散热过程。这些边界条件对模拟结果有着显著的影响。不同的边界条件会导致铸件表面的热量传递方式和速率不同，从而影响整个铸件的温度分布和凝固进程。如果边界条件设置不合理，可能会导致模拟结果与实际情况相差较大，无法准确预测铸件在凝固过程中可能出现的缺陷，如缩孔、缩松等。因此，在进行数值模拟时，需要根据实际的铸造工艺和环境条件，准确地确定边界条件，以确保模拟结果的可靠性。3.1.3凝固潜热的处理在铸件凝固过程中，凝固潜热的释放是一个重要的物理现象，它对温度场的分布和凝固进程有着显著影响。为了准确模拟凝固过程，需要合理处理凝固潜热。常见的处理方法有焓法、温度回升法等，它们各有优缺点。焓法是基于热焓的概念来处理凝固潜热。热焓H与温度T和固相分数f_s相关，其表达式为H=\int_{T_0}^T\rhoc(T)dT+f_sL，其中T_0是参考温度，\rho是密度，c(T)是温度相关的比热容，L是相变潜热。在数值计算中，通过求解热焓的变化来考虑凝固潜热的影响。焓法的优点是物理概念清晰，能够自然地将凝固潜热纳入能量方程的求解过程中，便于进行数值计算。它适用于各种凝固过程的模拟，尤其是对于结晶温度区间较宽的合金，能够较为准确地描述凝固过程中的温度变化。焓法也存在一些缺点，它需要预先知道材料的热焓与温度、固相分数之间的关系，这在实际应用中可能需要通过实验或其他方法来确定，增加了计算的复杂性。而且，在数值计算过程中，由于热焓的计算涉及积分运算，可能会引入一定的数值误差。温度回升法是将合金凝固时的潜热释放量折算成合金材质热容量所对应的温度变化量。其基本原理是，当合金凝固释放潜热时，假设这部分潜热全部用于使合金自身温度回升，通过计算温度的回升值来考虑凝固潜热的影响。对于结晶温度区间很窄的合金，温度回升法具有物理意义明确、计算简单的优点。在数值计算中，不需要复杂的积分运算，能够快速得到计算结果。但是，温度回升法也有局限性，它只适用于结晶温度区间较窄的合金，对于结晶温度区间较宽的合金，由于凝固过程中温度变化较为复杂，该方法难以准确描述凝固潜热的释放和温度场的变化，可能会导致较大的计算误差。在实际应用中，应根据合金的特性和计算精度的要求选择合适的凝固潜热处理方法。对于结晶温度区间较宽的合金，焓法通常是更好的选择，能够更准确地模拟凝固过程；而对于结晶温度区间较窄的合金，温度回升法因其计算简单、物理意义明确，可以在保证一定精度的前提下提高计算效率。在一些复杂的铸造过程中，也可以考虑将多种方法结合使用，以充分发挥各自的优势，提高模拟结果的准确性。3.1.4瞬态导热控制方程的有限元解法瞬态导热控制方程的有限元解法是求解铸件凝固过程温度场的重要方法之一，它通过将连续的求解区域离散化为有限个单元，将控制方程转化为代数方程组进行求解，从而得到温度场的数值解。在进行有限元求解时，首先需要对计算区域进行单元划分。根据铸件的几何形状和尺寸，利用专业的网格划分软件（如ANSYSMeshing等）将铸件和铸型划分为有限个单元，常用的单元类型有四面体单元、六面体单元等。划分单元时，需要遵循一定的原则，以保证计算精度和效率。在铸件的关键部位，如薄壁处、拐角处等，由于温度变化梯度较大，需要加密网格，以准确捕捉温度场的变化；而在温度变化较为平缓的区域，可以适当稀疏网格，以减少计算量。同时，单元的形状和大小应尽量均匀，避免出现过大或过小的单元，以免影响计算结果的准确性。形函数构造是有限元解法的关键步骤之一。对于每个单元，通过构造合适的形函数来近似表示单元内的温度分布。形函数是一种插值函数，它根据单元节点的温度值来确定单元内任意点的温度。不同的单元类型具有不同的形函数，例如，对于线性四面体单元，通常采用线性插值形函数；对于高阶单元，则需要采用更复杂的插值函数。形函数的构造应满足一定的条件，如在节点处取值为1，在其他节点处取值为0，以保证温度场的连续性和协调性。在完成单元划分和形函数构造后，将瞬态导热控制方程在每个单元上进行离散化处理。通过加权余量法或变分原理，将控制方程转化为单元节点温度的代数方程组。以加权余量法为例，假设控制方程为L(T)=0，其中L是微分算子，T是温度。选择一组权函数W_i，对控制方程在每个单元上进行积分，得到\int_{\Omega_e}W_iL(T)d\Omega=0，其中\Omega_e是单元的体积。将单元内的温度T用形函数和节点温度表示为T=\sum_{j=1}^nN_jT_j，其中N_j是形函数，T_j是节点温度，n是单元节点数。将其代入积分式中，经过一系列的数学推导和运算，得到关于节点温度T_j的代数方程组。对所有单元的代数方程组进行组装，形成全局的代数方程组。通过求解全局代数方程组，得到每个节点在不同时刻的温度值，从而得到整个计算区域的温度场分布。在求解过程中，常用的方法有直接解法和迭代解法。直接解法如高斯消去法、LU分解法等，适用于规模较小的方程组；对于大规模的方程组，迭代解法如雅可比迭代法、高斯-赛德尔迭代法、共轭梯度法等更为有效。在实际应用中，需要根据方程组的规模和特点选择合适的求解方法，以提高计算效率和精度。3.1.5三维瞬态温度场有限元法/有限体积法计算流程三维瞬态温度场的计算方法主要有有限元法和有限体积法，它们的计算流程各有特点，下面分别进行介绍并对比。有限元法的计算流程如下：几何建模：使用专业的三维建模软件（如SolidWorks、Pro/E等）或ANSYS自带的建模模块，根据铸件和铸型的实际尺寸和形状，创建精确的三维几何模型。在建模过程中，需要准确描述铸件的复杂结构和细节特征，以及铸型的形状和尺寸，确保模型能够真实反映实际的铸造情况。网格划分：将创建好的三维几何模型导入到网格划分软件（如ANSYSMeshing）中，根据铸件和铸型的几何形状和温度变化特点，选择合适的网格类型和尺寸进行划分。对于复杂的铸件，可能需要采用混合网格划分技术，在关键部位使用细密的网格，以提高计算精度；在其他部位使用相对稀疏的网格，以减少计算量。划分后的网格应保证质量良好，避免出现畸形单元，确保计算结果的准确性。材料参数定义：根据铸件和铸型所使用的材料，在ANSYS中定义相应的材料参数，包括密度、比热容、热导率等。这些材料参数通常是温度的函数，需要准确输入材料在不同温度下的参数值，以反映材料性能随温度的变化。对于一些特殊材料，还需要考虑其热膨胀系数、相变潜热等参数。边界条件设定：根据实际的铸造工艺和环境条件，在模型的边界上设定合适的边界条件。如前所述，常见的边界条件有第一类边界条件（给定表面温度）、第二类边界条件（给定表面热流密度）和第三类边界条件（给定表面与周围介质的换热系数和介质温度）。准确设定边界条件对于模拟结果的准确性至关重要，需要根据实际情况进行合理的选择和设置。初始条件设定：确定铸件在初始时刻的温度分布，通常假设铸件在初始时刻温度均匀，将初始温度值输入到ANSYS中作为初始条件。对于一些特殊的铸造工艺，如预热铸造，需要根据实际的预热情况设定初始温度分布。求解计算：在ANSYS中选择合适的求解器和求解设置，启动计算。求解器会根据设定的边界条件、初始条件和材料参数，对瞬态导热控制方程进行求解，得到每个时间步长下铸件和铸型内各节点的温度值。在计算过程中，需要密切关注计算的收敛情况，确保计算结果的可靠性。结果后处理：计算完成后，使用ANSYS的后处理模块对计算结果进行分析和可视化处理。可以绘制温度场分布云图、温度随时间变化曲线等，直观地展示铸件在凝固过程中的温度变化情况。通过后处理，可以提取关键部位的温度数据，分析温度场的分布规律和变化趋势，为评估铸件的质量和优化铸造工艺提供依据。有限体积法的计算流程如下：几何建模与网格划分：与有限元法类似，首先创建铸件和铸型的三维几何模型，并进行网格划分。不同的是，有限体积法将计算区域划分为一系列不重叠的控制体积，每个控制体积围绕一个网格节点。网格划分时，需要保证控制体积的形状和大小合理，以准确地离散控制方程。控制方程离散：将三维瞬态导热控制方程在每个控制体积上进行积分，通过对控制体积内的温度、热流等物理量进行离散化处理，得到关于控制体积节点温度的离散方程。在离散过程中，需要采用合适的插值方法和数值格式，以保证离散方程的准确性和稳定性。材料参数与边界条件处理：定义材料参数和设定边界条件的过程与有限元法相似，但在有限体积法中，边界条件的处理方式略有不同。对于不同类型的边界条件，需要根据控制体积的特点进行相应的处理，以确保边界条件能够准确地反映在离散方程中。初始条件设定：与有限元法相同，设定铸件在初始时刻的温度分布作为初始条件。迭代求解：由于离散方程通常是非线性的，需要采用迭代方法进行求解。常用的迭代方法有SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations）及其改进算法等。在迭代过程中，不断更新控制体积节点的温度值，直到满足收敛条件为止。结果处理：计算收敛后，对计算结果进行处理和分析。可以通过绘制温度场分布云图、温度随时间变化曲线等方式，直观地展示铸件在凝固过程中的温度变化情况。与有限元法类似，有限体积法的结果后处理也可以提取关键部位的温度数据，为后续的分析和决策提供支持。两种方法的计算过程和特点对比如下：计算过程：有限元法基于变分原理，通过对单元内的温度进行插值逼近，将控制方程转化为代数方程组进行求解；而有限体积法基于守恒原理，通过对控制体积内的物理量进行积分，得到离散方程进行求解。有限元法在处理复杂几何形状时具有优势，能够通过灵活的网格划分适应各种形状的铸件；有限体积法在处理守恒性问题上表现出色，能够准确地保证热量在每个控制体积内的守恒。特点：有限元法的计算精度较高，尤其是在处理复杂边界条件和材料特性时，能够通过合理的网格划分和形函数构造得到较为准确的结果；但其计算量较大，对计算机硬件要求较高，且在处理大变形和移动边界问题时相对复杂。有限体积法的计算效率较高，离散方程形式简单，易于编程实现，并且在处理流动和传热问题时具有较好的稳定性；但在处理复杂几何形状时，网格划分的灵活性相对较差，对于一些复杂的铸件结构，可能需要进行复杂的网格处理。在实际应用中，应根据铸件的具体特点、计算精度要求和计算机硬件条件等因素，选择合适的计算方法。对于形状复杂、对计算精度要求较高的铸件，有限元法可能更为合适；而对于一些对计算效率要求较高、形状相对简单的铸件，有限体积法可能是更好的选择。在一些情况下，也可以将四、Benchmark标准试件铸造过程流场与温度场耦合计算4.1伯明翰大学Sirrel.B标准验证实验简介伯明翰大学Sirrel.B标准验证实验旨在为铸造过程数值模拟提供一个标准的实验参考，以评估不同数值模拟方法和软件的准确性和可靠性。该实验聚焦于铸造过程中的充型和凝固环节，通过精确测量和记录实验数据，为数值模拟结果的验证提供了关键依据。实验装置主要包括铸型、浇注系统和数据采集设备。铸型采用特定的材料和结构，以确保实验条件的一致性和可重复性。浇注系统设计合理，能够精确控制液态金属的浇注速度和流量，保证充型过程的稳定性。数据采集设备则用于实时监测液态金属在充型和凝固过程中的温度变化、流速分布等关键参数。实验过程严格按照既定的操作流程进行。首先，将铸型预热至一定温度，以减少液态金属与铸型之间的温差，避免因热冲击导致的缺陷。然后，将高温的液态金属以设定的浇注速度通过浇注系统注入铸型型腔中。在充型过程中，利用高速摄像机和温度传感器等设备，同步记录液态金属的流动形态和温度分布。待充型完成后，继续监测铸件在凝固过程中的温度变化，直至铸件完全凝固。实验数据的获取采用了多种先进的测量技术。温度数据通过高精度的热电偶传感器进行测量，这些传感器被精确地布置在铸型和铸件的关键位置，能够实时捕捉温度的动态变化。流速数据则通过粒子图像测速（PIV）技术获得，该技术利用激光照射液态金属中的示踪粒子，通过拍摄粒子的运动轨迹来计算液态金属的流速。在数据处理方面，首先对采集到的原始数据进行筛选和滤波，去除异常值和噪声干扰，以提高数据的准确性和可靠性。然后，运用统计学方法对处理后的数据进行分析，计算出各参数的平均值、标准差等统计特征，以便与数值模拟结果进行对比分析。例如，通过对不同位置的温度数据进行统计分析，可以得到铸件在充型和凝固过程中的平均温度、温度梯度等信息，从而更全面地了解铸造过程的热行为。4.2充型过程流场与温度场计算4.2.1计算模型在ANSYS中建立充型过程计算模型时，首先利用专业的三维建模软件（如SolidWorks、Pro/E等），根据伯明翰大学Sirrel.B标准验证实验中铸型和浇注系统的实际尺寸和形状，创建精确的三维几何模型。在建模过程中，充分考虑铸型的结构细节、浇注系统的管道尺寸和弯曲角度等因素，确保模型能够真实反映实验的几何特征。将创建好的三维几何模型导入到ANSYSICEMCFD中进行网格划分。采用非结构化四面体网格对计算区域进行离散，在铸型型腔、浇注系统等关键部位，根据液态金属的流动特性和温度变化情况，适当加密网格，以提高计算精度，准确捕捉液态金属在这些区域的流动细节和温度变化。例如，在浇注系统的入口和出口处，以及铸型型腔的薄壁部位，网格尺寸设置得较小，以更好地模拟液态金属的高速流动和温度梯度变化；而在温度变化较为平缓的区域，网格尺寸则适当增大，以减少计算量，提高计算效率。在ANSYSFluent中进行材料参数设置。根据实验中所使用的液态金属和铸型材料，准确输入材料的密度、比热容、热导率、粘度等参数。这些材料参数通常是温度的函数，需要根据材料的热物性数据手册，输入不同温度下的参数值，以反映材料性能随温度的变化。例如，液态金属的粘度会随着温度的降低而增加，在设置材料参数时，需要考虑这一特性，确保模拟结果的准确性。4.2.2计算过程及结果对比在ANSYSFluent中，设置好计算模型和材料参数后，进行充型过程的计算。选择瞬态求解器，以模拟液态金属在充型过程中的非稳态流动和温度变化。设置合适的时间步长，根据液态金属的流速和计算区域的大小，通过多次试算，确定时间步长为[具体时间步长值]，以保证计算的稳定性和准确性。同时，开启能量方程和紊流模型（选择k-ε双方程模型），考虑液态金属在流动过程中的能量传递和紊流效应。在计算过程中，密切关注计算的收敛情况。通过监测残差曲线，确保各控制方程的残差在设定的收敛标准内。当所有控制方程的残差都收敛到足够小的值时，认为计算达到收敛，得到稳定的计算结果。将模拟结果与实验数据进行对比，分析模拟的准确性和误差来源。在流速分布方面，对比模拟得到的液态金属在浇注系统和铸型型腔中的流速与实验测量的流速。从对比结果来看，在浇注系统的直管段，模拟流速与实验流速较为接近，误差在[X]%以内，这表明模拟能够较好地捕捉液态金属在直管段的流动特性；然而，在浇注系统的弯管段和铸型型腔的复杂部位，模拟流速与实验流速存在一定偏差，误差在[X]%-[X]%之间，这可能是由于弯管段的流动较为复杂，存在二次流和漩涡等现象，而模拟中采用的紊流模型对这些复杂流动的描述存在一定局限性，以及网格划分在这些部位的精度可能不够，未能完全准确地捕捉到流动细节。在温度分布方面，对比模拟得到的铸件和铸型在充型过程中的温度与实验测量的温度。在充型初期，由于液态金属与铸型之间的热交换较为剧烈，模拟温度与实验温度的误差相对较大，在[X]℃左右，这可能是由于在模拟中对初始条件的设定存在一定误差，以及边界条件的处理不够精确，未能完全考虑到铸型与周围环境之间的辐射换热等因素；随着充型过程的进行，模拟温度与实验温度逐渐接近，在充型后期，误差减小到[X]℃以内，这说明模拟能够较好地反映液态金属在充型后期的温度变化趋势。模拟误差的来源主要包括以下几个方面：一是模型简化，在建立计算模型时，为了便于计算，对一些复杂的物理现象和几何细节进行了简化，这可能导致模拟结果与实际情况存在偏差；二是参数不确定性，材料参数的准确性对模拟结果有着重要影响，然而，实际材料的性能可能存在一定的波动，且一些参数的测量也存在误差，这会给模拟带来不确定性；三是数值算法的局限性，虽然有限体积法在流场和温度场模拟中具有较好的性能，但数值离散过程中仍然会引入一定的误差，尤其是在处理复杂流动和温度变化时，数值算法的局限性可能更加明显。4.3凝固过程温度场计算4.3.1边界条件及相变设置在凝固过程中，边界条件的设定对温度场计算有着重要影响。在ANSYSFluent中，将铸件与铸型的接触表面设置为第三类边界条件，即给定铸件表面与铸型之间的换热系数和铸型的初始温度。根据实验条件和相关文献资料，确定换热系数为[具体换热系数值]，铸型的初始温度为[具体温度值]。这是因为在实际铸造过程中，铸件与铸型之间主要通过对流和传导进行热交换，第三类边界条件能够较好地描述这种热交换情况。通过设置合适的换热系数，可以准确地模拟铸件在凝固过程中向铸型传递热量的过程，从而得到较为准确的温度场分布。对于铸件的外表面，考虑与周围空气的热交换，同样设置为第三类边界条件，给定铸件表面与空气之间的换热系数和空气的温度。根据实验环境和相关经验数据，确定铸件表面与空气之间的换热系数为[具体换热系数值]，空气温度为[具体温度值]。在凝固过程中，铸件外表面会向周围空气散热，通过设置这一边界条件，可以模拟铸件与空气之间的自然对流换热和辐射换热，使模拟结果更符合实际情况。相变设置是凝固过程模拟的关键环节。在ANSYSFluent中，采用焓-孔隙率法来处理凝固相变问题。该方法通过引入液相分数来描述材料的相变状态，当液相分数为1时，表示材料处于完全液态；当液相分数为0时，表示材料处于完全固态；在0-1之间时，表示材料处于固液两相区。在设置相变参数时，根据铸件材料的特性，确定相变温度范围为[具体相变温度范围]，相变潜热为[具体相变潜热值]。这些参数的准确设置对于模拟凝固过程中潜热的释放和吸收至关重要，直接影响温度场的计算结果。通过合理设置相变参数，可以准确地模拟铸件在凝固过程中由于相变潜热的作用而导致的温度变化，从而更真实地反映凝固过程的物理现象。4.3.2求解结果及对照经过在ANSYSFluent中的计算，得到了凝固过程温度场的求解结果。从模拟结果可以清晰地看到，在凝固初期，铸件的温度较高，随着时间的推移，热量逐渐从铸件表面向铸型和周围空气传递，铸件温度逐渐降低。在铸件的厚壁部位，由于热量积聚较多，凝固速度相对较慢，温度下降较为平缓；而在薄壁部位，热量散失较快，凝固速度相对较快，温度下降较为迅速。在固液两相区，由于相变潜热的释放，温度下降趋势减缓，出现了温度平台现象。将模拟结果与实验结果进行对照，分析模拟的可靠性。在凝固时间方面，模拟得到的铸件完全凝固时间为[具体时间值]，与实验测量的凝固时间[具体时间值]相比，误差在[X]%以内，这表明模拟能够较为准确地预测铸件的凝固时间，验证了模拟方法的可靠性。在温度分布方面，对比模拟得到的铸件不同位置在凝固过程中的温度与实验测量的温度。在铸件的关键部位，如厚壁与薄壁的交界处，模拟温度与实验温度的误差在[X]℃以内，能够较好地反映实际温度分布情况；然而，在一些局部区域，由于实验测量误差、材料特性的不均匀性以及模拟中对一些复杂物理现象的简化等因素，模拟温度与实验温度存在一定偏差，误差在[X]℃-[X]℃之间。总体而言，模拟结果与实验结果在趋势上基本一致，能够为铸造工艺的优化提供有价值的参考。4.4小结通过对Benchmark标准试件铸造过程流场与温度场耦合计算的研究，得出以下结论：在充型过程中，基于ANSYS建立的计算模型能够较好地模拟液态金属的流动和温度变化情况，模拟结果与实验数据在趋势上基本一致，但在一些复杂部位和初期阶段存在一定误差，主要原因包括模型简化、参数不确定性和数值算法的局限性等。在凝固过程中，合理设置边界条件和相变参数，采用焓-孔隙率法能够较为准确地模拟铸件的温度场分布和凝固时间，模拟结果与实验结果的对比验证了模拟方法的可靠性。综合来看，基于ANSYS的铸造过程流场与温度场耦合计算方法具有一定的准确性和有效性，能够为铸造工艺的研究和优化提供重要的技术支持，但仍需要进一步改进和完善，以提高模拟的精度和可靠性。五、基于ANSYSMechanical铸造过程流场、温度场、应力场计算5.1充型过程流场与温度场耦合计算5.1.1计算模型及材料参数在ANSYSMechanical中，建立基于实际铸件和铸型的充型过程计算模型。首先，利用三维建模软件（如SolidWorks、Pro/E等），根据铸件的设计图纸和铸型的结构特点，精确创建三维几何模型，确保模型的尺寸和形状与实际情况相符。将创建好的几何模型导入ANSYSMechanical中，利用其自带的网格划分工具或专业的网格划分软件（如ANSYSICEMCFD）进行网格划分。考虑到充型过程中液态金属的流动特性和温度变化情况，采用非结构化四面体网格对计算区域进行离散，在铸件的关键部位（如薄壁、拐角等）以及铸型与液态金属接触的区域，适当加密网格，以提高计算精度，准确捕捉液态金属在这些区域的流动细节和温度变化。在材料参数设置方面，充分考虑材料的热物理性质和流变特性。对于液态金属，根据其成分和工艺要求，设置密度、比热容、热导率、粘度等参数。例如，对于铝合金，其密度约为2700kg/m³，比热容在900-1000J/(kg・K)之间，热导率在150-200W/(m・K)左右，粘度则与温度和成分密切相关，在充型过程中会随着温度的降低而逐渐增大。这些参数的准确设置对于模拟液态金属的流动和传热过程至关重要，直接影响模拟结果的准确性。同时，考虑到材料的流变特性，引入合适的流变模型，如幂律模型或Bingham模型，以描述液态金属在不同剪切速率下的粘度变化。对于铸型材料，设置其密度、比热容、热导率等参数，这些参数会影响铸型与液态金属之间的热交换，进而影响充型过程中的温度分布。例如，砂型铸型的密度通常在1500-1800kg/m³之间，比热容在800-1000J/(kg・K)左右，热导率在0.5-1.0W/(m・K)之间，通过准确设置这些参数，可以更真实地模拟铸型对液态金属的冷却作用。5.1.2边界条件与初始条件确定充型过程的边界条件和初始条件是进行数值模拟的关键步骤。边界条件包括速度入口、压力出口、初始温度等，这些条件的设置依据实际的铸造工艺和物理过程。在浇口处设置速度入口边界条件，根据浇注系统的设计和实际浇注速度，确定液态金属的入口速度。例如，在重力浇注过程中，液态金属的入口速度可以根据伯努利方程计算得出，公式为v=\sqrt{2gh}，其中v是入口速度，g是重力加速度，h是浇口与液态金属液面的高度差。在铸型的排气孔或出口处设置压力出口边界条件，假设出口处的压力为大气压力，以模拟液态金属在充型过程中的压力变化和流动状态。在初始条件设置方面，假设铸件和铸型在充型开始时的温度均匀分布，将液态金属的浇注温度作为初始温度，铸型的初始温度根据实际铸造工艺确定。例如，在砂型铸造中，铸型通常在浇注前预热到一定温度，以减少液态金属与铸型之间的温差，避免因热冲击导致铸件产生缺陷。将预热后的铸型温度作为初始温度输入到模拟模型中，同时将液态金属的浇注温度也准确设置，以确保初始条件的准确性，为后续的模拟计算提供可靠的基础。5.1.3充型过程求解结果及分析通过在ANSYSMechanical中对充型过程进行求解计算，得到了流场与温度场的求解结果。从流场结果可以清晰地观察到金属液的流动形态。在充型初期，液态金属从浇口进入铸型型腔，由于惯性作用，呈现出较为集中的射流形态，流速较高。随着充型的进行，液态金属逐渐扩散，与铸型壁面接触后，流速逐渐降低，在铸型型腔的拐角和狭窄部位，由于流动阻力的增加，液态金属的流速明显减小，甚至可能出现漩涡等复杂流动现象。通过分析流速矢量图和流线图，可以直观地了解液态金属在铸型型腔内的流动路径和速度分布情况，为评估充型过程的合理性提供依据。在温度场方面，模拟结果展示了金属液在充型过程中的温度分布和传热规律。随着液态金属的流动，热量不断从高温的液态金属传递到低温的铸型中。在充型初期，液态金属的温度较高，铸型的温度较低，两者之间存在较大的温度梯度，导致热量传递速率较快。随着充型的进行，铸型被逐渐加热，温度升高，液态金属与铸型之间的温度梯度逐渐减小，热量传递速率也随之降低。在铸型的薄壁部位和与液态金属接触面积较大的区域，热量传递较快，温度变化较为明显；而在铸型的厚壁部位和远离液态金属流动路径的区域，热量传递较慢，温度变化相对较小。通过分析温度场分布云图和温度随时间变化曲线，可以深入了解充型过程中的传热特性，为预测铸件的凝固过程和质量缺陷提供重要信息。5.2凝固过程温度场与应力场耦合计算5.2.1材料参数在凝固过程中，材料参数的准确设置对于模拟结果的可靠性至关重要。除了考虑材料在充型过程中的热物理性质外，还需要特别关注材料的热膨胀系数、弹性模量和屈服强度等随温度的变化。热膨胀系数描述了材料在温度变化时的体积变化特性，它对铸件在凝固过程中的收缩和变形有着重要影响。一般来说，材料的热膨胀系数随着温度的升高而增大，在凝固过程中，由于铸件内部温度不均匀，不同部位的热膨胀程度不同，会产生热应力，进而导致铸件变形或产生裂纹。例如，对于铝合金材料，其热膨胀系数在20-25×10⁻⁶/℃之间，在凝固过程中，需要准确考虑其热膨胀特性，以预测铸件的变形情况。弹性模量反映了材料抵抗弹性变形的能力，在凝固过程中，随着温度的降低，材料的弹性模量会逐渐增大。准确设置弹性模量随温度的变化关系，能够更准确地模拟铸件在凝固过程中的应力应变状态。屈服强度则是材料开始发生塑性变形的临界应力，当铸件内部的应力超过屈服强度时，会产生塑性变形。在凝固过程中，由于热应力的作用，铸件某些部位可能会发生塑性变形，影响铸件的尺寸精度和质量。因此，需要根据材料的特性和实际凝固过程，合理设置屈服强度随温度的变化参数，以准确预测铸件的塑性变形情况。5.2.2计算结果及分析通过在ANSYSMechanical中进行凝固过程温度场与应力场的耦合计算，得到了丰富的计算结果。从温度场结果来看，在凝固初期，铸件的温度较高，随着时间的推移，热量逐渐从铸件表面向周围环境散失，温度逐渐降低。在铸件的厚壁部位，由于热量积聚较多，凝固速度相对较慢，温度下降较为平缓；而在薄壁部位，热量散失较快，凝固速度相对较快，温度下降较为迅速。在固液两相区，由于相变潜热的释放，温度下降趋势减缓，出现了温度平台现象。通过分析温度场分布云图和温度随时间变化曲线，可以清晰地了解铸件在凝固过程中的温度变化规律，为评估铸件的凝固质量提供重要依据。在应力应变状态方面，模拟结果显示，在凝固过程中，由于铸件各部位的温度变化不均匀，会产生热应力。在铸件的表面和薄壁部位，由于温度下降较快，收缩受到内部较热部位的约束，会产生拉应力；而在铸件的内部和厚壁部位，由于温度下降较慢，受到外部收缩的影响，会产生压应力。当热应力超过材料的屈服强度时，铸件会发生塑性变形，导致尺寸精度下降。通过分析应力场分布云图和应变分布云图，可以直观地了解铸件内部的应力应变状态，确定可能出现变形和裂纹的区域。例如，在铸件的拐角处和壁厚突变处，由于应力集中现象较为明显，容易产生裂纹，通过模拟结果可以提前发现这些潜在的问题，为优化铸造工艺提供指导。5.3缺陷预测5.3.1卷气和浇不足缺陷预测在铸造过程中，卷气和浇不足是常见的缺陷类型，准确预测这些缺陷对于提高铸件质量至关重要。卷气缺陷通常是由于液态金属在充型过程中卷入气体而形成的，它会降低铸件的力学性能和致密性。预测卷气缺陷的方法主要基于流速、压力和气体含量等判据。当液态金属的流速过高时，容易产生紊流，导致气体卷入。通过模拟充型过程中的流速分布，当流速超过一定阈值时，认为该区域存在卷气风险。例如，在某一铸造工艺中，当液态金属的流速超过[具体流速阈值]时，该区域出现卷气的可能性较大。同时，压力分布也会影响卷气的产生，在压力较低的区域，气体更容易析出并被液态金属卷入。通过分析模拟结果中的压力场分布，确定压力较低的区域，结合流速情况，评估卷气缺陷的可能性。此外，考虑液态金属中的初始气体含量以及气体在充型过程中的溶解和析出行为，综合判断卷气缺陷的产生位置和程度。浇不足缺陷则是由于液态金属未能完全填充铸型型腔而导致的，它会使铸件尺寸不完整，影响其使用性能。预测浇不足缺陷主要依据液态金属的充型能力和流动阻力。在模拟过程中，通过跟踪液态金属的流动前沿，观察其是否能够在规定时间内到达铸型型腔的各个部位。如果液态金属在流动过程中受到较大的阻力，如铸型型腔的狭窄通道、复杂结构等，导致流速降低过快，可能无法完全填充型腔，从而产生浇不足缺陷。通过分析模拟结果中的流速分布和流动时间，确定可能出现浇不足的区域。例如，在某一复杂铸件的模拟中，发现铸型型腔的某一薄壁部位，由于液态金属流动阻力较大，在充型后期流速明显降低