Tykhonov well-posedness of ﬁxed point problems in contact mechanics

We consider a ﬁxed point problem S u = u where S : C ( R + ; X ) → C ( R + ; X ) is an almost history-dependent operator. First, we recall the unique solvability of the problem. Then, we introduce the concept of Tykhonov triple, provide several relevant examples, and prove the corresponding well-posedness results for the considered ﬁxed point problem. This allows us to deduce various consequences which illustrate the stability of the solution with respect to perturbations of the operator S . Our results provide mathematical tools in the analysis of a large number of history-dependent problems which arise in solid and contact mechanics. To give some examples, we consider two mathematical models which describe the equilibrium of a viscoelastic body in frictionless contact with a foundation. We state the mechanical problems, list the assumptions on the data, and derive their associated ﬁxed point formulation. Then, we illustrate the use of our abstract results in order to deduce the continuous dependence of the solution with respect to the data and parameters. We also provide the corresponding mechanical interpretations.


Introduction
The concept of well-posedness in the sense of Tykhonov was introduced in [26] for a minimization problem. It is based on two main ingredients: the existence and uniqueness of solution to the problem and the convergence to it of any approximating sequence. The well-posedness in the sense of Tykhonov (well-posedness, for short) has been extended for various optimization problems, see [2,3,9,12] for instance. In addition, it was introduced in the study of variational and hemivariational inequalities in [13,14] and [6], respectively. References in the field include [5,8,22,27,28], where general results concerning the well-posedness of various classes of inequalities can be found. A Tykhonov well-posedness result for a fixed point problem can be found in [18]. A general framework which unifies the view on well-posedness of abstract problems in metric spaces was recently considered in [23,25]. The approach used in [25] was based on the concept of Tykhonov triple T = (I, , C) where I is a set of parameters, represents a family approximating sets, and C is a set of sequences of elements of I. History-dependent and almost history-dependent operators represent a special class of nonlinear operators defined on a space of continuous functions. They arise both in functional analysis, solid mechanics, and contact mechanics. In contact mechanics such operators could be involved either in the constitutive law, which describes the material behavior, or in the interface boundary conditions. Such kind of operators have a number of relevant properties, including fixed point properties and, for this reason, they have been intensively studied in the recent years. A reference in the field is the book [21]. Existence and uniqueness results on variational and hemivariational inequalities with history-dependent operators can be found in [20,21]. General results on their numerical analysis could be found in [7].
Processes of contact between a deformable solid and a foundation are ubiquitous, and they can be found in many industrial settings, in transportation, in various scientific experimental settings, and in everyday life. This is the reason for the very large engineering literature dedicated to the modeling, numerical approximations, and computer simulations of such processes. On the other hand, the mathematical theory of contact mechanics has expanded substantially in the last decades and is quickly maturing. A sample of references are, e.g., the books [1, 4, 10, 15-17, 20, 21]. The concepts of Tykhonov triple and Tykhonov well-posedness have been used in [11,24], in the study of quasistatic contact problems with viscoelastic or viscoplastic materials.
In this current paper we deal with the well-posedness of a class of fixed problems with emphasis to the study of history-dependent models of contact. To this end we need the notation that we introduce in what follows. First, we use N and R + to represent the set of positive integers and the set of real nonnegative numbers, i.e., N = {1, 2, 3, . . .} and R + = [0, +∞). Unless stated otherwise, X will be a real Banach space equipped with the norm · X . We also denote by C(R + ; X) the space of continuous functions on R + with values in X and use notation 0 X and 0 X for the zero element of the spaces X and C(R + ; X), respectively. Recall that C(R + ; X) can be organized in a canonical way as a Fréchet space, i.e., as a complete metric space in which the corresponding topology is induced by a countable family of seminorms. The convergence of a sequence {v n } to an element v can be described as follows: We also mention that, unless stated otherwise, all the limits below are considered as n → ∞, even if we do not mention it explicitly.
With these notation, the fixed point problem we study in this paper is stated as follows.
Problem P Given an operator S : C(R + ; X) → C(R + ; X), find a function u ∈ C(R + ; X) such that Note that here and below in this paper we use the shorthand notation Su(t) to represent the value of the function Su at the point t, i.e., Su(t) = (Su)(t). Moreover, for an element u 0 ∈ X, we still write u 0 for the constant function t → u 0 , for all t ∈ R + . Therefore, notation Su 0 used below defines an element of C(R + ; X).
Our aim in this paper is twofold. The first one is to deduce stability results for Problem P with respect to perturbations of the operator S. To this end we use a strategy based on the well-posedness of the fixed point problem P with respect to various Tykhonov triples. This represents, at the best of our knowledge, the first trait of novelty of our work. The second aim is to illustrate how these abstract results can be applied in the study of mathematical models of contact. To this end we consider two models of contact for which, in contrast with various references in the literature, we provide a variational formulation governed by fixed point problems. Using such formulation allows us to deduce some convergence results, which represents the second trait of novelty of this work.
The rest of the paper is structured as follows. In Sect. 2 we introduce the concept of well-posedness for the fixed point problem (2), then we prove its well-posedness with respect to four relevant Tykhonov triples. The proofs are based on various estimates, the properties of almost history-dependent operators, and the Gronwall argument. Next, we use these well-posedness results in Sect. 3. There, we prove various convergence results which show the stability of the solution with respect to perturbations of the operator S. In Sect. 4 we present our first example of contact problem. It concerns a model of quasistatic frictionless contact between a viscoelastic body and a deformable foundation. In Sect. 5, we present our second example which, in contrast, models the quasistatic contact between a viscoelastic body with a rigid foundation covered by a layer of elastic material, say asperities. We end this paper with Sect. 6 in which we present some concluding remarks.

Tykhonov well-posedness
We start with an existence and uniqueness result in the study of Problem P. To this end, we recall the following definition. Definition 2.1 An operator S : C(R + ; X) → C(R + ; X) is called an almost historydependent (a.h.d.) operator if for any m ∈ N there exist l m ∈ [0, 1) and L m > 0 such that An operator S : C(R + ; X) → C(R + ; X) is called a history-dependent (h.d.) operator if inequality (3) holds with l m = 0 for all m ∈ N.
In order to avoid any confusion, we underline that in Definition 2.1 and everywhere in this paper the superscript m does not represent a power.
The unique solvability of Problem P is guaranteed by the following existence and uniqueness result.

Theorem 2.2
Assume that S is an almost history-dependent operator. Then Problem P has a unique solution u ∈ C(R + ; X).
A proof of Theorem 2.2 can be found in [19] and [21, p. 42], based on the Banach fixed point principle. There, various examples of a.h.d. and h.d. operators are provided and general properties of these classes of operators are proved.
We now introduce the concepts of Tykhonov triple and well-posedness of the fixed point problem (2) with respect to a given Tykhonov triple. These concepts have been introduced in [25] in the framework of abstract problems in metric spaces. Here we specify them to the particular setting of Problem P and, to this end, we denote by X the set of nonempty subsets of the space C(R + ; X). Moreover, for any set J, we use the notation R(J) for the set of sequences with elements in J.

Definition 2.3
A Tykhonov triple is a mathematical object of the form T = (I, , C) where I is a given nonempty set, : I → X is a multivalued mapping, and C is a nonempty subset of the set R(I).

Definition 2.4
Given a Tykhonov triple T = (I, , C), a sequence {u n } ⊂ C(R + ; X) is called a T -approximating sequence if there exists a sequence {ω n } ∈ C such that u n ∈ (ω n ) for each n ∈ N. Definition 2.5 Given a Tykhonov triple T = (I, , C), Problem P is said to be well-posed with T (or T -well-posed, for short) if it has a unique solution and every T -approximating sequence converges in the space C(R + ; X) to its solution.
We now complete the definition above with the following comments. First, we refer to I as the set of parameters. In the examples we consider below in this paper we shall take either I = R + or I = R(R + ). In the case when I = R + we use the notation θ for an element of I and θ = {θ n } n (or θ = {θ n }, for short) for an element of R(I). In the case when I = R(R + ) we shall use the notation θ = {θ m } m for an element of I and θ = {θ n } n (or θ = {θ n }) with θ n = {θ m n } m for an element of R(I). Second, refer to the family of sets { (ω)} ω∈I as the family of approximating sets and, moreover, we say that C defines the criterion of convergence. Note that approximating sequences always exist since, by assumption, C = ∅ and for any sequence {ω n } ∈ C and any n ∈ N, the set (ω n ) is not empty.
Next, we underline that the concept of well-posedness above is not an intrinsic concept for Problem P. Indeed, it is defined through the convergence of the approximating sequences which, in turn, depends on the Tykhonov triple T we consider. For this reason, we use the terminology "T -approximating sequence" and "well-posedness with T " or "T -well-posedness", for short.
Assume now that S is an almost history-dependent operator and, therefore, Problem P has a unique solution u ∈ C(R + ; X). Denote Moreover, for any Tykhonov triple T , denote by R T the set of T -approximating sequences, i.e., Then it is easy to see that the well-posedness of a Problem P with respect to the Tykhonov triple T is equivalent with the inclusion R T ⊂ R P , i.e., We now introduce the following definition.

Definition 2.6
Given two Tykhonov triples T = (I, , C) and T = (I , , C ), we say that: (a) T and T are equivalent (and we write It is easy to see that "∼" represents an equivalence relation on the set of Tykhonov triples while "≤" defines a relation of order on the same set. Moreover, using (6) we deduce that the following statements hold.
If T ≤ T and P is T -well-posed, then P is T -well-posed, too.
We now construct four relevant triples in the study of Problem P.
Example 2.7 Assume that S is an a.h.d. operator and take T 1 = (I 1 , 1 , C 1 ) where for all θ ∈ I 1 , Note that for each θ ∈ I 1 the fixed point u obtained in Theorem 2.2 belongs to 1 (θ ), which shows that 1 (θ ) = ∅. Therefore, according to Definition 2.3, T 1 is a Tykhonov triple.
The following result is useful in order to compare the Tykhonov triples T 3 and T 4 .

Lemma 2.11
Assume that S is an a.h.d. operator, and let {u n } be a T 4 -approximating sequence. Then, for each m ∈ N, there exists Z m > 0 such that Proof Using Definition 2.4, (19), and (20), we deduce that there exists a sequence {θ n } n with θ n = {θ m n } m ∈ S(R + ) for any n ∈ N such that θ m n → 0 as n → ∞, ∀m ∈ N and, for any m, n ∈ N, the following inequality holds: Let m ∈ N, t ∈ [0, m] and n ∈ N. We write then we use inequalities (23) and (3) to find that Next, (22) and inequality l m < 1 imply that for n large enough we can assume that θ m n ≤ 1-l m 2 and, with notation inequality (24) yields We now use the Gronwall argument to find that 1-l m t and, using convergence (22) combined with inequality t ≤ m, we conclude the proof.
We now remark that, obviously, In addition, using Lemma 2.11 it is easy to see that R T 4 ⊂ R T 3 and, therefore, R T 3 = R T 4 . We now use Definition 2.6 to see that We are now in a position to introduce our main result in this section.

Theorem 2.12
Assume that S is an almost history-dependent operator. Then Problem P is well-posed with respect to the Tykhonov triples T 1 , T 2 , T 3 , and T 4 in Examples 2.7, 2.8, 2.9, and 2.10, respectively.
Proof First, we recall that the existence of a unique solution to problem P, needed for the well-posedness of Problem P with any Tykhonov triple, follows from Theorem 2.2.
Consider now the Tykhonov triple T 3 in Example 2.9, and let {u n } be a T 3 -approximating sequence. Then, using Definition 2.4, we deduce that there exists a sequence {θ n } n with θ n = {θ m n } m ∈ R(R + ) for all n ∈ N such that (22) holds and, moreover, Let n ∈ N, m ∈ N and let t ∈ [0, m]. We use (2) and write then we use inequalities (28) and (3) to find that Next, arguments similar to those used in the proof of Lemma 2.11 show that, for n large enough, inequality (29) yields and, therefore, We now combine inequality (31) and convergence (22) to see that and, since m ∈ N is arbitrary, using (1) we deduce that Convergence (32) combined with Definition 2.5 implies that Problem P is well-posed with respect to the Tykhonov triple T 3 in Example 2.9. The well-posedness of Problem P with the Tykhonov triples T 1 , T 2 , and T 4 is now a direct consequence of relations (27), (7), and (8).

Stability results
In this section we study the stability of the solution of Problem P with respect to perturbations of the operator S. To this end, for each n ∈ N, we consider the following fixed point problem.
Problem P n Given an operator S n : C(R + ; X) → C(R + ; X), find a function u n ∈ C(R + ; X) such that S n u n (t) = u n (t) for all t ∈ R + .
We assume that for each n ∈ N the operator S n is an almost history-dependent operator, i.e., For each n, m ∈ N, there exist l m n ∈ [0, 1) and L m n > 0 such that Then, using Theorem 2.2, it follows that Problem P n has a unique solution u n ∈ C(R + ; X) for each n ∈ N. Our aim in what follows is to state various conditions which link the operators S n and S and guarantee that the solution of Problem P n converges to the solution of Problem P, that is, To this end we use the following two step strategy.
(i) First, we impose conditions which guarantees that the sequence {u n } is a T i -approximating sequence for some i ∈ {1, 2, 3}.
(ii) Second, we use the T i -well-posedness of Problem P, obtained in Theorem 2.12, to deduce that convergence (34) holds.
The conditions we consider are the following.
We have the following convergence result. Proof (a) Let n ∈ N and t ∈ R + . We use (33) and (36)(a) to write Su n (t)u n (t) X ≤ Su n (t) -S n u n (t) X ≤ θ n . This inequality combined with assumption (36)(b) and definition (9)- (11) shows that the sequence {u n } is a T 1 -approximating sequence. Convergence (34) is now a consequence of the T 1 -well-posedness of Problem P, guaranteed by Theorem 2.12.
(b) Let n ∈ N and t ∈ R + . We use (33) and (37)(a) to write Su n (t)u n (t) X ≤ Su n (t) -S n u n (t) X ≤ θ n u n (t) X + 1 . This inequality combined with assumption (37)(b) and definition (12)- (14) show that in this case the sequence {u n } is a T 2 -approximating sequence. Convergence (34) is now a consequence of the T 2 -well-posedness of Problem P, guaranteed by Theorem 2.12.
(c) We first prove that for each m ∈ N there exists U m > 0 such that Indeed, let m ∈ N, n ∈ N and let t ∈ [0, m]. We use (33) and write then we use inequalities (38)(a) and (3) to find that Therefore, combining this inequality with notation (25), we obtain Next, (38)(b) and the inequality l m < 1 imply that for n large enough we can assume that θ m n ≤ 1-l m 2 < 1 and, therefore, (40) yields We now use the Gronwall argument to find that and, using inequality t ≤ m, we conclude that (39) holds with 1-l m m .

We now use (33) and (38)(a) to write
Su n (t)u n (t) X = Su n (t) -S n u n (t) X ≤ θ n 1 + u n (t) X + for each m, n ∈ N. Then inequality (42) implies that u n ∈ 3 (ω n ) where, recall, 3 (ω) is the set defined by (16) for each ω = {ω m } m ∈ R(R + ). On the other hand, assumptions (38)(b) and definition (43) imply that ω m n → 0 as n → ∞ for each m ∈ N. This shows that ω n = {ω m n } m ∈ C 3 where C 3 represents the set given by (17). It follows from the above that the sequence {u n } is a T 3 -approximating sequence for Problem P. We now use Theorem 2.12 and Definition 2.5 to deduce convergence (34), which concludes the proof.
We now remark that assumption (37) does not guarantee that the sequence {u n } is a T 1approximating sequence. Moreover, assumption (38) does not guarantee that the sequence {u n } is a T 2 -approximating sequence. An evidence of these statements is provided by the two examples below.
Example 3.2 Let X = R and let S n , S : C(R + ; X) → C(R + ; X) be the operators defined by for all u ∈ C(R + ; X), t ∈ R + , and n ∈ N. Then, it is easy to see that S n and S are almost history-dependent operators which satisfy condition (37). Moreover, with this choice, the solution of Problem P n is given by and the solution of Problem P is Arguing by contradiction, assume that the sequence {u n } is a T 1 -approximation sequence. Then there exists {θ n } ⊂ R(R + ) such that θ n → 0 and |Su nu n | ≤ θ n ∀n ∈ N, t ∈ R + . (47) Using now (44)-(47) we deduce that 1 n(n + 1) e n+1 n t + 1 n + 1 ≤ θ n ∀n ∈ N, t ∈ R + .
We now take t = n 2 in the previous inequality, then we pass to the limit as n → ∞ and arrive at a contradiction. We conclude from here that {u n } is not a T 1 -approximation sequence. Nevertheless, using (45) and (46) it is easy to see that convergence (34) holds. for all u ∈ C(R + ; X), t ∈ R + , and n ∈ N. Then it is easy to see that S n and S are historydependent operators which satisfy condition (38). Moreover, with this choice, the solution of Problem P n is given by u n (t) = e -t n ∀n ∈ N, t ∈ R + (49) and the solution of Problem P is Arguing by contradiction, assume that the sequence {u n } is a T 2 -approximation sequence. Then there exists {θ n } ⊂ R(R + ) such that θ n → 0 and |Su nu n | ≤ θ n u n (t) + 1 ∀n ∈ N, t ∈ R + . (51) Using now (48)-(51) we deduce that 1e -t n ≤ θ n e -t n + 1 ∀n ∈ N, t ∈ R + .
We take t = n in the previous inequality, then we pass to the limit as n → ∞ and arrive at a contradiction. We conclude from here that {u n } is not a T 2 approximation sequence. Nevertheless, convergence (34) holds.
It follows from Example 3.2 that, under condition (37), the T 1 -well-posedness of Problem P cannot be used to provide convergence (34) in the framework of the two step strategy (35). Similarly, Example 3.3 shows that, under condition (38), the T 2 -well-posedness of Problem P cannot be used to provide convergence (34) in the framework of the same strategy. We conclude from the above that the choice of the Tykhonov triple plays a crucial role in proving convergence result for Problem P.

A contact problem with normal compliance
In this section we apply Theorem 3.1 in the study of a quasistatic frictionless contact problem with normal compliance. The physical setting is the following. A viscoelastic body occupies a bounded domain ⊂ R d (d = 1, 2, 3) with a Lipschitz continuous boundary , divided into three measurable disjoint parts 1 , 2 , and 3 such that meas( 1 ) > 0. The body is fixed on 1 , is acted upon by a given surface traction on 2 , and could arrive in contact with a deformable foundation on 3 . Moreover, there is a gap between the body and the foundation, the mechanical process is quasistatic, the contact is frictionless, and the time interval of interest is R + = [0, +∞). We denote by ν the unit outward normal to and by S d the space of second order symmetric tensors on R d . Then the classical formulation of the contact problem we study in this section is the following.
Problem Q Find a displacement field u : × R + → R d and a stress field σ : × R + → S d such that, for any t ∈ R + , the following equalities hold: Div We now provide a short description of the equations and boundary conditions in Problem Q and send the reader to [20,21] for more details and mechanical interpretations. First, equation (52) represents the viscoelastic constitutive law of the material in which ε(u) denotes the linearized strain tensor, E is a fourth order elasticity tensor, and F is the relaxation tensor. Equation (53) represents the equilibrium equation in which f 0 denotes the density of body forces, (54) is the displacement boundary condition, and (55) is the traction boundary condition in which f 2 denotes the density of surface tractions. Condition (56) is the contact condition with normal compliance in which p is a given prescribed function, g represents the gap function, and k is a positive stiffness coefficient. Moreover, u ν and σ ν denote the normal components of the displacement and the stress field, respectively. Finally, condition (56) is the frictionless condition in which σ τ represents the tangential shear.
In the variational analysis of Problem Q we denote by "·", · , and 0 the inner product, the Euclidean norm, and the zero element of the spaces R d and S d , respectively. We use the standard notation for the Sobolev and Lebesgue spaces associated to and and, for an element v ∈ H 1 ( ) d , we write v for the trace γ v ∈ L 2 ( ) d of v on . Moreover, we denote by v ν and v τ the normal and tangential components of v on the boundary, given by v ν = v · ν and v τ = vv ν ν, respectively. We also use the following spaces: The spaces V and Q are real Hilbert spaces endowed with the canonical inner products where, for any v ∈ V , ε(v) represents the symmetric part of the gradient of V . The associated norms on these spaces are denoted by · V and · Q , respectively. The space Q ∞ is a real Banach space with the norm Moreover, the following inequality holds: We now list the assumption on the data of the contact problem Q. The elasticity tensor E is symmetric and positively definite, i.e., it satisfies the conditions The relaxation tensor F and the normal compliance function p satisfy the following conditions: for all r 1 , r 2 ∈ R, a.e. x ∈ 3 .
(c) p(·, r) is measurable on 3 , for all r ∈ R. (61) We also assume that the densities of the body forces and surface tractions, the stiffness coefficient and the gap function are such that f 0 ∈ C I; L 2 ( ) d . (62) We now use assumption (59) to endow the space V with the inner product and the associated norm · E . Then, using inequalities (59)(b) and (58), it follows that which shows that the norms · V and · E are equivalent norms on the space V . Therefore, V is a Hilbert space endowed with the inner product (·, ·) E , too. Moreover, it follows from the Sobolev trace argument and the continuity of the embedding V ⊂ L 2 ( ) d that there exist some constants c 0 > 0, d 2 > 0, and e 0 > 0 such that for all v ∈ V . Finally, we assume that the following smallness condition holds: where c 0 and L p are the positive constants in (68) and (61), respectively. We now turn to construct a fixed point weak formulation of Problem Q. Let v ∈ V and t ∈ R + . We use standard arguments based on the Green formula to see that if (u, σ ) is a smooth solution of Problem Q then We combine this equality with the constitutive law (52) to see that Next, we use the Riesz representation theorem to define the operator S : C(R + ; V ) → C(R + ; V ) by the equality for all u ∈ C(R + ; V ), t ∈ R + , and v ∈ V . We now combine equalities (66), (71), and (72) to deduce the following fixed point weak formulation of Problem Q in terms of displacements.
Problem Q V Find a displacement field u ∈ C(R + ; V ) such that Su(t) = u(t) for any t ∈ R + .
The unique solvability of Problem Q V is provided by the following existence and uniqueness result. Proof Let u, v ∈ C(R + ; X), m ∈ N, t ∈ [0, m], and w ∈ W . We use the definition (72) to see that and, therefore, assumption (61)(b) and inequality (58) yield We now use the trace inequality (68) to deduce that Finally, we use (67), then we take w = Su(t) -Su(t) in the resulting inequality to find that Inequality (73)  We now prove a continuous dependence result of the solution with respect to the data. To this end we consider the sequences {F n }, {f 0n }, {f 2n }, {k n } such that, for each n ∈ N, the following conditions hold: With these data we define the operator S n : C(R + ; V ) → C(R + ; V ) by the equality for all u ∈ C(R + ; V ), t ∈ R + , and v ∈ V . Then we consider the following variational problem.
Problem Q V n Find a displacement field u n ∈ C(R + ; V ) such that S n u n (t) = u n (t) for any t ∈ R + .
Using Theorem 4.1 it follows that Problem Q n has a unique solution for each n ∈ N. Consider now the following assumptions: f 0n → f 0 in C I; L 2 ( ) d as n → ∞, We have the following convergence result.
It follows now from (90) and (91) that condition (38) is satisfied. Theorem 4.2 is now a direct consequence of Theorem 3.1(c).
We end this section with the following comments. (1) We refer to the solution u of Problem Q V as the weak solution for the contact Problem Q. We conclude from the above that Theorem 4.1 provides the unique weak solvability of this contact problem and Theorem 4.2 shows the continuous dependence of the weak solution with respect to the data.
(2) In addition to the mathematical interest in convergence (84), it is important from the mechanical point of view since it shows that small perturbation on the relaxation tensor F , the density of body forces f 0 , the density of traction forces f 2 , and the stiffness coefficient k lead to small perturbation of the weak solution of the contact problem Q.

A contact problem with unilateral constraints
For the contact problem we consider in this section, we keep the physical setting and the notation in Sect. 4. Nevertheless, we now assume that the foundation is made by a rigid obstacle covered by a layer of deformable material, say asperities. Moreover, there is no gap between the body and the foundation, and we model the material's behavior with a different viscoelastic constitutive law. The classical formulation of the problem is the following.
Problem M Find a displacement field u : × R + → R d and a stress field σ : × R + → S d such thaṫ σ (t) = Eε u(t) + β(σ (t) -G ε u(t) in , Div σ (t) + f 0 (t) = 0 in , u ν (t) ≤ g, σ ν (t) + p(u ν (t)) ≤ 0, (u ν (t)g)(σ ν (t) + p(u ν (t)) = 0 on 3 , σ τ (t) = 0 on 3 (97) for any t ∈ R + and, moreover, The equations and boundary conditions in the statement of Problem M have a similar meaning as those used in the statement of Problem Q. The novelty arises in the fact that now we use the constitutive law (92), in which β represents a viscosity coefficient, G is a nonlinear relaxation operator, and the dot above denotes the derivative with respect to the time. In addition, we use a contact condition with unilateral constraint (96). Here g represents the thickness of the deformable layer. Finally, note that the functions u 0 and σ 0 in (98) are the initial displacement and the initial stress, respectively. Details on similar viscoelastic and viscoplastic contact models can be found in [20,21].