固定接触界面切向静弹性刚度问题研究

根据两球体单峰同时受法向、切向载荷时微滑切向应力的分布以及MB模型,给出了界面的总切向接触静弹性条件刚度、总条件法向载荷的解析解。将切向接触静弹性条件刚度的解析解嵌入到有限元软件中,获得了整机的理论模态。以一款八四七厂华中工学院XHK5140型自动换刀计算机数控立式镗铣床上的结合部为研究对象,通过实验对解析解进行了定量验证。研究结果表明：在理论振型与实验振型一致的条件下,界面模型的相对误差在-19．2％～16．8％之间。     

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

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

基于泛形理论的粗糙表面法向接触刚度研究

基于泛形理论和赫兹接触理论,通过泛形海岛分布描述粗糙表面的形貌从而建立结合面弹性接触模型,求解结合面的法向接触刚度。假设粗糙接触表面微凸体的高度满足高斯分布,通过赫兹接触理论建立单个微凸体的微观接触模型,利用粗糙表面的泛形复杂度D与面积度量尺码的最小下确界as确定表面形貌。泛形复杂度反映微凸体在粗糙表面上占据的空间大小程度;度量尺码下确界是接触过程中的最小接触面积。通过泛形复杂度和面积最小下确界推导出粗糙表面法向接触刚度的解析表达式。数值算例结果表明：在相同的面积尺度区间内,粗糙表面的法向接触刚度随着泛形复杂度的增加而增加,反之减小。当接触表面的泛形复杂度不变时,粗糙表面的法向接触刚度随最小下确界的减小而增大。泛形复杂度的较大时,最小下确界的变化对于接触刚度的影响更加明显。通过与已发表的文献结果对比分析之后,发现采用泛形海岛模型所得到的结合面接触刚度与文献中实验结果吻合较好。     

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

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

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

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

高速铁道车辆轮对振型函数及约束模态分析

本文将高速车辆轮对按弹性梁简化，解析推导了弹性梁在有弹性约束条件下的振型函数及频率方程．解释了弹性梁的低阶弹性振动频率和振型与边界约束刚度有关，同时引用弹性梁频率方程的数值解及轮对模态试验结果予以证实．     

Mindlin 矩形微板的热弹性阻尼解析解

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.      

ANLYTICAL SOLUTION OF THERMOELASTIC DAMPING IN RECTANGULAR MINDLIN MICRO PLATES 1)

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.     

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

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

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

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

Contact vibration analysis of the functionally graded material coated half-space under a rigid spherical punch

This paper studies the contact vibration problem of an elastic half-space coated with functionally graded materials (FGMs) subject to a rigid spherical punch. A static force superimposing a dynamic time-harmonic force acts on the rigid spherical punch. Firstly, we give the static contact problem of FGMs by a least-square fitting approach. Next, the dynamic contact pressure is solved by employing the perturbation method. Lastly, the dynamic contact stiffness with different dynamic contact displacement conditions is derived for the FGM coated half-space. The effects of the gradient index, coating thickness, internal friction, and punch radius on the dynamic contact stiffness factor are discussed in detail.     

A finite element parallel algorithm of the parametric variational principle for elastic contact problems

Basetl on the finite element solution of the parametric varialional principle of elastic con/del problem, a corresponding parallel algorithm has been created bv utilizing the specialities of parallel computer and the architecture of concurrent processing in this paper. In this algorithm. the parallelisms have heen realized in the processes of creation and assembly of stiffness matrix, of the static condensation, of the solution of stresses and in many other aspects. The programme of this algorithm has been realized on ELXSI-6400 parallel computer of Xi'an Jiaotong University. The results of computation show that the computational time can be saved efficiently and it is an effective parallel algorithm for the analyses of contact problems.     

DYNAMIC CONTACT STIFFNESS OF VIBRATING RIGID SPHERE CONTACTING SEMI-INFINITE TRANSVERSELY ISOTROPIC VISCOELASTIC SOLID

Dynamic contact stiffness at the interface between a vibrating rigid sphere and a semi-infinite transversely isotropic viscoelastic solid is investigated. An oscillating force superimposed onto a static compressive force in the vertical direction excites the vibration of a rigid sphere, which causes variable contact radius and contact pressure distribution in the contact region. The assumption of a sufficiently small oscillating force yields a dynamic contact-pressure distribution of a constant contact radius, which gives dynamic contact stiffness at the interface between the rigid sphere and the semi-infinite solid. Numerical calculations show the influence of vibration frequency of the sphere, and elastic constants of the transversely isotropic solid on dynamic contact stiffness, which benefits quantitative evaluation of elastic constants and orientation of single hexagonal grains by resonance-frequency shifts of the oscillator in resonance ultrasound microscopy.     

Ultrasonic assessment of rough surface contact between solids from elastoplastic loading-unloading hysteresis cycle

An ultrasonic method to characterize the elastoplastic contact between two rough surfaces is presented. Ultrasonic experiments are performed on three different interfaces formed by aluminum surfaces with different levels of roughness. The frequency-dependent ultrasonic reflection coefficient from the interface is measured during loading and unloading cycles as a function of pressure, from which the ultrasonic interfacial contact stiffness is reconstructed by the least-squares inversion procedure. It is shown that one should distinguish between the ultrasonic (dynamic) interfacial stiffness and static interfacial stiffness for rough surfaces in elastoplastic contact (they are identical for purely elastic contact). It is shown that ultrasonic stiffness is associated with local unloading stiffness. An elastoplastic micromechanical model is used to describe the plasticity-induced hysteresis in the ultrasonically measured interfacial stiffness during loading-unloading cycles. The topographic parameters of the interface contact are reconstructed by matching the model-predicted results with the experimentally determined ultrasonic stiffness. Using these parameters the real area of contact, which is not directly measurable, is predicted during loading-unloading cycles using the model.     

Frictional Hertzian contact problems under cyclic loading using static reduction

In this study we investigate an axisymmetric Hertzian contact problem of a rigid sphere pressing into an elastic half-space under cyclic loading. A numerical solution is sought to obtain a steady state, which demands a large amount of computer memory and computing speed. To achieve a tractable problem, the current numerical model uses a “static reduction” technique, which employs only a contact stiffness matrix rather than the entire stiffness of the problem and is more accurate than the approach used by most finite element codes. Investigation of the tendency of contact behavior in the transient and steady states confirms that a steady state exists, showing converged energy dissipation. The dependence of dissipation on load amplitude shows a power law of load amplitude less than 3, which may explain some deviations in the experimental findings.     

线接触弹性接触变形的解析算法

以一般光滑性体接触理论为基础,结合有限长弹性体接触的特点,求出线接触弹性接蟹变形的解析公式,并发现其解析解与数值解具有很好的一致性,所得公式可以对赫兹线接触理论加以补充,与经验公式相比,它能够确切反映材料、载荷以及曲率半径等对接触变形的影响,为工程中的精确计算提供了方便。     

摩擦接触弹塑性分析的数学规划法

本文提出摩擦接触弹塑性分析的数学规划法,是弹性接触问题解法的推广,使之能处理比例加载时的塑性效应。利用所导出接触问题总刚度阵的逆,结合参变量线性互补公式,仍可消去由接触单元而引入的惩罚因子。于是,可以获得一个弹塑性问题、弹性接触和弹塑性接触问题的有限元参数二次规划法模型。所选的几个实例说明了方法的实用性。     

Effective elastic modulus of film-on-substrate systems under normal and tangential contact

Load and depth sensing indentation methods have been widely used to characterize the mechanical properties of the thin film-substrate systems. The measurement accuracy critically depends on our knowledge of the effective elastic modulus of this heterogeneous system. In this work, based on the exact solution of the Green's function in Fourier space, we have derived an analytical relationship between the surface tractions and displacements, which depends on the ratio of the film thickness to contact size and the generalized Dundurs parameters that describe the modulus mismatch between the film and substrate materials. The use of the cumulative superposition method shows that the contact stiffness of any axisymmetric contact is the same as that of a flat-ended punch contact. Therefore, assuming a surface traction of the form of [1−(r/a)²]^−1/2 with radial coordinate r and contact size a, we can obtain an approximate representation of the effective elastic moduli, which agree extremely well with the finite element simulations for both normal and tangential contacts. Motivated by a recently developed multidimensional nanocontact system, we also explore the dependence of the ratio of tangential to normal contact stiffness on the ratio of film thickness to contact radius and the Dundurs parameters. The analytical representations of the correction factors in the relationship between the contact stiffness and effective modulus are derived at infinite friction conditions.