[发明专利]一种考虑三维分形的法向界面刚度预测方法

摘要

本发明涉及一种考虑三维分形的法向界面刚度预测方法，步骤为将描述二维分形曲线函数改进为模拟三维分形形貌的修正函数，将此函数描述的波峰与波谷幅值差表示为接触变形量；将两粗糙微凸体间的接触等效为一刚性平面和一等效微凸体间的接触，则等效微凸体和刚性平面间的实际接触面积；分别计算弹性变形阶段、弹塑性变形阶段的变形量；面积分布函数与实际接触面积计算；单个微凸体刚度和总界面刚度计算；弹性阶段刚度kn1和弹塑性阶段刚度kn2和总界面刚度的计算。本发明为精密机械界面间刚度的获得提供了一种简单易操作的预测方法，考虑了微凸体的弹塑性变形、接触间摩擦因数及三维分形分布的影响，得到的结果可为预测、控制界面动态特性提供技术参考。

主权项

一种考虑三维分形的法向界面刚度预测方法，其特征在于包括以下步骤：1)三维表面形貌模拟：将描述二维分形曲线函数改进为模拟三维分形形貌的修正函数，将此函数描述的波峰与波谷幅值差表示为接触变形量δ＝2GD‑2(lnγ)0.5(2r′)3‑D，其中，D为三维表面形貌分形维数，2<D<3，G为表面形貌的分形粗糙度，γ为频率密度参数，r′为微凸体截断半径；2)微凸体接触等效处理：将两粗糙微凸体间的接触等效为一刚性平面和一等效微凸体间的接触，则等效微凸体和刚性平面间的实际接触面积a＝πRδ；其中，δ为微凸体的接触变形量，R为微凸体的曲率半径；3)微凸体变形包括三个阶段，即弹性变形阶段、弹塑性变形阶段以及塑性变形阶段；分别计算弹性变形阶段、弹塑性变形阶段的变形量；4)面积分布函数与实际接触面积计算：整个接触界面上的面积分布函数为整个接触界面的实际接触面积为其中，D为三维表面形貌分形维数，2<D<3，al表示所有微凸体接触中最大的接触面积；a为实际接触面积；5)单个微凸体刚度和总界面刚度计算：单个微凸体刚度包含其弹性变形和弹塑性变形两个阶段，即弹性阶段刚度kn1和弹塑性阶段刚度kn2，kn1=dLedδ=4(2-0.5D)Eπ-0.5a0.53(1.5-0.5D);]]>kn2=dLepdδ=2.79λσyae-n(n+1)an-0.5+0.5D25.5-1.5DGD-2(lnγ)0.5π0.5D-1.5(1.5-0.5D),]]>总界面刚度为：Kn=2.79λσyae-n(n+1)(D-1)a10.5D-0.526.5-1.5DGD-2(lnγ)0.5π0.5D-1.5(1.5-0.5D)n(aen-apn)+2(2-0.5D)(D-1)3(1.5-0.5D)(1-0.5D)Eπ-0.5a10.5D-0.5(a11-0.5D-ae1-0.5D)]]>其中，Le为单个微凸体弹性阶段的载荷，ae为单个微凸体的临界弹性变形面积，al表示所有微凸体接触中最大的接触面积，δ为微凸体的接触变形量，D为三维表面形貌分形维数，2<D<3，Lep为单个微凸体弹塑性阶段的载荷，G为表面形貌的分形粗糙度，ap为临界塑性变形面积，σy表示相互接触材料中较软的屈服强度，λ为定义的系数，a为实际接触面积，n为材料硬度指数。