简介本资源是一份面向LS-DYNA与AUTODYN用户的专业级材料本构模型参数汇编适用于冲击动力学、侵彻仿真及弹体/破片毁伤效应建模等高阶工程仿真场景尤其适配装甲材料、穿甲弹体、战斗部破片等典型弹药结构的数值模拟需求。压缩包仅含1个HTML文件32KB内容为结构化整理的材料参数表涵盖多种常用金属、复合材料及高能材料的密度、屈服强度、硬化模型、失效准则等关键本构输入参数并附有AUTODYN与LS-DYNA双平台兼容性说明及典型MAT号索引便于工程师快速调用与校验。目前已有1733人学习下载内容源自工程实践验证可直接嵌入仿真前处理流程显著减少材料卡片搭建时间与参数试错成本是从事防护结构、武器效能仿真等方向研究人员与CAE工程师的实用速查工具。 做显式动力学仿真的朋友估计都经历过这种时刻翻遍整个材料库对着屏幕上的关键字卡片发呆不知道手头这个模型到底该选哪个材料号也不知道网上下载的那份“材料参数大全”能不能直接抄。尤其当你同时在用LS-DYNA和AUTODYN的时候两边材料库的组织方式、单位制、参数命名完全不一样稍不留神就得出一个看起来正常、实际全错的仿真结果。这篇内容围绕LS-DYNA和AUTODYN里的本构模型与材料参数展开讲讲我这些年攒下的材料模型选型思路、参数整理方法和踩坑记录希望对正在搞碰撞、冲击、爆炸、成形类仿真的朋友有点帮助。先说结论所谓“材料参数大全”真正值钱的不是那张参数表本身而是你拿到参数之后能不能判断“这组参数在什么条件下有效、适合什么模型、要怎么校验”。我会把常见材料模型按应用场景拆开讲再解释参数从哪来、怎么整理、怎么验证最后重点说说LS-DYNA和AUTODYN两套软件之间参数互通的坑。无论你是刚入门的新手还是已经跑过不少模型的老手这部分内容应该都能让你少走点弯路。1. 先搞清楚“材料模型”和“材料参数”到底是什么关系1.1 本构模型是数学框架材料参数是具体数据我经常看到有人把“材料模型”和“材料参数”混成一谈这其实是个容易出问题的误解。本构模型是描述材料力学行为的数学框架它决定了一组应力、应变、应变率、温度之间怎么互相影响而材料参数则是填进这个框架里的具体数字。比如Johnson-Cook模型有A、B、n、C、m这些参数A是参考应变率下的初始屈服应力B和n控制硬化行为C控制应变率敏感性m控制温度软化。这些参数离开了Johnson-Cook这个框架就没有意义同样换了框架之后材料参数必须重新标定不能直接平移。类比一下本构模型就像一辆车的底盘和悬挂结构材料参数就是具体调校数据。同一种钢材你可以用双线性弹塑性模型描述也可以用Johnson-Cook模型描述得到的仿真精度和适用范围完全不同。“哪个模型更高级就用哪个”这个思路在工程里不可取真正的选择依据是你关心的物理现象是什么、你要覆盖的应变率范围有多大、你手里有什么数据可以标定。1.2 显式求解器里的材料卡到底由哪几块组成在LS-DYNA里一张完整的材料卡通常分成几层密度和弹性参数如弹性模量、泊松比负责基础力学响应强度模型计算屈服应力和硬化状态方程EOS计算压力和体积变形关系失效准则决定单元什么时候删除或损伤开始。AUTODYN的经典材料模型则经常直接打包成“状态方程强度模型失效模型”三件套逻辑上更加清晰。这个分层关系直接决定了你整理参数时该怎么归档。我见过很多人只记了密度、弹性模量、屈服强度和几个硬化参数结果到了需要跑爆炸冲击、超高速碰撞这类高压工况时发现缺了状态方程参数只能临时找文献找到的EOS参数和强度模型参数还是不同来源、不同单位制的硬拼在一起结果自然没法看。所以一旦涉及高压、高应变率问题从源头就要按状态方程、强度、失效三块去收集参数。2. 常用本构模型按场景拆解选对模型才算入门2.1 金属材料从PLASTIC_KINEMATIC到Johnson-Cook金属材料是显式动力学里用得最多的一类。如果只做低速碰撞、成形分析且没有特别大的塑性变形用LS-DYNA里的*MAT_PLASTIC_KINEMATICMAT_003就足够了。这个模型本质是双线性弹塑性加Cowper-Symonds应变率修正关键参数有密度RO、弹性模量E、泊松比PR、初始屈服应力SIGY、切线模量ETAN和应变率参数C、P。这里有一个高频坑点ETAN填的是塑性切线模量也叫硬化模量不是弹性模量。假设某钢材弹性模量210GPa、屈服强度350MPa、抗拉强度600MPa、均匀延伸率10%估算切线模量大约在2.5GPa到3GPa左右很多人直接填了20GPa甚至200GPa结果模型刚度离谱偏大变形完全不对。正确做法是用真实应力应变曲线硬化段的斜率最好用实验数据拟合。如果问题涉及高应变率、高温或者大变形破坏Johnson-Cook模型MAT_JOHNSON_COOKMAT_015是业界标杆。它的表达式分成三项应变硬化项ABε^n、应变率项(1Clnε)、温度项(1-T*m)。注意A、B的单位和弹性模量一致n、C、m无量纲参考应变率通常取1/s或0.001/s要和标定条件对应。JC模型还带D1到D5五个断裂参数用来描述应力三轴度、应变率、温度共同作用下的失效应变。但在实际项目里Johnson-Cook模型也有它的毛病。它的应变率项是对数形式在10^4/s以下表现尚可到了10^5/s以上外推风险很大。另外它假设温度软化是单一幂次形式对某些铝合金、钛合金并不完全适用。这时候可以考虑LS-DYNA里的*MAT_MODIFIED_JOHNSON_COOKMAT_107或者分段表格化的JC变体直接把实验得到的应力应变曲线作为表格输入同时保留应变率和温度插值精度会明显提升。2.2 混凝土与岩土材料别指望一组参数吃遍所有标号混凝土材料的难点在于抗压/抗拉不对称、压力相关、应变率相关、损伤软化这些特性远不是弹塑性模型能描述的。工程上常见的有LS-DYNA的*MAT_CONCRETE_DAMAGE_REL3MAT_072R3、*MAT_CSCMMAT_159、*MAT_JOHNSON_HOLMQUIST_CONCRETEJHC和AUTODYN里常用的RHT模型。MAT_072R3是KC模型的一个简化版好处是只需要输入密度、泊松比、无侧限抗压强度和几个单位转换系数LS-DYNA可以自动生成剩余参数。它生成默认参数的过程依赖一个内部公式库适用于普通强度混凝土对高强混凝土、纤维混凝土或者含骨料特殊混凝土得谨慎。CSCM模型自带帽盖屈服面和损伤参数对隧道、地基、混凝土结构受撞问题比较合适但它的参数输入格式经常让新手摸不着头脑很多参数是可选的不填时用默认值问题是默认值未必适合你的单位制。如果你接触过AUTODYN会发现RHT模型在AUTODYN材料库里已经预置了很多岩石、混凝土的参数使用非常方便。但要留意这些预置参数不少来自于公开文献标定可能基于特定实验条件直接用于你的工况时必须做敏感性分析。混凝土这类材料每次仿真的核心问题是你的模型标定过哪个抗压强度等级C30和C50的参数差很远用C50的参数模拟C30构件应力峰值可能偏差30%以上。岩土材料则常用Mohr-Coulomb、Drucker-Prager以及LS-DYNA里的*MAT_FHWA_SOILMAT_147这类模型涉及黏聚力、内摩擦角、含水率、孔隙水压力等参数。这些参数一般来自土工实验没有实验数据时不要拍脑袋填否则连边坡稳定的定性趋势都可能算反。2.3 橡胶、泡沫与超弹性材料曲线比参数表更值钱橡胶、密封件、缓冲泡沫这类材料在汽车碰撞、包装跌落、防护结构中很常见。它们的核心特点是高度非线性弹性甚至在大变形下表现出体积近似不可压缩。LS-DYNA里的Mooney-Rivlin模型MAT_027和Ogden模型MAT_077适合橡胶类低密度泡沫用MAT_LOW_DENSITY_FOAMMAT_063可压碎泡沫用MAT_CRUSHABLE_FOAMMAT_026。我的建议是这类材料不要只抄几个弹性常数最好拿到完整的单轴压缩应力应变曲线以表格曲线方式输入。泡沫材料的迟滞、能量吸收、残余变形行为都体现在曲线形状里光靠几个材料号根本覆盖不了。低密度泡沫还要注意设置合适的阻尼系数DAMP否则卸载阶段可能出现剧烈的数值振荡这个系数本身没有实验意义就是数值稳定参数但调不好轻则沙漏能失控重则负体积报错。2.4 状态方程高压冲击问题里的另一半关键很多人做低速碰撞时从头到尾不碰状态方程但一旦涉及爆炸、高速冲击、水压冲击、航天器防护这类问题EOS就是决定成败的那一半。LS-DYNA常用的*EOS_GRUNEISEN适用于金属和水的冲击绝热关系*EOS_LINEAR_POLYNOMIAL适用于气体和一般流体*EOS_JWL用于炸药爆轰产物。Gruneisen EOS的核心参数包括C、S1、S2、S3和Gamma0。C是材料在零压力下的体声速S1等是冲击波速度与粒子速度关系式的斜率系数Gamma0是Gruneisen系数。以水为例C大约在1480m/s量级但在不同单位制下要换算成对应单位。很多人在低速流固耦合里也强行引入Gruneisen EOS其实没有必要线性多项式EOS反而更稳反过来爆炸冲击问题只用弹性模量而不填EOS那就完全没法计算冲击波传播。3. 参数从哪来如何积累一份真正能用的“参数大全”3.1 官方材料库和文档是第一站但别迷信LS-DYNA的LS-PrePost自带材料参数库打开界面可以看到按金属、非金属、混凝土、泡沫等分类的预置材料直接生成关键字卡片很方便。这些数据作为起点完全够用尤其是帮助你快速跑通模型。但官方库的局限也很明显它提供的多是典型值缺少标定条件和适用范围说明你要在仿真报告中交代参数来源时会发现没法溯源。AUTODYN的材料库在这方面做得更系统。AUTODYN把材料分成金属、非金属、炸药、水、空气、岩土等大类每个材料条目下挂载了EOS、强度模型、失效模型三组参数并且标注了数据来源文献。我在实际工作中经常用AUTODYN材料库作为参数交叉验证的基准尤其是同一种材料在LS-DYNA里参数怎么填不确定时去AUTODYN查一下对应量纲下的数值能很快发现量级错误。3.2 文献、数据库和实验标定怎么配合材料手册和公开数据库是第二站比如ASM手册、军工材料数据库、各大学公开的论文补充材料。网上流传的各种“LS-DYNA常用材料参数Excel表”不是不能用但使用前必须做三步核查第一步确认单位制第二步确认模型名第三步确认标定条件如温度、应变率、热处理状态。很多时候同一款铝合金挤压态和T6态的参数完全不同直接套用等于白算。如果项目对精度要求高最终归宿一定是实验标定。常规流程是先做准静态拉伸实验获得弹性模量、屈服强度、抗拉强度和硬化曲线再做不同应变率下的动态实验万能试验机做低速SHPB做中高应变率拟合Johnson-Cook的C值有条件的话做不同温度下的拉伸实验拟合m值最后用Hopkinson杆或Taylor撞击实验做模型验证。这个过程看着复杂但对高精度结构冲击仿真来说是绕不开的。我见过一个反面案例有人直接从文献里抄了一组Ti-6Al-4V的JC参数算出的弹道极限比实验低了30%。原因就是文献里的参数是基于退火态钛合金标定的而实验用的是锻造态屈服强度差了一大截。参数跟材料状态强相关这句话请务必写在你的参数管理表头里。3.3 把参数整理成带溯源信息的数据库真正的“参数大全”不是一页密密麻麻的数据而是一个可以持续积累、可溯源、可复用的数据结构。我的做法是用一张总表管理常用材料每个材料条目至少包含编号、材料牌号、材料状态、本构模型名、单位制、参数列表、数据来源、标定条件、适用范围、验证状态、备注这十一个字段。单位制和来源这两个字段绝对不能省因为这是后续排查问题最快的抓手。如果你长期在同一套单位制下工作还可以把整理好的材料卡保存成K文件片段用时直接include不用每次重新填卡。我自己维护了一个“通用材料卡”文件夹里面每个K文件都有命名规范比如“steel_q235_mat3_mm-ms-g.k”找起来一目了然。这个习惯第一次养成时花不了多少时间但是后面节省的返工时间数以十倍计。4. LS-DYNA和AUTODYN参数互通单位换算与格式差异4.1 单位制换算最容易翻车的环节我接触过不少同时使用LS-DYNA和AUTODYN的团队最典型的错误就是单位制混用。LS-DYNA本身不锁单位制只要保证封闭一致就行但LS-PrePost的交互式和某些转换工具默认按特定单位制解释容易让人产生“软件里填的就是标准值”的错觉。我通常把工作中最常用的一套单位制定义为mm-ms-g体系长度mm时间ms质量g应力MPa密度g/mm³速度mm/ms数值上等于m/s。在这个体系下钢材密度填7.85e-3弹性模量填210000屈服强度填350JC参数里A、B也按MPa填。如果从AUTODYN换到LS-DYNA或者反过来需要注意AUTODYN的经典参数常以米-千克-秒-帕m-kg-s-Pa为单位给出那么从AUTODYN里导出的模量2.1e11对应的是210GPa填到LS-DYNA的mm-ms-g体系必须改成210000直接填2.1e11会导致刚度大了100万倍模型完全锁死。反之从LS-DYNA的mm-ms-g体系换到m-kg-s-Pa则要乘以1e6。再补充一套可能遇到的cm-g-μs体系长度cm时间μs质量g密度g/cm³应力单位是Mbar1Mbar100GPa。这套单位制在爆炸冲击文献里很常见JWL状态方程参数经常用这套单位给出。因此当你在文献里看到JWL的A、B参数是“3.7e11”这类大数字时先确认单位是Pa还是Mbar再决定要不要除以1000。我把常用单位制的换算关系整理成表方便直接对照单位制长度时间质量密度单位应力/压力单位典型场景mm-ms-gmmmsgg/mm³MPa结构碰撞、汽车安全m-s-kgmskgkg/m³Pa通用流固耦合、AUTODYN默认cm-μs-gcmμsgg/cm³Mbar100GPa炸药爆轰、冲击波文献mm-ms-kgmmmskgkg/mm³GPa某些大型冲击模型这个表建议截图保存或者贴在工位旁边每次填材料卡之前看一眼成本几乎为零收益非常大。4.2 模型对应的映射关系同参数不同卡片写法AUTODYN的材料模型逻辑是“EOSStrengthFailure”三件套LS-DYNA则是通过独立的材料关键字、状态方程关键字、失效关键字组合实现。虽然底层物理模型可能是一样的但卡片写法完全不同。我总结过一套常用映射关系AUTODYN里的Linear EOS对应LS-DYNA的EOS_LINEAR_POLYNOMIALShock EOS对应EOS_GRUNEISENJWL EOS对应EOS_JWLAUTODYN里的Johnson-Cook强度模型对应LS-DYNA的MAT_JOHNSON_COOKAUTODYN里的Piecewise JC对应LS-DYNA的MAT_MODIFIED_JOHNSON_COOK或者带表格的JC变体AUTODYN里的Von Mises强度模型对应LS-DYNA的MAT_PLASTIC_KINEMATIC或者*MAT_ELASTIC_PLASTIC_HYDRO。失效模型这块差异更大。AUTODYN里失效条件通常是材料模型自带的参数比如最大拉伸压力、最大主应力、累计损伤等LS-DYNA则经常用MAT_ADD_EROSION给任何材料卡追加失效准则。很多人做LS-DYNA仿真时发现材料没有删除第一反应是材料卡填错了其实很可能是忘了加MAT_ADD_EROSION或者失效准则是基于压力而实际应力状态里压力没有达到阈值。4.3 互通时的三类典型坑第一类坑是压力符号约定。AUTODYN的很多材料模型里压力为正代表压缩而LS-DYNA的某些历史和应力输出里正负号惯例可能不一样尤其是在EOS卡片里填参数时要严格看手册里对压力的定义。我不止一次见过把Gruneisen EOS的C值填成负数导致模型一开始就爆炸式膨胀。第二类坑是失效模型参数的默认单位。*MAT_ADD_EROSION里压力单位必须和当前单位制匹配比如mm-ms-g体系下填的是MPa很多人从AUTODYN的Pa单位制切换过来后忘记转换等于失效阈值放大了一百万倍永远不失效。第三类坑是材料模型编号在不同LS-DYNA版本里可能有变体。同一个Johnson-Cook物理模型老版本用MAT_015新版本可能推荐MAT_015的增强版或者MAT_107。版本升级后老K文件里某些卡片会自动映射到新编号但参数含义未必完全一致建议升级后重跑一个标定算例做回归对比。5. 拿到参数后怎么验证单单元测试与常见坑位5.1 单单元测试五分钟找出参数低级错误无论参数是网上抄来的、文献里查来的还是实验标定的我的原则是不经过单单元测试验证的参数绝不上大模型。单单元测试思路特别简单建一个单单元模型施加单轴拉伸或压缩提取应力应变曲线和输入的理论曲线或者实验曲线对比。在LS-DYNA里单单元可以用一个8节点六面体单元一端固支另一端给恒定速度速度值选得足够小以避免应力波振荡干扰然后输出单元的应力应变。运行时间通常几秒钟K文件也就几十行。在AUTODYN里也可以建一个单单元网格做类似验证。这个测试能抓出大部分低级错误密度数量级不对、弹性模量差了1000倍、屈服强度填成了抗拉强度、曲线单位错位等。我发现一个特别高效的判断方式把仿真输出的名义应力-名义应变曲线跟输入的理论曲线叠在一起看。如果曲线斜率、屈服点、硬化趋势基本吻合说明参数基本可信如果曲线形态完全不对优先检查单位制和卡片填写。千万不要指望大模型能替你发现参数错误大模型里各种效应耦合在一起参数错了也常常能跑完只是结果没人敢用。5.2 动态参数验证SHPB数据的处理不能偷懒如果标定的是Johnson-Cook这类含应变率效应的模型单单元静态拉伸还不够最好用SHPB实验数据做动态验证。SHPB实验得到的是入射波、反射波、透射波三组波形需要经过二波法或三波法处理才能换算成材料的应力应变曲线。这里有一个新人常踩的坑直接把原始波形信号当成应力应变数据塞进参数拟合程序或者从论文里截取一张应力应变图再用软件抠数据点但忽略了原始实验的应变率并不是恒定值。正确做法是把SHPB数据按应变率分组每组取平均应变率下的代表性应力应变曲线再做JC参数拟合。拟合时还要注意应变范围SHPB在高应变率下试样容易出现绝热温升导致流动应力下降这个下降其实是温度效应不是材料本身的应变率软化。要想把C和m拟合得准最好结合不同温度和不同应变率的多批次数据做联合标定而不是单条曲线硬拟合。5.3 常见报错与异常结果速查我把这些年遇到过的高频问题整理成了一张速查表按症状、原因、排查方向排列症状最常见原因排查方向模型刚度异常大几乎不变形弹性模量单位制错误多填了6个数量级核对弹性模量数值是否是当前单位制下的正确值单元过早产生负体积材料失效删除不及时或者EOS缺失/参数不合理检查状态方程参数增大失效准则触发可能或细化网格沙漏能超过内能10%使用了单点积分单元且没有沙漏控制改为全积分单元或增加沙漏控制系数检查泡沫材料阻尼材料一直不出现失效删除忘了加*MAT_ADD_EROSION或者失效阈值单位错误检查失效关键字确认单位制一致自动生成的混凝土参数结果异常混凝土标号填错或单位制不匹配确认无侧限抗压强度、单位制、模型版本JC模型在高应变率下应力离谱C值拟合范围外推或者应变率单位不一致检查参考应变率、限制外推范围5.4 大模型前的最后一道检查能量平衡视图经验丰富的老手跑大模型前都会养成一个习惯大模型算完第一步先看能量平衡。总能量曲线应该基本守恒动能、内能、沙漏能、滑移界面能各自有合理比例。如果沙漏能占比猛涨说明网格或单元算法有问题如果总能量曲线持续上升可能意味着接触或材料模型注入了额外能量。这道检查不花什么时间但对结果可信度的判断非常关键。材料参数相关的错误很多都会在能量平衡视图里暴露出来。比如单位制错误会导致密度异常进而导致惯性效应异常动能在初始阶段就明显不对状态方程参数错误则可能导致内能异常升高。能量平衡永远比云图颜色更诚实我建议每个项目都保留这一步。6. 关于材料参数管理我最后想分享的一些习惯6.1 用“版本管理”的思路维护材料库材料参数不是一次整理完就一劳永逸的它会随着你接触的新模型、新工况、新实验不断更新。我建议像管理代码一样管理你的材料参数库至少要做到每次参数更新都记录变更原因每条数据保留低版本存档确认已经过验证的参数打上“已验证”标记。我自己在本地维护了一组文件夹一份是“原始文献参数”一份是“可用K文件”一份是“验证算例”。每当我要把文献参数转成K文件时先把原始出处和单位制记入表格然后在验证算例文件夹里建一个对应的单单元测试跑通之后才移入“可用K文件”。这个过程听起来繁琐但能让你在未来半年后回头复查时快速定位“当时为什么这么填”。6.2 最后分享一个小技巧参数表里加一列“置信度”我习惯在材料参数表里加一列“置信度”分为高、中、低三档。高档参数来自实验标定且经过验证算例复核中档参数来自权威材料库且量级核对无误低档参数来自网络表格且未经验证。这个习惯在项目紧急、需要快速决策时特别有用看到低置信度参数你会在方案评审里主动说明风险看到高档参数你就可以更有底气地坚持仿真结论。仿真这个行当真正拉开差距的往往不是软件操作技巧而是对参数来源和适用边界的把控。材料参数这东西抄对了一张表只能让你跑通模型搞懂了参数背后的物理含义和标定过程你才算真正把仿真结果用明白。希望这篇关于LS-DYNA和AUTODYN材料参数与模型选择的内容能帮你把自己的材料库建得更扎实一点。本文还有配套的精品资源点击获取 LS-DYNA与AUTODYN材料参数选型、本构模型及单位换算实战指南
简介本资源是一份面向LS-DYNA与AUTODYN用户的专业级材料本构模型参数汇编适用于冲击动力学、侵彻仿真及弹体/破片毁伤效应建模等高阶工程仿真场景尤其适配装甲材料、穿甲弹体、战斗部破片等典型弹药结构的数值模拟需求。压缩包仅含1个HTML文件32KB内容为结构化整理的材料参数表涵盖多种常用金属、复合材料及高能材料的密度、屈服强度、硬化模型、失效准则等关键本构输入参数并附有AUTODYN与LS-DYNA双平台兼容性说明及典型MAT号索引便于工程师快速调用与校验。目前已有1733人学习下载内容源自工程实践验证可直接嵌入仿真前处理流程显著减少材料卡片搭建时间与参数试错成本是从事防护结构、武器效能仿真等方向研究人员与CAE工程师的实用速查工具。 做显式动力学仿真的朋友估计都经历过这种时刻翻遍整个材料库对着屏幕上的关键字卡片发呆不知道手头这个模型到底该选哪个材料号也不知道网上下载的那份“材料参数大全”能不能直接抄。尤其当你同时在用LS-DYNA和AUTODYN的时候两边材料库的组织方式、单位制、参数命名完全不一样稍不留神就得出一个看起来正常、实际全错的仿真结果。这篇内容围绕LS-DYNA和AUTODYN里的本构模型与材料参数展开讲讲我这些年攒下的材料模型选型思路、参数整理方法和踩坑记录希望对正在搞碰撞、冲击、爆炸、成形类仿真的朋友有点帮助。先说结论所谓“材料参数大全”真正值钱的不是那张参数表本身而是你拿到参数之后能不能判断“这组参数在什么条件下有效、适合什么模型、要怎么校验”。我会把常见材料模型按应用场景拆开讲再解释参数从哪来、怎么整理、怎么验证最后重点说说LS-DYNA和AUTODYN两套软件之间参数互通的坑。无论你是刚入门的新手还是已经跑过不少模型的老手这部分内容应该都能让你少走点弯路。1. 先搞清楚“材料模型”和“材料参数”到底是什么关系1.1 本构模型是数学框架材料参数是具体数据我经常看到有人把“材料模型”和“材料参数”混成一谈这其实是个容易出问题的误解。本构模型是描述材料力学行为的数学框架它决定了一组应力、应变、应变率、温度之间怎么互相影响而材料参数则是填进这个框架里的具体数字。比如Johnson-Cook模型有A、B、n、C、m这些参数A是参考应变率下的初始屈服应力B和n控制硬化行为C控制应变率敏感性m控制温度软化。这些参数离开了Johnson-Cook这个框架就没有意义同样换了框架之后材料参数必须重新标定不能直接平移。类比一下本构模型就像一辆车的底盘和悬挂结构材料参数就是具体调校数据。同一种钢材你可以用双线性弹塑性模型描述也可以用Johnson-Cook模型描述得到的仿真精度和适用范围完全不同。“哪个模型更高级就用哪个”这个思路在工程里不可取真正的选择依据是你关心的物理现象是什么、你要覆盖的应变率范围有多大、你手里有什么数据可以标定。1.2 显式求解器里的材料卡到底由哪几块组成在LS-DYNA里一张完整的材料卡通常分成几层密度和弹性参数如弹性模量、泊松比负责基础力学响应强度模型计算屈服应力和硬化状态方程EOS计算压力和体积变形关系失效准则决定单元什么时候删除或损伤开始。AUTODYN的经典材料模型则经常直接打包成“状态方程强度模型失效模型”三件套逻辑上更加清晰。这个分层关系直接决定了你整理参数时该怎么归档。我见过很多人只记了密度、弹性模量、屈服强度和几个硬化参数结果到了需要跑爆炸冲击、超高速碰撞这类高压工况时发现缺了状态方程参数只能临时找文献找到的EOS参数和强度模型参数还是不同来源、不同单位制的硬拼在一起结果自然没法看。所以一旦涉及高压、高应变率问题从源头就要按状态方程、强度、失效三块去收集参数。2. 常用本构模型按场景拆解选对模型才算入门2.1 金属材料从PLASTIC_KINEMATIC到Johnson-Cook金属材料是显式动力学里用得最多的一类。如果只做低速碰撞、成形分析且没有特别大的塑性变形用LS-DYNA里的*MAT_PLASTIC_KINEMATICMAT_003就足够了。这个模型本质是双线性弹塑性加Cowper-Symonds应变率修正关键参数有密度RO、弹性模量E、泊松比PR、初始屈服应力SIGY、切线模量ETAN和应变率参数C、P。这里有一个高频坑点ETAN填的是塑性切线模量也叫硬化模量不是弹性模量。假设某钢材弹性模量210GPa、屈服强度350MPa、抗拉强度600MPa、均匀延伸率10%估算切线模量大约在2.5GPa到3GPa左右很多人直接填了20GPa甚至200GPa结果模型刚度离谱偏大变形完全不对。正确做法是用真实应力应变曲线硬化段的斜率最好用实验数据拟合。如果问题涉及高应变率、高温或者大变形破坏Johnson-Cook模型MAT_JOHNSON_COOKMAT_015是业界标杆。它的表达式分成三项应变硬化项ABε^n、应变率项(1Clnε)、温度项(1-T*m)。注意A、B的单位和弹性模量一致n、C、m无量纲参考应变率通常取1/s或0.001/s要和标定条件对应。JC模型还带D1到D5五个断裂参数用来描述应力三轴度、应变率、温度共同作用下的失效应变。但在实际项目里Johnson-Cook模型也有它的毛病。它的应变率项是对数形式在10^4/s以下表现尚可到了10^5/s以上外推风险很大。另外它假设温度软化是单一幂次形式对某些铝合金、钛合金并不完全适用。这时候可以考虑LS-DYNA里的*MAT_MODIFIED_JOHNSON_COOKMAT_107或者分段表格化的JC变体直接把实验得到的应力应变曲线作为表格输入同时保留应变率和温度插值精度会明显提升。2.2 混凝土与岩土材料别指望一组参数吃遍所有标号混凝土材料的难点在于抗压/抗拉不对称、压力相关、应变率相关、损伤软化这些特性远不是弹塑性模型能描述的。工程上常见的有LS-DYNA的*MAT_CONCRETE_DAMAGE_REL3MAT_072R3、*MAT_CSCMMAT_159、*MAT_JOHNSON_HOLMQUIST_CONCRETEJHC和AUTODYN里常用的RHT模型。MAT_072R3是KC模型的一个简化版好处是只需要输入密度、泊松比、无侧限抗压强度和几个单位转换系数LS-DYNA可以自动生成剩余参数。它生成默认参数的过程依赖一个内部公式库适用于普通强度混凝土对高强混凝土、纤维混凝土或者含骨料特殊混凝土得谨慎。CSCM模型自带帽盖屈服面和损伤参数对隧道、地基、混凝土结构受撞问题比较合适但它的参数输入格式经常让新手摸不着头脑很多参数是可选的不填时用默认值问题是默认值未必适合你的单位制。如果你接触过AUTODYN会发现RHT模型在AUTODYN材料库里已经预置了很多岩石、混凝土的参数使用非常方便。但要留意这些预置参数不少来自于公开文献标定可能基于特定实验条件直接用于你的工况时必须做敏感性分析。混凝土这类材料每次仿真的核心问题是你的模型标定过哪个抗压强度等级C30和C50的参数差很远用C50的参数模拟C30构件应力峰值可能偏差30%以上。岩土材料则常用Mohr-Coulomb、Drucker-Prager以及LS-DYNA里的*MAT_FHWA_SOILMAT_147这类模型涉及黏聚力、内摩擦角、含水率、孔隙水压力等参数。这些参数一般来自土工实验没有实验数据时不要拍脑袋填否则连边坡稳定的定性趋势都可能算反。2.3 橡胶、泡沫与超弹性材料曲线比参数表更值钱橡胶、密封件、缓冲泡沫这类材料在汽车碰撞、包装跌落、防护结构中很常见。它们的核心特点是高度非线性弹性甚至在大变形下表现出体积近似不可压缩。LS-DYNA里的Mooney-Rivlin模型MAT_027和Ogden模型MAT_077适合橡胶类低密度泡沫用MAT_LOW_DENSITY_FOAMMAT_063可压碎泡沫用MAT_CRUSHABLE_FOAMMAT_026。我的建议是这类材料不要只抄几个弹性常数最好拿到完整的单轴压缩应力应变曲线以表格曲线方式输入。泡沫材料的迟滞、能量吸收、残余变形行为都体现在曲线形状里光靠几个材料号根本覆盖不了。低密度泡沫还要注意设置合适的阻尼系数DAMP否则卸载阶段可能出现剧烈的数值振荡这个系数本身没有实验意义就是数值稳定参数但调不好轻则沙漏能失控重则负体积报错。2.4 状态方程高压冲击问题里的另一半关键很多人做低速碰撞时从头到尾不碰状态方程但一旦涉及爆炸、高速冲击、水压冲击、航天器防护这类问题EOS就是决定成败的那一半。LS-DYNA常用的*EOS_GRUNEISEN适用于金属和水的冲击绝热关系*EOS_LINEAR_POLYNOMIAL适用于气体和一般流体*EOS_JWL用于炸药爆轰产物。Gruneisen EOS的核心参数包括C、S1、S2、S3和Gamma0。C是材料在零压力下的体声速S1等是冲击波速度与粒子速度关系式的斜率系数Gamma0是Gruneisen系数。以水为例C大约在1480m/s量级但在不同单位制下要换算成对应单位。很多人在低速流固耦合里也强行引入Gruneisen EOS其实没有必要线性多项式EOS反而更稳反过来爆炸冲击问题只用弹性模量而不填EOS那就完全没法计算冲击波传播。3. 参数从哪来如何积累一份真正能用的“参数大全”3.1 官方材料库和文档是第一站但别迷信LS-DYNA的LS-PrePost自带材料参数库打开界面可以看到按金属、非金属、混凝土、泡沫等分类的预置材料直接生成关键字卡片很方便。这些数据作为起点完全够用尤其是帮助你快速跑通模型。但官方库的局限也很明显它提供的多是典型值缺少标定条件和适用范围说明你要在仿真报告中交代参数来源时会发现没法溯源。AUTODYN的材料库在这方面做得更系统。AUTODYN把材料分成金属、非金属、炸药、水、空气、岩土等大类每个材料条目下挂载了EOS、强度模型、失效模型三组参数并且标注了数据来源文献。我在实际工作中经常用AUTODYN材料库作为参数交叉验证的基准尤其是同一种材料在LS-DYNA里参数怎么填不确定时去AUTODYN查一下对应量纲下的数值能很快发现量级错误。3.2 文献、数据库和实验标定怎么配合材料手册和公开数据库是第二站比如ASM手册、军工材料数据库、各大学公开的论文补充材料。网上流传的各种“LS-DYNA常用材料参数Excel表”不是不能用但使用前必须做三步核查第一步确认单位制第二步确认模型名第三步确认标定条件如温度、应变率、热处理状态。很多时候同一款铝合金挤压态和T6态的参数完全不同直接套用等于白算。如果项目对精度要求高最终归宿一定是实验标定。常规流程是先做准静态拉伸实验获得弹性模量、屈服强度、抗拉强度和硬化曲线再做不同应变率下的动态实验万能试验机做低速SHPB做中高应变率拟合Johnson-Cook的C值有条件的话做不同温度下的拉伸实验拟合m值最后用Hopkinson杆或Taylor撞击实验做模型验证。这个过程看着复杂但对高精度结构冲击仿真来说是绕不开的。我见过一个反面案例有人直接从文献里抄了一组Ti-6Al-4V的JC参数算出的弹道极限比实验低了30%。原因就是文献里的参数是基于退火态钛合金标定的而实验用的是锻造态屈服强度差了一大截。参数跟材料状态强相关这句话请务必写在你的参数管理表头里。3.3 把参数整理成带溯源信息的数据库真正的“参数大全”不是一页密密麻麻的数据而是一个可以持续积累、可溯源、可复用的数据结构。我的做法是用一张总表管理常用材料每个材料条目至少包含编号、材料牌号、材料状态、本构模型名、单位制、参数列表、数据来源、标定条件、适用范围、验证状态、备注这十一个字段。单位制和来源这两个字段绝对不能省因为这是后续排查问题最快的抓手。如果你长期在同一套单位制下工作还可以把整理好的材料卡保存成K文件片段用时直接include不用每次重新填卡。我自己维护了一个“通用材料卡”文件夹里面每个K文件都有命名规范比如“steel_q235_mat3_mm-ms-g.k”找起来一目了然。这个习惯第一次养成时花不了多少时间但是后面节省的返工时间数以十倍计。4. LS-DYNA和AUTODYN参数互通单位换算与格式差异4.1 单位制换算最容易翻车的环节我接触过不少同时使用LS-DYNA和AUTODYN的团队最典型的错误就是单位制混用。LS-DYNA本身不锁单位制只要保证封闭一致就行但LS-PrePost的交互式和某些转换工具默认按特定单位制解释容易让人产生“软件里填的就是标准值”的错觉。我通常把工作中最常用的一套单位制定义为mm-ms-g体系长度mm时间ms质量g应力MPa密度g/mm³速度mm/ms数值上等于m/s。在这个体系下钢材密度填7.85e-3弹性模量填210000屈服强度填350JC参数里A、B也按MPa填。如果从AUTODYN换到LS-DYNA或者反过来需要注意AUTODYN的经典参数常以米-千克-秒-帕m-kg-s-Pa为单位给出那么从AUTODYN里导出的模量2.1e11对应的是210GPa填到LS-DYNA的mm-ms-g体系必须改成210000直接填2.1e11会导致刚度大了100万倍模型完全锁死。反之从LS-DYNA的mm-ms-g体系换到m-kg-s-Pa则要乘以1e6。再补充一套可能遇到的cm-g-μs体系长度cm时间μs质量g密度g/cm³应力单位是Mbar1Mbar100GPa。这套单位制在爆炸冲击文献里很常见JWL状态方程参数经常用这套单位给出。因此当你在文献里看到JWL的A、B参数是“3.7e11”这类大数字时先确认单位是Pa还是Mbar再决定要不要除以1000。我把常用单位制的换算关系整理成表方便直接对照单位制长度时间质量密度单位应力/压力单位典型场景mm-ms-gmmmsgg/mm³MPa结构碰撞、汽车安全m-s-kgmskgkg/m³Pa通用流固耦合、AUTODYN默认cm-μs-gcmμsgg/cm³Mbar100GPa炸药爆轰、冲击波文献mm-ms-kgmmmskgkg/mm³GPa某些大型冲击模型这个表建议截图保存或者贴在工位旁边每次填材料卡之前看一眼成本几乎为零收益非常大。4.2 模型对应的映射关系同参数不同卡片写法AUTODYN的材料模型逻辑是“EOSStrengthFailure”三件套LS-DYNA则是通过独立的材料关键字、状态方程关键字、失效关键字组合实现。虽然底层物理模型可能是一样的但卡片写法完全不同。我总结过一套常用映射关系AUTODYN里的Linear EOS对应LS-DYNA的EOS_LINEAR_POLYNOMIALShock EOS对应EOS_GRUNEISENJWL EOS对应EOS_JWLAUTODYN里的Johnson-Cook强度模型对应LS-DYNA的MAT_JOHNSON_COOKAUTODYN里的Piecewise JC对应LS-DYNA的MAT_MODIFIED_JOHNSON_COOK或者带表格的JC变体AUTODYN里的Von Mises强度模型对应LS-DYNA的MAT_PLASTIC_KINEMATIC或者*MAT_ELASTIC_PLASTIC_HYDRO。失效模型这块差异更大。AUTODYN里失效条件通常是材料模型自带的参数比如最大拉伸压力、最大主应力、累计损伤等LS-DYNA则经常用MAT_ADD_EROSION给任何材料卡追加失效准则。很多人做LS-DYNA仿真时发现材料没有删除第一反应是材料卡填错了其实很可能是忘了加MAT_ADD_EROSION或者失效准则是基于压力而实际应力状态里压力没有达到阈值。4.3 互通时的三类典型坑第一类坑是压力符号约定。AUTODYN的很多材料模型里压力为正代表压缩而LS-DYNA的某些历史和应力输出里正负号惯例可能不一样尤其是在EOS卡片里填参数时要严格看手册里对压力的定义。我不止一次见过把Gruneisen EOS的C值填成负数导致模型一开始就爆炸式膨胀。第二类坑是失效模型参数的默认单位。*MAT_ADD_EROSION里压力单位必须和当前单位制匹配比如mm-ms-g体系下填的是MPa很多人从AUTODYN的Pa单位制切换过来后忘记转换等于失效阈值放大了一百万倍永远不失效。第三类坑是材料模型编号在不同LS-DYNA版本里可能有变体。同一个Johnson-Cook物理模型老版本用MAT_015新版本可能推荐MAT_015的增强版或者MAT_107。版本升级后老K文件里某些卡片会自动映射到新编号但参数含义未必完全一致建议升级后重跑一个标定算例做回归对比。5. 拿到参数后怎么验证单单元测试与常见坑位5.1 单单元测试五分钟找出参数低级错误无论参数是网上抄来的、文献里查来的还是实验标定的我的原则是不经过单单元测试验证的参数绝不上大模型。单单元测试思路特别简单建一个单单元模型施加单轴拉伸或压缩提取应力应变曲线和输入的理论曲线或者实验曲线对比。在LS-DYNA里单单元可以用一个8节点六面体单元一端固支另一端给恒定速度速度值选得足够小以避免应力波振荡干扰然后输出单元的应力应变。运行时间通常几秒钟K文件也就几十行。在AUTODYN里也可以建一个单单元网格做类似验证。这个测试能抓出大部分低级错误密度数量级不对、弹性模量差了1000倍、屈服强度填成了抗拉强度、曲线单位错位等。我发现一个特别高效的判断方式把仿真输出的名义应力-名义应变曲线跟输入的理论曲线叠在一起看。如果曲线斜率、屈服点、硬化趋势基本吻合说明参数基本可信如果曲线形态完全不对优先检查单位制和卡片填写。千万不要指望大模型能替你发现参数错误大模型里各种效应耦合在一起参数错了也常常能跑完只是结果没人敢用。5.2 动态参数验证SHPB数据的处理不能偷懒如果标定的是Johnson-Cook这类含应变率效应的模型单单元静态拉伸还不够最好用SHPB实验数据做动态验证。SHPB实验得到的是入射波、反射波、透射波三组波形需要经过二波法或三波法处理才能换算成材料的应力应变曲线。这里有一个新人常踩的坑直接把原始波形信号当成应力应变数据塞进参数拟合程序或者从论文里截取一张应力应变图再用软件抠数据点但忽略了原始实验的应变率并不是恒定值。正确做法是把SHPB数据按应变率分组每组取平均应变率下的代表性应力应变曲线再做JC参数拟合。拟合时还要注意应变范围SHPB在高应变率下试样容易出现绝热温升导致流动应力下降这个下降其实是温度效应不是材料本身的应变率软化。要想把C和m拟合得准最好结合不同温度和不同应变率的多批次数据做联合标定而不是单条曲线硬拟合。5.3 常见报错与异常结果速查我把这些年遇到过的高频问题整理成了一张速查表按症状、原因、排查方向排列症状最常见原因排查方向模型刚度异常大几乎不变形弹性模量单位制错误多填了6个数量级核对弹性模量数值是否是当前单位制下的正确值单元过早产生负体积材料失效删除不及时或者EOS缺失/参数不合理检查状态方程参数增大失效准则触发可能或细化网格沙漏能超过内能10%使用了单点积分单元且没有沙漏控制改为全积分单元或增加沙漏控制系数检查泡沫材料阻尼材料一直不出现失效删除忘了加*MAT_ADD_EROSION或者失效阈值单位错误检查失效关键字确认单位制一致自动生成的混凝土参数结果异常混凝土标号填错或单位制不匹配确认无侧限抗压强度、单位制、模型版本JC模型在高应变率下应力离谱C值拟合范围外推或者应变率单位不一致检查参考应变率、限制外推范围5.4 大模型前的最后一道检查能量平衡视图经验丰富的老手跑大模型前都会养成一个习惯大模型算完第一步先看能量平衡。总能量曲线应该基本守恒动能、内能、沙漏能、滑移界面能各自有合理比例。如果沙漏能占比猛涨说明网格或单元算法有问题如果总能量曲线持续上升可能意味着接触或材料模型注入了额外能量。这道检查不花什么时间但对结果可信度的判断非常关键。材料参数相关的错误很多都会在能量平衡视图里暴露出来。比如单位制错误会导致密度异常进而导致惯性效应异常动能在初始阶段就明显不对状态方程参数错误则可能导致内能异常升高。能量平衡永远比云图颜色更诚实我建议每个项目都保留这一步。6. 关于材料参数管理我最后想分享的一些习惯6.1 用“版本管理”的思路维护材料库材料参数不是一次整理完就一劳永逸的它会随着你接触的新模型、新工况、新实验不断更新。我建议像管理代码一样管理你的材料参数库至少要做到每次参数更新都记录变更原因每条数据保留低版本存档确认已经过验证的参数打上“已验证”标记。我自己在本地维护了一组文件夹一份是“原始文献参数”一份是“可用K文件”一份是“验证算例”。每当我要把文献参数转成K文件时先把原始出处和单位制记入表格然后在验证算例文件夹里建一个对应的单单元测试跑通之后才移入“可用K文件”。这个过程听起来繁琐但能让你在未来半年后回头复查时快速定位“当时为什么这么填”。6.2 最后分享一个小技巧参数表里加一列“置信度”我习惯在材料参数表里加一列“置信度”分为高、中、低三档。高档参数来自实验标定且经过验证算例复核中档参数来自权威材料库且量级核对无误低档参数来自网络表格且未经验证。这个习惯在项目紧急、需要快速决策时特别有用看到低置信度参数你会在方案评审里主动说明风险看到高档参数你就可以更有底气地坚持仿真结论。仿真这个行当真正拉开差距的往往不是软件操作技巧而是对参数来源和适用边界的把控。材料参数这东西抄对了一张表只能让你跑通模型搞懂了参数背后的物理含义和标定过程你才算真正把仿真结果用明白。希望这篇关于LS-DYNA和AUTODYN材料参数与模型选择的内容能帮你把自己的材料库建得更扎实一点。本文还有配套的精品资源点击获取