Next Article in Journal
A Comprehensive Multi-Strategy Enhanced Biogeography-Based Optimization Algorithm for High-Dimensional Optimization and Engineering Design Problems
Previous Article in Journal
Meeting of Chemical Wave Fronts Creates Chaos
Previous Article in Special Issue
The Integrability and Modification to an Auxiliary Function Method for Solving the Strain Wave Equation of a Flexible Rod with a Finite Deformation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence and Uniqueness of Weak Solutions to Frictionless-Antiplane Contact Problems

1
Laboratory of Differential Equations, Department of Mathematics, University of Constantine 1, Ain El Bey Road, Constantine P.O. Box 325, Algeria
2
Laboratory of Applied Mathematics and Modeling, Department of Mathematics, University of Constantine 1, Ain El Bey Road, Constantine P.O. Box 325, Algeria
3
Center for Research and Development in Mathematics and Applications (CIDMA), Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal
4
Faculty of Business and Communications, INTI International University, Persiaran Perdana BBN, Putra Nilai, Nilai 71800, Negeri Sembilan, Malaysia
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(3), 434; https://doi.org/10.3390/math12030434
Submission received: 17 November 2023 / Revised: 22 January 2024 / Accepted: 26 January 2024 / Published: 29 January 2024

Abstract

:
We investigate a quasi-static-antiplane contact problem, examining a thermo-electro-visco-elastic material with a friction law dependent on the slip rate, assuming that the foundation is electrically conductive. The mechanical problem is represented by a system of partial differential equations, and establishing its solution involves several key steps. Initially, we obtain a variational formulation of the model, which comprises three systems: a hemivariational inequality, an elliptic equation, and a parabolic equation. Subsequently, we demonstrate the existence of a unique weak solution to the model. The proof relies on various arguments, including those related to evolutionary inequalities, techniques for decoupling unknowns, and certain results from differential equations.

1. Introduction

Over the past few decades, phenomena related to the contact between deformable bodies have had a significant impact on society. The contact of brakes with wheels, a ball with the ground, and the application of a force to cutting objects, are just a few everyday instances, among many other examples. Consequently, research in this area is growing across various fields, particularly in engineering and mathematical literature [1,2].
A comprehensive study of mechanical problems involves mathematical modeling. The modeling of mechanical phenomena is determined by a set of hypotheses that influence the system of partial differential equations. These hypotheses cover aspects such as the nature of the mechanical process (static, quasi-static, or dynamic), the behavior of the material (electromechanical, thermomechanical, elastic, etc.), and the boundary conditions on the contact surface (friction, adhesion, etc.) [3,4].
In 1933, notable advancements in the mathematical and mechanical exploration of contact-mechanics problems took place. Signorini seems to have been the first to delve into this subject when he formulated the challenge of contact between a deformable body and a foundation. The solution to this problem eventually emerged in 1964, courtesy of Fichera, who employed emerging mathematical techniques. The initial findings concerning the existence and uniqueness of contact problems were established by Duvaut and Lions. Following this, a multitude of subsequent works have concentrated on resolving these variational problems [5,6].
Mathematics plays a pivotal role in the field of contact mechanics by providing diverse contact conditions. Contemporary research in this field also encompasses behavior laws that establish connections between mechanical and electrical effects, known as piezoelectricity [7,8]. This area is extensively explored in engineering structures, due to its temperature-dependent characteristics. Several articles have tackled the thermo-piezoelectric contact problem with friction [5,9,10,11,12,13], thermo-elasto-visco-plasticity [14,15,16], and thermo-visco-elasticity [17,18]. These studies are rooted in variations in constitutive laws and contact conditions [19,20,21,22,23,24,25].
Here we propose a novel mathematical model designed to address the frictional-antiplane contact problem between a thermo-piezoelectric body and a conductive foundation. The innovation within our model lies in the introduction of a novel frictional contact tailored for such materials. This model incorporates a slip-rate-dependent friction law and introduces modifications to the electrical and thermal conditions used to describe the contact. We specifically focus on the case of antiplane-shear deformation.
The manuscript is organized as follows. In Section 2, we recall some notations and definitions used in the sequel. In Section 3, we introduce our mathematical model for the quasi-static thermo-electro-visco-elastic-antiplane contact problem. Given some consistency assumptions about the data, the variational formulation of the model is then obtained in Section 4. Finally, in Section 5, we prove the existence and uniqueness result of the formulated problem.

2. Preliminaries

In this brief section, we review the fundamental notations and definitions that are employed throughout the work. For more details, we refer the interested reader to [25].
The space of second-order symmetric tensors on R 3 is denoted by S 3 . The inner product and the norm on R 3 are given as follows:
u h · v h = u i h v i h , v h = ( v h · v h ) 1 2 for all u h = ( u i h ) , v h = ( v i h ) R 3 , 1 i , j 3 .
The inner product on S 3 is defined by
σ h · τ h = σ i , j h τ i , j h , for all σ h = ( σ i , j h ) , τ h = τ i , j h S 3 , 1 i , j 3 ,
while its corresponding norm is given by
τ h = τ h · τ h 1 2 .
Next, we consider the following function spaces:
V h = { v h H 1 ( Ω h ) , v h = 0 on Γ 1 h } ,
W h = { ψ h H 1 ( Ω h ) , ψ h = 0 on Γ 1 h } ,
E h = { θ h H 1 Ω h , θ h = 0 on Γ 1 h Γ 2 h } .
It is an established fact that V h , W h , and E h are real Hilbert spaces with the inner products
( u h , v h ) V h = Ω h u h · v h d x , u h , v h V h ,
( φ h , ψ h ) W h = Ω h φ h · ψ h d x , φ h , ψ h W h ,
( θ h , ϕ h ) E h = Ω h θ h · ϕ h d x , θ h , ϕ h E h .
Furthermore, the associated norms
v h V h = v h L 2 ( Ω h ) 2 , v h V h ,
ψ h W h = ψ h L 2 ( Ω h ) 2 , ψ h W h ,
θ h E h = θ h L 2 ( Ω h ) 2 , θ h E h
on V h , W h , and E h are equivalent to the usual norm · H 1 ( Ω h ) . In the light of Sobolev’s trace theorem, we deduce that
c V h > 0 : v h L 2 ( Γ 3 h ) c V h v h V h , v h V h ,
c W h > 0 : v h L 2 ( Γ 3 h ) c W h ψ h W h , ψ W h ,
c E h > 0 : v h L 2 ( Γ 3 h ) c E h θ h E h , θ E h .
In a real Banach space ( X h , · X h ) , we denote the spaces of continuous and continuously differentiable functions on [ 0 , T ] into X h by C ( 0 , T , X h ) and C 1 ( 0 , T , X h ) , respectively. These spaces are equipped with their respective norms,
x h C ( 0 , T , X h ) = max t [ 0 , T ] x h X h
and
x h C 1 ( 0 , T , X h ) = max t [ 0 , T ] x h X h + max t [ 0 , T ] x ˙ h X h .
Below, we shall use the notations for the Lebesgue space L 2 ( 0 , T , X h ) as well as the Sobolev space W 1 , 2 ( 0 , T , X h ) . Recall that the norm on L 2 ( 0 , T , X h ) is given by
u h L 2 ( 0 , T , X h ) 2 = 0 T u h ( t ) X h 2 d t .
Moreover,
u h W 1 , 2 ( 0 , T , X h ) 2 = 0 T u h ( t ) X h 2 d t + 0 T u ˙ h ( t ) X h 2 d t
defines a norm on W 1 , 2 ( 0 , T , X h ) .

3. Formulation of the Problem

We begin by making explicit the physical setting under investigation. Let us consider a thermo-electro-visco-elastic body B h with friction. Assume that B h R 3 is a cylinder with generators parallel to the x 3 h -axes having a regular region in its cross-section Ω h in the x 1 h , x 2 h -plane, O x 1 h x 2 h x 3 h being a Cartesian-coordinate system. The cylinder is presumptively long enough for the end effects in the axial direction to be insignificant. Thus, B h = Ω h × ( , + ) . Let Ω h = Γ h . We suppose that Γ h is divided into three disjoint measurable parts: Γ 1 h , Γ 2 h , and Γ 3 h . One has a partition of Γ 1 h Γ 2 h into two open parts, Γ a h and Γ b h , such that the measures of Γ a h and Γ b h are positive. For T > 0 , we denote by [ 0 , T ] the time interval of interest. The B h body is affected to time-dependent forces f 0 h and has a volume with free electric charges of density q 0 h . The body is clamped on Γ 1 h × ( , + ) and, therefore, the displacement field disappears there. The surface tractions of density f 2 h act on Γ 2 h × ( , + ) . We also assume that the electrical potential vanishes on Γ a h × ( , + ) and that a surface electrical charge of density q 2 h is prescribed on Γ b h × ( , + ) . The body is in contact with a rigid foundation Γ 3 h × ( , + ) [26,27]. We assume that
q 0 h = q 0 h ( x 1 h , x 2 h , t ) : Ω h × [ 0 , T ] R ,
q 2 h = q 2 h ( x 1 h , x 2 h , t ) : Γ b h × [ 0 , T ] R ,
f 0 h = ( 0 , 0 , f 0 h ) , with f 0 h = f 0 h ( x 1 h , x 2 h , t ) : Ω h × [ 0 , T ] R ,
f 2 h = ( 0 , 0 , f 2 h ) , with f 2 h = f 2 h ( x 1 h , x 2 h , t ) : Γ 2 h × [ 0 , T ] R .
The electric charges (7) and (8) and the forces (9) and (10) are expected to give rise to electric charges and to deformations of the piezoelectric cylinder corresponding to an electric-potential field φ h and to a displacement u h , which are independent of x 3 h and have the form
φ h = φ h ( x 1 h , x 2 h , t ) : Ω h × [ 0 , T ] R ,
u h = u h ( 0 , 0 , u h ) , with u h = u h ( x 1 h , x 2 h , t ) : Ω h × [ 0 , T ] R .
We denote by θ h the temperature field, which is of the form
θ h = θ h ( x 1 h , x 2 h , t ) : Ω h × [ 0 , T ] R ,
by ε h ( u h ) = ( ε i , j h ( u h ) ) the infinitesimal strain tensor, by σ h = ( σ i j h ) the stress field, by E h ( φ h ) = ( E i h ( φ h ) ) the electric field, and by D h = ( D i h ) the electric displacement field, where
ε i , j h ( u h ) = 1 2 ( u i , j h + u j , i h ) ,
E i h ( φ h ) = φ , i h .
The material is modeled by a thermo-electro-visco-elastic constitutive law of the form
σ h = 2 α h ε h ( u ˙ h ) + ζ h t r ε h ( u h ˙ ) I + 2 μ h ε h ( u h ) + λ h t r ε ( u h ) I E h * E h ( φ h ) M e h θ h ,
D h = E h * ε h ( u h ) + β h E h ( φ h ) ,
where ζ h and α h are viscosity coefficients, λ h and μ h are the Lame coefficients, tr ε h ( u h ) = ε i i h ( u h ) , I is the unit tensor in R 3 , β h is the electric permittivity constant, E h represents a third-order piezoelectric tensor, and E h * is its transpose. We assume that
E h ε h = e h ( ε 13 h + ε 31 h ) e h ( ε 23 h + ε 32 h ) e h ( ε 33 h ) , ε h = ( ε i , j h ) S 3 ,
where e h is the piezoelectric coefficient. Additionally, M e h = ( M i j h ) represents the thermal expansion tensor, which takes the form
M e h = 0 0 M e 1 0 0 M e 2 M e 1 M e 2 0 .
Conversely, we assume that the coefficients M e h , α h , μ h , β h , and e h depend on x 1 h and x 2 h . However, they are independent on x 3 h . Owing to E h ε h · v h = ε h · E h * v h for all ε h S 3 and v h R 3 , it follows from that (18)
E h * v h = 0 0 e h v 1 h 0 0 e h v 2 h e h v 1 h e h v 2 h e h v 3 h , v h = ( v i h ) R 3 .
In the antiplane context, having in mind (11)–(13), and given the constitutive Equations (16) and (17) and equalities (18) and (20), we obtain the stress field and the electric displacement field as
σ h = 0 0 σ 13 h 0 0 σ 23 h σ 31 h σ 32 h 0 , D h = e u , 1 h β h φ , 1 h e u , 2 h β h φ , 2 h 0 ,
where
σ 13 h = σ 31 h = α h u ˙ , 1 h + μ h u , 1 h + e h φ , 1 h M e 1 θ h ,
σ 23 h = σ 32 h = α h u ˙ , 2 h + μ h u , 2 h + e h φ , 2 h M e 2 θ h .
We presume that the process is electrically static and mechanically quasi-static. Therefore, the equilibrium equations that govern it are given by
div σ h + f 0 h = 0 , D i , i h q 0 h = 0 in B h × ( 0 , T ) ,
where div σ h = ( σ i j , j h ) is the divergence of the tensor field σ h . Accordingly, taking into account (7), (9), (11)–(13), and (21), the equilibrium Equation (22) condenses into the following scalar equations:
div ( α h u ˙ h + μ h u h ) + div ( e h φ h ) M e h div ( θ h ) + f 0 h = 0 , in Ω h × ( 0 , T ) ,
div ( e h u h ) div ( β h φ h ) = q 0 h , in Ω h × ( 0 , T ) ,
with
M e h = M e 1 M e 2 0 .
In the sequel, we use the notation
div τ h = τ 1 , 1 h + τ 1 , 2 h for τ h = ( τ 1 h ( x 1 h , x 2 h , t ) , τ 2 h ( x 1 h , x 2 h , t ) ) ,
v h = ( v , 1 h , v , 2 h ) , ν h v h = v , 1 h ν 1 h + v , 2 h ν 2 h for v h = v h ( x 1 h , x 2 h , t ) .
Keeping in mind that the cylinder is clamped on Γ 1 h × ( , + ) , the electrical potential vanishes on Γ a h × ( , + ) , and from (7) and (12) we find that
u h = 0 on Γ 1 h × ( 0 , T ) ,
φ h = 0 on Γ a h × ( 0 , T ) .
Note that ν h is the unit normal on Γ h × ( , + ) , where
ν h = ( ν 1 h , ν 2 h , 0 ) , with ν i h = ν i h ( x 1 h , x 2 h ) : Γ h R , for i = 1 , 2 .
We denote by v ν h and v τ h the normal and tangential components, respectively, of v h on the boundary—that is,
v ν h = v h · ν h , v τ h = v h v ν h ν h .
Therefore, we denote by σ ν h and σ τ h the normal and the tangential components, respectively, of σ h on the boundary, meaning that
σ ν h = ( σ h ν h ) · ν h , σ τ h = σ h ν h σ ν h ν h .
From (21) and (28), we conclude that the Cauchy stress vector and the normal component of the electric displacement field are given by
σ h ν h = ( 0 , 0 , α h ν h u ˙ h + μ h ν h u h + ν h φ h M e h θ h · ν h ) , D h · ν h = e h ν h u h β h ν h φ h ,
respectively. Following this, we utilize the following notations:
ν h u h = u 1 h ν 1 h + u 2 h ν 2 h ,
ν h φ h = φ 1 h ν 1 h + φ 2 h ν 2 h .
Taking into account the traction boundary condition on Γ 2 h × ( , + ) and the electric condition on Γ b h × ( , + ) , it follows from (8), (10), and (31) that
α h ν h u ˙ h + μ h ν h u h + e h ν h φ h M e h θ h · ν h = f 2 h on Γ 2 h × ( 0 , T ) ,
e h ν h u h β h ν h φ h = q 2 h on Γ b h × ( 0 , T ) .
We now describe the frictional contact condition on Γ 3 h × ( , + ) . First, we remark that from (12), (28), and (29), we find that u ν = 0 , which indicates that the contact is bilateral. Thus, the contact is kept during the whole process. Now, using (12) and (28)–(31), we obtain
u τ h = ( 0 , 0 , u h ) ,
σ τ h = ( 0 , 0 , σ τ h ) ,
where
σ τ h = α h ν h u ˙ h + μ h ν h u h + e h ν h φ h M e h θ h · ν h .
We suppose that the friction is invariant concerning the x 3 h axis, being modeled with a slip-rate-dependent friction law, where the strict inequality is satisfied in the stick zone and the equality in the slip zone—that is,
| σ τ h | r ( | u ˙ τ h | ) , on Γ 3 h × [ 0 , T ] , σ τ h = r ( | u ˙ h | ) u ˙ τ h | u ˙ τ h | , on Γ 3 h × [ 0 , T ] .
Here, r : Γ 3 h × R R 3 is a given function, the friction bound, and u ˙ τ h is the tangential velocity on the contact boundary: see [28,29] for details. By (34) and (36), we ascertain that the conditions in (37) are given by
| α h ν h u ˙ h + μ h ν h u h + e h ν h φ h M e h θ h · ν h | r ( | u ˙ h | ) , α h ν h u ˙ h + μ h ν h u h + e h ν h φ h M e h θ h · ν h = r ( | u ˙ h | ) u ˙ h | u ˙ h | , o n Γ 3 h × [ 0 , T ] .
Finally, we present the initial conditions:
u h ( 0 ) = u 0 h , φ h ( 0 ) = φ 0 h , θ h ( 0 ) = θ 0 h in Ω h ,
where u 0 h , φ 0 h , and θ 0 h are given functions on Ω h . Putting this all together, we obtain the problem under investigation.
Problem 1.
Determine a displacement field u h : Ω h × 0 , T R , the electric field φ h : Ω h × 0 , T R , and a temperature field θ h : Ω h × 0 , T R + , such that
div ( α h u ˙ h + μ h u h ) + div ( e h φ h ) M e h div ( θ h ) + f 0 h = 0 , in Ω h × ( 0 , T ) ,
div ( e h u h ) div ( β h φ h ) = q 0 h , in Ω h × ( 0 , T ) ,
θ ˙ h div ( k θ h ) = M e h ( u ˙ h ) + p ( t ) , in Ω h × ( 0 , T ) ,
u h = 0 , on Γ 1 h × ( 0 , T ) ,
α h ν h u ˙ h + μ h ν h u h + e h ν h φ h M e h ν h θ h = f 2 h , on Γ 2 h × ( 0 , T ) ,
| α h ν h u ˙ h + μ h ν h u h + ν h φ h M e h θ h · ν h | r ( | u ˙ h | ) , α h ν h u ˙ h + μ h ν h u h + ν h φ h M e h θ h · ν h = r ( | u ˙ h | ) u ˙ h | u ˙ h | , on Γ 3 h × ( 0 , T ) ,
φ h = 0 , on Γ a h × ( 0 , T ) ,
e h ν h u h β h ν h φ h = q 2 h , on Γ b h × ( 0 , T ) ,
θ 0 h = 0 , on Γ 1 h Γ 2 h × ( 0 , T ) ,
K i , j h θ h x j h v i h = K e h ( θ h θ R h ) , on Γ 3 h × ( 0 , T ) ,
u h ( 0 ) = u 0 h , φ h ( 0 ) = φ 0 h , θ h ( 0 ) = θ 0 h , in Ω h × ( 0 , T ) .
It should be noted that the problem we consider here is different from that of [20]. Indeed, in [20] they studied electro-visco-elastic material, meaning its behavior was influenced solely by electrical and mechanical factors. The focus was purely on the electrical impact without thermal interference, represented by equation (2.13) of [20]. By contrast, in our paper we add the factor of temperature, which plays a significant role in altering electromechanical properties through heating and cooling. For instance, the electro-rheological frictional force may vary depending on the temperature, thereby impacting the rheological properties of the material. Therefore, we can say that the material’s behavior in our study is influenced by mechanical, electrical, and thermal factors, represented by Equation (16), followed by the associated conditions boundary in the Problem 1. On the other hand, the addition of the temperature factor also provides us with an opportunity to explore and search for a different and simpler mathematical-solution method (the decoupling-of-unknowns method). Additionally, there is a difference between the friction conditions. Indeed, in [20] the authors used Tresca’s friction law (equation (2.25) there), which is characterized by setting a maximum limit for slip. It relies solely on the maximum force without considering time rates. This was what the equation described there, observing that the friction bound was constant and had one condition, making it easier to solve the posed problem. By contrast, in our paper we use the slip-rate-dependent friction law (37): here, the friction level depends on the rate of slip changes over time. This implies that variations in the slip rate impact the friction level (the friction being higher when the rate of change is greater and lower when the rate of change is smaller). This behavior is described by Equation (37), where we represent friction bound by the variable r.

4. Variational Formulation

To study our Problem 1, we assume that the viscosity coefficient, the electric permittivity coefficient, the Lame coefficient, and the piezoelectric coefficient, satisfy
α h L ( Ω h ) and there exists α h * > 0 , such that α h ( x h ) α h * , a . e . , x h Ω h ,
β h L ( Ω h ) and there exists β h * > 0 , such that β h ( x h ) β h * , a . e . , x h Ω h ,
μ h L and μ h ( x h ) > 0 , a . e . , x h Ω h ,
e h L ,
respectively. We also assume that the thermal tensors M e h = M i j : Ω h S 3 satisfy
M i j = M j i L Ω h , 1 i , j 3 .
The thermal conductivity tensor K h ˜ = K i j h : Ω h S 3 satisfies
( i ) K i j h = K j i h L Ω h , 1 i , j 3 , ( ii ) there exists m K h > 0 , such that K h s · s m K h s 2 , a . e . , x Ω h , s R 3 .
The boundary thermal data satisfy
p h W 1 , 2 ( 0 , T ; L 2 ( Ω h ) ) , θ R h W 1 , 2 ( 0 , T ; L 2 ( Γ 3 h ) ) , K e h L Ω h , R + .
The forces, tractions, volume, and surface free charge densities have the following regularity:
f 0 h W 1 , 2 ( 0 , T , L 2 ( Ω h ) ) , f 2 h W 1 , 2 ( 0 , T , L 2 ( Γ 2 h ) ) ,
q 0 h W 1 , 2 ( 0 , T , L 2 ( Ω h ) ) ,
q 2 h W 1 , 2 ( 0 , T , L 2 ( Γ b h ) ) , q 2 h = 0 , a . e . , x h Γ b h ,
respectively. The friction bound satisfies
a ) r : Γ 3 h × R R + , b ) L r > 0 , such that | r ( x h , s 1 ) r ( x h , s 2 ) | L r | s 1 s 2 | , c ) x h r ( x h , s ) is Lebesgue-measurable on Γ 3 h , s R , d ) the mapping x h r ( x h , 0 ) belongs to L 2 ( Γ 3 h ) .
We define the functional j : v h × v h R by
j ( u ˙ h , v h ) = Γ 3 h r ( | u ˙ h | ) | v h | d a , v h V h .
Let η 1 , η 2 , v 1 , v 2 V . By using (61) and (62), we find that
j ( η 1 , v 2 h ) j ( η 1 , v 1 h ) + j ( η 2 , v 1 h ) j ( η 2 , v 2 h ) |   L r η 1 η 2 L 2 ( Γ 3 h ) v 1 v 2 L 2 ( Γ 3 h ) ,
and applying the norm mentioned in (4), one obtains
L r η 1 η 2 L 2 ( Γ 3 h ) v 1 v 2 L 2 ( Γ 3 h ) c V 2 L r η 1 η 2 V h v 1 h v 2 h V h .
Hence, we conclude that
j ( η 1 , v 2 h ) j ( η 1 , v 1 h ) + j ( η 2 , v 1 h ) j ( η 2 , v 2 h ) c V 2 L r η 1 η 2 V h v 1 h v 2 h V h .
The initial data verify
u 0 h V h , θ 0 h L 2 ( Ω h ) .
We use functions f h : [ 0 , T ] V h and q h : [ 0 , T ] W h as
( f h , v h ) V h = Ω h f 0 h v h d x h + Γ 2 h f 2 h v h d a , v h V h ,
( q h , ψ h ) V h = Γ b h q 2 h ψ h d x h + Ω h q 0 h ψ h d a , ψ h W h .
The definitions of f h and q h are based on Riesz’s representation theorem. Therefore, by using assumptions (65) and (66), we ascertain that the above integrals are well defined and that
f h W 1 , 2 ( 0 , T , V h ) ,
q h W 1 , 2 ( 0 , T , W h ) .
Then, the function P : [ 0 , T ] ( E h ) and the operators K h ˜ : E ( E h ) , M h ˜ : V h ( E h ) are defined by
P ( t ) , μ h ( E h ) × E h = Ω h p h μ h d x h + Γ 3 h K e h θ R h μ h d s ,
K h ˜ τ , μ h ( E h ) × E h = Γ 3 h K e h τ μ h d s + i , j = 1 d Ω h K i , j h μ h x j h μ h x i h d x ,
M h ˜ v h , μ h ( E h ) × E h = Ω h ( M e h v h ) μ h d x + Γ 3 h h τ ( | v τ h | ) μ h d s
for all v h V h , τ E h , μ h E h . Moreover, we define the bilinear forms a α h : V h × V h R , a μ h : V h × V h R , a e h : V h × W h R , a e h : W h × V h R , a β h : W h × W h R and a M h : E h × V h R by
a α h ( u h , v h ) = Ω h α h u h · v h d x ,
a μ h ( u h , v h ) = Ω h μ h u h · v h d x ,
a M h ( u h , v h ) = Ω h M e h u h · v h d x ,
a e h ( u h , φ h ) = Ω h e h u h · v h d x = a e h ( v h , u h ) ,
a β h ( φ h , ψ h ) = Ω h β h φ h · ψ h d x
for all u h , v h V h , φ h , ψ h W h . We note that by assumptions (51)–(54), the above integrals are well defined. Using the definition of norms (1)–(5), we see that the forms a α h , a μ h , a e h , and a β h , are continuous and that the forms a α h , a μ h , and a β h are symmetric. Furthermore, the form a α h is V-elliptic, i.e.,
a α h ( u h , v h )   α h L ( Ω h ) u h V h v h V h , u h , v h V h ,
a α h ( v h , v h ) α h * v h V h 2 v h V h .
Now, we can state the variational formulation of our Problem 1.
Problem 2
(Variational Problem). Determine a displacement field u h : [ 0 , T ] V h , an electric-potential field φ h : [ 0 , T ] W h , and a temperature field θ h : [ 0 , T ] E h , such that
a α h ( u ˙ h ( t ) , v h u ˙ h ( t ) ) + a μ h ( u h ( t ) , v h u ˙ h ( t ) ) + a e h ( φ ( t ) , v u ˙ h ( t ) ) + a M h ( θ h ( t ) , v h u ˙ h ( t ) )
+ j ( u ˙ h ( t ) , v h ) j ( u ˙ h ( t ) , u ˙ h ( t ) ) ( f h , v h u ˙ h ) V h , v h V h , t [ 0 , T ] ,
a β h ( φ h , ψ h ) a e ( u h , ψ h ) = ( q h , ψ ) W h , ψ h W h ,
θ ˙ h ( t ) + K h ˜ θ h ( t ) = M h ˜ u ˙ h ( t ) + P ( t ) , in E h ,
u h ( 0 ) = u 0 h , φ h ( 0 ) = φ 0 h , θ h ( 0 ) = θ 0 h .

5. Existence and Uniqueness Result

In this section, we use the variational formulation to prove the existence and uniqueness of a weak solution to Problem 1. For that, we make use of some auxiliary problems and lemmas.
Problem 3
(Auxiliary Problem). Determine a displacement field u h : [ 0 , T ] V h , such that
a h ( u h ( t ) , v h u ˙ h ( t ) ) + b h ( u ˙ h ( t ) , v h u ˙ h ( t ) ) +   j h ( u ˙ h ( t ) , v h ) j h ( u ˙ h ( t ) , u ˙ h ( t ) ) ( F ( t ) , v h u ˙ h ( t ) ) V h , v h V h , t [ 0 , T ] ,
u h ( 0 ) = u 0 h .
In order to study Auxiliary Problem 3, we assume that
a h : V h × V h R is a bilinear form and there exists M h > 0 , such that | a h ( u h , v h ) |   M h u h V h v h V h , u h , v h V h ;
b h : V h × V h is a bilinear symmetric form verifying : ( a ) There exists M h > 0 , such that | b h ( u h , v h ) | M h u h V h v h V h , u h , v h X h , ( b ) there exists m h > 0 , such that b h ( v h , v h ) m h v h V h 2 , v h V h ;
j h : V h × V h R satisfies : ( a ) For all η V h , j h ( η , · ) is convex and l . s . c . on V h , ( b ) there exists α h 0 , such that for all η 1 , η 2 , v 1 h , v 2 h V h we have | j h ( η 1 , v 2 h ) j h ( η 1 , v 1 h ) + j h ( η 2 , v 1 h ) j h ( η 2 h , v 2 h ) |   α h η 1 η 2 V h v 1 h v 2 h V h ;
u 0 h V h ;
F W 1 , 2 ( 0 , T , V h ) .
From this point onward, we denote a generic constant as c h > 0 .
Lemma 1.
Assume that (82)–(86) holds. If m h > α h , then there exists a unique solution u h W 1 , 2 ( 0 , T , V h ) to Auxiliary Problem 3.
Proof. 
By using (80), for any t 1 and t 2 [ 0 , T ] we find that
a h ( u h ( t 1 ) , v h u ˙ h ( t 1 ) ) +   b h ( u ˙ h ( t 1 ) , v h u ˙ h ( t 1 ) ) + j h ( u ˙ h ( t 1 ) , v h ) j h ( u ˙ h ( t 1 ) , u ˙ h ( t 1 ) ) ( F ( t 1 ) , v h u ˙ h ( t 1 ) ) V h , v h V h ,
a h ( u h ( t 2 ) , v h u ˙ h ( t 2 ) ) +   b h ( u ˙ h ( t 2 ) , v h u ˙ h ( t 2 ) ) + j h ( u ˙ ( t 2 ) , v ) j h ( u ˙ h ( t 2 ) , u ˙ h ( t 2 ) ) ( F ( t 2 ) , v h u ˙ h ( t 2 ) ) V h , v h V h .
We take v h = u ˙ h ( t 2 ) in (87) and v h = u ˙ h ( t 1 ) in (88) and, adding the two inequalities, we obtain
b h ( u ˙ h ( t 1 ) u ˙ h ( t 2 ) , u ˙ h ( t 1 ) u ˙ h ( t 2 ) ) + j ( u ˙ h ( t 1 ) , u ˙ h ( t 2 ) ) j h ( u ˙ h ( t 1 ) , u ˙ h ( t 2 ) ) + j h ( u ˙ h ( t 2 ) , u ˙ h ( t 1 ) ) j h ( u ˙ h ( t 2 ) , u ˙ h ( t 2 ) ) a h ( u h ( t 1 ) u h ( t 2 ) , u ˙ h ( t 2 ) u ˙ h ( t 1 ) ) + ( F ( t 1 ) F ( t 2 ) , u ˙ h ( t 2 ) u ˙ h ( t 1 ) ) .
We now use assumptions (82)–(84) to obtain
u ˙ h ( t 1 ) u ˙ h ( t 2 ) V h   c h ( u h ( t 1 ) u h ( t 2 ) V h + F ( t 1 ) F ( t 2 ) V h ) ,
where c h = max M h m h α h , 1 m h α h . This inequality, combined with the regularity u h C 1 ( 0 , t , V h ) , shows that u ˙ h : [ 0 , T ] V h is an absolutely continuous function. Furthermore, it should be noted that
u ¨ h ( t ) V h c h ( u ˙ h ( t ) V h + F ( t ) V h ) a . e . , t [ 0 , T ] .
We conclude that u h W 2 , 2 ( 0 , T , V h ) . □
Using the solution u h asserted by Lemma 1, we proceed by considering a second auxiliary problem.
Problem 4
(Auxiliary Problem). Determine a function θ η h : [ 0 , T ] E h , such that
θ ˙ η h ( t ) + K h ˜ θ η h ( t ) = P ( t ) + M h ˜ u ˙ η h ( t ) , θ η h ( 0 ) = θ 0 h , t [ 0 , T ] .
Lemma 2.
For all η C ( 0 , T , V h ) , there exists a unique solution to Auxiliary Problem 4 with
θ η h L 2 0 , T , E h C 0 , T , L 2 Ω h W 1 , 2 ( 0 , T , ( E h ) ) .
Moreover, there exists c h > 0 , such that, for all η 1 , η 2 C ( 0 , T , V h ) , we have
θ η 1 h ( t ) θ η 2 h ( t ) E h c h 0 t u ˙ η 1 h ( s ) u ˙ η 2 h ( s ) V h 2 d s , t [ 0 , T ] .
Proof. 
The existence and uniqueness result to (91) follows by considering the Gelfand evolution (see, e.g., [30]):
E h F h ( F h ) ( E h ) ,
and verifying that the operator K h ˜ : E h ( E h ) is strongly monotonic and linearly continuous. As inferred from the expression of the operator M h ˜ , we have
u ˙ η h ( t ) W 1 , 2 ( 0 , T ; V h ) M h ˜ u ˙ η h ( t ) W 1 , 2 ( 0 , T , F h ) ,
and
P ( t ) W 1 , 2 ( 0 , T , E h ) M h ˜ u ˙ η h ( t ) + P ( t ) W 1 , 2 ( 0 , T , ( E h ) ) .
For η 1 , η 2 C ( 0 , T , V h ) , we have
θ ˙ η 1 h ( t ) θ η 2 ˙ h ( t ) , θ η 1 h ( t ) θ η 2 h ( t ) ( E h ) × E h + K h ˜ θ η 1 h ( t ) K h ˜ θ η 2 h ( t ) , θ η 1 h ( t ) θ η 2 h ( t ) E h = M h ˜ u ˙ η 1 h ( t ) M h ˜ u ˙ η 2 h ( t ) , θ η 1 h ( t ) θ η 2 h ( t ) L 2 ( Ω h ) , t [ 0 , T ] ,
and, by integrating the last property over ( 0 , t ) and utilizing the Lipschitz continuity of M h ˜ : V h ( E h ) , as well as the strong monotonicity of K h ˜ , we deduce that (93) holds for t [ 0 , T ] . □
We are now ready to prove our main result: the existence of a unique solution to our Problem 1.
Theorem 1.
Assume that (51)–(68) are satisfied. Then, there exists Z 0 , which depends on Ω h , Γ 1 h , Γ 2 h , and Γ 3 h , such that if L r < Z 0 then there exists a unique solution ( u h , φ h , θ h ) to (VP), satisfying
u h W 2 , 2 ( 0 , T , V h ) ,
θ h L 2 0 , T , E h C 0 , T , L 2 Ω h W 1 , 2 ( 0 , T , ( E h ) ) ,
φ h W 1 , 2 ( 0 , T , W h ) .
Proof. 
From (77) we ascertain that
( β h φ h , ψ h ) W h ( e h u h , ψ h ) W h = ( q h , ψ h ) W h ,
while the use of (98) gives
β h φ h ( t ) = e h u h ( t ) + q h .
Hence, we deduce that
φ h ( t ) = e h β h u h ( t ) + q h β h .
From (78) and (79), we obtain the solution to
θ h ( t ) = 0 t M h ˜ e s t K h ˜ d x u ˙ h ( s ) d s + 0 t e s t K h ˜ d x P ( s ) d s + θ 0 h e 0 t K h ˜ d s .
Taking (99) and (100) and substituting them into (76), we obtain
a θ h ( u ˙ h ( t ) , v h u ˙ h ( t ) ) + a μ h ( u h ( t ) , v h u ˙ h ( t ) ) + a e h e h β h u h ( t ) + q h β h , v h u h ˙ ( t ) + a M h ( 0 t M h ˜ e s t K h ˜ d x u ˙ h ( s ) d s + 0 t e s t K h ˜ d x P ( s ) d s + θ 0 h e 0 t K h ˜ d s , v h u ˙ h ( t ) ) + j h ( u ˙ h ( t ) , v h ) j h ( u ˙ h ( t ) , u ˙ h ( t ) ) ( f h ( t ) , v h u ˙ h ) V h , v h V h , t [ 0 , T ] , u h ( 0 ) = u 0 h .
Next, we define the bilinear forms a h : V h × V h R and b h : V h × V h R as follows:
a h ( u h ( t ) , v h u ˙ h ( t ) ) = a μ h ( u h ( t ) , v h u ˙ h ( t ) ) + a e h e h β h u h ( t ) , v h u ˙ h ( t ) , +   a M h 0 t M h ˜ e s t K h ˜ d x u ˙ h ( s ) d s , v h u ˙ h ( t ) ,
b h ( u ˙ h ( t ) , v h u ˙ h ( t ) ) = a θ h ( u ˙ h ( t ) , v h u ˙ h ( t ) ) .
We also consider the function F : [ 0 , T ] V h , defined by
( F ( t ) , v h u ˙ h ) V h = ( f h ( t ) , v h u ˙ h ) V h a e h q h β h , v h u ˙ h ( t ) a M h 0 t e s t K h ˜ d x P ( s ) d s + θ 0 h e 0 t K h ˜ d s , v h u ˙ h ( t ) .
From (101)–(103), we obtain the following results:
  • the bilinear form a h and the initial data u 0 h satisfy conditions (82) and (85);
  • f h W 1 , 2 ( 0 , T , V h ) , q h W 1 , 2 ( 0 , T , W h ) , and P W 1 , 2 ( 0 , T , E h ) , combined with the definition of F ( · ) in (103), satisfy (86);
  • for all η V h , the functional j h ( η , · ) : V h R is a continuous seminorm on V h , thus satisfying the condition (84)-(a) (recall that j h also fulfills inequality (63), demonstrating that condition (84)-(b) holds with α h = c V L r );
  • from (74), it is evident that the bilinear form b h satisfies condition (84) with m h = α h * .
Now, we choose Z 0 = α h * c V 2 , which depends on Ω h , Γ 1 h , Γ 2 h , Γ 3 h , and α h . Then, if L r < Z 0 , we have m h > α h . Therefore, (95) and (96) are a direct consequence of Lemmas 1 and 2. For the proof of (97), let φ h : [ 0 , T ] W h be the electrical-potential field defined by (99). We have
u h W 2 , 2 ( 0 , T , V h ) q h W 1 , 2 ( 0 , T , W h ) φ h W 1 , 2 ( 0 , T , W h ) ,
which completes the proof. □

6. Conclusions

In this paper, we theoretically investigated a quasi-static-antiplane contact problem with a slip-rate-dependent friction law involving a thermo-electro-visco-elastic body. The slip-rate-dependent friction law was interesting to investigate because the friction level depends on the rate at which the slip changes over time. This implies that variations in the slip rate impact the friction level, with the friction force increasing as the slip-rate speed increases, and decreasing as the slip-rate speed decreases. This is more reasonable because it considers both the quasi-static equilibrium conditions and the localized dynamic interactions that may occur within the contact interface.
We used Green’s formula to derive the variational formulation of our problem. Following that, we demonstrated the existence and uniqueness of the weak solution, using various techniques, including the time-dependent variational equation, the variational-evolution equation, the differential equation, and the decoupling of unknowns. The advantages of the developed model include both the quasi-static equilibrium conditions and the localized dynamic interactions that may occur within the contact interface. This approach provides a more comprehensive and realistic representation of the material response under varying thermal and loading conditions. One of the drawbacks we faced was the time consumption and the difficulty of providing practical examples without resorting to numerical analysis and simulation.
We have numerous prospects ahead. In particular, we intend to conduct work on numerical analysis and simulation based on the results obtained for the studied problem.

Author Contributions

Conceptualization, B.F. and M.D.; validation, B.F., M.D. and D.F.M.T.; formal analysis, B.F., M.D. and D.F.M.T.; investigation, B.F., M.D. and D.F.M.T.; writing—original draft preparation, B.F., M.D. and D.F.M.T.; writing—review and editing, B.F., M.D. and D.F.M.T.; supervision, M.D. and D.F.M.T.; project administration, M.D. and D.F.M.T.; funding acquisition, B.F., M.D. and D.F.M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundação para a Ciência e a Tecnologia grant number UIDB/04106/2020 (https://doi.org/10.54499/UIDB/04106/2020, accessed on 25 January 2024).

Data Availability Statement

Data are contained within the article.

Acknowledgments

Fadlia is grateful for the financial support of University of Constantine 1, Algeria, and for a one-month visit to the R&D Unit CIDMA, Department of Mathematics, University of Aveiro. The hospitality of the host institution is here gratefully acknowledged. The authors are grateful to four anonymous referees for several constructive comments and remarks.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Dosaev, M.; Samsonov, V.; Bekmemetev, V. Comparison between 2D and 3D simulation of contact of two deformable axisymmetric bodies. Int. J. Nonlinear Sci. Numer. Simul. 2020, 21, 123–133. [Google Scholar] [CrossRef]
  2. Fernandez, F.; Puso, M.A.; Solberg, J.; Tortorelli, D.A. Topology optimization of multiple deformable bodies in contact with large deformations. Comput. Methods Appl. Mech. Eng. 2020, 371, 113288. [Google Scholar] [CrossRef]
  3. Biber, S.W.; Champneys, A.R.; Szalai, R. Analysis of point-contact models of the bounce of a hard spinning ball on a compliant frictional surface. IMA J. Appl. Math. 2023, 88, 498–523. [Google Scholar] [CrossRef]
  4. Shaat, M.; Gao, X.-L.; Li, K.; Littlefield, A.G. New analytical model for thermomechanical responses of multi-layered structures with imperfect interfaces. Acta Mech. 2023, 234, 5779–5818. [Google Scholar] [CrossRef]
  5. Duvaut, G.; Lions, J.-L. Les Inéquations en Mécanique et en Physique; Dunod: Paris, France, 1972. [Google Scholar]
  6. Lions, J.-L. Quelques Méthodes de Résolution des Problèmes aux Limites non Linéaires; Gauthier-Villars: Paris, France, 1969. [Google Scholar]
  7. Rakici, S.; Kim, J. A stabilized non-ordinary peridynamic model for linear piezoelectricity. Appl. Math. Model. 2024, 125, 514–538. [Google Scholar] [CrossRef]
  8. Xu, J.; Xu, X. On uniqueness of recovering coefficients from localized Dirichlet-to-Neumann map for piecewise homogeneous piezoelectricity. SIAM J. Math. Anal. 2023, 55, 571–602. [Google Scholar] [CrossRef]
  9. Benaissa, H.; Essoufi, E.-H.; Fakhar, R. Existence results for unilateral contact problem with friction of thermo-electro-elasticity. Appl. Math. Mech. 2015, 36, 911–926. [Google Scholar] [CrossRef]
  10. Benaissa, H.; Essoufi, E.-H.; Fakhar, R. Analysis of a Signorini problem with nonlocal friction in thermo-piezoelectricity. Glas. Mat. Ser. III 2016, 51, 391–411. [Google Scholar] [CrossRef]
  11. Ammar, T.H. Quasistatic contact problem between thermo-electroelastic bodies with long-term memory and adhesion. Malaya J. Mat. 2016, 4, 211–223. [Google Scholar] [CrossRef] [PubMed]
  12. Sládek, J.; Sládek, V.; Stanak, P. Analysis of thermo-piezoelectricity problems by meshless method. Acta Mech. Slovaca 2010, 14, 16–27. [Google Scholar] [CrossRef]
  13. Tiersten, H.F. On the nonlinear equations of thermoelectroelasticity. Internat. J. Engrg. Sci. 1971, 9, 587–604. [Google Scholar] [CrossRef]
  14. Ahmed, A.A.; Boutechebak, S. Analysis of a dynamic thermo-elastic-viscoplastic contact problem. Electron. J. Qual. Theory Differ. Equ. 2013, 2013, 17. [Google Scholar] [CrossRef]
  15. Merouani, A.; Messelmi, F. Dynamic evolution of damage in elastic-thermo-viscoplastic materials. Electron. Differ. Equ. 2010, 2010, 1–15. [Google Scholar]
  16. Sofonea, M. Functional Methods in Thermo-Elasto-Visco-Plasticity. Ph.D. Thesis, University of Bucharest, Bucharest, Romania, 1988. [Google Scholar]
  17. Adly, S.; Chau, O. On some dynamic thermal non clamped contact problems. Math. Program. 2013, 139, 5–26. [Google Scholar] [CrossRef]
  18. Andrews, K.T.; Shillor, M.; Wright, S.; Klarbring, A. A dynamic thermoviscoelastic contact problem with friction and wear. Internat. J. Engrg. Sci. 1997, 35, 1291–1309. [Google Scholar] [CrossRef]
  19. Dalah, M. Analysis of electro-viscoelastic antiplane contact problem with total slip rate dependent friction. Electron. J. Differ. Equ. 2009, 2009, 15. [Google Scholar]
  20. Dalah, M.; Sofonea, M. Antiplane frictional contact of electro-viscoelastic cylinders. Electron. J. Differ. Equ. 2007, 2007, 14. [Google Scholar]
  21. Matei, A.; Motreanu, V.V.; Sofonea, M. A quasistatic antiplane contact problem with slip dependent friction. Adv. Nonlinear Var. Inequal. 2001, 4, 1–21. [Google Scholar]
  22. Mindlin, R.D. Polarization gradient in elastic dielectrics. Int. J. Solids Structures 1968, 4, 637–642. [Google Scholar] [CrossRef]
  23. Rachid, C.; Nemira, L. Analysis of a thermo-viscoelastic antiplane contact problem with long-term memory. An. Ştiinţ. Univ. Al. I. Cuza Iaşi. Mat. (N.S.) 2021, 67, 97–111. [Google Scholar] [CrossRef]
  24. Sofonea, M.; Dalah, M.; Ayadi, A. Analysis of an antiplane electro-elastic contact problem. Adv. Math. Sci. Appl. 2007, 17, 385–400. [Google Scholar]
  25. Sofonea, M.; Matei, A. Variational Inequalities with Applications: A Study of Antiplane Frictional Contact Problems; Springer: New York, NY, USA, 2009. [Google Scholar]
  26. Horgan, C.O. Anti-plane shear deformations in linear and nonlinear solid mechanics. SIAM Rev. 1995, 37, 53–81. [Google Scholar] [CrossRef]
  27. Horgan, C.O.; Miller, K.L. Antiplane shear deformations for homogeneous and inhomogeneous anisotropic linearly elastic solids. Trans. ASME J. Appl. Mech. 1994, 61, 23–29. [Google Scholar] [CrossRef]
  28. Han, W.; Sofonea, M. Quasistatic Contact Problems in Viscoelasticity and Viscoplasticity; American Mathematical Society: Providence, RI, USA; RIInternational Press: Somerville, MA, USA, 2002. [Google Scholar]
  29. Shillor, M.; Sofonea, M.; Telega, J.J. Models and Analysis of Quasistatic Contact: Variational Methods; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  30. Zeidler, E. Nonlinear Functional Analysis and Its Applications; II/A; Springer: New York, NY, USA, 1990. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fadlia, B.; Dalah, M.; Torres, D.F.M. Existence and Uniqueness of Weak Solutions to Frictionless-Antiplane Contact Problems. Mathematics 2024, 12, 434. https://doi.org/10.3390/math12030434

AMA Style

Fadlia B, Dalah M, Torres DFM. Existence and Uniqueness of Weak Solutions to Frictionless-Antiplane Contact Problems. Mathematics. 2024; 12(3):434. https://doi.org/10.3390/math12030434

Chicago/Turabian Style

Fadlia, Besma, Mohamed Dalah, and Delfim F. M. Torres. 2024. "Existence and Uniqueness of Weak Solutions to Frictionless-Antiplane Contact Problems" Mathematics 12, no. 3: 434. https://doi.org/10.3390/math12030434

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop