润滑科技信息平台

专利摘要

一种考虑涂层影响的点接变形与力载关系的建立方法，该方法包括：(1)采用半解析法或者有限元法计算获得固体润滑摩擦副在相应接触载荷下的接触变形；(2)应用最小二乘拟合获得考虑固体润滑涂层影响接触副的接触载荷与接触变形关系式。该方法克服了现有基于Hertz接触理论无法考虑固体润滑涂层影响的局限性，使固体润滑点接触摩擦副接触变形的计算结果更符合实际情况，保证了接触变形计算结果的准确性，为有固体润滑涂层的滚动轴承和齿轮力学特性分析提供了有效的求解方法，使滚动轴承和齿轮等涉及涂层接触问题的模型更符合实际情况，其计算结果更可靠可信。 ......

基本信息

专利类型:
发明专利
申请/专利号:
CN202110176144.2
申请日期:
2021-02-06
专利申请人:
["天津职业技术师范大学","哈尔滨工业大学"]
分类号:
G06F30/23 ; G06F119/14
发明/设计人:
王廷剑，吴跃，张传伟，李臻，翟晗，战利伟，段睿，李东昊，王黎钦，古乐，郑德志，赵小力，赵阳
权利要求: 1.一种考虑涂层影响的点接触力载与变形关系的建立方法，其特征在于：所述方法具体步骤为：步骤一：将接触载荷力载范围0～Qmax均匀划分为n+1个载荷点，第i个载荷点接触载荷Qi为：Qi＝(i-1)ΔQ (1)其中：1≤i≤n+1，Qmax为接触载荷力载范围的最大值，ΔQ＝Qmax/n，为相邻载荷点接触力载增量；步骤二：采用基于半解析法或有限元法的接触分析模型求解固体润滑摩擦副在各个点接触载荷Qi下的接触变形δi；基于半解析法的接触分析模型计算在接触载荷Qi下的接触变形δi，其具体过程如下：对于三维点接触问题，选取一个矩形区域Ω＝[-2aH,2aH]×[-2bH,2bH]作为计算域，并将其划分成Nx×Ny个大小相等的表面矩形单元，并以矩形单元的几何中心作为该单元的节点，Ω为分别以4aH，4bH为长和宽的矩形计算域，aH和bH分别为不考虑涂层情况下接触副赫兹椭圆点接触椭圆的长半轴和短半轴，Nx为x方向的单元数目，Ny为y方向的单元数目；作用在各表面单元上的接触压力近似为均匀分布的，并以单元的节点处的法向位移作为该单元的法向位移，那么固体润滑涂层体系弹性半空间与椭球体弹性点接触模型写为：式中：p(m,n)-作用在(m,n)表面单元上的接触压力，g(m,n)-接触外载荷Qi作用后(m,n)表面单元节点处的法向间距，Ωc为两接触体接触表面发生真实接触的区域，Δx-表面单元x方向的大小，Δy-表面单元y方向的大小，g0(m,n)-接触外载荷Qi作用前(m,n)表面单元节点处的方向间距，uz(m,n)-两接触表面在(m,n)表面单元节点处的法向位移之和，-接触外载荷作用下两接触体的法向趋近量；对于弹性接触问题，根据线性叠加原理，接触表面法向位移由接触压力导致的法向位移影响因子矩阵与接触压力矩阵的离散卷积运算获得，因此两接触表面在(m,n)表面单元节点处的法向位移之和uz(m,n)为：式中：Kκ(m-r,n-s)为由表面法向压力导致的接触表面法向位移的影响因子，其物理意义为(r,s)单元上的单位均布压力在(m,n)单元节点中心处引起的法向位移；p(r,s)指(r,s)表面单元上的压力；下标符号κ的取值1和2表示Kκ(m-r,n-s)分别为椭球体和涂层体系半空间表面法向压力导致的接触表面法向位移的影响因子；固体润滑摩擦副的接触变形δi为(Nx/2,Ny/2)表面单元的位移之和uz(Nx/2,Ny/2)，所以接触变形δi为：步骤三：基于非线性模型，应用最小二乘法拟合获得固体润滑三维点接触摩擦副接触载荷与接触变形的非线性关系：Q＝f(δ) (2)。2.根据权利要求1所述的一种考虑涂层影响的点接触力载与变形关系的建立方法，其特征在于：步骤三中，所述非线性模型为多项式或幂函数。3.根据权利要求2所述的一种考虑涂层影响的点接触力载与变形关系的建立方法，其特征在于：步骤三中，基于幂函数非线性模型应用最小二乘法拟合获得三维点接触摩擦副的接触载荷与接触变形的幂函数关系式，具体如下：①确定待拟合数据点：由步骤一、步骤二获得的一系列载荷点接触载荷下的固体润滑接触副的接触变形构建一组数据点：(δi,Qi),(i＝1,2,…,n+1) (6)式中，δi为接触载荷点Qi作用下的接触变形，i为划分后接触载荷点的下标；②对幂函数形式的固体润滑摩擦副的接触变形与接触载荷的非线性关系式：进行线性化处理如下：ln Q＝ln c1+c2lnδ (8)式中：Q为接触载荷，δ为接触变形，c1和c2为待定参数；令y＝ln Q，b＝ln c1，a＝c2，x＝lnδ，式(7)转化为y＝ax+b (9-a)yi＝ln Qi (9-b)xi＝lnδi (9-c)式中：xi为接触载荷Qi取对数，yi为接触变形δi取对数，a和b为待定参数；③最小二乘拟合：根据最小二乘拟合原理，待定参数a和b的取值应使各数据点的值yi与拟合关系曲线(9-a)的预测值axi+b的偏差平方和最小，即因此待定参数a和b的取值为④接触变形关系待定参数为：c2＝a (11)c1＝eb (12)。