新的柔性结合部法向接触刚度和接触阻尼方程

以修正分形几何学理论和赫兹法向接触力学方程为基础,推导出了柔性结合部法向接触刚度与阻尼方程。假设峰元顶端的曲率半径为变量,提出了一种全新的求导函数而非偏导函数的求解方法,建立了单峰元与平面接触的法向接触刚度方程。数值模拟表明:峰元承担的法向弹性载荷与其顶端的变形量之间符合非线性幂函数凹弧关系;降低表面粗糙度或增加法向接触载荷都将增大实际接触面积;当表面粗糙轮廓分形维数在较小范围内时,实际接触面积随着表面粗糙轮廓分形维数的增加而增大,而当表面粗糙轮廓分形维数在较大范围内时,实际接触面积随着表面粗糙轮廓分形维数的增加而变小;降低表面粗糙度或增加表面粗糙轮廓分形维数与法向接触载荷皆将增加法向接触刚度;法向接触阻尼随着表面粗糙轮廓分形维数的增加先减小后增大;当表面粗糙轮廓分形维数小于临界值时,法向接触阻尼随着分形粗糙度的增大而增大,而当表面粗糙轮廓分形维数超过转折点时,法向接触阻尼随着分形粗糙度的增大而减小;当法向接触载荷增大时,法向接触阻尼略微减小。     

结合面切向接触等效黏性阻尼分形模型

基于MB接触分形理论、结合面切向接触阻尼耗能机理以及阻尼损耗因子的定义,建立了结合面切向接触等效黏性阻尼的分形模型及其损耗因子模型。所建模型表明,结合面切向接触等效黏性阻尼与结合面法向接触载荷、摩擦系数、材料塑性指数、结合面上的切向动态载荷幅值与法向接触载荷之比(简称切法向载荷比)、结合面分形维数以及分形粗糙度参数之间具有复杂的非线性关系,而结合面切向接触阻尼损耗因子与结合面分形维数和分形粗糙度参数无关,仅与切法向载荷比和摩擦系数有关。模型的仿真结果表明,结合面切向接触阻尼损耗因子随着切法向载荷比的增大而增大,随结合面摩擦系数的增大而减小;结合面切向接触等效黏性阻尼随着结合面法向接触载荷、摩擦系数、材料塑性指数的增大而增大,随着结合面分形粗糙度的增大而减小;结合面切向接触等效黏性阻尼随结合面分形维数的变化规律较为复杂,先随着分形维数的增大而增大,在分形维数值1.65附近出现最大值,而后随着分形维数的增大而减小。     

固定结合面切向接触阻尼分形模型

基于三维接触分形理论和固定结合面切向接触阻尼耗能机理,建立了考虑弹塑性接触变形机制的固定结合面切向接触阻尼分形模型。模型仿真结果表明:切向接触阻尼会随着法向接触载荷的增大而增大;当表面分形维数小于等于2.5时,切向接触阻尼会随着分形维数的增大而增大;当表面分形维数大于2.5时,切向接触阻尼会随着分形维数的增大而减小。模型仿真结果也揭示了分形粗糙度尺度参数、塑性指数和摩擦因数对切向接触阻尼的影响。     

平面结合面切向接触阻尼分形模型及其仿真

基于接触分形理论和结合面接触阻尼耗能机理,以及球体与平面接触时的阻尼耗能理论,建立了平面结合面切向接触阻尼的分形模型,通过数值仿真直观揭示了平面结合面切向接触阻尼耗能与结合面法向载荷以及结合面表面粗糙轮廓分形维数之间的非线性关系.仿真结果表明:平面结合面切向接触阻尼的耗能随着结合面法向载荷的增大而减小;当结合面表面粗糙轮廓分形维数小于等于1.2时,平面结合面切向接触阻尼耗能随结合面表面粗糙轮廓分形维数的增大而增大;当结合面表面粗糙轮廓分形维数大于1.2时,平面结合面切向接触阻尼耗能随结合面表面粗糙轮廓分形维数的增大而减小.仿真结果验证了模型的有效性.     

考虑分布域扩展因子的结合面法向接触阻尼建模

为了准确计入微凸体接触面积分布对结合面特性的影响,提出了考虑分布域扩展因子的结合面法向接触阻尼分形模型。以分形理论为基础,采用Majumda-Bhushan修正模型,引入微接触大小分布域扩展因子,分析了结合面法向总载荷与接触面积的关系。利用阻尼耗能机理,提出了将结合面法向接触特性等效为弹簧和黏性阻尼器的动力学系统,推导出结合面法向接触阻尼损耗因子,进而建立结合面法向接触阻尼模型,并进行了归一化处理。仿真结果表明:结合面归一化法向接触阻尼随着分形维数的增大而增大,随着应变指数的增大而减小;归一化法向接触载荷和归一化分形粗糙度参数对结合面归一化法向接触阻尼的影响趋势均与分形维数所处范围有关。线轨滑台上进行的试验表明:试验模态参数与理论模态分析结果一致,所提模型能够准确地描述结合面的动态特性。     

考虑硬度变化的结合面法向刚度建模

基于分形理论及M-B模型,引入微接触点域扩展因子,综合考虑微凸体弹性接触变形、弹塑性接触变形和完全塑性变形,进而考虑微凸体弹塑性变形阶段硬度随其几何形貌的改变而变化,建立对应的结合面法向接触刚度模型。通过软件仿真发现:考虑微凸体硬度随其几何形貌改变后,无量纲接触载荷较将硬度视为定值时要小,且随着分形维数的增大,二者差异在逐渐增大;考虑微凸体硬度随微凸体几何形状改变而变化后,结合面无量纲法向接触刚度相较将硬度视为定值时大;无量纲法向接触刚度随着无量纲接触载荷、分形维数和塑性指数的增大而增大,但随着无量纲特征分形粗糙度的增大而减小。     

应用改进分形几何理论的结合部切向刚度模型

针对现有分形接触理论对2个机械部件粗糙表面相互接触形成的结合部的切向接触刚度分形模型存在违反赫兹法向接触力学的缺陷,以改进分形几何理论为基础、在严格应用赫兹法向接触力学的基础上,推导出结合部总切向接触静弹性条件刚度、总条件法向载荷的分析解。数值仿真表明:结合部的切向接触静弹性刚度随着总法向载荷的增加基本上呈线性增加的态势,随着表面轮廓分形维数的增加而增大,随着分形粗糙度的减小而增大;在恒定法向载荷作用下,最初作用于结合部的切向载荷使得切向接触静弹性刚度最大,该刚度随着切向载荷的增加而减小,随着静摩擦系数的增加而增大;随着法向载荷的增加,法向接触静弹性刚度的增量加大。该结果可为进一步研究粗糙表面的分形特性提供参考。     

基于微凸体连续变形理论的结合面切向刚度模型

提出了一种包含微凸体接触过程中弹塑性过渡阶段弹性的粗糙表面切向接触刚度模型。现有粗糙表面切向接触刚度模型中只考虑微凸体完全弹性接触阶段的弹性,未考虑弹塑性过渡阶段弹性。根据微凸体变形过程中应满足的连续性条件,微凸体变形过程可以按如下划分:完全弹性阶段、弹塑性过渡阶段(包含三个不同阶段)以及完全塑性阶段,依据分形理论和Hertz接触理论,建立了计入微凸体弹塑性过渡阶段弹性的结合面切向接触刚度模型。分析与试验结果表明:考虑微凸体的弹塑性过渡阶段的切向接触刚度相比与不考虑弹塑性过渡阶段更加的符合实际接触情况。新模型的切向接触刚度随着法向载荷和较软材料的屈服强度与两材料的复合弹性模量的比值的增大而增大,还由于分形粗糙度与切向载荷的增大,该切向接触刚度会减小,并且伴随着分形维数的增大,切向接触刚度先增大后减小,即接触刚度存在一个峰值。此外,当较软材料的屈服强度与两材料的复合弹性模量的比值增大时,刚度峰值所对应的分形维数在减小。     

三维分形固定结合面法向接触刚度的研究

基于传统的M-B模型,在考虑微凸体弹塑性变形的基础上,应用更能符合结合面实际表面形貌的修正的W-M函数,建立了三维分形结合面法向刚度模型.通过建立的刚度模型研究了分形尺度参数和分形维数对法向接触刚度的影响,并对出现的一些现象进行了分析.随着法向载荷与材料特性参数σ_y/E的增加,结合面的法向刚度也增加;但随着分形尺度参数的增加,结合面的法向刚度却减小;随着分形维数的增加,当2.1≤D≤2.6时,结合面法向刚度增加,而当2.6≤D≤2.9时,结合面的法向刚度却减小;在考虑弹塑性变形的情况下,三维分形结合面的法向刚度要小于二维分形结合面的法向刚度.     

结合面法向动态参数的分形模型

为了对结合面法向动态参数进行正确的理论计算,以分形接触理论为基础,建立了结合面法向动态参数的理论分形模型,揭示了接触刚度和接触阻尼与作用在粗糙表面上的法向载荷、粗糙表面材料性能常数,以及分形参数等因素之间的复杂关系,并对该模型进行了数值仿真.仿真结果表明:结合面的实际接触面积仅占名义接触面积的一小部分,降低表面粗糙度或增加法向载荷都将增大实际接触面积;接触刚度和接触阻尼与分形参数之间具有较强的非线性关系,而法向栽荷对接触刚度的影响较为明显,当栽荷增加时,刚度值也随之增大,但对接触阻尼的影响可以忽略.     

MB模型的修正及应用

由于MB模型存在3个缺陷,工程上会引入无法接受的误差.本文修正了MB模型的后两个缺陷.针对第2个缺陷,提出了实际临界接触面积、最小有效分形维数、量纲一的最大有效分形粗糙度参数、最小有效材料特性等术语.根据MB模型和提出的计算术语,对两弹塑性接触粗糙表面的接触行为及界面的静摩擦因数进行了研究.研究结果表明：界面的静摩擦因数首先随分形维数的增加而增加,然后随分形维数的增加而减小;界面的静摩擦因数随分形粗糙度参数的增加而减小,但随材料特性的增加而增加,也随总法向载荷的增加而增加;当分形维数较小或分形粗糙度参数较大或材料特性较小时,静摩擦因数-量纲一的总法向载荷曲线为凸弧.     

两弹塑性非赫兹接触粗糙表面温升的分形模型

针对赫兹接触理论存在的3个缺陷,考虑表面粗糙度和塑性变形,适当处理接触物体交界面处的摩擦,将赫兹接触理论以更符合实际的方式推广到滑动接触．采用球形微凸体的赫兹接触理论和MB修正模型,对微接触点的温升进行了分析,得到了低速滑动区域内的分形区域实际接触面积温升的补充累积概率分布函数的封闭形式表达式．分析结果表明：分形区域的最大温升随滑动速度增大而线性增大,非零域随滑动速度增大而扩展．对于固定的量纲一分形粗糙度参数,最大温升随分形维数增大而减小;对于固定的分形维数,最大温升随量纲一分形粗糙度参数增大而增大．温升的补充累积概率分布函数随滑动速度增大而增大,随分形维数增大或量纲一分形粗糙度参数减小而减小．平均温升为最大温升的0．4023倍,温升的标准差为最大温升的0．24倍．     

考虑微凸体弹塑性过渡变形机制的结合面法向接触刚度分形模型

摘要： 基于粗糙表面微凸体变形的连续性和光滑性原理,研究了在法向载荷逐渐增加时的粗糙表面单个微凸体弹塑性过渡变形机制,提出了考虑弹塑性过渡变形机制的结合面微凸体微观接触模型,建立了法向接触载荷和法向接触刚度的数学模型;基于分形几何理论,建立了结合面法向接触刚度的分形模型,并对结合面法向接触刚度进行仿真计算.结果表明：在较小的塑性指数条件下,法向接触载荷与法向接触刚度近似呈线性关系;在较大的塑性指数条件下,法向接触载荷与法向接触刚度呈非线性关系;法向接触刚度随着分形维数和法向接触载荷的增加而增大,随着无量纲分形特征尺度系数的增大而减小;所得结合面法向接触刚度的仿真计算值与铣削加工和磨削加工条件下的实验值较吻合.     

虑隧道表面特性的硬岩全断面掘进装备撑靴接触界面刚度分析

采用修正的各向同性三维分形表面的W M函数构造隧道掘进机掘进后的隧道粗糙表面;考虑法向支撑载荷作用和岩石在压缩载荷下的失效机制,建立了掘进机撑靴与岩石粗糙表面法向接触刚度模型.研究变载荷工况下,不同粗糙隧道表面特性与掘进机撑靴接触界面刚度特性变化规律.结果表明：岩石的失效对接触刚度特性影响明显,相同的载荷下,岩石的失效会带来接触刚度的减小;当岩石弹性模量相同时,接触刚度随着硬度的增大而增大;当岩石硬度相同时,接触刚度随着弹性模量的增加而减小;随着外部载荷的增加,粗糙表面的接触刚度随之增加,而表面粗糙度的增大会引起界面接触刚度相应减小.     

考虑三维结合部形貌的静摩擦因数非线性分形模型

考虑三维结合部形貌的W-M函数,推导了结合部静摩擦因数非线性分形理论模型.数值模拟了考虑三维形貌的结合部静摩擦因数与法向载荷P、分形维数D、分形尺度系数G的关系,以及在二维分形和三维分形模型中的差异.结果表明:结合部静摩擦因数与法向载荷成单调递增关系,与分形尺度系数成单调递减关系.当D小于2.5时,静摩擦因数随分形维数的增大而增大;当D大于2.5时,静摩擦因数随分形维数的增大而减小;三维分形下的静摩擦因数小于二维分形下的静摩擦因数.     

结合面静摩擦因数三维分形模型

摩擦磨损导致部分能源的消耗,且主要发生在接触区域,其接触特性依赖于接触表面之间的微观弹塑性变形.本文在三维接触分形理论的基础上,考虑微凸体的完全弹性和完全塑性阶段的变形机制,建立了结合面静摩擦因数的三维分形模型,推导了静摩擦因数f的解析解.通过仿真研究了分形维数D、无量纲分形特征尺度参数G＊以及无量纲法向总载荷P＊等因素对静摩擦因数的影响规律.结果显示,结合面静摩擦因数f先随着分形维数D的增大而增大,然后随着分形维数D的增大而减小;随着无量纲法向总载荷P＊的增大而增大;随着G＊的增加而减小.     

两弹性接触粗糙快速滑动表面温升的分形模型

采用球形微凸体的赫兹接触理论和MB模型,对微接触点的温升进行了分析,得到了快速滑动区域内的分形区域实际接触面积的温升概率分布密度、温升补充累积概率分布函数的封闭式表达式.分析结果说明:量纲一特定滑动区域的实际接触面积随量纲一最大Jaeger参数增加而单调减小.量纲一最大温升随分形维数增加而减小,但随分形粗糙度参数增加而增加.量纲一温升随分形维数增加而增加.当分形维数为1.5时,实际接触面积的温升概率分布密度等于在一个弹性微接触点面积上的温升概率分布密度,基于正八边形面积的近似解适当接近精确解.温升的补充累积概率分布函数随分形维数、滑动速度和量纲一分形粗糙度参数增加而增加.     

基于分形接触理论的结合面法向接触参数预估

基于分形接触理论,建立了结合面法向接触参数的分形预估模型,通过粗糙表面材料性能参数、法向载荷及粗糙表面的分形参数来预估法向接触刚度和接触阻尼,并对其变化规律进行数值仿真.结果表明：结合面之间的接触处于弹性变形和与塑性变形共存的状态,且小接触点面积的微凸体发生塑性变形,而大接触点面积的微凸体发生弹性变形;法向接触参数与分形参数之间具有较强的非线性关系;同时,法向接触刚度随法向载荷的增大而逐渐增加,但法向载荷对结合面的法向阻尼特性影响较小;仿真结果中极值点的存在,为结合面接触参数的优化设计提供了依据.     

考虑弹塑性变形机制的结合面静摩擦因数模型

在三维接触分形理论的基础上,从理论上建立了一种考虑微凸体的完全弹性、弹-塑性和完全塑性三种变形机制的结合面静摩擦因数三维分形模型。该模型反映了结合面法向载荷、分形维数以及分形特征尺度参数对结合面静摩擦因数的影响。仿真结果表明,静摩擦因数f随着无量纲法向接触总载荷P*的增大而增大;分形维数D的增大,先增大后减小;无量纲分形特征尺度参数G*的增大而减小;考虑弹-塑性变形机制的结合面静摩擦因数小于不考虑弹-塑性变形机制的结合面静摩擦因数。