非线性分析：粘弹性材料分析_（4）.本构关系及非线性特性.docxVIP

PAGE1

PAGE1

本构关系及非线性特性

1.本构关系的基本概念

1.1本构关系的定义

本构关系是指材料在受到外力作用时，其应力与应变之间的关系。在固体力学中，本构关系是描述材料行为的重要工具，通过它我们可以预测材料在不同加载条件下的响应。对于线性弹性材料，本构关系通常由胡克定律描述，即应力与应变成正比。然而，对于粘弹性材料，本构关系则更加复杂，因为它不仅依赖于当前的应变状态，还与时间有关。

1.2粘弹性材料的本构模型

粘弹性材料的本构模型可以分为两类：线性粘弹性模型和非线性粘弹性模型。

1.2.1线性粘弹性模型

线性粘弹性模型假设材料的应力与应变之间存在线性关系，并且这种关系随时间变化。常见的线性粘弹性模型包括：

Maxwell模型：由一个弹性元件和一个粘性元件串联组成。Maxwell模型可以描述应变随时间的松弛行为。

Kelvin-Voigt模型：由一个弹性元件和一个粘性元件并联组成。Kelvin-Voigt模型可以描述应力随时间的蠕变行为。

标准线性固体模型：由一个Maxwell模型和一个弹性元件并联组成。这种模型可以同时描述应变松弛和应力蠕变。

1.2.2非线性粘弹性模型

非线性粘弹性模型考虑了材料在大变形或高应变率下的非线性行为。常见的非线性粘弹性模型包括：

Burgers模型：由两个Maxwell模型并联组成，可以描述更复杂的松弛和蠕变行为。

Mullins模型：用于描述橡胶材料在大应变下的非线性行为。

Prony级数模型：通过一系列指数函数来拟合材料的松弛行为。

1.3本构关系的数学表示

本构关系可以用数学方程来表示。对于线性粘弹性材料，常用的数学表示方法有：

松弛模量：描述材料在恒定应变条件下的应力随时间的变化。

蠕变柔量：描述材料在恒定应力条件下的应变随时间的变化。

对于非线性粘弹性材料，本构关系通常更加复杂，需要考虑应变率、温度等因素的影响。

2.粘弹性材料的本构方程

2.1线性粘弹性模型的本构方程

2.1.1Maxwell模型

Maxwell模型的本构方程可以表示为：

σ

其中：-σt是应力-εt是应变-E是弹性模量-

2.1.2Kelvin-Voigt模型

Kelvin-Voigt模型的本构方程可以表示为：

σ

其中：-σt是应力-εt是应变-E是弹性模量-

2.1.3标准线性固体模型

标准线性固体模型的本构方程可以表示为：

σ

其中：-E1和E2是弹性模量-

2.2非线性粘弹性模型的本构方程

2.2.1Burgers模型

Burgers模型的本构方程可以表示为：

σ

其中：-E1和E2是弹性模量-Gt是松弛模量-

2.2.2Mullins模型

Mullins模型的本构方程可以表示为：

σ

其中：-Eε和η

2.2.3Prony级数模型

Prony级数模型的本构方程可以表示为：

G

其中：-G∞是无穷时间的松弛模量-Gi是第i个松弛模量-τi是第

2.3本构方程的数值解法

2.3.1显式时间积分法

显式时间积分法是一种常用的数值解法，适用于求解瞬态问题。以标准线性固体模型为例，其显式时间积分法可以表示为：

σ

2.3.2隐式时间积分法

隐式时间积分法适用于求解稳态问题。以Burgers模型为例，其隐式时间积分法可以表示为：

σ

2.4本构方程的软件实现

2.4.1Python实现

以下是一个使用Python实现标准线性固体模型的显式时间积分法的例子：

importnumpyasnp

importmatplotlib.pyplotasplt

#定义材料参数

E1=1000#弹性模量

E2=500#弹性模量

eta=100#粘度

#定义时间步长和总时间

dt=0.1

t_max=100

#初始化应力和应变数组

t=np.arange(0,t_max,dt)

stress=np.zeros_like(t)

strain=np.zeros_like(t)

#应用恒定应变

strain+=0.01

#计算应力

foriinrange(1,len(t)):

stress[i]=stress[i-1]+E1*(strain[i]-strain[i-1])+E2*(strain[i]-(1/E2)*np.trapz(np.exp(-(t[i]-t[:i])/eta),t[:i]))

#绘制应力-时间曲线

plt.plot(t,stress,label=Stress)

plt.xlabel(Time(s))

plt.ylabel(