Abaqus非线弹性材料模型

最新推荐文章于 2024-09-01 10:43:04 发布

heartwind2016

最新推荐文章于 2024-09-01 10:43:04 发布

阅读量926

点赞数 26

文章标签：学习方法

本文链接：https://blog.csdn.net/heartwind2007/article/details/140907885

版权

Abaqus中超弹性材料的应力应变关系是非线性的，分析步中必须打开几何非线性。几何非线性会增加计算复杂度。为了不使用超弹性材料，在Abaqus中找到一个porous elasticity模型。这个模型可以给出随静水压力变化的弹性模量和泊松比。此处把abaqus帮助文档翻译并记录一下，方便以后查找。

翻译内容：

一、孔隙弹性模型：

（1）用于研究材料的弹性随压力变化而变化的行为；

（2）包含logarithmic模型和power law模型两种；

（3）可以具有依赖于温度和其他自定义场变量的属性。

二、对数孔隙弹性模型（logarithmic porous elasticity model）

对数孔隙弹性模型适用于小弹性应变（一般小于5%）。它是一个非线性、各向同性的弹性模型。该模型中，静水压力是体积应变的指数函数。该模型允许零或非零的弹性拉伸极限应力。

1. 体积行为的定义

通常，多孔介质材料体积行为的弹性部分可用与静水压力的对数成比例的体积弹性应变来准确描述：

$\frac{\kappa}{1+e_0}\rm{ln}\it(\frac{p_0+p_t^{el}}{p+p_t^{el}})=J^{el}-\rm1$

其中， $\kappa$ 是“对数体积模量（logarithmic bulk modulus）”； $e_0$ 是初始孔隙比； $p$ 是静水压力，定义式为

$p=-\frac 13\rm trace({\bold{\sigma}})=-\frac 13(\sigma_{11}+\sigma_{22}+\sigma_{33});$

$p_0$ 是初始静水压力； $J^{el}$ 是当前构型与初始构型体积比的弹性部分； $p_t^{el}$ 是材料的“弹性拉伸强度（elastic tensile strenght, or tensile limit）”（当 $p\rightarrow-p_t^{el}$ 时 $J^{el}\rightarrow\infty$ ）。

图1 多孔弹性体积行为

在Input文件中使用：

使用下面的三行选项来定义孔隙弹性材料：

*POROUS ELASTIC, TYPE=LOGARITHMIC,SHEAR=G or POISSON to define κ and pelt

*INITIAL CONDITIONS, TYPE=STRESS to define p0

*INITIAL CONDITIONS, TYPE=RATIO to define e0

在Abaqus/CAE中使用：

Property module: material editor: Mechanical>Elasticity>Porous Elastic: Type: Logarithmic
Load module: Create Predefined Field: Step: Initial: choose Mechanical for the Category and Stress for the Types for Selected Step
Load module: Create Predefined Field: Step: Initial: choose Other for the Category and Void ratio for the Types for Selected Step

2. 剪切行为的定义

孔隙材料的偏弹性行为可以用两种方式中的任何一种来定义。

(1) 通过定义剪切弹性模量来描述

给定剪切弹性模量， $G$ ，偏应力， $\bold S$ ，与总弹性应变的偏应变部分 $\bold{e}\it^{el}$ 相关，

$\bold{S}=2G\bold{e}^{el}$

在这种情况下，剪切行为不受材料压实的影响。

在Input文件中使用：

使用下面的三行选项来定义孔隙弹性材料：

*POROUS ELASTIC, TYPE=LOGARITHMIC,SHEAR=G

在Abaqus/CAE中使用：

Property module: material editor: Mechanical>Elasticity>Porous Elastic: Type: Logarithmic: Shear: G

(2) 通过定义泊松比来描述

定义泊松比 $\nu$ ,瞬时剪切模量可由瞬时体积模量和泊松比来定义

$G=\frac{2(1-2\nu)(1+e_0)}{2(1+\nu)\kappa}(p+p_t^{el})\rm{exp}(\it\varepsilon_{vol}^{el})$

其中， $\varepsilon_{vol}^{el}=\rm{ln}\it{J^{el}}$ 是弹性体积改变的的对数测量。在这种情况下

$d\bold{S}=2Gd\bold{e}^{el}$

因此，弹性剪切刚度随材料的压实而增大。积分该公式可给出总应力和总弹性应变的关系。

Input File Usage

*POROUS ELASTIC, TYPE=LOGARITHMIC, SHEAR=POISSON

Abaqus/CAE Usage

Property module: material editor: MechanicalElasticityPorous Elastic: Type: Logarithmic: Shear: Poisson

二、幂律孔隙弹性模型

在幂律孔隙弹性模型中，杨氏模量与静水压力的幂相关：

$E=E_{ref}[\frac{p+p_0}{p_{ref}+p_0}]^n,for(p>0);E=fE_{ref},for(p\le0)$

其中， $E_{ref}$ 是静水压力为 $p_{ref}$ 是的参考杨氏模量， $p_0$ 和 $n$ 是材料常数。 $f$ 用来定义等效拉伸压力状态下的弹性模量（是个常数？），由下式计算：

$f=[\frac{p_0}{p_{ref}+p_0}]^n.$

泊松比是一个静水压力的指数函数

$\nu=\nu_0+(\nu_\infty-\nu_0)(1-e^{-mp}),for(p>0);\nu=\nu_0,for(p\le0)$

其中， $\nu_\infty$ 为p趋于无穷时的泊松比， $\nu_0$ 时p=0时的泊松比，m是个材料常数。

下面给一个幂律孔隙弹性模型参数的例子。

$E_{ref}=10\rm MPa$

$p_{ref}=2\rm MPa$

$p_0=4\rm MPa$

$n=1.2$

由以上常数，计算出的 $f=0.6147$ 。

$\nu_\infty=0.49$ （不是孔隙材料，只用abaqus的非线弹性属性）

$\nu_0=0.47$

$m=0.5$

由以上常数计算出的弹性模量随静水压力的变化用下表简单表示。

p(MPa)	E(MPa)	nu
0	6.147	0.47
10	27.642	0.4899
20	52.780	0.4900
30	80.167	0.4900

heartwind2016

关注

26
点赞
踩
22

收藏

觉得还不错? 一键收藏
0
评论
Abaqus非线弹性材料模型

通常，多孔介质材料体积行为的弹性部分可用与静水压力的对数成比例的体积弹性应变来准确描述：其中，是“对数体积模量（logarithmic bulk modulus）”；是初始孔隙比；是静水压力，定义式为是初始静水压力；是当前构型与初始构型体积比的弹性部分；是材料的“弹性拉伸强度（elastic tensile strenght, or tensile limit）”（当时图1 多孔弹性体积行为使用下面的三行选项来定义孔隙弹性材料：MechanicalElasticityType。
复制链接

扫一扫