Variational and Numerical Analysis of a Static Thermo-Electro-Elastic Problem with Friction

Baiz, Othman; Benaissa, Hicham; El Moutawakil, Driss; Fakhar, Rachid

doi:https://doi.org/10.1155/2018/8346891

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusion Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2018 | Article ID 8346891 | https://doi.org/10.1155/2018/8346891

Variational and Numerical Analysis of a Static Thermo-Electro-Elastic Problem with Friction

Othman Baiz,¹Hicham Benaissa,²Driss El Moutawakil,¹and Rachid Fakhar²

Academic Editor: Francesco Marotti de Sciarra

Received12 Aug 2017

Revised21 Nov 2017

Accepted27 Dec 2017

Published28 Jan 2018

Abstract

We consider a mathematical model which describes a static frictional contact between a piezoelectric body and a thermally conductive obstacle. The constitutive law is supposed to be thermo-electro-elastic and the contact is modeled with normal compliance and a version of Coulomb’s friction law. We derive a variational formulation of the problem and we prove the existence and uniqueness of its solution. The proof is based on some results of elliptic variational inequalities and fixed point arguments. Furthermore, a finite element approximation and a priori error estimates are obtained.

1. Introduction

Contact problems involving thermopiezoelectricity arise when there is an interaction between the mechanical, electrical, and thermal properties of the considered material. This kind of problems has a considerable interest in various fields of modern industries. Indeed, the thermal effects such as thermal deformations and pyroelectric effects are important to many smart ceramic materials. Thus, for some materials, it may be impossible to predict the electromechanical behavior without taking into account the thermal effects; see for example, [1–4]. There is a considerable interest in contact problems involving piezoelectric materials when thermal effects are considered. These so-called thermopiezoelectric materials can operate effectively as distributed sensors and actuators for controlling structural response. In sensor applications, mechanically or thermally induced disturbances can be determined from measurement of the induced electric potential difference (direct piezoelectric effect), whereas in actuator applications deformation or stress can be controlled through the introduction of an appropriate electric potential difference (converse piezoelectric effect). Among the numerous applications of these materials, we cite accelerometers, microphones, ultrasonic transducers, and so on. A good example is the use of this class of materials as sensors and actuators in microelectromechanical systems, for instance, the piezoelectric accelerometer which triggers an airbag in ten thousandths of a second during an accident.

Some general models and their analysis have been established for elastic-piezoelectric materials [5–8], for thermoelastic materials [9–11], and for electro/thermoviscoelastic bodies [12, 13]. The mathematical treatment of static contact process for thermopiezoelectric materials is recent. The reason lies in the considerable difficulties of the nonlinear evolutionary inequalities modeling the static contact problems present in the variational analysis. Existence and uniqueness results in the study of static contact problems can be found, for instance, in [14–16].

The present paper is a continuation of this kind of models. It deals with a new and nonstandard mathematical model, which is in a form of a system coupling a nonlinear variational inequality for the displacement field and two nonlinear variational equations for the electric potential and the temperature field. The boundary conditions on the contact surface used in this paper is described with a general normal compliance law and the associated Coulomb’s friction law, taking into account the electrical and thermal conductivity of the foundation. This work serves two purposes. The first purpose is to provide the variational analysis of the mechanical problems and to show the existence of a unique solution to each model. The second is to study a discrete scheme based on finite element method for numerical solutions and to establish the unique solvability of the scheme and derive an optimal order error estimate.

The rest of the paper is structured as follows. In Section 2, we present the model of our frictional thermopiezoelectric contact problem and we derive its variational formulation, given as a coupled system for the displacement field, the electric potential, and the temperature fields. In Section 3, we provide the assumptions on the data and we state the main existence and uniqueness theorem together with its proof. In Section 4, we present the main error estimate results for the finite element approximation of the problem.

2. The Mathematical Model

We consider a thermopiezoelectric body that occupies a bounded domain in () with a Lipschitz continuous boundary . In the sequel, we decompose into three open and disjoint parts , , and such that and on the one hand and we consider a partition of into two sets and with disjoint, relatively open sets and such that on the other hand.

We assume that the body is fixed on where the displacement field vanishes and it is subject to a volume forces of density , a volume electric charges of density , and a heat source term per unit volume in . We assume also that a surface tractions of density and a strength of heat source act on and a surface electric charges of density acts on . Finally, we suppose that the electrical potential vanishes in and the temperature is maintained constant at the part of the boundary, set to be . Over the contact surface , the body comes in a frictional contact with a thermally conductive foundation.

Here and below, the indices , , , and run between and , the summation convention over repeated indices is adopted, and the index that follows a comma indicates a partial derivative with respect to the corresponding component of the spatial variable; for example, . In the sequel, let be the space of second order symmetric tensors on . The canonical inner products and norms on and are given by

To present the mathematical model which describes the physical setting above, we denote by , , , , , and the displacement field, the stress tensor, the electric potential field, the electric displacement field, the temperature field, and the heat flux vector, respectively. These are functions which depend on the spatial variable . Nevertheless, in what follows we do not indicate explicitly the dependence of these quantities on ; that is, we write instead of . Also, denotes the linearized strain tensor and is the electric field.

We assume that the process is static, then the equations of stress equilibrium, the equation of the quasistationary electric field, and the heat conduction equation arewhere and are the divergence operator for tensor and vector field. The material is assumed to be thermopiezoelectric and satisfies the following constitutive laws [1, 4]:where is the elasticity tensor, is the third-order piezoelectric tensor, represents the electric permittivity tensor, is the thermal expansion, and denotes the pyroelectric tensor. Moreover, is the transpose of and it satisfiesFor the heat flux, we adopt the following Fourier-type law:where denotes the thermal conductivity tensor. We use classical decomposition in the normal and the tangential components of the displacement and of the stress on ; that is,where is the outward unit normal vector on and the physical setting to complete our model with the following boundary conditions:We model the frictional contact on with the following reduced normal compliance condition:where , , , and are material interface parameters and represents the penetration approach; for more detail see [17–19]. Here (15) is the normal compliance power law and (16) is a variant of Coulomb’s friction law. Furthermore, the thermoelectric contact is described with the following regularized conditions (see [13, 14]):where is the thermal conductance function, supposed to be zero for and positive otherwise, nondecreasing, and Lipschitz continuous, and is the foundation temperature. The truncate function is defined bywhere is a large positive constant. At last, we note that condition (17) describes the fact that the foundation is supposed to be a perfect electrical insulator. Under all these conditions, the classical formulation of our problem is as follows.

Problem (). Find a displacement field , an electric potential , and a temperature field such that (2)–(7) and (9)–(18) hold.

Note that once the triplet which solve is known, then the stress tensor , the electric displacement field , and the heat flux can be obtained from (5) and (7). In order to derive the variational formulation of the problem , we need the following Hilbert spaces:endowed with the following inner products:Next, we note that, by the Sobolev trace theorem, we can define the trace of a function on such that if . For simplicity, for an element we still denote by its trace on . Let and , then the trace operator is a linear continuous operator; that is, there exists a positive constant , depending only on , such thatKeeping in mind condition (9), we introduce the closed subspace of given bySince , Korn’s inequality holds; there exists a positive constant which depends only on and such that (see, e.g., [7, 20])We define over the space , the following inner product, and its associated norm:It comes from (24) that the norms and are equivalent on and therefore is a real Hilbert space. Finally, note that from (22) and (24), we deduce that there exists a positive constant depending on , , and such thatAccording to the boundary conditions (11) and (13), the electric potential field and the temperature field are, respectively, to be found in the closed subspace and of given bySince , the Friedrichs-Poincaré inequality holds and thus there exists a positive constant depending only on and such thatWe introduce on the inner product and its associated norm . It follows from (28) that the norms and are equivalent on and therefore, the space is a real Hilbert. Since , we can prove in an analogous way that the norm associated with the inner product is equivalent to the usual norm on and thus, is a real Hilbert space. Using the Sobolev trace theorem, we get that there exists a positive constant depending only on , , and such thatand a positive constant which depends only on , , and such that

In the study of the mechanical problem , the following assumptions will be needed:(h₁)The elasticity operator , the electric permittivity tensor , and the thermal conductivity tensor satisfy the usual properties of symmetry, boundedness, and ellipticity and there exist positive constants , , and such that(h₂)The piezoelectric tensor , the thermal expansion tensor , and the pyroelectric tensor satisfy(h₃)The thermal conductance satisfies the following conditions:(h₄)The function is a Lipschitz function on for all ; that is,(h₅)The forces, the traction, the charges densities, and the strength of the heat source satisfy(h₆)The foundation temperature satisfies(h₇)The friction bound function and the coefficient of friction satisfy(h₈)The material interface parameters satisfyUsing the Riesz representation theorem, we can define , , and as follows:We consider the functionals and defined byIt follows from (h₃) and (h₅)–(h₈) that the integrals above are well defined. We note that if and are sufficiently regular, the following Greens formulas hold:Using the previous Greens formulas, it is straightforward to see that if are sufficiently regular functions which satisfy (2)–(7) and (9)–(18), thenfor all , , and . We use (5), (7), and the notation to obtain the following variational formulation of our problem.

Problem . Find a displacement field an electric potential and a temperature field such that

3. Existence and Uniqueness of Weak Solution

The main existence and uniqueness result in the study of the problem are as follows.

Theorem 1. Assume (h₁)–(h₃) and (h₅)–(h₈) hold. Then one has the following:(1)The problem has at least one solution.(2)Under (h₄), there exists a constant such that if then the problem has a unique solution.Here the norms of the tensors and are given by

The proof will be carried out in several steps and it is based on arguments of variational inequalities and fixed point techniques. We assume in what follows that (h₁)–(h₃) and (h₅)–(h₈) hold and for every where , ≥ , we defineFor all supposed to be known, we consider the following auxiliary problem.

Problem . Find the elements , , and such thatIn the study of this problem, the two following results will be needed.

Lemma 2. Let be a positive real number. For all , one hasand consequently for all , one has

Proof. This lemma can be obtained by examining the two cases and .

Theorem 3. Let be open bounded set of with a Lipschitz boundary and . The trace operator satisfies the following results:(1)If , the map is compact for any . Then, there exists such that(2)If , the map is compact for any and then, there exists such that

Proof. For the proof of the trace Theorem 3, we can refer to [21, 22].

Lemma 4. Under the assumptions (h₁)–(h₃) and (h₅)–(h₈), the problem has a unique solution which depends Lipschitz continuously on .

Proof. We use Riesz’s representation theorem to define the element of and the operator such thatThus equation (55) will beUsing the assumptions of (see (h₁)), we can deduce that is a linear symmetric and positive definite operator. Hence, is linear continuous and invertible operator on and let denote its inverse. Thus, by the Lax-Milgram theorem, we get that problem (62) has a unique solutionIt follows from the properties of the operators , , and that and are linear continuous operators. Moreover, we apply Riesz’s representation theorem to define the element and the operators and as follows:Let be the adjoint operator of and then, it comes from (66) thatNext, we replace (64), (65), and (66) in (54) to obtainKeeping in mind the properties of (see (h₁)), we get that is a linear symmetric and positive definite operator. Hence, the operator is invertible and let be its inverse. We haveIt comes from (63), (69) that inequality (53) is equivalent to find such thatThe variational problem above is equivalent to the following minimization problem:where the functional is defined as follows:We consider the function whose derivative isApplying Lemma 2, we get for all and of thatHence the function is convex. In plus, it follows from the strict convexity of that the functional is strictly convex on . Moreover, since is coercive and , we deduce that is coercive. Consequently, the minimization problem has a unique solution . Therefore, keeping in mind (63) and (69), we conclude that the variational problem has a unique solution of
For the second part of Lemma 4, let and be two given elements of the reflexive space such that and . We consider , , the unique solution of the problems , , respectively. Then, the variational inequality (53) leads toWe take in the first inequality and in the second to getand by adding the two induced inequalities, we obtainMoreover, the definition of and Lemma 2 imply thatThen, we haveIt follows from (26) and the definition of thatIn addition, the variational equation (54) leads toAfter taking in (81) and in (81), we add the induced equationsUsing (79), (80), (83), (30), and (6), the strong monotonicity of , and the ellipticity of and after some algebra, we find that there exists a constant such thatFrom (60), (51), and (30), we obtainand from (61), (62), (85), and the ellipticity of , we find that there exists such thatFinally, we combine (84) and (87) to find that there exists a constant such thatHence the second part of Lemma 4 is proved.

Remark 5. The second part of Lemma 4 implies that the function where the triplet is the solution of is a continuous function from to .

Lemma 6. If the triplet is a solution of problem , then there exists a positive constant such that

Proof. Taking in the variational equation (55), we getUsing (43), (30), and the ellipticity of the operator , we findMoreover, if we take in (53) and in (54), we haveKeeping in mind (6), the ellipticity of and , the positivity of and , and the properties of and , we deduce that there exists a constant such thatWe combine the two inequalities (91) and (93) to obtainwhich finishes the proof.

Remark 7. Using the same argument as in the proof of in Lemma 6, we have that if the triplet is a solution of problem , thenwhere is the same constant as in (94).

In this step, we consider the operator defined byNow, we shall prove that the operator has a unique fixed point . For that, we need to introduce the two closed convex subsets of where and will be defined below.

Lemma 8. For a specified values of and , the operator has at least one fixed point.

Proof. Let . We have and and thenOn another hand, it follows from the definition of thatUsing (h₃), (h₈), and Theorem 3, we deduce that there exists a constant such thatSince , it becomes from Lemma 6 that there exists a constant such thatIf we choose and , we getHence is an operator of into itself. Since is a nonempty, convex, and closed subset of the reflexive space , then is weakly compact. Using the continuity of the functions and and Remark 5, we deduce that is a continuous operator. Hence, by Schauder’s fixed point theorem the operator has a fixed point.

Proof of Theorem 1.
Existence. Let be the fixed point of the operator obtained in Lemma 8. We denote by the solution of the problem for . It follows from the definition of and that is a solution of which concludes the proof of the existence part.
Uniqueness. Here, we will show that there exists a positive constant such that has a unique solution if . We consider and two solutions of . From (46), we haveAfter taking in the first inequality and in the second and adding the resulting inequalities, it follows from the positivity of thatPlus the variational equation (47) leads toWe take in the first equation and in the second to obtainHence, the addition of (104) and (106) implies thatMoreover, we use the variational equation (48) to deduceWe substitute by and we subtract the two induced equations to obtainLet us consider . Then, we havewhereTaking in mind (26), (30), and the assumptions (h₃) and (h₄), we getHence, it comes from (109), (110), and (112) and the ellipticity of thatRecalling , , we deduceFurthermore, the definition of implies thatUsing the following mathematical inequalities (see [23])Then, we obtainKeeping in mind Theorem 3 and Remark 7, we obtainwhere .
Using (107), (114), and (118), after some algebra it follows that there exists such thatwhereChoose . Then, if holds, we will conclude that , , and which leads to the uniqueness part of Theorem 1.

4. Discrete Approximation

This section deals with the discrete approximation of the problem . We assume that the conditions (h₁)–(h₈) hold, then the problem has a unique solution × . Let be a family of regular triangulations of the polygonal domain such thatHere and below is a discretization parameter. We define the following finite dimensional subspaces , , and which approximate, respectively, the spaces , , and bywhere denotes the space of polynomials of a degree lower or equal to one on . Then, the discrete approximation of the problem is as follows.

Problem (). Find the displacement , the electric potential , and the temperature such thatUnder the assumptions of Theorem 1 and with the same arguments, we can prove that the discrete problem has a unique solution in . Now, we proceed to derive some error estimates for the discrete solution. In the sequel, denotes positive constants which are independent of the discretization parameter .

Theorem 9. Assume the conditions of Theorem 1 hold. One has the following error estimate:

Proof. Taking in the second equation of the problem and subtracting the obtained equation to the second equation of the problem , we obtain We substitute by to find after simplificationsTake in the first inequality of the problem to deduceBy using , it follows from the first inequality of thatFrom Lemma 2, we know that . Then, keeping in mind t (128) and (6), the addition of (129) and (130) giveswhereNext, we rewrite (48) for and we subtract it from (125) to obtainWe replace by . Hence, we haveBy adding the two results (134) and (131), we deduce the following inequality:where and are given byNow, we begin to estimate the quantities , , , and . We haveWe recall that there exists a positive constant independent of such thatand similar to Remark 7, we can find that . ThusFinally, remembering (95), we find that there exists a constant such thatUsing the same process as in the proof of (118), we can deducewhere is a positive constant. Moreover, we havewhere is a positive constant. Finally, we havewhere is a positive constant. Keeping in mind the following -inequality , it follows from (140)–(143) that if , then there exists a positive constant such thatand that finishes the proof of the Theorem 9.

As a result of the previous theorem, the following corollary is about the convergence order error estimates for the fully discrete approximations with the previous subspaces , , and . Let be, as usual, the interpolation operator . Using the standard finite element interpolation error estimates, we have the following approximations:

Corollary 10. Assume that conditions of Theorem 9 hold. Under the regularity conditionsthere exists positive constant , independent of , such that

This corollary gives an estimation of the numerical errors of the problem and its proof is based on the above approximation properties of the finite element spaces , , and .

Now, we investigate the particular case where and . Indeed, under some conditions on and , we can prove the following result.

Theorem 11. Let and . One assumes that, for some , one hasThen, the following error estimate can be obtained:

Proof. Since is zero, the functional becomesThen, the variational inequalities (46) and (123) can be writtenWe take in (151) to obtainand from (152), we concludeAdding (153) and (154), we deduceWe recall that and , thenRemembering (128) and (6), we replace by its value. HenceIt follows from the sum of (134) and (157) that, for all of , we haveTo simplify the calculations, let us consider the following quantities:Taking in the integral term of (159), it comes from (26), (30), and the assumptions – thatUsing the Green formula and recalling (2), (9), (10), (15), and (16), we obtainand since , we haveHenceMoreover, taking in the first integral of (161) and using the Lemma 2, we haveand then we concludeHowever, we take in the integral term of (162) and we use conditions ()– and (26), (30) to findwhere . If we suppose withthenKeeping in mind and (158), (159), (160), (161), (162), (163), (166), (168), (169), and (171) we apply several times the -inequality to conclude that there exists a constant such that, for all of , we have

5. Conclusion

In this work, we presented a model for the static process of frictional contact between a piezoelectric body and an electrically thermally conductive foundation. The constitutive relation of the material is assumed to be thermo-electro-elastic. The contact was modeled with the normal compliance condition and the associated Coulomb’s friction law, including the electrical and thermal conductivity conditions. The existence of the unique weak solution for the problem was established by using arguments from the theory of variational inequalities and a fixed point theorem. A discrete scheme by finite element method was used to approach the problem and an optimal order error estimate was derived. A numerical validation of the convergence result included in this method will be provided in a forthcoming paper.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

M. Aouadi, “Generalized thermo-piezoelectric problems with temperature-dependent properties,” International Journal of Solids and Structures, vol. 43, no. 21, pp. 6347–6358, 2006.
View at: Publisher Site | Google Scholar
D. S. Chandrasekharaiah, “A generalized linear thermoelasticity theory for piezoelectric media,” Acta Mechanica, vol. 71, no. 1-4, pp. 39–49, 1988.
View at: Publisher Site | Google Scholar
W. Nowacki, “Some general theorems of thermo-piezoelectricity,” Journal of Thermal Stresses, vol. 1, no. 2, pp. 171–182, 1978.
View at: Publisher Site | Google Scholar
H. F. Tiersten, “On the nonlinear equations of thermoelectroelasticity,” International Journal of Engineering Science, vol. 9, pp. 587–604, 1971.
View at: Publisher Site | Google Scholar | MathSciNet
E.-H. Essoufi, E.-H. Benkhira, and R. Fakhar, “Analysis and numerical approximation of an electro-elastic frictional contact problem,” Advances in Applied Mathematics and Mechanics, vol. 2, no. 3, pp. 355–378, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
P. Bisegna, F. Lebon, and F. Maceri, “The unilateral frictional contact of a piezoelectric body with a rigid support,” in Contact Mechanics, vol. 103, pp. 347–354, Kluwer, Dordrecht, The Netherlands, 2002.
View at: Publisher Site | Google Scholar | MathSciNet
G. Duvaut and J.-L. Lions, Inequalities in Mechanics and Physics, Springer, Berlin, Germany, 1976.
View at: MathSciNet
M. Sofonea and A. Matei, Mathematical Models in Contact Mechanics, Lecture Note, Series 398, Cambridge University Press, Cambridge, UK, 2012.
View at: Publisher Site
L.-E. Andersson, A. Klarbring, J. R. Barber, and M. Ciavarella, “On the existence and uniqueness of steady state solutions in thermoelastic contact with frictional heating,” Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 461, no. 2057, pp. 1261–1282, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
G. Duvaut, “Free boundary problem connected with thermoelasticity and unilateral contact,” in Free boundary problems, pp. 217–236, 1980.
View at: Google Scholar | MathSciNet
M. A. Ezzat and E. S. Awad, “Constitutive relations, uniqueness of solution, and thermal shock application in the linear theory of micropolar generalized thermoelasticity involving two temperatures,” Journal of Thermal Stresses, vol. 33, no. 3, pp. 226–250, 2010.
View at: Publisher Site | Google Scholar
A. S. El-Karamany and M. A. Ezzat, “On the boundary integral formulation of thermo-viscoelasticity theory,” International Journal of Engineering Science, vol. 40, no. 17, pp. 1943–1956, 2002.
View at: Publisher Site | Google Scholar | MathSciNet
Z. Lerguet, M. Shillor, and M. Sofonea, “A frictional contact problem for an electro-viscoelastic body,” Electronic Journal of Differential Equations, vol. 170, pp. 1–16, 2007.
View at: Google Scholar | MathSciNet
H. Benaissa, E.-H. Essoufi, and R. Fakhar, “Existence results for unilateral contact problem with friction of thermo-electro-elasticity,” Applied Mathematics and Mechanics-English Edition, vol. 36, no. 7, pp. 911–926, 2015.
View at: Publisher Site | Google Scholar | MathSciNet
H. Benaissa, E. L.-H. Essoufi, and R. Fakhar, “Analysis of a Signorini problem with nonlocal friction in thermo-piezoelectricity,” Glasnik Matematički, vol. 51, no. 2, pp. 391–411, 2016.
View at: Publisher Site | Google Scholar
H. Benaissa, E.-H. Essoufi, and R. Fakhar, “Variational analysis of thermo-piezoelectric contact problem with friction,” Journal of Advanced Research in Applied Mathematics, vol. 7, no. 2, pp. 52–75, 2015.
View at: Publisher Site | Google Scholar
A. Klarbring, A. Mikelic, and M. Shillor, “On friction problems with normal compliance,” Nonlinear Analysis: Theory, Methods & Applications, vol. 13, no. 8, pp. 935–955, 1989.
View at: Publisher Site | Google Scholar | MathSciNet
J. A. C. Martins and J. T. Oden, “Existence and uniqueness results for dynamic contact problems with nonlineaire normal and friction interface laws,” Nonlineaire Analysis. Theory, Methods & Applications, vol. 11, no. 3, pp. 407–428, 1986.
View at: Google Scholar
P. Hild, “Two results on solution uniqueness and multiplicity for the linear elastic friction problem with normal compliance,” Nonlinear Analysis: Theory, Methods & Applications, vol. 71, no. 11, pp. 5560–5571, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
J. Nečas and I. Hlaváček, Mathematical Theory of Elastic and Elasto-Plastic Bodies: An Introduction, Elsevier Scientific, Amsterdam, The Netherlands, 1981.
View at: MathSciNet
R. A. Adams, Sobolev Spaces, Academic Press, NY, USA, 1975.
View at: MathSciNet
S. a. Migórski, A. Ochal, and M. Sofonea, Nonlinear inclusions and hemivariational inequalities: Models and analysis of contact problems, Springer, NY, USA, 2014.
View at: MathSciNet
G. H. Hardy, J. E. Littlewood, and G. Pólya, Inequalities, Cambridge University Press, London, UK, 2nd edition, 1952.
View at: MathSciNet

Copyright

Copyright © 2018 Othman Baiz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

854

Downloads

823

Citations