带损伤弹性反问题的数值分析

考虑一类由椭圆性方程和热传导方程共同来刻画的准静态弹性模型,通过给定观测值来反演边界的牵引力.首先构造一个凸目标泛函,并引入Tikhonov正则化方法,使之极小化得到一个稳定的近似解.再用有限元离散求解,导出误差估计.最后,用数值例子说明算法的可行性和有效性.     

Numerical analysis of a bilateral frictional contact problem for linearly elastic materials

We consider a mathematical model which describes the contactbetween a linearly elastic body and an obstacle, the so-calledfoundation. The process is quasistatic and the contact is bilateral,i.e. there is no loss of contact during the process. The frictionis modelled with Tresca's law. The variational formulation ofthe problem is a nonlinear evolutionary inequality for the displacementfield which has a unique solution under certain assumptionson the given data. We study spatially semi-discrete and fullydiscrete schemes for the problem with finite-difference discretizationin time and finite-element discretization in space. The numericalschemes have unique solutions. We show the convergence of thescheme under the basic solution regularity. Under appropriateregularity assumptions on the solution, we derive optimal ordererror estimates. Finally, we present numerical results in thestudy of two-dimensional test problems.     

Numerical analysis and simulations of a quasistatic frictional contact problem with damage in viscoelasticity

We consider a quasistatic problem which models the bilateral contact between a viscoelastic body and a foundation, taking into account the damage and the friction. The damage which results from tension or compression is then involved in the constitutive law and it is modelled using a nonlinear parabolic inclusion. The variational problem is formulated as a coupled system of evolutionary equations for which we state the existence of a unique solution. Then, we introduce a fully discrete scheme using the finite element method to approximate the spatial variable and the Euler scheme to discretize the time derivatives. Error estimates are derived and, under suitable regularity hypotheses, the convergence of the numerical scheme obtained. Finally, a numerical algorithm and results are presented for some two-dimensional examples.     

The three-dimensional contact problem with friction for an elastic wedge

Solutions of three-dimensional boundary-value problems of the theory of elasticity are given for a wedge, on one face of which a concentrated shearing force is applied, parallel to its edge, while the other face is stress-free or is in a state of rigid or sliding clamping. The solutions are obtained using the method of integral transformations and the technique of reducing the boundary-value problem of the theory of elasticity to a Hilbert problem, as generalized by Vekua (functional equations with a shift of the argument when there are integral terms). Using these and previously obtained equations, quasi-static contact problems of the motion of a punch with friction at an arbitrary angle to the edge of the wedge are considered. In a similar way the contact area can move to the edge of a tooth in Novikov toothed gears. The method of non-linear boundary integral equations is used to investigate contact problems with an unknown contact area.     

The periodic contact problem for an elastic half-space

A method of solving the periodic contact problem for a system of indentors of arbitrary shape and an elastic half-space is proposed. Different versions of the arrangement of the indentors, at one and at several levels, are considered. The results are used to analyse the effect of the parameters of the microgeometry of the characteristics of a discrete contact and the stressed state of solids possessing regular microrelief.     

A frictionless contact problem for elastic–viscoplastic materials with normal compliance: Numerical analysis and computational experiments

Summary. In this paper we consider a frictionless contact problem between an elastic–viscoplastic body and an obstacle. The process is assumed to be quasistatic and the contact is modeled with normal compliance. We present a variational formulation of the problem and prove the existence and uniqueness of the weak solution, using strongly monotone operators arguments and Banach's fixed point theorem. We also study the numerical approach to the problem using spatially semi-discrete and fully discrete finite elements schemes with implicit and explicit discretization in time. We show the existence of the unique solution for each of the schemes and derive error estimates on the approximate solutions. Finally, we present some numerical results involving examples in one, two and three dimensions. Received May 20, 2000 / Revised version received January 8, 2001 / Published online June 7, 2001     

Numerical analysis of dynamic thermoviscoelastic contact with damage of a rod

We describe and analyse a model for a problem of thermoviscoelasticdynamic contact which allows for the general evolution of thematerial damage. The effects on the mechanical properties ofthe material due to crack expansion are described by a damagefield, which measures the decrease in the load-bearing capacityof the material. The damage process is assumed to be reversibleand the microcracks which develop as a result of tension orcompression may grow or disappear. The geometric setting isthat of a 1D rod which may contact a deformable obstacle. Thecontact is modelled by the normal compliance condition and thestress–strain constitutive equation is of Kelvin–Voigttype. The model consists of a coupled system of energy–elasticityequations together with a non-linear parabolic inclusion forthe damage field. The existence of a local weak solution isestablished using penalization, a finite element algorithm forthe solution is constructed and analysed and the results ofnumerical simulations based on this algorithm are presented.The simulations illustrate how the size of the normal compliancecoefficients, the damage rate coefficients and the applied forceaffect the character of the evolution of the damage. In particular,cycles of bonding and debonding can occur.     

A contact problem for an elastic plate with a thin rigid inclusion

Under study is an equilibrium problem for a plate under the influence of external forces. The plate is assumed to have a thin rigid inclusion that reaches the boundary at the zero angle and partially contacts a rigid body. On the inclusion face, there is a delamination. We consider the complete Kirchhoff–Love model, where the unknown functions are the vertical and horizontal displacements of the middle surface points of the plate. We present differential and variational formulations of the problem and prove the existence and uniqueness of a solution.     

An electro-viscoelastic frictional contact problem with damage

We consider a quasistatic frictional contact problem between a piezoelectric body and a foundation. The contact is modeled with normal compliance and friction is modeled with a general version of Coulomb's law of dry friction; the process is quasistatic and the material's behavior is described by an electro-viscoelastic constitutive law with damage. We derive a variational formulation for the model which is in the form of a system involving the displacement field, the electric potential field, and the damage field. Then we provide the existence of a unique weak solution to the model. The proof is based on arguments of evolutionary variational inequalities and fixed point.     

Numerical analysis of a non-clamped dynamic thermoviscoelastic contact problem

In this work, we analyze a non-clamped dynamic viscoelastic contact problem involving thermal effect. The friction law is described by a non-monotone relation between the tangential stress and the tangential velocity. This leads to a system of second-order inclusion for displacement and a parabolic equation for temperature. We provide a fully discrete approximation of the problem and find optimal error estimates without any smallness assumption on the data. The theoretical result is illustrated by numerical simulations.     

Numerical analysis of a contact problem including bone remodeling

In this work, a contact problem between an elastic body and a deformable obstacle is numerically studied. The bone remodeling of the material is also taken into account in the model and the contact is modeled using the normal compliance contact condition. The variational problem is written as a nonlinear variational equation for the displacement field, coupled with a first-order ordinary differential equation to describe the physiological process of bone remodeling. An existence and uniqueness result of weak solutions is stated. Then, fully discrete approximations are introduced based on the finite element method to approximate the spatial variable and an Euler scheme to discretize the time derivatives. Error estimates are obtained, from which the linear convergence of the algorithm is derived under suitable regularity conditions. Finally, some 2D numerical results are presented to demonstrate the behavior of the solution.     

A contact problem for a viscoelastic plate and an elastic beam

Under consideration is the problem of contact of a viscoelastic plate with an elastic beam. To characterize the viscoelastic deformation of the plate, the hereditary integrals are used. The differential formulation of the problem with the conditions in the form of a system of equalities and inequalities in the domain of possible contact is presented, and its equivalence to a variational inequality is proved. The unique solvability of the problem is proved as well as the existence of the time derivative of the solution. A limit problem is also considered as the bending rigidity of the plate tends to infinity.     

The axisymmetric contact problem for an elastic layer loaded with a deformable cover plate

The axisymmetric problem of the contact interaction of an elastic cover plate with an elastic layer, loaded at infinity with a uniform stretching force, directed parallel to the boundaries of the layer, is considered. The cover plate resists stretching but does not resist bending. The contact shearing stress under the cover plate, the displacement of the points of the cover plate and the deformation distortion coefficient of the elastic layer are determined.     

Numerical analysis of a contact problem between two thermoviscoelastic beams

We propose and analyze in this article a finite element approximation, based on a penalty formulation, to a quasi‐static unilateral contact problem between two thermoviscoelastic beams. An error bound is given and some numerical experiments discussed. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 644–661, 2011     

Numerical analysis of a frictional contact problem for viscoelastic materials with long-term memory

We consider a mathematical model which describes the bilateral contact between a deformable body and an obstacle. The process is quasistatic, the material is assumed to be viscoelastic with long memory and the friction is modeled with Tresca’s law. The problem has a unique weak solution. Here we study spatially semi-discrete and fully discrete schemes using finite differences and finite elements. We show the convergence of the schemes under the basic solution regularity and we derive order error estimates. Finally, we present an algorithm for the numerical realization and simulations for a two-dimensional test problem.     

An approximate solution of the axisymmetric contact problem for an elastic sphere

An axisymmetric, fractionally non-linear contact problem for an elastic sphere with a priori unknown boundary of the contactarea is considered. An integral equation for determining the density of the contact pressures is constructed taking account of the shear displacements of the boundary points of the elastic body. An approximate solution, which refines the equations of Hertz' theory, is constructed in the case of a small contact area.     

Shape sensitivity analysis of an elastic contact problem: Convergence of the Nitsche based finite element approximation

In a recent work, we introduced a finite element approximation for the shape optimization of an elastic structure in sliding contact with a rigid foundation where the contact condition (Signorini’s condition) is approximated by Nitsche’s method and the shape gradient is obtained via the adjoint state method. The motivation of this work is to propose an a priori convergence analysis of the numerical approximation of the variables of the shape gradient (displacement and adjoint state) and to show some numerical results in agreement with the theoretical ones. The main difficulty comes from the non-differentiability of the contact condition in the classical sense which requires the notion of conical differentiability.     

Numerical solution of an equilibrium problem for an elastic body with a thin delaminated rigid inclusion

Under consideration is a 2D-problem of elasticity theory for a body with a thin rigid inclusion. It is assumed that there is a delamination crack between the rigid inclusion and the elastic matrix. At the crack faces, the boundary conditions are set in the form of inequalities providing mutual nonpenetration of the crack faces. Some numerical method is proposed for solving the problem, based on domain decomposition and the Uzawa algorithm for solving variational inequalities.We give an example of numerical calculation by the finite element method.