实际结合面的法向接触刚度多尺度计算方法.docxVIP

实际结合面的法向接触刚度多尺度计算方法 孙 伟，黄 信，孙志勇，黄伟强 （大连理工大学 机械工程学院，辽宁 大连 116024） 摘 要：针对结合面接触刚度解析计算与实际接触状态差异较大，基于有限元的大平面 接触刚度计算由于计算量大难以实现的问题，提出一种结合面法向接触刚度多尺度计算方法。该方法在应用测量仪器获取结合面平面度、波纹度及粗糙度形貌特征的基础上，应用小波分解技术获取真实形貌特征，并基于有限元微观接触分析建立局部接触刚度与压强之间的关系，基于有限元宏观接 触分析获取结合面压强分布， 宏 - 微观结合计算获得结合面法 向接触刚度。该方法具有能够真实反映结合面接触状态、计 算效率高等特点，法向接触刚度试验结果表明，该方法是一 种计算结合面法向接触刚度的有效方法。 关键词：结合面； 表面形貌；接触刚度；多尺度；有限元 1 引言 机械装配体 中广泛涉及平面或曲面接触，结合面的接触刚度对机械系统 的静动特性有着重要的影响 [1] 。零件表面质量受加工过程中 诸多因素影响，不可能是理想光滑表面，而是由不同尺寸和 形状的凸峰和凹谷组成的，其形貌特征按尺度从小到大分别 对应于粗糙度、波纹度和形状误差。准确计算结合面接触刚 度，需要考虑界面不同尺度形貌特征及相应的接触特性。 对 于粗糙表面接触特性，学者们从理论上进行了大量的研究。 文献 [2] 率先提出了 GW 统计接触模型，文献 [3] 提出了 CEB 弹塑性接触模型， 文献 [4] 提出了 MB 分形接触模型。 近年来学者们改进了 GW 、CEB 及 MB 模型，提出一系列粗糙表面 弹塑性接触问题的解析模型 [5-6] ，取得了不少成果。但现有解析模型中针对粗糙表面微凸体分布规律及形状的假设、简化不可避免地导致计算结果与实际接触状态存在一定差异。 此外，解析模型主要针对微观形貌接触问题， 而对形状误差、 波纹度等宏观形貌特征的影响缺乏细致研究。 随着计算机技术的发展，有限元法在结合面接触分析中得到广泛应用。文献 [7] 针对分形粗糙表面与刚性光滑平面之间的弹塑性接 触问题进行了有限元分析。 文献 [8] 基于弹塑性理论对具有粗糙表面的长方微元体进行有限元接触分析，给出了根据受力 和变形关系计算粗糙表面接触刚度的方法。文献 [9] 通过 3D 粗糙表面的数字化表征方法构建两粗糙表面接触的有限元 模型并分析了粗糙表面接触特性。文献 [10] 用真实粗糙表面 形貌数据进行了有限元仿真，模型与实际情况更加吻合。有限元方法虽然克服了解析模型的一些局限性，但要准确描述表面形貌，就必须对模型特别是接触界面进行细致的单元划分，计算时间、内存等限制了计算规模，通常仅适用于微小接触面的计算，而微小接触面不能完整反映机床结合面、拉 杆转子接触界面等实际表面的形貌特征。 针对以上研究现状，基于实际表面的形貌数据，应用小波分解技术获取真实 形貌特征，进行微观和宏观两种尺度的有限元弹塑性接触分析，通过提取微观尺度的局部法向接触刚度和宏观尺度的接触压强分布，最终获取包含实际表面完整形貌信息的结合面 法向接触刚度，并进行了实验验证。 2 法向接触刚度多尺 度计算方法 机械零件表面是由粗糙度、波纹度和形状误差三种不同尺度的形貌特征综合影响产生的结果。实际表面的接触实际上是由许多局部微观接触构成，决定宏观接触特性的形貌特征主要是波纹度及形状误差，而决定微观接触特性 的形貌特征则主要是粗糙度形貌，如图 1 所示。结合面总接 触刚度可以认为是单位面积的局部接触刚度在接触面内的 积分，而一般认为是接触压强的函数，则有： 图1实际表 面接触模型 Fig.1Contact Model of Actual Interface  提出的结合面法向接 触刚度计算方法，将粗糙度形貌特征归为微观尺度特征，将 波纹度及形状误差形貌特征归为宏观尺度特征，进行两种尺 度的有限元接触分析。微观尺度计算作为接触微元为宏观尺 度传递局部接触刚度细节信息，而宏观尺度的计算为微观尺 度提供接触压强边界条件，多尺度间递进，最终获取基于实 际表面形貌的结合面法向接触刚度，主要步骤如下： （ 1） 用表面轮廓仪测试获取微观、宏观两种尺度的结合面形貌， 并用小波分解进行形貌特征提取； （ 2）将形貌数据以命令 流导入 ANSYS 中，用直接建模方法建立有限元接触模型， 进行微、宏观两种尺度的有限元弹塑性接触分析； （3）提 取微观尺度接触体的一系列平均接触压强 p 及相应的法向变 形δ ，以下的幂函数形式 [1] 进行拟合。 式中： c、m—待定 系数。则单位面积的局部法向接触刚度为： 式中 （4）提 取宏观尺度接触体的结合面各节点的压强，对于离散的节 点，若节点在 XY 平面内是均匀分布的， 则根据式 （ 1）可求 出单位面积结合面的总法向接触刚