Frictional elastic contact with periodic loading

Quasi-static frictional contact problems for bodies of fairly general profile that can be represented as half planes can be solved using an extension of the methods of Ciavarella and Jäger. Here we consider the tangential traction distributions developed when such systems are subjected to loading that varies periodically in time. It is shown that the system reaches a steady state after the first loading cycle. In this state, part of the contact area (the permanent stick zone) experiences no further slip, whereas other points may experience periods of stick, slip and/or separation. We demonstrate that the extent of the permanent stick zone depends only on the periodic loading cycle and is independent of the initial conditions or of any initial transient loading phase. The exact traction distribution in this zone does depend on these factors, but the resultant of these tractions at any instant in the cycle does not. The tractions and slip velocities at all points outside the permanent stick zone are also independent of initial conditions, confirming an earlier conjecture that the frictional energy dissipation per cycle in such systems depends only on the periodic loading cycle. We also show that these parameters remain unchanged if the loading cycle is changed by a time-independent tangential force, provided this is not so large as to precipitate a period of gross slip (sliding).     

Modeling of fault slip in earthquakes

The total energy release and the length of interfacial crack growth associated with an earthquake are analyzed based on a horizontal one-dimensional semi-infinite elastic model of fault of uniform rectangular cross section between two rigid rocks. Initial interfacial cracks of finite length are assumed to exist on the top and bottom surfaces of the fault. A compressive force is applied at the end of the fault and uniform frictional forces are exerted on the crack surfaces. When the compressive force reaches a critical value, fault slip against the surrounding rock occurs and the frictional force changes from a static value to a dynamic value within a short time interval. The change in frictional force generates longitudinal wave propagation in the fault. As a result, a large amount of energy release takes place. The total amount of energy release and the length of crack growth can be analyzed based on the balance of energy.     

Effect of phase on the frictional dissipation in systems subjected to harmonically varying loads

When contacting elastic systems are subjected to periodic loading, frictional slip occurs resulting in energy dissipation. Here we investigate the effect of the relative phase of harmonically varying tangential and normal loads on the frictional dissipation in a very simple uncoupled frictional system. We demonstrate that this effect is substantial when the system experiences periods of separation, but more modest when contact is continuous. The maximum dissipation occurs when the normal and tangential loads are approximately π/2 out of phase.     

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

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

阶跃载荷作用下弹塑性悬臂梁的变形机制与能量耗散

基于大挠度动力控制方程,应用有限差分离散求解,研究了阶跃载荷作用下弹塑性悬臂梁的动力行为。通过对动力响应早期内力、变形以及能量分布规律的分析,考察了悬臂梁的弹塑性响应模式和变形机制,并与已有的刚塑性分析进行了系统的比较。数值计算表明,阶跃载荷的不同幅值使得梁的响应模式存在较大差异,弹塑性分析肯定了刚塑性理论在处理中载情形的准确性,同时也指出了其在处理低载和高载情形时的缺陷。通过与小变形理论计算结果的比较,指出了考虑大变形效应的必要性,为今后的大变形刚塑性动力分析提出了建议。     

Frictional complete contacts between elastically similar bodies subject to normal and shear load

The problem of trapezium shaped punches pressed into a frictional, elastically similar half-plane and subject to sequential normal and shear loading is studied. The initial mix of stick, slip and separation regions is described, together with the steady state response when the shearing force is cycled. Conditions for full stick are established.     

Breakage mechanics—Part II: Modelling granular materials

The compression of granular materials has been traditionally modelled with the limitations of classical elasto-plasticity. The energy was implicitly assumed to dissipate from the frictional interaction of particles. However, the fact that brittle granular materials crush suggests that energy must also be dissipated from the fracturing of the grains, as in fracture mechanics. The concept of breakage as a thermomechanical internal variable was introduced in Part I [Einav, I., 2006. Breakage mechanics—Part I: theory. J. Mech. Phys. Solids 00,000-000] to describe the fracturing mechanisms. The theory allows to treat ideal theoretical materials that undergo dissipation purely from breakage with no other mechanism allowed for the energy consumption. However, as accounted for in elasto-plasticity, dissipation must also occur from the frictional rearrangement of grains. The combination of the two dissipative mechanisms of breakage and plasticity must therefore be investigated, as we do in this paper. Those two mechanisms are generally coupled, in the sense that one inevitably appears when the other develops. Plastic dissipation emerges as a by-product of breakage dissipation because after grains crush, local rearrangement must occur. This scenario may be termed an ‘active breakage mechanism’, and typifies compression deformations. In shear the plastic dissipation is dominant but breakage appears inevitably from grains abrasion. This scenario may be called a ‘passive breakage mechanism’. Based on the coupling assumption, models are developed for granular materials. In particular, we show that in compression isotropic hardening of sands may appear without involving plastic strains, i.e., independent of frictional dissipation. This interpretation of hardening is different from the one used in classical critical state soil mechanics. However, frictional dissipation leads to plastic straining that are necessary for the models to be predictive in unloading.     

A model for frictional slip in woven fabrics

A continuum model for frictional slip of the yarns of a plain-weave fabric is presented. The model is based on the assumption that the weave is composed of two families of continuously distributed yarns constrained at all times to occupy a common evolving surface in three-dimensional space. The two families may slide relative to one another on the surface, subject to their respective equations of motion, fiber constitutive equations, and frictional slip rules. The theory is intended for the quantitative analysis of deformation, slip and energy dissipation during a ballistic impact event. To cite this article: B. Nadler, D.J. Steigmann, C. R. Mecanique 331 (2003).     

微液滴在PDMS软基体表面的动态摩擦行为研究

通过固液界面摩擦力测试装置研究了微液滴在PDMS软基体表面运动时的动态摩擦学行为,并对微液滴体积、滑动速度及软基体力学性能对固液界面动态摩擦行为的影响进行了分析. 结果表明:微液滴在软基体表面运动时表现出最大静摩擦力和动态摩擦力. 最大静摩擦力与微液滴黏度和速度梯度呈正比,动态摩擦力与微液滴体积、滑动速度和基体力学性能有关. 随着微液滴体积的增加,三相接触线长度增加,动态摩擦力增加;随着相对滑动速度增加,三相接触线长度及接触角滞后增加,动态摩擦力增加;随着软基体弹性模量降低,固液界面黏附力增加,固液界面运动能量耗散增加,动态摩擦力增加. 研究结果可为PDMS软基体表面微液滴的精确驱动和运动参数优化提供理论指导,也可进一步丰富固液界面摩擦理论.      

A Dynamic-induced Direct-shear Model for Dynamic Triggering of Frictional Slip on Simulated Granular Gouges

This study presents a dynamic-induced direct-shear model to investigate the dynamic triggering of frictional slip on simulated granular gouges. An incident P-wave is generated as a shear load and a normal stress is constantly applied on the gouge layer. The shear stress accumulates in the incident stage and the frictional slip occurs in the slip stage without the effect of the reflected wave. The experimental results show a non-uniform shear stress distribution along the gouge layer, which may be induced by a shear load induced torque and by normal stress vibration along the layer. The shear stress at the trailing edge strongly affects the frictional slip along the P-wave loading direction, while the rebound stress at the leading edge propagates along the opposite direction. The frictional slip is triggered when the maximum shear stress at the trailing edge reaches a critical value. The normal stress influences the maximum shear stress at the trailing edge, the maximum slip displacement and the slip velocity. The advantages and the limitations of this model are discussed at the end.     

Frictional energy dissipation in materials containing cracks

Kachanov's simplified model of microcrack interaction is applied to an investigation of the behaviour of a cracked body under predominantly compressive periodic loading, so that the cracks experience periods of closure and slip, with frictional dissipation. The model is shown to be equivalent to a discrete elastic frictional system with each crack representing one node. Theorems and algorithms from such systems are applied to determine the conditions under which the system shakes down to a state with no slip and hence no energy dissipation in friction. For conditions not too far beyond the shakedown state, the dissipation is significantly affected by the initial conditions, but with larger oscillating loads, it becomes a unique and increasing function of load amplitude. The effect of crack interaction is assessed by comparison with an uncoupled model, for which the dissipation is obtained as a summation of closed form expressions over the crack population. For small numbers of cracks, the results are significantly dependent on the randomly chosen crack locations and sizes, but with larger populations, a statistically significant decrease in dissipation is observed with increasing interaction terms.     

Transition waves in bistable structures. I. Delocalization of damage

We consider chains of dimensionless masses connected by breakable bistable links. A non-monotonic piecewise linear constitutive relation for each link consists of two stable branches separated by a gap of zero resistance. Mechanically, this model can be envisioned as a ”twin-element” structure which consists of two links (rods or strands) of different lengths joined by the ends. The longer link does not resist to the loading until the shorter link breaks. We call this construction the waiting link structure. We show that the chain of such strongly non-linear elements has an increased in-the-large stability under extension in comparison with a conventional chain, and can absorb a large amount of energy. This is achieved by two reasons. One is an increase of dissipation in the form of high-frequency waves transferring the mechanical energy to heat; this is a manifestation of the inner instabilities of the bonds. The other is delocalization of the damage of the chain. The increased stability is a consequence of the distribution of a partial damage over a large volume of the body instead of its localization, as in the case of a single neck formation in a conventional chain. We optimize parameters of the structure in order to improve its resistance to a slow loading and show that it can be increased significantly by delocalizing a damage process. In particular, we show that the dissipation is a function of the gap between the stable branches and find an optimal gap corresponding to maximum energy consumption under quasi-static extension. The results of numerical simulations of the dynamic behavior of bistable chains show that these chains can withstand without breaking the force which is several times larger than the force sustained by a conventional chain. The formulation and results are also related to the modelling of compressive destruction of a porous material or a frame construction which can be described by a two-branched diagram with a large gap between the branches. We also consider an extension of the model to multi-link chain that could imitate plastic behavior of material.     

Peristaltic Hemodynamic Flow of Couple-Stress Fluids Through a Porous Medium with Slip Effect

The present investigation deals with a theoretical study of the peristaltic hemodynamic flow of couple-stress fluids through a porous medium under the influence of wall slip condition. This study is motivated towards the physiological flow of blood in the micro-circulatory system, by taking account of the particle size effect. Reynolds number is small enough and the wavelength to diameter ratio is large enough to negate inertial effects. Analytical solutions for axial velocity, pressure gradient, frictional force, stream function and mechanical efficiency are obtained. Effects of different physical parameters reflecting couple-stress parameter, permeability parameter, slip parameter, as well as amplitude ratio on pumping characteristics and frictional force, streamlines pattern and trapping of peristaltic flow pattern are studied with particular emphasis. The computational results are presented in graphical form. This study puts forward an important observation that pressure reduces by increasing the magnitude of couple-stress parameter, permeability parameter, slip parameter, whereas it enhances by increasing the amplitude ratio.     

Analysis of Nonequilibrium Dissociation Kinetics in a Boundary Layer Using Various Reaction Models

Values of the nonequilibrium macroscopic reaction rate for a nonisothermal boundary layer of a monatomic diluent gas are calculated using a number of models for thermal dissociation of diatomic molecules — anharmonic Morse oscillators. Analysis is performed for conditions where the diffusive transfer of excited molecules has a significant effect on the population of their upper vibrational levels, which does not only amount to change in vibrational temperature. Under the joint influence of diffusive transfer of molecules, vibrational exchanges, and reactions involving vibrationally excited particles, the local vibrational distribution functions are substantially nonequilibrium. The kinetic models considered take into account the possible contribution of the energy of molecular translational and rotational degrees of freedom to the energy required to overcome the reaction threshold. The effect of multiquantum vibrational—translational exchanges on the distribution of dissociating molecules in their upper vibrational levels is taken into account approximately.     

Physics-based modeling for partial slip behavior of spherical contacts

A physics-based modeling approach for partial slip behavior of a spherical contact is proposed. In this approach, elastic and elastic–plastic normal preload and preload-dependent friction coefficient models are integrated into the Cattaneo–Mindlin partial slip solution. Partial slip responses to cyclic tangential loading (fretting loops) obtained by this approach are favorably compared with experiments and finite element results from the literature. In addition to load-deformation curves, tangential stiffness of the contact and energy dissipation per fretting cycle predictions of the models are also provided. Finally, the critical assumptions of elastically similar bodies, smooth contact surface and negligible adhesion, and limitations of this physics-based modeling approach are discussed.     

Impact comminution of solids due to local kinetic energy of high shear strain rate: I. Continuum theory and turbulence analogy

The modeling of high velocity impact into brittle or quasibrittle solids is hampered by the unavailability of a constitutive model capturing the effects of material comminution into very fine particles. The present objective is to develop such a model, usable in finite element programs. The comminution at very high strain rates can dissipate a large portion of the kinetic energy of an impacting missile. The spatial derivative of the energy dissipated by comminution gives a force resisting the penetration, which is superposed on the nodal forces obtained from the static constitutive model in a finite element program. The present theory is inspired partly by Grady's model for expansive comminution due to explosion inside a hollow sphere, and partly by analogy with turbulence. In high velocity turbulent flow, the energy dissipation rate gets enhanced by the formation of micro-vortices (eddies) which dissipate energy by viscous shear stress. Similarly, here it is assumed that the energy dissipation at fast deformation of a confined solid gets enhanced by the release of kinetic energy of the motion associated with a high-rate shear strain of forming particles. For simplicity, the shape of these particles in the plane of maximum shear rate is considered to be regular hexagons. The particle sizes are assumed to be distributed according to the Schuhmann power law. The condition that the rate of release of the local kinetic energy must be equal to the interface fracture energy yields a relation between the particle size, the shear strain rate, the fracture energy and the mass density. As one experimental justification, the present theory agrees with Grady's empirical observation that, in impact events, the average particle size is proportional to the (−2/3) power of the shear strain rate. The main characteristic of the comminution process is a dimensionless number B_a (Eq. (37)) representing the ratio of the local kinetic energy of shear strain rate to the maximum possible strain energy that can be stored in the same volume of material. It is shown that the kinetic energy release is proportional to the (2/3)-power of the shear strain rate, and that the dynamic comminution creates an apparent material viscosity inversely proportional to the (1/3)-power of that rate. After comminution, the interface fracture energy takes the role of interface friction, and it is pointed out that if the friction depends on the slip rate the aforementioned exponents would change. The effect of dynamic comminution can simply be taken into account by introducing the apparent viscosity into the material constitutive model, which is what is implemented in the paper that follows.     

Elastic beam on a rigid frictional foundation under monotonic and cyclic loading

An elastic beam resting on a frictional foundation and loaded by the concentrated force or moment applied at its tip is considered. The evolution of slip zones along the beam is discussed for both monotonie and cyclic loading. It is shown that an infinite number of slip zones develop and their propagation satisfies in some cases a self-similarity property. Transient hysteretic effects under cyclic loading are discussed. The closed form analytical solution is presented for the elastic friction model in the case of monotonie loading.     

Nonlinear multimodal model for TLD of irregular tank geometry and small fluid depth

Tuned liquid dampers (TLDs) utilize sloshing fluid to absorb and dissipate structural vibrational energy. TLDs of irregular or complex tank geometry may be required in practice to avoid tank interference with fixed structural or mechanical components. The literature offers few analytical models to predict the response of this type of TLD, particularly when the fluid depth is small. In this paper, a multimodal model is developed utilizing a Boussinesq-type modal theory which is valid for small TLD fluid depths. The Bateman–Luke variational principle is employed to develop a system of coupled nonlinear ordinary differential equations which describe the fluid response when the tank is subjected to base excitation. Energy dissipation is incorporated into the model from the inclusion of damping screens. The fluid model is used to describe the response of a 2D structure–TLD system when the structure is subjected to external loading and the TLD tank geometry is irregular.Shake table experiments are conducted on a rectangular and chamfered tank subjected to unidirectional base excitation. Comparisons of the experimental and predicted sloshing forces and energy dissipation per cycle indicate that the model is able to predict the fluid response at fluid depth ratios greater than h/L=0.10. Next, structure–TLD system tests are conducted and it is found that the model can predict the structural and TLD responses. The simulated and experimental results show that the TLD tank transfers energy between orthogonal structural sway modes.     

Direct simulation of multiphase flows with modeling of dynamic interface contact angle

We describe a modeling technique for dynamic contact angle between a phase interface and a solid wall using a generalized Navier boundary condition in the context of a front-tracking-based multiphase method. The contact line motion is determined by the generalized Navier slip boundary condition in order to eliminate the infinite shear stress at the contact line. Applying this slip boundary condition only to the interface movement with various slip ratios shows good agreement with experimental results compared to allowing full fluid slip along the solid surface. The interface slip model performs well on grid convergence tests using both the slip ratio and slip length models. A detailed energy analysis was performed to identify changes in kinetic, surface, and potential energies as well as viscous and contact line dissipation with time. A friction coefficient for contact line dissipation was obtained based on the other computed energy terms. Each energy term and the friction coefficient were compared for different grid resolutions. The effect of varying the slip ratio as well as the contact angle distribution versus contact line speed was analyzed. The behavior of drop impact on a solid wall with different advancing and receding angles was investigated. Finally, the proposed dynamic contact model was extended to three dimensions for large-scale parallel calculations. The impact of a droplet on a solid cylinder was simulated to demonstrate the capabilities of the proposing formulation on general solid structures. Widely different contact angles were tested and showed distinctive characteristic behavior clearly.     

Tractive performance of a driven rigid wheel on soft ground based on the analysis of soil-wheel interaction

A new analytical method has been presented to predict the tractive performance of a rigid wheel running on soft ground. The resultant stress of the normal stress and the shear resistance applied around the peripherical contact part of the rigid wheel should be calculated by use of the dynamic pressure-sinkage curve measured from the plate loading and unloading test, considering the rolling locus of the wheel in the direction of the external resultant force of the effective driving force and the axle load. The effective driving force could be calculated as the difference of the driving force, i.e. the integration of shear resistance and the locomotion resistance calculated from the total amount of sinkage. As a result, the analytical relations between the driving force, the effective driving force and the slip ratio, the amount of sinkage and the slip ratio, the amount of eccentricity of resultant force and the slip ratio, and the entry angle, the exit angle and the slip ratio could be verified experimentally.