短纤维／橡胶复合材料数值模拟.pdfVIP

现代橡胶技术 2011年第37卷 短纤维／橡胶复合材料 的数值模拟 吴贻珍 编译 摘 要：应变能密度函数的超弹性体本构模型有各种形式，既有以右 Cauchy—Green张量不变量，又有以主伸长为基础， 均假设材料为各向同性。对于像传动V带由纤维增强构成的超弹性材料 ，大多具有横向各向同性特征，则应变能量函数有 5 个基本不变量，除3个通用基本不变量外，还附加第 4和第 5不变量。本文提出了一个纤维增强材料超弹性体多项式应变能 函数本构模型并确定了其材料常数。 关键词：横向各向同性；超弹性体；纤维增强；各向异性；本构方程；数值模拟 1 前言 的宽度方向(如图 1的箭头方 向)。定向后 的短纤 橡胶材料具有非线性弹性性能，在有限元分 维／橡胶复合材料 (SFRP)表现出很强的各 向异 析(FEM)中一般作为超弹性体处理。按照 Fung 性 ，即在短纤维取 向方 向的拉伸刚性非常大，而且 有关超弹性体的定义，所谓超弹性材料，如果它具 从拉伸到压缩的切 向刚度还表现出明显的过渡 有每单位质量的应变能函数w，则这是一个与在 区。 自然均质无应力状态下有关的应变张量解析函 数。这样每单位质量应变能变化率等于应力所做 的功的速率。换言之，超弹性体是一种含有应变 能函数的材料。这种弹性体应力可通过不同的与 应变有关的应变能密度函数来表达。常用的应变 能密度函数模型有基于应变不变量的门尼里夫林 (Mooney—Rivlin)模型，基于主拉伸比的Ogden模 型等。其他还有以分子统计学为基础的Gent模 型和Arruda—Boyce模型等，后两种模型由于易于 通过简单的单轴 向试验就能确立，近来也常常被 使用。以上这些模型都是假设材料不可压缩的， 即各向同性的，其本构方程通常被用于有限元分 图1 齿形切边V带 (箭头为纤维方向) 析。而在 自然界，许多材料 (如哺乳动物的心肌 等)不仅有不可压缩性和非线性，而且还具强烈的 本文提出了一种多项式应变能函数，可以很 各向异性的特征，以上这些模型并不适用，须用各 好处理从拉伸到压缩的切向刚度之间的差异性 ， 向异性本构方程 。 确定了纤维增强材料的材料常数 。最后还探讨 了 除此之外，在许多机械工程领域，如广泛应用 模拟传动带工作状态的弯曲试验，以图示方式 比 于汽车、OA机器和机床等传动的传动带 (见图 较了数值模拟与实验结果之间的误差。 1)，其结构复杂，除橡胶外还含有织物和芯线，也 具有各向异性的特征 。V带是被锲人皮带轮 中使 2 短纤维增强超弹性体材料的本构方程 用 ，其在宽度方 向上会受到很大的压缩力 ，因此在 ¨u给出了基于假设各 向同性变形不变量各 其橡胶部分常常加入定向短纤维，纤维方 向为带 向同性的切 向模量的本构方程 ；Ogden给出了基 第 1期 吴贻珍 ．短纤维 ／橡胶复合材料 的数值模拟 25 于拉伸 比的本构方程 ；Ishikawa等给出基于变形不 ha=Fa0 (13) 变量各 向异性切 向模量的本构方程 。在这一节 中 拉伸 比A与初始方 向a。、变形方向a的关系 首先简单介绍一下基于不变量的横 向各向异性的 为 ： 本构方程 。 14=aoCa0=Fa0Fa0=AaAa=A aa=A 2．1 Green—Lagrange变形张量