We consider two mathematical models which describe the frictionless process of contact between a rate-type viscoplastic body and a foundation. The contact is modelled with normal compliance and memory term such that penetration is not restricted in the first problem, but is restricted with unilateral constraint in the second one.
For each problem, we derive a variational formulation in terms of displacements, which is in a form of a history-dependent variational equation and a history-dependent variational inequality. Then we prove the unique weak solvability of each model. Next, we prove the convergence of the weak solution of the first problem and the weak solution of the second problem, as the stiffness coefficient of the foundation converges to infinity.
Finally, we provide numerical simulations which illustrate this convergence result.
Authors
Mikael Barboteu
(Laboratoire de Mathématiques et Physique, Université de Perpignan)
Flavius Patrulescu
(Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy)
Ahmad Ramadan
(Laboratoire de Mathématiques et Physique, Université de Perpignan)
Mircea Sofonea
(Laboratoire de Mathématiques et Physique, Université de Perpignan)
M. Barboteu, F. Pătrulescu, A. Ramadan, M. Sofonea, History-dependent contact models for viscoplastic materials, IMA J. Appl. Math., 79 (2014) no. 6, pp. 1180-1200.
[1] P. Alart and A. Curnier, A mixed formulation for frictional contact problems prone to Newton like solution methods, Computer Methods in Applied Mechanics 26 and Engineering 92 (1991), 353–375.
[2] M. Barboteu, A. Matei, M. Sofonea, Analysis of quasistatic viscoplastic contact problems with normal compliance, Quarterly Journal of Mechanics and Applied Mathematics, 65 (2012), 555–579.
[3] N. Cristescu and I. Suliciu, Viscoplasticity, Martinus Nijhoff Publishers, Editura Tehnica, Bucharest, 1982.
[4] W. Han and M. Sofonea, Quasistatic Contact Problems in Viscoelasticity and Viscoplasticity, Studies in Advanced Mathematics 30, American Mathematical Society–International Press, 2002.
[5] I.R. Ionescu and M. Sofonea, Functional and Numerical Methods in Viscoplasticity, Oxford University Press, Oxford, 1993.
[6] H. B. Khenous, P. Laborde and Y. Renard, On the discretization of contact problems in elastodynamics, Lecture Notes in Applied Computational Mechanics 27 (2006), 31–38.
[7] H. B. Khenous, J. Pommier, Y. Renard, Hybrid discretization of the Signorini problem with Coulomb friction. Theoretical aspects and comparison of some numerical solvers, Applied Numerical Mathematics 56 (2006), 163–192.
[8] A. Klarbring, A. Mikeliˇc and M. Shillor, Frictional contact problems with normal compliance, Int. J. Engng. Sci. 26 (1988), 811–832.
[9] A. Klarbring, A. Mikeliˇc and M. Shillor, On friction problems with normal compliance, Nonlinear Analysis 13 (1989), 935–955.
[10] T. Laursen, Computational Contact and Impact Mechanics, Springer, Berlin, 2002.
[11] J. A. C. Martins and J. T. Oden, Existence and uniqueness results for dynamic contact problems with nonlinear normal andfriction interface laws, Nonlinear Analysis : Theory, Methods and Applications 11 (1987), 407–428.
[12] J.T. Oden and J.A.C. Martins, Models and computational methods for dynamic friction phenomena, Computer Methods in Applied Mechanics and Engineering 52 (1985), 527–634.
[13] M. Shillor, M. Sofonea and J.J. Telega, Models and Analysis of Quasistatic Contact, Lecture Notes in Physics, 655, Springer, Berlin, 2004.
[14] M. Sofonea, C. Avramescu and A. Matei, A Fixed point result with applications in the study of viscoplastic frictionless contact problems, Communications on Pure and Applied Analysis 7 (2008), 645–658.
[15] M. Sofonea and A. Matei, History-dependent quasivariational inequalities arising in Contact Mechanics, European Journal of Applied Mathematics, 22 (2011), 471–491.
[16] M. Sofonea and F. P˘atrulescu, Analysis of a history-dependent frictionless contact problem, Mathematics and Mechanics of Solids, 18, no. 4 (2013), 409–430.
[17] P. Wriggers, Computational Contact Mechanics, Wiley, Chichester, 2002.
Paper (preprint) in HTML form
History-dependent Contact Models for Viscoplastic Materials
M. Barboteu 1, F. Pătrulescu 2, A. Ramadan 1 and M. Sofonea 1 1 Laboratoire de Mathématiques et Physique
Université de Perpignan, 52 Avenue de Paul Alduy, 66860 Perpignan, France
2 Tiberiu Popoviciu Institute of Numerical Analysis
P.O. Box 68-1, 400110 Cluj-Napoca, Romania
Abstract
We consider two quasistatic contact problems which describe the contact between a viscoplastic body and an obstacle, the so-called foundation. The contact is frictionless and is modelled with normal compliance and memory term of such a type that the penetration is not restricted in the first problem, but is restricted with unilateral constraint, in the second one. For each problem we derive a variational formulation, then we prove its unique solvability. Next, we prove the convergence of the weak solution of the first problem to the weak solution of the second problem, as the stiffness coefficient of the foundation converges to infinity. And, finally, we provide numerical simulations which illustrate the convergence result.
1 Introduction
Phenomena of contact between deformable bodies or between deformable and rigid bodies abound in industry and everyday life. Contact of braking pads with wheels, tires with roads, pistons with skirts are just a few simple examples. Common industrial processes such as metal forming, metal extrusion, involve contact evolutions and, for this reason, a considerable effort has been developed in their modelling, mathematical analysis and numerical solution. Owing to their inherent complexity, contact phenomena lead to nonlinear and nonsmooth mathematical problems.
The aim of this paper is to study two frictionless contact problems for rate-type viscoplastic materials, within the framework of the Mathematical Theory of Contact Mechanics. We model the behavior of the material with a constitutive law of the form
(1.1)
where denotes the displacement field, represents the stress and is the linearized strain tensor. Here is a fourth order tensor which describes the elastic properties of the material and is a constitutive function which describes its viscoplastic behavior. In (1.1) and everywhere in this paper the dot above a variable represents derivative with respect to the time variable .
Various results, examples and mechanical interpretations in the study of viscoplastic materials of the form (1.1) can be found in [3, 5] and references therein. Displacement-traction boundary value problems with such materials were considered in [5], both in the dynamic and quasistatic case. Quasistatic frictionless and frictional contact problems for materials of the form (1.1) were studied in various papers, see 4 and 13 for a survey. There, various models of contact were stated and they variational analysis, including existence and uniqueness results, was provided. The numerical analysis of the corresponding models can be found in [4] and the references therein. In all the above papers the process of contact was studied in a finite interval of time.
In [14] a quasistatic frictionless contact problem for viscoplastic materials of the form (1.1) was considered. The process was assumed to be quasistatic and the contact was modelled by using the normal compliance condition with infinite penetration. The unique solvability of the solution was obtained by using a fixed point argument. In contrast, in [2] was considered a problem with normal compliance and finite penetration and was proved the unique solvability of the models using new arguments on history-dependent variational inequalities presented in [15]. The present paper represents a continuation of [2] and we consider the problem with normal compliance, finite penetration and memory term. The same contact condition for viscoelastic materials was used in [16]. Also, we state and prove the convergence of the solution of the problem with infinite penetration to the solution of the problem with finite penetration as the stiffness coefficient converges to infinity. And, finally, we provide numerical simulations which illustrate this convergence.
The rest of the paper is structured as follows. In Section 2 we introduce the notations we shall use as well as some preliminary material. In Section 3 we present the classical formulation of the two contact problems. In Section 4 we list the assumptions on the data and derive the variational formulation of the problems. Then we state and prove the unique weak solvability of each model. In Section 5 we state and prove a converge result, Theorem 5.1. Next, in Section 6 we present the numerical solution of the contact problem with normal compliance restricted by unilateral constraints and memory term. And, finally, in Section 7 some numerical simulations are presented including a numerical validation of the convergence result.
2 Notations and preliminaries
Everywhere in this paper we use the notation for the set of positive integers and will represent the set of non negative real numbers, i.e. . We denote by the space of second order symmetric tensors on or, equivalently, the space of symmetric matrices of order . The inner product and norm on and are defined by
Let be a bounded domain with a Lipschitz continuous boundary and let be a measurable part of such that meas . We use the notation for a typical point in and we denote by the outward unit normal at . Here and below the indices run between 1 and and, unless stated otherwise, the summation convention over repeated indices is used. An index that follows a comma represents the partial derivative with respect to the corresponding component of the spatial variable, e.g. . We consider the spaces
These are real Hilbert spaces endowed with the inner products
and the associated norms and , respectively. Here represents the deformation operator given by
Completeness of the space ( ) follows from the assumption meas , which allows the use of Korn’s inequality.
For an element we still write for the trace of and we denote by and the normal and tangential components of on given by . Let be a measurable part of . Then, by the Sobolev trace theorem, there exists a positive constant which depends on and such that
(2.1)
Also, for a regular stress function we use the notation and for the normal and the tangential traces, i.e. and . Moreover, we recall that the divergence operator is defined by the equality and, finally, the following Green’s formula holds:
(2.2)
For a normed space ( ) we use the notation for the space of continuously functions defined on with values on , and for the space of continuous differentiable functions defined on with values on . For a subset we still use the symbols and for the set of continuous and continuously differentiable functions defined on with values on , respectively.
Let be a real Hilbert space with inner product and associated norm . Assume given a set , the operators and a function such that:
(2.3)
(2.4)
We proceed with the following existence and uniqueness result in the study of nonlinear equations involving monotone operators.
Theorem 2.1 Let be a Hilbert space and let be a strongly monotone Lipschitz continuous operator. Then, for each there exists a unique element such that .
The following result, proved in [15, will be used in Section 4 of this paper.
Theorem 2.2 Assume that (2.3)-(2.6) hold. Then there exists a unique function such that for all , the inequality below holds:
(2.7)
We have the following consequence of Theorem 2.2.
Corollary 2.3 Let be a Hilbert space and assume that (2.3)-(2.6) hold. Then there exists a unique function such that
(2.8)
Following the terminology introduced in [15] we refer to (2.7) as a history-dependent quasivariational inequality. To avoid any confusion, we note that here and below the notation and are short hand notation for and , i.e. and , for all .
3 The models
In this section we present the two problems which describe the frictionless contact process and present the assumption on the data. The physical setting is as follows. A viscoplastic body occupies a bounded domain with a Lipschitz continuous boundary , divided into three measurable parts and , such that meas . The body is subject to the action of body forces of density . We also assume that it is fixed on and surface tractions of density act on . On , the body is in frictionless contact with a deformable obstacle, the so-called foundation. We assume that the problem is quasistatic, and we study the contact process in the interval of time .
In the first problem, unlike in [2], the contact is modelled with normal compliance and memory term in such a way that the penetration is not limited. Under these conditions, the classical formulation of the problem is the following.
Problem . Find a displacement field and a stress field : such that
Here and below, in order to simplify the notation, we do not indicate explicitly the dependence of various functions on the variables or . Equation (3.1) represents the viscoplastic constitutive law of the material introduced in Section 1 and equation (3.2) is the equilibrium equation. Conditions (3.3) and (3.4) are the displacement and traction boundary conditions, respectively, and condition (3.6) shows that the tangential stress on the contact surface, denoted , vanishes. We use it here since we assume that the contact process is frictionless. Finally, (3.7) represents the initial
conditions in which and denote the initial displacement and the initial stress field, respectively.
The function is Lipschitz continuous, increasing and vanishes for a negative argument, i.e.
In the second problem the contact is again modelled with normal compliance and memory term but in such a way that the penetration is limited and associated to a unilateral constraint. The classical formulation of the problem is the following.
Problem . Find a displacement field and a stress field : such that
(3.9)
(3.10)
(3.11)
(3.12)
(3.13)
(3.14)
(3.15)
Here is given and is a function which satisfies (3.8). Conditions (3.9)(3.12) and (3.14) - (3.15) have the same interpretation as in the contact problem .
We now present the new contact condition (3.13), condition (3.5) can be presented using similar arguments. It can be derived in the following way. First, we assume that the penetration is limited by the bound and, therefore, at each time moment , the normal displacement satisfies the inequality
(3.16)
Next, we assume that the normal stress has an additive decomposition of the form
(3.17)
in which the functions and describe the deformability, the rigidity and the memory properties of the foundation, at each . Also, we assume that the function satisfies the normal compliance contact condition
(3.18)
Condition (3.18) combined with assumption (3.8) shows that when there is separation between the body and the obstacle (i.e. when ), then the reaction of the foundation vanishes (since ); also, when there is penetration (i.e. when 0 ), then the reaction of the foundation is towards the body (since ) and it is increasing with the penetration (since is an increasing function). Finally, we note that in this condition the penetration is not restricted and the normal stress is uniquely determined by the normal displacement.
Condition (3.18) was first introduced in [11, 12] in the study of dynamic contact problems with elastic and viscoelastic materials. The term normal compliance for this condition was first used in (8, 9. A first example of normal compliance function which satisfies condition (3.8) is
(3.19)
where and is a positive constant. In this case condition (3.18) shows that the reaction of the foundation is proportional to the penetration and, therefore, (3.8), (3.18) model the contact with a linearly elastic foundation. A second example of normal compliance function which satisfies condition (3.8) is given by
where is a positive coefficient related to the wear and hardness of the surface and, again, . In this case the contact condition (3.18) means that when the penetration is too large, i.e. when it exceeds , the obstacle backs off and offers no additional resistance to the penetration. We conclude that in this case the foundation has an elastic-perfectly plastic behavior.
The part of the normal stress satisfies the Signorini condition in the form with a gap function, i.e.
(3.20)
And, finally, the function satisfies the memory condition
(3.21)
in which represents a given function, the so-called surface memory function. Contact conditions of the form (3.21) have a simple physical interpretation if there are no cycles
of contact and separation during the time interval of interest. For instance, assume in what follows that is a positive function. Moreover, assume that in the time interval there is only penetration (i.e. for all ). Then (3.21) shows that the reaction of the foundation at is towards the body (since ). Also, if in the time interval there is separation (i.e. for all ) then there is no reaction at the moment (since ).
Now, assume a situation in which is positive in time interval and negative on the time interval . Then, following (3.21) we have
since the integral on the remaining interval vanishes. Assume, in addition, that the support of the function is included in the interval with . Two possibilities arise. First, if it follows that for all and (3.21) shows the normal stress vanishes. Second, if (3.21) implies that i.e. a residual pression exists at the moment on the body’s surface. We interpret this as a memory effect in which the foundation prevents the separation, moves towards the body and exerts a pression on a short interval of time of length . Various other mechanical interpretation of the condition (3.21) could be obtained if is assumed to be a negative function.
We combine equalities (3.17), (3.18) and (3.21) to see that
(3.22)
Then we substitute equality (3.22) in (3.20) and use inequality (3.16) to obtain the contact condition (3.13).
4 Existence and uniqueness results
In this section we list the assumptions on the data, derive the variational formulations of the problems and and then we state and prove their unique weak solvability. To this end we assume that the elasticity tensor and the constitutive function satisfy the following conditions.
The surface memory function satisfies
(4.3)
We also assume that the body forces and the surface tractions have the regularity
(4.4)
In the study of Problem we assume that the initial data satisfy
(4.5)
and, finally, in the study of Problem we assume that
(4.6)
where denotes the set of admissible displacements defined by
(4.7)
In the rest of the section we denote by a positive generic constant that may depend on time and whose value may change from line to line. Also, we use the symbol " " to denote the weak convergence in the Hilbert space .
We turn now to the variational formulation of the problems and . To this end, we use Riesz’s representation Theorem to define the operator and the function by equalities
(4.8)
(4.9)
Assume in what follows that ( ) are sufficiently regular functions which satisfy (3.1)-(3.7) and let be given. We integrate equation (3.1) with the initial conditions (3.7) to obtain
(4.10)
Then we use the Green formula (2.2), the equilibrium equation (3.2), the boundary conditions (3.3)-(3.6) and notation (4.8)-(4.9) to see that
(4.11)
We present the following existence and uniqueness result proved in 2 .
Lemma 4.1 Assume that (4.2) and (4.5) hold. Then, for each function there exists a unique function such that
(4.12)
Moreover, the operator satisfies the following property: for every there exists such that
(4.13)
Using the operator defined in Lemma 4.1 we deduce that (4.10) and (4.11) are equivalent with
We use again Riesz’s representation Theorem to define the operator by equality
(4.14)
and we obtain the following variational formulation of the Problem .
Problem . Find a displacement field and a stress field , such that
(4.15)
(4.16)
hold, for all .
In the study of the problem we have the following existence and uniqueness result.
Theorem 4.2 Assume that (3.8) and (4.1)-(4.5) hold. Then, Problem has a unique solution, which satisfies
(4.17)
Proof. We define the operator by equality
(4.18)
With this notation we consider the problem of finding a function such that, for all , the following equality holds
(4.19)
To solve (4.19) we employ Corollary 2.3 with . We use (4.1), (3.8) and (2.1) to see that the operator verifies condition (2.4), i.e. it is strongly monotone and Lipschitz continuous. In addition, using (4.3) we note that the operator satisfies condition (2.5) with
(4.20)
for every .
Finally, using (4.4) and (4.9) we deduce that has the regularity expressed in (2.6). It follows now from Corollary 2.3 that there exists a unique function which solves the equality (4.19), for any .
Based on the results above we deduce the existence of a unique function which satisfies (4.16) for any . Let be defined by (4.15). Then it follows that the couple ( ) is the unique couple of functions with regularity (4.17) which satisfies (4.15)-(4.16).
Assume now that ( ) are sufficiently regular functions which satisfy (3.9)-(3.15) and, again, let be given. Then, using similar arguments as above we obtain the following variational formulation of Problem .
Problem . Find a displacement field and a stress field , such that
(4.21)
(4.22)
hold, for all .
In the study of the problem we have the following existence and uniqueness result.
Theorem 4.3 Assume that (3.8), (4.1) -(4.4) and (4.6) hold. Then, Problem has a unique solution, which satisfies
(4.23)
Proof. We use the Theorem 2.2 with and and arguments similar to those used in the proof of Theorem 4.2.
5 A convergence result
Everywhere in this section we assume that the function satisfies condition (3.8) and let be a function which satisfies
Let and consider the function defined by
Using assumption (5.1) it follows that the function satisfies condition (3.8), i.e.
This allows us to consider the operator defined by
(5.4)
and, moreover, we note that is a monotone, Lipschitz continuous operator.
With these notation, we consider the following contact problem.
Problem . Find a displacement field and a stress field such that
(5.5)
(5.6)
(5.7)
(5.8)
(5.9)
(5.10)
(5.11)
The equations and boundary conditions in problem (5.5)-(5.11) have a similar interpretations as those in problem (3.9)-(3.15). The difference arises in the fact that here we replace the contact condition with normal compliance, unilateral constraint and memory term (3.13) with the contact condition with normal compliance and memory term (5.9). In this condition represents a penalization parameter which may be interpreted as a deformability of the foundation, and then is the surface stiffness coefficient. Indeed, when is smaller the reaction force of the foundation to penetration is larger and so the same force will result in a smaller penetration, which means that the foundation is less deformable. When is larger the reaction force of the foundation to penetration is smaller, and so the foundation is less stiff and more deformable.
Note that here and below is the normal component of the displacement field and represent the normal and tangential components of the stress tensor , respectively.
Assume now that (4.1)-(4.4) and (4.5) hold. Using arguments similar as in Section 4 for contact problem we obtain the following variational formulation for Problem .
Problem . Find a displacement field and a stress field , such that
(5.12)
(5.13)
hold, for all .
It follows from Theorem 4.2 that Problem has a unique solution ( ) which satisfies (4.17). Finally, it follows from Theorem 4.3 that Problem has a
unique solution ( ) which satisfies (4.23). The behavior of the solution ( ) as is given in the following result.
Theorem 5.1 Assume that (3.8), (4.1)-(4.4), (4.6) and (5.1) hold. Then, the solution ( ) of Problem converges to the solution ( ) of Problem , that is
(5.14)
as , for all .
In addition to the mathematical interest in the result above, this result is important from the mechanical point of view, since it shows that the weak solution of the viscoplastic contact problem with normal compliance, memory term and finite penetration may be approached as closely as we wish by the solution of the viscoplastic contact problem with normal compliance, memory term and infinite penetration, with a sufficiently small deformability coefficient.
The proof of Theorem 5.1 is carried out in several steps.
Let . In the first step we consider the auxiliary problem of finding a displacement field such that, for all ,
(5.15)
This problem is an intermediate problem between (5.13) and (4.22), since here is known, taken from the problem with finite penetration .
We have the following existence and uniqueness result.
Lemma 5.2 There exists a unique function which satisfies (5.15), for all .
Proof. We define the operator and the function by equalities
(5.16)
(5.17)
and note that (4.3), (4.4), (4.9), (4.12) and (4.14) yield
(5.18)
Let . Based on (5.16)-(5.17), it is easy to see that the variational equation (5.15) is equivalent with the nonlinear equation
(5.19)
Next, as in the proof of Theorem 4.2 by (4.1) and the properties of operator it follows that is a strongly monotone and Lipschitz continuous operator. And, Theorem 2.1 implies the existence of a unique solution for the nonlinear equation (5.19), which concludes the proof.
We proceed with the following convergence result.
Lemma 5.3 As ,
for all .
Proof. Let . We take in (5.15) to obtain
(5.20)
On the other hand, the properties (5.3) of the function yield
(5.21)
We combine (5.20), (5.21) and use (4.1) (c) to obtain that
(5.22)
This inequality shows that the sequence is bounded. Hence, there exists a subsequence of the sequence , still denoted , and an element such that
(5.23)
It follows from (5.20) that
and, since is a bounded sequence in , we deduce that
This implies that
and, since , it follows that
(5.24)
We consider now the measurable subsets of defined by
(5.25)
Clearly, both and depend on and but, for simplicity, we do not indicate explicitly this dependence. We use (5.24) to write
and, since
we obtain
Thus, taking into account that for , by the monotonicity of the function we can write
Therefore, we deduce that
(5.26)
We use now the definitions (5.2) and (5.25) to see that, a.e. on , we have
Consequently, the inequality (5.26) yields
(5.27)
Next, we consider the function defined by
and we note that by (5.1) it follows that is a continuous increasing function and, moreover,
(5.28)
We use (5.27), equality a.e. on and (5.25) to deduce that
where denotes the positive part of . Therefore, passing to the limit as , by using (5.23) as well as compactness of the trace operator we find that
Since the integrand is positive a.e. on , the last inequality yields
and, using (5.28) and definition (4.7) we conclude that
(5.29)
Next, we test in (5.15) with , where , to obtain
(5.30)
Since we have a.e. on . Taking into account this equality and the monotonicity of the function we have
and, therefore, by using (5.4) we obtain
(5.31)
Then, using (5.31) and (5.30) we find that
(5.32)
We pass to the lower limit in (5.32) and use (5.23) to obtain
(5.33)
Next, we take in (4.22) and in (5.33). Then, adding the resulting inequalities we find that
This inequality combined with (4.1) implies that
It follows from here that the whole sequence is weakly convergent to the element , which concludes the proof.
We proceed with the following strong convergence result.
Lemma 5.4 As ,
for all .
Proof. Let . Using (4.1) we write
Next, we take in (5.32) to obtain
and, therefore, combining the above inequalities we find that
We pass to the upper limit in this inequality and use Lemma 5.3 to conclude the proof of the lemma.
We are now in position to provide the proof of Theorem 5.1.
Proof. Let and let be such that . Let also . Then, testing with in (5.15) and (5.13), we have
We subtract the previous equalities and use the monotonicity of the operator to deduce that
and, therefore,
(5.34)
We use (5.34) to find that
where is given by (4.20). It follows from here that
and, using a Gronwall argument, we obtain
(5.35)
Note that for all and, therefore, (5.35) yields
(5.36)
On the other hand, by estimate (5.22), Lemma 5.4 and Lebesgue’s convergence Theorem it follows that
(5.37)
We use now (5.36), (5.37) and Lemma 5.4 to see that
(5.38)
Next, by (4.21) and (5.12) we obtain
and, using (4.1), (4.13), (4.3) and (4.20) it follows that
We use again the convergence (5.38) and Lebesque’s Theorem to find that
(5.39)
Theorem 5.1 is now a consequence of the convergences (5.38) and (5.39).
6 Numerical solutions
This section is devoted to the numerical solution of the contact problems presented in Section 3, including the numerical validation of the convergence result in Theorem 5.1. In order to avoid repetitions, we restrict ourselves to present details only on
the numerical approach of Problem , which is based on penalization and the augmented Lagrangean method. To this end we introduce a new variational formulation of Problem , more convenient for the numerical solution.
An adapted variational formulation. We consider the space with its usual norm and denote by and the dual of and the duality pairing mapping, respectively. We also consider the function and the operators defined by
(6.1)
where represents the indicator function of the set .
We note that, for all , condition (3.13) is equivalent to the subdifferential inclusion
(6.2)
where denotes the subdifferential of . This inclusion suggests to introduce a new unknown of the problem, the Lagrange multiplier, which represents the normal stress on the contact surface. Thus, proceeding in a standard way and using the inclusion (6.2) we obtain the following variational formulation of Problem , in terms of three unknown fields.
Problem . Find a displacement field , a stress field and a Lagrange multiplier such that, for all ,
(6.3)
(6.4)
(6.5)
Fully-discrete approximation. Let be the time step size and define
where is a sufficiently large integer. Below, for a continuous function with values in a function space, we use the notation , for . Assume that is a polyhedral domain. Moreover, let be a regular family of triangular finite element partitions of that are compatible with the boundary decomposition
, i.e. if one side of an element has more than one point on , then the side lies entirely in or . The space is approximated by the finite dimensional space of continuous and piecewise affine functions, that is,
(6.6)
where represents the space of polynomials of degree less or equal to one in . The space is approximated by the finite element space of piecewise constants, denoted . For any , we denote by its finite element projection onto , that is
We also consider the discrete space related to the discretization of the Lagrange multiplier . see [6, 7] for considerations about the discretization step.
Let and be the finite element approximations of and , respectively. Then, we consider the following fully discrete numerical approximation of Problem .
Problem . Find a discrete displacement field , a discrete stress field and a discrete Lagrange multiplier such that, for all ,
(6.7)
(6.8)
(6.9)
Note that the sum in (6.7) corresponds to the approximation of the integral in (6.3) by using a rectangle method (Top-left corner approximation) for the time integration.
Furthermore, in (6.9) we propose to approximate the operator by using a trapezoidale rule for the time integral which appears in (6.1). The approximated operator is defined as follows:
(6.10)
Here and below we use the short-hand notation .
Numerical method. In the rest of this subsection, to simplify the notation, we skip the dependence of various variables with respect to the discretization parameters , and , i.e., for example, we write instead of .
For the numerical treatment of the condition (6.9), we use the penalized method for the compliance contact combined with the augmented Lagrangean approach for the unilateral condition. To this end, we consider additional fictitious nodes for the Lagrange multiplier in the initial mesh. The construction of these nodes depends on the contact element used for the geometrical discretization of the interface . In the case of the numerical example presented below, the discretization is based on "node-to-rigid" contact element, which is composed by one node of and one Lagrange multiplier node. This contact interface discretization is characterized by a finite dimensional subspace . Let be the total number of nodes and denote by the basis functions used to define the space for . Moreover, let represent the number of nodes on the interface and let be the shape functions of the finite element space , for , i.e.
Usually, if a finite element method is used for the displacement, then a finite element method is considered for the multipliers. The expression of functions and is given by
(6.11)
where represents the value of the corresponding functions at the -th node of . Also, denotes the value of the function at the -th node of the contact element of the discretized contact interface. More details about this discretization step can be found in [1, 6, 7, 17].
It can be shown that the numerical approach of Problem is governed at each time step by a system of nonlinear equations of the form
(6.12)
where the functions and are defined below. Here the unknowns are the discrete displacement field and the Lagrange multiplier generalized vector , defined by
(6.13)
where represents the value of the corresponding function at the -th node of . Also, denotes the value of the corresponding function at the -th node of the contact element of the discretized contact interface. The generalized elastic term is defined by , where is the zero element of denotes the term given by
is defined by (6.11) and is related to by the discrete constitutive law (6.7). The contact operator , which allows to deal with the contact condition (6.9) is defined by
(6.14)
Here are derivable functions such that on and is the function defined by
and represents the gradient operator with respect the variable ; also, denotes the augmented Lagrangean functional given by
(6.15)
where is a positive penalty coefficient.
The solution of the nonlinear system (6.12) is based on a generalized Newton method, which permits to treat simultaneously the two unknowns and . Nevertheless, to keep this paper in a reasonable length, we skip the presentation of the numerical algorithm and we pass in what follows to description of the numerical example. Details on this kind of algorithms can be found in [1, 10, 17.
7 Numerical simulations
Physical setting of the numerical example. For the numerical simulations we consider the physical setting depicted in Figure 1. There, with . The domain represents the cross section of a three-dimensional deformable body subjected to the action of tractions in such a way that a plane stress hypothesis is assumed. On the part the body is clamped and, therefore, the displacement field vanishes there. Vertical tractions act on the part of the boundary and the part is traction free. No body forces are assumed to act on the body during the process. The body is in frictionless contact with an obstacle on the part of the boundary. For the discretization we use 7935 elastic finite elements and 129 contact elements.
The total number of degrees of freedom is equal to 8326 and we take a time step equal to .
We model the material’s behavior with a constitutive law of the form (1.1) in which elasticity tensor satisfies
(7.1)
where is the Young modulus, the Poisson ratio of the material and denotes the Kronecker symbol.
Moreover, in order to facilitate the numerical implementation, we assume that , where the tensor satisfies
(7.2)
For the computation below we use the following data:
Numerical results. The main purpose of this part consists to present a numerical validation of the theoretical convergence result obtained in Theorem 5.1. Our results are presented in Figures 2-6 and are described in what follows.
First, the deformed configuration as well as the contact interface forces at are plotted in Figure 2, which corresponds to the numerical solution of problem .
In order to compare the deformed mesh related to Problem with those obtained for the numerical solution of problem , we plotted in Figures 3 and 4 , respectively, the deformed configurations for the numerical solution of problems with memory term (in which the function ) and without memory term ( ), respectively. Then, in Figures 3 and 4, we note that the penetration of the contact nodes is no longer restricted by unilateral constraint and exceed the limit g. Moreover, the absence of the memory term leads to larger penetrations in the foundation.
In Figure 5 we present the evolution of the convergence of the discrete solution of the problem to the discrete solution of the problem as the deformability of the foundation tends to zero. More precisely, we plot 4 deformed meshes and the associated contact forces for 4 values of which represents here the stiffness
of the foundation after the limit is reached. One can see that for all the contact nodes are in strong penetration, whereas for two-third of the nodes slightly exceed the limit and will come into a unilateral contact.
For the convergence result, we denote by and the discrete solution of the contact problems and , respectively, for a given . The numerical estimations of the difference
at the time , for various values of the coefficient , are presented in Figure 6. It results from here that this difference converges to zero as tends towards infinity. We conclude that our results in Figure 6 represent a numerical validation of the theoretical convergence result obtained in Theorem 5.1.
References
[1] P. Alart and A. Curnier, A mixed formulation for frictional contact problems prone to Newton like solution methods, Computer Methods in Applied Mechanics
and Engineering 92 (1991), 353-375.
[2] M. Barboteu, A. Matei, M. Sofonea, Analysis of quasistatic viscoplastic contact problems with normal compliance, Quarterly Journal of Mechanics and Applied Mathematics, 65 (2012), 555-579.
[3] N. Cristescu and I. Suliciu, Viscoplasticity, Martinus Nijhoff Publishers, Editura Tehnica, Bucharest, 1982.
[4] W. Han and M. Sofonea, Quasistatic Contact Problems in Viscoelasticity and Viscoplasticity, Studies in Advanced Mathematics 30, American Mathematical Society-International Press, 2002.
[5] I.R. Ionescu and M. Sofonea, Functional and Numerical Methods in Viscoplasticity, Oxford University Press, Oxford, 1993.
[6] H. B. Khenous, P. Laborde and Y. Renard, On the discretization of contact problems in elastodynamics, Lecture Notes in Applied Computational Mechanics 27 (2006), 31-38.
[7] H. B. Khenous, J. Pommier, Y. Renard, Hybrid discretization of the Signorini problem with Coulomb friction. Theoretical aspects and comparison of some numerical solvers, Applied Numerical Mathematics 56 (2006), 163-192.
[8] A. Klarbring, A. Mikelič and M. Shillor, Frictional contact problems with normal compliance, Int. J. Engng. Sci. 26 (1988), 811-832.
[9] A. Klarbring, A. Mikelič and M. Shillor, On friction problems with normal compliance, Nonlinear Analysis 13 (1989), 935-955.
[10] T. Laursen, Computational Contact and Impact Mechanics, Springer, Berlin, 2002.
[11] J. A. C. Martins and J. T. Oden, Existence and uniqueness results for dynamic contact problems with nonlinear normal and friction interface laws, Nonlinear Analysis : Theory, Methods and Applications 11 (1987), 407-428.
[12] J.T. Oden and J.A.C. Martins, Models and computational methods for dynamic friction phenomena, Computer Methods in Applied Mechanics and Engineering 52 (1985), 527-634.
[13] M. Shillor, M. Sofonea and J.J. Telega, Models and Analysis of Quasistatic Contact, Lecture Notes in Physics, 655, Springer, Berlin, 2004.
[14] M. Sofonea, C. Avramescu and A. Matei, A Fixed point result with applications in the study of viscoplastic frictionless contact problems, Communications on Pure and Applied Analysis 7 (2008), 645-658.
[15] M. Sofonea and A. Matei, History-dependent quasivariational inequalities arising in Contact Mechanics, European Journal of Applied Mathematics, 22 (2011), 471-491.
[16] M. Sofonea and F. Pătrulescu, Analysis of a history-dependent frictionless contact problem, Mathematics and Mechanics of Solids, 18, no. 4 (2013), 409-430.
[17] P. Wriggers, Computational Contact Mechanics, Wiley, Chichester, 2002.