Existence and numerical approximation of a solution to frictional contact problem for electro-elastic materials

In this paper,a frictional contact problem between an electro-elastic body and an electrically conductive foundation is studied.The contact is modeled by normal compliance with finite penetration and a version of Coulomb’s law of dry friction in which the coefficient of friction depends on the slip.In addition,the effects of the electrical conductivity of the foundation are taken into account.This model leads to a coupled system of the quasi-variational inequality of the elliptic type for the displacement and the nonlinear variational equation for the electric potential.The existence of a weak solution is proved by using an abstract result for elliptic variational inequalities and a fixed point argument.Then,a finite element approximation of the problem is presented.Under some regularity conditions,an optimal order error estimate of the approximate solution is derived.Finally,a successive iteration technique is used to solve the problem numerically and a convergence result is established.

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.

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.

A NEW METHOD FOR SOLUTION OF 3D ELASTIC-PLASTIC FRICTIONAL CONTACT PROBLEMS

The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.

Propagation of harmonic waves through micro gap with consideration of frictional contact

Transmission of elastic waves through a micro gap between two solids with consideration of frictional contact is investigated. By using the Fourier analysis technique and the corrective solution method, the nonlinear boundary problem is reduced to a set of algebraic equations. Numerical results exhibit the locations and extents of separation, slip, and stick zones, the interface tractions, and the energy partition. The effects of gap width, frictional coefficients, and the incident angle on the wave transmission are discussed in detail. The results show that higher harmonics are generated due to the local contact/slip at the interface.

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-functionbased 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.

Numerical Analysis of a Sliding Frictional Contact Problem with Normal Compliance and Unilateral Contact

This paper represents a continuation of [1] and [2]. Here, we consider the numerical analysis of a non-trivial frictional contact problem in a form of a system of evolution nonlinear partial differential equations. The model describes the equilibrium of a viscoelastic body in sliding contact with a moving foundation. The contact is modeled with a multivalued normal compliance condition with memory term restricted by a unilateral constraint and is associated with a sliding version of Coulomb’s law of dry friction. After a description of the model and some assumptions, we derive a variational formulation of the problem, which consists of a system coupling a variational inequality for the displacement field and a nonlinear equation for the stress field. Then, we introduce a fully discrete scheme for the numerical approximation of the sliding contact problem. Under certain solution regularity assumptions, we derive an optimal order error estimate and we provide numerical validation of this result by considering some numerical simulations in the study of a two-dimensional problem.

FRICTIONAL CONTACT MULTIPOLE-BEM AND 3-D ANALYSIS OF SCREWPAIRS

The 3-D traction field in the pressure screw-pair of a 3 500 heavy and mediumplate mill press down system is successfully calculated by applying the 3-D frictional contactmultipole-BEM and the corresponding program that has been developed. The computing results show themedium diameter orientation is unreliable, especially under the interference of an outer forcecouple. Under such working conditions, the circumferential traction distribution on the screw teethis extremely uneven, which is the main reason for the destruction and short life time ofscrew-pairs. When utilizing the same precision (the relative tolerance is 10X10^(-5)), themultipole-BEM uses almost the same CPU time as used by the FEM, but the needed computer memory sizeis only one eightieth of that needed by the FEM (10 MB vs. 800 MB). The multipole-BEM is well suitedfor computing large-scale engineering problems.

Frictional contact problem for elastic-viscoplastic materials with thermal effect

A dynamic contact problem for elastic-viscoplastic materials with thermal effects is investigated. The contact is bilateral, and the friction is modeled with Tresca＇s friction law with heat exchange. A variational formulation of the model is derived, and the existence of a unique weak solution is proved. The proofs are based on the classical result of nonlinear first order evolution inequalities, the equations with monotone operators, and the fixed point arguments. Finally, the continuous dependence of the solution on the friction yield limit is studied.

F0RMULATION OF BOUNDARY ELEMENT-LINEAR COMPLEMENTARY EQUATION FOR THE FRICTIONAL ELASTIC CONTACT PROBLEMS

Boundary element-linear complementary equations are formulated to solve elastic contact problems with Coulomb frictions. It is also a new attempt to solve free boundary problems in solid mechanics by means of boundary element-mathematical programming techniques.

EFFECT OF FRICTIONAL FORCE ON GEARING CONTACT FATIGUE STRENGTH

The model for computing frictional coefficient between two teeth faces at the state of mixed elastohydrodynamic lubrication is established. And then more than 80 sets of numerical calculations and six sets of disc fatigue tests are completed. The results show that when the film thickness ratio λ 〈1.6, frictional coefficient μ is drastically decreased as λ. rises; Thereafter it decreases smoothly until λ=4.5. When λ〉4.5, however, it goes up again with λ, which indicates that the excessive film thickness ratio will deteriorate gearing contact fatigue strength. At the end, the formulae for determining the frictional coefficients are formed.

Influence of contact curvature on frictional energy dissipation under varying tangential loads

This study investigates the effect of contact surface curvatures on the friction response under varying tangential loadings using a finite element(FE)model.The results showed that the geometry of the surface influences the contact force at the interface and reduces the friction effect through an unsteady distribution of the contact force.The relationship between the friction effect,excitation,and contact surface shape was also examined,revealing a linear inverse relationship between the friction and curvature.The findings provide a comprehensive understanding of the frictional interactions between elastic bodies and highlight the role of curvature as a design parameter for regulating the friction effect.

CONTACT FRICTION ANALYSIS AND STRESS OSCILLATION SUPPRESSION WITH A SIMPLE INTERFACE ELEMENT

A simple interface element for analyzing contact friction problems is developed. Taking nodal displacements and contact stresses as unknowns, this element can simulate frictional slippage, decoupling and re-bonding of two bodies initially mating or having gaps at a common interface. The method is based on the Finite Element Method and load incremental theory. The geometric and static constraint conditions on contact surfaces are treated as additional conditions and are included in stiffness equations. This simple element has the advantages of easy implementation into standard finite element programs and fast speed for convergence as well as high accuracy for stress distribution in interface. Undesirable stress oscillations are also investigated whenever large stress gradients exist over the contact surfaces. Exact integration or the conventional Gauss integration scheme used to evaluate the interpolation function matrix of the interface element is found to be the source of the oscillations. Eigenmode analysis demonstrates that the stress behavior of an interface element can be improved by using the Newton-Cotes integration scheme. Finally, the test example of a strip footing problem is presented.

Analysis of composite material interface crack face contact and friction effects using a new node-pairs contact algorithm

A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.

ON LARGE DEFORMATION UNILATERAL CONTACT PROBLEM WITH FRICTION(Ⅰ)—INCREMENTAL VARIATIONAL EQUATIONS

Based on the theory and technique of nonlinear geometric field theory of continuum, a more general incremental variational equation for elastic and plastic large deformation in co-moving coordinate is established in this paper. An expression for two and three-ditnensicnal continua is derived, and the incremental variational equation for large deformation of changing boundary contact and the variational inequality in rate form tire obtained, which provides the theoretical basis for the computation of elastic-plastic large deformation contact problem with friction.

SOLVING CONTACT PROBLEM WITH FRICTION BY A NEW FAST BOUNDARY ELEMENT METHOD

The formulation of boundary element method for handling contact problems with friction and the technique for high speed contact analysis are presented. This formulation is based on the idea of modifying the length of contact elements without altering the total number of elements. The high precision of solution and high speed analysis are verified according to the results of conventional method and analysis method.

Multiphysical contact of functionally graded magneto-electroelastic material with imperfect bonding interface

Three-dimensional(3D)frictional contact model of functionally graded magneto-electro-elastic(FGMEE)material with a conducting spherical punch under electromagnetic fields is presented.Two types of imperfect bonding interface of layers,dislocation-like interface and force-like interface,are considered.Frequency response functions(FRFs)for multilayered MEE material with imperfect interface subjected to unit mechanical,electric,and magnetic loads are derived.The FRFs are used with the semi-analytical method(SAM)to solve present multiphysical contact problem.The present model is verified by comparing with literatures and the finite element method(FEM)and used to study the contact problem of FGMEE film imperfectly bonded on homogenous MEE half-space under electromagnetic fields.Parametric studies are conducted to reveal the effects of imperfect interfaces and also film properties including gradient index and thickness.

ANALYSIS OF AN ELASTO-PIEZOELECTRIC SYSTEM OF HEMIVARIATIONAL INEQUALITIES WITH THERMAL EFFECTS

In this paper we prove the existence and uniqueness of a weak solution for a dynamic electo-viscoetastic problem that describes a contact between a body and a foundation. We assume the body is made from thermoviscoelastic material and consider nonmonotone boundary conditions for the contact. We use recent results from the theory of hemivariational inequalities and the fixed point theory.

MESHFREE METHOD BASED ON POINT COLLOCATION FOR METAL FORMING SIMULATION

A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM （finite element method） and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.

ANALYSIS OF DYNAMIC INSTABILITY OF INTERFACIAL SLIP WAVES BASED ON THE SURFACE IMPEDANCE TENSOR

A new method relying on the Stroh formulism and the theory of the surface impedance tensor was developed to investigate the dynamic instability of interfacial slip waves.The concept of the surface impedance tensor was extended to the case where the wave speed is of a complex value,and the boundary conditions at the frictionally contacting interface were expressed by the surface impedance tensor.Then the boundary value problem was transformed to searching for zeroes of a complex polynomial in the unit circle.As an example,the steady frictional sliding of an elastic half-space in contact with a rigid flat surface was considered in details.A quartic complex characteristic equation was derived and its solution behavior in the unit circle was discussed.An explicit expression for the instability condition of the interfacial slip waves was presented.