橡胶密封元件轴对称结构在ABAQUS中的数值模拟---HOFF.doc

橡胶密封元件轴对称结构在ABAQUS中的数值模拟 和强 王宇 火进 （北京霍夫技术服务有限公司 北京 100083） 摘 要：运用ABAQUS/Standard模拟了橡胶密封元件装配过程中的受力特性；实现了用数值方法解决包含超弹性材料，轴对称结构，摩擦接触以及多步骤分析的问题；为复杂工程问题的数值解法提供了参考价值。 关键词：橡胶密封 超弹性 轴对称 多步骤分析 引言 橡胶是在使用温度下处于高弹态的高分子材料，变形中表现出很强的几何物理非线性，与其他材料最基本的区别是其弹性模量低，具有很高的伸缩性和储能能力，因此，被广泛的用来制成密封，减振，防护等用品。传统的橡胶元件的设计方法是根据试验和试验所取得的规律进行，造成产品开发周期长，成本高，复杂结构或工况试验难以进行等后果。近年来，计算机仿真技术和现代非线性理论的快速发展为设计人员提供了良好的基础。本文借助有限元理论和分析软件ABAQUS对汽车转向轴橡胶密封圈的安装过程进行了研究。 油封结构简介 目前中国已经标准化的产品有内包骨架油封（B型），外露、半外露骨架油封（W型）， 装配式油封（Z型），有副唇内包骨架油封（FB型），有副唇外包骨架油封（FW型），有副唇装配式油封（FZ型）以及流体动力型旋转油封等等。其中外露、半外露骨架油封具有油封定位准确，同轴性好，安装性能高，摩擦升热小，导热性能好和材质消耗小等优点。 油封结构包括结构形式和结构参数两部分。结构形式（油封结构和几何形状）是由配合及安装要求，密封介质及轴的旋转方向等使用要求来决定；而结构参数的内容包括工作情况参数，性能参数，胶料性能参数和剖面结构参数等。 橡胶变形特点及超弹性本构模型 2．1 橡胶变形特点 橡胶材料在承受拉压过程中，体积的变化量很小，可以忽略不计，但是作为超弹性体，其最明显也是最重要的物理特性是在较小的外力作用下就能产生很大的变形，其伸长率可达100％～1000％。 2．2 超弹性本构模型 弹性材料的变形过程是可逆的，如无其他不可逆伴随，单纯的弹性变形过程的熵产率为零，也就是单位质量的热力学能等于单位质量的应变能，对??等温过程，单位质量的自由能便是单位质量的应变能，存在应变能的材料称为超弹性材料，因此橡胶材料属于超弹性材料。一般认为，橡胶材料是各项同性不可压缩的超弹性体，它的物理属性主要通过应变能函数表达。一旦确定了应变势能函数W的形式，柯西应力张量就可以由下式给出： （3－1） 式中：单位张量；为由于不可压缩假设引入的水静压力；为的不变量。 在大多数连续介质力学处理橡胶弹性时，都把橡胶材料的变形看成是各项同性的超弹性的均匀变形，这样应变能密度函数就可以表示成主伸率或者是变形张量的三个不变量的函数。在诸多的应变势能模型中，本文采用Ogden模型，其表达式如下： （3－2） 为应变势能；为主伸长率；指数不一定是整数； 必须有； Ogden模型是有限元分析中较常用的以伸长率表示的应变势能模型，它的优点在于，简单直观，且能够用合理的拟合系数较精确地表达橡胶材料的力学特性。 ABAQUS应用应变势能来表达超弹性材料的应力－应变关系，而不是用杨氏模量和泊松比。 3，有限元模型建立 3．1超弹性 ABAQUS模拟超弹性材料时，提出如下假设： （1）材料行为是弹性； （2）材料行为是各项同性； (3) 考虑几何非线性效应。 当模拟超弹性材料时，用户可以提供弹性势能表达式所需要的参数，也可以提供模拟材料特性的试验参数。ABAQUS可以接受试验数据并能够应用最小二乘法拟合出材料参数。为了得到最佳材料模型，使计算更加精确和稳定，本论文提供了单轴拉压，等双轴拉压以及平面拉压（纯剪）的材料实验数据输入到ABAQUS中。有限元模型如图3－1所示： 图3－1： 橡胶密封圈有限元模型（及其局部放大图） 3．2 轴对称性 由于模型的几何形状和载荷均是轴对称的，所以为了减小节省计算时间，降低对存储的需求，允许使用更为精细的网格，这里采用轴对称单元进行模拟，只需要模拟通过构件的一个平面：每个单元代表一个完整的3600圆环。模型如图3－2所示： 图3－2 轴对称结构密封圈模型 3．3 摩擦接触 有限元分析中，接触条件是一种特殊类型的不连续约束，它允许力从一部分传递到另一部分。ABAQUS/Standard的接触模拟可以是基于表面（surface）或基于接触单元（contact element）的，设置步骤：（1）创建接触表面（2）创建接触对（3）定义控制接触面之间相互作用的本构关系。当表面发生接触时，在接触面之间一般传递切向力以及法向力。这样在分析中就要考虑阻止表面之间的相对滑动。本文采用库仑摩擦模型，用摩擦系