Nonlinear shear wave interaction at a frictional interface: energy dissipation and generation of harmonics

Analytical and numerical modeling of the nonlinear interaction of shear wave with a frictional interface is presented. The system studied is composed of two homogeneous and isotropic elastic solids, brought into frictional contact by remote normal compression. A shear wave, either time harmonic or a narrow band pulse, is incident normal to the interface and propagates through the contact. Two friction laws are considered and the influence on interface behavior is investigated: Coulomb's law with a constant friction coefficient and a slip-weakening friction law which involves static and dynamic friction coefficients. The relationship between the nonlinear harmonics and the dissipated energy, and the dependence on the contact dynamics (friction law, sliding, and tangential stress) and on the normal contact stress are examined in detail. The analytical and numerical results indicate universal type laws for the amplitude of the higher harmonics and for the dissipated energy, properly non-dimensionalized in terms of the pre-stress, the friction coefficient and the incident amplitude. The results suggest that measurements of higher harmonics can be used to quantify friction and dissipation effects of a sliding interface.     

Micro/macro-slip damping in beams with frictional contact interface

Friction in contact interfaces of assembled structures is the prime source of nonlinearity and energy dissipation. Determination of the dissipated energy in an assembled structure requires accurate modeling of joint interfaces in stick, micro-slip and macro-slip states. The present paper proposes an analytical model to evaluate frictional energy loss in surface-to-surface contacts. The goal is to develop a continuous contact model capable of predicting the dynamics of friction interface and dissipation energy due to partial slips. To achieve this goal, the governing equations of a frictional contact interface are derived for two distinct contact states of stick and partial slip. A solution procedure to determine stick–slip transition under single-harmonic excitations is derived. The analytical model is verified using experimental vibration test responses performed on a free-frictionally supported beam under lateral loading. The theoretical and experimental responses are compared and the results show good agreements between the two sets of responses.     

Effect of substrate dimensions on dynamics of ultrasonic consolidation

An FEM model is developed for a fundamental study of the time-dependent mechanical behavior of the substrate and its dimensions on ultrasonic consolidation. The simulation shows that for a given vibration condition, the amplitude of contact friction stress and displacement stabilizes to a saturated state after certain number of ultrasonic cycles. With the increased substrate height, the amplitude of contact frictional stress decreases, while that of contact interface displacement increases. The reason for the decrease in the frictional stress at the contact interface for certain substrate heights is the complicated wave interference occurring in the substrate. An analytical wave model has been built to validate the FEM model. A specific substrate geometry (height:width = 1.0) generates a minimum frictional strain state at the interface as a result of wave superposition. Such minimum strain state is believed to have produced the “lack of bonding” defect for the geometry. The energy density and transfer coefficient at the contact interface with different substrate heights is used as an indicator to correlate with the bond formation in ultrasonic bonding.     

A model for the characterization of friction contacts in turbine blades

Stresses produced by the forced vibrations can lead to a significant reduction of the life of turbo engine blades. To predict the vibration amplitudes of this components an accurate dynamic analysis is necessary. The forced response calculation of these dynamic systems is strongly affected by the presence of the contact interfaces (i.e., underplatform dampers, shrouds, root joints). Different contact models are available in literature. These models make use of contact parameters, contact stiffness and friction coefficient to evaluate the damping and stiffness related to the contact interfaces. In this paper a model is proposed to characterize friction contact of non-spherical contact geometries obeying the Coulomb friction law with constant friction coefficient and constant normal load. The hysteresis curves of the oscillating tangential contact forces vs. relative tangential displacements and the dissipated energy at the contact are obtained for different contact geometries. The developed model is suitable to be implemented in numerical solvers for the calculation of the forced response of turbine blades with embedded friction contacts.     

Modal Amplitude Stability Analysis and its application to brake squeal

In the present study, a new approach is proposed to predict the occurrence of squeal in brake systems. This strategy, called Modal Amplitude Stability Analysis (MASA), is based on the calculation of the first harmonic state-space system of nonlinear original equations using a specific linearization of the nonlinear contact forces at the frictional interfaces. An estimation of the occurrence and generation of increasing self-excited vibration is proposed on the basis of monitoring and the evolution of the real parts of the dynamic system considered as a function of modal amplitudes.     

A dynamic Lagrangian frequency-time method for the vibration of dry-friction-damped systems

A new frequency-time domain procedure, the dynamic Lagrangian mixed frequency-time method (DLFT), is proposed to calculate the non-linear steady state response to periodic excitation of structural systems subject to dry friction damping. In this formulation, the dynamic Lagrangians are defined as the non-linear contact forces obtained from the equations of motion in the frequency domain, with the adjunction of a penalization on the difference between the interface displacements calculate by the non-linear solver in the frequency domain and those calculated in the time domain from the non-linear contact forces, thus accounting for Coulomb friction and non-penetration conditions. The dynamic Lagrangians allow one to solve for the non-linear forces between two points in contact without using artifacts such as springs. The new DLFT method is thus particularly well suited to handling finite element models of structures in frictional contact, as it does not require a special model for the contact interface. Dynamic Lagrangians are also better suited to frequency-domain friction problems than the traditional time-domain method of augmented Lagrangians. Furthermore, a reduction of the non-linear system to relative interface displacements is introduced to decrease the computation time. The DLFT method is validated for a beam in contact with a flexible dry friction element connected to ground, for frictional constraints that feature two-dimensional relative motion. Results are also obtained for a large-scale structural system with a large number of one-dimensional dry-friction dampers. The DLFT method is shown to be accurate and fast, and it does not suffer from convergence problems, at least in the examples studied.     

Stability and transient dynamics of a propeller–shaft system as induced by nonlinear friction acting on bearing–shaft contact interface

This paper investigates the friction-induced instability and the resulting self-excited vibration of a propeller–shaft system supported by water-lubricated rubber bearing. The system under consideration is modeled with an analytical approach by involving the nonlinear interaction among torsional vibrations of the continuous shaft, tangential vibrations of the rubber bearing and the nonlinear friction acting on the bearing–shaft contact interface. A degenerative two-degree-of-freedom analytical model is also reasonably developed to characterize system dynamics. The stability and vibrational characteristics are then determined by the complex eigenvalues analysis together with the quantitative analysis based on the method of multiple scales. A parametric study is conducted to clarify the roles of friction parameters and different vibration modes on instabilities; both the graphic and analytical expressions of instability boundaries are obtained. To capture the nature of self-excited vibrations and validate the stability analysis, the nonlinear formulations are numerically solved to calculate the transient dynamics in time and frequency domains. Analytical and numerical results reveal that the nonlinear coupling significantly affects the system responses and the bearing vibration plays a dominant role in the dynamic behavior of the present system.     

Wetting of phospholipid membranes on hydrophilic surfaces - Concepts towards self-healing membranes

We report on the wetting behavior of phospholipid membranes on solid surfaces immersed in aqueous solution. Using fluorescence microscopy, the spreading velocity of fluid bilayers advancing from a lipid source is investigated. The kinetic spreading coefficient was measured as a function of temperature for pure DMPC membranes and as a function of charge density and cholesterol content for binary membranes. A theoretical model for the membrane flow is presented, which takes into account the liquid crystalline bilayer architecture of the lipid membrane. The spreading power results from the membrane-solid VdW interaction and is dissipated in hydrodynamic shear flow as well as by inter-monolayer friction within the bilayer. The frictional drag causes a dynamic tension gradient in the spreading membrane, which is manifested by a single exponential decay of the fluorescence intensity profile along the spreading direction. Obstacles are shown to act as pinning centers deforming the advancing line interface. However, no depinning was observed, since the centers are circumflown without abrupt relaxation. Received 6 November 1998     

A method to determine dynamic loads on spur gear teeth and on bearings

A method is described which can be used to calculate dynamic gear tooth force and bearing forces. The model includes elastic bearings. The gear mesh stiffness and the path of contact are determined using the deformations of the gears and the bearings. This gives contact outside the plane-of-action and a time-varying working pressure angle. In a numerical example it is found that the only important vibration mode for the gear contact is the one where the gear tooth deformation is dominant. The bearing force variation, however, will be much more affected by the other vibration modes. The influence of the friction force is also studied. The friction has no dynamic influence on the gear contact force or on the bearing force in the gear mesh line-of-action direction. On the other hand, the changing of sliding directions in the pitch point is a source for critical oscillations of the bearings in the gear tooth frictional direction. These bearing force oscillations in the frictional direction appear unaffected by the dynamic response along the gear mesh line-of-action direction.     

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses.In the finite element model,the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method,and the dynamical interaction between both crack faces is modeled using a contact-impact theory.In the simulations,the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis,especially,the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail.Meanwhile,the related ultrasonic infrared images are also obtained experimentally,and the theoretical simulation results are in agreement with that of the experiments.The results show that,by using the theoretical method,a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.     

Mechanism analysis and improved model for stick-slip friction behavior considering stress distribution variation of interface

Studying the evolution of interface contact state, revealing the "black box" behavior in interface friction and establishing a more accurate friction model are of great significance to improve the prediction accuracy of mechanical system performance. Based on the principle of total reflection, a visual analysis technology of interface contact behavior is proposed. Considering the dynamic variation of stress distribution in interface contact, we analyze the nonlinear characteristics of contact parameters in different stages of stick-slip process using the above-mentioned experimental technology. Then, we find that the tangential stiffness of the interface is not a fixed value during the stick-slip process and the stress distribution variation is one of the important factors affecting the tangential stiffness of interface. Based on the previous experimental results, we present an improved stick-slip friction model, considering the change of tangential stiffness and friction coefficient caused by the stress distribution variation. This improved model can characterize the variation characteristics of contact parameters in different stages of stick-slip process, whose simulation results are in good agreement with the experimental data. This research may be valuable for improving the prediction accuracy of mechanical system performance.     

Local friction at a sliding interface between an elastomer and a rigid spherical probe

This paper reports on spatially resolved measurements of the shear stress distribution at a frictional interface between a flat rubber substrate and a glass lens. Silicone rubber specimens marked close to their surface by a colored pattern have been prepared in order to measure the surface displacement field induced by the steady-state friction of the spherical probe. The deconvolution of this displacement field then provides the actual shear stress distribution at the contact interface. When a smooth glass lens is used, a nearly constant shear stress is achieved within the contact. On the other hand, a bell-shaped shear stress distribution is obtained with rough lenses. These first results suggest that simple notions of real contact area and constant interface shear stress cannot account for the observed changes in local friction when roughness is varied.     

粗糙接触界面超声非线性效应的概率模型

提出描述粗糙接触界面超声非线性效应的概率模型,利用分段均匀概率函数描述粗糙接触界面劲度系数变化,结合表面粗糙峰的几何分布特征,得到界面粗糙度和两侧表面相对运动对声波的非线性调制作用。实验观测了铝合金材料的粗糙接触界面的高次谐波现象和阈值现象,进一步分析了归一化非线性参数与界面加载压力、粗糙度之间的关系。实验测量结果与概率模型的理论预测一致,证明了该模型的正确性,为利用超声非线性效应评价粗糙接触界面提供了理论依据。      

On linear modeling of interface damping in vibrating structures

Dissipation of mechanical vibration energy at contact interfaces in a structure, commonly referred to as interface damping, is an important source of vibration damping in built-up structures and its modeling is the focus of the present study. The approach taken uses interface forces which are linearly dependent on the relative vibration displacements at the contact interfaces.The main objective is to demonstrate a straightforward technique for simulation of interface damping in built-up structures using FE modeling and simple, distributed, damping forces localized to interfaces where the damping occurs.As an illustration of the resulting damping the dissipated power is used for evaluation purposes. This is calculated from surface integrals over the contact interfaces and allows for explicit assessment of the effect of simulated interface forces for different cases and frequencies. The resulting loss factor at resonance is explicitly evaluated and, using linear simulations, it is demonstrated that high damping levels may arise even though the displacement differences between contacting surfaces at damped interfaces may be very small.     

Transverse shear oscillator investigation of boundary lubrication in weakly adhered films

We investigate the boundary lubrication in weakly adhered molecularly thin films deposited between a sphere and a plane, below the sliding threshold. The shear contact stiffness and interfacial dissipation at the micrometer scale are determined with a high-frequency quartz oscillator. Two distinct behaviors are found as a function of the shear oscillation: a linear viscoelastic response at low amplitude and a nonlinear frictional microslip at high amplitude. A friction model is proposed to analyze the data, which allows evaluating the shear strength, the friction coefficient, and the interfacial viscosity at different solid interfaces under low load.     

Dynamics of precursors to frictional sliding

We measure the spatial and temporal behavior of the true contact area A along a rough spatially extended interface between two blocks in frictional contact. Upon the application of shear the onset of motion is preceded by a discrete sequence of cracklike precursors, which are initiated at shear levels that are well below the threshold for static friction. These precursors arrest well before traversing the entire interface. They systematically increase in length with the applied shear force and significantly redistribute the true contact area along the interface. Thus, when frictional sliding occurs, the initially uniform contact area along the interface has already evolved to one that is highly nonuniform in space.     

Amplitude reduction of parametric resonance by dynamic vibration absorber based on quadratic nonlinear coupling

The paper proposes an amplitude reduction method for parametric resonance with a new type of dynamic vibration absorber utilizing quadratic nonlinear coupling. A main system with asymmetric nonlinear restoring force and harmonic excitation causes parametric resonance in the system. In contrast with autoparametric vibration absorber, the natural frequency of the vibration absorber is tuned to be in the neighborhood of twice that of the main system. For such a vibration absorber, we investigate the effect on the amplitude reduction for a parametrically excited main system. Analytical results using the method of multiple scales show that the amplitude of parametric resonance is reduced by the effect of the vibration absorber. The experimental results by a simple apparatus indicate that the parametric resonance is stabilized by the effects of both vibration absorber and Coulomb friction of the main system. Moreover, numerical results considering the Coulomb friction of the main system show that the amplitude of parametric resonance becomes close to zero by the proposed vibration absorber.     

Influence of the normal force and contact geometry on the static force of friction of an oscillating sample

The paper is devoted to an experimental and theoretical investigation of the static friction force between a rapidly oscillating sample and a steel plate. The static frictional force is studied experimentally as function of the oscillating amplitude, the normal force and the contact geometry. A simplest model of tangent contact with a constant friction coefficient is proposed and shows a good agreement with experiment. The static friction force is proved to be a universal function of the ratio of the oscillation amplitude, the indentation depth and to the friction coefficient.     

计入固液界面作用的润滑热力学模型与分析

摩擦与润滑过程是典型的能量耗散过程, 在机理上与非平衡热力学中的熵增、耗散结构等理论颇有相似之处. 通过热力学分析可以对一些典型的摩擦磨损过程做出合理的机理揭示与推测. 本文利用热力学理论对典型的润滑过程进行了建模分析. 采用分离压模型表征和计入了微尺度下的固液界面作用, 揭示分析了润滑热力学模型与润滑状态Stribeck曲线的联系. 从分析计算结果来看, 润滑Stribeck曲线的摩擦系数最低点与系统热力学上的熵增率最低点具有相当好的对应关系, 而润滑状态从弹流润滑向薄膜润滑的转变过程, 可以用耗散结构理论加以机理解释. 文中的热力学模型和方法能够有效地体现出润滑过程中多物理要素跨尺度非线性耦合的作用, 对实际工程与实验有着重要的指导作用.     

Squeal and chatter phenomena generated in a mountain bike disc brake

This paper examines squeal and chatter phenomena generated experimentally in mountain bike disc brakes. There are two kinds of frictional self-excited vibrations in the bike disc brakes, called squeal with frequency of 1 kHz and chatter with frequency of 500 Hz. In order to reproduce the squeal and chatter, a bench test apparatus using an actual bike was set up to determine the associated frequency characteristics experimentally. The results show the frequencies to be independent of pad temperature and disc rotating speed. Squeal is shown to be in-plane vibration in the direction of the disc surface which is caused by the frictional characteristics having negative slope with respect to the relative velocity in the vibrating system, which includes brake unit, spokes and hub. Chatter is generated within a limited high temperature region. Again, it is frictional vibration in which the squeal and out-of-plane vibration of the disc due to Coulomb friction combine through the internal resonance relation between in-plane and out-of-plane nonlinear vibration caused by the temperature increase of the disc during braking.