ANSYS粘弹体分析分析和总结.docx

ANSYS中粘弹材质属性参数输入和分析 1

ANSYS中表征粘弹性属性问题 1

Prony级数形式 1

Maxwell形式 3

建模与载荷条件 5

模型设计 5

有限元建模 5

理论解析解计算式 6

有限元数值解与结果比较 6

Plane183，Prony级数方式 6

1.4.5 算例结论 10

ANSYS中粘弹材质属性参数输入和分析

ANSYS中表征粘弹性属性问题

粘弹性材料的应力响应包括弹性部分和粘性部分，在载荷作用下弹性部分是即时响应的，而粘性部分需要经过一段时间才能表现出来。一般的，应力函数是由积分形式给出的，在小应变理论下，各向同性的粘弹性本构方程可以写成如下形式：

???t2G?t???ded??I?tK?t???d?d? (0.1)

0 d? 0 d?

其中

?＝Cauchy应力

G?t?＝为剪切松弛核函数K?t?＝为体积松弛核函数

e＝为应变偏量部分（剪切变形）

?＝为应变体积部分（体积变形）

t＝当前时间

?＝过去时间

I＝为单位张量。

该式是根据松弛条件本构方程(0.1)，通过将一点的应变分解为应变球张量（体积变形）和应变斜张量（剪切变形）两部分，推导而得的。这里不再敖述，可参考相关文献等。

ANSYS中描述粘弹性积分核函数G?t?和K?t?参数表示方式主要有两种，一种是广义 Maxwell单元

（VISCO88和VISCO89）所采用的Maxwell形式，一种是结构单元（如Plane183，Plane182等）所采用的Prony级数形式。实际上，这两种表示方式是一致的，只是具体数学表达式有一点点不同。

Prony级数形式

用Prony级数表示粘弹性属性的基本形式为：

G?t??G

?nGG

? t?

exp?? ?

(0.2)

??? i ?Gi?1 i

?

?

K?t??K

?nKK

? t?

exp?? ?

(0.3)

??? i ?Ki?1 i

?

?

其中，G 和G是剪切模量，K

和K是体积模量，?G和?K

是各Prony级数分量的松弛时间。再定义下面

??i

相对模量

?G?G G

? i i i

(0.4)

i i 0

?K?K K

i i 0

(0.5)

其中，G，K

0 0

分别为粘弹性材质（固体推进剂）的瞬态模量，并定义式如下：

G ?G?t?0??G

0 ?

K ?K?t?0??K

?n G

Gi

G

i?1

?nKK

(0.6)

(0.7)

0 ??i

i?1

在ANSYS中，Prony级数的阶数n 和n 可以不必相同，当然其中的松弛时间?G和?K

也不必相同。

G K i i

对于粘弹性问题，粘弹体的泊松比一般是取为时间的函数????t?。不过有时情况允许也可近似设为常

数，这时根据弹性常数关系就有：

G?t??

E?t?

2?1???

E?t?

(0.8)

K?t??3?1?2??

其中，E?t?为松弛模量，由实验来确定。E?t?,G?t?,K?t?的相应系数比相同。

这样就可以将G?t?和K?t?统一于E?t?形式。若我们将松弛模量表示为Prony级数形式，即：

E?t??E

??n

? t?

Eexp?? ?

(0.9)

??i

i?1

? ?i?

于是，G?t?和K?t?中有，n?n

?n ，?

??G

??K，?

??G??K。类似于G、K

，我们也同样定

义瞬态松弛模量E：

0

G K i i i

i i i 0 0

E ?E?t?0??E

0 ?

?n E

Gi

G

i?1

(0.10)

这样，由(0.8)可得

E

G ? ? 0 ?

0 21??

E

(0.11)

K ? ? 0 ?

0 31?2?

要注意的是，ANSYS中对Prony级数的支持项数不能超过6项，即n?6。这确实是一个遗憾。

另外，

――――――――――――――――――――――――――――――――――――――――――

TheviscoelasticityinputforSHELL181,PLANE182,PLANE183,SOLID185,SOLID186,SOLID187,SOLSH190,SHELL208,andSHELL209consists