Instability of thermoelastic contact for two half-planes sliding out-of-plane with contact resistance and frictional heating

Thermoelastic contact is known to show instabilities when the heat transmitted across the interface depends on the pressure, either because of a pressure-dependent thermal contact resistance R(p) or because of frictional heating due to the product of friction coefficient, speed, and pressure, fVp. Recently, the combined effect of pressure-dependent thermal contact resistance and frictional heating has been studied in the context of simple rod models or for a more realistic elastic conducting half-plane sliding against a rigid perfect conductor “wall”. Because R(p) introduces a non-linearity even in full contact, the “critical speed” for the uniform pressure solution to be unstable depends not just on material properties, and geometry, but also on the heat flux and on pressure.Here, the case of two different elastic and conducting half-planes is studied, and frictional heating is shown to produce significant effects on the stability boundaries with respect to the Zhang and Barber (J. Appl. Mech. 57 (1990) 365) corresponding case with no sliding. In particular, frictional heating makes instability possible for a larger range of prescribed temperature drop at the interface including, at sufficiently high speeds, the region of opposite sign of that giving instability in the corresponding static case. The effect of frictional heating is particularly relevant for one material combinations of the Zhang and Barber (J. Appl. Mech. 57 (1990) 365) classification (denominated class b here), as above a certain critical speed, the system is unstable regardless of temperature drop at the interface.Finally, if the system has a prescribed heat flow into one of the materials, the results are similar, except that frictional heating may also become a stabilizing effect, if the resistance function and the material properties satisfy a certain condition.     

Homogenization in non-linear dynamics due to frictional contact

This work is devoted to a study of the classical homogenization process and its influence on the behavior of a composite under non-linear dynamic loading due to contact and friction. First, the general problem of convergence of numerical models subjected to dynamic contact with friction loading is addressed. The use of a regularized friction law allows obtaining good convergence of such models. This study shows that for a dynamic contact with friction loading, the classical homogenization process, coupled with an homogenization of the frictional contact, enables replacing the entire heterogeneous model by a homogenized one. The dynamic part of the frictional contact must be homogenized by modifying the dynamic parameter of the friction law. Modification of the dynamic parameter of the friction law is function of the type and regime of instability. A calculation of a homogenized friction coefficient is presented in view to homogenizing the static part of the frictional contact when the friction coefficient is not constant over the contact surface. Finally matrix and heterogeneities stresses in the heterogeneous models are identified by using the relocalization process and a frictional contact dynamic analysis of a homogeneous model.     

接触问题的非线性互补－接触柔度法

基于势能互补原理，不引入额外松弛变量，对库仑摩擦定律未做预先线性化近似，提出了三维弹性摩擦接触问题的非线性互补－接触柔度法，收敛性和收敛速率得到了严格理论保证。数值实验结果表明方法可靠、有效、低存储量、高精度     

Response of frictional contact problems in thermo-rheologically complex structures

This paper presents a comprehensive computational model for predicting the nonlinear response of frictional viscoelastic contact systems under thermo-mechanical loading and experience geometrical nonlinearity. The nonlinear viscoelastic constitutive model is expressed by an integral form of a creep function, whose elastic and time-dependent properties change with stresses and temperatures. The thermo-viscoelastic behavior of the contacting bodies is assumed to follow a class of thermo-rheologically complex materials. An incremental-recursive formula for solving the nonlinear viscoelastic integral equation is derived. Such formula necessitates data storage only from the previous time step. The contact problem as a variational inequality constrained model is handled using the Lagrange multiplier method for exact satisfaction of the inequality contact constraints. A local nonlinear friction law is adopted to model friction at the contact interface. The material and geometrical nonlinearities are modeled in the framework of the total Lagrangian formulation. The developed model is verified using available benchmarks. The effectiveness and accuracy of the developed computational model is validated by solving two thermo-mechanical contact problems with different natures. Moreover, obtained results show that the mechanical properties and the class of thermo-rheological behavior of the contacting bodies as well as the coefficient of friction have significant effects on the contact response of nonlinear thermo-viscoelastic materials.     

NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS

A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.     

Analysis of FRP composites under frictional contact conditions

This work studies numerically the tribological behavior of fiber-reinforced plastics (FRP) under different frictional contact conditions, using a boundary element methodology. The formulation uses the Boundary Element Method (BEM) with an explicit approach for fundamental solutions evaluation, for computing the elastic influence coefficients. To enforce the contact constraints on the potential contact zone: Signorini’s contact conditions and an orthotropic law of friction, contact operators over the augmented Lagrangian are considered in the formulation. The methodology and the proposed algorithm are applied to study two types of glass FRP and two types of carbon FRP, with the same fiber volume fraction, under frictional contact. In these studies, it can be observed how the fiber orientation and sliding orientation affect the normal and tangential contact compliance, as well as the contact traction distribution. Furthermore, the formulation considers a micromechanics model for FRP that allows also to study the influence of fiber volume fraction on normal and tangential contact compliances.     

Two-dimensional sliding frictional contact of functionally graded materials

A multi-layered model for sliding frictional contact analysis of functionally graded materials (FGMs) with arbitrarily varying shear modulus under plane strain-state deformation has been developed. Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into several sub-layers and in each sub-layers the shear modulus is assumed to be a linear function while the Poisson's ratio is assumed to be a constant. In the contact area, it is assumed that the friction is one of Coulomb type. With this model the fundamental solutions for concentrated forces acting perpendicular and parallel to the FGMs layer surface are obtained. Then the sliding frictional contact problem of a functionally graded coated half-space is investigated. The transfer matrix method and Fourier integral transform technique are employed to cast the problem to a Cauchy singular integral equation. The contact stresses and contact area are calculated for various moving stamps by solving the equations numerically. The results show that appropriate gradual variation of the shear modulus can significantly alter the stresses in the contact zone.     

Instability and stochastic analyses of a pad-on-disc frictional system in moving interactions

In this paper, we establish a new nonlinear equation which is called the two-mode Korteweg–de Vries–Burgers equation (TMKdV–BE). The new equation describes the propagation of two different wave modes simultaneously. First, we introduce a new Cole–Hopf transformation needed for the modified simplified bilinear method to find necessary conditions on the solution of TMKdV–BE. Also, a finite series in terms of tanh–coth function is presented as an alternative method to study the solution of the equation. Finally, the reliability of the obtained results is discussed in the last section.     

Solving frictional contact problems by two aggregate-function-based algorithms

Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-function-based algorithms for solving complementarity problems are proposed. One is called the self-adjusting interior point algorithm, the other is called the aggregate function smoothing algorithm. Numerical experiment shows the efficiency of the proposed two algorithms. The project supported by the National Natural Science Foundation of China(10225212, 50178016, 10302007), the National Key Basic Research Special Foundation and the Ministry of Education of China The English text was polished by Ron Marshall.     

Optimal prestress of structures with frictional unilateral contact interfaces

Summary The problem of optimal prestress stabilization of elastic structures with frictional contact interfaces subject to static loads is studied in this paper. A linear elastic structure with given unilateral contact at frictional interfaces is considered. The prestressing control is modelled by the pin-load method. The static problem is formulated as a nonsymmetric variational inequality. The goal of the optimal control design is closing of the unilateral contact joints as well as minimization of the friction induced slips with a minimum effort. The resulting optimal control problem is nonsmooth and nonconvex, as it concerns the control of structures governed by variational inequalities. Appropriate techniques of nonsmooth analysis are used for its numerical solution. Effective computer realization and integration into existing finite element software is facilitated by appropriate static condensation techniques, which are outlined in the paper. Numerical examples illustrate the theory.     

The thermoelastic Aldo contact model with frictional heating

In the study of the essential features of thermoelastic contact, Comninou and Dundurs (J. Therm. Stresses 3 (1980) 427) devised a simplified model, the so-called “Aldo model”, where the full 3D body is replaced by a large number of thin rods normal to the interface and insulated between each other, and the system was further reduced to 2 rods by Barber's Conjecture (ASME J. Appl. Mech. 48 (1981) 555). They studied in particular the case of heat flux at the interface driven by temperature differences of the bodies, and opposed by a contact resistance, finding possible multiple and history dependent solutions, depending on the imposed temperature differences.The Aldo model is here extended to include the presence of frictional heating. It is found that the number of solutions of the problem is still always odd, and Barber's graphical construction and the stability analysis of the previous case with no frictional heating can be extended. For any given imposed temperature difference, a critical speed is found for which the uniform pressure solution becomes non-unique and/or unstable. For one direction of the temperature difference, the uniform pressure solution is non-unique before it becomes unstable. When multiple solutions occur, outermost solutions (those involving only one rod in contact) are always stable.A full numerical analysis has been performed to explore the transient behaviour of the system, in the case of two rods of different size. In the general case of N rods, Barber's conjecture is shown to hold since there can only be two stable states for all the rods, and the reduction to two rods is always possible, a posteriori.     

三维接触问题的非光滑算法

给出了一种非光滑算法直接用于求解三维摩擦接触问题的不可微非线性互补模型,不再对模型进行光滑化处理,使算法更加简单。文中对非光滑算法的收敛性给出了严格的数学证明,数值实验表明该算法列式简单,但与光滑化算法同样有效。     

Block pivot methods for solving frictional contact problems

Based on elementary group theory, the block pivot methods for solving two-dimensional elastic frictional contact problems are presented in this paper. It is proved that the algorithms converge within a finite number of steps when the friction coefficient is “relative small”. Unlike most mathematical programming methods for contact problems, the block pivot methods permit multiple exchanges of basic and nonbasic variables. The project supported by the National Natural Science Foundation of China     

3D frictional contact of anisotropic solids using BEM

A numerical method to study three-dimensional (3D) contact problems in solids with anisotropic elastic behavior is developed in this work. This formulation is based on the Boundary Element Method (BEM) for computing the elastic influence coefficients and on projection functions over the augmented Lagrangian for contact restrictions fulfillment. The constitutive equations of the potential contact zone are Signorini’s contact conditions and Coulomb’s law of friction. The formulation uses a recently introduced explicit approach for fundamental solutions evaluation, which are valid for general anisotropic behavior meanwhile mathematical degeneracies are allowed. The accuracy and robustness of the proposed method is illustrated by solving some examples previously presented in the literature. This approach is further applied to study the influence of solids anisotropy on the contact problem.     

Slip-shakedown analysis of a system of circular beams in frictional contact

In automotive components, the cumulative microslip phenomenon is often observed for engine assemblies. This phenomenon results in an accumulation of the relative slips in a preferred tangential direction on the contact interface of two solids under cyclic loadings. A significant relative displacement may occur and leads to the assembly failure. In particular, a global rotation of the bearing shell may result from this mechanism of cumulated slips in conrod big end systems. To discuss this rotation problem, a model of two circular beams in frictional contact and submitted to a periodical rotating load is considered here. The aim is to give some simplified estimates of the critical rotation load based on a slip-shakedown analysis. The discussion holds for Tresca friction and can be extended to Coulomb friction under the assumption of small coupling. The static and kinematic slip-shakedown approaches are discussed. The obtained analytical results are shown to be in agreement with the finite element computations.     

Interaction of thermal contact resistance and frictional heating in thermoelastic instability

Thermoelastic contact problems can posess non-unique and/or unstable steady-state solutions if there is frictional heating or if there is a pressure-dependent thermal contact resistance at the interface. These two effects have been extensively studied in isolation, but their possible interaction has never been investigated. In this paper, we consider an idealized problem in which a thermoelastic rod slides against a rigid plane with both frictional heating and a contact resistance. For sufficiently low sliding speeds, the results are qualitatively similar to those with no sliding. In particular, there is always an odd number of steady-state solutions; if the steady-state is unique it is stable and if it is non-unique, stable and unstable solutions alternate, with the outlying solutions being stable. However, we identify a sliding speed V₀ above which the number of steady states is always even (including zero, implying possible non-existence of a steady-state) and again stable and unstable states alternate. A parallel numerical study shows that for V>V₀ there are some initial conditions from which the contact pressure grows without limit in time, whereas for V<V₀ the system will always tend to one of the stable steady states.     

The least-squares meshfree method for rigid-plasticity with frictional contact

The least-squares meshfree method (LSMFM) for rigid-plasticity based on J₂-flow rule and infinitesimal theory is proposed. In the least-squares formulation the squared residuals of the constitutive and equilibrium equations are minimized. Those residuals are represented in a form of first-order differential system using the velocity and stress components as nodal unknowns and thus the proposed formulation is a mixed-type method. Also the penalty scheme for the enforcement of the boundary and frictional contact conditions is devised and the reshaping of nodal supports is introduced to avoid the difficulties due to the severe local deformation near the contact interface. The proposed method does not require any structure of extrinsic cells for the construction of shape functions, the treatment of incompressibility, the integration of variational formulation and the reconstruction of approximation. Through some numerical examples of metal forming processes, the validity and effectiveness of the method are discussed.     

The equilibrium state of a structure subject to frictional contact

A structure in frictional contact subject to static loads has not, in general, a unique static equilibrium state. This is because the state, displacements and contact forces, depend on the load history of the structure.In cases where the exact load history is not known it would be of interest to find a state that is in some sense likely and define this as the equilibrium state. In this paper, it is assumed that the state with the smallest potential energy is the most likely one. The implication of this definition of likely state is analysed and shows that the resulting problem basically can be seen as a generalization of the frictionless contact problem to structures where no frictionless state is possible, i.e. structures where non-zero friction forces are necessary to satisfy force equilibrium.The results of several numerical experiments are given. The structures in the experiments are trusses and structures modelled by the finite element method. Both a sequential quadratic programming method and an enumeration method are used to solve the likely-state problem.     

General contact boundary conditions and the analysis of frictional systems

General factional systems, i.e. solid bodies interacting by contact forces, are investigated and variational formulations are presented. Following recent ideas, general contact boundary conditions are formulated using the notions of subdifferentials and generalized gradients.