１ １４ 液压与气动 ２０１３年第５期 ］Ｊ ＤＯＩ：１０．１１８３２／ｊ．ｉｓｓｎ．１０００－４８５８．２０１３．０５．０３２ 旋转轴唇形密封圉的有限元分析与仿真 桑建兵，邢素芳，刘宝会，周婧，卢宸华 Ｆｉｎｉｔｅ Ｅｌｅｍｅｎｔ Ａｎａｌｙｓｉｓ ａｎｄ Ｓｉｍｕｌａｔｉｏｎ ｏｎ Ｒｏｔａｒｙ Ｓｈａｆｔ Ｌｉｐ Ｓｅａｌ ＳＡＮＧ Ｊｉａｎ－ｂｉｎｇ，ＸＩＮ Ｓｕ—ｆａｎｇ，ＬＩＵ Ｂａｏ—ｈｕｉ，ＺＨＯＵ Ｊｉｎｇ，ＬＵ Ｃｈｅｎ—ｈｕａ （河北工业大学机械工程学院，天津３００１３０） 摘要：利用橡胶类材料单轴拉伸所得到的应力应变曲线，通过数据拟合确定了Ｍｏｏｎｅｙ．Ｒｉｖｉｌｉｎ模型的 材料参数。建立了旋转轴唇形密封圈的有限元计算模型，模拟了弹簧圈预紧、过盈配合以及在不同的介质压 力对旋转轴唇形密封圈密封性能的影响，得到了橡胶密封圈的范・米塞斯应力的分布规律以及唇口处接触 应力的分布曲线。研究结果表明：随着介质压力的增加，唇口处范・米塞斯应力随之增加。唇口处接触应力 的分布近似为二次抛物线，接触应力的最大值出现唇口尖端处。随着介质压力的增加，唇口处接触应力的应 力峰值和接触宽度明显增加，且峰值均大于介质压力，满足旋转轴唇形密封圈的密封条件。 关键词：旋转轴唇形密封圈；Ｍｏｏｎｅｙ—Ｒｉｖｉｌｉｎ模型；非线性有限元分析；仿真 中图分类号：ＴＨ１３７；ＴＢ４２文献标志码：Ｂ文章编号：１０００４８５８（２０１３）０５－０１ｌ４＿（）４ 引言 用，其工作的主要机理是通过油封与轴的接触来防止 橡胶类材料在航天航空、机械以及汽车等各领域 润滑油或其他介质的泄漏。接触宽度和密封界面上接 有着非常广泛的应用 。由于在变形过程中具有材料 非线性、几何非线性以及接触非线性的典型特征，使得 收稿日期：２０１２－１１－０１ 橡胶类材料的有限元分析成为一个高度非线性问题。 基金项目：河北省自然科学基金项目（Ｅ２０１１２０２１１５）；天津市 旋转轴唇形密封圈，又称油封，因其成本低廉、结 自然科学基金项目（１２ＪＣＹＢＪＣ１９６００） 作者简介：桑建兵（１９７４一），男，河北邢台人，副教授，博士， 构简单以及安装方便等优点，在工业中得到了广泛应 主要从事橡胶类材料的有限变形理论及数值模拟的研究工作。 一５０． ｃｌｅｌＤ＝１８５８３． 崔勇，崔松彪．踝足矫形器对４０例脑卒中后偏瘫患者下 ［９］郭忠武，王广志，丁辉，等．正常青年人步态运动学参数的 肢运动功能的影响［Ｊ］．延边大学医学学报，２００８，（１）： 研究［Ｊ］．中国康复理论与实践，２００２，８（９）：５３２—５３３，５４４． ３７—３８． ［１Ｏ］ 肖亮子，韩建海，赵书尚，等．基于气动人工肌肉并联驱 Ｄａｎｉｅｌ Ｐ．Ｆｅｒｒｉｓ，Ｋｅｉｔｈ Ｅ．Ｇｏｒｄｏｎ，Ｇｒｅｇｏｒｙ Ｓ．Ｓａｗｉｃｋｉ．Ａｎ 动手腕康复训练器［Ｊ］．液压与气动，２００７，（６）：６４—６６． Ｉｍｐｒｏｖｅｄ Ｐｏｗｅｒｅｄ Ａｎｋｌｅ－－ｆｏｏｔ Ｏｒｔｈｏｓｉｓ Ｕｓｉｎｇ Ｐｒｏｐｏｒｔｉｏｎａｌ Ｍｙ－－ ［１１］李继磊，韩建海，赵书尚，等．气动人工肌肉驱动关节 ｏｅｌｅｃｔｒｉｃ Ｃｏｎｔｒｏｌ［Ｊ］．Ｇａｉｔ＆Ｐｏｓｔｕｒｅ，２００６，（２３）：４２５—４２８． ＰＩＤ位置控制［Ｊ］．河南科技大学学报（自然科学版）， Ｓ．Ｑ．Ｘｉｅ，Ｐ．Ｋ．Ｊａｍｗａ１．Ａｎ ｉｔｅｒａｔｉｖｅ Ｆｕｚｚｙ Ｃｏｎｔｒｏｌｌｅｒ ｆｏｒ ２００８，（３）：６３—６６． Ｐｎｅｕｍａｔｉｃ Ｍｕｓｃｌｅ Ｄｒｉｖｅｎ Ｒｅｈａｂｉｌｉｔａｔｉｏｎ Ｒｏｂｏｔ［Ｊ］．Ｅｘｐｅｒｔ ［１２］ 隋立明，赵铁，张立勋．气动肌肉驱动关节的位置控制策 Ｓｙｓｔｅｍｓ ｗｉｔｈ Ａｐｐｌｉｃａｔｉｏｎｓ，２０１１，（３８）：８１２８—８１３７． 略研究［Ｊ］．液压与气动，２００８，（６）：１１—１３． 余伟正．步行康复用踝关节系统的研究［Ｄ］．上海：上海 ［１３］ＧＢ １００００—８８，中国成年人人体尺寸［Ｓ］． 大学。２００９． ［１４］ＳＭＣ（中国）有限公司．现代实用气动技术［Ｍ］．北京：机 张今瑜，董玉红，王岚，等．国家自然科学基金．下肢康复 械工业出版社，２００４．１． 训练机器人关键技术研究［ＥＢ／ＯＬ］．２００６—８—３０．ｈｔ． ［１５］ 张贵兰，耿德旭，刘晓敏，等．伸长型气动人工肌肉轴向 ｔｐ：／／ｇｏｎｇｘｕｅ．ｃｎ／ｊｉｄｉａｎｘｕｅｙｕａｎ／ＳｈｏｗＡｒｔｉｃｌｅ．ａｓｐ？Ａｒｔｉ— 力学特性分析［Ｊ］．机械工程师，２０１０，（１２）：７—８． ２０１３年第５期 液压与气动 ｌ １５ 触压应力的分布是影响油封密封性能的重要参数，但 要得到其精确的理论解析解是十分困难的。随着数值 计算方法、材料学以及大型有限元分析软件的发展，使 得利用非线性有限元软件对油封在安装和使用中的高 度非线性接触问题进行研究成为可能，并取得了一系 列研究成果。 的物理参数，可以在理论上阐明材料变形的物理含义， 但通过实验或数据来拟合和确定这些参数是非常困难 的。物理参数越多，本构模型的形式越复杂，在实际应 用时就越困难，这通常是这类本构理论模型难以在实 际中进行应用的重要原因之一。为了简化，取Ｒｉｖｌｉｎ 模型的第一项，满足此应变能函数的材料为Ｎｅｏ． 本文基于非线性有限元软件ＭＳＣ．Ｍａｒｃ，并结合 单轴拉伸试验，确定了橡胶类材料Ｍｏｏｎｅｙ—Ｒｉｖｉｌｉｎ模型 Ｈｏｏｋｅａｎ材料，而取Ｒｉｖｌｉｎ模型的线性组合，满足此应 变能函数的材料为Ｍｏｏｎｅｙ．Ｒｉｖｌｉｎ材料 Ｊ，其本构模型 的本构参数。利用所得到的本构参数对油封的密封性 能进行了数值模拟，研究了介质压力对接触应力、接触 宽度的影响，为进一步对油封力学性能的研究及优化 提供了理论上的技术参考。 １ 橡胶类材料的本构关系 在常温情况下，橡胶类材料为各向同性超弹性材 料，其应变能函数可表示为： Ｗ＝Ｗ（， ，，２， ） （１） 其中：，ｌ、，２和，３为Ｃａｕｃｈｙ应变张量的三个不变量。 ， ＝Ａ ＋Ａ；＋Ａ； ，２＝Ａ Ａ；＋Ａ；Ａ；＋Ａ；Ａ ，３＝Ａ２１＾２２Ａ； （２） 其中：Ａ 、Ａ：和Ａ 为主延伸率。 对于不可压缩的橡胶类材料，厶＝１，则Ｃａｕｃｈｙ应 力张量可表示为： 一ｐ，＋２ 一２ （３） 其中：Ｐ表示为静水压力；Ｂ为左Ｃａｕｃｈｙ．Ｇｒｅｅｎ变形 张量。 令 ＝ ， ＝ ，则ＣａｕｃｈＹ应力分量的表 达式可表示为： ＝２Ａ ＷＩ一２Ａ ，２一Ｐ （４） 其中：Ｐ为静水压力。 国内外学者从连续介质力学、热力学、相变唯象理 论、统计物理学等多学科、多领域出发，建立了许多描 述橡胶类材料的物理行为的本构模型。从唯象角度， １９４８年Ｒｉｖｌｉｎ对于各向同性超弹性材料提出了应变 能函数模型为 Ｊ： ＝ｉＪ＝Ｏ ∑ｃｊｊ（， 一３） （，２—３） （５） 其中：Ｃｉｉ为材料常数，， 和，２分别为左Ｃｕａｃｈｙ—Ｇｒｅｅｎ 变形张量的第一和第二不变量。在本构模型中含有许 多物理参数，使得模型过于复杂。研究表明，引入较多 的表达式如下： Ｗ＝Ｃ１０（， 一３）＋ 。（，２—３） （６） １９７２年，Ｏｇｄｅｎ＿４ 提出了主伸长比Ａｉ作为应变能 函数自变量新的应变能函数： ＝ ［ｊ－￣ － （ｘ７”＿＋Ａ；“＋Ａ； ）】－３ ＋４．５Ｋ（ 一１） （７） 其中： ｉ和 ｊ为材料常数， 为初始体积模量，＇，为体 积率。 Ｊ＝Ａ１Ａ２Ａ ３ （８） 对于不可压缩的橡胶材料，Ｊ＝１且最后一项消 失，此时应变能函数可以改写为： Ⅳ ，． Ｗ＝∑ ［（Ａ ”＋Ａ；“＋Ａ； ）一３ ３ （９） 二阶Ｍｏｏｎｅｙ—Ｒｉｖｉｌｉｎ模型能更好地反映橡胶类材 料的受力状态，且材料参数容易确定，在单轴拉伸、纯 剪切和等比双轴拉伸实验中都得到了验证，是目前准 确性比较高的的不可压缩弹性体的材料模型。因此本 论文采用了二阶Ｍｏｏｎｅｙ—Ｒｉｖｉｌｉｎ模型，利用单轴拉伸所 得到橡胶类材料的应力应变曲线，拟合了橡胶类材料 的本构参数，并将拟合结果应用于工程实例中的油封 大变形接触变形分析，为油封的设计与优化提供有益 的参考依据。 ２单轴拉伸实验及材料参数的确定 采用哑铃型标准试件，按照国家标准采用电子拉 伸实验机对橡胶类材料进行单轴拉伸实验，所得到应 力应变曲线如图１所示。 应变 图１单轴拉伸应力应变曲线 １１６ 液压与气动 ３有限元分析模型 ２０１３年第５期 对于单轴拉伸， 。＝ ， ＝ ，＝０，考虑橡胶类 材料为不可压缩材料，因此可取主延伸率Ａ。＝Ａ，Ａ： 旋转轴唇形密封圈由于其结构、外载荷及约束条 ：Ａ ＝１／√Ａ。则应变不变量分别为： ｆ，１＝Ａ２＋２Ａ一 件关于其中心旋转轴对称，因此进行有限元分析的时 候可以将旋转轴唇形密封圈简化为一轴对称问题，选 择其横截面进行有限元计算＿６］。图３为油封的装配 示意图。为了减小误差，橡胶部分采用四节点的四边 形赫曼单元，金属骨架以及弹簧圈均采用四节点四边 形全积分单元，其有限元网格如图４所示。为了防止 ｛【 厶＝２Ａ＋Ａ一 ：（１０） １ 由表达式（４）可得： ｆ １＝一Ｐ＋２Ａ２Ｗ１—２Ａ （１１） 油封在弹簧圈预紧力、介质压力以及过盈配合的作用下 Ｌ ２＝ ３＝一Ｐ＋２Ａ １—２Ａ 由于 ２＝ ３：０，由此可得： Ｐ＝２Ａ～Ｗｌ一２Ａ （１２） 则拉伸应力 表示为： ＝２（Ａ 一Ａ－１）（ ＋Ａ＿１ ｒ２） （１３） 由Ｍｏｏｎｅｙ—Ｒｉｖｌｉｎ本构本构模型（６）可知： ＝ ＝Ｃｌｏ ＝ ＝Ｃｏ。 （１４） 由此可得： ＝２（Ａ 一Ａ ）（Ｃ１０＋Ａ～Ｃ０１） （１５） 实验表明单轴拉伸、单轴压缩、双轴拉伸以及平面 剪切等状态下的应力应变关系，能较好的反映橡胶类 材料的本构关系 Ｊ。为了得到Ｍｏｏｎｅｙ—Ｒｉｖｉｌｉｎ橡胶类 材料模型的材料参数，采用标准试件以及标准的实验 方法，对橡胶胶类材料进行了单轴拉伸实验，得到了在 单轴拉伸模式下的应力应变关系数据曲线，将数据导 入非线性有限元软件ＭＳＣ．Ｍａｒｃ中，与程序提供的单 轴拉伸等标准曲线进行拟合，得到的拟合曲线如图２ 所示。 雳一图 从图２中可以看出，实验所得到单轴拉伸曲线与 Ｍｏｏｎｅｙ—Ｒｉｖｉｌｉｎ本构模型的标准曲线吻合的较好，因此 Ｍｏｏｎｅｙ—Ｒｉｖｉｌｉｎ模型能够较好地描述橡胶类材料的变 形行为。 经过回归分析及数据处理后所得到橡胶类材料的 材料参数如下： Ｃ１０＝０．１４５ ＭＰａ，Ｃｏｌ＝０．５５８ ＭＰａ （１６） 有限元网格出现畸变，有限元分析过程中采用网格自适 应以及网格重划分的方法，以减小有限元的计算误差。 图 图４油封的有限元模型 本模型采用直接约束法来处理有限元分析中的接 触问题。通过追踪物体的运动轨迹，一旦探测出发生 接触，便将接触所需的约束和节点力作为边界条件直 接施加在产生接触的节点上。由于相对于橡胶来说轴 的变形很小，因此在进行有限元接触分析的过程中将 轴定义为刚体， 定义了橡胶与骨架的接触，橡胶与弹簧 圈的接触以及橡胶与轴的接触。 在加载过程中，首先施加弹簧圈的预紧力，然后利 用刚体的强制位移，将旋转轴的的竖直方向的位移作 为油封与旋转轴之间的过盈配合量，以此来模拟油封的 安装过程，使油封处于一定的压缩状态 。然后在油封 的右侧逐步施加介质压力，使其达到最终的工作状态。 ４接触问题的有限元分析 利用ＭＳＣ．Ｍａｒｃ软件对不同介质压力下油封的受 力和变形进行分析。图５为介质压力为０．０３ ＭＰａ时 油封整体的范．米塞斯应力（Ｖｏｎ—Ｍｉｓｅｓ）分布云图。从 应力云图上可以看出，最大应力主要集中在金属骨架 及弹簧圈上，可以得出在介质压力作用下，油封的应力 主要集中作用在金属骨架以及弹簧圈上，而对于橡胶
