橡胶有限元分析之研究_魏泳涛.pdfVIP

第 1 卷 第5 期 四川联 合大学学报( 工程科学版) Vol. 1 No. 5 1997 年9 月 JOU RNAL OF SICHUAN UN ION U NIVERSITY (ENGINEERING SCIENCE EDIT ION) Sept. 1997 魏泳涛 于建华 ( 四川联合大学土木工程及应用力学系 成都61 65) 本文研究了橡胶材料 限元应力分析的两种途径; 罚 限单元法及混合插 值 限单元法。分别将罚函数及静水压力引入材料应变能密度函数, 对橡胶类不可 压缩超弹性材料进行了 限元分析; 数值结果和解析解符合得较好, 表明了算法的 效性。 不可压缩超弹性材料; 罚函数; 静水压力; 限单元法 O242 1 - 橡胶( 天然橡胶或硫化橡胶) 从本质上说是一种粘弹性材料, 其力学持性与时间及温 度 关, 即呈现出蠕变、松弛、老化等现象。然而, 当时间较短且温度变化不大时, 则可认 为橡胶材料是完全弹性的: 即当载荷被撤走后, 橡胶材料完全恢复到它未变形的状态。大 量实验阐明了橡胶的弹性源于其热力学过程中的熵, 因此其弹性较之其它弹性体是完全 不同的。在本文中, 我们根据连续介质力学唯象的观点来讨论橡胶的力学特性。这种方 法根据大量实验事实对橡胶的力学特性作如下假定: 1) 橡胶在变形过程中存在着自由 能函数, 即等温条件下的应变能密度函数W, 它可以表达为变形状态的函数, 满足此条件 的材料称为超弹性材料; 而此函数的形式及其中所包含的常数则应实验确定; 2) 橡胶在 变形中体积变化极小, 因而可认为是不可压缩的; 3) 橡胶为各向同性材料。 文献[ 1] 列出了常用的应变能密度函数的形式, 其中Mooney - Rivilin 理论较好地描 述了橡胶类不可压缩超弹性材料在大变形下的力学特性, 且给出如下形式的应变能密度 函数: W = C ( I - 3) + C ( I - 3) 且 I = 1 ( 1) 1 1 2 2 3 其中C 、C 为材料常数, I 、I 和I 分别是Cauchy - Green 变形张量的三个不变量。 1 2 1 2 3 不可压缩超弹材料用物质描述的应力- 应变关系为[ 2] 收稿日期 1996 - 9 - 17 魏泳涛 等 橡胶 限元分析之研究 79 X X W i j S = 2 - P (2) ij C x x ij l l 其中, Sij 是第二类Piola- Kirchhoff 应力张量, P 是一坐标函数, Cij 为右Cauchy - Green 张量, X 及x 分别为Langrange 和Euler 坐标。将式( 1) 代入式(2) , : i i I I X X 1