Next Article in Journal
Exact Traveling Wave Solutions of the Local Fractional Bidirectional Propagation System Equations
Next Article in Special Issue
Unique Solvability of the Initial-Value Problem for Fractional Functional Differential Equations—Pantograph-Type Model
Previous Article in Journal
A Novel Analytical LRPSM for Solving Nonlinear Systems of FPDEs
Previous Article in Special Issue
Qualitative Properties of Positive Solutions of a Kind for Fractional Pantograph Problems using Technique Fixed Point Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Three-Field Variational Formulation for a Frictional Contact Problem with Prescribed Normal Stress

by
Andaluzia Cristina Matei
Department of Mathematics, University of Craiova, A.I. Cuza 13, 200585 Craiova, Romania
Submission received: 26 September 2022 / Revised: 19 October 2022 / Accepted: 23 October 2022 / Published: 4 November 2022
(This article belongs to the Special Issue Initial and Boundary Value Problems for Differential Equations)

Abstract

:
In the present work, we address a nonlinear boundary value problem that models frictional contact with prescribed normal stress between a deformable body and a foundation. The body is nonlinearly elastic, the constitutive law being a subdifferential inclusion. We deliver a three-field variational formulation by means of a new variational approach governed by the theory of bipotentials combined with a Lagrange-multipliers technique. In this new approach, the unknown of the mechanical model is a triple consisting of the displacement field, a Lagrange multiplier related to the friction force and the Cauchy stress tensor. We obtain existence, uniqueness, boundedness and convergence results.

1. Introduction

Everywhere around us we can see deformable bodies in interaction. However, even though very common, the contact phenomenon is not a trivial one; in addition, in engineering, handling the interactions between deformable bodies and obstacles is very important and requires advanced applied mathematics. The contact phenomenon can be mathematically modeled by means of boundary value problems governed by partial differential equations. Actually, the topic is very complex, involving continuum mechanics, differential equations, function spaces, calculus of variations, nonlinear analysis, control theory and numerical analysis. The importance and the abundance of the applications of contact problems in the real world has motivated a large number of scientists to investigate this kind of model. It is worth underlining that, due to their complexity, contact models do not have classical solutions. Thus, the variational methods play a crucial role in the qualitative and quantitative analysis.
In the present paper, we focus on the well-posedness and approximation results addressing a stationary frictional contact problem with prescribed normal stress, for materials governed by a multi-valued elastic operator. Using the bipotential theory, we deliver a variational formulation of the mechanical model in a form of a variational system consisting of three inequalities. Placing us in an appropriate functional setting governed by Lebesgue and Sobolev spaces for vector functions including fractional spaces on the boundary, we apply the saddle point theory and a minimization technique in order to prove the existence of at least one solution. We also pay attention to the uniqueness of the solution. Firstly, we draw attention to a partial uniqueness result related to the uniqueness in the first component. Subsequently, we discuss a global uniqueness result, not for the original problem, but for a perturbed version of it. After we investigate the boundedness of the solution of the perturbed problem, we prove a convergence result allowing an approximation of a weak solution of the contact model under consideration. The present study can be seen as a continuation of [1]; there, a two-field variational formulation for the same model was delivered; with the well-posedness of the model being studied under a more restrictive hypothesis for the prescribed normal stress. In [1], the weak solution is a pair consisting of the displacement vector and the Cauchy stress tensor. In the present study, the weak solution is a triple by considering a Lagrange multiplier related to the friction force as a component of the weak solution, in addition to the displacement vector and the Cauchy stress tensor. From the mathematical point of view, the new study is more complex. Besides the theory of bipotentials and the minimization techniques, the present study requires a saddle point technique, fractional Lebesgue and Sobolev spaces, weak topologies and a convergence of the Mosco type. It is worth emphasizing that the variational approach we propose leads to a new class of variational problems governed by Lagrange multipliers. From the mechanical point of view, the new variational approach we propose is important because it makes possible an estimation of the friction force, even a numerical computation, after passing from the qualitative study to the quantitative analysis in a future investigation.
Let us specify some helpful references for background knowledge: for the mechanics of deformable solids/contact mechanics see, e.g., [2,3,4,5,6,7,8,9]; for bipotential theory, we refer to, e.g., [10,11,12,13,14,15]; for the saddle point theory see, e.g., [16,17,18,19,20]; and for the functional spaces the reader can consult, e.g., [21,22,23,24,25,26,27]. For recent results related to the topic of the present work, we refer, for instance, to [28,29,30]. However, in order to increase the clarity of the exposure, we specify herein some mathematical tools that will play a crucial role.
Let ( X , ( · , · ) X , · X ) be a Hilbert space.
Definition 1.
A bipotential is a function B : X × X ( , ] with the following three properties:
  • ( i ) B is convex and lower semicontinuous in each argument;
  • ( i i ) For each x , y X , we have B ( x , y ) ( x , y ) X ;
  • ( i i i ) For each x , y X , y B ( · , y ) ( x ) x B ( x , · ) ( y ) B ( x , y ) = ( x , y ) X .
Recall that the Fenchel conjugate of a functional ϕ : X ( , ] is the functional
ϕ * : X ( , ] , ϕ * ( x * ) = sup x X x * , x ϕ ( x ) .
Always, the Fenchel conjugate is a convex and lower semicontinuous functional, see, e.g., [22].
Theorem 1.
Let ϕ : X ( , ] be a proper, convex, lower semicontinuous functional. Then:
  • ( i ) For each x , y X , we have ϕ ( x ) + ϕ * ( y ) ( x , y ) X ;
  • ( i i ) For each x , y X , y ϕ ( x ) x ϕ * ( y ) ϕ ( x ) + ϕ * ( y ) = ( x , y ) X .
Notice that if ϕ is a proper, convex, lower semicontinuous functional, then its Fenchel conjugate has all these properties too, see, e.g., [22]. For some details and additional elements in the convex analysis, we refer to, e.g., [18,31,32,33,34].
In addition, we shall need the following theorem.
Theorem 2
(See, e.g., [35]). Let ( X , · X ) be a reflexive Banach space and let K X be a nonempty, convex, closed, unbounded subset of X. Suppose φ : K R is coercive, convex and lower semicontinuous. Then, φ is bounded from below on K and attains its infimum in K. If φ is strictly convex, then φ has a unique minimizer.
We recall that φ : K R is coercive if, for all u K , we have φ ( u )   a s u X .
Theorem 3.
Let ( X , ( · , · ) X , · X ) , ( Y , ( · , · ) Y , · Y ) be two Hilbert spaces and let A X , B Y be nonempty, closed, convex subsets. Assume that a bifunctional L : A × B R satisfies the following conditions
v L ( v , μ ) i s   c o n v e x   a n d   l o w e r   s e m i c o n t i n u o u s f o r   a l l μ B ,
μ L ( v , μ ) i s   c o n c a v e   a n d   u p p e r   s e m i c o n t i n u o u s f o r   a l l v A .
Moreover, we assume that
A is s   bounded or lim v X , v A L ( v , μ 0 ) =   for   some   μ 0 B
and
B is   bounded   or lim μ Y , μ B inf v A L ( v , μ ) = .
Then:
(a)
The bifunctional L ( · , · ) has at least one saddle point;
(b)
The set A 0 × B 0 of the saddle points of L is convex, where A 0 A and B 0 B ;
(c)
If v L ( v , μ ) i s   s t r i c t l y   c o n v e x f o r   a l l μ B , then A 0 contains at most one point;
(d)
If μ L ( v , μ ) i s   s t r i c t l y   c o n c a v e f o r   a l l v A , then B 0 contains at most one point.
For a proof of ( a ) , see [18] (p. 176). For a proof of ( b ) , ( c ) , ( d ) , see [18] (p. 169).
The rest of the paper has the following structure. In Section 2 we provide the functional setting we use. In Section 3 we describe the mechanical model and deliver its three-field variational formulation. In Section 4 we obtain existence, uniqueness, boundedness and convergence results. The last section provides some conclusions and final comments.

2. Functional Setting

Let Ω R N ( N > 1 ) be a bounded domain with a regular enough boundary denoted by Γ . We use standard notation for L p ( Ω ) , W 1 , p ( Ω ) , L p ( Γ ) , H 1 ( Ω ) , see, e.g., [21,22,23,24,25,26,27].
Let 2 p < . Recall that, according to the trace theorem—see, e.g., [6] (p. 34)—there exists a unique linear continuous operator γ : W 1 , p ( Ω ) L p ( Γ ) such that
(a)
γ u = u | Γ if u C 1 ( Ω ¯ ) ;
(b)
γ u L p ( Γ ) c t r u W 1 , p ( Ω ) with c t r = c t r ( p , Ω ) > 0 ;
(c)
If 1 < p < , then γ ( W 1 , p ( Ω ) ) = W 1 1 / p , p ( Γ ) ;
(d)
If 1 < p < N , then γ : W 1 , p ( Ω ) L r ( Γ ) is compact for any r such that 1 r < N p p N p .
The function γ u is called the trace of the scalar function u on Γ and the operator γ : W 1 , p ( Ω ) L p ( Γ ) is called the trace operator. The trace operator is neither an injection, nor a surjection from W 1 , p ( Ω ) to L p ( Γ ) . However, according to, e.g., Theorem 1.5.1.2 in [23], there exists a unique linear continuous and surjective operator from W 1 , p ( Ω ) to γ ( W 1 , p ( Ω ) ) = W 1 1 / p , p ( Γ ) . Therefore, there are some d ¯ = d ¯ ( p , Ω ) , such that
γ u W 1 1 p , p ( Γ ) d ¯ u W 1 , p ( Ω ) for   all u W 1 , p ( Ω ) .
Furthermore, there exists a map called the right inverse of γ , denoted here by ,
: W 1 1 / p , p ( Γ ) W 1 , p ( Ω ) , γ ( ( ζ ) ) = ζ for   all ζ W 1 1 / p , p ( Γ ) ;
see, e.g., [36]. Notice that the right inverse operator is a linear and continuous map.
The space
W 1 1 / p , p ( Γ ) = { u L p ( Γ ) | v W 1 , p ( Ω )   such   that u = γ v   a . e . on   Γ }
is a normed space endowed with the following norm
u W 1 1 / p , p ( Γ ) = u L p ( Γ ) p + Γ Γ | u ( x ) u ( y ) | p x y N + p 2 d s ( x ) d s ( y ) 1 / p ,
see, e.g., [26], (p. 43). Actually, the space ( W 1 1 / p , p ( Γ ) , · W 1 1 / p , p ( Γ ) ) is a Banach space, see, e.g., [24], (p. 332).
If p = 2 , γ ( H 1 ( Ω ) ) = H 1 / 2 ( Γ ) . The space H 1 / 2 ( Γ ) is a Hilbert space endowed with the inner product
( v , w ) H 1 / 2 ( Γ ) = ( v , w ) L 2 ( Γ ) + Γ Γ ( v ( x ) v ( y ) ) ( w ( x ) w ( y ) ) x y N d Γ ( x ) d Γ ( y )
and the corresponding Sobolev–Slobodeckij norm
v H 1 / 2 ( Γ ) = v L 2 ( Γ ) 2 + Γ Γ ( v ( x ) v ( y ) ) 2 x y N d Γ ( x ) d Γ ( y ) 1 / 2 .
Notice that
u L 2 ( Γ ) u H 1 / 2 ( Γ )   for   all   u H 1 / 2 ( Γ ) ,
so H 1 / 2 ( Γ ) is continuously embedded in L 2 ( Γ ) . The space H 1 / 2 ( Γ ) is called the image of H 1 ( Ω ) by the trace operator γ .
To proceed, we consider N = 3 . The fields in R 3 will be typeset in boldface. The inner product and the Euclidean norm on R 3 will be denoted by · and · , respectively.
Let us introduce the vector spaces
L 2 ( Ω ; R 3 ) = { w = ( w i ) | w i L 2 ( Ω ) , 1 i 3 } ; H 1 ( Ω ; R 3 ) = { w = ( w i ) | w i H 1 ( Ω ) , 1 i 3 } ; H 1 / 2 ( Γ ; R 3 ) = γ ( H 1 ( Ω ; R 3 ) ) = { w = ( w i ) | w i H 1 / 2 ( Γ ) , 1 i 3 } .
The spaces L 2 ( Ω ; R 3 ) and H 1 ( Ω ; R 3 ) are Hilbert spaces endowed with the inner products
( u , v ) L 2 ( Ω ; R 3 ) = i = 1 3 Ω u i v i d x , ( u , v ) H 1 ( Ω ; R 3 ) = i = 1 3 ( u i , v i ) H 1 ( Ω ) ,
and the associated norms · L 2 ( Ω ; R 3 ) = ( · , · ) L 2 ( Ω ; R 3 ) and · H 1 ( Ω ; R 3 ) = ( · , · ) H 1 ( Ω ; R 3 ) , respectively. The space H 1 / 2 ( Γ ; R 3 ) is also a Hilbert space endowed with the inner product
( v , w ) H 1 / 2 ( Γ ; R 3 ) = i = 1 3 ( v i , w i ) H 1 / 2 ( Γ ) = i = 1 3 ( v i , w i ) L 2 ( Γ ) + i = 1 3 Γ Γ ( v i ( x ) v i ( y ) ) ( w i ( x ) w i ( y ) ) x y 3 d Γ ( x ) d Γ ( y )
and the corresponding Sobolev–Slobodeckij norm
w H 1 / 2 ( Γ ; R 3 ) = w L 2 ( Γ ; R 3 ) 2 + Γ Γ w ( x ) w ( y ) 2 x y 3 d Γ ( x ) d Γ ( y ) 1 / 2 .
We easily observe that
u L 2 ( Γ ; R 3 ) u H 1 / 2 ( Γ ; R 3 )   for   all   u H 1 / 2 ( Γ ; R 3 ) ,
so H 1 / 2 ( Γ ; R 3 ) is continuously embedded in L 2 ( Γ ; R 3 ) .
The trace operator for vector functions γ : H 1 ( Ω ; R 3 ) L 2 ( Γ ; R 3 ) is a linear, continuous and compact operator, but it is neither an injection nor a surjection; see, e.g., [37].
γ u L 2 ( Γ ; R 3 ) c ¯ t r u H 1 ( Ω ; R 3 ) for   all     u H 1 ( Ω ; R 3 ) ( c ¯ t r > 0 ) .
Notice that
γ : H 1 ( Ω ; R 3 ) H 1 / 2 ( Γ ; R 3 )
is a linear, continuous operator but it is not a compact operator. Let us point out that there are some c ˜ > 0 such that
γ u H 1 / 2 ( Γ ; R 3 ) c ˜ u H 1 ( Ω ; R 3 ) for   all   u H 1 ( Ω ; R 3 ) .
Furthermore, there exists a linear, continuous operator : H 1 / 2 ( Γ ; R 3 ) H 1 ( Ω ; R 3 ) such that
γ ( ( ξ ) ) = ξ for or all ξ H 1 / 2 ( Γ ; R 3 ) .
The operator is called the right inverse of the trace operator γ .
Let S 3 be the space of second-order symmetric tensors on R 3 . Every field in S 3 is typeset in boldface. By : and · S 3 we denote the inner product and the Euclidean norm on S 3 .
We introduce now two tensor Lebesgue spaces, as follows.
L 2 ( Ω ; S 3 ) = { μ = ( μ i j ) : μ i j L 2 ( Ω ) for   all   i , j { 1 , 2 , 3 } } ; L s 2 ( Ω ; S 3 ) = { μ = ( μ i j ) : μ i j L 2 ( Ω ) , μ i j = μ j i for   all   i , j { 1 , 2 , 3 } } .
The space L 2 ( Ω ; S 3 ) is a Hilbert space endowed with the inner product
( μ , τ ) L 2 ( Ω ; S 3 ) = Ω i , j = 1 3 μ i j ( x ) τ i j ( x ) d x ;
the space L s 2 ( Ω ; S 3 ) is a Hilbert space endowed with the inner product
( μ , τ ) L s 2 ( Ω ; S 3 ) = ( μ , τ ) L 2 ( Ω ; S 3 ) .
Next, we introduce the following space,
H 1 = { u L 2 ( Ω ; R 3 ) | ε ( u ) L s 2 ( Ω ; S 3 ) } ,
where
ε : H 1 ( Ω ; R 3 ) L 2 ( Ω ; S 3 )
is the linear operator
ε ( u ) = ( ε i j ( u ) ) , ε i j ( u ) = 1 2 ( u i , j + u j , i ) ;
the index that follows a comma indicates a weak partial derivative with respect to the corresponding component of the independent variable.
The space H 1 is a real Hilbert space endowed with the inner product
( u , v ) H 1 = ( u , v ) L 2 ( Ω ; R 3 ) + ( ε ( u ) , ε ( v ) ) L s 2 ( Ω ; S 3 ) .
The associated norm on the space H 1 is denoted by · H 1 . According to, e.g., [38], H 1 = H 1 ( Ω ; R 3 ) algebraically and the norms · H 1 and · H 1 ( Ω ; R 3 ) are equivalent.
Notice that ε is a linear and continuous operator from H 1 ( Ω ; R 3 ) to L s 2 ( Ω ; S 3 ) .
Let Γ 1 be a measurable part of Γ with positive surface measure.
We consider the space
V = { v H 1 | γ v = 0   almost   everywhere   on Γ 1 } .
This is a closed subspace of H 1 ( Ω ; R 3 ) , so ( V , ( · , · ) H 1 ( Ω ; R 3 ) , · H 1 ( Ω ; R 3 ) ) is a Hilbert space.
Let us recall Korn’s inequality: there exists c K = c K ( Ω , Γ 1 ) > 0 such that
ε ( v ) L s 2 ( Ω ; S 3 ) c K v H 1 for   all   v V ;
see, e.g., [38,39]. Using this inequality, it can be proved that the space V is a Hilbert space endowed with the following inner product,
( · , · ) V : V × V R ; ( u , v ) V = ( ε ( u ) , ε ( v ) ) L s 2 ( Ω ; S 3 ) for   all   u , v V ,
and the corresponding norm
v V = ε ( v ) L s 2 ( Ω ; S 3 ) for   all   v V .
Notice that, since the norms · H 1 and · H 1 ( Ω ; R 3 ) are equivalent, then there exists c ˜ K = c ˜ K ( Ω , Γ 1 ) > 0 such that
v V c ˜ K v H 1 ( Ω ; R 3 ) for   all   v V .
We proceed by introducing a closed subspace of H 1 / 2 ( Γ ; R 3 ) as follows:
γ ( V ) = { v ˜ H 1 / 2 ( Γ ; R 3 ) | v ˜ = γ v   for   some   v V } .
According to [40], the space γ ( V ) is a closed subspace of H Γ = H 1 / 2 ( Γ ; R 3 ) . Thus, γ ( V ) is a Hilbert space endowed with the inner product
( · , · ) γ ( V ) : γ ( V ) × γ ( V ) R , ( ζ , ϕ ) γ ( V ) = ( ζ , ϕ ) H 1 / 2 ( Γ ; R 3 ) f o r   a l l   ζ , ϕ γ ( V ) .
We note that
γ ( ( γ v ) ) = γ v for   all   v V .
As
( ζ ) V for   all   ζ γ ( V ) ,
we can introduce an operator as follows
R : γ ( V ) V , R ( ζ ) = ( ζ ) .
The operator R is a linear and continuous operator. Hence, there are some c ¯ > 0 such that
R v ˜ V c ¯ v ˜ H 1 / 2 ( Γ ; R 3 ) for   all   v ˜ γ ( V ) .

3. The Model and Its Three-Field Variational Formulation

The physical setting is as follows: a deformable body occupies a bounded domain Ω R 3 with smooth enough boundary Γ , partitioned in three measurable parts Γ 1 , Γ 2 and Γ 3 with positive surface measures. The body is clamped on Γ 1 , body forces of density f 0 act on Ω , surface tractions of density f 2 act on Γ 2 while on Γ 3 the body is in frictional contact with a foundation.
According to this physical setting, we state the following boundary value problem.
Problem 1.
Find u : Ω ¯ R 3 and σ : Ω ¯ S 3 , such that
Div σ ( x ) + f 0 ( x ) = 0   i n   Ω ,
σ ( x ) ω ( ε ( u ) ( x ) )   i n   Ω ,
u ( x ) = 0   o n   Γ 1 ,
σ ( x ) ν ( x ) = f 2 ( x )   o n   Γ 2 ,
σ ν ( x ) = F ( x ) o n   Γ 3 ,
σ τ ( x ) k ( x ) | σ ν ( x ) | , σ τ ( x ) = k ( x ) | σ ν ( x ) | u τ ( x ) u τ ( x ) if u τ ( x ) 0 on Γ 3 .
As usual, u = ( u i ) denotes the displacement field, ε = ε ( u ) = ( ε i j ( u ) ) denotes the infinitesimal strain tensor, σ = ( σ i j ) denotes the Cauchy stress tensor, ω is a constitutive map, ν stands for the unit outward normal to Γ , the normal and the tangential components of the displacement vector on the boundary are defined by the formulas u ν = u · ν , u τ = u u ν ν and the normal and the tangential components of the Cauchy vector on the boundary are defined by σ ν = ( σ ν ) · ν , σ τ = σ ν σ ν ν .
Due to the condition
σ ν = F ,
according to the engineering literature, the frictional contact model we treat is a bilateral frictional contact problem. In this context, we have to mention that the bilateral frictional contact phenomenon can be found in many components of mechanical equipment. Thus, many real-world examples can be envisaged. The mathematical and the engineering literature contains relevant applications of bilateral frictional contact models; see, e.g., [41,42]. Referring to the behavior of the materials, for significant examples of nonlinearly elastic constitutive laws described by means of subdifferential inclusions for various constitutive maps ω , see, e.g., [38] and the references therein. For the convenience of the reader, we indicate here an example of such a constitutive map:
ω : S 3 R , ω ( τ ) = 1 2 A τ : τ + k 2 τ P K τ S 3 2 ,
where A : S 3 S 3 , A = ( A i j k l ) , A i j k l = λ 0 δ i j δ k l + μ 0 ( δ i k δ j l + δ i l δ j k ) , 1 i , j , k , l 3 , with λ 0 , μ 0 and k small enough positive material coefficients, and P K : S 3 K denotes the projection operator on the closed and convex set K S 3 which contains 0 S 3 .
In order to study Problem 1, we make the following assumptions.
Assumption 1.
ω : S 3 R is a convex and lower semicontinuous functional. In addition, there exist α , β such that 1 > β α > 0 and β ε S 3 2 ω ( ε ) α ε S 3 2 f o r   a l l   ε S 3 .
Assumption 2.
f 0 L 2 ( Ω ; R 3 ) a n d   f 2 L 2 ( Γ 2 ; R 3 ) .
Assumption 3.
F L 2 δ 1 δ ( Γ 3 ) ( 0 < δ < 1 ) and F ( x ) 0   a . e . x Γ 3 .
Assumption 4.
The coefficient of friction satisfies k L ( Γ 3 ) a n d k ( x ) 0 a . e . x Γ 3 .
Notice that the example in (14) fulfills Assumption 1.
Let ( u , σ ) be a pair of smooth-enough functions that verify Problem 1. Using Green’s formula
( σ , ε ( v ) ) L 2 ( Ω ; S 3 ) + ( Div σ , v ) L 2 ( Ω ; R 3 ) = Γ σ ( x ) ν ( x ) · γ v ( x ) d Γ for   all   v H 1 ( Ω ; R 3 ) ,
(see, e.g., [3], (p. 145)), by taking into account (8), (10) and (11) we obtain, for all v V
( σ , ε ( v ) ) L 2 ( Ω ; S 3 ) = ( f 0 , v ) L 2 ( Ω ; R 3 ) + Γ 2 f 2 ( x ) · γ v ( x ) d Γ + Γ 3 σ ( x ) ν ( x ) · γ v ( x ) d Γ ,
where V is the space defined in (3).
As
σ ( x ) ν ( x ) · γ v ( x ) = σ ν ( x ) v ν ( x ) + σ τ ( x ) · v τ ( x )   a . e . on   Γ ,
then by (15) we infer that
( σ , ε ( v ) ) L 2 ( Ω ; S 3 ) = ( f 0 , v ) L 2 ( Ω ; R 3 ) + Γ 2 f 2 ( x ) · v ( x ) d Γ + Γ 3 σ τ ( x ) · v τ ( x ) d Γ Γ 3 F ( x ) v ν ( x ) d Γ .
Herein and everywhere below, v ν ( x ) = γ v ( x ) · ν ( x ) and v τ ( x ) = γ v ( x ) v ν ( x ) ν ( x ) a.e. x Γ . By taking into consideration H3, H4, using the trace theorem for N = 3 and p = 2 and the Hölder’s inequality, it follows that F ( · ) v ν ( · ) L 1 ( Γ 3 ) for all v V .
Since
V v ( f 0 , v ) L 2 ( Ω ; R 3 ) + Γ 2 f 2 ( x ) · γ v ( x ) d Γ Γ 3 F ( x ) v ν ( x ) d Γ R
is a linear and continuous map, according to Riesz’s representation theorem, there exists a unique element f V such that
( f , v ) V = ( f 0 , v ) L 2 ( Ω ; R 3 ) + Γ 2 f 2 ( x ) · γ v ( x ) d Γ Γ 3 F ( x ) v ν ( x ) d Γ for   all   v V .
Therefore,
( σ , ε ( v ) ) L 2 ( Ω ; S 3 ) = ( f , v ) V + Γ 3 σ τ ( x ) · v τ ( x ) d Γ .
Let D be the dual of the space γ ( V ) defined in (6).
We define λ D such that
λ , v ˜ = Γ 3 σ τ ( x ) · v ˜ τ ( x ) d Γ , for   all   v ˜ γ ( V ) ,
where · , · denotes the duality pairing between D and γ ( V ) and
v ˜ τ ( x ) = v ˜ ( x ) v ˜ ν ( x ) ν ( x ) , v ˜ ν ( x ) = v ˜ ( x ) · ν ( x ) a . e . x Γ 3 .
Furthermore, we define a form c ( · , · ) as follows,
c : V × D R , c ( v , μ ) = μ , γ v , for   all   v V , μ D .
Important properties of the form c ( · , · ) are given by the following lemma.
Lemma 1.
The form c ( · , · ) is bilinear. In addition,
  • c ( · , · ) is continuous of rank M c > 0 i.e.,
    | c ( v , ζ ) | M c v V ζ D   f o r   a l l   v V , ζ D ,
  • c ( · , · ) verifies the inf-sup property:
    inf ζ D , ζ 0 D sup v V , v 0 V c ( v , ζ ) v V ζ D α c   f o r   s o m e   α c > 0 .
Proof. 
The bilinearity of c ( · , · ) is obvious keeping in mind the linearity of the trace operator γ .
Let v V and ζ D be arbitrarily given.
| c ( v , ζ ) | ζ D γ v H 1 / 2 ( Γ ; R 3 ) c ˜ ζ D v H 1 ( Ω ; R 3 ) c ˜ c ˜ K ζ D v V .
We can consider M c = c ˜ c ˜ K where c ˜ > 0 appears in (2) and c ˜ K > 0 appears in (5).
Next, we prove that c ( · , · ) verifies the inf-sup property, or equivalently, there are some α c >0 such that ζ D 1 α c sup v V , v 0 V c ( v , ζ ) v V for all ζ D , ζ 0 D .
Indeed, let ζ D , ζ 0 D .
ζ D = sup w ˜ γ ( V ) , w ˜ 0 γ ( V ) < ζ , w ˜ > w ˜ H 1 / 2 ( Γ ; R 3 ) .
As ζ , w ˜ = ζ , γ ( w ˜ ) = ζ , γ ( R w ˜ ) ) , then
ζ D = sup w ˜ γ ( V ) , w ˜ 0 γ ( V ) ζ , γ ( R w ˜ ) ) w ˜ H 1 / 2 ( Γ ; R 3 ) = sup w ˜ γ ( V ) , w ˜ 0 γ ( V ) c ( R w ˜ , ζ ) w ˜ H 1 / 2 ( Γ ; R 3 ) = sup w ˜ γ ( V ) , R w ˜ 0 V c ( R w ˜ , ζ ) w ˜ H 1 / 2 ( Γ ; R 3 ) .
Using (7), we can write
ζ D c ¯ sup w ˜ γ ( V ) , R w ˜ 0 V c ( R w ˜ , ζ ) R w ˜ V c ¯ sup v V , v 0 V c ( v , ζ ) v V .
Thus, we can take α c = 1 c ¯ .
Using the hypotheses H3, H4, as γ v L 2 δ ( Γ ; R 3 ) and v τ ( x ) γ v ( x ) a.e. x Γ , we deduce that k ( · ) F ( · ) v τ ( · ) L 1 ( Γ 3 ) for all v V using the generalized Hölder’s inequality.
Let us introduce the following subset of D ,
Λ = ζ D : ζ , v ˜ Γ 3 k ( x ) F ( x ) v ˜ τ ( x ) d Γ f o r   a l l   v ˜ γ ( V ) ;
see (18) for the definition of v ˜ τ .
It is easy to observe that λ Λ by using (17) and (20).
Lemma 2.
The set Λ is a closed convex bounded subset of D that contains 0 D .
Proof. 
It is easy to observe that 0 D Λ . In addition, the convexity can be easily obtained by means of the definition of the convex sets.
Let ( ζ n ) Λ be a convergent sequence,
ζ n ζ i n   D   a s   n .
We have to prove that ζ Λ . As the sequence ( ζ n ) converges strongly to ζ , then ζ n ζ in D as n and then
ζ n * ζ   i n   D   a s   n .
Using the definition of the weak* convergence, we infer that
ζ n , v ˜ ζ , v ˜ as   n for   all   v ˜ γ ( V ) .
Since ( ζ n ) Λ , then
ζ n , v ˜ Γ 3 k ( x ) F ( x ) v ˜ τ ( x ) d Γ for   all   v ˜ γ ( V ) .
Passing to the limit n in the previous inequality, we conclude that ζ Λ . As a result, Λ is a closed set.
Let ζ Λ be arbitrarily fixed.
ζ D = sup v ˜ γ ( V ) , v ˜ 0 γ ( V ) ζ , v ˜ v ˜ H 1 / 2 ( Γ ; R 3 ) .
Let v ˜ γ ( V ) , v ˜ 0 γ ( V ) .
ζ , v ˜ v ˜ H 1 / 2 ( Γ ; R 3 ) Γ 3 k ( x ) F ( x ) v ˜ τ ( x ) d Γ v ˜ H 1 / 2 ( Γ ; R 3 ) Γ 3 k ( x ) F ( x ) v ˜ ( x ) d Γ v ˜ H 1 / 2 ( Γ ; R 3 ) k L ( Γ 3 ) F L 2 δ 1 δ ( Γ 3 ) v ˜ L 2 δ ( Γ ; R 3 ) v ˜ H 1 / 2 ( Γ ; R 3 ) k L ( Γ 3 ) F L 2 δ 1 δ ( Γ 3 ) c ( m e a s ( Γ ) , δ ) v ˜ L 2 ( Γ ; R 3 ) v ˜ H 1 / 2 ( Γ ; R 3 ) .
As
v ˜ H 1 / 2 ( Γ ; R 3 ) v ˜ L 2 ( Γ ; R 3 ) ,
see (1), then
ζ , v ˜ v ˜ H 1 / 2 ( Γ ; R 3 ) k L ( Γ 3 ) F L 2 δ 1 δ ( Γ 3 ) c ( m e a s ( Γ ) , δ ) .
As a result,
ζ D = sup v ˜ γ ( V ) , v ˜ 0 γ ( V ) ζ , v ˜ v ˜ H 1 / 2 ( Γ ; R 3 ) k L ( Γ 3 ) F L 2 δ 1 δ ( Γ 3 ) c ( m e a s ( Γ ) , δ ) .
Therefore, Λ is a bounded set. □
Next, we claim that
c ( u , λ ) = Γ 3 k ( x ) F ( x ) u τ ( x ) d Γ .
Indeed, since c ( u , λ ) = λ , γ u = Γ 3 σ τ ( x ) · u τ ( x ) d Γ , in order to justify (21), we have to prove that
σ τ ( x ) · u τ ( x ) = k ( x ) F ( x ) u τ ( x ) a . e . on   Γ 3 .
To prove (22), we consider x Γ 3 arbitrarily fixed.
  • If u τ ( x ) = 0 then 0 = 0 and, thus, (22) holds true.
  • If u τ ( x ) 0 then
    σ τ ( x ) · u τ ( x ) = k ( x ) F ( x ) u τ ( x ) · u τ ( x ) u τ ( x ) = k ( x ) F ( x ) u τ ( x ) .
    Hence, (22) holds true in this situation too.
As a result, (21) is fulfilled.
On the other hand, by (20),
c ( u , ζ ) Γ 3 k ( x ) F ( x ) u τ ( x ) d Γ for   all   ζ Λ .
Therefore, by (21) and (23), we obtain
c ( u , ζ λ ) 0 for   all   ζ Λ .
Let ω * be the Fenchel conjugate of the constitutive function ω ,
ω * : S 3 ( , ] , ω * ( τ ) = sup ξ S 3 { τ : ξ ω ( ξ ) } .
This is a proper, convex and lower semicontinuous functional.
Notice that
( 1 β ) τ 2 ω * ( τ ) 1 4 α τ 2 for   all   τ S 3 ,
where α , β are the constants in the hypothesis H1; see [43] for a proof of (25).
Due to (25),
ω * ( τ ( · ) ) L 1 ( Ω ) for   all   τ L 2 ( Ω ; S 3 ) ,
and, due to the hypothesis H1,
ω ( τ ( · ) ) L 1 ( Ω ) for   all   τ L 2 ( Ω ; S 3 ) .
Thus,
ω ( ε ( v ( · ) ) ) L 1 ( Ω ) for   all   v V .
Let us introduce the bifunctional B : S 3 × S 3 R ,
B ( τ , μ ) = ω ( τ ) + ω * ( μ ) for   all   τ , μ S 3 .
Using Theorem 1 for X = S 3 and ϕ = ω , we deduce that B ( · , · ) is a bipotential in the sense of Definition 1.
By (27) and (26) we are lead to
B ( ε ( v ( · ) ) , τ ( · ) ) L 1 ( Ω ) for   all   v V , τ L s 2 ( Ω ; S 3 ) .
This regularity allows us to define a form b ( · , · ) as follows,
b : V × L s 2 ( Ω ; S 3 ) R b ( v , μ ) = Ω B ( ε ( v ( x ) ) , μ ( x ) ) d x .
Since σ ( x ) ω ( ε ( u ) ( x ) ) in Ω and B ( · , · ) defined in (28) is a bipotential, then a.e. x Ω ,
B ( ε ( u ) ( x ) , σ ( x ) ) = σ ( x ) : ε ( u ) ( x ) , B ( ε ( v ) ( x ) , μ ( x ) ) μ ( x ) : ε ( v ) ( x ) .
After integration over Ω , we obtain
b ( u , σ ) = ( σ , ε ( u ) ) L s 2 ( Ω ; S 3 ) ,
b ( v , μ ) ( μ , ε ( v ) ) L s 2 ( Ω ; S 3 ) for   all   v V , μ L s 2 ( Ω ; S 3 ) .
In particular, setting in (31) v = u and μ = σ , respectively, the following inequalities can be written:
b ( u , μ ) ( μ , ε ( u ) ) L s 2 ( Ω ; S 3 ) for   all   μ L s 2 ( Ω ; S 3 ) ;
b ( v , σ ) ( σ , ε ( v ) ) L s 2 ( Ω ; S 3 ) for   all   v V .
Hence, by (32) and (30), we arrive at
b ( u , μ ) b ( u , σ ) ( μ σ , ε ( u ) ) L s 2 ( Ω ; S 3 ) .
Consider now a variable subset of L s 2 ( Ω ; S 3 ) defined as follows: given ζ Λ ,
Σ ( ζ ) = { μ L s 2 ( Ω ; S 3 ) : ( μ , ε ( v ) ) L s 2 ( Ω ; S 3 ) + c ( v , ζ ) = ( f , v ) V for   all   v V } .
Lemma 3.
Let ζ Λ . The set Σ ( ζ ) is a nonempty, convex, closed, unbounded subset of L s 2 ( Ω ; S 3 ) .
Proof. 
Let ζ Λ and let us define
φ ζ : V R , φ ζ ( v ) = c ( v , ζ ) .
As c ( · , · ) is a bilinear form, then the linearity of φ ζ is obvious. Let us prove its continuity. In order to do this, we have to prove that there exists K > 0 such that
| φ ζ ( v ) | K v V for   all   v V .
Indeed, let v V .
| φ ζ ( v ) | = | c ( v , ζ | c ˜ c ˜ K ζ D v V .
We can set K = c ˜ c ˜ K , where c ˜ is the constant in (2) and c ˜ K is the constant in (5). As a result, φ ζ is a linear and continuous map. Then, due to Riesz’s representation theorem, there exists a unique G ζ V such that
c ( v , ζ ) = ( G ζ , v ) V for   all   v V .
Hence,
( f , v ) V c ( ζ , v ) = ( f , v ) V ( G ζ , v ) V = ( f G ζ , v ) V ,
and keeping in mind the definition of the inner product on V, (4), actually we can write
( f , v ) V c ( ζ , v ) = ( ε ( f G ζ ) , ε ( v ) ) L s 2 ( Ω ; S 3 ) .
Let us take μ = ε ( f G ζ ) to conclude that Σ ( ζ ) is a nonempty subset of L s 2 ( Ω ; S 3 ) . The convexity can be easily proved by using the definition of the convex sets.
In order to prove that Σ ( ζ ) is a closed subset of L s 2 ( Ω ; S 3 ) , we consider ( μ n ) Σ ( ζ ) such that μ n μ in L s 2 ( Ω ; S 3 ) as n and we prove that μ Σ ( ζ ) . As
( μ n , ε ( v ) ) L s 2 ( Ω ; S 3 ) + c ( v , ζ ) = ( f , v ) V for   all   v V
and
( μ n , ε ( v ) ) L s 2 ( Ω ; S 3 ) ( μ , ε ( v ) ) L s 2 ( Ω ; S 3 )   as   n   for   each   v V
then, passing to the limit n in (35), we obtain
( μ , ε ( v ) ) L s 2 ( Ω ; S 3 ) + c ( v , ζ ) = ( f , v ) V for   all   v V .
Hence, μ Σ ( ζ ) and, thus, Σ ( ζ ) is a closed set.
Finally, let us prove that Σ ( ζ ) is an unbounded subset. Indeed, there exist at least one sequence ( μ n ) n Λ such that μ n L s 2 ( Ω ; S 3 ) as n . We can construct such a sequence as follows: we take τ ε ( V ) and for each positive integer n we define
μ n = ε ( f G ζ ) + n τ ,
where ε ( V ) = { τ L s 2 ( Ω ; S 3 ) | ( τ , σ ) L s 2 ( Ω ; S 3 ) = 0 for   all   σ ε ( V ) } , with ε ( V ) = { ε ( v ) | v V } . Recall that ε ( V ) and ε ( V ) are closed subspaces of the space L s 2 ( Ω ; S 3 ) and L s 2 ( Ω ; S 3 ) = ε ( V ) ε ( V ) , see, e.g., Theorem 1.16 in [8]. As f G ζ V and τ ε ( V ) , by using (36) it follows that μ n L s 2 ( Ω ; S 3 ) for all positive integer n . Moreover,
( μ n , ε ( v ) ) L s 2 ( Ω ; S 3 ) + ζ , γ v = ( ε ( f ) , ε ( v ) ) L s 2 ( Ω ; S 3 ) .
Taking into consideration (4) we can write
( μ n , ε ( v ) ) L s 2 ( Ω ; S 3 ) + ζ , γ v = ( f , v ) V .
Thus, ( μ n ) n Σ ( ζ ) .
Furthermore, μ n L s 2 ( Ω ; S 3 ) as n . Indeed,
μ n L s 2 ( Ω ; S 3 ) = n ε ( f ) ε ( G ζ ) n + τ L s 2 ( Ω ; S 3 )
and
ε ( f ) ε ( G ζ ) n + τ
is a bounded sequence. Therefore, the subset Σ ( ζ ) is unbounded. □
We observe that σ Σ ( λ ) .
Let μ Σ ( λ ) . Then, ( μ , ε ( u ) ) L s 2 ( Ω ; S 3 ) = ( ε ( f ) ε ( G λ ) , u ) V = ( σ , ε ( u ) ) L s 2 ( Ω ; S 3 ) . Hence,
( μ σ , ε ( u ) ) L s 2 ( Ω ; S 3 ) = 0 .
By (34), we obtain
b ( u , μ ) b ( u , σ ) 0 for   all   μ Σ ( λ ) .
On the other hand, by (33) and (30) combined with (16), (17) and (19), we can write
b ( v , σ ) b ( u , σ ) ( σ , ε ( v ) ε ( u ) ) L s 2 ( Ω ; S 3 ) = ( f , v u ) V λ , γ v γ u = ( f , v u ) V c ( v u , λ ) .
As a consequence,
b ( v , σ ) b ( u , σ ) + c ( v u , λ ) ( f , v u ) V for   all   v V .
Keeping in mind (38), (24) and (37), we are led to the following variational problem.
Problem 2.
Find u V , λ Λ , σ Σ ( λ ) such that
b ( v , σ ) b ( u , σ ) + c ( v u , λ ) ( f , v u ) V f o r   a l l   v V c ( u , ζ λ ) 0 f o r   a l l   ζ Λ b ( u , μ ) b ( u , σ ) 0 f o r   a l l   μ Σ ( λ ) .
Each solution ( u , λ , σ ) V × Λ × Σ ( λ ) of Problem 2 is called a weak solution of Problem 1.
Keeping in mind (29) and (28), the bifunctional b ( · , · ) can be written by means of two functionals J ( · ) and J * ( · ) as follows,
b ( v , μ ) = J ( v ) + J * ( μ ) ,
where
J : V R , J ( v ) = Ω ω ( ε ( v ) ( x ) ) d x
and
J * : L s 2 ( Ω ; S 3 ) R , J * ( μ ) = Ω ω * ( μ ( x ) ) d x .
Then, we observe that
b ( v , σ ) b ( u , σ ) = J ( v ) J ( u ) , b ( u , μ ) b ( u , σ ) = J * ( μ ) J * ( σ ) .
In consequence, Problem 2 can be equivalently written as follows.
Problem 3.
Find u V , λ Λ , σ Σ ( λ ) such that
J ( v ) J ( u ) + c ( v u , λ ) ( f , v u ) V f o r   a l l   v V
c ( u , ζ λ ) 0 f o r   a l l   ζ Λ J * ( μ ) J * ( σ ) 0 f o r   a l l   μ Σ ( λ ) .
The weak solvability of Problem 1 will be studied by means of the variational formulation stated in Problem 3.

4. Well-Posedness and a Convergence Result

This section is devoted to the solvability of Problem 3.
Theorem 4.
Assume that the hypotheses H1–H4 hold true. Then, Problem 3 has at least one solution ( u , λ , σ ) V × Λ × Σ ( λ ) . If ω is, in addition, strictly convex, there is uniqueness in the first component.
Proof. 
Let us introduce the following bifunctional:
L : V × Λ R , L ( v , ζ ) = J ( v ) ( f , v ) V + c ( v , ζ ) .
A pair ( u , λ ) V × Λ verifies (39) and (40) if and only if it is a saddle point of the bifunctional L ( · , · ) , i.e.,
L ( u , ζ ) L ( u , λ ) L ( v , λ ) , for   all   v V , for   all   ζ Λ .
Indeed, (40) is equivalent with the first inequality in the chain above. On the other hand, using the definition of L ( · , · ) , by a similar technique with that used in [44], we deduce that (39) is equivalent with the second inequality in (41).
Keeping in mind H1, we immediately conclude that J ( · ) is a convex lower semicontinuous functional such that
J ( v ) α v V 2 for   all   v V .
Furthermore, according to Lemma 1, c ( · , · ) is a bilinear continuous form. Therefore, L ( · , · ) fulfills the conditions in Theorem 3. Moreover, according to Lemma 2, Λ is a closed convex bounded set which contains 0 D . In addition, we note that
L ( v , 0 D ) = J ( v ) + c ( v , 0 D ) ( f , v ) X α v V 2 f V v V ,
which allows us to write
lim v V L ( v , 0 D ) = .
Therefore, applying Theorem 3, we conclude that the functional L ( · , · ) has at least one saddle point ( u * , λ * ) V × Λ .
Let us introduce the set Σ ( λ * ) . According to Lemma 3, this is a nonempty closed convex unbounded subset of L s 2 ( Ω ; S 3 ) . Moreover, due to (25), we infer that J * ( · ) is convex, lower semicontinuous and
J * ( μ ) ( 1 β ) μ L s 2 ( Ω ; S 3 ) 2 .
Using Theorem 2, we conclude that J * ( · ) has a unique minimizer σ * on Σ ( λ * ) . Consequently, the triple ( u * , λ * , σ * ) V × Λ × Σ ( λ * ) is a solution of Problem 3.
Finally, if ω is, in addition, strictly convex, using Theorem 3 we immediately conclude that Problem 3 has at least one solution which is unique in its first component. □
In order to obtain a global uniqueness result, one option could be to perturb Problem 3 as follows.
Problem 4.
Let ϵ > 0 . Find u ϵ V , λ ϵ Λ , σ ϵ Σ ( λ ϵ ) such that
J ( v ) J ( u ϵ ) + c ( v u ϵ , λ ϵ ) ( f , v u ϵ ) V f o r   a l l   v V
c ( u ϵ , ζ λ ϵ ) ϵ ( ζ D 2 λ ϵ D 2 ) 0 f o r   a l l   ζ Λ
J * ( μ ) J * ( σ ϵ ) 0 f o r   a l l   μ Σ ( λ ϵ ) .
Let us introduce the following perturbed bifunctional.
L ϵ : V × Λ R , L ( v , ζ ) = J ( v ) ( f , v ) V + c ( v , ζ ) ϵ ζ D 2 ( ϵ > 0 ) .
Using similar techniques with those used in [44], it can be verified that a pair ( u ϵ , λ ϵ ) V × λ verifies (42) and (43) if and only if it is a saddle point of the bifunctional L ϵ ( · , · ) , i.e.,
L ϵ ( u ϵ , ζ ) L ϵ ( u ϵ , λ ϵ ) L ϵ ( v , λ ϵ ) , for   all   v V , ζ Λ .
We observe that L ϵ ( · , · ) is strictly concave in the second argument.
Theorem 5.
Let ϵ > 0 . Assume H1–H4 hold true. In addition, we assume that ω and ω * are strictly convex functionals. Then, Problem 4 has a unique solution ( u ϵ , λ ϵ , σ ϵ ) V × Λ × Σ ( λ ϵ ) . Moreover, if ω is Lipschitz continuous of rank L then
u ϵ V f V α ;
λ ϵ D f V + L α c ;
σ ϵ L s 2 ( Ω ; S 3 ) 1 4 α ( 1 β ) f V + M c ( f V + L ) α c .
Proof. 
Let ϵ > 0 . We observe that L ϵ ( · , · ) is strictly convex in the first argument and strictly concave in the second argument. We apply Theorem 3 in order to conclude that there exists a unique ( u ϵ * , λ ϵ * ) V × Λ such that (42) and (43) are fulfilled. Let Σ ( λ ϵ * ) . As J * ( · ) is strictly convex, lower semicontinuous and coercive we deduce that it has a unique minimizer σ ϵ * Σ ( λ ϵ * ) . The triple ( u ϵ * , λ ϵ * , σ ϵ * ) V × Λ × Σ ( λ ϵ * ) is the unique solution of Problem 4.
In order to prove the boundedness of the solution, let us set v = 0 V in (42) and ζ = 0 D in (43). With this choice, (42) and (43) lead us to
J ( u ϵ ) ( f , u ϵ ) V ,
since, due to H1, J ( 0 V ) = 0 . Additionally, from this,
α u ϵ V 2 f V u ϵ V .
Hence, (45) holds true.
By (42), setting v = u ϵ w with w V arbitrarily chosen, we immediately obtain
c ( w , λ ϵ ) ( f , w ) V + J ( u ϵ w ) J ( u ϵ ) .
As ω is Lipschitz continuous of rank L, J is Lipschitz continuous of the same rank L . Therefore,
c ( w , λ ϵ ) ( f V + L ) w V for   all   w V .
Let w 0 V . Then
sup w V , w 0 V c ( w , λ ϵ ) w V f V + L .
By using the inf-sup property of the form c ( · , · ) , we can write
α c λ ϵ D f V + L .
As a result, we obtain (46).
According to (44) and keeping in mind (25),
( 1 β ) σ ϵ L s 2 ( Ω ; S 3 ) 2 J * ( σ ϵ ) J * ( μ ) 1 4 α μ L s 2 ( Ω ; S 3 ) 2 for   all   μ Σ ( λ ϵ ) .
Let G λ ϵ be the unique element of V such that
c ( v , λ ϵ ) = ( G λ ϵ , v ) V for   all   v V .
Using G λ ϵ we define
μ = ε ( f ) ε ( G λ ϵ ) .
Then,
μ L s 2 ( Ω ; S 3 ) ε ( f ) L s 2 ( Ω ; S 3 ) + ε ( G λ ϵ ) L s 2 ( Ω ; S 3 ) = f V + G λ ϵ V .
As
G λ ϵ V = sup v V , v 0 V ( G λ ϵ , v ) V v V = sup v V , v 0 V c ( v , λ ϵ ) v V M c λ ϵ D ,
then
μ L s 2 ( Ω ; S 3 ) f V + M c λ ϵ D .
and from this, keeping in mind (46), we deduce that
μ L s 2 ( Ω ; S 3 ) f V + M c ( f V + L ) α c .
Using now (48), we can write
σ ϵ L s 2 ( Ω ; S 3 ) 1 4 α ( 1 β ) μ L s 2 ( Ω ; S 3 ) .
Combining this last inequality with the previous one, we immediately obtain (47). □
With these preliminaries, an approximation result can be obtained. Everywhere below, we keep the assumptions H1, H2, H3, H4. In addition, we assume that ω and ω * are strictly convex. In addition, we assume that ω is Lipschitz continuous of rank L and ω * is upper semicontinuous.
Theorem 6
(A convergence result). Let ϵ > 0 and let ( u ϵ , λ ϵ , σ ϵ ) be the unique solution of Problem 4. Then, passing eventually to a subsequence ( u ϵ , λ ϵ , σ ϵ ) of the sequence ( u ϵ , λ ϵ , σ ϵ ) , the subsequence ( u ϵ , λ ϵ , σ ϵ ) is weakly convergent to a solution of Problem 3 as ϵ 0 .
Proof. 
Let ϵ > 0 and let ( u ϵ , λ ϵ , σ ϵ ) be the unique solution of Problem 4. Let G λ ϵ be the unique element of V such that c ( v , λ ϵ ) = ( G λ ϵ , v ) for all v V . As
G λ ϵ V = sup v V , v 0 V ( G λ ϵ , v ) V v V = sup v V , v 0 V c ( v , λ ϵ ) v V M c λ ϵ D M c ( f V + L ) α c ,
then, G λ ϵ is a bounded sequence. This conclusion together with the boundedness results (45)–(47) allow us to assert that ( u ϵ , λ ϵ , σ ϵ , G ϵ ) is a bounded sequence. Thus, passing eventually to a subsequence ( u ϵ , λ ϵ , σ ϵ , G λ ϵ ) of the sequence ( u ϵ , λ ϵ , σ ϵ , G λ ϵ ) , we can write
u ϵ u *   in   V   as   ϵ 0 ; λ ϵ λ *   in   D   as   ϵ 0 ; σ ϵ σ *   in   L s 2 ( Ω ; S 3 )   as   ϵ 0 ; G λ ϵ G *   in   V   as   ϵ 0 .
To proceed, we prove that ( u * , λ * , σ * ) is a solution of Problem 3.
Passing to the limit ϵ 0 in (42) and (43), we observe that ( u * , λ * ) V × D verifies (39) and the inequality in (40). Moreover, Λ being a convex and closed set, it is weakly closed too. Therefore, λ * Λ . Thus, (40) is verified.
Next, we claim that
G λ ϵ G λ *   in   V   as   ϵ 0 ,
where G λ * is the unique element of V such that c ( v , λ * ) = ( G λ * , v ) for all v V . Indeed, (49) implies that λ ϵ * λ * . Therefore,
λ ϵ , v ˜ λ * , v ˜   as   ϵ 0 for   all   v ˜ γ ( V ) .
Let v V be arbitrarily fixed. As λ ϵ , γ v λ * , γ v as ϵ for all v V , then
c ( v , λ ϵ ) c ( v , λ * )   as   ϵ 0 for   all   v V .
and from this,
( G λ ϵ , v ) V ( G λ * , v ) V   as   ϵ 0 for   all   v V .
Hence, (50) holds true.
Next, by (50) we can prove that
ε ( G λ ϵ ) ε ( G λ * )   in   L s 2 ( Ω ; S 3 )   as   ϵ 0 .
Indeed, ε being a continuous operator, it is also a lower semicontinuous and upper semicontinuous operator. On the other hand, ε is a linear operator so it is also a convex operator. Therefore, ε is a weakly lower semicontinuous operator and a weakly upper semicontinuous operator as well. Hence, keeping in mind (50), we can write
lim sup ε ( G λ ϵ ) ε ( G λ * ) lim inf ε ( G λ ϵ )   as   ϵ 0 .
As a result, (51) is true.
Let μ Σ ( λ * ) . Then, there exists μ 0 ε ( V ) such that
μ = ε ( f ) ε ( G λ * ) + μ 0 .
We claim that there exists ( μ ϵ ) L s 2 ( Ω ; S 3 ) such that
μ ϵ Σ ( λ ϵ ) , μ ϵ μ i n L s 2 ( Ω ; S 3 )   as   ϵ 0 .
Indeed, let us define
μ ϵ = ε ( f ) ε ( G λ ϵ ) + μ 0 .
Clearly, μ ϵ Σ ( λ ϵ ) .
Due to (51), we immediately observe that (52) holds true.
Passing to the superior limit ϵ 0 in
J * ( μ ϵ ) J * ( σ ϵ ) 0 ,
we obtain
J * ( μ ) J * ( σ * ) 0 .
As μ was arbitrarily chosen in Σ ( λ * ) , we conclude that
J * ( μ ) J * ( σ * ) 0 for   all   μ Σ ( λ * ) .
It remains to justify that
σ * Σ ( λ * ) .
Indeed, as
σ ϵ σ *   as   ϵ 0 ,
λ ϵ λ *   as   ϵ 0 ,
keeping in mind the definition of Σ ( λ ϵ ) and Σ ( λ * ) , we obtain the conclusion passing to the limit ϵ 0 in
( σ ϵ , ε ( v ) ) L 2 ( Ω ; S 3 ) + c ( v , λ ϵ ) = ( f , v ) V .

5. Conclusions and Final Comments

In the present paper, we address a frictional contact model with prescribed normal stress. We deliver a new weak formulation that is a three-field variational formulation governed by a bipotential related to the constitutive function ω and a Lagrange multiplier related to the friction force σ τ . We establish existence, uniqueness, boundedness and convergence results. Theorem 6 indicates us that the unique solution of the perturbed Problem 4 helps us to approximate a weak solution of Problem 1. Delivering uniqueness results by omitting a perturbation technique is left open.
The advantage of the approach we propose is twofold. On the one hand, the new approach allows the inclusion of the friction force in the unknown in addition to the displacement field and the Cauchy stress tensor. On the other hand, the qualitative analysis we perform allows moving on to the quantitative analysis in order to efficiently approximate the triple weak solutions.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the referees for their careful evaluation and valuable suggestions.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Matei, A. A variational approach via bipotentials for a class of frictional contact problems. Acta Appl. Math. 2014, 134, 45–59. [Google Scholar] [CrossRef]
  2. Doghri, I. Mechanics of Deformable Solids; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
  3. Han, W.; Sofonea, M. Quasistatic Contact Problems in Viscoelasticity and Viscoplasticity; Studies in Advanced Mathematics; American Mathematical Society, International Press: Providence, RI, USA, 2002. [Google Scholar]
  4. Kikuchi, N.; Oden, J.T. Contact Problems in Elasticity: A Study of Variational Inequalities and Finite Element Methods; SIAM: Philadelphia, PA, USA, 1988. [Google Scholar]
  5. Laursen, T. Computational Contact and Impact Mechanics: Fundamentals of Modeling Interfacial Phenomena in Nonlinear Finite Element Analysis; Springer Science and Business Media: New York, NY, USA, 2013. [Google Scholar]
  6. Migorski, S.; Ochal, A.; Sofonea, M. Nonlinear Inclusions and Hemivariational Inequalities. Models and Analysis of Contact Problems; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  7. Sofonea, M.; Matei, A. Variational Inequalities with Applications. A Study of Antiplane Frictional Contact Problems; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  8. Sofonea, M.; Matei, A. Mathematical Models in Contact Mechanics; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  9. Wriggers, P. Computational Contact Mechanics, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  10. Buliga, M.; de Saxcé, G.; Vallxexe, C. Bipotentials for Non-Monotone Multivalued Operators: Fundamental Results and Applications. Acta Appl. Math. 2010, 110, 955–972. [Google Scholar] [CrossRef] [Green Version]
  11. Buliga, M.; de Saxcé, G.; Vallxexe, C. Non-maximal cyclically monotone graphs and construction of a bipotential for the Coulomb’s dry friction law. J. Convex Anal. 2010, 17, 81–94. [Google Scholar]
  12. Buliga, M.; de Saxcé, G.; Vallxexe, C. Existence and construction of bipotentials for graphs of multivalued laws. J. Convex Anal. 2008, 15, 87–104. [Google Scholar]
  13. de Saxcé, G. Une généralisation de l’inégalité de Fenchel et ses applications aux lois constitutives. C. R. Acad. Sci. Paris 1992, 314, 125–129. [Google Scholar]
  14. de Saxcé, G. The bipotential method, a new variational and numerical treatment of the dissipative laws of materials. In Proceedings of the 10th International Conference on Mathematical and Computer Modelling and Scientific Computing, Boston, MA, USA, 5–8 July 1995. [Google Scholar]
  15. de Saxcé, G.; Feng, Z.-Q. The Bipotential Method: A Constructive Approach to Design the Complete Contact Law with Friction and Improved Numerical Algorithms. Mathl. Comput. Model. 1998, 28, 225–245. [Google Scholar] [CrossRef]
  16. Braess, D. Finite Elements. Theory, Solvers, and Applications in Solid Mechanics; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar] [CrossRef]
  17. Brezzi, F.; Fortin, M. Mixed and Hybrid Finite Element Methods; Springer: New York, NY, USA, 1991. [Google Scholar]
  18. Ekeland, I.; Témam, R. Convex analysis and variational problems. In Studies in Mathematics and Its Applications; North-Holland Publishing Company: Amsterdam, The Netherlands, 1976; Volume 1. [Google Scholar]
  19. Haslinger, J.; Hlaváček, I.; Nečas, J. Numerical methods for unilateral problems in solid mechanics. In Handbook of Numerical Analysis; Ciarlet, P.G., Lions, J.L., Eds.; North-Holland: Amsterdam, The Netherlands, 1996; Volume IV, pp. 313–485. [Google Scholar]
  20. Hüeber, S.; Matei, A.; Wohlmuth, B. Efficient algorithms for problems with friction. SIAM J. Sci. Comput. 2007, 29, 70–92. [Google Scholar] [CrossRef]
  21. Adams, R.A. Sobolev Spaces; Academic Press: New York, NY, USA, 1975. [Google Scholar]
  22. Brézis, H. Functional Analysis, Sobolev Spaces and Partial Differential Equations; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  23. Grisvard, P. Elliptic Problems in Nonsmooth Domains; Pitman: London, UK, 1985. [Google Scholar]
  24. Kufner, A.; John, O.; Fučik, S. Function Spaces; Noordhoff International Publishing: Groningen, The Netherlands, 1977. [Google Scholar]
  25. Megginson, R.E. An Introduction to Banach Space Theory (Graduate Texts in Mathematics); Springer: New York, NY, USA, 1998; Volume 183. [Google Scholar]
  26. Monk, P. Numerical Mathematics and Scientific Computation. Finite Element Methods for Maxwell’s Equations; Oxford University Press: Oxford, UK, 2003. [Google Scholar]
  27. Lions, J.-L.; Magenes, E. Non-Homogeneous Boundary Value Problems and Applications I; Springer: New York, NY, USA, 1972. [Google Scholar]
  28. Han, W.; Matei, A. Minimax Principles for Elliptic Mixed Hemivariational-Variational Inequalities. Nonlinear Anal. Real World Appl. 2022, 64, 103448. [Google Scholar] [CrossRef]
  29. Matei, A.; Osiceanu, M. Two-Field Weak Solutions for a Class of Contact Models. Mathematics 2022, 10, 369. [Google Scholar] [CrossRef]
  30. Matei, A.; Osiceanu, M. Two-field variational formulations for a class of nonlinear mechanical models. Math. Mech. Solids 2022, 27, 2532–2547. [Google Scholar] [CrossRef]
  31. Goebel, R.; Rockafellar, R.T. Local strong convexity and local Lipschitz continuity of the gradient of convex functions. J. Convex Anal. 2008, 15, 263–270. [Google Scholar]
  32. Niculescu, C.P.; Persson, L.-E. Convex Functions and Their Applications. A Contemporary Approach; CMS Books in Mathematics; Springer: New York, NY, USA, 2006; Volume 23. [Google Scholar]
  33. Phelps, R. Convex Functions, Monotone Operators and Differentiability, 2nd ed.; Lecture Notes in Math; Springer: Berlin/Heidelberg, Germany, 1993; Volume 1364. [Google Scholar]
  34. Rockafellar, R.T. Convex Analysis; Princeton University Press: Princeton, NJ, USA, 1997. [Google Scholar]
  35. Struwe, M. Variational Methods: Applications to Nonlinear Partial Differential Equations and Hamiltonian Systems; Springer: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  36. Nečas, J. Direct Methods in the Theory of Elliptic Equations; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  37. Arendt, W.; Kreuter, M. Mapping theorems for Sobolev-spaces of vector-valued functions. Stud. Math. 2018, 240, 275–299. [Google Scholar] [CrossRef] [Green Version]
  38. Sofonea, M. Problèmes Mathématiques en Elasticité et Viscoplasticité; Cours de DEA de Mathématiques Appliquées à l’Université Blaise Pascal (Clermind Ferrand II), Laboratoire de Mathématiques Appliqées et URA-CNRS: Aubière, France, 1993. [Google Scholar]
  39. Nečas, J.; Hlaváček, I. Mathematical Theory of Elastic and Elastico-Plastic Bodies: An Introduction; Elsevier Scientific Publishing Company: Amsterdam, The Netherlands; Oxford, UK; New York, NY, USA, 1981; Volume 79. [Google Scholar]
  40. Matei, A.; Ciurcea, R. Contact problems for nonlinearly elastic materials: Weak solvability involving dual Lagrange multipliers. ANZIAM J. 2010, 52, 160–178. [Google Scholar] [CrossRef]
  41. Duvaut, G.; Lions, J.-L. Inequalities in Mechanics and Physics; Springer: Berlin/Heidelberg, Germany, 1976. [Google Scholar]
  42. Panagiotopoulos, P.D. Inequality Problems in Mechanics and Applications; Birkhäuser: Boston, MA, USA, 1985. [Google Scholar]
  43. Matei, A. A variational approach via bipotentials for unilateral contact problems. J. Math. Anal. Appl. 2013, 397, 371–380. [Google Scholar] [CrossRef]
  44. Cojocaru, M.C.; Matei, A. On a class of saddle point problems and convergence results. Math. Model. Anal. 2020, 25, 608–621. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Matei, A.C. A Three-Field Variational Formulation for a Frictional Contact Problem with Prescribed Normal Stress. Fractal Fract. 2022, 6, 651. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6110651

AMA Style

Matei AC. A Three-Field Variational Formulation for a Frictional Contact Problem with Prescribed Normal Stress. Fractal and Fractional. 2022; 6(11):651. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6110651

Chicago/Turabian Style

Matei, Andaluzia Cristina. 2022. "A Three-Field Variational Formulation for a Frictional Contact Problem with Prescribed Normal Stress" Fractal and Fractional 6, no. 11: 651. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6110651

Article Metrics

Back to TopTop