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.     

Finite element approximation to a contact problem for a nonlinear thermoviscoelastic beam

In this paper we propose and analyze a finite element method to the solution of a quasi-static contact problem between a nonlinear beam and a rigid obstacle. Error estimates and energy decay are obtained and some numerical simulations described.     

Damageable contact between an elastic body and a rigid foundation

In this work, the contact problem between an elastic body and a rigid obstacle is studied, including the development of material damage which results from internal compression or tension. The variational problem is formulated as a first-kind variational inequality for the displacements coupled with a parabolic partial differential equation for the damage field. The existence of a unique local weak solution is stated. Then, a fully discrete scheme is introduced using the finite element method to approximate the spatial variable and an Euler scheme to discretize the time derivatives. Error estimates are derived on the approximate solutions, from which the linear convergence of the algorithm is deduced under suitable regularity conditions. Finally, three two-dimensional numerical simulations are performed to demonstrate the accuracy and the behaviour of the scheme.     

Numerical analysis of two frictionless elastic-piezoelectric contact problems

In this work, we consider two frictionless contact problems between an elastic-piezoelectric body and an obstacle. The linear elastic-piezoelectric constitutive law is employed to model the piezoelectric material and either the Signorini condition (if the obstacle is rigid) or the normal compliance condition (if the obstacle is deformable) are used to model the contact. The variational formulations are derived in a form of a coupled system for the displacement and electric potential fields. An existence and uniqueness result is recalled. Then, a discrete scheme is introduced based on the finite element method to approximate the spatial variable. Error estimates are derived on the approximate solutions and, as a consequence, the linear convergence of the algorithm is deduced under suitable regularity conditions. Finally, some two-dimensional examples are presented to demonstrate the performance of the algorithm.     

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.     

Dynamic frictionless contact of a nonlinear beam with two stops

In this work, we consider mathematical and numerical approaches to a dynamic contact problem with a highly nonlinear beam, the so-called Gao beam. Its left end is rigidly attached to a supporting device, whereas the other end is constrained to move between two perfectly rigid stops. Thus, the Signorini contact conditions are imposed to its right end and are interpreted as a pair of complementarity conditions. We formulate a time discretization based on a truncated variational formulation. We prove the convergence of numerical trajectories and also derive a new form of energy balance. A fully discrete numerical scheme is implemented to present numerical results.     

A quasi-static Signorini contact problem for a thermoviscoelastic beam

This paper is concerned with the existence, uniqueness and numerical solution of a system of equations modelling the evolution of a quasi-static thermoviscoelastic beam that may be in contact with two rigid obstacles. A finite element approximation is proposed and analysed and some numerical results are given. Work partially supported by the Brazilian institution CNPq.     

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     

Existence for a thermoviscoelastic beam model of brakes

The existence of a weak solution to a model for the dynamic thermomechanical behavior of a viscoelastic beam, which is in frictional contact with a rigid rotating wheel, is established. The model describes a simple braking system in which a rotating wheel comes to a stop as a result of the frictional traction generated by the beam. The classical model consists of a system of coupled equations for the beam temperature and displacement, the wear of the beam's contacting end, the wheel temperature and its angular velocity. The weak formulation is an abstract differential inclusion involving set-valued pseudomonotone operators, The existence is proved by using recent results for such operators. Uniqueness is shown to hold when the wheel's angular velocity and temperature are known.     

Dynamic thermoviscoelastic thermistor problem with contact and nonmonotone friction

Abstract

The paper studies the evolution of the thermomechanical and electric state of a thermoviscoelastic thermistor that is in frictional contact with a reactive foundation. The mechanical process is dynamic, while the electric process is quasistatic. Friction is modeled with a nonmonotone relation between the tangential traction and tangential velocity. Frictional heat generation is taken into account and so is the strong dependence of the electric conductivity on the temperature. The mathematical model for the process is in the form of a system that consists of dynamic hyperbolic subdifferential inclusion for the mechanical state coupled with a nonlinear parabolic equation for the temperature and an elliptic equation for the electric potential. The paper establishes the existence of a weak solution to the problem by using time delays, a priori estimates and a convergence method.     

Dynamic vibrations of a damageable viscoelastic beam in contact with two stops

A model for the material damage, due to dynamic vibrations of a Kelvin‐Voigt viscoelastic beam whose tip is constrained to move between two stops, is presented and numerically analyzed. The contact of the free tip with the stops is described by the normal compliance condition. The evolution of damage of the beam's material, which measures the reduction of its load carrying capacity, is modeled with a parabolic inclusion. The existence of the unique local solution is stated. A numerical algorithm is presented, in which spatially it is approximated by finite elements, and the time derivatives are discretized with the Euler scheme. Error estimates are derived for sufficiently regular solutions, and four numerical simulations are shown. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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 porous thermoviscoelastic mixture problem: numerical analysis and computational experiments

In this paper, we study, from the numerical point of view, a porous thermoviscoelastic mixture problem. The mechanical problem is written as a linear coupled system of two hyperbolic partial differential equations for the porosities and a parabolic partial differential equation for the temperature field. An existence and uniqueness result and an energy decay property are stated. Then, fully discrete approximations are introduced by using the finite element method to approximate the spatial variable and the backward Euler scheme to discretize the time derivatives. A priori error estimates are proved from which, under suitable regularity conditions, the linear convergence of the algorithm is derived. Finally, some numerical simulations are presented to demonstrate the accuracy of the approximations in an academical one-dimensional example and the behaviour of the solutions in one- and two-dimensional problems.     

A dynamic contact problem in thermoviscoelasticity with two temperatures

This work is concerned with the study of a one-dimensional dynamic contact problem arising in thermoviscoelasticity with two temperatures. The existence and uniqueness of a solution to the continuous problem is established using the Faedo–Galerkin method. A finite element approximation is proposed, a convergence result given and some numerical simulations described.     

A convergence result in the study of bone remodeling contact problems

We consider the approximation of a bone remodeling model with the Signorini contact conditions by a contact problem with normal compliant obstacle, when the obstacle's deformability coefficient converges to zero (that is, the obstacle's stiffness tends to infinity). The variational problem is a coupled system composed of a nonlinear variational equation (in the case of normal compliance contact conditions) or a variational inequality (for the case of Signorini's contact conditions), for the mechanical displacement field, and a first-order ordinary differential equation for the bone remodeling function. A theoretical result, which states the convergence of the contact problem with normal compliance contact law to the Signorini problem, is then proved. Finally, some numerical simulations, involving examples in one and two dimensions, are reported to show this convergence behaviour.     

Global existence and exponential stability for a contact problem between two thermoelastic beams

In this paper we analyze a dynamic unilateral contact problem between two thermoelastic beams. We establish the existence of a weak global-in-time solution, by a penalization method. Moreover, we study the asymptotic behavior of such a solution proving that the energy associated to the system decays exponentially to zero, as time goes to infinity.     

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.     

A coupling of FEM-BEM for a kind of Signorini contact problem

In this paper, we consider a kind of coupled nonlinear problem with Signorini contact conditions. To solve this problem, we discuss a new coupling of finite element and boundary element by adding an auxiliary circle. We first derive an asymptotic error estimate of the approximation to the coupled FEM-BEM variational inequality. Then we design an iterative method for solving the coupled system, in which only three standard subproblems without involving any boundary integral equation are solved. It will be shown that the convergence speed of this iteration method is independent of the mesh size.     

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

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

Existence of a solution for a Signorini contact problem for Maxwell-Norton materials

The aim of this article is to study the quasistatic evolutionof a Maxwell–Norton three-dimensional viscoelastic solidwith contact constraints. After introducing the appropiate functionalframework, we will discretize the problem in time using an implicitscheme whose resultant variational inequality is well posed.By using monotonicity arguments together with compensated compactnesstechniques, we will prove that the corresponding discrete solutionconverges to a solution of the continuous problem.