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.     

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.     

Numerical analysis of a quasi-static contact problem for a thermoviscoelastic beam

In this paper we revisit a quasi-static contact problem of a thermoviscoelastic beam between two rigid obstacles which was recently studied in [1]. The variational problem leads to a coupled system, composed of an elliptic variational inequality for the vertical displacement and a linear variational equation for the temperature field. Then, its numerical resolution is considered, based on the finite element method to approximate the spatial variable and the implicit Euler scheme to discretize the time derivatives. Error estimates are proved from which, under adequate regularity conditions, the linear convergence is derived. Finally, some numerical simulations are presented to show the accuracy of the algorithm and the behavior of the solution.     

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.     

Solvability of a dynamic contact problem between a piezoelectric body and a conductive foundation

We consider a mathematical model which describes the dynamic process of contact between a piezoelectric body and an electrically conductive foundation. We model the material’s behavior with a nonlinear electro-viscoelastic constitutive law; the contact is frictionless and is described with the normal compliance condition and a regularized electrical conductivity condition. We derive a variational formulation for the problem and then, under a smallness assumption on the data, we prove the existence of a unique weak solution to the model. We also investigate the behavior of the solution with respect the electric data on the contact surface and prove a continuous dependence result. Then, we introduce a fully discrete scheme, based on the finite element method to approximate the spatial variable and the backward Euler scheme to discretize the time derivatives. We treat the contact by using a penalized approach and a version of Newton’s method. We implement this scheme in a numerical code and, in order to verify its accuracy, we present numerical simulations in the study of two-dimensional test problems. These simulations provide a numerical validation of our continuous dependence result and illustrate the effects of the conductivity of the foundation, as well.     

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

考虑一类由椭圆性方程和热传导方程共同来刻画的准静态弹性模型,通过给定观测值来反演边界的牵引力.首先构造一个凸目标泛函,并引入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.     

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.     

Fully discrete finite element approximations of the forced Fisher equation

We consider fully discrete finite element approximations of the forced Fisher equation that models the dynamics of gene selection/migration for a diploid population with two available alleles in a multidimensional habitat and in the presence of an artificially introduced genotype. Finite element methods are used to effect spatial discretization and a nonstandard backward Euler method is used for the time discretization. Error estimates for the fully discrete approximations are derived by applying the Brezzi-Rappaz-Raviart theory for the approximation of a class of nonlinear problems. The approximation schemes and error estimates are applicable under weaker regularity hypotheses than those that are typically assumed in the literature. The algorithms and analyses, although presented in the concrete setting of the forced Fisher equation, also apply to a wide class of semilinear parabolic partial differential equations.     

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.     

Mathematical and numerical approaches to a one-dimensional dynamic thermoviscoelastic contact problem

In this paper, we propose numerical schemes for solving a nonlinear system which consists of a coupled partial differential equations and two conditions, called normal compliance contact condition and Barber’s heat exchange condition. The convergence of numerical trajectories is shown by using a time discretization and passing the limit of the time step size. The uniqueness of the weak solution is proved as well. We derive the extensive form of an energy balance which will be a criterion to examine numerical stability. An example is provided to present and discuss numerical results.     

Numerical Analysis of a Problem Involving a Viscoelastic Body with Double Porosity

We study from a numerical point of view a multidimensional problem involving a viscoelastic body with two porous structures. The mechanical problem leads to a linear system of three coupled hyperbolic partial differential equations. Its corresponding variational formulation gives rise to three coupled parabolic linear equations. An existence and uniqueness result, and an energy decay property, are recalled. Then, fully discrete approximations are introduced using the finite element method and the implicit Euler scheme. A discrete stability property and a priori error estimates are proved, from which the linear convergence of the algorithm is derived under suitable additional regularity conditions. Finally, some numerical simulations are performed in one and two dimensions to show the accuracy of the approximation and the behaviour of the solution.     

Deterministic particle method approximation of a contact inhibition cross-diffusion problem

We use a deterministic particle method to produce numerical approximations to the solutions of an evolution cross-diffusion problem for two populations.According to the values of the diffusion parameters related to the intra- and inter-population repulsion intensities, the system may be classified in terms of an associated matrix. When the matrix is definite positive, the problem is well posed and the finite element approximation produces convergent approximations to the exact solution.A particularly important case arises when the matrix is only positive semi-definite and the initial data are segregated: the contact inhibition problem. In this case, the solutions may be discontinuous and hence the (conforming) finite element approximation may exhibit instabilities in the neighborhood of the discontinuity.In this article we deduce the particle method approximation to the general cross-diffusion problem and apply it to the contact inhibition problem. We then provide some numerical experiments comparing the results produced by the finite element and the particle method discretizations.     

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 a dual-phase-lag model involving two temperatures

In this paper, we numerically analyse a phase-lag model with two temperatures which arises in the heat conduction theory. The model is written as a linear partial differential equation of third order in time. The variational formulation, written in terms of the thermal acceleration, leads to a linear variational equation, for which we recall an existence and uniqueness result and an energy decay property. Then, using the finite element method to approximate the spatial variable and the implicit Euler scheme to discretize the time derivatives, fully discrete approximations are introduced. A discrete stability property is proved, and a priori error estimates are obtained, from which the linear convergence of the approximation is derived. Finally, some one-dimensional numerical simulations are described to demonstrate the accuracy of the approximation and the behaviour of the solution.     

Analysis and numerical solution of a piezoelectric frictional contact problem

We consider a mathematical model which describes the frictional contact between an electro-elastic–visco-plastic body and a conductive foundation. The contact is modelled with normal compliance and a version of Coulomb’s law of dry friction, in which the stiffness and the friction coefficients depend on the electric potential. We derive a variational formulation of the problem and we prove an existence and uniqueness result. The proof is based on a recent existence and uniqueness result on history-dependent quasivariational inequalities obtained in [15]. Then we introduce a fully discrete scheme for solving the problem and, under certain solution regularity assumptions, we derive an optimal order error estimate. Finally, we present some numerical results in the study of a two-dimensional test problem which describes the process of contact in a microelectromechanical switch.     

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     

阻尼Sine-Gordon方程的H1-Galerkin混合元方法数值解

利用H1-Galerkin混合有限元方法讨论阻尼Sine-Gordon方程,得到一维情况下半离散和全离散格式的最优阶误差估计,并且推广应用到二维和三维情况,而且不用验证LBB相容性条件.