基于ANSYS有限元的微动疲劳接触特征分析.pdfVIP

第 11卷 第 5期 重庆科技学院学报 自然科学版 2009年 l0月 基~ANSYS有限元的微动疲劳接触特征分析 孙 明俊 刘 玉金 河南理工大学万方科技学院，焦作 454003 摘 要 ：接触压力是影响构件微动疲劳寿命 的关键 因素之一 。微动疲劳实验 中，桥式微动装置可 以合理控制接触压 力，以便分析接触压力对构件疲劳寿命的影响规律 。基于数值分析法，建立桥式微动装置 的有限元模型来研究接触 特征 ；考虑微动桥 自身弯曲的影响，分析 了接触副上节点的应力应变分布规律，提 出了特征接触 区域 的判断方法。 关键词 ：微动疲劳 ；数值分析 ；接触区域 ；接触压力 中图分类号：TH123 文献标识码 ：A 文章编号：1673—1980 2009 05—0148—04 机械系统和工程设备服役于交变载荷作用下 ， 分析 了微动疲劳损伤机制 ，探讨 了磨屑形成与演变 可导致构件接触表面之间相对微小幅度的滑动 ．这 过程及其对微动疲劳行为的影响，考察了摩擦系数 种重复的滑动称为微动。同时交变载荷和微动能促 随时间的变化规律 使疲劳裂纹早期萌生和加速扩展 。最后导致构件在 但在此装置中，微动桥和试样间的接触条件不 大大低于材料疲劳极限，甚至在低于材料弹性极限 好描述，特别是当微动桥本身有弯曲时这一问题尤 时失效 ，这一现象被称为微动疲劳 ，它是指构件在有 为突出。上述研究均没有考虑微动桥本身的弯曲对 微动磨损时的疲劳强度和疲劳寿命问题 1J【。影响构件 微动接触的影响。要使分析结果更符合实验结果，必 微动疲劳寿命 的因素很多，接触压力P是影响构件微 须考虑微动桥 自身弯曲对微动接触区域的影响。文 动疲劳寿命的关键因素之一嘲。在实验过程 中 通过 中在部分滑移状态条件下，提 出了接触表面节点的 改变P来研究其对构件微动疲劳强度的影响规律。 应力应变分析方法和接触区域的判断方法。在考虑 为简易地控制P的变化 ．桥式微动试验装置常用于微 微动桥本身弯曲时 ，为基于桥式微动装置实验的数 动疲劳试验研究I一。图l为试验装置结构示意图。 值分析提供了有效的方法 。 1 计算模型建立 模型材料45钢 ，弹性模量E 210GPa，v 0．3。由 于模型和载荷关于 ．y轴对称 ，可以取其1／4模型进 行分析。由于1／4模型的对称截面在其法 向是不动 的．所以在对称截面上节点法 向位移约束为0。模型 选用plane42单元 。接触面为conta172．目标 面为 图1 微动疲劳 实验结构示意 图 targe169单元；为了节省计算时间和计算机资源，又 微动桥两足对称地夹在试样上 ，试样上的循环 要满足接触计算要求 ．所以在接触面上及附近网格 应变引起微动桥足与试样间产生相对运动。通过此 较细，离接触面较远部位 网格较粗 。对于接触模型， 装置，微动疲劳试验可以在普通疲劳试验机上进行。 接触面节点与 目标面节点最好成对出现 ，形成接触 刘道新等人5[1运用桥式微动试验装置研究了不 同的 对 。网格模型共有节点l217个 ，单元4498个 ，接触 微动桥足圆角半径对微动疲劳中磨损状况的影响规 面上单元长度0．1mm，共40个单元 ，40对接触节点， 律。基于运用桥式微动试验装置进行的试验 ．卫 中 41个节点。为了模拟垂直载荷 ，在微动桥平行y轴对 山等人6】『通过观察微动 区的磨损特征和截面形貌 ， 称截面节点上施加集 中力 ．在试样右表面上施加轴 收稿 日期 ：2009—03—40 作者简介 ：孙明俊 1980一 ，男，山东青岛人，硕士，河南理工大学万方科技学院讲师，研究方 向为力学。 · 148· IuI l ＆ ^J【 III三 E惦0+ 0叫+ F0∞+ 0眦+ E0叭+ 0∞+ EO叭一 孙明俊 ．刘玉金 ： ff-ANSYS有限元的微动疲劳接触特征分析 向载荷 二维微动疲劳有限元网格、边界条件 对称 约束 和载荷 Jp 一2880N，Ocyli,．一480MPa 如图2所 示 为了研究两接触表面相对滑移距离 ，这里用相 同的单元长度来划分接触面，并使每个节点上下成 对出现在同一位置 ，形成接触对 。由于模型为接触 非线性有限元计算 ，选择大变形静态分析 ，100个载 荷步 ．其收敛情况如图3所示。 T丽 4-1l5l721~129．013l 8i _444l 一 I1．99 ． 图4 试样接触部分y方向应 力云图 l 一f 、 ＼ ＼ 、 、 ＼ 图2 计算模型网格、边界条件及载荷图 ＼ Tin1e：100 、、、～ DIST ／ | 图5 试样接触表面上节点沿接触路径 l，方向应 力曲线 、 ＼ I^^^。^ 0 、 ＼ 柙 ￥ — 4l899～ ＼ — 83．804