[发明专利]一种基于亚塑性理论的机械关节接触面本构模型创建方法

说明书

本发明实施例公开了一种基于亚塑性理论的机械关节接触面本构模型创建方法，涉及机械结构动力学以及金属摩擦学领域。该模型视关节界面为一种虚拟的材料，并建立在Kolymbas提出的亚塑性框架和连续介质假设的基础上。首先，根据金属接触面的特征，简化了一般的亚塑性公式。然后，通过考虑初始刚度比、法向接触定律和库仑破坏面，将通用公式中的待定系数与实际中广泛使用的物理参数关联起来。推导出的本构模型具有简单的数学公式，且仅包含四个物理参数。该模型可有效表征接触界面的粘着、微观滑移和宏观滑移行为以及各方向间的耦合行为，并具有压力敏感性。因此，所提模型可以更精确地预测任意多轴加载下机械关节接触面的行为。

技术领域

本发明属于机械结构动力学以及金属摩擦学领域，具体涉及一种基于亚塑性理论的机械关节接触面本构模型创建方法。

背景技术

机械结构通常通过紧固件连接组装，如螺栓连接、铆接连接和过盈配合连接。对于单个部件，有限元分析可以准确预测响应；然而，对于装配结构来说，由于无法准确描述关节非线性行为，要准确预测响应还存在较大的挑战。在实验观测中，经常可以看到由关节非线性引起的频率偏移、幅变阻尼以及超谐波共振等特征行为。这些特性的先验知识对于设计过程中避免运行过程中的意外共振至关重要。因此，开发一种能够精确描述机械关节非线性行为的模型具有相当大的实际需求。近几十年来，机械关节的本构模型受到了广泛的关注，这些模型可以分为基于物理驱动的模型以及现象学模型两类。前者基于微凸体接触理论，并由微凸体相互作用的统计总和导出接触力，其参数依赖于粗糙表面的力学性能和统计参数。对于现象学模型，关节接触面行为由一参数化的数学公式描述，通过对实验数据的拟合可获得其参数。其中，Iwan模型因其再现接触面微滑移行为的能力而广受欢迎。在Iwan模型基础上已由许多学者提出了大量的现象学模型，其对接触面行为的不同特征有不同的侧重点。这些模型涉及的层面包括粘着、微滑、宏滑、拍打、微冲击、切向与法向运动的耦合效应、部分卸载时的局部环以及大载荷下的钉扎现象。然而，考虑的行为越多数学公式也变得更加复杂，且材料参数的数量随着模型的复杂程度而显著增加，使得响应计算和参数识别变得困难。在实际应用中，最好采用一种能捕捉主要特征、数学公式简单、参数少的本构模型。此外，属于这两类的本构模型大多是1维或2维模型，不能考虑接触面中三个运动分量之间的耦合和相互作用。因此，开发一个能够考虑不同运动部件之间耦合的3维本构模型是必要的。

发明内容

基于上述动机，本文基于亚塑性框架建立了的机械关节接触面三维本构模型创建方法，推导出的模型属于现象学模型。

为达到上述目的，本发明的实施采用如下技术方案：

本发明实施例提供一种基于亚塑性理论的机械关节接触面本构模型，其包括：

S1、根据亚塑性理论得到满足设定限制条件且具有滞回行为的接触面材料的亚塑性本构方程的一般形式；

S2、不考虑接触面材料的面内应力及应变的条件下，将亚塑性本构方程的一般形式写为矩阵形式；

S3、根据机械关节接触面的特征，由其法向力学行为的规律确定亚塑性本构方程中的与法向行为相关的系数；

S4、根据机械关节接触面的特征，由其切向力学行为的规律确定亚塑性本构方程中的与切向行为相关的系数；

S5、将由切法向行为确定的系数代回亚塑性本构方程最终得到描述机械关节接触面材料的本构模型。

步骤S1根据亚塑性理论得到满足设定限制条件且具有滞回行为的接触面材料的亚塑性本构方程一般形式的具体方法如下：

将机械关节的接触面区域建模为具有厚度b，b→0，的连续介质，根据亚塑性理论，其本构模型通过满足如下关系的张量函数H进行定义

其中，σ和ε分别为接触面材料的应力和应变，且顶标·表示物理量的时间导数；