基于JC模型的umat子程序实例-Hopkinson压杆的隐式求解器数值仿真实现

众所周知，Abaqus中自带的JC模型只能应用于显示求解器（Abaqus/Explicit），但是隐式求解器(Abaqus/Standard)通常具有更高的精度，而且随着研究的进步JC模型也出现了新的修正形式，因此我们有必要基于Umat实现传统的JC模型，以便有需要的朋友在此基础上进一步创新，避免重复“造车”！本文基于JC模型的Umat子程序实现hopkinson压杆的数值模拟，本文的一些参数设置和程序借鉴了卢剑锋老师的硕士论文《冲击载荷作用下材料和结构力学行为有限元模拟》，且没有考虑温度软化。

装配图如下

JC模型的方程形式如下

试样的材料参数如下

入力杆和出力杆使用线弹性材料，弹性模量和泊松比分别为200GPa和0.3。

在入力杆的一端施加应力脉冲

入力杆某点和出力杆某点上的应力波时程曲线如下

建议大家按照前面提到的论文自己独立把程序走一遍，也欢迎大家下载本次的*.cae文件和*.for文件，后期打算基于JC模型的umat子程序做一个侵彻的案例，感谢持续关注！