3D Finite Element-Based Multiphysics Simulation of a Shape Memory Alloy Hybrid Composite Module

这是一份关于《基于 3D 有限元的形状记忆合金混合复合材料模块多物理场仿真》论文的详细技术总结。

1. 研究背景与问题 (Problem)

形状记忆合金混合复合材料（SMAHC）通过将形状记忆合金（SMA）嵌入聚合物基体中，能够实现可逆且连续的形变，广泛应用于智能蒙皮、自适应后缘、机器人抓手等领域。尽管已有多种解析模型和有限元（FE）模型用于预测 SMAHC 的瞬态响应，但现有方法仍存在以下局限性：

• 耦合不足：许多模型未能充分整合电 - 热 - 力（Electro-thermomechanical）的多物理场耦合。
• 验证缺失：缺乏全面的实验验证，特别是针对复杂的三维几何结构和边界条件（如环境温度、外部机械载荷）。
• 初始化困难：难以在仿真中准确模拟 SMA 丝在室温下的初始微观结构状态（如部分去孪晶的马氏体状态），这直接影响其收缩潜力的预测。
• 缺乏通用工具：目前缺乏一种能够同时考虑完整电 - 热 - 力耦合、边界条件及三维几何细节的经过充分验证的仿真方法。

2. 方法论 (Methodology)

本研究提出了一种基于 ANSYS LS-DYNA 的耦合多物理场三维有限元（3D FE）仿真方法，具体包括：

• 多物理场求解器：集成了机械、热和电磁（EM）求解器。
  • 电磁场：采用简化的 EM 求解器（仅计算电阻，忽略感应效应），通过欧姆热作为热源与热求解器耦合。
  • 热场：模拟焦耳加热、对流冷却及环境温度影响。
  • 力学场：模拟 SMA 丝的相变收缩及复合材料的弯曲变形。
• 本构模型：采用了由 Kelly 等人提出的微观力学本构模型（参考文献 [25]），该模型基于亥姆霍兹自由能最小化原理，能够描述马氏体与奥氏体相变的成核、饱和及生长动力学。
  • 模型考虑了马氏体体积分数（$\lambda$）和名义马氏体应变张量（$\varepsilon_m$）作为内部状态变量。
• 初始状态预设（关键创新）：
  • 由于商业软件的材料卡无法直接指定初始相态（通常默认为完全奥氏体），研究引入了一个预拉伸步骤。
  • 在正式仿真前，先对缩放的 SMA 丝进行机械加载和拉伸，使其在 23°C 下从初始奥氏体转变为孪晶马氏体，再部分转变为去孪晶马氏体。这为后续仿真提供了符合物理实际的初始状态。
• 几何模型：
  • 基于商用 SMAHC 执行器（CompActive GmbH, Type A3950）构建。
  • 包含 SMA 丝、弹性基底、软中间层和锚固结构。
  • 使用 32,360 个六面体实体单元进行网格划分，利用对称性简化模型。
• 验证方案：
  • 将仿真结果与一维交错方案模型（SSM）（基于 Python 的完全耦合模型）进行对比。
  • 与实验数据（不同电流、环境温度、机械载荷下的偏转数据）进行对比。

3. 主要贡献 (Key Contributions)

1. 首次实现全耦合 3D 仿真：在 ANSYS LS-DYNA 中成功实现了 SMAHC 执行器的机械、热、电磁全耦合三维仿真，能够模拟复杂的电 - 热 - 力响应。
2. 物理驱动的初始化流程：提出了一种通过预拉伸步骤来模拟 SMA 丝初始去孪晶状态的方法，解决了商业软件无法直接定义初始微观相态的难题。
3. 微观与宏观的结合：将微观尺度的相变动力学（马氏体/奥氏体转变）直接应用于宏观三维有限元模型中，能够捕捉伪塑性、超弹性和单程/双程形状记忆效应。
4. 全面的验证体系：不仅对比了实验数据，还引入了高精度的 SSM 数值模型作为中间参考，评估了 3D 方法的预测能力。

4. 研究结果 (Results)

• 温度分布：仿真能够重现 SMA 丝沿长度方向的不均匀加热（中心快于两端），与实验观测一致。
• 加热与偏转：
  • 定性一致性：仿真结果在定性上与实验数据吻合良好，能够复现执行器偏转随温度变化的特征滞后环（Hysteresis）。
  • 定量偏差：预测的偏转量级正确，但在某些情况下（如 2A 电流）略低于实验值，且处于 95% 实验置信区间边缘。
  • 冷却差异：仿真中的冷却速度明显快于实验（实验中存在较慢的冷却过程），导致仿真中执行器复位时间早于实验。这在高环境温度下尤为明显。
• 载荷影响：
  • 在低载荷（0.1 N - 19.6 N）下，仿真与实验偏差较小。
  • 在高载荷下，仿真模型表现出比实际执行器更“硬”的特性（偏转量偏小），且伪塑性变形不如实验显著。
  • 对于较薄的弹性层（Type A/B），仿真未能完全复现实验中因层间距减小而导致的更大偏转，推测原因是软中间层的线性弹性建模及其刚度被高估。
• 与 SSM 对比：总体而言，一维 SSM 模型在定量上与实验数据的吻合度略优于 3D FE 模型，但 3D 模型提供了 SSM 无法实现的几何变形可视化和复杂边界条件处理能力。

5. 意义与展望 (Significance & Outlook)

• 工程价值：该研究证明了基于 3D 有限元的多物理场方法在设计和预测复杂 SMAHC 系统方面的可行性。尽管存在定量偏差，但其能够捕捉关键物理现象（如滞后、相变动力学），为计算机辅助工程（CAE）设计提供了有力工具。
• 局限性：目前的模型在软中间层的非线性行为、冷却过程的热惯性以及高载荷下的伪塑性变形方面仍需改进。
• 未来工作：
  • 利用模型的空间分辨率计算马氏体体积分数分布。
  • 改进弹性中间层的材料模型（从线性弹性改为非线性）。
  • 进行材料属性和几何参数的敏感性分析，以优化模型精度。

总结：该论文成功开发并验证了一种用于 SMAHC 执行器的 3D 多物理场有限元仿真框架。虽然在高载荷和冷却阶段的定量预测上仍有提升空间，但该模型在捕捉复杂物理机制和几何效应方面展现了巨大潜力，为下一代智能复合材料执行器的设计奠定了坚实基础。