【OptiStruct要领】材料非线性模型

首先介绍下OptiStruct中弹塑性材料模型。

弹塑性材料在加载的初始阶段即弹性阶段，应力—应变关系基本呈线性，斜率即弹性模量。但当应力大于屈服应力时，材料进入塑形，若此后继续加载，应力-应变不呈线性关系，一些情况下仍可近似为线性，但斜率与弹性阶段不同。卸载时，卸载曲线和弹性段曲线斜率相同，但完全卸载后，材料将保留永久的塑形变形。这里由于屈服点和比例极限相差较小，故而OptiStruct中设定两者是同一点。OptiStruct中弹塑性材料的应力应变曲线为：

屈服准则用于预测材料在一定应力下是处于弹性或者塑性状态，即判别材料是否有塑性应变产生。OptiStruct中用到的屈服准则为Mises屈服准则。Mises屈服准则认为，当材料在复杂应力状态下的形状改变能达到单项拉伸屈服时的形状改变能时，材料开始屈服。Mises屈服准则是除了土壤和脆性材料外典型的屈服准则。

‒ 应力状态在屈服面之内为弹性状态

‒ 应力状态在屈服面上为塑性状态

‒ 应力不会落在屈服面之外

前面提到一般工程材料当达到一定的应力状态后，即出现塑性流动，呈现材料非线性。流动准则是指当材料发生屈服时，塑性应变增量的方向。OptiStruct中流动准则假定塑性势函数与屈服函数一致，塑性变形增量方向总是沿塑性势法线方向。

当材料发生屈服卸载后，重新加载，材料将重新屈服，但材料屈服极限提高了，这种现象成为冷作硬化，此时的屈服成为后继屈服。当材料发生硬化后，屈服准则也将发生改变。

OptiStruct中考虑的硬化模型为，各向同性硬化、随动硬化以及混合硬化。在MATS1中控制参数为HR。

在进行弹塑性分析时，我们选择进行增量/全量的有限元分析方法。全量分析应变包括弹性和塑性应变，而增量分析中应变仅为塑性应变部分。在OptiStruct中，我们可以通过参数MATS1—TYPSTRN来进行控制。那么对应不同的分析方法，应力应变曲线表格的输入MATS1-TID也是不一样的。

➡全量分析(TYPSTRN=0)应力应变曲线TABLES1输入：

起点：(X1 = 0, Y1 = 0)；第二点(X2, Y2)=初始屈服点，且斜率=弹性模量（E）

➡增量分析(TYPSTRN=1)应力应变曲线TABLES1输入：

起点：(X1=0, Y1 )

（具体卡片关键字详见HWSOLVER HELP—MATS1）

OptiStruct中的弹塑性分析可以在积分点、高斯点、角点输出塑性应变、应力。

等效塑性应变

只要有应变输出要求，都会输出等效塑性应变

塑性应变分量

               STRAIN（PLASTIC）=YES

节点塑性应变的输出

               GPSTRAIN(PLASTIC) = YES

等效塑性应力

               STRESS(VON)=YES

但是需要注意一些限制

➡只能输出实体单元的节点塑性应变

➡只支持静力分析

➡只能输出到H3D文件

这里需要注意的是，如果加载过程是承接上一步加载，而不是若干个独立的工况，则需要在GLOBAL_CASE_CONTROL这里设置CNTNLSUB,YES

这是一个材料拉伸实验的例子。如下图所示试件，一端固支，一端收+x方向拉力。考察其塑性应变。

其中应力应变曲线由TABLES1表格输入。起点为屈服点。

位移：

应力：

以及spc force的实验及仿真曲线的对比。蓝色为实验曲线，红色为仿真结果。

除了弹塑性材料本构外，OptiStruct还可以定义超弹性材料。

超弹性材料实质就是一种非线性弹性材料，在较大的变形下也保持弹性特性。超弹性材料的行为比金属材料行为要复杂得多，该模型可用于分析各种橡胶类、泡沫类材料，这种材料能够承受大应变和大位移，但体积改变极小（超弹性材料一般为不可压或近似不可压缩材料），比如吹气球，又如橡胶按钮受压缩等，这种分析需用到大应变理论。

超弹性材料应力应变关系通过应变能密度函数来描述

S = second Piola_Kirchhoff stress tensor

E = Lagrangian strain tensor

U= strain energy function per unit undeformed volume

C = right Cauchy-Green deformation tensor

在实际应用中应变能密度一般采用体积应变及主方向上的伸长率来描述

让我们先看下超弹性中的几个术语吧~

➡伸长率：变形后长度除以变形前长度

➡体积应变：变形后体积除以变形前体积；对于不可压缩的超弹性材料=1

➡应变能密度函数：应变能密度函数一般通过3个不变量来描述，其中3个不变量又可用3个主方向的伸长率来表示

对于不可压缩超弹性材料，其应变能密度可通过偏应变张量不变量及体积应变来表示。

为了完整描述超弹性材料模型，理论上需要6种纯应变状态的力学实验：单轴拉伸、单轴压缩、等双轴压缩、等双轴拉伸、平面拉伸、平面压缩。

由于超弹性材料的不可压缩特性，有三组实验理论上存在应力状态等效关系。

因此，超弹性材料的3中基本力学实验为：单轴拉伸（等效的等双轴压缩）、单轴压缩（或等效的等双轴拉伸）、剪切（平面拉伸或压缩）。

在OptiStruct中，超弹性材料模型有：Mooney-Rivlin多项式形式(MOONEY)，简化的多项式形式(RPOLY)，Physical Mooney-Rivlin 模型(MOOR), Neo-Hookean 模型(NEOH)，Yeoh 模型(YEOH)，Arruda-Boyce 模型(Model=ABOYCE)，Ogden Material 模型(Model=OGDEN)

方法1

直接输入简单载荷下的实验曲线，TAB1单轴拉伸，TAB2等双轴拉伸，TAB4纯剪，TABD体积应变下的压力曲线，OptiStruct会自动拟合多项式参数进行求解。

方法2

直接输入多项式的各个参数

介绍完理论部分，按照国际惯例，接下来就是动手环节啦~

这个模型我们将进行橡胶材料加载卸载的变形分析，模型对象如下图所示