机械结合面切向接触阻尼计算模型

针对两粗糙表面在法向力和切向力共同作用下相互接触时结合面切向阻尼的问题进行了研究。首先,根据KE模型对单个微凸体在弹性、弹塑性、塑性变形阶段的切向接触行为进行了分析,获得了微凸体在三个变形阶段的黏滑特性;然后,基于GW统计模型建立了一种在微凸体法向弹性、弹塑性和塑性变形机制基础上,考虑微凸体黏滑摩擦行为的机械结合面切向接触阻尼统计模型;最后,分别讨论了机械结合面的法向预载荷、切向激振频率和切向动态位移幅值对机械结合面切向阻尼的影响。研究表明:结合面切向接触阻尼系数随着结合面法向载荷的增大而增大,随着切向激振频率和切向动态位移幅值的增大而减小;在高频率、大幅值下,结合面切向接触阻尼系数几乎与动态位移幅值和激振频率无关。为了验证模型的准确性,构建了动态切向力作用下的结合面切向阻尼试验,其试验结果与理论仿真变化规律与量级基本一致,从而证明了本文所提出的切向阻尼模型的有效性。      

固-液结合面法向动态接触刚度及阻尼的理论模型与实验分析

固-液接触状态广泛存在于机床核心单元关键零部件的接触运动副中,精确获得固-液结合面法向接触刚度及阻尼参数是高档数控机床产品在研发阶段就存在的一个关键理论与技术问题,并且仍然尚未根本解决。固-液结合面在介观层面上表现为两个粗糙表面的接触,在微观层面上表现为微凸体之间的接触,并在中/重载荷作用下微凸体可能会发生弹性/弹塑性/塑性变形。为了揭示静动态外载荷对固-液结合面接触刚度及阻尼的影响,分别基于GW模型、KKE模型和AF模型对接触微凸体弹性/弹塑性/塑性变形展开研究,并结合流体动力润滑REYNOLDS方程,建立了考虑接触微凸体弹性/弹塑性/塑性变形的固-液结合面接触刚度及阻尼模型。并对其进行实验验证,结果表明：随着预载荷的增大固-液结合面法向动态接触刚度表现出先减小后增大的规律,当接触载荷小于某阈值时动态接触刚度较大,反之静态接触刚度较大;法向动态接触刚度随着法向相对位移幅值的增大而增大,在低载荷时呈线性规律,而高载荷呈非线性规律;法向动态接触刚度随激振频率增大呈线性增大,且载荷越大线性斜率越小。对于法向接触阻尼,随着法向相对位移幅值和接触载荷增大呈非线性增大,随着激振频率增大几乎不变。精确获得固-液结合面法向接触刚度和阻尼及其关键因素的影响规律,对机械系统的分析、设计、优化以及静、动态性能控制都具有重要的理论意义。     

弹塑性微凸体侧向接触相互作用能耗

传统的结合面研究多基于光滑刚性平面与等效粗糙表面接触假设,忽略了结合面上微凸体侧向接触及相邻微凸体之间的相互作用,这导致理论模型与实际结合面存在较大出入.针对承受法向静、动态力的机械结合面,从微观上研究了微凸体侧向接触及相互作用的接触能耗.将法向静、动态力分解为法向分力和切向分力,获取弹性/弹塑性/塑性阶段考虑微凸体侧接触及相互作用的加、卸载法向分力-变形和切向分力-位移的关系.通过力的合成定理,从而获取加、卸载法向合力与总变形之间的关系,由于法向分力产生的塑性变形及切向分力产生的摩擦,导致加载、卸载法向合力-总变形曲线存在迟滞回线.通过对一个加、卸载周期内的法向合力-总变形曲线积分,获得一个周期的微凸体接触能耗,包括应变能耗及摩擦能耗.仿真分析表明:微凸体在3个阶段的能耗均随变形的增大而非线性增大.微凸体侧向接触角度越大,能耗越大,且在弹性阶段最为明显.在弹性阶段,仅存在侧向的摩擦能耗,故结合面在低载荷作用下必须采用双粗糙表面假设.在塑性阶段,由于微凸体接触能耗为应变能耗,且接触角对其能耗影响甚微,故结合面在大载荷作用下可采用单平面假设对其进行研究.相对于KE和Etsion模型,本文提出的模型与Bartier的实验结果更吻合.     

考虑粗糙结合面弹塑性接触的黏滑摩擦建模

提出一种同时考虑粗糙面上微凸体弹性变形和塑性接触的切向黏滑摩擦建模方法。采用Hertz弹性理论和Mindlin解描述弹性接触微凸体的切向载荷和相对变形的关系;采用AF(Abbott-Firstone)塑性理论和Fujimoto模型描述塑性接触微凸体切向载荷和相对变形的关系。再利用GW(Greenwood-Williamson)模型统计分析方法建立粗糙表面切向载荷和相对变形之间的关系。将模型与仅考虑微凸体弹性接触情况的模型进行对比,并研究了不同塑性指数对切向载荷和相对变形关系的影响。结果表明：与完全弹性接触模型相比,本文模型引入了塑性接触理论,能够更好地描述粗糙表面切向载荷和相对变形关系,并且考虑不同接触条件下弹性变形微凸体和塑性变形微凸体对切向接触载荷的贡献,在微滑移阶段,主要由弹性接触变形影响,而在进入宏观滑移阶段之后,切向行为主要由塑性变形影响。界面切向载荷由黏着和滑移接触作用共同决定,随着切向变形的增加,滑移接触力逐渐增加,而黏着接触力先增加后减少,反映了界面由微滑移逐渐向宏滑移演化的过程。随着塑性指数的增加,粗糙面上发生塑性接触的微凸体数目逐渐增加,切向黏滑行为主要受到塑性接触特征的控制。     

基于微凸体接触压力分布的粗糙表面法向接触建模

提出一种考虑微凸体弹塑性接触变形影响的粗糙表面法向接触力学模型。采用有限元模拟微凸体弹塑性接触过程,分析不同塑性屈服条件对微凸体接触载荷和实际接触面积的影响。再根据微凸体接触面上压力分布的变化规律,将微凸体的接触状态分为完全弹性接触阶段、弹塑性接触阶段、完全塑性接触阶段。分析接触面压力变化规律对微凸体法向接触载荷-变形的影响,再利用GW模型的数理统计分析的方法得到粗糙表面的法向接触载荷。将本文提出的模型与完全弹性模型、CEB模型、ZMC模型、KE模型、JG模型进行对比,并且研究了塑性指数对粗糙面接触载荷-平均高度距离的影响。结果表明,本文提出的模型能够更好地描述微凸体法向接触载荷与接触变形的变化趋势,模型预测粗糙表面法向载荷与ZMC、KE模型具有较好的一致性;粗糙面接触载荷随着平均接触距离增加而减少,随着塑性指数的增加,不同模型预测的法向接触载荷差异逐渐增大。     

依据各向异性分形理论的固定结合面椭圆弹塑性法向接触刚度建模及仿真分析

依据各向异性分形几何理论,将结合面微凸体拓展为椭球体,结合微凸体接触点面积大小分布函数以及概率论相关理论获得关于椭圆形接触点接触面积以及离心率的二维联合分布密度函数,应用椭球体的赫兹接触理论,进而建立了依据各向异性分形理论的结合面椭圆弹塑性法向接触刚度模型,并采用MATLAB软件对影响结合面法向刚度的相关因素进行了数值仿真及结果分析。结果表明：结合面椭圆形接触点离心率分布情形对结合面总刚度具有明显影响,结合面总刚度随形状参数 的增大而增大,却随着形状参数 的增大而减小;结合面法向刚度随法向载荷的增大而增大;在同一载荷作用下,结合面法向接触刚度随塑性指数增大而增大,随分形粗糙度的增大而减小,但是随分形维数的增大而先增后减。该模型对模型优化进而提高计算精度提供了一定的理论依据。     

基于三维各向异性分形理论的固定机械结合面法向接触刚度建模研究

摘要：将结合面微凸体拓展为椭球体,基于KE有限元模型,类比球形微凸体在弹塑性接触变形阶段法向载荷、接触面积以及变形量之间的关系,根据椭球体的弹性接触理论,采用代入法得到了表征椭球形微凸体弹塑性接触变形机制的对应关系式。结合三维各向异性分形几何理论得到了结合面接触点离心率的有效区间为[0,0.7374],在此基础上,假设结合面接触点离心率分布在有效区间内,在该区间上服从指数分布,且与接触点面积分布相互独立,根据概率论以及接触点的面积大小分布函数,获得了关于结合面接触点面积与离心率的二维联合分布密度函数,进而建立了包含椭球形微凸体完全弹性、弹塑性以及完全塑性三种变形机制的结合面法向接触刚度分形模型。所建模型理论刚度与实验数据的对比结果,表明了模型的正确性及有效性,能较好的预测固定结合面在轻载状态下的法向接触刚度。     

介观双粗糙弹塑性流体动力润滑界面法向接触刚度模型

弹性流体动力润滑状态通常出现在机械高副零部件的点/线接触部位,如齿轮、轴承和蜗轮蜗杆等. 宏观上点/线接触在介观层面表现为两粗糙表面的接触,在微观层面上则又表现为微凸体间的接触. 由于在中/重载荷作用下,粗糙表面上的微凸体发生接触后会产生弹塑性/塑性变形,从而使得两粗糙表面的弹流润滑接触转变为弹塑性流体动力润滑接触. 此外,界面的接触刚度决定了机械装备的整机刚度. 为了精确获得弹性流体动力润滑状态下界面法向接触刚度及其主要影响因素,基于界面的法向接触刚度由固体接触刚度和润滑油膜刚度两部分构成的思想,根据固体弹塑性理论和流体动力学理论,分别对界面间微凸体侧接触及部分膜流体动力润滑进行分析,从微观入手揭示双粗糙表面弹塑性流体动力润滑接触机理,进而建立考虑微凸体侧接触弹塑性变形的流体动力润滑界面法向接触刚度模型. 通过仿真分析,揭示了法向载荷、卷吸速度、表面粗糙度及润滑介质特性等因素对润滑界面法向接触刚度的影响规律. 研究表明:在相同速度、粗糙度及润滑油黏度的工况下,固体接触刚度和油膜接触刚度均随着法向接触载荷的增加呈非线性增大;在相同载荷、速度及润滑油黏度的工况下,接触表面粗糙度越大,表面形貌对于润滑状态的影响较强,固体接触刚度占界面总刚度的主要部分,界面主要由固体承载;在相同载荷、粗糙度及润滑油黏度工况下,随着卷吸速度的增大,固体接触刚度逐渐减小,油膜刚度占界面总刚度的主要部分;在相同载荷、粗糙度及速度工况下,随着润滑油黏度的增大,油膜刚度基本保持不变,固体接触刚度基本不受润滑油黏度的影响. 通过理论建模准确获得单位面积弹塑性流体动力润滑结合面法向接触刚度,对改善机械装备动态性能、提高机械装备的可靠性具有重要的理论和实际意义.      

考虑摩擦系数和微凸体相互作用的粗糙表面接触热导分形模型

基于三维分形理论，建立了考虑摩擦系数和微凸体相互作用的粗糙表面接触热导分形模型，并且考虑了微凸体的弹性变形、弹塑性变形和完全塑性变形. 通过该模型，分析了摩擦系数、分形维数、分形粗糙度和接触载荷对热接触热导的影响. 研究结果表明:接触热导随着摩擦系数和分形粗糙度的增大而减小，随着分形维数和接触载荷的增大而增大. 该研究为开展接合面的热传递提供了一定的理论基础.      

包含几何误差的机械结合面法向刚度研究

对包含几何误差的机械结合面进行离散化,离散后的微表面的基准平面高度满足结合面几何误差分布.每个微表面内,微凸体的高度只受粗糙度的影响.基于接触理论建立了微表面的法向刚度模型,通过对微表面模型集成获得了结合面的法向刚度模型.通过对所建模型的数值仿真,揭示了结合面法向刚度与间隙的非线性关系,几何误差的幅值和波长对法向刚度的影响以及非线性刚度对结合面振动特性的影响.计算结果表明:法向刚度随着间隙的减少而迅速增加,几何误差会导致结合面宏观上的局部接触和应力集中;在相同干涉量下,法向刚度随着几何误差幅值的增加而增加,但与结合面的波长没有关系;非线性刚度会导致结合面固有频率的下降和振动位移的不对称.     

分形粗糙表面弹塑性接触力学模型

为了建立更为精准的粗糙表面接触模型,基于分形理论提出一种修正的弹塑性接触模型,该模型同时考虑了微凸体的弹性、弹塑性、完全塑性三种变形状态;建立了临界第一弹塑性接触面积、临界第二弹塑性接触面积与临界弹性接触面积之间的关系,分别推导出三种变形阶段的接触载荷与接触面积之间的关系式;结合海洋岛屿面积分布密度函数,获得了接触表面上总接触载荷与真实接触面积之间的表达式。计算结果表明:单个微凸体的临界弹性接触面积不但与材料属性和分形参数有关,也与微凸体的尺寸有关;分形参数一定时,随着接触载荷的增大,真实接触面积也在增大,弹性接触面积与真实接触面积的比值逐渐减小。该模型的建立为进一步研究粗糙表面的摩擦、磨损、润滑提供了理论依据。     

基于分形理论的双渐开线齿轮接触应力研究

综合考虑接触面粗糙度、材料特性等因素对齿轮接触应力的影响，基于分形理论和经典Hertz接触理论建立双渐开线齿轮分形接触模型. 该模型中，影响载荷和实际接触面积的主要因素包括分形维数、粗糙度幅值和材料特性参数. 理论分析表明:分形维数一定时，真实接触面积随着载荷的增大而增大；载荷一定时，接触面积随着粗糙度幅值的增大而减小；随着材料特性参数值的增加，在一定程度上加强了软材料轮齿承载能力，同时会使得微凸体由弹性变形到塑性变形的临界面积减小. 对比分形接触模型和有限元模型两种计算双渐开线齿轮轮齿接触应力方法，结果证明了分形接触模型计算双渐开线齿轮接触应力的有效性.      

基于Iwan模型的接合面切向响应建模

针对工程中常见预紧力作用下的搭接接头,研究其在小幅切向位移激励时的切向位移响应问题,为此提出一种新的基于实际表面形貌和材料性能参数的滑移力密度分布函数.应用该分布函数得到搭接接头切向响应本构模型,并获得单位加载周期内的迟滞曲线和能量耗散值, 通过与已出版的实验结果相对比,发现得到的模拟值与实验结果吻合, 证明该模型的合理性.在此基础上利用该分布函数研究了接合面切向位移与切向力、切向接触刚度及能量耗散之间的关系,结果表明: 建立的模型能很好地描述接合面间切向力与切向位移之间的关系,临界滑移力函数开始迅速上升, 到达最大值后迅速收敛到零;切线力与切向位移之间表现出非线性特性, 随着切向位移的增大,切向接触刚度表现出"软化"现象;初始切向刚度与法向载荷、粗糙度参数及塑性指数有关, 对于确定的接触表面,法向力越大, 初始切向刚度越大; 初始切向刚度同样也随着塑性指数的增大而增大.      

基于Iwan模型的接合面切向响应建模

针对工程中常见预紧力作用下的搭接接头,研究其在小幅切向位移激励时的切向位移响应问题,为此提出一种新的基于实际表面形貌和材料性能参数的滑移力密度分布函数.应用该分布函数得到搭接接头切向响应本构模型,并获得单位加载周期内的迟滞曲线和能量耗散值,通过与已出版的实验结果相对比,发现得到的模拟值与实验结果吻合,证明该模型的合理性.在此基础上利用该分布函数研究了接合面切向位移与切向力、切向接触刚度及能量耗散之间的关系,结果表明:建立的模型能很好地描述接合面间切向力与切向位移之间的关系,临界滑移力函数开始迅速上升,到达最大值后迅速收敛到零;切线力与切向位移之间表现出非线性特性,随着切向位移的增大,切向接触刚度表现出"软化"现象;初始切向刚度与法向载荷、粗糙度参数及塑性指数有关,对于确定的接触表面,法向力越大,初始切向刚度越大;初始切向刚度同样也随着塑性指数的增大而增大.     

接触式机械密封混合润滑状态摩擦演化规律研究

接触式机械密封在运转中主要处于混合润滑状态，为探究其混合润滑状态下摩擦机理，结合粗糙面弹塑性模型，求解考虑密封端面粗糙度效应的雷诺方程，探究了转速和介质压力等工况条件对密封混合润滑状态密封摩擦参数的影响，推导了密封声发射波能量公式，将密封混合润滑状态分为磨损期和稳定期，应用1.5维谱理论提取密封特征频率，探究混合润滑状态摩擦演化规律.研究结果表明摩擦参数影响密封声发射波能量幅值且密封摩擦形式随混合润滑状态发生变化：在磨损期，微凸体接触特征频率幅值较大，密封端面间摩擦以微凸体接触为主；在稳定期，流体膜黏性剪切摩擦特征频率幅值变大，黏性剪切效应增强，密封端面只存在局部微凸体接触.所得结论对接触式机械密封混合润滑状态摩擦机理的研究具有一定的理论指导意义.     

粗糙表面的弹塑性接触研究

建立了综合载荷作用下粗糙表面弹塑性接触的确定性模型,考虑了微凸峰接触的弹塑性变形阶段,数值求解得到实际接触面积、压力分布和微凸峰塑性形变.分析了实际接触面积与法向载荷的关系,并研究了点接触的椭圆参数对上述关系的影响.建立了结点增长模型,分析了结点增长与滑动摩擦系数的关系以及滑动摩擦系数随椭圆参数的变化.结果表明:随着法向载荷增大,实际接触面积与法向载荷之间的非线性关系愈加显著;椭圆参数越大,实际接触面积越小,选择较小的椭圆参数可降低平均接触压力;结点增长的速率随滑动摩擦系数增大而增大;表面剪切作用力使最大Mises应力值升高,疲劳裂纹的发生源向表面靠近;重载时选择较小的滚动轴承沟曲率半径系数有利于减小摩擦功耗.     

计及摩擦的球面弹塑性分形模型与应用

为了研究球面接触微凸体的弹塑性变化规律,在M-B分形模型的基础上考虑基底直径的实际变化情况,引入摩擦力修正因子,并基于修正的球面接触基底直径密度分布函数,最终建立了考虑摩擦因素的球面弹塑性分形模型;并以两球体接触为研究对象,分析了模型中的相关变量对接触强度的影响规律。仿真分析结果表明:当载荷不变时,分形维数的增加会导致真实接触面积先保持不变后增大再减小;当分形维数在1.35～1.7区间时,接触面积上升的速度最快;增大摩擦系数或降低材料特性参数,有利于提高接触承载能力;粗糙度幅值参数对真实接触面积和弹性接触面积占比的影响大略一致;当粗糙度幅值增加时,真实接触面积和弹性接触面积占比都随之降低。     

基于分形理论的齿面法向接触阻尼计算模型

为了建立更为准确的微观尺度下齿面法向阻尼分形模型,以M-B模型为基础,利用存在摩擦时的弹塑性变形临界接触面积计算理论,并结合前期获得的两圆柱体接触面积分布公式,推导了考虑摩擦因素的齿面法向接触阻尼分形模型。通过MATLAB仿真,获得了模型中主要参数对分形接触模型影响的预测分析。结果表明:齿面法向接触阻尼在齿轮啮合过程中呈周期性变化,且随着法向载荷的增加而上升;摩擦系数与法向接触阻尼成正比关系;随着分形维数的增加,法向接触阻尼先下降后上升再下降;法向接触阻尼随着材料特性参数的变大而减小,随着粗糙度幅值参数的增加而下降;齿轮从初始啮合到结束啮合的过程中,法向接触阻尼由大变小;内接触的法向接触阻尼大于外接触。该模型的研究为后续进行齿轮动力学分析提供一定的理论参考。     

橡胶弹性支座动态承载特性及其影响因素分析

基于复刚度理论,采用正交实验法,对风力发电机橡胶弹性支座进行动态性能实验,研究载荷频率、载荷幅值、预载荷诸因素对弹性支座动态性能的影响规律. 实验结果表明,橡胶材料的动态性能与载荷频率和幅值具有显著的相关性. 频率增加,橡胶弹性元件的弹性刚度和阻尼因子明显增大;幅值增加,弹性刚度减小,阻尼因子却逐渐增大. 对实验数据进行方差分析,表明载荷频率对动态性能的影响最大,其次是载荷幅值,预加载因素影响非常有限.     

粗糙表面接触力学问题的重新分析

为了克服基于统计学参数的接触模型的尺度依赖性以及现有接触分形模型推导过程中初始轮廓表征受控于接触面积或取样长度的不足,基于粗糙表面轮廓分形维数D、尺度系数G和最大微凸体轮廓基底尺寸l,建立了新的粗糙表面接触分形模型,探讨了微凸体变形机制、粗糙表面的真实接触面积和接触载荷的关系,揭示了接触界面的孔隙率和真实接触面积随端面形貌、表面接触压力等参数变化的规律,给出了不同形貌界面被压实的最大变形量.结果表明:微凸体变形从弹性变形开始,并随着平均接触压力pm的增大逐步向弹塑性变形和完全塑性变形转变;接触界面的初始孔隙率φ_0随D的增大而增大,压实孔隙所需要的最大变形量δ也随之增大;接触压力p_c增大,孔隙率φ减小,并随着D的增大和G减小,φ快速减小,直至填实,变为零;D较小时,G的增大对真实接触面积的增大影响较小;D较大时,G的增大对真实接触面积的增大作用明显.研究成果为端面摩擦副的润滑与密封设计提供了理论基础.