BOUNDARY ELEMENT ANALYSIS OF CONTACT PROBLEMS USING ARTIFICIAL BOUNDARY NODE APPROACH

An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.

Two kinds of contact problems in three-dimensional icosahedral quasicrystals

Two kinds of contact problems, i.e., the frictional contact problem and the adhesive contact problem, in three-dimensional （3D） icosahedral quasicrystals are dis- cussed by a complex variable function method. For the frictional contact problem, the contact stress exhibits power singularities at the edge of the contact zone. For the adhe- sive contact problem, the contact stress exhibits oscillatory singularities at the edge of the contact zone. The numerical examples show that for the two kinds of contact problems, the contact stress exhibits singularities, and reaches the maximum value at the edge of the contact zone. The phonon-phason coupling constant has a significant effect on the contact stress intensity, while has little impact on the contact stress distribution regu- lation. The results are consistent with those of the classical elastic materials when the phonon-phason coupling constant is 0. For the adhesive contact problem, the indentation force has positive correlation with the contact displacement, but the phonon-phason cou- pling constant impact is barely perceptible. The validity of the conclusions is verified.

EFFECTIVE NUMERICAL METHODS FOR ELASTO-PLASTIC CONTACT PROBLEMS WITH FRICTION

Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.

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.

A NEW STATIC MECHANICS MODEL TO SOLVE CONTACT PROBLEMS IN MECHANICAL SYSTEMS

The coupling of the local contact problems between the components and the deformation of the components in the mechanical system were discovered.A series of coordinate systems have been founded to describe the mechanical system with the contact problems.The method of isolating the boundary of contact body from others has been used to describe the constraint between the contacting points.A more generalized static mechanics model of the mechanical system with the contact problems has been founded through the principle of virtual work.As an application,the model was used to study the multi_teeth engagement problems in the inner meshed planet gear systems.The stress distribution of contact gears was got.A test has verified that the static contact model and the computational method are right.

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.

A MULTIGRID SEMISMOOTH NEWTON METHOD FOR SEMILINEAR CONTACT PROBLEMS

This paper develops and analyzes multigrid semismooth Newton methods for a class of inequality-constrained optimization problems in function space which are motivated by and include linear elastic contact problems of Signorini type. We show that after a suitable Moreau-Yosida type regularization of the problem superlinear local convergence is obtained for a class of semismooth Newton methods. In addition, estimates for the order of tile error introduced by the regularization are derived. The main part of the paper is devoted to the analysis of a multilevel preconditioner for the semismooth Newton system. We prove a rigorous bound for the contraction rate of the multigrid cycle which is robust with respect to sufficiently small regularization parameters and the number of grid levels. Moreover, it applies to adaptively refined grids. The paper concludes with numerical results.

Two Kinds of Contact Problems in Dodecagonal Quasicrystals of Point Group 12 mm

In this paper,two kinds of contact problems in 2-D dodecagonal quasicrystals were discussed using the complex variable function method：one is the finite frictional contact problem,the other is the adhesive contact problem.The analytic expressions of contact stresses in the phonon and phason fields were obtained for a flat rigid punch,which showed that：（1） for the finite frictional contact problem,the contact stress exhibited power-type singularities at the edge of the contact zone;（2） for the adhesive contact problem,the contact stress exhibited oscillatory singularities at the edge of the contact zone.The distribution regulation of contact stress under punch was illustrated;and the low friction property of quasicrystals was verified graphically.

Book Review “Fabrikant,V.I.Contact and Crack Problems in Linear Theory of Elasticity,Bentham Science Publishers,Sharjah,UAE(2010)”

As suggested by the title, this extensive book is concerned with crack and contact prob- lems in linear elasticity. However, in general, it is intended for a wide audience ranging from engineers to mathematical physicists. Indeed, numerous problems of both academic and tech- nological interest in electro-magnetics, acoustics, solid and fluid dynamics, etc. are actually related to each other and governed by the same mixed boundary value problems from a unified mathematical standpoint

CONTACT PROBLEMS AND DUAL VARIATIONAL INEQUALITY OF 2－D ELASTOPLASTIC BEAM THEORY

In Order to study the frictional contact problems of the elastoplastic beam theory,an extended two-dimensional beam model is established, and a second order nonlinear equilibrium problem with both internal and external complementarity conditions is proposed. The external complementarity condition provides the free boundary condition. while the internal complemententarity condition gives the interface of the elastic and plastic regions. We prove that this bicomplementarity problem is equivalent to a nonlinear variational inequality The dual variational inequality is also developed.It is shown that the dual variational inequality is much easier than the primalvariational problem. Application to limit analysis is illustrated.

Reduced projection augmented Lagrange bi-conjugate gradient method for contact and impact problems

Based on the numerical governing formulation and non-linear complementary conditions of contact and impact problems, a reduced projection augmented Lagrange bi- conjugate gradient method is proposed for contact and impact problems by translating non-linear complementary conditions into equivalent formulation of non-linear program- ming. For contact-impact problems, a larger time-step can be adopted arriving at numer- ical convergence compared with penalty method. By establishment of the impact-contact formulations which are equivalent with original non-linear complementary conditions, a reduced projection augmented Lagrange bi-conjugate gradient method is deduced to im- prove precision and efficiency of numerical solutions. A numerical example shows that the algorithm we suggested is valid and exact.

Analytical solution of axisymmetric contact problem about indentation of a circular indenter into a soft functionally graded elastic layer

The paper addresses a contact problem of the theory of elasticity,i.e.,the penetration of a circular indenter with a flat base into a soft functionally graded elastic layer.The elastic properties of a functionally graded layer arbitrarily vary with depth,and the foundation is assumed to be elastic,yet much harder than a layer.Approximated analytical solution is constructed,and it is shown that the solutions are asymptotically exact both for large and small values of characteristic dimensionless geometrical parameter of the problem.Numerical examples are analyzed for the cases of monotonic and nonmonotonic variations of elastic properties.Numerical results for the case of homogeneous layer are compared with the results for nondeformable foundation.

AXISYMMETRIC CONTACT PROBLEM OF CUBIC QUASICRYSTALLINE MATERIALS

The axisymmetric elasticity theory of cubic quasicrystal was developed in Ref. [1]. The axisymmetric elasticity problem of cubic quasicrystal is reduced to a single higher-order partial differential equation by introducing a displacement function, based on which, the exact analytic solutions for the elastic field of an axisymmetric contact problem of cubic quasicrystalline materials are obtained for universal contact stress or contact displacement. The result shows that if the contact stress has order - 1/2 singularity on the edge of the contact domain, die contact displacement is a constant in the contact domain. Conversely, if the contact displacement is a constant, the contact stress must have order - 1/2 singularity on the edge of die contact domain.

BOUNDARY ELEMENT METHOD FOR MOVING AND ROLLING CONTACT OF 2D ELASTIC BODIES WITH DEFECTS

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.

Numerical modeling of large deformation and nonlinear frictional contact of excavation boundary of deep soft rock tunnel

Roadways excavated in soft rocks at great depth are difficult to be maintained due to large deformation of surrounding rocks, which greatly influences the safety and efficiency of deep resources exploitation. During the excavation process of a deep soft rock tunnel, the rock wall may be compacted due to large deformation. In this paper, the technique to address this problem by a two-dimensional (2D) finite element software, large deformation engineering analyses software (LDEAS 1.0), is provided. By using the Lagrange multiplier method, the kinematic constraint of non-penetrating condition and static constraint of Coulomb friction are introduced to the governing equations in the form of incremental displacement. The numerical example demonstrates the efficiency of this technology. Deformations of a transportation tunnel in inclined soft rock strata at the depth of 1 000 m in Qishan coal mine and a tunnel excavated to three different depths are analyzed by two models, i.e. the additive decomposition model and polar decomposition model. It can be found that the deformation of the transportation tunnel is asymmetrical due to the inclination of rock strata. For extremely soft rock, large deformation can converge only for the additive decomposition model. The deformation of surrounding rocks increases with the increase in the tunnel depth for both models. At the same depth, the deformation calculated by the additive decomposition model is smaller than that by the polar decomposition model.

Fundamental formulation for frictional contact with graded materials

In the paper, we develop the fundamental solutions for a graded half-plane subjected to concentrated forces acting perpendicularly and parallel to the surface. In the solutions, Young’s modulus is assumed to vary in the form of E(y)=E0eαy and Poisson’s ratio is assumed to be constant. On the basis of the fundamental solutions, the singular integral equations are formulated for the unknown traction distributions with Green’s function method. From the fundamental integral equations, a series of integral equations for special cases may be deduced corresponding to practical contact situations. The validity of the fundamental solutions and the integral equations is demonstrated with the degenerate solutions and two typical numerical examples.

THE EXACT INTEGRAL EQUATION OF HERTZ’S CONTACT PROBLEM

This paper presents the exact integral equation of Hertz's contact problem, which is obtained by taking into account the horizontal displacement of points in the contacted surfaces due to pressure.

UNIQUENESS FOR THE SOLUTIONS OF ELASTIC THIN PLATES AND SHALLOW SHELLS AS WELL AS THEIR CONTACT PROBLEM WITH HALF SPACE

The uniqueness for the solutions mentioned in the subject is proved by using the uniqueness of the solution for the internal boundary problem of Laplace and bi-Laplace equations of the first kind as well as of the second.

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.