评估固体火箭羽流中氧化铝颗粒冲击下延性合金侵蚀率的计算方法

《Computational Particle Mechanics》：Computational approach for evaluating the erosion rate of ductile alloys impacted by alumina particles in solid rocket plumes

　　

本研究介绍了一个全面的计算框架，用于预测典型固体火箭发动机羽流中高速氧化铝颗粒作用下延性合金的稳态侵蚀率。该框架结合了有限元法（FEM）模拟与Johnson-Cook本构和损伤模型，并采用周期性边界条件，能够在重复颗粒冲击下准确预测稳态侵蚀率。这一能力超越了先前仅限于单次冲击或刚性颗粒假设的模型。周期性边界条件的新颖应用允许系统地探索稳态侵蚀率，同时最小化计算成本。研究人员直接比较了固态和熔融氧化铝颗粒，揭示了颗粒相强烈影响侵蚀行为。虽然两种相在法向冲击下产生相当的侵蚀，但只有固态颗粒在倾斜角度下增强侵蚀。此外，证明了颗粒形状效应，立方体颗粒产生的侵蚀率比球形颗粒高出27%。该框架还突出了材料特定的抗侵蚀性，表明4340钢在高速氧化铝颗粒冲击下优于Ti-6Al-4V。预测的侵蚀趋势与现有实验数据高度一致，支持了计算方法的有效性。本研究为固体火箭发动机羽流中的侵蚀现象提供了有价值的见解，并建立了一个预测延性合金侵蚀率的框架，在航空航天工程、涡轮机械和石油天然气行业具有潜在应用。

论文解读文章

**研究背景与问题**

现代火箭操作强调可重复使用性和寿命，固体火箭发动机产生的高温、高速、含颗粒羽流对发射台材料构成严重侵蚀威胁。侵蚀主要通过热通量和颗粒冲击两种机制实现，其中颗粒冲击在许多情况下占据主导地位。尽管已有大量关于单次颗粒冲击的实验和计算研究，但这些研究在预测多次冲击导致的大规模机械侵蚀时存在局限，因为先前冲击造成的损伤会显著影响后续侵蚀。现有理论模型多局限于刚性球形颗粒法向冲击的理想配置，难以处理角度、形状和状态等复杂因素。实验研究虽提供了有价值见解，但速度范围有限（如80-300 m/s），远低于火箭羽流中的典型值（可达1050 m/s）。因此，需要开发一种能够模拟多次颗粒冲击并预测稳态侵蚀率的计算方法，以量化实际火箭羽流条件下的侵蚀行为。

**研究内容与意义**

本研究旨在通过三维有限元法（FEM）模拟，量化与火箭羽流相关的侵蚀率。研究人员系统研究了氧化铝颗粒速度（400-900 m/s）、入射角（0-90°）、形状（球形和立方体）和状态（固态和液态）对4340钢和Ti-6Al-4V靶板侵蚀的影响。目标是确定速度和角度的影响、比较颗粒形状和状态的作用，以及对比不同靶材的响应。该研究发表在《Computational Particle Mechanics》。研究意义在于提供了一个经过验证的计算框架，用于预测延性合金在固体火箭发动机羽流中的稳态侵蚀率，可应用于航空航天工程、涡轮机械和石油天然气行业，为保护发射系统、延长寿命提供设计依据。

**关键技术方法**

研究人员采用有限元法（FEM）模拟，结合Johnson-Cook本构和损伤模型，并应用周期性边界条件以模拟无限大靶板。使用IMPETUS Afea^?软件在图形处理器（GPU）上运行，以加速计算。固态氧化铝颗粒采用线弹性模型，液态氧化铝颗粒采用包含密度、线性体积模量和粘度的流体模型。靶材（4340钢和Ti-6Al-4V）的弹塑性响应通过Johnson-Cook模型描述，材料参数来自文献（如Johnson-Cook原始测试及后续研究）。颗粒尺寸为50 μm（球形）或等体积的立方体（边长40.3 μm），颗粒冲击位置在靶面随机分布。通过单元尺寸研究、网格分辨率研究和孵化质量（incubation mass）研究，确保模拟达到稳态侵蚀。验证案例将FEM模拟结果与Hutchings等人（文献[46]）的单次钢球冲击实验对比，误差在实验误差范围内（±26%），验证了模型准确性。

**研究结果**

**4.1 本构模型选择**
通过比较Zerilli-Armstrong模型、Arrhenius型模型和IMPETUS内置“MAT_METAL”模型，研究人员最终选择Johnson-Cook本构和损伤模型，因其在文献中广泛用于高速冲击问题，且能有效捕捉应变硬化、挤出、低周疲劳和熔化等侵蚀机制，但不包括粘着材料转移。

**4.2 单元尺寸确定**
通过模拟不同表面面积（2a×2b至5a×5b）的单元胞，发现4a×4b×3d（a为损伤椭圆长轴，b为半短轴，d为颗粒直径）的尺寸使侵蚀率差异小于1%，因此确定为标准单元尺寸。对于不同入射角，单元尺寸根据损伤椭圆调整（见表3）。

**4.3 网格分辨率**
对靶板网格（元素边长2.0-3.33 μm）和颗粒网格密度（球形3-7，立方体3-9）进行独立研究。结果表明，5 μm立方体元素、球形颗粒网格密度6（固态）、5（液态）以及立方体颗粒7×7×7网格可确保侵蚀率变化小于2-4%，满足精度要求。

**4.4 状态方程（EOS）的影响**
对比包含Mie-Grüneisen EOS与忽略EOS的模拟，在900 m/s法向冲击下，侵蚀率差异小于3%。由于冲击速度（≤900 m/s）远低于目标材料弹性波速（>3 km/s），EOS影响可忽略，故在后续模拟中未包含EOS。

**4.5 侵蚀率预测**
所有模拟结果以相对侵蚀率（归一化于900 m/s固态球形颗粒对钢的侵蚀率889.1 mg/g）表示，展示于图7。

**4.5.1 入射角的影响**
对于固态球形颗粒，4340钢和Ti-6Al-4V靶的侵蚀率在30-70°入射角出现峰值，约50°时最大。立方体颗粒在50°时侵蚀率比法向高24%（1.24倍）。这一结果与Oka等人实验一致，侵蚀率峰值位置取决于靶材硬度，且由重复塑性变形和切割作用共同导致。

**4.5.2 颗粒形状的影响**
立方体颗粒比球形颗粒侵蚀率高27%（1194 mg/g vs 943.8 mg/g，在50°入射角）。单次冲击模拟显示，立方体颗粒的不同取向（平面、棱边、角）侵蚀率分别比球形高73%、12%和31%，且坑深和坑宽增加。这表明棱角增强穿透，侵蚀率取决于具体几何和取向。

**4.5.3 法向冲击下颗粒状态的影响**
在法向冲击下，液态和固态球形颗粒的侵蚀率几乎相同（差异<1%）。液态颗粒虽单次造成更少材料去除，但损伤面积更大，多次冲击后稳态侵蚀率与固态相当。这一结果仅通过多次冲击模拟获得。

**4.5.4 液态颗粒倾斜冲击**
液态颗粒侵蚀率随入射角减小而降低，遵循Baker等人提出的关系：E_θ = K(V sinθ - V_c)ⁿ¹。本研究对液态氧化铝颗粒在900 m/s下拟合得到V_c=1 m/s，n1=3.45。固态与液态颗粒在掠射角下的侵蚀行为差异显著，表明颗粒状态在倾斜冲击中至关重要。

**4.5.5 速度效应**
固态球形颗粒法向冲击时，900 m/s侵蚀率是400 m/s的3.70倍，符合幂律关系E₉₀ ∝ Vⁿ²，指数n2=1.6。这与实验报道的2-3指数范围基本一致，表明速度对侵蚀率有显著影响。

**4.5.6 材料抗侵蚀性能**
4340钢的抗侵蚀性远优于Ti-6Al-4V。在900 m/s法向冲击下，Ti-6Al-4V的侵蚀率是4340钢的2.14倍（1907 mg/g vs 889.1 mg/g）。该差异无法通过断裂韧性和硬度等常见材料参数直接预测，需通过计算或实验确定。

**4.5.7 验证案例**
模拟Hutchings等人的单次9.5 mm钢球以30°入射角冲击钢靶的实验（速度310 m/s）。FEM模拟预测侵蚀率为7.56 mg/g，实验经验公式预测9.77 mg/g，偏差23%，在实验误差±26%范围内。模拟坑长10.7 mm（实验11.8 mm），坑宽7.0 mm（实验9.0 mm），偏差分别为9%和22%，验证了模型的准确性。

**总结讨论与结论翻译**

**讨论总结**
本研究通过FEM框架成功预测了多次颗粒冲击下的稳态侵蚀率，揭示了入射角、颗粒形状、颗粒状态和靶材类型对侵蚀率的影响。关键发现包括：法向冲击下固态和液态颗粒侵蚀率相同；固态颗粒在50°角侵蚀增强；立方体颗粒侵蚀率高于球形；液态颗粒侵蚀率随角度正弦下降；侵蚀率随速度呈幂律（指数1.6）；4340钢抗侵蚀性优于Ti-6Al-4V；EOS影响可忽略。验证案例表明模拟结果与实验一致。该框架为火箭羽流侵蚀研究提供了有效工具，但也指出未来需耦合热-机械侵蚀、考虑颗粒碎片化及实际羽流环境中的复杂相互作用。

**研究结论翻译**
在计算固体火箭发动机羽流侵蚀率的具体背景下，本研究使用有限元框架与Johnson-Cook本构和损伤模型及周期性边界条件，量化了感兴趣颗粒速度、入射角、形状和相态的侵蚀率。显著发现包括：在900 m/s法向冲击下，液态和固态球形颗粒导致相同侵蚀率，表明对于直接羽流冲击，氧化铝颗粒相态在侵蚀计算中相对不重要；与法向冲击相比，固态颗粒在50°角冲击可增强侵蚀高达24%，而液态颗粒侵蚀率随入射角正弦函数减小，因此对于掠射冲击，氧化铝状态是侵蚀计算的关键因素；在900 m/s冲击下，随机取向的立方体颗粒比球形颗粒侵蚀率增加高达27%，由于火箭羽流中氧化铝本质上是棱角形的，理解相对于理想球形颗粒的增强至关重要；单个立方体颗粒的特定冲击取向对侵蚀有显著影响，增强范围从仅12%到高达73%相对于单个球形颗粒；侵蚀率随速度呈幂律缩放，固态球形颗粒在本研究中指数为1.6；4340钢具有优于Ti-6Al-4V的抗侵蚀性，钛合金侵蚀量约高2.1倍，在设计中需平衡热响应和密度等其他因素；对于≤900 m/s的颗粒冲击，在FEM模型中包含EOS导致模拟侵蚀率变化小于3%，表明其在开发材料模型时重要性较小。这些结果得益于一种新颖的FEM框架，该框架集成了Johnson-Cook材料建模与周期性边界条件，使得在重复颗粒冲击下预测稳态侵蚀率成为可能。这种创新方法超越了先前单次冲击模型，提供了一个经过验证的计算途径，用于量化与火箭发动机羽流及相关高速颗粒环境相关的侵蚀机制。