2. SMA本构模型：从一维到三维的工程视角
各位同学，今天我们来聊聊SMA的本构模型。说实话，这部分内容我当年刚接触时也觉得头大——一堆公式、一堆变量，到底怎么用？但做了十几年机器人关节设计，我慢慢发现，这些模型其实就是我们和材料“对话”的语言。
先给大家一个整体框架。SMA的本构模型，说白了就是描述“温度、应力、应变、马氏体体积分数”这四个量之间关系的数学工具。我习惯把它们分成两类：一维模型和三维模型。一维模型适合细丝、弹簧这类简单结构；三维模型则用于复杂形状的构件。
2.1 一维本构模型：Brinson模型与Liang-Rogers模型
先讲一维模型。为什么先讲它？因为简单、直观，而且我们做机器人关节时，SMA丝、SMA弹簧这些一维构件用得最多。我早期做的一个手指关节，用的就是0.5mm的SMA丝，一维模型完全够用。
2.1.1 Liang-Rogers模型
Liang-Rogers模型是1987年提出的，算是SMA本构模型的“老祖宗”。它的核心思想是：把马氏体体积分数ξ作为内变量，用余弦函数来描述相变过程。
公式长这样：
σ - σ₀ = E(ξ)(ε - ε₀) + Ω(ξ)(ξ - ξ₀) + Θ(T - T₀)其中：
- E(ξ)：弹性模量，随马氏体含量变化
- Ω(ξ)：相变模量，描述相变引起的应变
- Θ：热弹性系数
- ξ：马氏体体积分数（0~1）
嗯，这里要注意：Liang-Rogers模型假设相变过程中，马氏体体积分数ξ和温度T、应力σ之间满足余弦关系。说白了，就是用一个平滑的曲线来描述“从奥氏体变成马氏体”这个过程。
2.1.2 Brinson模型
Brinson模型是Liang-Rogers模型的“升级版”。它把马氏体分成了两种：温度诱发马氏体（ξ_T）和应力诱发马氏体（ξ_S）。为什么要分？因为这两种马氏体的行为不一样。
你想想看：温度降低时形成的马氏体，和受力变形时形成的马氏体，恢复特性完全不同。Brinson模型抓住了这个关键点。
核心公式：
σ = E(ξ)(ε - ε_L · ξ_S) + Θ(T - T₀)注意这里：只有应力诱发马氏体ξ_S贡献了相变应变ε_L。温度诱发马氏体ξ_T不贡献应变，但它会影响弹性模量。
| 对比项 | Liang-Rogers模型 | Brinson模型 |
|---|---|---|
| 马氏体分类 | 单一变量ξ | ξ = ξ_T + ξ_S |
| 相变应变描述 | 通过Ω(ξ)间接描述 | 直接使用ε_L · ξ_S |
| 适用应力范围 | 低应力（<100MPa） | 中高应力（<500MPa） |
| 计算复杂度 | 较低 | 中等 |
2.2 三维本构模型简介
好了，一维模型讲完了。但现实中的SMA构件不总是细丝或弹簧。比如我做的一个SMA驱动的柔性手腕，用的是SMA薄片，受力是三维的——有拉伸、有弯曲、还有扭转。这时候一维模型就不够用了。
三维本构模型的核心思想是：把一维的标量关系，推广到张量形式。说白了，就是把σ、ε这些量变成3×3的矩阵（张量），然后描述它们之间的关系。
目前主流的三维模型有：
- Lagoudas模型：基于热力学势函数，考虑相变和重取向
- Auricchio模型：基于塑性力学框架，适合有限元实现
- Panico-Brinson模型：Brinson模型的三维推广
三维模型的通用形式：
σ_ij = C_ijkl(ξ) · (ε_kl - ε^tr_kl - ε^th_kl)其中：
- σ_ij：应力张量
- C_ijkl：弹性刚度张量（随ξ变化）
- ε^tr_kl：相变应变张量
- ε^th_kl：热应变张量
2.3 模型选择指南
讲了这么多，到底怎么选？我给大家一个简单的判断标准：
- 构件形状：细丝/弹簧 → 一维模型；薄片/块体 → 三维模型
- 应力水平：低应力 → Liang-Rogers；中高应力 → Brinson
- 计算资源：手算/简单编程 → 一维；有限元仿真 → 三维
- 精度要求：初步估算 → 一维；详细设计 → 三维
最后说一句：模型终究是工具，不是真理。我见过太多人死磕公式，却忽略了实验验证。记住，再好的模型，也要用实验数据来校准。这是我做了十几年SMA设计最深的体会。