Contact Problem in Decagonal Two-Dimensional Quasicrystal

As a new structure of solid matter quasicrystal brings profound new ideas to the traditional condensed matter physics, its elastic equations are more complicated than that of traditional crystal. A contact problem of decagonal two? dimensional quasicrystal material under the action of a rigid flat die is solved satisfactorily by introducing displacement function and using Fourier analysis and dual integral equations theory, and the analytical expressions of stress and displacement fields of the contact problem are achieved. The results show that if the contact displacement is a constant in the contact zone, the vertical contact stress has order -1/2 singularity on the edge of contact zone, which provides the important mechanics parameter for contact deformation of the quasicrystal.

On existence and uniqueness of the solution of elastoplastic contact problems

Contact problems and elastoplastic problems are unified and described by the variational inequality formulation, in which the constraints of the constitutional relations for elastoplastic materials and the contact conditions are relaxed totally. First, the coerciveness of the functional is proved. Then the uniqueness of the solution of variational inequality for the elastoplastic contact problems is demonstrated. The existence of the solution is also demonstrated according to the sufficient conditions for the solution of the elliptic variational inequality. A mathematical foundation is developed for the variational extremum principle of elastoplastic contact problems. The developed variational extremum forms can give an effective and strict mathematical modeling to solve contact problems with mathematical programming.

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.

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.

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.

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.

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.

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.

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.

A Perturbation Result for Dynamical Contact Problems

This paper is intended to be a first step towards the continuous dependence of dynamical contact problems on the initial data as well as the uniqueness of a solution. Moreover,it provides the basis for a proof of the convergence of popular time integration schemes as the Newmark method.We study a frictionless dynamical contact problem between both linearly elastic and viscoelastic bodies which is formulated via the Signorini contact conditions.For viscoelastic materials fulfilling the Kelvin-Voigt constitutive law,we find a characterization of the class of problems which satisfy a perturbation result in a non-trivial mix of norms in function space.This characterization is given in the form of a stability condition on the contact stresses at the contact boundaries.Furthermore,we present perturbation results for two well-established approximations of the classical Signorini condition:The Signorini condition formulated in velocities and the model of normal compliance,both satisfying even a sharper version of our stability condition.

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

Nano-Contact Problem with Surface Effects on Triangle Distribution Loading

This work presents a theoretical study of contact problem. The Fourier integral transform method based on the surface elasticity theory is adopted to derive the fundamental solution for the contact problem with surface effects, in which both the surface tension and the surface elasticity are considered. As a special case, the deformation induced by a triangle distribution force is discussed in detail. The results are compared with those of the classical contact problem. At nano-scale, the contributions of the surface tension and the surface elasticity to the stress and displacement are not equal at the contact surface. The surface tension plays a major role to the normal stress, whereas the shear stress is mainly affected by the surface elasticity. In addition, the hardness of material depends strongly on the surface effects. This study is helpful to characterize and measure the mechanical properties of soft materials through nanoindentation.

REDUCED AUGMENTED LAGRANGIAN BI-CONJUGATE GRADIENT METHOD FOR IMPACT-CONTACT PROBLEMS

To avoid the numerical oscillation of the penalty method andnon-compatibility with ex- plicit operators of conventional Lagrangemultiplier methods used in transient contact problems to en- forcesurface contact conditions, a new approach to enforcing surfacecontact constraints for the tran- sient nonlinear finite elementproblems, referred to as 'the reduced augmented Lagrangianbi-conjugate gradient method (ALCG)', is developed in this paper.Based on the nonlinear constrained optimization theory and iscompatible with the explicit time integration scheme, this approachcan also be used in implicit scheme naturally.

Study of adaptive finite element techniques for contact problem in elastic bodies

The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, then the Gauss - Newton iterations are used which allow the nonlinear problem to be transformed into a sequence of linear sub- problems then easily solved. In addition, the algorithm can be applied into the simulation of de -bonding of fiber - reinforced composites.

Mathematical Description and Finite Element Equation of 3D Coupled Thermo-elastic Contact Problem

Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-contact coupled problem have been listed. On this foundation, the finite element equation and definite solution condition of contact heat transfer have been given out. Based on virtual work principle and contact element technology, the finite element equation of 3D elastic contact system has been deduced under the effect of thermal stress. The pseudo load brought by contact gap have been introduced into this equation in order to reflect the contact state change. During iteration, once contact rigidity matrix is formed, it won’t change, which will make calculation reduce greatly.

Approximation of thermoelasticity contact problem with nonmonotone friction

The paper presents the formulation and approximation of a static thermoelasticity problem that describes bilateral frictional contact between a deformable body and a rigid foundation. The friction is in the form of a nonmonotone and multivalued law. The coupling effect of the problem is neglected. Therefore, the thermic part of the problem is considered independently on the elasticity problem. For the displacement vector, we formulate one substationary problem for a non-convex, locally Lipschitz continuous functional representing the total potential energy of the body. All problems formulated in the paper are approximated with the finite element method.

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.

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.