基于Abaqus/Explicit的钛基功能梯度复合材料动态冲击模拟方法：参数分析、网格策略、数值伪影抑制与CAE建模

《Materials》：A Dynamic Impact Simulation Method for Titanium-Based Functionally Graded Composites Based on Abaqus/Explicit: Parametric Analysis, Mesh Strategy, Numerical Artifact Mitigation, and CAE Modeling

　　

均匀TC4钛合金存在强度-塑性权衡问题，在高应变率冲击下表现出不足的抗侵彻能力。TiCp增强的功能梯度钛基复合材料（FGTMCs）具有“外硬内韧”的梯度结构，是很有前景的轻质装甲材料。使用Abaqus/Explicit有限元模拟对FGTMCs进行动态冲击分析，可以捕捉瞬态变形和损伤演化，同时能够快速评估不同材料的抗冲击性能；然而，该领域的研究仍然很少。研究人员对关键模拟参数进行了系统参数分析，包括TC4钛合金的Johnson-Cook（J-C）本构和损伤模型参数标定、网格划分策略优化、沙漏控制方案、模型尺寸和边界条件，以抑制“鬼网格”数值伪影。研究人员设计了TC4和三种典型钛基复合材料的材料属性分配，以及构建功能梯度材料模型的方法。比较了两种弹靶匹配配置（小弹/薄靶 vs. 大弹/厚靶），并确定了700 m/s小口径钨弹撞击50 mm TC4靶板的最优模型。参数分析表明，损伤参数D₄=0.1显著提高了数值稳定性，而沿侵彻路径进一步细化的分级网格策略平衡了计算精度和效率。使用优化后的材料系统，模拟结果再现了钛合金在冲击下的三阶段损伤演化——成坑、塑性侵彻和背面崩落——为梯度装甲材料的结构优化提供了可靠的数值支持。

**论文解读：基于Abaqus/Explicit的钛基功能梯度复合材料动态冲击模拟方法研究**

**研究背景与问题**
钛合金（尤其是Ti-6Al-4V，即TC4）因其低密度、高比强度、优异耐腐蚀性和良好弹道性能，在装甲防护领域备受关注。然而，均匀TC4存在强度-塑性权衡问题，在高应变率冲击下抗侵彻能力不足。通过引入TiC颗粒增强相，构建具有“外硬内韧”梯度结构的功能梯度钛基复合材料（FGTMCs），可有效提升抗冲击性能。但实验弹道测试成本高、耗时长，且难以捕捉弹丸侵彻过程中的瞬态变形和失效机制，因此数值模拟成为重要补充。Abaqus/Explicit有限元软件能够模拟高速冲击事件，但关键参数（如网格密度、沙漏控制、分析步设置）的选择对精度和效率影响显著，尤其是减缩积分单元在高应变率下易产生“鬼网格”零能量变形模式，严重破坏模拟可靠性。目前，针对钛合金及FGTMCs的动态冲击模拟参数化研究尚属空白。为此，研究人员开展了系统性参数分析，旨在建立可靠的有限元建模框架，优化Johnson-Cook（J-C）本构与损伤参数、网格策略及沙漏控制方案，并比较不同弹靶配置，为梯度装甲材料设计提供数值工具。该论文发表在《Materials》。

**主要关键技术方法**
研究人员采用Abaqus/Explicit（Abaqus/CAE 2016）进行动态冲击模拟。使用C3D8R三维八节点线性减缩积分单元离散靶板，并启用增强沙漏控制抑制零能量模式。通过Johnson-Cook（J-C）本构模型和损伤准则描述TC4及复合材料的力学行为，基于文献数据标定参数。采用分层分区方法构建FGTMC模型：沿厚度方向将靶板分为多层，每层赋予不同材料属性，相邻层通过共享节点实现位移连续，模拟冶金结合。参数分析包括：对比三组损伤参数集（重点分析D₄的影响）、评估网格尺寸敏感性（分级网格：侵彻路径区域1.5–2 mm，外围15 mm）、比较两种弹靶配置（小弹/薄靶：φ21 mm×56 mm钨弹、50 mm TC4靶；大弹/厚靶：φ26 mm×312 mm长杆弹、500 mm靶）。优化后选定700 m/s小口径钨弹撞击50 mm TC4靶板为标准条件。

**研究结果**

**3.1. 模型标定与参数优化**
- **3.1.1 Johnson-Cook损伤参数D₄敏感性分析**：通过对比三组参数集（Set 1: D₄=0.01；Set 2: D₄=0.1；Set 3: 调整D₂、D₃、D₅），发现在1200 m/s小弹冲击下，Set 1出现严重沙漏“鬼网格”畸变；Set 2（D₄=0.1）显著改善数值稳定性，再现了成坑、塑性侵彻、完全穿孔三阶段；Set 3进一步抑制了伪影。残速分别为1015、920、825 m/s，表明D₄增大可提高高应变率下断裂应变率敏感性，及时删除过度畸变单元，从而抑制零能量模式。
- **3.1.2 网格敏感性与收敛性**：在大弹/厚靶配置中，粗网格（15 mm）在500 m/s和600 m/s时均未预测穿孔，而细网格（5 mm局部细化）预测了穿孔，表明粗网格严重低估局部应变集中。采用沿侵彻路径局部细化的分级网格策略，可在精度与效率间取得平衡。
- **3.1.3 沙漏单元形成与抑制**：大弹/厚靶配置中，侵彻路径区域出现持续沙漏模式，即使增强沙漏控制也难以完全消除。通过将靶厚从500 mm减至200 mm可移除未受应力波影响的区域，残速偏差<4%，计算时间减少约70%，但失去了厚靶模拟的原始意义。小弹/薄靶配置通过组合措施（增强沙漏控制、D₄=0.1、适当细化网格）有效控制了沙漏伪影，因此被推荐为参数分析和机理研究的标准模型。

**3.2. 材料系统与梯度模型模拟验证**
- **3.2.1 TC4与典型钛基复合材料属性分配**：基于混合法则和文献数据，设计了TC4及三种TMCs（TiCp含量不同）的J-C本构参数，如表4所示。
- **3.2.2 材料系统参数化与梯度模型开发**：采用分层分区方法，将靶板沿厚度分为多层，每层赋予逐步变化的材料属性，相邻层共享节点，无需定义宏观界面。以三层梯度材料为例，成功模拟了“外硬内韧”梯度响应。

**3.3. 侵彻过程分析**
- **3.3.1 三阶段侵彻机制**：模拟结果清晰再现了FGTMC靶板在高速弹丸冲击下的三阶段损伤演化：Stage I（成坑）——弹丸撞击表面产生局部高压，形成冲击坑，需高硬度表面钝化弹头；Stage II（塑性侵彻）——弹丸持续压缩产生绝热剪切带，应力波向背面扩散，需过渡层兼具强度与韧性；Stage III（崩落/穿孔）——反射应力波导致背面崩落，需高韧性背层吸收能量。
- **3.3.2 速度阈值与侵彻深度**：在20–1200 m/s范围内进行系统模拟，发现20 m/s时仅弹性变形；80 m/s时表面塑性变形开始（Stage I）；200 m/s时出现侵彻特征（Stage II）；300 m/s时Stage II完全发展；700 m/s时进入Stage III但未完全穿孔，且无沙漏畸变；1200 m/s时完全穿孔。综合考虑数值稳定性与物理现象丰富性，选定700 m/s为最优模拟条件。
- **3.3.3 抗冲击FGTMC设计准则**：基于模拟结果，提出设计准则：表面需高硬度（高TiC含量）以钝化弹头；中间层硬度逐渐降低、韧性增加，以缓解界面应力集中；背层需高断裂韧性与延性，防止灾难性崩落。平滑梯度可减少内部应力波反射，而陡峭梯度易导致层间脱粘。

**总结讨论与结论**
讨论部分指出，小弹/薄靶模型在数值稳定性、计算效率和机理清晰度上优于大弹/厚靶模型，适用于参数筛选和材料优化；大弹/厚靶模型虽更贴近实际装甲场景，但数值困难较多，建议仅在最终工程验证中使用。建立的仿真框架可快速评估梯度设计参数，为定制抗冲击FGTMC结构提供实用指南。

**研究结论翻译**：本研究基于Abaqus/Explicit对FGTMCs的动态冲击模拟方法进行了系统研究，得出以下结论：(1)建立了TC4钛合金动态冲击的标准化模拟参数体系。通过系统参数分析，优化后的Johnson-Cook损伤参数D₄=0.1有效抑制了高应变率下中心变形区的沙漏“鬼网格”数值畸变，显著提高了计算稳定性。沿侵彻路径采用1.5–2 mm细化的分级网格，在捕捉应力-应变梯度与计算效率之间取得平衡。确定小口径钨弹（φ21 mm×56 mm，300 g）撞击50 mm TC4靶板为机理研究的最优弹靶配置，其数值稳定性和计算经济性明显优于大弹/厚靶方案。(2)构建了涵盖均匀和梯度材料系统的动态冲击模拟模型。设计了TC4和三种典型TMCs的材料属性参数体系，并通过多区域截面定义方法建立了功能梯度钛基复合材料（FGTMCs）的有限元模型，实现了沿厚度方向材料属性的梯度赋值，为梯度装甲结构的多参数优化提供了可扩展的模拟平台。(3)重现了三阶段侵彻损伤演化规律，阐明了梯度结构的抗冲击机理。模拟结果清晰呈现了成坑（要求表面高硬度）、侵彻（要求过渡层匹配强度与韧性）和崩落（要求背层高韧性）三阶段损伤演化过程，验证了“外硬内韧”梯度构型与弹丸多阶段侵彻载荷特征的物理匹配关系。所建立的模拟方法可有效指导梯度装甲材料关键结构参数的优化设计，为开发轻质高抗冲击梯度复合材料提供可靠的数值支持。