⚡ electrical engineering

A Finite Element Process Modeling Framework for Predicting Residual Stresses and Microcracking in Ceramic Matrix Composites

本文提出了一种基于物理的有限元框架，该框架将固化和热解阶段联系起来，以预测陶瓷基复合材料中的残余应力、收缩和微裂纹，证明了其在模拟聚合物浸渍和热解工艺过程中缺陷形成方面的能力。

原作者： Michael N. Olaya, Marianna Maiaru

发布于 2026-06-25

📖 1 分钟阅读☕ 轻松阅读

原作者： Michael N. Olaya, Marianna Maiaru

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，你正试图烘焙一种非常特别、超强韧的陶瓷蛋糕。但不同于使用面粉和鸡蛋，你使用的是一种在加热后会变成岩石状的液态塑料。这就是科学家制造**陶瓷基复合材料（CMCs）**的方法，这种材料被用于像高超音速喷气机这样需要承受极端高温的设备中。

问题在于，这种“烘焙”过程非常棘手。随着塑料变成岩石，它会发生收缩并发生化学变化。这会在材料内部产生看不见的内应力（就像被拉得太紧的橡皮筋），从而导致内部形成微小的裂纹和孔洞（空隙）。如果这些裂纹形成，最终的产品就会变得脆弱且毫无用处。

这篇论文介绍了一个计算机模拟框架——一个数字化的“虚拟烤箱”——它可以帮助工程师在实际制造零件之前，精准预测裂纹会在何时何地发生。

以下是作者的“虚拟烤箱”是如何运作的，通过简单的步骤进行拆解：

1. 两阶段“烘焙”过程

作者意识到，你不能直接跳到高温阶段。你必须分两个不同的阶段来模拟这个过程：

第一阶段：“凝固”阶段（固化）： 想象一下液态塑料就像明胶。首先，你必须让它凝固成固体果冻。计算机模拟了这个“凝固”过程（称为 C-staging）。在此期间，材料开始变硬，但还没有明显的收缩。计算机也会检查是否仅因材料硬化就产生了微小的裂纹。
第二阶段：“岩石化”阶段（热解）： 一旦果冻凝固，你就调高温度将其变为岩石。这正是奇迹（以及危险）发生的地方。材料会失去重量（气体逸出）、显著收缩，并硬化成陶瓷。计算机会追踪材料如何收缩，以及在如此剧烈的加热过程中内部压力是如何积聚的。

2. “数字孪生”微观结构

为了使模拟更加真实，作者并没有仅仅观察一块材料。他们构建了材料微观结构的数字孪生。

纤维： 把它们想象成蛋糕内部的“骨架”或“钢筋”。他们使用了碳纤维（AS4）的模型，这种纤维保持强韧且几乎不发生变化。
基体： 这是“填充物”或“蛋糕糊”（SMP-10），它包裹着纤维。这是会发生收缩和开裂的部分。
模拟过程： 他们创建了一个包含 20 根纤维并被这种“蛋糕糊”包围的微型重复数字方块，并在计算机中观察加热时的变化。

3. “裂纹检测器”

他们的框架中最巧妙的部分是一个渐进损伤模型。

想象材料有一个“断裂点”。当计算机模拟加热过程时，它会不断检查材料上的应力。
如果应力过高（就像拉扯橡皮筋直到它断裂一样），计算机就会让模拟中的那个微小点“断裂”。
一旦某个点断裂，计算机就会更新该处的材料属性，使其在该处变得更弱，从而改变应力向邻近区域传播的方式。这使得模拟能够展示出微裂纹网络是如何形成的，就像现实生活中一样。

4. 他们的发现（结果）

由于他们并没有关于这种特定材料的所有完美现实世界数据，因此他们利用了来自其他研究的最佳可用信息，并做了逻辑性的推测（例如使用一种类似的材料作为“替代品”）。

他们的数字烤箱揭示了以下内容：

大幅收缩： 当材料从塑料转变为陶瓷时，它收缩了约 8%。对于一个固体物体来说，这可是一个巨大的变化！
应力峰值： 最大的危险并不发生在烘焙的最末端。最高的应力（以及裂纹开始产生的地方）发生在早期，即温度刚刚开始上升的时候（大约在 300–400°C 左右）。
裂纹： 由于“蛋糕糊”在收缩的同时，“骨架”（纤维）保持不变的大小，导致“蛋糕糊”被拉开了。在模拟结束时，约 85% 的基体材料出现了某种程度的裂纹。
结果： 最终的材料比完美的、无裂纹的陶瓷件要弱得多，也更不坚硬。

为什么这很重要

作者并不是说他们的计算机模型现在就是完美的；他们承认由于获取现实世界的数据非常困难，他们不得不通过假设来填补一些空白。然而，他们建立了一个基础。

把这个框架看作是制造领域的飞行模拟器。正如飞行员使用模拟器在不坠毁真实飞机的情况下学习如何在风暴中驾驶飞机一样，工程师现在可以使用这个工具来测试不同的加热配方。他们可以询问：“如果我们加热得慢一点会怎样？”或者“如果我们改变纤维的排列方式会怎样？”以便在计算机内部测试是否能阻止裂纹形成，从而在进行昂贵的现实实验之前节省成本。

简而言之： 他们构建了一个数字工具，用于模拟超强韧陶瓷的“烘焙”过程，证明了该过程会导致材料显著收缩和开裂，并证明了我们需要通过精确控制加热过程来解决这些缺陷。

技术摘要：一种用于预测陶瓷基复合材料残余应力和微裂纹的有限元过程建模框架

问题陈述
陶瓷基复合材料（CMCs）因其在超过 1200–1600 °C 高温下卓越的热阻性和结构完整性，成为下一代高超音速应用的关键材料。然而，主要的制造方法——聚合物浸渍热解法（PIP）存在显著的缺陷形成问题，特别是微裂纹和孔隙。这些缺陷使得实现完全致密化需要多次、高成本的再浸渍循环。其根本原因在于复杂的耦合热化学反应，以及在热解过程中先驱体聚合物（PCP）内力学性质的同时演变。这些过程产生了难以通过实验量化的巨大残余应力，因为在极端加工温度下，反应动力学、质量损失、致密化、体积收缩以及气体传输过程是同时发生的。因此，迫切需要物理驱动的计算框架来模拟这些耦合机制并预测缺陷形成。

研究方法
本研究引入了一种基于物理的有限元（FE）框架，旨在进行 CMCs 的数字化制造，特别侧重于微观尺度模拟。该方法采用了一种创新的两步工艺，将固化阶段（C-staging）直接与随后的热解过程相联系。

控制物理机制： 该框架求解三个耦合方程：瞬时线性弹性本构方程、带有放热源项的瞬时热传导方程，以及用于固化和热解的动力学模型。总应变被分解为机械应变、热应变和收缩应变。收缩量在固化阶段被建模为固化度（degree of cure）的线性函数，而在热解阶段则作为体积产率（由质量损失和致密化导出）的函数。
损伤建模： 引入了基于 Bažant 和 Oh 的裂纹带理论的渐进损伤模型，以捕捉原位微裂纹。损伤起始遵循最大主应力准则，后峰行为通过具有应变软化的线性牵引-分离规律进行建模。标量损伤参数（ $D$ ）用于退化柔度张量。
实现方式： 控制方程通过用户自定义子程序（UMAT 和 UMATHT）在 Abaqus 有限元环境中实现。该框架采用两阶段分析：固化分析建立初始状态（残余应力、损伤和材料属性），随后将其作为热解分析的初始条件导入。
材料建模： 由于针对特定体系（SMP-10 基体与 AS4 碳纤维）的实验数据有限，本研究构建了一个替代材料描述。固化特性源自类比系统（酚醛 SC-1008 和环氧 RIM R135-H1366），而热解动力学和收缩特性基于 Key 等人建立的 SMP-10 模型。热解过程中的材料属性被假设为随陶瓷密度变化的函数而演变。
几何结构： 该框架在包含 20 根随机分布的 AS4 碳纤维（体积分数为 45%）于 SMP-10 基体中的代表性体积单元（RVE）上进行演示，并采用周期性边界条件。

核心贡献

统一的两步式框架： 本文提出了一个直接连接固化与热解的框架，允许将固化过程中积累的残余应力和损伤状态作为热解模拟的初始条件。这与以往通常将这两个阶段分开处理或侧重于未显式模拟纤维的介观尺度研究不同。
集成渐进损伤模型： 引入渐进损伤模型使得能够模拟整个加工阶段中的原位微裂纹，从而捕捉应力演变与材料退化之间的相互作用。
物理一致性的假设： 在缺乏完整材料数据的情况下，作者开发了一种稳健的 SMP-10 属性代理模型，通过从密度变化和经过验证的动力学模型中导出属性演变，确保了物理一致性。
微观尺度分辨率： 该模型在微观尺度上显式解析了纤维-基体界面，从而能够观察到在较大尺度模型中通常会被均质化的应力集中点和损伤局部化现象。

结果
对 CMC 微观结构的模拟揭示了以下结果：

固化阶段： 固化过程中产生的残余应力相对较低（升温阶段峰值约为 7 MPa 拉应力，冷却至室温后约为 65 MPa）。损伤起始极小，阶段结束时体积平均损伤参数（ $D$ ）保持在 0.015 以下。在纤维位置较近处观察到了应力集中。
热解阶段： 显著的基体微裂纹是由热解诱导的收缩驱动的。损伤在热升温阶段早期（约 285°C，与快速分解的开始时间吻合）即开始发生。体积平均损伤参数急剧上升，在温度达到最终 1200°C 目标值的一半之前，该参数已超过 0.7。
应力与收缩： 峰值拉应力出现在热升温阶段早期（约 300–400°C），达到约 55 MPa，随后由于渐进损伤限制了进一步的应力积累。模拟预测最终体积收缩率约为 8%。
属性演变： 由于累积损伤，基体的最终弹性模量和抗拉强度较理论上完全致密 SiC 的本体属性显著退化（分别为 ~160 GPa 和 1.4 GPa）。所得微观结构表现出的破碎形态与文献观察一致。

意义与主张
作者将这项工作定位为一种用于链接固化与热解以预测 CMC 缺陷形成的计算框架的原理验证（proof-of-concept）。该研究声称建立了一个基础，用于实现“预测性、过程驱动的建模”，从而可以系统地研究加工条件（热循环、PIP 参数）如何影响残余应力、收缩和微裂纹。

论文也谦虚地指出，由于材料数据有限，目前的结果依赖于代理描述和假设。然而，作者断言该框架提供了一个理解缺陷形成耦合机制的必要平台。作者总结道，随着更准确的材料表征数据被获取以及更多物理机制（如热解气体传输）被纳入，该方法具有指导设计制造策略以提高 CMC 加工后性能的潜力。这项工作并不声称能立即解决所有的优化问题，而是提供了开展此类研究必不可少的计算工具集。