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.

Frictional contact algorithm study on the numerical simulation of large deformations in deep soft rock tunnels

There exist three types of nonlinear problems in large deformation processes of deep softrock engineering, i.e., nonlin- earity caused by material, geometrical and contact boundary. In this paper, the numerical method to tackle the nonlinear eontact and large deformation problem in A Software on Large Deformation Analysis for Soft Rock Engineering at Great Depth was presented. In the software, based on Lagrange multiplier method and Coulomb friction law, kinematic constraints on contact boundaries were introduced in functional function, and the finite element equations was established for two incremental large deformation analyses models, polar decomposition model and additive decomposition model. For every incremental loading step, by searching for the contact points in the potential contact interfaces （the excavation boundaries）, the Lagrange multipliers, i.e., contact forces are cal- culated iteratively by Gauss-Seidel method, and justified through satisfy the inequalities of static constraint on contact boundaries. With the software, large deformation and frictional contact of a transport roadway were analyzed numerically by the two models. The numerical examples demonstrated the efficiency of the method used in the software.

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.

On the correction of the boundary deficiency in SPH for the frictional contact simulation

Smoothed particle hydrodynamics(SPH)is a mesh-free method which is powerful for large deformation computation of soils.However,the algorithm for the simulation of frictional contact which is common in geotechnical engineering is still quite immature due to the boundary deficiency.In this study,the cause of boundary deficiency in the SPH simulation for frictional contact is analysed.Then,based on mathematical derivation,the method to correct boundary deficiency related to frictional contact is discussed theoretically,where the frictional contact algorithm is established by dividing the computational domain into several subdomains according to the existing contact boundaries and by using contact forces as bridges of these subdomains to fulfil problem solving,and the value of correction coefficient is obtained by comparing the SPH outcome of the contact particles with that calculated through Newton’s second law of motion.At the same time,from the perspective of numerical computation,an optimized value for the correction coefficient is proposed,and a thorough investigation is performed on the cubic spline kernel function and quintic spline kernel function,whose correction coefficients are found to be 2.0 and[2.0,2.16],respectively.Finally,numerical tests are carried out to verify the proposed method.The outcome of the study is helpful to providing theoretical support for the research of frictional contact simulation within the framework of SPH.

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.

Failure of Semi-infinite Beams of Variable Thickness and Curvature on Elastic Foundation under Contact Loading

In early winter it is usual, in cold regions, that ice features approach offshore structures, like offshore platforms, impacting them, in a slow process of constant deformation build up. Interaction follows, in many cases, up to the point where ice-failure caused by bending fracture takes place. This supposes very large contact forces that the structure has to resist. Therefore, quantification of these efforts is of vital importance to the structural design of platforms. In several designs, these platforms are constructed with inclined walls so as to cause ice to fail in a flex-compression mode. In such a case the ice feature is analyzed as a beam constituted of a linear elastic material in brittle state with constant ice thickness. The simplification renders the problem solvable in a close form. However, this hypothesis goes against field observations. Marine currents action, wind and the sequence of contacts among features lead to thickness variations. Here this factor is addressed in the construction of a model, for harmonic forms of variation of thickness profile, and the accompanying curvature variations, whose solution determines field variables used to address the failure question. Due to the deformation dependency of the loading, a numerical scheme for the two-point boundary value problem in the semi-infinite space is developed. Failure pressures are computed based on a Rankine locus of failure. Variations of the order of 20% in the failure loads, as compared to the uniform beam model, are observed.

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.