- Research
- Open access
- Published:
Dirichlet boundary value problem for differential equations involving dry friction
Boundary Value Problems volume 2015, Article number: 106 (2015)
Abstract
Sufficient conditions in terms of growth restrictions are given for the solvability of the Dirichlet boundary value problem to forced nonlinear differential equations involving the combination of viscous and dry frictions. Explicit estimates of solutions and their derivatives allow us to restrict ourselves to a sufficiently large neighbourhood of the origin, when formulating these effective conditions. In this way, the behaviour of nonlinearities outside of this neighbourhood can be quite arbitrary. In order to get optimal solvability criteria, the problems with one-term and complete linear differential operators will be treated separately by means of various Green’s functions. The obtained results are compared with some of their analogies of the other authors.
1 Introduction
The Dirichlet problem is one of the mostly studied boundary value problems for differential equations. For ordinary differential equations, the first results due to Hamel, Hammerstein and Lichtenstein were obtained by variational methods (see [1], where also more recent results established in this way are systematically described). After publishing the Schauder fixed point theorem in 1930, topological methods started to be alternatively applied also for vector equations (see e.g. [2–13]). It is well known from the 1930s that, according to the Scorza Dragoni theorem for vector second-order ordinary differential equations with bounded continuous right-hand sides (r.h.s.), the Dirichlet problem is always solvable (see [12, 13]). Thus, for instance as a very particular case, the Dirichlet problem for a forced mathematical pendulum equation, i.e.
where b, \(x_{0}\), \(x_{T}\) and \(T>0\) are real constants, admits a solution, for any b and \(p \in C([ 0, T ], \mathbb{R})\). In fact, \(p \colon[ 0, T ] \rightarrow\mathbb{R}\) can be only Lebesgue integrable, i.e. measurable and \(\int_{0}^{T} |p(t)| \,\mathrm{d}t < \infty\), and, for the same goal, the pendulum equation can involve the viscous damping term, i.e.
(see e.g. [6]).
On the other hand, the Dirichlet problem for a forced linear oscillator with a viscous friction, i.e.
has a unique solution, provided again \(\int_{0}^{T}|p(t)| \,\mathrm{d}t<\infty\), but still the homogeneous problem, namely
should have only a trivial solution which is not always the case; otherwise, the problem is in resonance.
In the presence of a dry friction, the notion of a Carathéodory solution, i.e. the one with an absolutely continuous derivative, is insufficient. The appropriate notion is a Filippov solution which is a Carathéodory solution, but of a differential inclusion with a Filippov regularized right-hand side (see e.g. [14, 15]). For the history and phenomenology of dry friction problems in general, see e.g. [16–20].
For the combination of viscous and dry frictions, it means to consider the problem
when dealing with a Filippov solution of a forced pendulum equation, resp. the problem
where
when dealing with a Carathéodory solution of a forced pendulum inclusion.
In the case of a forced ‘linear’Footnote 1 oscillator, the related problems read as follows:
resp.
Applying the Kakutani-Ky Fan fixed point theorem (see e.g. Theorem II.8.4. on pp.168-169 in [21]), Lasota and Opial (see Theorem 3 in [22]) formulated in 1965 a theorem for a first-order vector problem. This theorem solves problem (2) resp. (1), for any a, b, c, provided only \(\int_{0}^{T}|p(t)| \,\mathrm{d}t < \infty\), as well as problem (5) resp. (4), provided again \(\int_{0}^{T}|p(t)| \,\mathrm{d}t < \infty\), and additionally requiring either \(|b| T\) to be sufficiently small or a trivial solvability of a homogeneous problem
On the other hand, the estimates of solutions and derivatives are not indicated explicitly in [22].
There are also further related results concerning problems (2) and (5) obtained mainly by means of degree arguments (see e.g. [23–34]), but none of them allows us to get such simple criteria as those in [22].
For instance, according to Theorem 6.1 in [25] as well as Theorem 2 in [33], where the combination of sign and growth restriction was employed, both problems (2) and (5) are solvable only in the lack of a viscous friction, i.e. only for \(a=0\), provided \(\int_{0}^{T} |p(t)|\,\mathrm{d}t <\infty\) for (2), and \(\int_{0}^{T} |p(t)|\,\mathrm{d}t <\infty\) together with \(b \leq0\) for (5). Moreover, the estimates of solutions are again not indicated explicitly in [25, 33].
Similarly, in Theorem 4.1 in [26], Theorem 4.1 in [29] and Theorem 3.1 in [30], the viscous friction term \(a x'\) cannot be involved in (2) and (5) but, under more restrictive assumptions than those in [25, 33], the solution estimates are available there.
Using the Hartman-type conditions (called along the lines of the classical paper [8]), as in Corollary 4.1 in [23], the related sufficient conditions for the solvability of (2) are again more restrictive than those in Theorem 3 in [22], namely \(b<0\), p is continuous on \([ 0, T]\), and
where
but a can be this time different from zero, i.e. \(a \neq0\).
Those for the solvability of (5) hold, according to Corollary 4.1 in [23], without (6), i.e. with \(b<0\) and p is continuous on \([0, T ]\) (\(\Rightarrow|p(t)| \leq P\), \(t \in[ 0, T]\)), only.
In the absence of a viscous friction, i.e. when \(a=0\), the forcing term p can be, according to Corollary 4.3 in [23] and Theorem 3 in [30], Lebesgue measurable and essentially bounded, for both problems (2) and (5). For \(b=0\), however, it should be \(c=0\), \(p(t) \equiv0\) and \(k_{3}=0\) or \(a=0\) in order Corollaries 4.1 and 4.3 in [23] and Theorem 3 in [30] to be applied to (2) and (5) which reduces to a trivial solvability of (2) and (5).
On the other hand, under the above assumptions with \(b<0\), the solution \(x(\cdot)\) estimates for (2) and (5) can be this time expressed explicitly as
and
respectively.
If further possibly discontinuous nonlinearities or multivalued maps are implemented into the right-hand sides of given differential equations or inclusions whose growth has, for instance, a superlinear character sufficiently far from the origin, then Theorem 3 in [22] does not any longer apply. Having, however, to our disposal explicit estimates of solutions like (7) or (8) and their derivatives, we can formulate criteria such that the implemented new terms can behave in an arbitrary way outside of the domains characterized by these estimates. In this way, all the results under our consideration can be naturally extended.
This will be therefore our main aim of the present paper. Of course, for obtaining the explicit estimates of solutions and their derivatives, we should study, unlike in [22], exclusively the second-order Dirichlet boundary value problems. For the sake of brevity we will prove only one main theorem for the scalar problem. Nevertheless, since all the related proofs of all the cases under our consideration are quite analogous and differ just by the technical details, we can also present the related solvability criteria and some solutions estimates.
On the other hand, we will discuss in detail the advantages and disadvantages of the usage of one-term vs. complete linear differential operator for the associated Green’s functions.
Hence, our paper will be organized as follows. After the preliminaries, the scalar multivalued Dirichlet problem of the form
will be considered. We will distinguish between the one-term linear differential operator \(L x:= x''\) (then the remaining terms \(ax'+bx\) will be considered in Section 3 as a part of a multivalued perturbation of the r.h.s. of the inclusion in (9)) and the complete linear differential operator \(L x:= x''+ax'+bx\) in Section 4, where particular vector analogies of some obtained results are also indicated.
2 Preliminaries
Let X, Y be two metric spaces. We say that φ is a multivalued mapping from X to Y (written \(\varphi\colon X \multimap Y\)) if, for every \(x \in X\), a nonempty closed subset \(\varphi(x)\) of Y is given. We can associate with φ its graph \(\Gamma_{\varphi}\), the subset of \(X \times Y\), defined by
A multivalued mapping \(\varphi\colon X \multimap Y\) is called upper semicontinuous (u.s.c.) if, for each open \(U \subset Y\), the set \(\{x \in X \mid\varphi(x) \subset U \}\) is open in X.
Lemma 1
(cf., e.g., Proposition I.3.15 in [14])
If \(\varphi\colon X \multimap Y\) is u.s.c., then the graph \(\Gamma _{\varphi}\) is a closed subset of \(X \times Y\).
Lemma 2
(cf., e.g., Proposition I.3.16 in [14])
Assume that \(\varphi\colon X \multimap Y\) is a multivalued mapping such that \(\varphi(X) \subset K\), where \(K \subset Y\) is a compact set, and the graph \(\Gamma_{\varphi}\) of φ is closed. Then φ is u.s.c.
Let Y be a separable metric space and \((\Omega, \mathcal{U}, \mu)\) be a measurable space, i.e. a set Ω equipped with a σ-algebra \(\mathcal{U}\) of subsets and a countably additive measure μ on \(\mathcal{U}\). For our needs here, Ω will be a bounded domain in \(\mathbb{R}^{k}\), equipped with the Lebesgue measure. A multivalued mapping \(\varphi\colon\Omega\multimap Y\) is called measurable if \(\{\omega\in\Omega\mid\varphi(\omega ) \subset V \} \in\mathcal{U}\), for each open set \(V \subset Y\).
A multivalued mapping \(F \colon J \times\mathbb{R}^{n} \multimap \mathbb{R}^{m}\) with compact and convex values, where \(J \subset\mathbb {R}\) is a compact interval, is an upper-Carathéodory (shortly, u-Carathéodory) mapping if it satisfies
-
(i)
\(t \multimap F(t,x)\) is measurable, for every \(x \in \mathbb{R}^{n}\),
-
(ii)
\(x \multimap F(t,x)\) is u.s.c., for almost all (a.a.) \(t\in J\),
-
(iii)
\(|y| \leq r(t) (1+|x|)\), for every \((t,x) \in J \times \mathbb{R}^{n}\), and every \(y\in F(t,x)\), where \(r\colon J \rightarrow [0, \infty)\) is a Lebesgue integrable function.
Let \(F \colon X \multimap Y\) be a multivalued mapping and \(f \colon X \rightarrow Y\) be a single-valued mapping. We say that f is a selection of F (written \(f \subset F\)) if \(f(x) \in F(x)\), for every \(x \in X\).
We will employ the following definitions and statements.
Proposition 1
(Castaing representation, cf. Theorem III.7 in [35])
Let X be a separable metric space, \((\Omega, \mathcal{U}, \mu)\) a measurable space, and φ a measurable multivalued mapping from Ω to complete subsets of X. Then there exists a sequence \(\{ f_{n}\}\) of measurable selections of φ such that, for every x,
where the bar denotes the closure in X.
Lemma 3
(cf., e.g., Lemma 7.1 in [36])
Let \(F\colon J\times\mathbb{R}^{n} \multimap\mathbb{R}^{n}\) be a u-Carathéodory multivalued mapping. Then the composition \(F(t, q(t))\) admits, for every \(q \in C(J, \mathbb{R}^{n})\), a single-valued measurable selection.
If \(X \subset Y\) and \(\varphi\colon X \multimap Y\), then a point \(x \in X\) is called a fixed point of φ if \(x \in\varphi (x)\). We set \(\operatorname{Fix} (\varphi) := \{ x \in X \mid x \in\varphi(x) \}\).
Definition 1
Let X and Y be subsets of normed linear spaces and \(\varphi\colon X \multimap Y\) be a multivalued mapping. If Y is convex, then φ is called a Kakutani map, provided φ is u.s.c. with (nonempty) compact, convex values.
Let X and Y be metric spaces. A multivalued mapping \(\varphi\colon X \multimap Y\) is called compact if its image \(\varphi(X)= \bigcup\{\varphi(x) \mid x \in X \}\) is contained in a compact subset of Y.
The following statement, which we state in the form of proposition, is usually called the Kakutani-Ky Fan-type fixed point theorem.
Proposition 2
(cf., e.g., Theorem II.8.4 in [21])
Let C be a convex (not necessarily closed) subset of a normed linear space, and let \(\varphi\colon C \multimap C\) be a compact Kakutani map. Then φ has a fixed point.
Definition 2
Let \(\varphi\colon J\times\mathbb{R}^{n} \multimap\mathbb{R}^{n}\) be a convex-valued, bounded and measurable mapping. Then the mapping
is called the Filippov-like regularization of φ, where μ denotes Lebesgue measure on \(\mathbb{R}^{n}\), \(N \subset\mathbb {R}^{n}\), and \(\overline{\operatorname{conv}}\) denotes the closed convex hull of a set. \(B((t,x),\delta)\) is an open ball of radius \(\delta>0\), centered at a point \((t,x)\).
Observe that Φ is a bounded u.s.c. mapping with (nonempty) compact and convex values.
Let \(P\colon J \multimap\mathbb{R}^{n}\) be a multivalued mapping having a nonempty set of Lebesgue integrable selections. The Aumann integral of P is defined as follows (cf. [37]):
Now, we shall consider the vector boundary value problem (BVP) in \(\mathbb{R}^{n}\):
where L is a linear ordinary differential operator of order 2 and \(f\colon[a, b] \rightarrow\mathbb{R}^{n}\) is a Lebesgue measurable function in \(a \leq t \leq b\), with homogeneous boundary conditions
Proposition 3
(cf., e.g., Theorem 3.2.1 in [38])
If the homogeneous BVP associated with (10), (11), i.e. the one with \(f(t)\equiv0\), has only a trivial solution, then the BVP (10), (11) has exactly one solution which is given by
where \(G \colon[a, b] \times[a, b] \rightarrow\mathbb{R}\) is the Green’s function related to the homogeneous BVP associated with (10), (11).
It will also be convenient to recall the following lemmas.
Lemma 4
(cf., e.g., Theorem 0.3.4 in [39])
Assume that the sequence of absolutely continuous functions \(x_{k}\colon J \rightarrow\mathbb{R}^{n}\), where J is a compact interval, satisfies the following conditions:
-
the set \(\{x_{k}(t) \mid k \in\mathbb{N}\}\) is bounded \(\forall t \in J\),
-
there is a Lebesgue integrable function \(\alpha\colon J \rightarrow\mathbb{R}\) such that
$$\bigl\vert x'_{k}(t)\bigr\vert \leq\alpha(t), \quad \textit{for almost all }t \in J, \forall k \in\mathbb{N}. $$
Then there exists a subsequence \(\{x_{k}\}\) (denoted just the same) convergent to an absolutely continuous function \(x\colon J \rightarrow \mathbb{R}^{n}\) in the following sense:
-
(i)
\(\{x_{k}\}\) uniformly converges to x,
-
(ii)
\(\{x'_{k}\}\) weakly converges in \(L^{1}(J,\mathbb{R}^{n})\) to \(x'\).
Lemma 5
(cf. [40], p.88)
Let \(J \subset\mathbb{R}\) be a compact interval and \(F \colon J \times\mathbb{R}^{n} \multimap\mathbb{R}^{n}\) be a u-Carathéodory mapping. Let \(N_{F} \colon C(J, \mathbb{R}^{n}) \multimap L^{1}(J, \mathbb {R}^{n})\) be the Nemytskiı̌ operator defined as follows:
for every \(x \in C(J, \mathbb{R}^{n})\). Then, if the sequences \(\{x_{i}\} \subset C(J, \mathbb{R}^{n})\) and \(\{f_{i}\} \subset L^{1}(J, \mathbb{R}^{n})\), \(f_{i} \in N_{F}(x_{i})\), \(i \in\mathbb{N}\), are such that \(x_{i} \rightarrow x\) in \(C(J, \mathbb{R}^{n})\) and \(f_{i} \rightarrow f\) weakly in \(L^{1}(J, \mathbb{R}^{n})\), then \(f \in N_{F}(x)\).
3 Main theorem
Let us consider at first the scalar Dirichlet problem for differential equations involving a dry friction (4), where a, b, c, \(x_{0}\), \(x_{T}\) and \(T >0\) are real constants and \(p\colon J \rightarrow \mathbb{R}\), \(J=[0,T]\), is a Lebesgue integrable function.
Since the function \(\operatorname{sgn}(\cdot)\) is discontinuous in the spatial variable, problem (4) need not have a Carathéodory solution, i.e. a function \(x \colon J \rightarrow\mathbb {R}\) with an absolutely continuous derivative, satisfying (4), for almost all \(t \in J\). Therefore, we need another notion of an appropriate solution, namely the one in the sense of Filippov. For this goal, we use the concept of the Filippov-like regularization (see [15]) of spatially discontinuous maps. More precisely, applying Definition 2 to the right-hand side involving spatial discontinuities, we can speak about a solution in the sense of Filippov of the original problem, provided it is a Carathéodory solution of a multivalued problem with a Filippov-like regularized right-hand side.
In our situation, the discontinuous function to be regularized is the function signum. On the basis of the Filippov-like regularization of \(\operatorname{sgn}(\cdot)\), we obtain the multivalued mapping Signum defined in (3), i.e. \(\operatorname{Sgn}(\cdot)\).
One can readily check that the Signum mapping is u.s.c. with compact and convex values. Hence, after the described Filippov-like regularization, problem (4) with a discontinuous function \(\operatorname{sgn}(\cdot)\) becomes multivalued, i.e. (5).
By a Filippov solution of (4), we understand a function \(x(\cdot)\colon J \rightarrow\mathbb{R}\) with absolutely continuous derivative, satisfying problem (5), almost everywhere on J.
For our needs, it will be convenient to consider still the problem involving more multivalued terms, namely (9), where \(F_{1} \colon\mathbb{R}\multimap\mathbb{R}\), \(F_{2} \colon\mathbb {R}\multimap\mathbb{R}\) are u.s.c. maps with compact, convex values and \(P \colon J \multimap\mathbb{R}\) is an Aumann integrable mapping, \(J=[0,T]\).
Making the change of variables \(y(t) = x(t) - v(t)\), where \(v(t) = \frac{x_{T}-x_{0}}{T} t + x_{0}\), we can immediately see that \(x(t)\) is a Carathéodory solution of (9) if and only if \(y(t)\) satisfies
One can also easily check that \(|v(t)| \leq k_{2} := \max\{|x_{0}|, |x_{T}|\} \), for all \(t\in J\), and \(|v'(t)|= k_{3} :=\frac{|x_{T}-x_{0}|}{T}\), for all \(t \in J\).
Our aim at this moment is to prove the existence and localization theorems for problem (9). Since the desired sufficient conditions will tendentiously take the form of growth restrictions, let us start with the problem involving the truncated maps of \(F_{1}\), \(F_{2}\), i.e.
where
Obviously, \(F_{1} \vert_{\overline{B}_{D}} \colon\overline{B}_{D} \multimap\mathbb{R}\) and \(F_{2} \vert_{\overline{B}_{D}} \colon \overline{B}_{D} \multimap\mathbb{R}\) are u.s.c. maps with compact and convex values, where \(\overline{B}_{D} := \{ z \in\mathbb{R}\mid|z| \leq D\}\) is the closed ball, constant D is such that \(D = D_{0}+k_{2}+k_{3}\), where \(D_{0}>0\) is a suitable constant which will be specified later, and \(k_{2}\), \(k_{3}\) are defined above. Thus, the same is true for \(F_{1}^{*}, F_{2}^{*} \colon\mathbb{R}\multimap\mathbb{R}\) defined as above.
Observe that, in this way, \(F_{1}^{*}(z+v(t)) = F_{1}(z+v(t))\), for \(|z| \leq D_{0}\), and \(F_{2}^{*}(z+v'(t)) = F_{2}(z+v'(t))\), for \(|z| \leq D_{0}\).
Hence, let us find sufficient conditions for the solvability of (13). We distinguish two cases in order to separate formally a linear differential operator and a multivalued perturbation.
In this section, we will deal with the problem formally written in the following way:
where \(F(t,y(t),y'(t)):=P(t)+F_{1}^{*}(y(t)+v(t))+F_{2}^{*}(y'(t)+v'(t))-a (y'(t)+v'(t)) -b (y(t)+v(t))-c \operatorname{Sgn} (y'(t)+v'(t) )\), i.e. \(F\colon J \times\mathbb {R}\times\mathbb{R}\multimap\mathbb{R}\) is a u-Carathéodory multivalued mapping and \(L y(t) := y''(t)\).
For a solvability of the multivalued nonlinear problem (14), we use the Schauder linearization device. Thus, we parametrize the right-hand side (r.h.s.) F in order to have a one-parameter family of linear problems. Let
be the set of candidate solutions, where \(\|u\|_{C^{1}}:=\sup\{ |u(t)|+|u'(t)|, t \in J\}\) and \(D_{0}>0\) is a suitable constant.
Then, for each \(q \in Q\), we get a fully linearized problem
where \(F_{q}(t)=P(t)+F_{1}^{*}(q(t)+v(t))+F_{2}^{*}(q'(t)+v'(t))-a q'(t)-a v'(t) -b q(t)-b v(t)-c \operatorname{Sgn} (q'(t)+v'(t) )\). Obviously, \(F_{q} \colon J\multimap\mathbb{R}\) is, for every \(q \in Q\), an Aumann integrable function of t.
From Lemma 3, the existence follows of at least one measurable selection \(f_{q} \subset F_{q}\) of the multivalued composition \(F_{q}(t)\). Thus, we can consider the single-valued linear Dirichlet problem
The homogeneous problem associated with the problem (16), i.e.
has only a trivial solution. Hence, it follows from the Fredholm alternative (see Proposition 3) that problem (16) has a unique Carathéodory solution \(y(\cdot)\), which takes the form
where \(G_{1} \) is the Green’s function of problem (17), i.e.
One can easily check that \(\vert G_{1}(t, s) \vert \leq \frac{T}{4}\), for all \(t, s \in[0, T]\).
In fact, since the multivalued composition \(F_{q}(t)\) is, for every \(q \in Q\), obviously measurable, according to Proposition 1, we can even write
where \(\{f_{n, q}(t) \subset F_{q}(t) \}_{n \in\mathbb {N}}\) is a sequence of measurable selections of \(F_{q}\) and the integral is understood in the sense of Aumann.
Observe that because of
we also have
where the integral is again understood in the sense of Aumann. Moreover, \(\vert \frac{\partial G_{1}}{\partial t} (t, s)\vert \leq1\), for all \(t, s \in[0, T]\).
Hence, denoting by \(\varphi \colon Q \multimap C^{1}(J, \mathbb{R})\) the solution operator of (15), where
instead of the existence of the solution of differential problem (15), we can equivalently investigate the existence of a fixed point of the multivalued operator φ. For this purpose, we will apply the Kakutani-Ky Fan-type fixed point theorem (Proposition 2). In this way, we will prove the following theorem.
Theorem 1
Let a, b, c, \(x_{0}\), \(x_{T}\) and \(T>0\) be real constants such that
Assume that \(P \colon J \multimap\mathbb{R}\) is an Aumann integrable multivalued mapping and \(F_{1} \vert_{\overline{B}_{D}} \colon\overline{B}_{D} \multimap \mathbb{R}\), \(F_{2} \vert_{\overline{B}_{D}} \colon\overline{B}_{D} \multimap\mathbb{R}\) are u.s.c. multivalued mappings with convex and compact values, where \(D=D_{0}+k_{2}+k_{3}\), \(D_{0}>0\) is still supposed to be a suitable constant such that
where
and
Then problem (9) admits a solution \(x(\cdot)\) such that
Proof
In order to check all the assumptions of Proposition 2, we will proceed in four steps.
(i) Since problem (16) is uniquely solvable, the set \(\varphi(Q)\) is nonempty.
(ii) Let us prove that the set \(\varphi(Q)\), i.e. the set of solutions of (15), is relatively compact. According to the well known Arzelá-Ascoli lemma, the set of solutions is relatively compact in \(C^{1}(J, \mathbb{R})\) if and only if it is uniformly bounded and equi-continuous, both in the \(C^{1}\)-norm.
(a) Let us show that the set of solutions of (15) is uniformly bounded in \(C^{1}(J, \mathbb{R})\). Let \(u(\cdot)\) be a solution of (15) and \(f_{q} \subset F_{q}\) be a measurable selection of \(F_{q}\). We could see that such a measurable selection exists and that \(u(\cdot)\) takes the form
It is obvious that
where \(v(t)= \frac{x_{T}-x_{0}}{T} t +x_{0}\) and \(v'(t)=\frac{x_{T}-x_{0}}{T}\).
By means of Lemma 3 and in view of (18), (19), for any \(t \in J\), we obtain the following estimate:
where \(f_{1}^{*} \subset F_{1}^{*}\), \(f_{2}^{*} \subset F_{2}^{*}\), \(p \subset P\) and \(\operatorname{Sgn}_{\mathrm{Sel}} \subset\operatorname{Sgn}\) are the related measurable selections.
Since this estimate holds in the same way for all \(q \in Q\), this already means that the solutions \(u(\cdot)\) of (15), i.e. the set \(\varphi(Q)\), are uniformly bounded in the \(C^{1}\)-norm.
Moreover, according to (20), (22), if there exists a positive constant \(D_{0}\) such that \(D_{0} \geq\Delta_{1}(D_{0})\), then the set \(\varphi(Q)\) satisfies \(\varphi(Q) \subset Q\).
(b) Now, let us show that the elements of the set \(\varphi(Q)\), i.e. the solutions \(u(\cdot)\) of (15), are equi-continuous in the \(C^{1}\)-norm. For any \(t_{1}, t_{2} \in J\) with \(t_{1}< t_{2}\), we have by means of Lemma 3 and in view of (19)
Furthermore, we have still
Therefore, the solutions \(u(\cdot)\) of (15) are equi-continuous in the \(C^{1}\)-norm. Summing up (a) and (b), the elements of the set \(\varphi(Q)\) are relatively compact in the \(C^{1}\)-norm, as claimed.
(iii) We will show that the operator φ is u.s.c. In view of Lemma 2, and since φ was shown to be compact, it is sufficient to show that the graph \(\Gamma_{\varphi}\) is closed. Let \(\{(q_{k}, u_{k})\} \subset\Gamma_{\varphi}\) be a sequence such that \(\{(q_{k}, q'_{k}, u_{k})\} \rightarrow(q, q', u)\), where \(q \in Q\). For all \(k \in\mathbb{N}\) and a.a. \(t \in J\), the sequence \(\{u_{k}'\}\) is bounded and \(|u_{k}''(t)|\leq |p(t)|+M_{1}(D_{0})+M_{2}(D_{0})+k_{1} D_{0}+|c|+|a| k_{3}+|b| k_{2}\), for a.a. \(t \in J\). The sequence \(\{w_{k}:=u'_{k}\}\) satisfies all the assumptions of Lemma 4.
Thus, applying Lemma 4 to the sequence \(\{w_{k}:=u'_{k}\}\), we find that there exists a subsequence of \(\{u_{k}'\}\), for the sake of simplicity denoted in the same way as the sequence, which converges uniformly to \(u'\) on J and such that \(\{u''_{k}\}\) converges weakly to \(u''\) in \(L^{1}(J,\mathbb{R})\).
If we set \(z_{k}:=(u_{k}, w_{k})\), then \(z'_{k}= (u_{k}', w_{k}')=(u_{k}', u_{k}'')\rightarrow(u', u'')\), weakly in \(L^{1}(J, \mathbb{R}^{2})\). Let us consider the system
where \(z_{k}'(t) = (u_{k}'(t), w_{k}'(t))\) and \(H(t, q_{k}(t), q_{k}'(t)) = (w_{k}, F_{q_{k}} (t))\).
Applying Lemma 5, for \(f_{i}:=z'_{k}\), \(f := (u', u'')\), \(x_{i}:=(q_{k}, q'_{k})\), it follows that
i.e. \(u''(t) \in F_{q}(t)\), for a.a. \(t \in J\).
The set \(\varphi(Q)\) is relatively compact and graph \(\Gamma_{\varphi }\) is closed. Therefore, the mapping φ is u.s.c., compact and, in particular, with compact values.
(iv) Finally, we will show that the mapping φ has convex values.
Let \(u_{1}\), \(u_{2}\) be two distinct solutions of problem (15) associated with measurable selections \(f_{1,q}, f_{2,q} \subset F_{q}\). Then, for all \(t \in J\), we have
Let \(\lambda\in[ 0, 1 ]\) be arbitrary. Since F is convex-valued, the same is true for \(F_{q}\), and subsequently
must be a measurable selection of \(F_{q}\), i.e. \(f_{q} \subset F_{q}\). We have
Thus, \(u(\cdot)\) must also be a solution of (15), by which the mapping φ has convex values.
After all, applying Proposition 2, we obtain the existence of a fixed point of the multivalued mapping φ which represents a solution of problem (13). However, because of the definitions of \(F_{1}^{*}\), \(F_{2}^{*}\), such a solution must be a solution of problem (12) as well as (9). □
Example 1
As an illustrative example of a Dirichlet problem satisfying rather implicit conditions (20)-(22), we can consider (9), where
In this way, conditions (21), (22) take, for \(D_{0}=\frac{1}{4}\), the form:
Taking, furthermore, \(\mathcal{P}\leq\frac{\sqrt{2}-1}{16}\) and \(|c|=\frac{1}{32}\), we get for \(F_{1}(z):=f_{1}z^{m}\) and \(F_{2}(z):=f_{2}z^{n}\) the inequality:
Hence, if \(|f_{1}|,|f_{2}|\leq\frac{1}{16}\) and \(m,n\geq1\), then we arrive at (21), (22), because for \(|f_{1}|=|f_{2}|=\frac{1}{16}\) and \(m=n\), we have
which holds for all \(m\geq1\).
After all, under these assumptions, problem (9) admits a solution \(x(\cdot)\) such that
In fact, the same is true for (9), where \(F_{1}\), \(F_{2}\) are as above for \(z\in[-\frac{1}{2},\frac{1}{2}]\) and can be arbitrary outside of the interval \([-\frac{1}{2},\frac{1}{2}]\).
Remark 1
If, additionally, a multivalued mapping P is essentially bounded, then condition (20) can be replaced, for the same conclusion for (9), by
and (21) with \(\Delta_{1}(D_{0})\) in (22) can be replaced by \(\widetilde{D}_{0} \geq\Delta'_{1}(\widetilde{D}_{0}) \), where
\(\mathcal{P}_{1}:=\operatorname{ess}\sup_{t \in J} | P(t)|\), because
Thus, the same conclusion holds for (9) with D replaced by \(\widetilde{D}=\widetilde{D}_{0}+k_{2}+k_{3}\).
Remark 2
If, additionally, a multivalued mapping \(F_{2}\) is bounded, i.e. \(M_{2} := \max_{z \in\mathbb{R}} |F_{2}(z)|\), and \(a=0\), we can still improve the estimate for a solution of (9). Condition (20) can then be replaced by \(\frac{4}{T^{2}} > |b|\), provided there still exists a positive constant \(D_{1}\) such that
Under these assumptions, problem (9) has a solution \(x(\cdot )\) such that
where \(M_{1}(D_{1}):=\max_{|z|\leq D_{1}+k_{2}}|F_{1}(z)|\), \(k_{2}:=\max\{|x_{0}|, |x_{T}|\}\), \(k_{3}:=\frac{|x_{T}-x_{0}|}{T}\).
Similarly if, additionally, a multivalued mapping \(F_{1}\) is bounded, i.e. \(M_{1}:= \max_{z \in\mathbb{R}} |F_{1}(z)|\), and \(b=0\), then condition (20) can be replaced by \(\frac{1}{T}>|a|\), provided there still exists a positive constant \(D_{2}\) such that
where \(M_{2}(D_{2}):=\max_{|z|\leq D_{2}+k_{3}}|F_{2}(z)|\). Here, problem (9) has a solution \(x(\cdot)\) such that
4 Discussion of further possibilities
Now, we will make comments on the existence of a solution of the problem (9), but this time considered in the following form:
where \(F(t, x(t), x'(t)):=P(t)+F_{1}(x(t))+F_{2}(x'(t))-c\operatorname{Sgn} x'(t) \) is a u-Carathéodory multivalued mapping and \(L x(t):= x''(t) +a x'(t)+bx(t)\). This allows us to consider different Green’s functions associated with the complete linear differential operator.
In this way, the related solvability criteria are as follows:
-
(1)
For \(a^{2}-4b>0\): \(D_{0} \geq\Delta_{2}(D_{0})\), where
$$\begin{aligned} \Delta_{2}(D_{0}) :=& \frac{e^{ (\lambda_{1}-\lambda_{2} )T}}{\sqrt{a^{2}-4 b}} \bigl[1+\vert \lambda_{1}\vert +|\lambda_{2}| \bigr] \\ &{}\times \bigl[\mathcal{P}+T \bigl(M_{1}(D_{0})+M_{2}(D_{0})+|c|+|a| k_{3}+|b| k_{2} \bigr) \bigr], \end{aligned}$$\(\lambda_{1}=\frac{-a+\sqrt{a^{2}-4b}}{2}\), \(\lambda_{2}=\frac{-a-\sqrt {a^{2}-4b}}{2}\).
-
(2)
For \(a^{2}-4b=0\): \(D_{0} \geq\Delta_{3}(D_{0})\), where
$$\begin{aligned} \begin{aligned} \Delta_{3}(D_{0}) :={} &e^{\frac{|a|}{2} T} \biggl[1+ \frac{T}{4} \bigl(1+2\vert a\vert \bigr) \biggr] \\ &{}\times \bigl[\mathcal{P}+T \bigl(M_{1}(D_{0})+M_{2}(D_{0})+|c|+|a |k_{3}+|b| k_{2} \bigr) \bigr]. \end{aligned} \end{aligned}$$ -
(3)
For \(a^{2}-4b<0\) and \(T \neq \frac{2\pi k}{\sqrt{4b-a^{2}}}\), \(k \in \mathbb{N}\): \(D_{0} \geq\Delta_{4}(D_{0})\), where
$$\begin{aligned} \Delta_{4}(D) :=& \frac{e^{\frac{|a|}{2} T} [2+|a|+\sqrt {4b-a^{2}} ]}{ \sqrt{4b-a^{2}}\vert \sin(\frac{T}{2}\sqrt{4b-a^{2}})\vert } \\ &{}\times \bigl[\mathcal{P}+T \bigl(M_{1}(D_{0})+M_{2}(D_{0})+|c|+|a| k_{3}+|b| k_{2} \bigr) \bigr]. \end{aligned}$$
In all the cases (1), (2), (3), the problem (9) admits a solution \(x(\cdot)\) such that \(\max_{t \in J}\{|x(t)|+|x'(t)|\} \leq D_{0}+k_{2}+k_{3}\), where \(k_{2}:= \max\{|x_{0}|, |x_{T}|\}\), \(k_{3}:=\frac{|x_{T}-x_{0}|}{T}\).
In view of Remark 2 and the estimates above, we can immediately formulate the following corollary for a Lebesgue integrable selection \(p \subset P\) and \(F_{1}(x) := d\sin x\).
Corollary 1
Let b, c, d, \(x_{0}\), \(x_{T}\) and \(T>0\) be real constants, and \(p \colon[ 0, T] \rightarrow\mathbb{R}\) be a Lebesgue integrable mapping, where \(\mathcal{P} := \int_{0}^{T} |p(t)|\, \mathrm{d}t<\infty\). Then there exists a Filippov solution \(x(\cdot)\) of problem
such that, for \(0< b \neq (\frac{k\pi}{T} )^{2}\), \(k \in \mathbb{N}\),
For \(b<0\) and \(b=0\), we obtain the respective estimates
and
where \(k_{2}=\max\{|x_{0}|, |x_{T}|\}\).
Remark 3
Observe that, for \(b = 0\), the conclusion of Corollary 1 reduces to the Dirichlet problem for a forced pendulum with a dry friction while, for \(d = 0\), the result concerns the Dirichlet problem for a forced ‘linear’ oscillator with a dry friction. Moreover, the non-resonance condition \(b \neq (\frac{k \pi}{T} )^{2}\), \(k \in\mathbb{N}\), distinguished in Corollary 1 in three cases (\(b>0\), \(b<0\), \(b=0\)), is evidently optimal, when comparing it with all its analogies discussed above (cf. (20), Remark 1 and Remark 2).
Finally, we can also consider a particular case with \(c=0\) of the Dirichlet problem (9) in a vector form in \(\mathbb{R}^{n}\), i.e. let a, b (\(c=0\)) and \(T>0\) remain real constants, but \(x_{0}, x_{T} \in\mathbb{R}^{n}\), \(F_{1} \colon\mathbb{R}^{n} \multimap \mathbb{R}^{n}\), \(F_{2} \colon\mathbb{R}^{n} \multimap\mathbb{R}^{n}\) be u.s.c. with convex and compact values and \(P \colon J \multimap \mathbb{R}^{n}\) be an Aumann (component-wise)-integrable mapping. In this way, dry friction terms can be involved, after their Filippov’s regularization, in some components of \(F_{2}\).
All the above statements can then be reformulated for the vector case in the appropriate norm \(\|\cdot\|\), when the truncated vector analogies of \(F_{1}^{*}\), \(F_{2}^{*}\) in (13) are defined as follows (\(i=1, 2\)):
Observe that in this way \(F_{i}^{*}\) is obviously u.s.c. with nonempty, compact, convex values, provided \(F_{i} \vert_{\overline{B}_{D}}\) is so, where \(\overline{B}_{D}:=\{z \in\mathbb{R}^{n} \mid\|z\| \leq D\}\). Moreover, \(M_{i}(D):=\max_{z \in\overline{B}_{D}}\|F_{i}(z)\|=\max_{z \in \mathbb{R}^{n}}\|F_{i}^{*}(z)\|\), \(i=1, 2\).
Since \(\overline{B}_{D}\) is an absolute retract, there exists a continuous retraction \(r \colon\mathbb{R}^{n} \rightarrow\overline {B}_{D}\) such that \(r \vert_{\overline{B}_{D}} = \mathrm{id}\vert _{\overline{B}_{D}}\) (see e.g. Chapter I.2 in [14]). Therefore, we can alternatively replace \(F_{i}^{*}(z)\) by \(F_{i}(r(z))\), where (\(i= 1, 2\))
because the compositions \(F_{i}(r(z))\), \(i=1, 2\), must also be u.s.c. maps with convex, compact values with the same properties (see e.g. Chapter I.3.3 in [14]).
Notes
Because of the terms \(c\operatorname{sgn}x'(t)\) resp. \(c\operatorname{Sgn} x'(t)\), the related equations resp. inclusions are no longer linear but, for the sake of simplicity, we reserve for them this name.
References
Mawhin, J: Problèmes de Dirichlet variationnels non linéaires. Partie 1 des comptes rendus du cours d’ été OTAN ‘Variational Methods in Nonlinear Problems’. Presses University Montreal, Montreal (1987)
Andres, J, Sanchez, L: A note on vector boundary value problems. Nonlinear Differ. Equ., TMA 3(1-2), 49-58 (1997)
Bernfeld, SR, Lakshmikantham, V: An Introduction to Nonlinear Boundary Value Problems. Academic Press, New York (1974)
Erbe, LH, Knobloch, HW: Boundary value problems for systems of second order differential equations. Proc. R. Soc. Edinb. A 101(1-2), 61-76 (1985)
Fabry, C, Habets, P: The Picard boundary value problem for nonlinear second order vector differential equations. J. Differ. Equ. 42(2), 186-198 (1981)
Granas, A, Guenther, RB, Lee, JW: Nonlinear Boundary Value Problems for Ordinary Differential Equations. Dissertationes Math., vol. 244. PWN, Warsaw (1985)
Habets, P, Schmitt, K: Nonlinear boundary value problems for system of differential equations. Arch. Math. (Basel) 40(1), 441-446 (1983)
Hartman, P: On boundary value problems for systems of ordinary nonlinear second order differential equations. Trans. Am. Math. Soc. 96(3), 493-509 (1960)
Knobloch, HW, Schmitt, K: Nonlinear boundary value problems for systems of differential equations. Proc. R. Soc. Edinb. A 78(1-2), 139-159 (1977)
Mawhin, J: Boundary value problems for nonlinear second-order vector differential equations. J. Differ. Equ. 16(2), 257-269 (1974)
O’Regan, D: Boundary value problems for second and higher order differential equations. Proc. Am. Math. Soc. 113(3), 761-775 (1991)
Scorza-Dragoni, G: Sul problema dei valori ai limiti per i sistemi di equazioni differenziali del secondo ordine. Boll. Unione Mat. Ital. 14, 225-230 (1935)
Scorza-Dragoni, G: Sui sistemi di equazioni integrali non lineari. Rend. Semin. Mat. Univ. Padova 7, 1-35 (1936)
Andres, J, Górniewicz, L: Topological Fixed Point Principles for Boundary Value Problems. Kluwer Academic, Dordrecht (2003)
Filippov, AF: Differential Equations with Discontinuous Righthand Sides. Kluwer Academic, Dordrecht (1988)
Kunze, M: Non-Smooth Dynamical Systems. Lecture Notes in Mathematics, vol. 1744. Springer, Berlin (2000)
Monteiro Marques, MDP: Differential Inclusions in Nonsmooth Mechanical Problems: Shocks and Dry Friction. Birkhäuser, Basel (1993)
Popp, K, Hinrichs, N, Oestreich, M: Dynamical behaviour of a friction oscillator with simultaneous and external excitation. Sādhanā 20(2-4), 627-654 (1995)
Stewart, DE: Rigid-body dynamics with friction and impact. SIAM Rev. 42(1), 3-39 (2000)
Zhuravlev, VP: On the history of the dry friction law. Mech. Solids 48(4), 364-369 (2013)
Dugundji, J, Granas, A: Fixed Point Theory. Springer, Berlin (2003)
Lasota, A, Opial, Z: An application of the Kakutani-Ky Fan theorem in the theory of ordinary differential equations. Bull. Acad. Pol. Sci., Sér. Sci. Math. Astron. Phys. 13(11-12), 781-786 (1965)
Andres, J, Malaguti, L, Pavlačková, M: Hartman-type conditions for multivalued Dirichlet problem in abstract spaces. Discrete Contin. Dyn. Syst., Ser. A (to appear)
Cubiotti, P, Yao, JC: Two-point problem for vector differential inclusions with discontinuous right-hand side. Appl. Anal. 93(9), 1811-1823 (2014)
Granas, A, Guenther, RB, Lee, JW: Some general existence principles in the Carathéodory theory of nonlinear differential systems. J. Math. Pures Appl. 70(2), 153-196 (1991)
Kožušníková, M: A bounding functions approach to multivalued Dirichlet problem. Atti Semin. Mat. Fis. Univ. Modena Reggio Emilia 55, 1-19 (2007)
Lee, JW, O’Regan, D: Existence of solutions to some initial value, two point and multipoint boundary value problems with discontinuous nonlinearities. Appl. Anal. 33(1-2), 57-77 (1989)
Nistri, P: Positive solutions of non-linear eigenvalue problem with discontinuous non-linearity. Proc. R. Soc. Edinb. A 83(1-2), 133-145 (1979)
Pavlačková, M: A bound sets technique for Dirichlet problem with an upper-Carathéodory right-hand side. Acta Univ. Palacki. Olomuc., Fac. Rerum Nat., Math. 49(2), 95-106 (2010)
Pavlačková, M: A Scorza-Dragoni approach to Dirichlet problem with an upper-Carathéodory right-hand side. Topol. Methods Nonlinear Anal. 44(1), 239-247 (2014)
Pruszko, T: Some Applications of the Topological Degree Theory to Multi-Valued Boundary Value Problems. Dissertationes Math., vol. 229. PWN, Warsaw (1984)
Stuart, CA: Differential equations with discontinuous non-linearities. Arch. Ration. Mech. Anal. 63(1), 59-75 (1976)
Šenkyřík, M, Guenther, RB: Boundary value problems with discontinuities in the spatial variable. J. Math. Anal. Appl. 193(1), 296-305 (1995)
Zuev, AV: On the Dirichlet problem for a second-order ordinary differential equation with discontinuous right-hand side. Differ. Equ. 42(3), 340-346 (2006)
Castaing, C, Valadier, M: Convex Analysis and Measurable Multifunctions. Springer, Berlin (1977)
Appell, J, De Pascale, E, Nguyen, HT, Zabreiko, PP: Multi-Valued Superpositions. Dissertationes Math., vol. 345. PWN, Warsaw (1995)
Aumann, RJ: Integrals of set-valued functions. J. Math. Anal. Appl. 12(1), 1-12 (1965)
Stakgold, I: Green’s Functions and Boundary Value Problems. Wiley, New York (1979)
Aubin, J-P, Cellina, A: Differential Inclusions. Springer, Berlin (1984)
Vrabie, II: Compactness Methods for Nonlinear Evolutions. Longman, Harlow (1990)
Acknowledgements
The first author was supported by the grant No. 14-06958S ‘Singularities and impulses in boundary value problems for nonlinear ordinary differential equations’ of the Grant Agency of the Czech Republic. The second author was supported by the grant IGA Mathematical Models, IGA_PrF_2014028 and by the grant No. 14-06958S ‘Singularities and impulses in boundary value problems for nonlinear ordinary differential equations’ of the Grant Agency of the Czech Republic.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally in this article. They read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Andres, J., Machů, H. Dirichlet boundary value problem for differential equations involving dry friction. Bound Value Probl 2015, 106 (2015). https://doi.org/10.1186/s13661-015-0371-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-015-0371-z