两种型面轮轨滚动接触蠕滑率和摩擦功

采用数值分析方法详细分析了2种型面轮对与轨道在滚动接触过程中的接触几何、蠕滑率和摩擦功。在摩擦功分析计算中利用了Kalker三维弹性体非赫兹滚动接触理论,并考虑了轮轨的结构变形的影响。分析计算表明,在轮对运动相同的情况下,磨耗型轮对与轨道之间的蠕滑率和摩擦功与锥型轮对轨道之间的蠕滑率和摩擦功相差甚大。数值分析结果表明,现行铁路中正在推广使用的磨耗型轮对并不能减少轮轨接触表面之间的磨损,磨损型车轮的型面与钢轨尺寸的匹配还有待进一步改进。     

机械系统摩擦动力学研究进展

摩擦环节对机械系统动力学行为有重要、甚至可能是关键的影响. 深入研究摩擦及摩擦动力学特点,对于解决机械系统中摩擦带来的不利影响, 发挥其有利作用, 是非常重要的. 本文介绍和评述了机械系统摩擦动力学的研究进展, 包括常用的摩擦模型及其特性, 摩擦系统自激振动、强迫振动和摩擦振动控制等. 除理论研究方法之外, 重点讨论了制动噪声振动、摩擦耗能和多领域摩擦振动控制等摩擦动力学的应用研究.     

原子尺度摩擦行为的分子动力学模拟

应用大规模分子动力学方法,模拟了具有原子级光滑和原子级粗糙形貌的刚性球形探头与弹性平面基体的干摩擦行为,研究了无/有粘附条件下的载荷与摩擦力、载荷与真实接触面积,以及摩擦力与真实接触面积之间的关系,对纳米尺度下的摩擦行为规律进行了分析。几种系统的真实接触面积－载荷关系都与相应的连续力学接触模型定性的一致,它们分别是Hertz光滑表面接触模型、Greenwood-Williamson粗糙表面接触模型和Maugis-Dugdale粘着接触模型。无论是由光滑表面还是粗糙表面构成的摩擦系统,在无粘附条件下摩擦力与载荷成正比,而摩擦力与真实接触面积之间没有一个简单的关系;在粘附条件下摩擦力与真实接触面积成正比,而摩擦力与载荷之间表现为Maugis-Dugdale模型预测的亚线性关系。我们的研究表明,当表面作用从无粘附到粘附时,控制摩擦力的决定因素从载荷转变为接触面积,摩擦行为从载荷控制摩擦转变为粘着控制摩擦。     

Influence of rolling friction on single spout fluidized bed simulation

In this paper we study the effect of rolling friction on the dynamics in a single spout fluidized bed using Discrete Element Method (DEM) coupled to Computational Fluid Dynamics (CFD). In a first step we neglect rolling friction and show that the results delivered by the open source CFD–DEM framework applied in this study agree with previous simulations documented in literature. In a second step we include a rolling friction sub-model in order to investigate the effect of particle non-sphericity. The influence of particle–particle as well as particle–wall rolling friction on the flow in single spout fluidized bed is studied separately. Adequate rolling friction model parameters are obtained using first principle DEM simulations and data from literature. Finally, we demonstrate the importance of correct modelling of rolling friction for coupled CFD–DEM simulations of spout fluidized beds. We show that simulation results can be improved significantly when applying a rolling friction model, and that experimental data from literature obtained with Positron Emission Particle Tracking (PEPT) technique can be satisfactorily reproduced.     

Rigid body impact with moment of rolling friction

In this paper, the impact between a rigid pendulum and rough surfaces is studied. The rolling friction moment and the coefficient of rolling friction are introduced, and an improved mathematical model of the planar impact with friction is presented. The influence of the moment of rolling friction on the energy dissipated by friction during the impact is analyzed. For a simple pendulum, using the energetic coefficient of restitution, more energy is dissipated for larger values of the coefficient of kinetic friction and contact radius, and for smaller values of the length of the beam. For a double pendulum using the kinematic coefficient of restitution, some energetically inconsistent results can be solved for some values of the coefficient of rolling friction.     

轮轨滚动摩擦温升分析

利用有限元法,考虑轮轨间非稳态热传导、与环境的热对流以及热辐射的影响,建立了轮轨滚动接触热耦合计算模型来模拟轮轨滚滑摩擦温升;在模拟轮轨纯滑动条件下,计算分析了由磨损引起的滑动接触斑的尺寸增大对轮轨温度场的影响;在模拟轮轨接触斑部分滑动工况时,针对不同蠕滑率、摩擦系数以及轴重对轮轨温度场的影响进行了相应的计算分析.结果表明:接触斑材料的磨损速度只影响磨损过程中的温度场分布,其稳态温度场分布基本一致;热载荷随着纵向载荷、蠕滑率以及摩擦系数的增大而增大,进而影响轮轨滚动接触热疲劳.     

非公度接触石墨烯层间摩擦对压入深度的依赖性

近年来, 以石墨烯为代表的层状二维纳米材料的摩擦力学行为受到广泛关注, 许多新的纳尺度摩擦现象、规律及机理被陆续报道, 推动纳米摩擦学取得了重要进展. 然而, 由于纳米级摩擦十分复杂, 在建立摩擦力与影响因素之间的直接关联方面依然进展非常缓慢. 本文利用分子动力学模拟方法, 研究了衬底支撑石墨烯基底与石墨烯滑片之间的摩擦行为, 着力考察了非公度接触情况下的摩擦规律. 结果表明, 石墨烯滑片和基底之间的摩擦力与压入深度直接相关, 说明压入深度可作为纳尺度摩擦力的重要度量指标. 特别地, 法向载荷和衬底刚度对石墨烯摩擦的影响，都可通过压入深度归一化处理. 该结果对理解二维材料表面弹性影响的摩擦规律具有重要的理论意义.     

Chaotic behavior of a single spherical gas bubble surrounded by a Giesekus liquid: A numerical study

In the present work, nonlinear oscillations of a spherical, acoustically driven gas bubble in a Giesekus liquid are examined numerically. A novel approach based on the Gauss–Laguerre quadrature (GLQ) method is implemented to solve the integro-differential equation governing bubble dynamics in a Giesekus liquid. It is shown that, using this robust method, numerical results could be obtained at very high amplitudes and frequencies typical of ultrasound applications. The GLQ method also enabled obtaining results at very high Deborah and Reynolds numbers over prolonged dimensionless times not reported previously. Based on the results obtained in this work, it is concluded that the GLQ method is well suited for bubble dynamics studies in viscoelastic liquids. It is also concluded that the extensional-flow behavior of the liquid surrounding the bubble (as represented by the mobility factor in the Giesekus model) has a strong effect on the chaotic behavior of the bubble, and this is particularly so at high Deborah numbers, high amplitudes and/or high frequencies of the acoustic field. A period-doubling bifurcation structure is predicted to occur for certain values of the mobility factor.     

非公度接触石墨烯层间摩擦对压入深度的依赖性

近年来, 以石墨烯为代表的层状二维纳米材料的摩擦力学行为受到广泛关注, 许多新的纳尺度摩擦现象、规律及机理被陆续报道, 推动纳米摩擦学取得了重要进展. 然而, 由于纳米级摩擦十分复杂, 在建立摩擦力与影响因素之间的直接关联方面依然进展非常缓慢. 本文利用分子动力学模拟方法, 研究了衬底支撑石墨烯基底与石墨烯滑片之间的摩擦行为, 着力考察了非公度接触情况下的摩擦规律. 结果表明, 石墨烯滑片和基底之间的摩擦力与压入深度直接相关, 说明压入深度可作为纳尺度摩擦力的重要度量指标. 特别地, 法向载荷和衬底刚度对石墨烯摩擦的影响,都可通过压入深度归一化处理. 该结果对理解二维材料表面弹性影响的摩擦规律具有重要的理论意义.     

An Approximate Analysis of Dry-Friction-Induced Stick-Slip Vibrations by a Smoothing Procedure

This paper deals with a systematic procedure to find both stable and unstable periodic stick-slip vibrations of autonomous dynamic systems with dry friction. In this procedure, the discontinuous friction forces are approximated by smooth functions. Using the simple shooting method with a stiff-ODE solver, in combination with a path following algorithm, branches of periodic solutions are computed for a changing design variable. For testing purposes, both 1 and 2-DOF autonomous block-on-belt models and a 1-DOF autonomous drill string model from literature are investigated. Comparison of the results shows that the smoothing procedure accurately describes the behavior of the discontinuous systems. The proposed procedure can also easily be applied to more complex MDOF models, as well as to nonautonomous dynamic systems.     

A hybrid scheme for simulation of planar rigid bodies with impacts and friction using impact mappings

This paper introduces numerical techniques necessary for the implementation of impact maps derived from an energetic impact law for rigid-body impacts with friction at isolated contact points. In particular the work focuses on methodologies for long-term simulation with behaviours such as dynamic transitions and chatter. The methods are based on hybrid event-driven numerical solvers for ordinary differential equations together with system states to deal with the transitions. A slender rod impacting a periodically oscillating surface is used as an example to illustrate implementation and methods. The numerical scheme for the rod system is used to show how symmetry can play an important role in the presence of friction for long-term dynamics. This will show that surface oscillations with low frequencies tend to lead to periodic motions of the rod that are independent of friction. For higher frequencies however the periodic solutions are not that common and irregular motion ensues.     

Dynamics of a homogeneous ball on a horizontal plane with sliding,spinning, and rolling friction taken into account

We analyze the dynamics of a homogeneous ball on a horizontal plane with friction of all kinds, namely, sliding, spinning, and rolling friction, taken into account. The qualitative-analytic study of the ball dynamics is supplemented with numerical experiments.     

Analytical approximations for stick-slip vibration amplitudes

The classical “mass-on-moving-belt” model for describing friction-induced vibrations is considered, with a friction law describing friction forces that first decreases and then increases smoothly with relative interface speed. Approximate analytical expressions are derived for the conditions, the amplitudes, and the base frequencies of friction-induced stick-slip and pure-slip oscillations. For stick-slip oscillations, this is accomplished by using perturbation analysis for the finite time interval of the stick phase, which is linked to the subsequent slip phase through conditions of continuity and periodicity. The results are illustrated and tested by time-series, phase plots and amplitude response diagrams, which compare very favorably with results obtained by numerical simulation of the equation of motion, as long as the difference in static and kinetic friction is not too large.     

车轮钢滚动剥离摩擦磨损特性研究

在NENE-2型摩擦磨损试验机上利用往复滚动试验装置研究了不同滚滑状态下车轮钢的剥离摩擦磨损特性和碳含量对车轮钢滚动剥离磨损性能的影响.结果表明:在不同滚滑状态下摩擦副之间的摩擦力不同,平面试样的表面磨痕形貌随着不同的切向摩擦力而明显不同,随着切向摩擦力的增大滚动磨损机制亦发生改变,剥离磨损加剧且磨损深度变大,当相对滑动量增大到一定程度后,磨损表现为明显的剥层机制;碳含量对车轮钢的滚动磨损表面磨痕形貌影响显著,碳含量低时磨痕以犁沟为主,碳含量高时剥离磨损发生的概率增加.     

车辆纵向非光滑多体动力学建模与数值算法研究

在建立车辆纵向多体系统的动力学模型中, 将车身与车轮视为刚体, 两者通过减振器链接; 将传动系统视为一个圆盘通过扭簧和阻尼器与驱动轮连接; 将车轮与路面间的接触力简化为法向约束力、摩擦力和滚阻力偶,其中摩擦力的模型采用库仑干摩擦模型. 采用笛卡尔坐标作为该系统的广义坐标用于描述该系统的位形, 给出系统单双边的约束方程, 应用第一类拉格朗日方法建立了系统的动力学方程. 由于摩擦与滚阻的非光滑性, 使得该系统动力学方程不连续. 为便于计算, 建立了车轮与路面接触点的相对切向加速度与摩擦力余量的互补条件、车轮角加速度与滚阻力偶余量的互补条件, 以及车轮轮心法向加速度与路面法向约束力的互补条件. 将接触—分离、黏滞—滑移的判断问题转化成线性互补问题的求解, 并给出了具有约束稳定化的基于事件驱动法的数值计算方法. 最后, 应用该方法对车辆纵向多体系统进行了仿真, 分析了输出扭矩、摩擦及滚阻系数对其动力学行为的影响.      

Dynamics of mechanical systems involving impact and friction using an efficient contact detection algorithm

In this work, some new results are presented on the dynamics of a class of multibody mechanical systems, involving contact and friction. The main contribution refers to the development of a systematic, accurate and efficient method for detecting contact among the components of a system of solid bodies. For some simple geometries, this task is achieved by employing analytical means. For systems possessing components with complex geometric shapes a more involved numerical methodology is developed. In both cases, once a potential contact point is detected, the common tangent plane and normal vector are located and the penetration depth is calculated, leading to determination of the force arising between the contacting bodies. This information is then passed to a solver, providing the full dynamic response of the system. The validity and numerical efficiency of the methodology developed is first demonstrated by considering a number of examples with relatively small geometric complexity but large traditional value and interesting dynamic response. Some new results are obtained and presented on the dynamics of these systems. Finally, the same methodology is also tested in a more complicated and demanding mechanical application.     

On the dynamics of a rolling ball actuated by internal point masses

The motion of a rolling ball actuated by internal point masses that move inside the ball’s frame of reference is considered. The equations of motion are derived by applying Euler–Poincaré’s symmetry reduction method in concert with Lagrange–d’Alembert’s principle, which accounts for the presence of the nonholonomic rolling constraint. As a particular example, we consider the case when the masses move along internal rails, or trajectories, of arbitrary shape and fixed within the ball’s frame of reference. Our system of equations can treat most possible methods of actuating the rolling ball with internal moving masses encountered in the literature, such as circular motion of the masses mimicking swinging pendula or straight line motion of the masses mimicking magnets sliding inside linear tubes embedded within a solenoid. Moreover, our method can model arbitrary rail shapes and an arbitrary number of rails such as several ellipses and/or figure eights, which may be important for future designs of rolling ball robots. For further analytical study, we also reduce the system to a single differential equation when the motion is planar, that is, considering the motion of the rolling disk actuated by internal point masses, in which case we show that the results obtained from the variational derivation coincide with those obtained from Newton’s second law. Finally, the equations of motion are solved numerically, illustrating a wealth of complex behaviors exhibited by the system’s dynamics. Our results are relevant to the dynamics of nonholonomic systems containing internal degrees of freedom and to further studies of control of such systems actuated by internal masses.     

On a mechanical lens

In this paper, we consider the dynamics of a heavy homogeneous ball moving under the influence of dry friction on a fixed horizontal plane. We assume the ball to slide without rolling. We demonstrate that the plane may be divided into two regions, each characterized by a distinct coefficient of friction, so that balls with equal initial linear and angular velocity will converge upon the same point from different initial locations along a certain segment. We construct the boundary between the two regions explicitly and discuss possible applications to real physical systems.     

Angle of repose in the numerical modeling of ballast particles focusing on particle-dependent specifications: Parametric study

The discrete element method (DEM) is widely used in the realistic simulation of the shapes of particles. Researchers have considered the simplification of particle shapes owing to the high computational cost of such simulation. In this regard, the modeling of calibrated particles is a major challenge owing to the simultaneous effects of particle properties. The angle-of-repose test is a standard test method used to calibrate the bulk behavior of simulated particles. In the present study, the hollow-cylinder (slump) test was modeled for the verification of discrete element simulations. In this regard, a sensitivity analysis was conducted for all effective parameters, namely the static friction, rolling friction, restitution coefficient, sphericity, roundness, particle size distribution, and number of ballast particles. The results indicate that the rolling friction, roundness, number of particles, and size of particles are the most important parameters in the determination of the angle of repose (AOR). For particles in the range of ballast (20–60 mm), the effect of the number of particles on the angle of repose is reduced when the number is greater than 426. Additionally, it is concluded that angular particles can be replaced with sub-angular particles (R ≈ 0.2–0.45) with a higher rolling friction coefficient (μ_r > 0.14).