Abaqus UMAT example

简介：
Abaqus UMAT 案例解析：Johnson-Cook 本构模型 改写说明Johnson-Cook本构模型作为一种在高速冲击及塑性成型等领域具有广泛应用的本构关系式，其主要功能在于精确表征材料在高温和高速状态下流动应力场的特性。该模型的具体表达式为：$$\sigma = A + B \varepsilon^n [1 + C\ln(\dot{\varepsilon}^*)] [1 - (T^*/T_m)^m]$$其中各项参数包括：等效应力σ、应变ε、无量纲化的速率值$\dot{\varepsilon}^*$、无量纲温度$T^*$以及材料常数A、B、n、C和m。具体而言，σ表示等效应力状态，ε代表总的应变程度，$\dot{\varepsilon}^*$为无量纲化的速率值，$T^*$则表示相对于临界温度的比例关系。而参数A和B分别对应材料的基体强度和强化度，n是应变 hardening指数，C用于描述温度敏感性的影响范围，m则是硬化随动量变化的调节因子。三、UMAT子程序解析通过分析给定的代码片段，我们能够深入理解其内部工作原理：UMAT的主要输入参数为： - 应力状态信息：通过`STRESS`获取。 - 材料状态信息：由`STATEV`提供。 - 应力对弹性应变的影响程度：体现于`DDSDDE`中。 - 应变量密度值：表示为`sse`。 - 塑性功密度值：记录在`spsd`里。 - 热能密度信息：包含在`s cd`内。 - 恢复力大小：通过`rpl`获得。 - 应力随时间变化率：反映在`ddsddt`中。 - 恢复力对弹性应变的变化程度：由`drplde`描述。 - 恢复力对时间的响应速度：体现在`drplt dt`. - 总变形信息：通过`stran`呈现。 - 增量变形数据：记录在dstrain中. - 当前时间步的时间值：指定为time. - 时间步增量大小：设定在dtime里. - 环境温度状况：由temp给出. - 温度变化趋势：通过dtemp反映。 - 预定义变量设置：使用predef实现. - 定义变量的改变量：记录为dpre. - 材料名称标识：以MATERL`指定. - 直接应力分量独立数量：设定在Ndi`里. - 独立剪切应力分量数目：由NsHR确定。 - 总体应力分量总数：通过Ntens给出。 - 材料状态变量个数：由nstatv决定。 - 材料属性数据集合：存储于props中。 - 属性数量设定：以Nprops`完成. - 单元坐标信息：记录在coords`里. - 旋转矩阵数值：通过Drot`获得。 - 预测新时间步值：由Pnewdt`提供。 - 单元类型分类：指定为Celement`. - 变形梯度初始状态：由dfrgcd0给出. - 变形梯度增量数据：记录在dfrgcd1里. - 元素编号标识：通过Noel确定. - 积分点编号信息：由Npt指明。 - 层号数值设定：以Kslay完成。 - 层内积分点数目：设定为Kspt值。 - 时间步编号指定：使用kstep确定. - 增量编号赋值：通过kinc实现.材料属性及模型参数 Johnson-Cook 模型的相关参数在 PROPS 数组中定义： - 表示为杨氏模量的值由 PROPS(1) 给出。 - 泊松比设置在 PROPS(2) 位置。 - 弹性体应变能转换效率存储于 PROPS(3) 中。 - 模型中的参数 A 值定义在 PROPS(4) 处。 - B 参数对应于 PROPS(5) 的值。 - 非线性硬化指数设置为 PROPS(6)。 - C 参数则由 PROPS(7) 给出。 - 最后，幂指数 n 存储在 PROPS(8) 中。通过杨氏模量和泊松比计算出弹性模量，其中剪切模量(G)和体积模量(K)，从而求得弹性刚度矩阵。通过弹性刚度矩阵计算弹性应变所导致的应力增量变化，并进一步确定当前的应力状态。首先，通过 Johnson-Cook 模型求解塑性应变的增量；接着，在塑性流动准则的基础上进行等效应力和塑性应变的求解或确定。 四、代码实现细节 在给定的代码片段中，可观察到具体实施流程如下： 1. 当 `NDI` 值不等于3时，应显示错误提示信息，因该UMAT仅在单元拥有三个直接应力分量的情况下有效。 2. 基于提供的杨氏模量和泊松比参数，计算得出剪切模量(G)及体积模量(K)，随后构建完整的弹性刚度矩阵。 3. 通过结合弹性刚度矩阵与增量应变张量，可精确计算出新的应力状态。 4. 可计算出当前的塑性应变增量值，并相应地更新整体的塑性应变状态以反映最新变化。以Johnson-Cook模型为例，在ABAQUS软件中UMAT的应用过程具有重要的指导意义。该方法涵盖从材料参数设定、应力-应变关系建立直至塑性行为仿真四个关键环节，为材料力学研究者与工程实践者提供了完整的建模思路和操作指导。