从工程应力应变曲线到仿真材料卡片：一位CAE工程师的实战笔记

国高材高分子材料产业创新中心

2026年5月15日 16:45

材料卡片是仿真分析的"基因"，决定了有限元计算结果的精度上限。

在碰撞仿真、NVH分析、产品可靠性评估等场景中，材料参数设置的准确性直接影响仿真的可信度。然而，实验室提供的原始材料曲线与仿真软件所需的有效应力应变曲线之间，存在一道需要跨越的转化鸿沟。本文基于实战经验，系统梳理从材料曲线获取到仿真材料卡片生成的完整流程，供从事CAE工作的工程师参考。

一、工程应力应变曲线

1.1 材料的关键参数

开展有限元分析前，必须明确材料的几项基础参数，这些参数构成了材料卡片的骨架。

弹性模量（杨氏模量）是工程应力应变曲线屈服段的斜率，即应力与应变的比值。金属材料通常为210000 MPa或20600 MPa，塑料材料约为2350 MPa。这一参数直接决定了结构在弹性阶段的刚度表现。

图1 带引伸计拉伸测试

泊松比是材料在单向受拉或受压时，横向正应变与轴向正应变的比值，用于反映材料的横向变形特性。金属材料泊松比通常取0.34，塑料材料约为0.39。密度是质量与体积的比值，在碰撞仿真和NVH分析中尤为重要——不同单位制模型中，密度参数容易出现数量级错误，导致分析结果严重失真。

屈服强度是材料从弹性变形进入塑性变形的临界点。拉伸过程中，材料在屈服点之前仅产生弹性变形；过了屈服点则进入塑性阶段，产生永久不可恢复的变形。塑料材料由于韧性较差，拉伸试验中基本没有明显的屈服阶段，工程设计中常以产生0.2%残余应变时的应力作为条件屈服极限。

抗拉强度是材料应力值的极限点，超过此值材料即被判定破坏失效。断裂延伸率则是抗拉强度所对应的应变值，塑性应变值超过断裂延伸率时，材料同样被视为失效。

图2 应力应变曲线

1.2 获取途径

工程应力应变曲线的获取主要有三种途径，各有优劣。

第一种方式是向材料供应商直接索取，这是最理想的信息来源，尤其对于成熟牌号的商业材料，供应商通常能提供完整的测试报告。

第二种方式是委托第三方实验室进行拉伸试验，这种方法获得的数据最为准确可靠，但成本较高，适用于对仿真精度要求极高的关键零部件。

第三种方式是通过专业工具从已发表的技术文献或网络资源中"白嫖"曲线数据，再利用数字化工具提取坐标点，这种方式成本最低但数据质量参差不齐，仅推荐用于项目前期的快速可行性分析阶段。

二、工程曲线→真实曲线→有效曲线

2.1 为什么必须转换？

实验室万能试验机直接输出的拉伸曲线称为工程应力应变曲线，其定义方式为：工程应力 = 力 / 原始截面积；工程应变 = 伸长量 / 原始标距长度。这种表达方式假设样条在整个拉伸过程中截面积不变，与实际情况存在偏差。

CAE仿真软件（以LS-Dyna为例）使用的则是有效应力应变曲线，这条曲线需要满足两个条件：一是真实反映材料在大变形阶段的应力-应变关系；二是曲线形态必须单调递增，以便于数值计算。因此，从工程曲线到有效曲线需要经过两次数学转换。

2.2 第一次转换：工程曲线→真实曲线

真实应力与工程应力的转换公式为：

真实应变与工程应变的转换公式为：

这一转换的本质是引入瞬时截面积的概念。当材料被拉伸时，样条的截面积随着变形而减小，因此真实的应力值实际上高于按原始截面积计算的工程应力值。转换后的真实应力应变曲线已经呈现出单调递增的形态。

2.3 第二次转换：真实曲线→有效曲线

在塑性大变形分析中，有效应力应变曲线采用等效应力的概念进行计算。对于单轴拉伸情况，有效应力与真实应力之间存在以下关系：

经过这两次转换得到的有效应力应变曲线，才能真正作为LS-Dyna等仿真软件的输入数据使用。

三、聚合物的高非线性与粘塑性力学特征

与具有明确晶格滑移机制的金属不同，非晶态与半结晶态聚合物（如PC, ABS, PP）的变形源于高分子链段的滑移、解缠结与取向。这种微观机制使得聚合物在宏观上表现出极其复杂的力学特征：强烈的静水压敏感性（拉压屈服不对称，压缩屈服强度往往远高于拉伸）、显著的粘弹性/粘塑性耦合响应、极低应变率下的颈缩后冷拉（Cold Drawing）现象，以及伴随微裂纹（Crazing）与剪切带（Shear Banding）竞争的损伤演化。

在构建聚合物材料卡片时，传统的金属本构模型完全失效。工程界目前倾向于采用两类策略：

第一类是基于Drucker-Prager或Mohr-Coulomb这类原本用于岩土材料的屈服准则，通过引入静水压力项来修正拉压不对称性；

第二类则是采用专为聚合物开发的半解析模型，如SAMP-1（Semi-Analytical Model for Polymers）。SAMP-1模型允许用户直接输入单轴拉伸、单轴压缩、双轴拉伸及纯剪切四条不同应力状态下的屈服曲线，并根据加载路径自动插值构建动态的三维屈服面。

同时，该模型引入了非关联流动法则（Non-associated Flow Rule）以准确控制聚合物在塑性变形过程中的体积膨胀（即泊松比随塑性应变的变化），这对于精确预测塑料扣位的插拔失效与手机外壳的跌落开裂至关重要。

四、材料卡片应用逻辑

各大主流商业求解器（如LS-DYNA, Abaqus, PAM-CRASH）在底层动力学积分算法上殊途同归，但在材料卡片的关键字定义、输入逻辑与容错处理机制上存在显著差异。深入理解这些差异，是避免“垃圾输入，垃圾输出”（Garbage In, Garbage Out）的关键。

LS-DYNA：显式碰撞仿真的材料库开发经验

LS-DYNA以其在强非线性、大变形和瞬态冲击领域的绝对统治力而闻名。其内置材料库极其庞杂，包含逾三百种材料模型，每一种模型都针对特定的物理场景进行了底层代码级的极速优化。

*MAT_024 (MAT_PIECEWISE_LINEAR_PLASTICITY)： 这是整个汽车工业与消费电子防撞设计中使用频率最高的通用弹塑性卡片。它的魅力在于卓越的鲁棒性。工程师可以通过LCSS表格输入离散的真实应力-真实塑性应变曲线族来定义应变率效应。实践经验警告： LCSS表格的插值机制是基于自然对数插值的。如果输入的应变率曲线出现交叉（即高应变率下的应力低于低应变率下的应力），或者硬化曲线呈现负斜率（未激活损伤模块时），求解器的材料刚度矩阵将出现非正定，导致不可控的网格畸变。此外，必须通过外推确保表格覆盖到极高应变率（如10000 /s），以防求解器在局部高变形区发生错误的常数外推。
*MAT_081 (MAT_PLASTICITY_WITH_DAMAGE) 与 MAT_SAMP-1 (MAT_187)： 在涉及破坏仿真的塑料件（如保险杠、手机外壳）中，MAT_187是业界公认的最佳实践。SAMP-1模型内部的屈服面随拉压状态动态改变，这要求输入的数据必须在原点保持严格的相切连续性。如果单轴拉伸与纯剪切曲线在微小塑性应变处的模量不匹配，求解器会在极短时间内由于屈服面不封闭而崩溃。
*复合材料损伤 MAT_054 (MAT_ENHANCED_COMPOSITE_DAMAGE)： 针对碳纤维层合板（CFRP），MAT_054利用Chang-Chang失效准则分别判断基体与纤维的拉压破坏。由于复合材料的极度脆性，单元失效极易引发应力波的虚假反射。工程实践中，必须精细调节DFAIL（失效应变控制）与SOFT（软化系数控制）参数，同时强制约束单元的最小破坏时间步，以防止仿真因为局部高频振荡而中止。

Abaqus：从隐式非线性到用户子程序的深度定制

Abaqus采用极其模块化的*MATERIAL关键字树状结构，使得多物理场耦合特性的定义更加符合人类直觉。

金属可塑性与损伤演化： 在Abaqus中定义延性金属，往往通过组合*ELASTIC, *PLASTIC, *RATE DEPENDENCE, 以及延性损伤起始*DAMAGE INITIATION, CRITERION=DUCTILE。Abaqus的强大之处在于其引入了断裂能正则化机制*DAMAGE EVOLUTION, TYPE=ENERGY。这一机制彻底改变了传统材料卡片随网格尺寸变小而急剧变“脆”的网格敏感性缺陷，使得能量耗散成为一个相对客观的物理不变量。
UMAT / VUMAT 的二次开发： 当标准材料库无法覆盖新兴材料（如具有形状记忆效应的镍钛合金、相变诱发塑性的TRIP钢、或者超高周疲劳退化材料）时，最高阶的仿真工程师必须依赖Fortran或C++编写用户自定义材料子程序（UMAT用于Abaqus/Standard隐式求解，VUMAT用于Abaqus/Explicit显式求解）。编写UMAT的核心挑战并不在于屈服函数的编码，而在于推导高度复杂的“一致切线刚度矩阵”（Consistent Tangent Modulus）。隐式非线性求解严重依赖该矩阵进行牛顿迭代，如果切线刚度推导存在微小误差，将导致模型在屈服点附近彻底丧失二次收敛性（Quadratic Convergence），陷入无尽的迭代发散死循环。

材料卡片定制