New limit cycles of dry friction oscillators under harmonic load

We consider a system composed of two masses connected by linear springs. One of the masses is in contact with a driving belt moving at a constant velocity. Friction force, with Coulomb??s characteristics, acts between the mass and the belt. Moreover, the mass is also subjected to a harmonic external force. Several periodic orbits including stick phases and slip phases are obtained. In particular, the existence of periodic orbits including a part where the mass in contact with the belt moves in the same direction at a higher speed than the belt itself is proved. Non-sticking orbits are also found for a non-moving belt. We prove that this kind of solution is symmetric in space and in time.     

On the dynamics of the rigid body with a fixed point: periodic orbits and integrability

The aim of the present paper is to study the periodic orbits of a rigid body with a fixed point and quasi-spherical shape under the effect of a Newtonian force field given by different small potentials. For studying these periodic orbits, we shall use averaging theory. Moreover, we provide information on the $\mathcal{C}^{1}$ -integrability of these motions.     

基于分子动力学-格林函数法的微凸体接触数值分析

表面接触是摩擦的先决条件,其真实接触面积、压应力大小、空间分布等一直是接触力学关注的核心问题.采用分子动力学-格林函数法(GFMD)模拟粗糙面的接触过程,验证了其在大规模接触分析中的高效及准确性,同时探讨了由微球体组成的粗糙面的接触力学特性,并分析了分子尺度下的结果和传统力学模型计算结果的差异.结果表明,单个微凸体接触结果和分子动力学-格林函数法模拟所得非常接近,误差在5%以内.数值模拟发现,在微凸体高度符合高斯分布的情况下,接触面积和接触力成线性关系;在相同接触面积下,微凸体模型得出的接触力偏高,是上限值.微凸体模型没有考虑微凸体间的相互影响,实际是高估了弹性体的刚度;实际接触过程中微凸体相互影响,微凸体对临域形变影响尤其大,使接触区域更加离散.GFMD模型可以准确计算数十亿量级别分子、原子接触过程中真实接触面积及分布,为后续摩擦、滑移等分析提供可靠的参考.     

Adhesive contact of an elastic semi-infinite solid with a rigid rough surface: Strength of adhesion and contact instabilities

The effect of adhesion on the contact behavior of elastic rough surfaces is examined within the framework of the multi-asperity contact model of Greenwood and Williamson (1966), known as the GW model. Adhesive surface interaction is modeled by nonlinear springs with a force–displacement relation governed by the Lennard–Jones (LJ) potential. Constitutive models are presented for contact systems characterized by low and high Tabor parameters, exhibiting continuous (stable) and discontinuous (unstable) surface approach, respectively. Constitutive contact relations are obtained by integrating the force–distance relation derived from the LJ potential with a finite element analysis of single-asperity adhesive contact. These constitutive relations are then incorporated into the GW model, and the interfacial force and contact area of rough surfaces are numerically determined. The development of attractive and repulsive forces at the contact interface and the occurrence of instantaneous surface contact (jump-in instability) yield a three-stage evolution of the contact area. It is shown that the adhesion parameter introduced by Fuller and Tabor (1975) governs the strength of adhesion of contact systems with a high Tabor parameter, whereas the strength of adhesion of contact systems with a low Tabor parameter is characterized by a new adhesion parameter, defined as the ratio of the surface roughness to the equilibrium interatomic distance. Applicable ranges of aforementioned adhesion parameters are interpreted in terms of the effective surface separation, obtained as the sum of the effective distance range of the adhesion force and the elastic deformation induced by adhesion. Adhesive strength of rough surfaces in the entire range of the Tabor parameter is discussed in terms of a generalized adhesion parameter, defined as the ratio of the surface roughness to the effective surface separation.     

Elasto-plastic contact of rough surfaces: a mixed-lubrication model for the textured surface analysis

An improved asperity contact model for two rough surfaces with misalignment is presented in this study. The contact model is statistical and can account for the elasto-plastic behavior of interacting asperities. By combining the improved asperity contact model and the average flow Reynolds equation together, a mixed-lubrication model is developed to understand the effect of surface texturing. By comparing with the results of the purely elastic asperity contact model, it is found that the improved asperity contact model can predict the contact force and actual contact area more accurately, particularly under high load conditions. Moreover, comparing with the elasto-plastic model with an equivalent rough surface against a plane, the improved contact model can consider the effect of permitting misalignment of two rough surfaces. This is beneficial for analyzing the performance of the textured piston ring/liner system, especially when asperities contact and wear happen.     

粗糙表面法向接触刚度的分形模型

提出了以往有关粗糙表面法向接触刚度理论研究工作的缺陷与不足,并在一定的前提假设下,基于球体与平面的接触理论和粗糙表面的分形接触理论,从理论上给出了具有尺度独立性的粗糙表面法向接触刚度分形模型,并进行了数字仿真研究。     

The effect of oblateness in the perturbed restricted three-body problem

In this paper the restricted three-body problem is generalized in the sense that the effects of oblateness of the three participating bodies as well as the small perturbations in the Coriolis and centrifugal forces are considered. The existence of equilibrium points, their linear stability and the periodic orbits around these points are studied under these effects. It is found that the positions of the collinear points and y-coordinate of the triangular points are not affected by the small perturbations in the Coriolis force. While x-coordinate of the triangular points is neither affected by the small perturbations in the Coriolis force nor the oblateness of the third body. Furthermore, the critical mass value and the elements of periodic orbits around the equilibrium points such as the semi-major and the semi-minor axes, the angular frequencies and corresponding periods may change by all the parameters of oblateness as well as the small perturbations in the Coriolis and centrifugal forces. This model could be applicable to send satellite or place telescope in stable regions in space.     

KINETICS ANALYSIS OF NON-SMOOTH VIBRATION-DRIVEN SYSTEM WITH THREE-PHASE CONTROL 1)

可移动式机器人已成为机器人研究领域的重要分支,为实现其在狭小特殊环境中的运动, 学者们提出并研究了振动驱动移动系统.本文基于二维LuGre摩擦模型和拉格朗日方程,给出了一类振动驱动系统在各向同性摩擦环境中的动力学建模方法和数值算法.这类振动驱动系统结构简单且密封性好,依靠箱体与地面间的摩擦力实现自身的定向运动.该系统由一个外部箱体和两个内部质量块构成,两个质量块在箱体内的两个平行轨道上作三相振动驱动,箱体通过三个刚性支撑足与地面保持接触. 二维LuGre摩擦模型的利用,可有效避免库伦摩擦模型的不连续性给动力学方程的数值求解带来的困难,且可有效揭示该系统在运动过程中的黏滞-滑移切换现象. 数值仿真结果表明,通过调整其内部质量块的驱动参数,可实现箱体的直线平移、定轴转动和平面一般运动,且箱体在移动和转动过程中会出现擦滑、穿滑、回滑和不黏等4种现象; 另外,通过调节驱动参数, 不仅可以改变箱体移动和转动的快慢,还可以改变箱体形心运动轨迹的曲率半径.     

基于基底长度的粗糙表面分形接触模型的构建与分析

为建立更为精确的粗糙表面接触模型,根据微凸体变形特征、分形理论以及摩擦的作用,从微观角度基于基底长度建立了粗糙表面分形接触模型.通过与其他粗糙表面接触模型和实验数据的比较,验证了本文模型的正确与合理,并由数值仿真分析了分形维数、接触载荷与真实接触面积之间的相互关系.分析结果表明:特征尺度一定时,要维持一定的真实接触面积,分形维数越大,所需要的力也越大;分形维数与特征尺度一定时,随着载荷的增加,接触面积也在增加;特征尺度与接触力一定时,随着分形维数的增大,真实接触面积在减小.该模型的建立为进一步研究粗糙表面的摩擦、磨损与润滑提供了理论依据.     

三相驱动下非光滑振动驱动系统的动力学分析

可移动式机器人已成为机器人研究领域的重要分支,为实现其在狭小特殊环境中的运动, 学者们提出并研究了振动驱动移动系统.本文基于二维LuGre摩擦模型和拉格朗日方程,给出了一类振动驱动系统在各向同性摩擦环境中的动力学建模方法和数值算法.这类振动驱动系统结构简单且密封性好,依靠箱体与地面间的摩擦力实现自身的定向运动.该系统由一个外部箱体和两个内部质量块构成,两个质量块在箱体内的两个平行轨道上作三相振动驱动,箱体通过三个刚性支撑足与地面保持接触. 二维LuGre摩擦模型的利用,可有效避免库伦摩擦模型的不连续性给动力学方程的数值求解带来的困难,且可有效揭示该系统在运动过程中的黏滞-滑移切换现象. 数值仿真结果表明,通过调整其内部质量块的驱动参数,可实现箱体的直线平移、定轴转动和平面一般运动,且箱体在移动和转动过程中会出现擦滑、穿滑、回滑和不黏等4种现象; 另外,通过调节驱动参数, 不仅可以改变箱体移动和转动的快慢,还可以改变箱体形心运动轨迹的曲率半径.      

微重力下圆柱形贮箱内液体晃动的分岔现象

为了得到微重力下液体晃动的特性,晃动模态用具有水平静液面的液体晃动模态近似. 液体晃动取前5 个重要模态,用Lagrange 方程推导了横向力作用下液体晃动和航天器结构的无量纲的耦合动力学方程并进行数值模拟. 模拟了液体晃动模态随外力振幅和频率变化产生的分岔现象,并分析了系统参数,如Bond数、接触角、接触角迟滞、充液高度、频率比、质量比以及外力的周期形式和方向等对晃动模态的分岔的影响.     

Relative equilibria with clusters of small masses

Relative equilibria of the Newtoniann-body problem are studied in the limiting case where some of the masses tend to zero. The small masses may form clusters. Equations for these clusters are derived and analyzed. Stability conditions for the resulting periodic orbits of then-body problem are given.     

Hybrid dynamics of two coupled oscillators that can impact a fixed stop

We consider two linearly coupled masses, where one mass can have inelastic impacts with a fixed, rigid stop. This leads to the study of a two degree of freedom, piecewise linear, frictionless, unforced, constrained mechanical system. The system is governed by three types of dynamics: coupled harmonic oscillation, simple harmonic motion and discrete rebounds. Energy is dissipated discontinuously in discrete amounts, through impacts with the stop. We prove the existence of a non-zero measure set of orbits that lead to infinite impacts with the stop in a finite time. We show how to modify the mathematical model so that forward existence and uniqueness of solutions for all time is guaranteed. Existence of hybrid periodic orbits is shown. A geometrical interpretation of the dynamics based on action coordinates is used to visualize numerical simulation results for the asymptotic dynamics.     

Bifurcation analysis of a smoothed model of a forced impacting beam and comparison with an experiment

A piecewise-linear model with a single degree of freedom is derived from first principles for a driven vertical cantilever beam with a localized mass and symmetric stops. The aim is to show that this model constitutes a considerable step toward developing a vibro-impact model that is able to make qualitative and quantitative predictions of the observed dynamics. The resulting piecewise-linear dynamical system is smoothed by a switching function (nonlinear homotopy). For the chosen smoothing function, it is shown that the smoothing can induce bifurcations in certain parameter regimes. These induced bifurcations disappear when the transition of the switching is sufficiently and increasingly localized as the impact becomes harder. The bifurcation structure of the impact oscillator response is investigated via the one- and two-parameter continuation of periodic orbits in the driving frequency and/or forcing amplitude. The results are in good agreement with experimental measurements.     

Elucidating the escape dynamics of the four hill potential

This paper treads discontinuous bifurcation in piecewise smooth systems of Filippov type. These bifurcations occur when a fixed point or a periodic orbit crosses with the border between two regions of smooth behavior. A detailed analysis of generalization Poincaré map and monodromy matrix which are related shows that subfamily of system with invariant cone-like objects is foliated by periodic orbits and determines its stability. In addition, we introduce a theoretical framework for analyzing 3D perturbed nonlinear piecewise smooth systems and give necessary conditions so that different types of bifurcations occur. The analysis identifies criteria for the existence of a novel bifurcation based on sensitively the location of the return map. Moreover, the piecewise smooth Melnikov function and sufficient conditions of the existence of the periodic orbits for nonlinear perturbed system are explicitly obtained.     

Study of a transition in the qualitative behavior of a simple oscillator with Coulomb friction

A simple mass-spring system is submitted to a constant force in addition to a periodic perturbation of rectangular wave shape. It has been obtained in a previous study that the range of the period-amplitude plane of this perturbation, where the trajectories are sliding with no loss of contact, is divided into two parts, one in which there exist infinitely many equilibrium states and no periodic solutions, and another one where there exist periodic solutions and no equilibrium states. The present work focuses on the transition between these two parts. All along the transition line, there exists a single equilibrium state. Initial data out of equilibrium lead either to a periodic trajectory, or to a trajectory, which tends to the equilibrium or to a periodic solution, either in finite time or at infinity.     

Dynamics of a mass–spring–pendulum system with vastly different frequencies

We investigate the dynamics of a simple pendulum coupled to a horizontal mass?Cspring system. The spring is assumed to have a very large stiffness value such that the natural frequency of the mass?Cspring oscillator, when uncoupled from the pendulum, is an order of magnitude larger than that of the oscillations of the pendulum. The leading order dynamics of the autonomous coupled system is studied using the method of Direct Partition of Motion (DPM), in conjunction with a rescaling of fast time in a manner that is inspired by the WKB method. We particularly study the motions in which the amplitude of the motion of the harmonic oscillator is an order of magnitude smaller than that of the pendulum. In this regime, a pitchfork bifurcation of periodic orbits is found to occur for energy values larger that a critical value. The bifurcation gives rise to nonlocal periodic and quasi-periodic orbits in which the pendulum oscillates about an angle between zero and ??/2 from the down right position. The bifurcating periodic orbits are nonlinear normal modes of the coupled system and correspond to fixed points of a Poincare map. An approximate expression for the value of the new fixed points of the map is obtained. These formal analytic results are confirmed by comparison with numerical integration.     

滤波截止频率对粗糙点接触机理影响研究

使用基于快速傅里叶变换和共轭梯度法的接触程序,分析了滤波截止频率对粗糙点接触机理的影响.计算结果表明:粗糙表面形貌数据再处理时使用的滤波器滤波截止频率,对力-接近量关系有着明显的影响;滤波截止频率越小,接触面积越大,接触压力越小;滤波截止频率的改变会导致Von M ises应力的显著改变,这种改变在接触面上特别明显.因此,在实际研究粗糙表面的接触机理时,选择合适的滤波截止频率是非常必要的.     

A station-keeping strategy for collinear libration point orbits

Spacecrafts in periodic or quasi-periodic orbits near the collinear libration points are proved to be excellent platforms for scientific investigations of various phenomena. Since such periodic or quasi-periodic orbits are exponentially unstable, the station-keeping maneuver is needed. A station-keeping strategy which is found by an analytical method is presented to eradicate the dominant unstable component of the libration point trajectories. The inhibit force transforms the unstable component to a stable component. In this method, it is not necessary to determine a nominal orbit as a reference path.