In this paper, we take into consideration the mathematical analysis of time-dependent quasistatic processes involving the contact between a solid body and an extremely rigid structure, referred to as a foundation. It is assumed that the constitutive law is fractional long-memory viscoelastic. The contact is considered to be bilateral and is modeled around Tresca’s law. We establish the existence of the generalized solution’s result. The proof is supported by the surjectivity of the multivalued maximum monotone operator, Rothe’s semidiscretization method, and arguments for evolutionary variational inequality.

1. Introduction

It is well known that the empirical models containing the fractional derivative fit the experimental data more precisely than the model containing the integer derivative. The early application of the fractional derivative in viscoelasticity dates back to the beginning of this century. Germent and Baglet proposed the fractional Kelvin–Voigt model, and at the same time, Koeller obtained the fractional Maxwell model [1, 2]. The theoretical and numerical investigation has been done for fractional differential equations and inequalities in finite-dimensional spaces. On this subject, we cite, for instance [36]. The application of the variational framework, more precisely the variational inequalities, proves to be very useful in the study of engineering and mechanical problems. It was initiated by Panagiotopoulos; see monographs [7, 8], and several results have emerged; see [913]. In this paper, we are interested in studying the contact problem for a fractional viscoelastic law with long-term memory. The corresponding constitutive law iswhere , and represent the fractional viscoelasticity, elasticity operator, and relaxation tensor, and stand for the stress and the displacement fields, respectively, and is related to Caputo fractional derivative of order and .

The long-memory viscoelastic constitutive law has been the subject of numerous papers; see [1416]. The distinction of this study is the selection of contact with Tresca’s friction type, governed by a fractional constitutive law, and the use of Rothe’s approach and the maximum monotone operators’ surjectivity results to demonstrate the existence of a solution to a variational inequality. Actually, there are only a few contributions to the literature on the resolution of variational and hemivariational inequalities using Rothe’s approach; see, for example, [1821].

This article is structured according to the following outlines: in Section 2, we state the equilibrium equations modeling the process of viscoelastic material coming into contact with a foundation. We give some of the notations and the assumptions on data. Then, we obtain the variational formulation of the problem and present the main results concerning the existence of a weak solution. In Section 3, we first prove the existence of the discrete generalized solution to the discrete problem, in the second, we establish the a priori estimates, and in the last, we pass to the limit to obtain the existence of solutions to the variational problem.

2. Setting of the Problem

Before setting up the physical model, we give some definitions and preliminary results.

Definition 1. For any positive integer and , the Caputo fractional derivative and fractional integral of order of a given function , are respectively defined aswhere is the standard Gamma function.

Definition 2. Let be a real Banach space. An operator is said to be(a)est hemicontinuous if the real function is continuous from [0,1] in , for all ,.(b)bounded, if maps bounded sets of into bounded sets of .(c)is monotone if, .(d) is maximal monotone if, is monotone, and for all ,, such that then .(e)The generalized gradient (subdifferential) of a convex function at is a subset of the dual space given by(f)A sequence is weakly convergent to , noted , if for all , we havewhere represents the duality between and .

Theorem 1 (see [22]). Let is convex, lower semicontinuous, and . Then, is maximal monotone.

Theorem 2 (see [22]). Let the following hypotheses hold: the operator is maximal monotone.: the operator is monotone hemicontinuous and bounded.: is -coercive, meaning that a point and a number exist, that wayHence, is surjective.

2.1. Physical Setting

Suppose that in the reference configuration, the domain which is considered bounded with a smooth boundary . Three open, disconnected sections make up this boundary: , , and such that . Let time interval of interest, where .

The body is submitted to the action of body forces of density . It also submitted to mechanical constraints on the boundary. The body is supposed squeezed into , thus the displacement field vanishes over . On , there are the density of traction forces denoted . At some point, the body comes into contact over .

Let denote by the space of the second order of symmetric tensors and by and , respectively, the scalar product and the Euclidean norm in (resp in ). We adopt the summation convention over repeated indices. All indices , take values from 1 to .

Let , we shall use the notationwith and are, respectively, operators of deformation and divergence defined by

The space , , , and are Hilbert spaces endowed with the inner products given by

We use the notation for the trace of on . The normal and tangential components of on are represented by the symbols and , respectively, and are determined by the formulas , , , and . The trace of the regular stress field applied to is noted .

The physical model for the process is as follows:

Problem 1. Find a displacement field The boundary and initial conditions

For the sake of simplifying the notations in equations (10)–(16) and below, we do not show explicitly the dependence of all functions on, respectively, space and time variables and . The fractional viscoelastic constituency with memory law is given by (10), for , where is linear viscosity function, is linear elasticity, and the linear relaxation function. The expression (11) is the equilibrium equation, we suppose the process is quasistatic. Here, the boundary conditions on displacement and traction are represented by (12) and (13). Conditions (14) and (15) are the bilateral contact and the friction conditions, respectively, it is modeled by Tresca’s law with an imposed friction bound (hence and the nonseparation condition ). The equation given in (16) is the initial condition on displacement. To be able to give a variational formulation of the above problem, we give the following additional notations. We setequipped with the inner product for , . In view of , there exist, respectively, and , called Korn’s and Poincaré constants, such that

Let denote by the Bochner–Lebesgue function space. Since the space is reflexive, it remains that and its dual space are reflexive too. To study the problem and present our main result, we make the following assumptions: the viscosity operator , the elasticity one , and the relaxation tensor satisfy the conditions

The initial data and the density forces satisfy

3. Variational Formulation of the Problem

In this section, we establish the existence of a weak solution to Problem 1. First, we multiply the equilibrium equation by the test function, taking into account the initial and boundary conditions, and we obtain Problem 2 in the form of a variational inequality. Then, we reformulate the problem in operational form (Problem 3), which allows us to pass from a heavy scripture relating to tensorial products to a more convenient one. In order to apply the classical results of functional inclusions, we rewrite the last problem in the form of a nonlinear inclusion, thus giving Problem 4. Note that Problem 2, 3, and 4 are equivalent.

To obtain the variational formulation of Problem 1, we assume that the pair of functions and are sufficiently smooth. We recall that Green’s formula holds:

From the (11), the boundary conditions (12)–(15) and equality (23), we obtain for all


Let define , for all , expression (25) leads to the following variational formulation:

Problem 2. Find , such that

Let define the operators , , and by the expressionswhere is defined below. Note that since is convex, we have

The following problem can be established in light of all those considerations.

Problem 3. Find , such thatx

Let now suppose that be a solution to Problem 3 and put , for where is the Caputo fractional derivative. Note that when , then have a sense.

So, we may rewrite our problem as a nonlinear inclusion, driven by a fractional integral operator, and it is given by the problem below.

Problem 4. Find , such that

Using data from Problems 3 and 4, we reformulate the assumptions. :  is coercive, i.e., there exists a constant , such that : . there exists , such that, and . .

Let be fixed, , , and for all , be defined by , for . Consider the so-called Rothe problem, i.e., the discretized problem corresponding to Problem 4, given by as follows:

Problem 5. Find , such that , andwhere and for , and , are defined by

First, we shall show the existence of solution for a Problem 5.

Theorem 3. Suppose that conditions , , , and are verified. Hence, there exists, such that for all, there is at least one solution to Problem 5.

Proof. Given , we claim that there exist , such that (33) hold. From equality (34), we obtain elements . We have to show that there exists elements and such that inclusion (33) holds. Let denote byTo accomplish this, we will demonstrate that the multivalued operatorSuppositions , , and imply the linearity and boundedness of the operatorso, it is hemicontinuous, and satisfies the conditionThus, under the condition , the operatorFor the mapping , which it is, continuous and convex for all , then by Theorem 1, we obtain that is maximal monotone for all . The coercivity is so fulfilled if . Therefore, for all , where is given byWe may now use Theorem 2 to guarantee the surjectivity of the operator , for all . Thus, we conclude that there exists elements , satisfying the (33). This completes the proof of the theorem.
Note that is solution of (33), also means that there exists , such that

It is noteworthy that hereafter, we use to denote a generic positive constant independent of . The sequence of solutions to Rothe Problem 5 is estimated in the finding that follows:

Lemma 1. Let’s say that all four hypotheses , , , and , are fulfilled, then there existsindependent of , such that for all, the solutions of Problem 5 satisfywhere the constant here is independent of and .

Proof. We multiply equation (33) by taking and to getBy the expression (34) of and hypotheses and , we havewhereAccording to the definition of , we havewhere is Poincaré’s constant. This implies thatNow, we conclude that because of the coercivity of the operator and expressions (46), (47), and (50),