功能梯度SMA悬臂梁ABAQUS-UMAT有限元仿真

Kang Z. et al. | Mathematics 2024, 12, 282 | ABAQUS-UMAT二次开发

一、研究背景与工程问题

功能梯度形状记忆合金将FGM梯度材料的消除应力集中优势与SMA的超弹性和形状记忆效应结合，在航空航天、微机电系统和军事等高端领域具有广泛应用前景。然而，FG-SMA材料参数沿特定方向连续变化，传统均质SMA有限元方法无法直接适用，工程中缺乏有效的数值计算手段。Kang等人基于形状记忆因子概念构建了SMA增量本构方程，通过ABAQUS UMAT用户子程序二次开发，提出一种离散建模方法——沿梯度方向将梁分层/分段、每层材料参数取该位置坐标值、层间位移连续——从而实现了FG-SMA结构的有限元数值模拟。

图1 横向功能梯度SMA悬臂梁——沿厚度方向离散为10层，材料参数梯度变化

二、数值方法与模型验证

UMAT子程序采用增量迭代算法：整个加载过程划分为若干增量步，每步内材料参数保持不变。正相变过程通过等效应力与临界应力（σms、σmf，温度影响系数CM=13MPa/℃）的比较判定马氏体相变启停；逆相变以CA=8MPa/℃控制奥氏体相变边界。模型用C3D20R实体单元，环境温度60℃，NiTi基材料参数取EA=70GPa、EM=30GPa、εL=6.9%、泊松比0.43。横向功能梯度梁沿厚度方向分10层（梯度参数e=5），轴向功能梯度梁沿长度方向分20段。我们对均质梁首先进行了基准校核——加载3.0kN时无相变、3.5kN部分马氏体相变、4.0kN触发完全马氏体相变，自由端挠度分别为16.96/22.11/52.00mm，完全相变区域不在固定端而在距固定端约10mm处、最大深度3mm。

图2 SMA均质悬臂梁截面Mises应力分布——沿厚度方向对称，最大值在上/下表面固定端

图3 SMA均质悬臂梁挠度曲线——P=4.0kN时已激发完全马氏体相变，自由端最大挠度52mm

三、核心发现与工程指导

横向FG-SMA梁：应力沿厚度方向**非对称**分布，最大值在上表面固定端（梯度参数导致上表面材料刚度更大）；马氏体相变也从下表面固定端先触发（该层σms最小）；同等挠度下承载力为均质SMA梁的3倍——这是工程设计的核心优势。轴向FG-SMA梁：应力沿厚度方向对称但集中在固定端较均质梁更显著，马氏体相变起始位置距固定端仅7.5mm（均质梁为10mm），承载力介于横向FG梁与均质梁之间。三者排序：**横向FG梁 > 轴向FG梁 > 均质梁**——同等挠度下横向FG承载力约为轴向FG的2倍，这一规律为FG-SMA结构按需设计提供了关键数据支撑。这里我们认为，离散建模方法避免了繁琐的理论非线性求解，通过工程中常用的ABAQUS+UMAT路径即可模拟FG-SMA的相变与力学耦合行为，对实际工程应用具有直接参考价值。

图4 横向FG-SMA截面形状记忆因子分布——马氏体相变沿厚度方向非对称，下表面固定端先触发

图5 横向FG-SMA与均质梁载荷-位移对比——同等挠度下横向FG梁承载力为均质梁3倍

图6 轴向/横向FG-SMA与均质梁三者对比——横向FG梁承载力最强，轴向次之

四、我们提供的服务

① ABAQUS UMAT子程序定制开发——基于形状记忆因子或其他本构模型（Auricchio/Lagoudas/Brinson等），实现SMA超弹性、形状记忆效应的数值仿真 ② 功能梯度材料（FGM/FG-SMA）离散建模方案设计与有限元实现 ③ NiTi基/铜基/铁基SMA结构的相变-力学-热多场耦合分析 ④ 航空航天/生物医疗/微机电等领域SMA结构的工程仿真与技术咨询 ⑤ ABAQUS-UMAT/VUMAT/UEL二次开发技术培训与项目协作 欢迎通过公众号"320科技工作室"联系我们，获取更详细的技术方案和工程案例。

如有需要，欢迎通过公众号"320科技工作室"联系我们。