LS-DYNA材料模型实战解析：从选型到参数标定

LS-DYNA的材料模型，说它是整个软件的灵魂一点也不夸张。搞过几年冲击、碰撞、跌落仿真的人都有这种感觉：网格能画好、接触能设对，最后算出来结果像不像样，七八成的功劳或者锅都在材料本构上。这玩意儿不像边界条件那么直观，也不像求解器设置那样按几个按钮就完事，它需要你真正理解材料在什么状态下怎么变形、怎么失效，还得会用关键字卡片把这种物理行为翻译给求解器听。

我这些年断断续续攒了不少关于LS-DYNA材料本构的笔记，踩过的坑比走过的路还多。从最早拿MAT_001算啥都算不准，到后来给橡胶密封件标定超弹性参数，再到混凝土在爆炸载荷下开裂失效的模拟，每一步都离不开对材料模型本身的理解。这篇笔记不是教科书式的罗列，是我把平时翻手册、调参数、对比实验数据时最实用的那部分经验整理出来，希望能给正在研究材料卡片、对着Help文档一头雾水的朋友一些实在的帮助。不管是刚入门的仿真工程师，还是做过几个项目但总觉得材料参数没底的老手，这篇内容应该都有值得你记一笔的地方。

1. 先从材料库的整体地图说起

1.1 LS-DYNA材料库到底有多大

很多人第一次打开LS-DYNA的Material菜单时是懵的，二百多种材料模型，从最简单的线弹性到复杂的炸药爆轰、复合材料渐进失效、生物组织，什么都有。这不怪你，我刚入行的时候也是抓瞎，心想这不是逼人做选择题吗。

但你要真把这二百多种模型都过一遍，那是不现实的，也没必要。实际工程项目里，翻来覆去用到的核心模型就二三十种。LS-DYNA官方手册把材料模型按大类分成金属、塑料、橡胶/超弹性、泡沫、混凝土/岩土、复合材料、玻璃、炸药、黏弹性材料等若干类别。每个类别下面又有若干编号，比如金属类里有MAT_003塑性随动、MAT_024分段线性塑性、MAT_015 Johnson-Cook；混凝土类里有MAT_072 RHT、MAT_084 Winfrith、MAT_159 CSCM等等。

我的建议是别想着把所有材料模型都弄明白，先把你的主攻方向对应的那十几个模型吃透。做汽车碰撞的就盯着金属和塑料的应变率相关模型，做消费电子跌落的就重点关注泡沫和塑料类，做土木防护的就把混凝土和岩土的几个模型搞扎实。每一个模型，不仅要知道它的关键字卡片怎么写，还得知道它内部是什么屈服准则、什么强化模型、怎么处理应变率效应、怎么定义失效。这些东西不搞清楚，参数填得再漂亮，算出来也是一堆没有物理意义的数字。

1.2 几大类材料模型的选型逻辑

选材料模型，第一步不是找最先进的，而是找最匹配你需求的。我的经验是，先问自己三个问题：第一，这是个什么问题，准静态还是动态，要不要考虑应变率；第二，材料的变形量有多大，是小变形弹性问题还是大变形塑性流动；第三，要不要考虑失效，失效的判据是什么。

第一类问题决定你要不要用带率相关效应的模型。如果是冲击碰撞，材料在高速变形下表现出的应力会比准静态高很多，这时候还拿一个不考虑应变率的弹性模型去算，结果肯定不对，MAT_024和MAT_015这类带Cowper-Symonds或Johnson-Cook率项的本构就派上用场了。准静态结构分析则用简单一点的模型就够了，没必要非上复杂的。

第二类问题决定你要不要用超弹性或高精度塑性模型。橡胶、硅胶这类大变形材料要用Mooney-Rivlin、Ogden等超弹性模型，不能用小变形线弹性。金属材料如果变形量不大，MAT_001线弹性也不是不能用，但一旦进入塑性，它的应力状态就完全失真了。

第三类问题决定你要不要关注失效定义。LS-DYNA的很多金属模型都自带失效准则，比如MAT_015里用失效应变，MAT_024里用塑性应变失效，MAT_020里带Gurson-Tvergaard-Needleman损伤模型来模拟金属的韧性断裂。混凝土模型里还区分了拉伸损伤和压缩损伤，比如MAT_159的FEA混凝土模型。

说到底，选型就是把你对材料物理行为的理解翻译成关键字卡片上那几行参数。材料模型不是越高级越好，而是越贴合你的问题越好。高模型如果参数标不准，反而不如中等的模型配合扎实的参数来得可靠。

2. 最常用的几类材料模型实操解析

2.1 金属材料：从弹性到Johnson-Cook

金属材料是我做碰撞仿真时打交道最多的。先说最基础的MAT_001 *MAT_ELASTIC，这个模型就是最简单的线弹性，需要定义弹性模量E、泊松比NU和密度RO。做非常简单的校核计算、或者模拟不会进入塑性的结构件时用它，其他时候基本用不上。这里的坑在于，如果你不确定结构的受力是否超过屈服强度，直接用线弹性算出来的应力分布可能云图上好看，但数值已经超出了材料的物理承载极限，成了几何上的应力集中数字游戏。

真正做金属件冲击碰撞，我主力用的是MAT_024 *MAT_PIECEWISE_LINEAR_PLASTICITY和MAT_015 *MAT_JOHNSON_COOK。MAT_024的精髓在于可以用一条真实的应力-应变曲线来定义塑性段，同时支持Cowper-Symonds应变率模型。关键字卡片里要填材料密度、弹性模量、泊松比、屈服应力，然后用DEFINE_TABLE定义多根不同应变率下的应力-应变曲线。它的优势是直观，实验测到什么曲线就输入什么曲线，不用做太多数学拟合，工程上很好用。

MAT_015式Johnson-Cook模型则是把流动应力写成了一个关于塑性应变、应变率和温度的函数。公式是这样的：

$$\sigma_{y} = (A + B \varepsilon_{p}^{n})(1 + C \ln \dot{\varepsilon}^{*})(1 - T^{*m})$$

A是初始屈服应力，B和n是应变强化系数，C是应变率敏感系数，m是热软化指数。这个模型的好处是只要A、B、n、C这几个常数标定好了，整个应力状态都能用一个公式描述，不需要一大堆曲线。缺点是参数的获取有门槛，通常需要通过不同应变率下的拉伸或压缩实验拟合。我见过不少初学者不知道Johnson-Cook里的A和B都是真实应力而非工程应力，直接拿材料手册上的抗拉强度填进去，结果算出来的变形模式完全不靠谱。

3. 这里提醒一句：不管哪个金属模型，输入的应力应变曲线一定要是真实应力-真实应变。把工程应力应变数据直接给求解器等于给自己埋雷，后面会专门讲转换方法。

2.2 超弹性和黏弹性材料：橡胶密封圈的活路在哪

做跌落、密封、减振类项目时，橡胶和软胶材料是主角。这类材料的特点是可恢复的大变形，本构关系不再是线性的，而是基于应变能密度函数来描述的。

LS-DYNA里橡胶材料常用的有MAT_004 *MAT_ELASTIC_PLASTIC、MAT_027 *MAT_MOONEY-RIVLIN、MAT_077 *MAT_HYPERELASTIC_RUBBER、MAT_181 *MAT_SIMPLIFIED_RUBBER等。其中最经典的是Mooney-Rivlin模型，它把应变能密度表示成两个不变量，关键参数是C10和C01，两个剪切响应系数。做单轴拉伸、等双轴拉伸、平面剪切三种实验，把数据拟合一下，C10和C01就有了。

我个人的习惯是能用MAT_077就用MAT_077。它内部集成了多种超弹性本构，包括Mooney-Rivlin、Ogden、Yeoh等，配合实验数据可以自动拟合，而且对不可压缩行为的处理更好。很多新手容易忽略橡胶材料的体积响应，直接设置一个很大的体积模量，结果计算中会出现体积锁死或者极度小的步长。实际工程中橡胶处理成近似不可压缩，泊松比设在0.495附近，并用B-bar或Uzawa算法配合，能显著减少体积锁死。

有一点要注意，超弹性参数是有有效范围的。你用单轴拉伸拟合出来的C10、C01，换到压剪主导的工况时可能会失真。所以我做密封结构模拟时，通常会把单轴、等双轴、纯剪切三类实验数据一起用来拟合成套参数，这样得到的参数适用范围更宽。

对阻尼减振橡胶，还要考虑黏弹性，这时候要用MAT_076 *MAT_GENERAL_VISCOELASTIC。它用Prony级数描述松弛模量，简单理解就是用一系列弹簧-阻尼器串联并联来模拟橡胶的应力松弛和蠕变行为。这个模型需要频率依赖的动刚度数据或时间依赖的松弛数据，做起来更麻烦，但对减振类项目精度提升很大。

2.3 混凝土与岩土：失效形态怎么还原

混凝土和岩土材料，在LS-DYNA里的选择也不少。做爆破拆除、冲击防护这些强动力问题，我经常用的有MAT_072R3 *MAT_CONCRETE_DAMAGE_REL3、MAT_084 *MAT_WINFRITH_CONCRETE和MAT_159 *MAT_CSCM_CONCRETE。这三个是混凝土模拟的主流选择。

MAT_072R3是LS-DYNA里名气很大的一个，因为它自动生成参数：只需要给质量密度RO和单轴抗压强度，求解器就能基于这些输入生成完整的屈服面方程。这对没有太多实验数据的场景很友好，特别适合军方和防护工程领域做快速评估。但它的参数内部是压垮性标定的，如果你把抗压强度从40MPa改成80MPa，就出声了，因为失效面是缩放函数，不能做简单的直接改值处理，我遇到很多朋友在这个模型上翻车，原因都是随意改了单轴强度而没有调整其他参数。

MAT_159 *MAT_CSCM是基于连续面帽模型的混凝土本构，它模拟混凝土的硬化、软化、损伤和率效应，对拉伸开裂和压缩破碎的区分做得比较细。它同样可以只输入抗压强度来生成默认参数，但更方便用户根据实际实验调整。在穿透分析和混凝土靶板冲击模拟里，MAT_159的表现很不错，因为它考虑到了拉压不对称、静水压力对强度的影响和软化下降段。

岩土方面，莫尔-库仑、Drucker-Prager这些经典准则在LS-DYNA里有对应的MAT_005 *MAT_SOIL_AND_FOAM和MAT_016 *MAT_PSEUDO_TENSOR。纯土体模拟的关键是状态方程EOS配合压力相关屈服面，不能孤立看材料参数。说实话，岩土材料的模拟比金属更依赖经验和标定，不要指望有一个万能的材料模型解决所有土体问题。

混凝土和岩土模型最大的共同坑点是：网格尺寸对结果的影响非常大。软化段参数如果不做网格尺寸正则化，同一个材料在5mm和20mm网格下算出来的破坏模式可能判若两人。LS-DYNA的MAT_159内部自带了一定的网格正则化处理，这也是我偏好它的原因之一。

3. 参数标定、单位制与曲线处理，这三关过不了就别谈精度

3.1 单位制：最先要过的一关

每次看到有人在论坛上问“为什么我的模型应力大了1e6倍”，九成是单位制出的问题。LS-DYNA本身不强制单位，你只要保证输入参数的单位自洽就行。但LS-DYNA的官方材料参数很多用的是kg-mm-ms的单位制，这和很多工程师习惯的kg-m-s或者g-cm-us差别很大。

常用的三套单位制对比我直接列出来：

我最推荐的是mm-ms-kg这套，因为它对应的应力单位是GPa，数值上比较好读。很多材料手册上的参数就是按这个单位制来的，比如钢的弹性模量210GPa，填E=210就行。密度的话，钢材填7.85e-9（单位是Mg/mm³，等同于kg/mm³）。

单位制搞错，轻则结果数量级不对，重则求解器直接报出“zero or negative pressure”之类的错误。我的习惯是在每个模型文件的最前面加一个注释卡片，写清楚当前文件使用的单位制，免得隔几个月回头再看时忘了当初用什么单位建的模型。

3.2 工程应力应变曲线怎么转换成真实曲线

给MAT_024这类模型输入塑性曲线之前，必须先把实验测得的工程应力应变数据转成真实应力真实应变。公式其实很简单：

真实应力 = 工程应力 × (1 + 工程应变) 真实塑性应变 = ln(1 + 工程应变) - 真实应力 / 弹性模量

这里有个容易忽视的地方，输入的是真实塑性应变，而不是真实总应变。如果你直接拿真实总应变输进去，相当于重复计入了弹性部分，计算出来的屈服面从一开始就是歪的。我见过好几次这种情况，模型算出来的应力水平比实验低不少，排查到最后发现是塑性应变没有扣掉弹性段的贡献。

还有一个经验，输入的塑性曲线应该在应变率范围内尽量覆盖实际工况，尽量不要在0应变率下只给一条曲线。MAT_024支持定义多条不同应变率下的曲线，LS-DYNA会在两条曲线之间对数插值。如果只有一条准静态曲线，动态冲击时材料应变率硬化效应完全体现不出来，结果会偏软。实在没有高应变率实验数据，可以先用Cowper-Symonds参数做一个修正。

3.3 参数从哪里来：实验、文献与经验的优先级

材料的参数来源，我的优先级排序是这样：

第一是自有实验数据，这个最靠谱，但也很贵。冲击仿真需要的高应变率数据通常要做霍普金森杆实验，不是每个项目都做得起。第二是公开发表的文献数据，尤其是那些做了系统实验研究并给出完整参数表的论文，可以参考其测试条件和适用场景。第三是材料手册和LS-DYNA官方范例中的传感器库，LS-DYNA官方手册里面的材料参数很多来自汽车工业多年积累，比如各类钢的Johnson-Cook参数在文献里能找到多种版本，注意选用与你材料牌号及热处理状态一致的。

第四是我自己积累的项目数据，我平时做项目时会顺手把每次实验标定出来的参数记在一个个人材料库里，包括材料名称、来源、适用应变率区间、单位制、实测曲线文件路径。后续建模时直接调用，效率高得多，而且经过实际项目验证的参数比抄来的参数可靠得多。

这里分享一个教训：不要从别处的卡片文件里直接复制材料参数而不去核对。我有个同事，从百度文库下载了一个含MAT_015的K文件，里面A、B、C等系数都齐得很，他直接拿来做铝合金冲击仿真，结果算出来的变形和实验对比差将近40%，后来发现那个参数根本就是某种钢的Johnson-Cook数据，被上传的人误标为铝合金。做材料仿真，参数溯源和验证永远不能省。

4. 材料相关的常见报错与排查技巧

4.1 负体积（Negative Volume）怎么排查

做LS-DYNA经常遇到“negative volume in solid element”的报错。报错之后求解器可能继续跑也可能直接中断，不管怎样，这个模型已经处于不可信状态了。材料相关的原因占据负体积报错的大头，金属单元在极端压应力下屈服面收缩到零，橡胶泡沫单元发生过度压缩，都会产生负体积。

排查步骤我一般是这样：

第一步，先看报错出现在什么位置，是哪个Part，对应什么材料。有些负体积出现在应力集中区域，那很可能是网格问题而非材料问题，加密网格或重新划分就能解决。

第二步，看该位置的变形程度。如果变形量远超材料的实际物理极限，比如橡胶被压缩到原有厚度的十分之一以下，那说明材料参数里的体积响应不对。像泡沫材料MAT_063 *MAT_CRUSHABLE_FOAM，它的压强-体积应变曲线决定了材料在压溃过程中的承载能力，如果曲线输入不合理，材料会被压成负体积。

第三步，检查是否使用了正确的单元算法。超弹性材料、泡沫材料对单元算法非常敏感，推荐用全积分单元或带B-bar修正的算法，否则容易沙漏和锁定并存。我自己用MAT_181和MAT_063时，会优先选择1点积分ALE或全积分算法，再配合刚度沙漏控制。

第四步，排查接触参数。材料过软或过度变形时，默认接触刚度可能出现穿透，特别是软材料和硬材料接触的场景。这时调整SOFT参数为1或2，并对软材料Part单独设置接触刚度，能缓解大部分负体积问题。

4.2 沙漏能超标：材料性质可能背了锅

沙漏（Hourglass）是LS-DYNA单点积分单元的典型病，表现为网格出现锯齿状形变，模型的能量不守恒，沙漏能占比过大时结果就是废的。材料模型虽然不能完全决定沙漏，但确实有关系。

像橡胶、泡沫这样刚度较低的材料，在单点积分单元中抵抗沙漏的能力弱，比刚度大的金属更容易出现沙漏。陶瓷材料由于本身刚度响应快，如果设置全部单点积分单元，也很容易在冲击波传播阶段出现局部沙漏。

控制沙漏的办法我按优先级排序：优先用全积分单元，代价是计算量增大；其次用刚度形式沙漏控制，设置IHQ=2并给合适的QM参数；再者细化网格，从物理上消除过大变形引起的沙漏模式。

这里有一个很实在的经验：沙漏能占总内能的比例一般控制在5%以下才算健康，冲击问题可以稍微放宽到10%以内。如果超了，先看是不是网格太粗，再看是不是材料刚从低刚度切换到大刚度产生数值震荡。有些材料模型自带黏性阻尼或横向剪切修正，比如MAT_181内部对沙漏的抵抗就要比MAT_001好不少，这也是我偏爱用专业材料模型的原因。

4.3 压强、体积模量和不可压缩性的坑

不少新手在定义超弹性材料时，不知道体积模量该怎么填，就随便填了一个大数。这往往是灾难的开始。体积模量填得过大，材料会表现出近似刚性体积响应，时间步长被拉得非常小，计算变得极慢，甚至出现压强震荡。

实际上，超弹性材料有一个独立的体积响应定义方式。以MAT_077为例，它允许单独定义体积模量K，也可以用一个EOS来描述压力-体积变化关系。如果用了EOS，材料卡片里的体积模量就不能再过度设定了，二者配合不好会出现双体积响应的冲突。

我踩过一个印象很深的坑：做橡胶密封条压缩时，为了追求不可压缩性，把体积模量设成了默认钢材的数量级，结果计算时间暴增到原来的十倍以上，压缩力也比实验值大了一倍。问题根源就在于体积模量过大，导致材料在体积变形上过于“刚”，压缩时压不下去，应力全部堆在了接触面上。后来按橡胶的真实体积模量大约在2GPa左右设置（相对剪切模量大两到三个数量级），计算速度恢复正常，力值也对上了。

对金属材料，体积模量通常不是直接输入的，弹性模量和泊松比已经隐含了体积模量，不需要单独特意设置。只有谐波要格外检查：钢板在爆炸冲击下，压力波在材料内部传播，如果体积模量和密度换算出的声速与物理不符，波传递的时间序列就会出偏差，最终应力波到达时间和反射波形全乱。

5. 给新人的几条实战建议

5.1 从案例库开始，不要从零自建

LS-DYNA安装包里有大量的官方验证案例和示例K文件，很多材料模型都有对应的验证算例。我强烈建议新手不要自己凭空想一个工况去试材料，而是去找官方验证案例，把里面的材料参数提取出来，对照手册理解参数含义，再改造成自己的工况。

官方材料卡片的参数为什么可靠？因为它们和实验对照过，很多是关键标准实验的复现。比如LS-DYNA验证手册里对MAT_024、MAT_015、MAT_159都有专门案例，跑通这些案例你就能看到材料模型在不同受力状态下的标准表现，这个认知比看十篇教程都有用。

我在线上不是学院派的，当年入行时就是靠着反复跑官方案例，逐步摸索出每个参数对结果的影响规律。现在带新人，我让他们做的第一件事就是跑通材料案例库，不是用鼠标点两下GUI，而是手动打开K文件，逐行读卡片，理解每个参数对应的物理含义。

5.2 材料卡片的版本兼容性

LS-DYNA版本升级后，有些材料模型的关键字卡片格式会变化。比如老版本里MAT_024的某一位参数在R8系列里增加了新选项，或者默认值变了，如果你拿着老K文件直接提交到新版本求解器，可能不会有报错，但结果可能已经悄悄改变了。

这提醒我：保存K文件时，养成记录LS-DYNA版本号的习惯。我通常在K文件头部注释里写上“Generated with LS-DYNA R12.0 / tested on R11.0”，这样后续换版本时方便对照。另外，升级版本后最好用相同算例跑一遍，对比一下关键结果，确认没有偏离再进入正式项目。

LS-DYNA在不同版本中对材料模型的支持也有差别，个别新引入的模型如MAT_297 *MAT_SAMP-1等只在较新版本中可用，老版本里根本不识别。如果你的团队还在用老版本求解器，看到没见过的材料编号时，先确认版本支持再花时间研究。

5.3 材料模型的在线资源和社区

除了官方手册，LS-DYNA的学习资源还包括官方技术论文、每年的LS-DYNA用户大会论文集、以及各大论坛里工程师分享的实战经验。我在做材料标定时遇到各种稀奇古怪的报错，很多都能在社区里找到同类问题。只是要注意，社区里的参数和结论要谨慎使用，最好带着批判眼光验证一下，特别是那些“亲测有效”的帖子，没有说明工况和边界条件的参数，照搬很可能翻车。

LS-DYNA官方的材料模型指南文件，几乎每年更新，里面有每个材料模型的详细介绍、适用场景、参数说明甚至推荐单元算法。我常把这份PDF当词典查，比在Help里搜索高效得多。新版本发布时，我会重点关注新增了哪些材料模型、哪些模型更新了算法，这些往往反映了行业前沿需求。

说实话，做LS-DYNA材料这块，最大的门槛不是学会读关键字卡片，而是建立起参数和物理行为之间的映射感知。你得看到某组参数，就能大概想象出这条材料在冲击下会怎么表现；建模出了问题，能快速判断是材料参数的问题、网格的问题，还是接触算法的问题。这种感知没有捷径，就是一点点积累，多算、多对比、多复盘。

我平时做项目时有个习惯，每做一个项目就把材料标定过程中的心得、参数、曲线和遇到的问题都记在一个本子里，包括参数来源、实验数据文件放哪、试过哪些组合、为什么最终改成现在的值。时间长了，这个笔记本身就成了一个很大的资料库。回头看这篇文章里提到的一些思路，其实都是从这个本子里翻出来的。材料模型的积累没有终点，每个新项目都可能推翻你对某个参数的理解。保持记录，保持怀疑，保持对实验数据的敬畏，这条路才能越走越稳。

最后再分享一个我摸索出来挺管用的小技巧：拿到一组新的材料参数，先别急着跑完整大模型，单独建一个单单元模型或者一个小立方体模型，施加几种简单的载荷工况（单向拉伸、纯剪切、单轴压缩），观察单元在这些载荷下的响应是否合理。这一步相当于给材料模型做一次体检，花不了多少计算资源，却能帮你提前排除绝大部分参数错误。正是这样一次次的“体检”，才让复杂工况下的仿真结果有了底气。