橡胶类材料大变形本构关系和其有限元方法ij.pdfVIP

Vol. 20 　No. 4 　　 　　　　　　 　 　　　　　　　　　 　 第20卷 第4期 固 体 力 学 学 报 1999年 12 月 ACTA MECHANICA SOLIDA SINICA December 1999 橡胶类材料大变形本构关系及其有限元方法 危银涛 　　杨挺青 杜星文 (华中理工大学力学系 ,武汉 ,430074) ( 哈尔滨工业大学复合材料研究室 ,哈尔滨 ,150001) 摘　要 　讨论大变形拟不可压缩橡胶类材料的本构关系及有限元分析方法. 采用乘法分 解 ,将变形梯度表示成等容和体积变形两部分 ,在此基础上 ,推导了克希荷夫应力和格林应变表 示的 Yeoh 形式应变能橡胶类材料的本构关系及数值处理方法. 为处理不可压缩问题 ,采用三场 变分原理 ,其中静水压力 ,体积膨胀 , 以及位移均作为独立变量进行处理. 并指出该变分原理同 胡鹫津广义变分原理的联系. 变形采用相容等参插值 ,压力及体积膨胀采用低阶插值 ,推导了 详细的有限元列式. 最后给出了两个数值算例 ,结果表明了该方法的有效与可靠. 为这类材料的 精确的有限元分析打下了良好的基础. 关键词 　橡胶 ,Yeoh 模型 ,不可压缩性 ,广义变分原理 ,有限元法 1 　引言 橡胶类材料广泛应用于轮胎 、胶管、缓冲气囊、阻尼器等工程结构 ,这种材料有其独特的 力学性能 : 超弹性 ,不可压缩性 ,大变形等. 其本构关系的表征、确定及数值实现 ,和考虑不 可压缩性的大变形非线性弹性的有效数值方法都较为困难 ,本文主要讨论这两个问题. 在橡胶结构分析与设计中建立和发展适当的本构关系和分析方法是极其重要的。本文 首先给出了以第二类克希荷夫应力和格林应变表示的橡胶类材料本构关系 ,其中用到了变 形梯度的乘法分解。结合 Yeoh 模型的具体形式 ,给出了本构关系的具体数值实现。Yeoh 模型形式简单而又有足够的精度 ,而且仅由单轴拉伸实验即可确定其系数。因而近年来 Yeoh 模型较为流行[1～3 ] . 乘法分解可以将变形梯度分解为等容变形部分和体积变形部分 , 这对处理不可压缩问题非常有利。本文的一个主要工作是将 Yeoh 模型与变形梯度的乘法 分解结合 ,给出橡胶类材料的本构关系及详细的数值实现方式 ,为进一步的有限元分析打下 了良好的基础。本文的另一个主要工作是考虑不可压缩性的橡胶结构非线性有限元分析 , 其中考虑了橡胶的大变形非线性弹性。采用三类变量变分原理 ,将静水压力 ,体积变形和位 移场进行独立的处理 ,避免了体积闭锁现象。变形采用相容等参插值 ,压力及体积变形采用 低阶插值。并建立了胡鹫津广义变分原理与三类变量变分原理的联系。本文的推导均在 Lagrangian 描述下进行 ,本构关系 ,变分原理及 FEM 实现协调一致 ,避免了某些文献中的混 乱. 最后给出了两个算例分析 ,结果表明了本文方法的有效性和可靠性. ( ) 国家自然科学基金项目 资助. 收到第 1 稿 收到修改稿. © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 　　　　　　　　　　　　　　　　　固体力学学报 　　　　　　　　　　　　　　　1999 年 第 20 卷 2·82 · 2 　基本理论 变形梯度一般写为 x xi