1. Introduction
Variational–hemivariational inequalities represent a powerful mathematical tool in the study of nonlinear boundary value problems. Their study is motivated by various applications in Physics, Mechanics and Engineering Sciences, among others. In contrast with variational inequalities (which are governed by convex functions) and hemivariational inequalities (which are governed by nonsmooth locally Lipschitz functions which could be nonconvex), variational–hemivariational inequalities are governed by both convex and locally Lipschitz functions. As a consequence, they have both a convex and nonconvex structure and, therefore, their study is carried out by using arguments on both convex and nonsmooth analysis.
Introduced in the pioneering work of Panagiotopoulos [
1], the theory of variational–hemivariational inequalities grew up rapidly, as shown in [
2,
3,
4,
5] and the references therein. It includes existence, uniqueness and numerical approximation results, obtained in the study of different classes of inequalities, by using various methods and functional arguments. Reference in the field include [
6,
7,
8,
9,
10,
11]. Among the inequalities studied in these papers we distinguish the class of elliptic, the class of time-dependent and the class of evolutionary variational–hemivariational inequalities. A variational–hemivariational inequality is said to be elliptic if it does not involve the time variable; it is said to be time-dependent if both the data and the solution depend on time but no time-derivatives of the solution appear in its statement; finally, a variational–hemivariational inequality is said to be evolutionary if it is formulated in terms of the derivative of the unknown function.
A first example of elliptic variational-hemivariational inequality is the following: find
u such that
Here and below, H represents a real Hilbert space endowed with the inner product and the associated norm , K is a nonempty subset of H, is a nonlinear operator, and are given functions and, finally, . The function is assumed to be convex while the function g is assumed to be locally Lipschitz. Moreover, notation represents the Clarke directional derivative of g at the point u, in the direction v.
A second example of elliptic variational-hemivariational inequality is the following: find
u such that
Note that, in contrast with (
1), here the functions
and
j are defined on the product space
, that is,
and
. The function
is assumed to be convex with respect to the second variable,
j is assumed to be locally Lipschitz with respect to the second argument and notation
represents the Clarke directional derivative of
at the point
u, in the direction
v.
A special case of time-dependent variational-hemivariational inequalities is given by the so-called history-dependent variational-hemivariational inequalities. A typical example is the following: find a function
such that
Note that in (
3) and below in this paper
and
are operators defined on the space of continuous functions defined on
with values
H, denoted in what follows by
. Moreover, for any function
we use the shorthand notation
and
, to represent the value of the functions
and
at the point
, that is,
and
.
Inequality problems of the form (
1), (
2), arise in the study of mathematical models which describe the equilibrium of elastic body in frictionless and frictional contact with a foundation, respectively. Moreover, inequality problems of the form (
3) arise in the study of mathematical models of contact with elastic or viscoelastic materials, in which memory effects are taken into consideration, either in the constitutive law or in the contact conditions. References in the field are the books [
4,
5] as well as the survey article [
2]. Moreover, it is worth noting that variational–hemivariational inequalities arise in the study of complex fluids and history-dependent viscoelastic and elasto-viscoplastic models. A comprehensive reference in the field is the book [
12]. There, an introduction to the modeling of complex fluids is provided, up-to-date mathematical and numerical analysis of the corresponding equations can be found, together with several numerical algorithms for the approximation of the solutions. Furthermore, subdifferential operators have been used in [
13] in the study of various magnetorheological mixtures composed of a fluid and a solid continuum.
Existence and uniqueness results in the study of elliptic variational–hemivariational inequality have been obtained in many papers, under different assumptions on the data. For instance, a surjectivity result for pseudomonotone multivalued operators was used in [
14] in order to obtain the unique solvability of inequality (
1). There, the operator
A was assumed to be pseudomonotone and strongly monotone and the Clarke subdifferential of the function
j was assumed to satisfy a growth condition. The method used in [
14] can be used in the study of inequality (
2), as shown in [
5], for instance. Recently, problem (
1) was considered in [
15], under the assumption that
A a strongly monotone Lipschitz continuous operator and
is a continuous convex function. The unique solvability of the problem was obtained by using a minimization principle which avoids any pseudomonotonicity argument.
Motivated by the importance of the topic in both pure and applied mathematics, in this paper we introduce a new approach which allows us to prove existence, uniqueness and convergence results for variational–hemivariational inequalities in Hilbert spaces. The novelty of the results we present here arises in the fact that the approach we use is based on arguments of multivalued maximal operators in Hilbert spaces and fixed point. It can be used for various classes of elliptic or history-dependent variational–hemivariational inequalities. Nevertheless, for simplicity, we restrict ourselves to the study of inequalities (
1)–(
3), which represent three relevant examples. Our results are obtained under assumptions which are slightly different from those used in [
5,
14,
15] and, therefore, they complete the results obtained in these references. For instance, here
j is a bifunction, no growth assumption on its subdifferential is assumed and the smallness assumptions involving the constants
,
,
,
(related to the data
A,
,
j) are relaxed. Relaxing this assumption was possible by using the Browder–Godhe–Kirk fixed point argument instead of the classical Banach fixed point principle. For all these reasons we believe that our results contribute to a better knowledge of the structure of variational–hemivariational inequalities and, in addition, they open the way to the approach of the solution by using various iterative methods.
The rest of the manuscript is organized as follows. In
Section 2 we recall some preliminary material. In
Section 3 we consider a variational–hemivariational inequality which depends on two parameters. We use arguments of convex and nonsmooth analysis in order to prove that this inequality is governed by a maximal monotone operator. This allows us to obtain various properties for the resolvent of this operator, which have interest in their own. We use these properties in
Section 4,
Section 5 and
Section 6 in order to deduce existence and uniqueness results for elliptic and history-dependent variational–hemivariational inequalities of the form (
1), (
2) and (
3), respectively. In addition to the properties of the resolvent operator, our proofs are based on equivalence and fixed point arguments. We also introduce several iterative methods in solving these inequalities and deduce various convergence results. Finally, in
Section 7 we present some concluding remarks.
2. Preliminaries
The results we present in this section can be found in many books and surveys, including [
5,
16,
17,
18,
19]. For this reason we present them without proofs. Everywhere below
H represents a real Hilbert space endowed with the inner product
and the associated norm
. We use the symbols “→” and “⇀” to denote the strong and the weak convergence in the space
H and employ the notation
for the space
H equipped with the weak topology. The limits, lower limits and upper limits are considered as
, even if we do not mention it explicitly. Moreover, we use
for the interior of the set
, in the strong topology of
H. Finally, we denote by
the identity map of
H, by
the zero element of
H and by
the set of parts of
H. We start with the following definitions for single-valued operators.
Definition 1. The operator is said to be:
- (a)
demicontinuous if in H implies in H;
- (b)
strongly monotone if there exists constant such that - (c)
Lipschitz continuous if there exists constant such that
Definition 2. Let . The operator is said to be:
- (a)
nonexpansive on K if there exists a constant such that - (b)
a contraction if it is nenexpansive on K with constant .
In this paper, in addition to the well-known Banach contraction principle we shall use the following Browder–Godhe–Kirk fixed point Theorem, proved in [
20], (p. 55).
Theorem 1. Let K be a nonempty closed bounded convex subset of the Hilbert space H and let be a nonexpansive operator. Then A has at least one fixed point.
We now proceed with some results concerning multivalued operators defined on the space
H. To this end we recall that, given a multivalued operator
, its domain
, range
and graph
are the sets defined by
Definition 3. The operator is said to be:
- (a)
- (b)
relaxed monotone if there exists constant such that - (c)
maximal monotone if it is monotone and, for any , the following implication holds:
There is a close connection between the property of maximal monotonicity of T and the surjectivity property of the operator with . The fundamental result in this direction is the celebrated theorem of Minty that we recall below.
Theorem 2. Let be a maximal monotone operator and let . Then . Moreover, for any there exists a unique element such that .
Theorem 2 allows us to consider the resolvent operator
defined by
for any
. In other words,
is the inverse of the operator
, i.e.,
. Note that the resolvent operator exists for each
and is a single valued operator.
Next, we recall two sufficient conditions which guarantee the maximal monotonicity of a multivalued operator.
Proposition 1. Assume that is a monotone operator such that for every , the set is nonempty convex and weakly closed. Moreover, assume that for all , the graph of the mapping is closed in . Then the operator T is maximal monotone.
Proposition 2. Let be two maximal monotone operators such that . Then is a maximal monotone operator, too.
We now proceed with the definition and the properties of the Clarke subdifferential of locally Lipschitz functions.
Definition 4. The Clarke directional derivative of the locally Lipschitz function at the point in the direction is defined by The Clarke subdifferential of j is the multivalued operator defined by For the Clarke subdifferential and directional derivative we have the following properties.
Proposition 3. Let be a locally Lipschitz function. Then:
- (a)
is a nonempty convex and compact subset of H, for all ;
- (b)
the graph of the Clarke subdifferential is closed in topology;
- (c)
for all , one has
We now move to the properties of the subdifferential in the sense of convex analysis.
Definition 5. The subdifferential of a convex function is the multivalued operator defined by The following result represents an important property of the subdifferential of a convex function.
Proposition 4. Assume that is a convex lower semicontinuous function. Then the subdifferential operator is maximal monotone and .
A relevant example of function defined on
H with values on
is the indicator function defined by
where
. It is well known that if the subset
K is nonempty closed and convex, then the indicator function
is proper, convex and lower semicontinuous. Moreover,
. In the rest of this paper we shall use notation
for the subdifferential of the convex function
. Moreover, using Propositions 4 and 2 we deduce the following result.
Proposition 5. Assume that is a proper convex lower semicontinuous function and K is a nonempty closed convex subset of H. Then the operator is maximal monotone and, moreover, .
We now recall the following result proved in [
21].
Proposition 6. Let K be a nonempty closed convex subset of H, , a nonempty closed convex bounded subset of H and a proper convex lower semicontinuous function. Assume that for each there exists such that Then, there exists such that We end this section by recalling the notion of history-dependent operator. To this end, throughout this paper, for a normed space
we use the notation
for the space of continuous functions on
with values in
W. Recall that
can be organized in a canonical way as a complete metric space. The convergence of a sequence
to an element
v, in the space
, can be described as follows:
The next definition introduces two important classes of operators defined on spaces of continuous functions.
Definition 6. Let and be two normed spaces. An operator is said to be almost history-dependent if for any nonempty compact set there exist and such thatfor all and all . If, in particular, for any nonempty compact set , then is said to be a history-dependent operator. History-dependent and almost history-dependent operators arise in Functional Analysis, Solid Mechanics and Contact Mechanics, as well. General properties, examples and mechanical interpretations can be found in [
5]. In particular, the following fixed point property was proved in [
5], (p. 41).
Theorem 3. Let W be a Banach space and let be an almost history-dependent operator. Then Λ has a unique fixed point, i.e., there exists a unique element such that .
We shall use Theorem 3 in
Section 2 in the study of the history-dependent variational hemivariational inequality (
3).
3. A Parametric Variational–Hemivariational Inequality
In this section, in addition to the Hilbert space H, we assume that Y and Z are normed spaces endowed with the norms and , respectively. We also denote by the product of the spaces Y and Z. A typical point of X will be denoted by where and .
We now consider the following elliptic variational–hemivariational inequality: given the parameter
, find
u such that
In the study of this problem we consider the following assumptions.
Next, for any
we use the notation
and
for the functions defined by
Moreover, we introduce the operator
given by
and we recall that Propositions 3(a) and 5 guarantee that
. In addition, we have the following comment.
Remark 1. Assumption (10)(b)
implies thatfor any . Then, using Lemma 7 of [5] we deduce that the Clarke subdifferential of the function satisfies the relaxed monotonicity condition (4) with constant . We now state and prove various results related to the parametric variational-hemivariational inequality (
6).
Proposition 7. Assume (7)–(10) and let , , . Then u satisfies the inequalityif and only if Proof. Assume that
u satisfies the inequality (
16). We deduce from Proposition 3(c) that for each
there exists
such that
and, therefore,
Moreover, from Proposition 3(a), we obtain that the set
is a nonempty closed convex weakly compact subset of
H which implies that it is bounded, too. Hence, using Proposition 6 with
we see that there exists
which does not depend on
v, such that
Next, by the definition of the subdifferential of convex functions and inclusion
we have
This implies that
which shows that (
17) holds.
Conversely, assume that (
17) holds. Then, the definition (
15) of the operator
yields
Therefore, there exist
and
such that
Moreover, the definitions of the Clarke subdifferential and the subdifferential of a convex function imply that
,
for all
. Combining these inequalities with equality (
18) we deduce that
which shows that (
16) holds and concludes the proof. □
We now focus on the main property of the operator .
Proposition 8. Assume (7)–(11). Then, for any the operator is maximal monotone and, moreover, .
Proof. Let be fixed. The proof is structured in several steps, as follows.
Step (i) The operator is monotone. Indeed, assume that
,
. Then, there exist
and
such that
and, therefore,
We now use the inequalities
guaranteed by assumption (
8) and Remark 1, respectively, to see that
Therefore, assumption (
11) implies that the multivalued operator
is monotone.
Step (ii) The operator is maximal monotone. We start by proving that the mapping
has a closed graph in
. To this end let
and assume that
in
,
in
H as
and
for each
. Then,
and, since
, it is obvious to see that
Therefore, using (
19), assumption (
8) and the convergence
in
H, we deduce that
We now use the closedness of the graph of
in the product space
to see that
i.e.,
.
We conclude from above that the mapping has a closed graph in . Moreover, we use Proposition 3(a) to see that for any the set is a nonempty convex and weakly closed subset in H. The maximality of the monotone operator is now a consequence of Proposition 1.
Step (iii) The operator is maximal monotone. Indeed, Step (ii) and Proposition 3(a) show that the operator
is maximal monotone and
. Moreover, using (
7), (
9) and Proposition 5 we deduce that the operator
is maximal monotone and
. This implies that
. We now use Proposition 2 in order to deduce that the operator
is maximal monotone. Now, since (
15) shows that
, it follows that
is a maximal monotone operator. Moreover,
, which concludes the proof. □
Proposition 8 guarantees that, under assumptions (
7)–(
12), the operator
given by
is maximal monotone, too. Moreover,
. Therefore, for any
we are in a position to define its resolvent, denoted in what follows by
. We use (
5) to see that
and
for each
,
and
. We proceed with the following result.
Proposition 9. Assume (7)–(12), let , , , and let . Thenif and only if Proof. We use (
21), (
20) and (
15) to see that the following equivalences hold:
Proposition 9 is now a direct consequence of Proposition 7, used with the choice . □
We now take and obtain the following consequence of Proposition 9.
Corollary 1. Assume (7)–(12), let , and . Then u is a solution of the variational–hemivariational inequality (6) if and only if u is a fixed point of the operator , i.e., .
The following result represents a Lipschitz continuity result for the resolvent operator .
Proposition 10. Assume (7)–(11), let , , , and, for , let . Then, the following inequality holds: Proof. We use Proposition 9 to see that
Then, we take
in (
23),
in (
24) and add the resulting inequalities to find that
Next, we use the strong monotonicity of the operator
A, notation (
13) and (
14) and assumptions (
9)(b), (
10)(b) on the functions
and
j, respectively. In this way we deduce that
This inequality implies the bound (
22), which concludes the proof. □
We now consider the following additional assumptions.
We end this section with the following result concerning the parametric variational–hemivariational inequality (
6).
Theorem 4. Assume (7)–(10), (12). Then, the following statements hold.
(a) Under assumptions (25) and (26) the variational–hemivariational inequality (6) has at least one solution.
(b) Under assumption (27) the variational–hemivariational inequality (1) has a unique solution which depends Lipschitz continuously on f.
Proof. Let
,
,
and let
,
. Note that if (
26) or (
27) hold, than (
11) holds, too. Therefore, Proposition 10 implies that
(a) Assume that (
25), and (
26) hold. Then, using inequality (
28) we deduce that the operator
is non expansive. Therefore, we are in a position to use Theorem 1 to see that
has at least one fixed point. We now use Corollary 1 to deduce the solvability of the variational–hemivariational inequality (
1).
(b) Assume now (
27). Then, it follows from inequality (
28) that that the operator
is a contraction. Therefore, using the Banach fixed point principle we deduce that
has a unique fixed point. The unique solvability of the variational–hemivariational inequality (
1) is, again, a direct consequence of Corollary 1.
Finally, let
,
and let
,
denote the solution of inequality (
1) for
,
, respectively. Then
,
and, using Proposition 10, we deduce that
We now combine inequality (
29) with the smallness condition (
27) to deduce that the operator
is Lipschitz continuous, which concludes the proof. □
4. An First Elliptic Variational–Hemivariational Inequaliy
In this section, we study the solvability of the elliptic variational–hemivariational inequality (
1). To this end, in addition to assumptions (
7), (
8) and (
12), we consider the following assumptions on the functions
and
g:
Our main result in this section is the following
Theorem 5. Assume (7), (8), (12), (30) and (31). Then:
(a) Under assumptions (25) and (32) the variational–hemivariational inequality (1) has at least one solution.
(b) Under assumption (33) the variational–hemivariational inequality (1) has a unique solution which depends Lipschitz continuously on f.
Proof. Let
Y and
Z be arbitrary normed spaces and let
. For any
let
and
be the functions defined by
First, we see that the functions
and
j satisfy assumption (
9) and (
10) with
,
and
. Moreover, (
13), (
14) show that
and, therefore, with the notation above, inequality (
1) can be written in the equivalent form (
6). Based on this remark, we use in what follows the results in
Section 3.
(a) Assume that (
25) and (
32) hold and note that this implies that (
26) hold, too. Theorem 5(a) is now a direct consequence of Theorem 4(a).
(b) Assume now that (
33) and note that this implies that (
27) hold, too. Theorem 5(b) is now a direct consequence of Theorem 4(b). □
The proof of Theorem 5 shows that the solution of inequality (
1) is a fixed point for the resolvent operator
, defined for any
and
. This suggests us to consider several iterative methods in the solution of this inequality. Details on iterative methods for nonexpansive and contractions operators can be found in [
20,
22]. Here, we restrict ourselves to present only two examples. Note that, since in the particular case of Theorem 5 the operator
does not depend on
w, we shall denote it in what follows by
.
Example 1. (Picard iterations.) Under assumptions of Theorem 5
, let be arbitrary given and define the sequence by equality Using Proposition 9 and notation (34), (35), it is easy to see that equality (36) can be written, equivalently, as follows: We conclude that at each step of this iterative scheme we have to solve an elliptic variational-hemivariational inequality. Now, since the operator is a contraction, the sequence converges strongly in H to the fixed point of this operator and, therefore, to the solution u of inequality (1).
Example 2. (Krasnoselski iterations) Under assumptions of Theorem 5
, let , be arbitrary given and define the sequence by equality Then, since the operator is nonexpansive, it is well known that the sequence converge weakly in H to a fixed point of this operator. A proof of this result can be found in ([20], p. 61). Therefore, in H, where u is a solution of inequality (1). Consider now the particular case when . Then equality (37) becomesand, therefore, using notation , Proposition 9 implies that equality (38) leads to the following iterative scheme: given , the sequence is determined, recursively, by solving the systemfor all . 5. An Second Elliptic Variational–Hemivariational Inequaliy
In this section, we use the results in
Section 3 in the study the solvability of the elliptic variatioanal-hemivariational inequality (
2). To this end, in addition to assumptions (
7), (
8) and (
12) we consider the following assumptions.
Note that assumptions (
39) and (
40) represent a particular case of assumptions (
9) and (
10), respectively, obtained when
. Therefore, the preliminary results in
Section 3 can be used in the study of inequality (
2).
Our main result in this section is the following.
Theorem 6. Assume (7), (8), (12), (39), (40). Then:
(a) Under assumptions (25) and (41) the variational–hemivariational inequality (2) has at least one solution.
(b) Under assumption (42) the variational–hemivariational inequality (2) has a unique solution which depends Lipschitz continuously on f.
Proof. (a) Assume that (
25) and (
41) hold and recall the equivalence (
21). Let
and consider the operator
defined by
where
, for all
. Let
and let
,
. Then, using Proposition 10 with
,
,
for
, we deduce that
We now use assumption (
41) and inequality (
43) to see that the operator
is nonexpansive. Therefore, we are in a position to use Theorem 1 to see that
has at least one fixed point. We now use Corollary 1 to deduce the solvability of the variational–hemivariational inequality (
2).
(b) Assume now that (
42) hold. Then, inequality (
43) shows that the operator
is a contraction and, using the Banach principle we deduce that
has a unique fixed point. The unique solvability of the variational–hemivariational inequality (
2) is, again, a direct consequence of Corollary 1. Finally, the Lipschitz continuity of the mapping
follows from arguments similar to those used in the proof of Theorem 4(b). □
The proof of Theorem 6 shows that the solution of inequality (
2) is a fixed point for the resolvent operator
, for any
. This allows us to consider several iterative methods in the solution of this inequality. In order to avoid repetitions we restrct here to the Picard iterations.
Example 3. Under assumptions of Theorem 6
, let be arbitrary given and define the sequence by equalitywhere, recall, represents a short hand notation for the pair . Then, using Proposition 9, it is easy to see that equality (44) can be written, equivalently, as follows: Now, since the operator is a contraction, it follows that the sequence converge strongly in H to the solution u of inequality (2).