Skip to main content

Triple positive solutions for two classes of delayed nonlinear fractional FDEs with nonlinear integral boundary value conditions

Abstract

This paper is concerned with two classes of delayed nonlinear fractional functional differential equations (FDEs) with nonlinear Riemann-Stieltjes integral boundary value conditions. By employing the well-known Leggett-Williams fixed point theorem and a generalization of Leggett-Williams fixed point theorem, some new sufficient criteria are established to guarantee the existence of at least triple positive solutions. As applications, some interesting examples are presented to illustrate our main results.

1 Introduction

Fractional differential equations arise in many engineering and scientific disciplines as the mathematical modeling of systems and processes in the fields of physics, chemistry, aerodynamics, electrodynamics of complex medium, polymer rheology, Bode’s analysis of feedback amplifiers, capacitor theory, electrical circuits, electron-analytical chemistry, biology, control theory, fitting of experimental data, and so forth. Fractional differential equations also serve as an excellent tool for the description of hereditary properties of various materials and processes. For example, in physics, the traditional way to deal with the behavior of certain materials under the influence of external forces in mechanics is to use the laws of Hooke and Newton. If we are dealing with viscous liquids, then we can use Newton’s law \(\eta\varepsilon'(t)=\sigma(t)\), where \(\sigma(t)\) and \(\varepsilon(t)\) denote stress and strain at time t respectively, η is the so-called viscosity of the material. In view of some possible interpolation properties, it is natural for us to design the classical Newton’s law according to

$$ \eta D^{k}_{0^{+}}\varepsilon(t)=\sigma(t),\quad k\in(n-1,n), n \in\mathbb{N}, $$

which is called Nutting’s law [1]. As a consequence, the subject of fractional differential equations is gaining much importance and attention. Especially, the boundary value problems of fractional differential equations have been one of the hottest problems. There have been many papers focused on boundary value problems of fractional ordinary differential equations; see [1–16]. Moreover, the boundary value problems with Riemann-Stieltjes integral boundary condition arise in a variety of different areas of applied mathematics and physics (for more comments on Stieltjes integral boundary condition and its importance, we refer the reader to the papers by Webb and Infante [17, 18] and their other related works). For example, blood flow problems, chemical engineering, thermo-elasticity, underground water flow, population dynamics, and so on can be reduced to nonlocal integral boundary problems. By means of some well-known fixed point theorems, some papers deal with the existence and multiplicity of solutions or positive solutions for this type of boundary value problems involving fractional differential in the recent references (see [19–25]).

In the real world, the time-delay phenomenon exists commonly and is inevitable. Many changes and processes not only depend on the present status but also on the past status. Therefore, it is necessary to consider the time-delay effect in the mathematical modeling of fractional differential equations. To the best of our knowledge, there are rare papers dealing with the existence of positive solutions for fractional Riemann-Stieltjes integral BVPs with time-delays by the well-known Leggett-Williams fixed point theorem. Therefore, the main goal of this paper is to study the existence of multiple positive solutions of Riemann-Stieltjes integral boundary value problems (BVP for short) involving time-delays for two classes of nonlinear Caputo fractional differential equations (1.1) and (1.2) as follows:

$$ \left \{\textstyle\begin{array}{l} D_{0+}^{q}u(t)+f(t,u,u_{t})=0, \quad t\in I, 2< q\leq3, \\ \alpha u(0)-\beta u'(0)=g_{1}(\int_{0}^{1}u(s)\, dA_{1}(s)),\qquad u''(0)=0, \\ \gamma u(1)+\delta u'(1)=g_{2}(\int_{0}^{1}u(s)\, dA_{2}(s)), \\ u(s)=\phi(s), \quad s\in[-\tau,0]\triangleq J \end{array}\displaystyle \right . $$
(1.1)

and

$$ \left \{\textstyle\begin{array}{l} D_{0+}^{q}u(t)+g(t,u,u_{t},u')=0,\quad t\in I, 2< q\leq3, \\ \alpha u(0)-\beta u'(0)=\int_{0}^{1}h_{1}(s,u(s))\, dA_{1}(s), \qquad u''(0)=0, \\ \gamma u(1)+\delta u'(1)=\int_{0}^{1}h_{2}(s,u(s))\, dA_{2}(s), \\ u(s)=\phi(s),\quad s\in[-\tau,0]\triangleq J, \end{array}\displaystyle \right . $$
(1.2)

where \(I\triangleq[0,1]\), \(D_{0+}^{q}\) is the standard Caputo fractional derivative of fractional order q. \(f\in C([0,1]\times\mathbb{R}^{2},\mathbb{R}^{+})\), \(g_{i}\in C(\mathbb{R}^{+},\mathbb{R}^{+})\) (\(i=1,2\)), \(g\in C([0,1]\times\mathbb{R}^{3},\mathbb{R}^{+})\), \(h_{i}\in C([0,1]\times\mathbb{R},\mathbb{R}^{+})\) (\(i=1,2\)), \(\mathbb{R} = (-\infty,+\infty)\), \(\mathbb{R}^{+} =[0,+\infty)\). α, β, γ, δ are all nonnegative constants with \(\rho\triangleq\alpha\gamma+\alpha\delta+\beta\gamma>0\). \(\int_{0}^{1}u(s)\, dA_{i}(s)\) (\(i=1, 2\)) denotes the Riemann-Stieltjes integrals. \(A_{i}: [0,1]\rightarrow \mathbb{R}\) (\(i=1, 2\)) is the increasing function of bounded variation. \(\tau>0\) is the constant time-delay. \(\phi(t)\in C_{\tau}\) (\(C_{\tau}\) will be given in Section 3), \(u_{t}\in C_{\tau}\), \(u_{t}(\theta)=u(t+\theta)\), \(\theta\in[-\tau,0]\).

The rest of this paper is organized as follows. In Section 2, we introduce some definitions and lemmas to prove our main results. In Section 3, one sufficient condition is given by the well-known Leggett-Williams fixed point theorem to guarantee the existence multiple positive solutions for BVP (1.1). Applying a generalization of the Leggett-Williams fixed point theorem, we establish the existence of at least three positive solutions for BVP (1.2) in Section 4. As applications, some interesting examples are presented to illustrate the main results in Section 5.

2 Preliminaries and statements

For the convenience of the reader, we state some background materials from the theory of both fractional calculus and cones in Banach spaces. These definitions and properties can be found in the literature.

Definition 2.1

(see [26, 27])

The Riemann-Liouville fractional integral of order \(\alpha>0\) of a function \(f:(0,\infty)\rightarrow \mathbb{R}\) is given by

$$I_{0+}^{\alpha}f(t)=\frac{1}{\Gamma(\alpha)}\int_{0}^{t}(t-s)^{\alpha-1}f(s) \, ds $$

provided that the right-hand side is pointwise defined on \((0,\infty)\).

Definition 2.2

(see [26, 27])

The Caputo fractional derivative of order \(\alpha>0\) of a continuous function \(f:(0,\infty)\rightarrow \mathbb{R}\) is given by

$$D_{0+}^{\alpha}f(t)=\frac{1}{\Gamma(n-\alpha)}\int_{0}^{t} \frac {f^{(n)}(s)}{(t-s)^{\alpha-n+1}}\, ds, $$

where \(n-1<{\alpha}\leq n\), provided that the right-hand side is pointwise defined on \((0,\infty)\).

Lemma 2.1

(see [26])

Assume that \(u\in C(0,1)\cap L(0,1)\) with a Caputo fractional derivative of order \(\alpha>0\) that belongs to \(u\in C^{n}[0,1]\), then

$$I_{0+}^{\alpha}D_{0+}^{\alpha}u(t)=u(t)+c_{0}+c_{1}t+ \cdots+c_{n-1}t^{n-1} $$

for some \(c_{i}\in\mathbb{R}\) (\(i=0,1,\ldots,n-1\)), here n is the smallest integer greater than or equal to α.

Definition 2.3

Let X be a real Banach space. A nonempty closed convex set \(P\subset X\) is called a cone if it satisfies the following two conditions:

  1. (i)

    \(x\in X\), \(\lambda\geq0\) implies \(\lambda x\in P\);

  2. (ii)

    \(x\in P\), \(-x\in P\) implies \(x=0\).

Every cone \(P\subset X\) induces an ordering in X given by \(x\preceq y\) if and only if \(y-x\in P\).

Definition 2.4

The map ψ is said to be nonnegative continuous concave on a cone P of a real Banach space E provided that \(\psi:P \rightarrow [0,\infty)\) is continuous and for all \(u, v\in P\), \(\lambda\in[0,1]\) such that

$$ \psi\bigl(\lambda u+(1-\lambda)v\bigr)\geq\lambda\psi(u)+(1-\lambda)\psi(v). $$

Similarly, the map Ï– is said to be nonnegative continuous convex on a cone P of a real Banach space E provided that \(\varpi:P \rightarrow [0,\infty)\) is continuous and for all \(u, v\in P\), \(\lambda\in[0,1]\) such that

$$ \varpi\bigl(\lambda u+(1-\lambda)v\bigr)\leq \lambda\varpi(u)+(1-\lambda) \varpi(v). $$

Let P be a cone in a real Banach space E. For \(c>0\), \(0< b< d\), we define

$$ P_{c}=\bigl\{ x\in P:\Vert x\Vert < c\bigr\} , \qquad \overline{P}_{c}=\bigl\{ x\in P:\Vert x\Vert \leq c\bigr\} $$

and

$$ P(\psi,b,d)=\bigl\{ x\in P:b\leq\psi(x),\|x\|\leq d\bigr\} . $$

It is easy to see that \(P(\psi,b,d)\) is a convex and closed subset of P.

Lemma 2.2

(Leggett-Williams fixed point theorem, see [28])

Let P be a cone in a real Banach space E, ψ be a nonnegative continuous concave functional on P such that \(\psi(x)\leq\|x\|\) for \(x\in\overline{P}_{c}\). Suppose that \(A:\overline{P}_{c} \rightarrow \overline{P}_{c}\) is completely continuous and there exist \(0< a< b< d\leq c\) such that

  1. (i)

    \(\{x\in P(\psi,b,d):\psi(x)>b\}\neq\emptyset\) and \(\psi(Ax)>b\) for all \(x\in P(\psi,b,d)\);

  2. (ii)

    \(\|Ax\|< a\) for all \(x\in\overline{P}_{a}\);

  3. (iii)

    \(\psi(Ax)>b\) for all \(x\in P(\psi,b,c)\) with \(\|Ax\|>d\).

Then A has at least three fixed points \(x_{1}\), \(x_{2}\) and \(x_{3}\) satisfying

$$ \|x_{1}\|< a< \|x_{3}\|, \qquad \psi(x_{3})< b< \psi(x_{2}). $$

Next, we are prepared to state an important generalization of the Leggett-Williams fixed point theorem, which comes from Bai and Ge in [29].

Let ψ be a nonnegative continuous concave functional on P, and let ϖ and ω be nonnegative continuous convex functionals on P. For nonnegative real numbers r, a and l, we define the following convex sets:

$$\begin{aligned}& P(\varpi,r;\omega,l)=\bigl\{ u\in P: \varpi(u)< r, \omega(u)< l\bigr\} , \\& \bar{P}(\varpi,r;\omega,l)=\bigl\{ u\in P: \varpi(u)\leq r, \omega (u)\leq l\bigr\} , \\& P(\varpi,r;\omega,l;\psi,a)=\bigl\{ u\in P: \varpi(u)< r, \omega(u)< l, \psi(u)>a \bigr\} , \\& \bar{P}(\varpi,r;\omega,l;\psi,a)=\bigl\{ u\in P: \varpi(u)\leq r, \omega(u)\leq l, \psi(u)\geq a\bigr\} . \end{aligned}$$

Lemma 2.3

(see [29])

Let P be a cone in a real Banach space E. Assume that constants \(r_{1}\), b, d, \(r_{2}\), \(l_{1}\) and \(l_{2}\) satisfy \(0< r_{1}< b< d\leq r_{2}\) and \(0< l_{1}\leq l_{2}\). If there exist two nonnegative continuous convex functionals ϖ and ω on P and a nonnegative continuous concave functional ψ on P such that:

(A1):

there exists \(M>0\) such that \(\|u\|\leq M\max\{\varpi (u),\omega(u)\}\) for all \(u\in P\);

(A2):

\(P(\varpi,r; \omega,l)\neq\emptyset\) for any \(r>0\) and \(l>0\);

(A3):

\(\psi(u)\leq\varpi(u)\) for all \(u\in P(\varpi,r_{2}; \omega,l_{2})\);

and if \(A: P(\varpi,r_{2}; \omega,l_{2})\rightarrow P(\varpi,r_{2}; \omega,l_{2})\) is a completely continuous operator which satisfies

(B1):

\(\{u\in\bar{P}(\varpi,d;\omega,l_{2};\psi,b):\psi (u)>b\}\neq\emptyset\), \(\psi(Au)>b\) for all \(u\in\bar{P}(\varpi,d;\omega,l_{2};\psi,b)\);

(B2):

\(\varpi(Au)< r_{1}\), \(\omega(Au)< l_{1}\) for \(u\in\bar {P}(\varpi,r_{1};\omega,l_{1})\);

(B3):

\(\psi(Au)>b\) for \(u\in \bar{P}(\varpi,r_{2};\omega,l_{2};\psi,b)\) with \(\varpi(Au)>d\),

then A has at least three different fixed points \(u_{1}\), \(u_{2}\) and \(u_{3}\) in \(\bar{P}(\varpi,r_{2};\omega,l_{2})\) with

$$\begin{aligned}& u_{1}\in P(\varpi,r_{1};\omega,l_{1}), \qquad u_{2}\in\bigl\{ u\in \bar{P}(\varpi,d;\omega,l_{2};\psi,b): \psi(u)>b\bigr\} , \\& u_{3}\in\bar{P}(\varpi,r_{2};\omega,l_{2}) \setminus \bigl(\bar{P}(\varpi ,r_{2};\omega,l_{2};\psi,b)\cup \bar{P}(\varpi,r_{1};\omega,l_{1}) \bigr). \end{aligned}$$

3 Triple positive solutions for BVP (1.1)

In this section, we discuss the existence of multiple positive solutions for boundary value problem (1.1).

Let \(C_{\tau}\triangleq \{\varphi |\varphi :[-\tau,0]\rightarrow \mathbb{R}^{+}\mbox{ is continuous}\}\). Then \(C_{\tau}\) is the space with the norm \(\|\varphi\|_{[-\tau,0]}=\max_{\theta\in[-\tau,0]}|\varphi (\theta)|\) for all \(\varphi\in C_{\tau}\). \(C(I,\mathbb{R})\) represents the Banach space of continuous functions from I to \(\mathbb{R}\) with the norm \(\|u\|_{I}=\max_{t\in I}|u(t)|\), where \(I\triangleq[0,1]\).

Now let us consider the boundary value problem as follows:

$$ \left \{\textstyle\begin{array}{l} D_{0+}^{q}u(t)+y(t)=0,\quad t\in (0,1), 2< q\leq3, \\ \alpha u(0)-\beta u'(0)=g_{1}(\int_{0}^{1}u(s)\, dA_{1}(s)),\qquad u''(0)=0, \\ \gamma u(1)+\delta u'(1)=g_{2}(\int_{0}^{1}u(s)\, dA_{2}(s)). \end{array}\displaystyle \right . $$
(3.1)

Lemma 3.1

Assume that \(A_{i}: [0,1]\rightarrow \mathbb{R}\) (\(i=1,2\)) is a function of bounded variation, \(g_{i}\in C(\mathbb{R},\mathbb{R})\), \(\rho\triangleq\alpha\gamma+\alpha\delta+\beta\gamma\neq0\) and \(y\in C([0,1])\). Then \(u\in C(I,\mathbb{R})\) is a solution of the boundary value problem (3.1) if and only if \(u(t)\) is a solution of the following integral equation:

$$\begin{aligned} u(t) =&\int_{0}^{1}G(t,s)y(s)\, ds+ \frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s) \, dA_{1}(s) \biggr) \\ &{}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \, dA_{2}(s) \biggr), \end{aligned}$$
(3.2)

where

$$ G(t,s)=\left \{\textstyle\begin{array}{l@{\quad}l} \frac{-\rho(t-s)^{q-1}+(\alpha t+\beta)[\gamma(1-s)+(q-1)\delta] (1-s)^{q-2}}{\rho\Gamma(q)}, & 0\leq s\leq t\leq1, \\ \frac{(\alpha t+\beta)[\gamma(1-s)+(q-1)\delta](1-s)^{q-2}}{\rho \Gamma(q)}, & 0\leq t\leq s\leq1. \end{array}\displaystyle \right . $$
(3.3)

Proof

Applying Lemma 2.1, Eq. (3.1) can be translated into the following equivalent integral equation:

$$ u(t)=-I_{0+}^{q}y(t)+c_{0}+c_{1}t+c_{2}t^{2}=- \frac{1}{\Gamma(q)}\int_{0}^{t}(t-s)^{q-1}y(s) \, ds+c_{0}+c_{1}t+c_{2}t^{2}. $$
(3.4)

From (3.4), we obtain

$$ u'(t)=-\int_{0}^{t} \frac{(t-s)^{q-2}}{\Gamma(q-1)}y(s)\, ds+c_{1}+2c_{2}t, u''(t)=-\int_{0}^{t} \frac{(t-s)^{q-3}}{\Gamma(q-2)}y(s)\, ds+2c_{2}. $$
(3.5)

Condition \(u''(0)=0\) gives \(c_{2}=0\). By the second boundary value condition of problem (3.1), we have

$$ \gamma c_{0}+(\gamma+\delta)c_{1}=\gamma I_{0+}^{q}y(1)+\delta I_{0+}^{q-1}y(1)+g_{2} \biggl(\int_{0}^{1}u(s)\, dA_{2}(s) \biggr). $$
(3.6)

From (3.5) and the first boundary value condition of problem (3.1), we have

$$\begin{aligned}& c_{1}=\frac{\alpha}{\rho} \biggl[\gamma I_{0+}^{q}y(1)+\delta I_{0+}^{q-1}y(1)+g_{2} \biggl(\int_{0}^{1}u(s)\,dA_{2}(s) \biggr) -\frac{\gamma}{\alpha}g_{1} \biggl(\int_{0}^{1}u(s) \,dA_{1}(s) \biggr) \biggr], \end{aligned}$$
(3.7)
$$\begin{aligned}& c_{0}=\frac{\beta}{\rho} \biggl[\gamma I_{0+}^{q}y(1)+\delta I_{0+}^{q-1}y(1)+g_{2} \biggl(\int_{0}^{1}u(s)\,dA_{2}(s) \biggr) \biggr] +\frac{\gamma+\delta}{\rho}g_{1} \biggl(\int _{0}^{1}u(s)\,dA_{1}(s) \biggr). \end{aligned}$$
(3.8)

Substituting (3.7) and (3.8) into (3.4), we get

$$\begin{aligned} u(t) =&-\frac{1}{\Gamma(q)}\int_{0}^{t}(t-s)^{q-1}y(s) \,ds +\frac{\alpha}{\rho} \biggl[\gamma I^{q}y(1)+\delta I^{q-1}y(1) +g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr) \\ &{}-\frac{\gamma}{\alpha}g_{1} \biggl(\int_{0}^{1}u(s) \,dA_{1}(s) \biggr) \biggr]t +\frac{\beta}{\rho} \biggl[\gamma I^{q}y(1)+\delta I^{q-1}y(1)+g_{2} \biggl(\int _{0}^{1}u(s)\,dA_{2}(s) \biggr) \biggr] \\ &{}+\frac{\gamma+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s) \,dA_{1}(s) \biggr) \\ =&-\int_{0}^{t}\frac{(t-s)^{q-1}}{\Gamma(q)}y(s)\,ds+ \frac{(\alpha t+\beta)\gamma}{\rho} \biggl(\int_{0}^{t}+\int _{t}^{1} \biggr)\frac {(t-s)^{q-1}}{\Gamma(q)}y(s)\,ds \\ &{}+\frac{(\alpha t+\beta)\delta}{\rho} \biggl(\int_{0}^{t}+\int _{t}^{1} \biggr)\frac {(t-s)^{q-2}}{\Gamma(q-1)}y(s)\,ds + \frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s) \,dA_{1}(s) \biggr) \\ &{}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr) \\ =&\int_{0}^{1}G(t,s)y(s)\,ds + \frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s) \,dA_{1}(s) \biggr) +\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int _{0}^{1}u(s)\,dA_{2}(s) \biggr), \end{aligned}$$

where \(G(t,s)\) is defined by (3.3). This indicates that u is a solution of (3.2). Conversely, noting that the above derivations are reversible, we assert that if u is a solution of the integral equation (3.2), then u is also the solution of BVP (3.1). The proof is complete. □

Lemma 3.2

Assume \(\alpha, \beta, \gamma, \delta\in[0,\infty)\) with \(\rho\triangleq\alpha\gamma+\alpha\delta+\beta\gamma>0\), then the function \(G(t,s)\) defined by (3.3) has the following properties:

  1. (1)

    \(G(t,s)\) is continuous on \([0,1]\times[0,1]\);

  2. (2)

    \(G(t,s)>0\) for any \(t, s\in[0,1]\);

  3. (3)

    \(G(t,s)\leq G(s,s)\) for any \(t, s\in(0,1)\);

  4. (4)

    there exists a positive number λ such that \(G(t,s)\geq\lambda G(s,s)\) for any \(t, s\in(0,1)\), where \(\lambda\triangleq\frac{4\alpha\gamma\delta[(q-2)\alpha +(q-1)\beta]\times\min\{1,\beta\}}{ [(q-1)\alpha\delta+\alpha\gamma-\beta\gamma]^{2}+4\alpha\beta \gamma[(q-1)\delta+\gamma]}<1\);

  5. (5)

    \(\vert \frac{\partial G(t,s)}{\partial t}\vert \leq \Lambda(s)\triangleq \frac{(q-1)\rho(1-s)^{q-2}+\alpha[\gamma(1-s)+(q-1)\delta ](1-s)^{q-2}}{\rho\Gamma(q)}\) for all \(t,s\in[0,1]\).

The proof is similar to the proof of Lemmas 3.2 and 3.3 in [7] or Lemma 3.2 in [8], so we omit it here.

For each \(\phi\in C_{\tau}\) and \(u\in C(I,\mathbb{R})\), we define

$$ u_{t}(s, \phi)\triangleq \left \{\textstyle\begin{array}{l@{\quad}l} u(t+s), & t+s\geq0, \\ \phi(t+s), & t+s< 0, s\in J. \end{array}\displaystyle \right . $$

Obviously, \(u_{t}(\cdot, \phi)\in C(I,\mathbb{R})\). Thus we have the following assertion.

Proposition 3.1

A function \(u\in C(I,\mathbb{R})\) is a solution of BVP (1.1) if and only if, for certain \(\phi\in C_{\tau}\), u is a solution of the following BVP:

$$ \left \{\textstyle\begin{array}{l} D_{0+}^{q}u(t)+f(t,u,u_{t}(\cdot,\phi))=0,\quad t\in[0,1], 2< q\leq 3, \\ \alpha u(0)-\beta u'(0)=g_{1}(\int_{0}^{1}u(s)\,dA_{1}(s)), \qquad u''(0)=0, \\ \gamma u(1)+\delta u'(1)=g_{2}(\int_{0}^{1}u(s)\,dA_{2}(s)). \end{array}\displaystyle \right . $$

Therefore, by Lemma 3.1, we have \(u(t)\in C(I,\mathbb{R})\) is a solution of BVP (1.1) if and only if \(u(t)\) is a solution of the integral equation as follows:

$$\begin{aligned} u(t) =&\int_{0}^{1}G(t,s)f \bigl(s,u(s),u_{s}(\cdot,\phi)\bigr)\,ds+\frac{\gamma (1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s)\,dA_{1}(s) \biggr) \\ &{}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr), \quad t\in[0,1], \end{aligned}$$
(3.9)

where \(G(t,s)\) is defined by (3.3).

In order to study the existence of solution of (1.1), we define the operator \(A_{\phi}:C(I,\mathbb{R}) \rightarrow C(I,\mathbb{R})\) as

$$\begin{aligned} (A_{\phi}u) (t) \triangleq&\int_{0}^{1}G(t,s)f \bigl(s,u(s),u_{s}(\cdot,\phi )\bigr)\,ds+\frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s)\,dA_{1}(s) \biggr) \\ &{}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr), \quad t\in[0,1]. \end{aligned}$$
(3.10)

Then solving the solutions of BVP (1.1) reduces to solving the fixed points of the operator equation \(u=A_{\phi}u\), where \(A_{\phi}\) is given by (3.10). Thus, the fixed point of operator \(A_{\phi}\) coincides with the solution of BVP (1.1).

For the sake of convenience, we introduce some assumptions as follows:

(H1):

\(\alpha, \beta, \gamma, \delta\in[0,\infty )\) with \(\rho\triangleq\alpha\gamma+\alpha\delta+\beta\gamma>0\);

(H2):

\(f\in C([0,1]\times\mathbb{R}^{+}\times\mathbb{R}^{+},\mathbb{R}^{+})\) with \(f(t,0,0)\neq0\) for all \(t\in[0,1]\);

(H3):

\(g_{i}\in C(\mathbb{R}^{+}, \mathbb{R}^{+})\) (\(i=1,2\)); \(\phi\in C([-\tau,0],\mathbb{R}^{+})\);

(H4):

\(A_{i}: [0,1]\rightarrow \mathbb{R}\) (\(i=1, 2\)) is the increasing positive function of bounded variation.

Theorem 3.1

Assume that conditions (H1)-(H4) hold. Suppose further that there exist \(\theta\in(0,\frac{1}{2})\), \(\xi\in(0,1)\), \(\sigma\in(1,+\infty)\), \(\kappa\in[3,+\infty)\), \(l_{1}\in (0,\frac{\rho}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)} ]\), \(l_{2}\in (0,\frac{\rho}{3(\alpha+\beta)\int_{0}^{1} dA_{2}(s)} ]\) and some positive constants a, b, c with \(0< a< b<\mu^{2}c\) such that the following conditions

(H5):

\(|g_{i}(u)-g_{i}(v)|\leq l_{i}|u-v|\) with \(g_{i}(0)=0\) for \(u,v\in[0,+\infty)\), \(i=1,2\);

(H6):

\(f(t,u,v)<\frac{a}{\kappa\int_{0}^{1}G(s,s)\,ds}\) for \((t,u,v)\in[0,1]\times[0,a]\times[0,a]\);

(H7):

\(f(t,u,v)\geq\frac{b\sigma}{\mu\int_{\theta }^{1-\theta}G(\xi,s)\,ds}\) for \((t,u,v)\in[\theta,1-\theta]\times [b,\frac{b}{\mu^{2}}]\times[b,\frac{b}{\mu^{2}}]\);

(H8):

\(f(t,u,v)\leq\frac{c}{\kappa\int_{0}^{1}G(s,s)\,ds}\) for \((t,u,v)\in[0,1]\times[0,c]\times[0,c]\),

have also been fulfilled, where \(\mu\triangleq\min\{\lambda,\frac{\gamma\theta+\delta}{\gamma +\delta},\frac{\alpha\theta+\beta}{\alpha+\beta}\}\), \(\lambda=\frac{4\alpha\gamma\delta[(q-2)\alpha+(q-1)\beta]\times \min\{1,\beta\}}{ [(q-1)\alpha\delta+\alpha\gamma-\beta\gamma]^{2}+4\alpha\beta \gamma[(q-1)\delta+\gamma]}<1\). Then BVP (1.1) has at least three positive solutions \(u_{1}\), \(u_{2}\) and \(u_{3}\) satisfying

$$ \|u_{1}\|_{I}< a< \|u_{3}\|_{I}, \qquad \min_{\theta\leq t\leq1-\theta}u_{3}(t)< b< \min_{\theta\leq t\leq1-\theta}u_{2}(t). $$

Proof

Define a cone \(\mathcal{P}\) in \(C(I,\mathbb{R})\) as follows:

$$ \mathcal{P}\triangleq\Bigl\{ u(t)\in C(I,\mathbb{R}): u(t)\geq0, \min _{\theta\leq t\leq1-\theta}u(t)\geq\mu\|u\|_{I}\Bigr\} . $$

Let \(\psi:\mathcal{P}\rightarrow [0,+\infty)\) be the nonnegative continuous concave functional defined by

$$ \psi(u)=\min_{\theta\leq t\leq1-\theta}u(t), \quad u\in\mathcal{P}. $$

Evidently, for each \(u\in\mathcal{P}\), we have \(\psi(u)\leq\|u\|_{I}\).

Now we will prove the assertions of Theorem 3.1 through two steps.

Step 1. Take \(\mathcal{\mathcal{P}}_{c}=\{u\in\mathcal{\mathcal {P}}: \|u\|_{I}< c\}\), \(\mathcal{\mathcal{P}}_{a}=\{u\in\mathcal {\mathcal{P}}: \|u\|_{I}< a\}\), then \(\overline{\mathcal {P}}_{c}=\{u\in\mathcal{P}: \|u\|_{I}\leq c\}\), \(\overline{\mathcal{P}}_{a}=\{u\in\mathcal{P}: \|u\|_{I}\leq a\}\). Define an operator \(A_{\phi}: \overline{\mathcal{P}}_{c}\rightarrow \overline{\mathcal{P}}\) as (3.10). Now it is necessary to show that \(A_{\phi}:\overline{\mathcal{P}}_{c} \rightarrow \overline {\mathcal {P}}_{c}\) is completely continuous. In fact, for any \(u\in\overline{\mathcal{P}}_{c}\subset\mathcal{P}\) and any \(t\in{I}\), by Lemma 3.2, we have

$$\begin{aligned} \|A_{\phi}u\|_{I} =&\max_{t\in I}\bigl\vert (A_{\phi}u) (t)\bigr\vert \\ =& \max_{t\in I}\biggl\vert \int_{0}^{1}G(t,s)f \bigl(s,u(s),u_{s}(\cdot,\phi)\bigr)\,ds +\frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s)\,dA_{1}(s) \biggr) \\ &{}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr)\biggr\vert \\ \leq& \int_{0}^{1}G(s,s)f\bigl(s,u(s),u_{s}( \cdot,\phi)\bigr)\,ds +\frac{\gamma+\delta}{\rho}g_{1} \biggl(\int _{0}^{1}u(s)\,dA_{1}(s) \biggr) \\ &{}+\frac{\alpha+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr), \end{aligned}$$

which implies that

$$\begin{aligned}& \min_{\theta\leq t\leq1-\theta}(A_{\phi}u) (t) \\& \quad = \min_{\theta\leq t\leq1-\theta}\biggl\vert \int_{0}^{1}G(t,s)f \bigl(s,u(s),u_{s}(\cdot,\phi)\bigr)\,ds +\frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s)\,dA_{1}(s) \biggr) \\& \qquad {}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr)\biggr\vert \\& \quad \geq \int_{0}^{1}\min_{0\leq t\leq 1}G(t,s)f \bigl(s,u(s),u_{s}(\cdot,\phi)\bigr)\,ds +\min_{\theta\leq t\leq 1-\theta} \biggl[\frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s) \,dA_{1}(s) \biggr) \\& \qquad {}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr) \biggr] \\& \quad \geq \lambda\int_{0}^{1}G(s,s)f \bigl(s,u(s),u_{s}(\cdot,\phi)\bigr)\,ds +\frac{\gamma\theta+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s)\,dA_{1}(s) \biggr) \\& \qquad {}+\frac{\alpha\theta +\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr) \\& \quad \geq \min \biggl\{ \lambda,\frac{\gamma\theta+\delta }{\gamma+\delta},\frac{\alpha\theta +\beta}{\alpha+\beta} \biggr\} \biggl[\int_{0}^{1}G(s,s)f\bigl(s,u(s),u_{s}( \cdot,\phi)\bigr)\,ds \\& \qquad {}+\frac{\gamma+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s) \,dA_{1}(s) \biggr) +\frac{\alpha+\beta}{\rho}g_{2} \biggl(\int _{0}^{1}u(s)\,dA_{2}(s) \biggr) \biggr] = \mu\|A_{\phi}u\|_{I}. \end{aligned}$$
(3.11)

On the other hand, when \(u\in\overline{\mathcal{P}}_{c}\), then \(\|u\|_{I}\leq c\). Noting the assumptions of κ, \(l_{1}\) and \(l_{2}\), by applying Lemma 3.2 and conditions (H3) and (H8), we obtain

$$\begin{aligned} \|A_{\phi}u\|_{I} =&\max_{t\in I} \bigl\vert (A_{\phi}u) (t)\bigr\vert \\ =& \max_{t\in I}\biggl\vert \int_{0}^{1}G(t,s)f \bigl(s,u(s),u_{s}(\cdot ,\phi)\bigr)\,ds +\frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s)\,dA_{1}(s) \biggr) \\ &{}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr)\biggr\vert \\ \leq& \int_{0}^{1}G(s,s)f\bigl(s,u(s),u_{s}( \cdot,\phi)\bigr)\,ds +\frac{\gamma+\delta}{\rho}g_{1} \biggl(\int _{0}^{1}u(s)\,dA_{1}(s) \biggr) \\ &{} + \frac{\alpha+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr) \\ \leq& \frac{c}{\kappa\int_{0}^{1}G(s,s)\,ds}\int_{0}^{1}G(s,s)\,ds+ \frac{l_{1}(\gamma+\delta)\|u\|_{I}}{\rho }\int_{0}^{1} dA_{1}(s) \\ &{}+ \frac{l_{2}(\alpha+\beta)\|u\|_{I}}{\rho}\int_{0}^{1} dA_{2}(s) \\ \leq& \frac{c}{\kappa}+\frac{\rho}{3(\gamma+\delta )\int_{0}^{1} dA_{1}(s)}\times\frac{c(\gamma+\delta)}{\rho }\int _{0}^{1} dA_{1}(s) \\ &{}+\frac{\rho}{3(\alpha+\beta)\int_{0}^{1} dA_{2}(s)} \times\frac{c(\alpha+\beta)}{\rho}\int_{0}^{1} dA_{2}(s) \\ \leq& \frac{c}{\kappa}+\frac{c}{3}+\frac{c}{3}\leq \frac {c}{3}+\frac{c}{3}+\frac{c}{3}= c. \end{aligned}$$
(3.12)

By (3.11) and (3.12), we conclude that \(A_{\phi}(\overline{\mathcal{P}}_{c}) \subset\overline{\mathcal {P}}_{c}\), that is, \(A_{\phi}:\overline{\mathcal{P}}_{c} \rightarrow \overline {\mathcal{P}}_{c}\) is well defined. Similar to the arguments of (3.12), it is easy to show that \(A_{\phi}\) maps any bounded subset of \(\mathcal{P}\) into the bounded subset of \(\mathbb{R}\). So we omit it here. Thus, according to the Arzela-Ascoli theorem, we know that \(A_{\phi}: \overline{\mathcal{P}}_{c}\rightarrow \overline {\mathcal {P}}_{c}\) is completely continuous. Similarly, one can prove that \(A_{\phi}: \overline{\mathcal{P}}_{a}\rightarrow \overline{\mathcal{P}}_{a}\) defined as (3.10) is also completely continuous.

Step 2. In the following, we will verify conditions (i)-(iii) of Lemma 2.2. In fact, when \(u\in\overline{\mathcal{P}}_{a}\), according to assumption (H6), it is similar to (3.12) that \(A_{\phi}(\overline{\mathcal{P}}_{a})\subset\overline{\mathcal {P}}_{a}\), which implies that \(\|A_{\phi}u\|_{I}\leq a\) for all \(\|u\|_{I}\leq a\). Noticing that assumption (H6) is a strict inequality, we conclude that \(\|A_{\phi}u\|_{I}< a\) for all \(\|u\|_{I}\leq a\). This indicates that condition (ii) of Lemma 2.2 holds.

Next, we show that condition (i) of Lemma 2.2 is satisfied. Clearly, \(\{u\in\mathcal {P}(\psi,b,\frac{b}{\mu^{2}}):\psi(u) >b\}\neq\emptyset\). Moreover, if \(u\in\mathcal{P}(\psi,b,\frac{b}{\mu^{2}})\), then \(\psi(u)=\frac{b}{\mu^{2}}>b\), so \(b\leq u(t)\leq\frac{b}{\mu^{2}}\) for \(t\in[\theta,1-\theta]\). Thus, for \(t\in[\theta,1-\theta]\), from condition (H7), we have

$$\begin{aligned} (A_{\phi}u) (t) =&\int_{0}^{1}G(t,s)f \bigl(s,u(s),u_{s}(\cdot,\phi)\bigr)\,ds +\frac{\gamma(1-t)+\delta}{\rho}g_{1} \biggl(\int_{0}^{1}u(s)\,dA_{1}(s) \biggr) \\ &{}+\frac{\alpha t+\beta}{\rho}g_{2} \biggl(\int_{0}^{1}u(s) \,dA_{2}(s) \biggr) \\ \geq&\int_{\theta}^{1-\theta}G(t,s)f\bigl(s,u(s),u_{s}( \cdot,\phi)\bigr)\,ds \geq\frac{b\sigma}{\mu\int_{\theta}^{1-\theta}G(\xi,s)\,ds}\int_{\theta}^{1-\theta}G(t,s) \,ds. \end{aligned}$$
(3.13)

In view of (3.13) and the definition of ψ, we get

$$\begin{aligned} \psi(A_{\phi} u)&=\min_{\theta\leq t\leq1-\theta}(A_{\phi }u) (t) \geq \mu\|A_{\phi}\|_{I}=\max_{t\in I}\bigl\vert \mu(A_{\phi}u) (t)\bigr\vert \geq\mu(A_{\phi}u) (\xi ) \\ &\geq\frac{b\sigma\mu}{\mu\int_{\theta}^{1-\theta}G(\xi ,s)\,ds}\int_{\theta}^{1-\theta}G (\xi,s ) \,ds=\sigma b>b. \end{aligned}$$

Therefore, condition (i) of Lemma 2.2 is satisfied. Finally, we show that condition (iii) of Lemma 2.2 also holds. Indeed, assume that \(u\in\overline {\mathcal{P}}_{\frac{b}{\mu^{2}}}\) with \(\|A_{\phi}\|_{I}>\frac{b}{\mu^{2}}\), then by the definition of cone \(\mathcal{P}\), we have

$$ \psi(A_{\phi} u)=\min_{\theta\leq t\leq1-\theta}(A_{\phi }u) (t) \geq \mu\|A_{\phi}\|_{I}>\mu\frac{b}{\mu^{2}}= \frac{b}{\mu}>b. $$

Hence, according to Lemma 2.2, BVP (1.1) has at least three positive solutions \(u_{1}\), \(u_{2}\), \(u_{3}\) satisfying

$$ \|u_{1}\|_{I}< a< \|u_{3}\|_{I}, \qquad \min_{\theta\leq t\leq1-\theta}u_{3}(t)< b< \min_{\theta\leq t\leq1-\theta}u_{2}(t). $$

The proof is complete. □

4 Triple positive solutions for BVP (1.2)

In this section, by employing a generalization of the Leggett-Williams fixed point theorem, we investigate the existence of at least three positive solutions for boundary value problem (1.2).

Consider the following boundary value problem:

$$ \left \{\textstyle\begin{array}{l} D_{0+}^{q}u(t)+y(t)=0,\quad t\in (0,1), 2< q\leq3, \\ \alpha u(0)-\beta u'(0)=\int_{0}^{1}h_{1}(s,u(s))\,dA_{1}(s), \qquad u''(0)=0. \\ \gamma u(1)+\delta u'(1)=\int_{0}^{1}h_{2}(s,u(s))\,dA_{2}(s). \end{array}\displaystyle \right . $$
(4.1)

Lemma 4.1

Assume that \(A_{i}: [0,1]\rightarrow \mathbb{R}\) (\(i=1,2\)) is a function of bounded variation, \(h_{i}\in C([0,1]\times\mathbb{R},\mathbb{R})\), \(\rho\triangleq\alpha\gamma+\alpha\delta+\beta\gamma\neq0\) and \(y\in C([0,1])\). Then \(u\in C(I,\mathbb{R})\) is a solution of the boundary value problem (4.1) if and only if \(u(t)\) is a solution of the following integral equation:

$$\begin{aligned} u(t) =&\int_{0}^{1}G(t,s)y(s)\,ds+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s) \\ &{}+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s), \end{aligned}$$
(4.2)

where \(G(t,s)\) is defined by (3.3).

The proof is similar to the proof of Lemma 3.1, so we omit it here.

Let \(C_{\tau}\) and \(C(I,\mathbb{R})\) be defined as in Section 3. \(\mathbb{E}\triangleq C^{1}[0,1]=\{u: u, u'\in C(I,\mathbb{R})\}\). Then \(\mathbb{E}\) is a Banach space with respect to the norm

$$ \|u\|_{C^{1}}=\max\Bigl\{ \max_{t\in I}\bigl\vert u(t) \bigr\vert , \max_{t\in I}\bigl\vert u'(t)\bigr\vert \Bigr\} , $$

where \(I\triangleq[0,1]\).

Define

$$ P\triangleq\bigl\{ u\in\mathbb{E}: u(t)\geq0, u(t) \mbox{ is concave on } [0,1] \bigr\} . $$

Clearly, P is a cone.

For each \(\phi\in C_{\tau}\) and \(u\in\mathbb{E}\), we define

$$ u_{t}(s, \phi)\triangleq \left \{\textstyle\begin{array}{l@{\quad}l} u(t+s), & t+s\geq0, \\ \phi(t+s), & t+s< 0, s\in J. \end{array}\displaystyle \right . $$

Obviously, \(u_{t}(\cdot, \phi)\in\mathbb{E}\). Thus we have the following assertion.

Proposition 4.1

A function \(u\in\mathbb{E}\) is a solution of BVP (1.2) if and only if, for certain \(\phi\in C_{\tau}\), u is a solution of the following BVP:

$$ \left \{\textstyle\begin{array}{l} D_{0+}^{q}u(t)+g(t,u,u_{t}(\cdot,\phi),u')=0,\quad t\in[0,1], 2< q\leq 3, \\ \alpha u(0)-\beta u'(0)=\int_{0}^{1}h_{1}(s,u(s))\,dA_{1}(s), \qquad u''(0)=0, \\ \gamma u(1)+\delta u'(1)=\int_{0}^{1}h_{2}(s,u(s))\,dA_{2}(s). \end{array}\displaystyle \right . $$

Consequently, by Lemma 4.1, we have \(u(t)\in\mathbb{E}\) is a solution of BVPs (1.1) if and only if \(u(t)\) is a solution of the integral equation as follows:

$$\begin{aligned} u(t) =&\int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi),u'(s)\bigr)\,ds+ \frac{\gamma (1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s) \\ &{}+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s), \quad t\in[0,1], \end{aligned}$$
(4.3)

where \(G(t,s)\) is defined by (3.3).

By way of investigating the existence of solution of (1.2), we define an operator \(A_{\phi}:P\rightarrow \mathbb{E}\) by

$$\begin{aligned} (A_{\phi}u) (t) \triangleq&\int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi ),u'(s)\bigr)\,ds+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s) \\ &{}+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s),\quad t\in[0,1]. \end{aligned}$$
(4.4)

Then solving the solutions of BVP (1.2) reduces to solving the fixed points of the operator equation \(u=A_{\phi}u\), where \(A_{\phi}\) is given by (4.4). Thus, the fixed point of operator \(A_{\phi}\) coincides with the solution of BVP (1.2).

In this section, we assume that the following conditions are satisfied:

(G1):

\(\alpha, \beta, \gamma, \delta\in[0,\infty )\) with \(\rho\triangleq\alpha\gamma+\alpha\delta+\beta\gamma>0\);

(G2):

\(g\in C([0,1]\times\mathbb{R}^{+}\times\mathbb{R}^{+}\times\mathbb {R}^{+},\mathbb{R}^{+})\) with \(g(t,0,0,0)\neq0\) for all \(t\in[0,1]\);

(G3):

\(h_{i}\in C([0,1]\times\mathbb{R}^{+}, \mathbb{R}^{+})\) (\(i=1,2\)); \(\phi\in C([-\tau,0],\mathbb{R}^{+})\);

(G4):

\(A_{i}: [0,1]\rightarrow \mathbb{R}\) (\(i=1, 2\)) is the increasing positive function of bounded variation.

Lemma 4.2

Assume that (G1)-(G4) hold. Then, for \(u\in P\), we have:

  1. (i)

    \((A_{\phi}u)(t)\) is concave on \([0,1]\);

  2. (ii)

    \((A_{\phi}u)(t)\geq0\) for \(t\in[0,1]\).

Proof

(i) By the definition of \(A_{\phi}\) and \(G(t,s)\), for \(u\in P\), we have

$$ (A_{\phi}u)''(t)=-\int _{0}^{1}\frac{(t-s)^{q-3}}{\Gamma (q-2)}g\bigl(s,u(s),u_{s}( \cdot,\phi),u'(s)\bigr)\,ds\leq0, $$
(4.5)

so \((A_{\phi}u)'(t)\) is nonincreasing. This implies that \((A_{\phi}u)(t)\) is concave.

(ii) According to the nonnegativity of \(G(t,s)\) and \(g(s,u(s),u_{s}(\cdot,\phi),u'(s))\), we can verify that \((A_{\phi}u)(t)\geq0\) for \(t\in[0,1]\).

The proof is complete. □

Lemma 4.3

Suppose that (G1)-(G4) hold. Then \(A_{\phi}: P\rightarrow P\) is a completely continuous operator.

Proof

From Lemma 4.2 it follows that \(A_{\phi} : P\rightarrow P\) is well defined. Next, we show that \(A_{\phi}\) is completely continuous. To this end, we assume that r is a positive constant and \(u\in\overline{\Omega}_{r}=\{u\in P: \|u\|_{C^{1}}\leq r\}\). Note that the continuity of \(g(t,u(t),u_{t},u'(t))\) and \(h_{i}(t,u(t))\) (\(i=1,2\)) guarantees that there exist some constants \(M_{i}>0\) (\(i=1,2,3\)) such that \(g(t,u(t),u_{t},u'(t))\leq M_{1}\), \(h_{1}(t,u(t))\leq M_{2}\) and \(h_{2}(t,u(t))\leq M_{3}\) for all \(t\in[0,1]\). Therefore, by Lemma 3.2, we have

$$\begin{aligned}& \max_{t\in[0,1]}\bigl\vert (A_{\phi}u) (t)\bigr\vert \\& \quad =\max_{t\in[0,1]}(A_{\phi}u) (t) \\& \quad = \max_{t\in[0,1]} \biggl\{ \int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi ),u'(s)\bigr)\,ds+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s) \\& \qquad {}+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr\} \\& \quad \leq M_{1}\int_{0}^{1}G(s,s)\,ds+ \frac{(\gamma+\delta)M_{2}}{\rho}\int_{0}^{1} dA_{1}(s) +\frac{(\alpha+\beta)M_{3}}{\rho}\int_{0}^{1} dA_{2}(s) \end{aligned}$$

and

$$\begin{aligned}& \max_{t\in[0,1]}\bigl\vert (A_{\phi}u)'(t) \bigr\vert \\& \quad = \max_{t\in[0,1]}\biggl\vert \int_{0}^{1} \frac{\partial G(t,s)}{\partial t}g\bigl(s,u(s),u_{s}(\cdot,\phi),u'(s) \bigr)\,ds+\frac{-\gamma}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s) \\& \qquad {}+\frac{\alpha}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s)\biggr\vert \\& \quad \leq M_{1}\int_{0}^{1}\Lambda(s) \,ds+\frac{\gamma M_{2}}{\rho}\int_{0}^{1} dA_{1}(s) +\frac{\alpha M_{3}}{\rho}\int_{0}^{1} dA_{2}(s), \end{aligned}$$

which imply that \(A_{\phi}(\overline{\Omega}_{r})\) is uniformly bounded.

Next, we shall prove that \(A_{\phi}: P\rightarrow P\) is equicontinuous. Indeed, for any \(u\in\overline{\Omega}_{r}\), \(t_{1}, t_{2}\in[0,1]\), we have

$$\begin{aligned}& \bigl\vert (A_{\phi}u) (t_{2})-(A_{\phi}u) (t_{1})\bigr\vert \\& \quad = \biggl\vert \int_{0}^{1} \bigl[G(t_{2},s)-G(t_{1},s)\bigr]g\bigl(s,u(s),u_{s}( \cdot ,\phi),u'(s)\bigr)\,ds \\& \qquad {}+\frac{-\gamma(t_{2}-t_{1})}{\rho}\int _{0}^{1}h_{1}\bigl(s,u(s)\bigr) \,dA_{1}(s)+\frac{\alpha(t_{2}-t_{1})}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s)\biggr\vert \\& \quad \leq \int_{0}^{1}\biggl\vert \frac{\partial G(\xi,s)}{\partial t}\biggr\vert |t_{2}-t_{1}|g \bigl(s,u(s),u_{s}(\cdot,\phi),u'(s)\bigr)\,ds \\& \qquad {}+ \frac{\gamma |t_{2}-t_{1}|}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha|t_{2}-t_{1}|}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \\ & \quad \leq \biggl[M_{1}\int_{0}^{1} \Lambda(s)\,ds+\frac{M_{2}\gamma}{\rho}\int_{0}^{1} dA_{1}(s) +\frac{M_{3}\alpha}{\rho}\int_{0}^{1} dA_{2}(s) \biggr]|t_{2}-t_{1}| \\ & \quad \rightarrow 0,\quad \mbox{as } t_{2}\rightarrow t_{1} \end{aligned}$$
(4.6)

and

$$\begin{aligned}& \bigl\vert (A_{\phi}u)'(t_{2})-(A_{\phi}u)'(t_{1}) \bigr\vert \\ & \quad =\biggl\vert \int_{0}^{1} \biggl[ \frac{\partial G(t_{2},s)}{\partial t}-\frac{\partial G(t_{1},s)}{\partial t} \biggr]g\bigl(s,u(s),u_{s}(\cdot, \phi),u'(s)\bigr)\,ds\biggr\vert \\ & \quad \leq \frac{M_{1}}{\Gamma(q-1)}\int_{0}^{1} \bigl\vert (t_{2}-s)^{q-2}-(t_{1}-s)^{q-2}\bigr\vert \,ds\rightarrow 0, \quad \mbox{as } t_{2}\rightarrow t_{1}. \end{aligned}$$
(4.7)

Therefore, (4.6) and (4.7) imply that \(A_{\phi}\) is equicontinuous for all \(u\in\overline{\Omega}_{r}\). By applying the Arzela-Ascoli theorem, we can see that \(A_{\phi}(\overline{\Omega}_{r})\) is relatively compact. In view of Lebesgue’s dominated convergence theorem, it is clear that \(A_{\phi}\) is a continuous operator. Hence, \(A_{\phi}: P\rightarrow P\) is a completely continuous operator. The proof is complete. □

For \(u\in P\), we define

$$\varpi(u)=\max_{t\in[0,1]}\bigl\vert u(t)\bigr\vert , \qquad \omega(u)=\max_{t\in [0,1]}\bigl\vert u'(t)\bigr\vert ,\qquad \psi(u)=\min_{t\in[\vartheta,1]}u(t), $$

where \(\vartheta\in(0,1)\). It is easy to verify that \(\varpi, \omega: P\rightarrow [0,+\infty)\) are nonnegative continuous convex functionals with \(\|u\|_{C^{1}}=\max\{\varpi(u), \omega(u)\}\). \(\psi: P\rightarrow [0,+\infty)\) is a nonnegative concave functional. We have \(\psi(u)\leq\varpi(u)\) for \(u\in P\), this means that assumptions (A1)-(A3) in Lemma 2.3 hold.

Theorem 4.1

Assume that conditions (G1)-(G4) hold. If there exist constants \(r_{1}\), r, \(r_{2}\), \(l_{1}\) and \(l_{2}\) with \(0< r_{1}< r<\frac{r}{\vartheta}\leq r_{2}\), \(0< l_{1}\leq l_{2}\). Suppose further that g, \(h_{i}\) (\(i=1,2\)) satisfy the following conditions:

(G5):

\(g(t,u,v,w)\leq\min \{\frac{r_{2}}{3\int_{0}^{1}G(s,s)\,ds}, \frac{l_{2}}{3\int_{0}^{1}\Lambda(s)\,ds} \}\) for \((t,u,v,w)\in[0,1]\times[0,r_{2}]\times[0,r_{2}]\times[-l_{2}, l_{2}]\); \(h_{1}(t,u)\leq\frac{\rho\min\{r_{2}, l_{2}\}}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)}\), \(h_{2}(t,u)\leq\frac{\rho\min\{r_{2}, l_{2}\}}{3(\alpha+\beta)\int_{0}^{1} dA_{2}(s)}\), for \((t,u)\in[0,1]\times[0,r_{2}]\);

(G6):

\(g(t,u,v,w)>\frac{r}{\lambda\int_{0}^{1}G(s,s)\,ds}\) for \((t,u,v,w)\in[\vartheta,1] \times[r,\frac{r}{\vartheta}]\times[r,\frac{r}{\vartheta}]\times [-l_{2},l_{2}]\);

(G7):

\(g(t,u,v,w)<\min \{\frac{r_{1}}{3\int_{0}^{1}G(s,s)\,ds}, \frac{l_{1}}{3\int_{0}^{1}\Lambda(s)\,ds} \}\) for \((t,u,v,w)\in[0,1]\times[0,r_{1}]\times[0,r_{1}]\times[-l_{1}, l_{1}]\); \(h_{1}(t,u)<\frac{\rho\min\{r_{1}, l_{1}\}}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)}\), \(h_{2}(t,u)<\frac{\rho\min\{r_{1}, l_{1}\}}{3(\alpha+\beta)\int_{0}^{1} dA_{2}(s)}\), for \((t,u)\in[0,1]\times[0,r_{1}]\);

(G8):

\(\frac{\min \{\lambda, \frac{\delta}{\rho }, \frac{\alpha\vartheta+\beta}{\rho} \}}{ \max \{1, \frac{\gamma+\delta}{\rho}, \frac{\alpha+\beta }{\rho} \}} >\vartheta\),

where \(\lambda\triangleq\frac{4\alpha\gamma\delta[(q-2)\alpha +(q-1)\beta]\times\min\{1,\beta\}}{ [(q-1)\alpha\delta+\alpha\gamma-\beta\gamma]^{2}+4\alpha\beta \gamma[(q-1)\delta+\gamma]}<1\), \(\Lambda(s)\triangleq \frac{(q-1)\rho(1-s)^{q-2}+\alpha[\gamma(1-s)+(q-1)\delta ](1-s)^{q-2}}{\rho\Gamma(q)}\). Then BVP (1.2) has at least three nonnegative solutions \(u_{1}\), \(u_{2}\) and \(u_{3}\) satisfying

$$\begin{aligned} \begin{aligned} &\max_{t\in[0,1]}\bigl\{ u_{1}(t)\bigr\} < r_{1}, \qquad \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{1}(t) \bigr\vert \bigr\} < l_{1}, \\ &r< \min_{t\in[\vartheta,1]}\bigl\{ u_{2}(t)\bigr\} \leq\max _{t\in[0,1]}\bigl\{ u_{2}(t)\bigr\} \leq r_{2}, \qquad \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{2}(t) \bigr\vert \bigr\} < l_{2}, \\ &\min_{t\in[\vartheta,1]}\bigl\{ u_{3}(t)\bigr\} < r,\qquad r_{1}< \max_{t\in[0,1]}\bigl\{ u_{3}(t)\bigr\} < \frac{r}{\vartheta},\qquad l_{1}< \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{3}(t)\bigr\vert \bigr\} \leq l_{2}. \end{aligned} \end{aligned}$$

Proof

The boundary value problem (1.2) has a solution \(u=u(t)\) if and only if u solves the operator equation \(A_{\phi}u=u\). Thus, we set out to verify that the operator \(A_{\phi}\) satisfies Lemma 2.3, which will prove the existence of a fixed point of \(A_{\phi}\).

We first prove that if assumption (G5) is satisfied, then \(A_{\phi}:\bar{P}(\varpi,r_{2};\omega,l_{2})\rightarrow \bar {P}(\varpi ,r_{2}; \omega,l_{2})\) defined as (4.4). In fact, let \(u\in\bar{P}(\varpi,r_{2};\omega,l_{2})\), then

$$\varpi(u)=\max_{t\in[0,1]}\bigl\vert u(t)\bigr\vert \leq r_{2}, \qquad \omega(u)=\max_{t\in[0,1]}\bigl\vert u'(t)\bigr\vert \leq l_{2}, $$

and assumption (G5) implies

$$g\bigl(t,u(t),u_{t},u'(t)\bigr)\leq\min \biggl\{ \frac{r_{2}}{3\int_{0}^{1}G(s,s)\,ds}, \frac{l_{2}}{3\int_{0}^{1}\Lambda(s)\,ds} \biggr\} , \quad t\in[0,1]. $$

For all \(u\in P\), we have \(A_{\phi}u\in P\), therefore,

$$\begin{aligned} \varpi(A_{\phi}u) =&\max_{t\in[0,1]}\bigl\vert (A_{\phi}u) (t)\bigr\vert =\max_{t\in [0,1]}(A_{\phi}u) (t) \\ =&\max_{t\in[0,1]} \biggl\{ \int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi ),u'(s)\bigr)\,ds \\ &{}+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr\} \\ \leq&\frac{r_{2}}{3\int_{0}^{1}G(s,s)\,ds}\times\int_{0}^{1}G(s,s) \,ds+\frac {\gamma+\delta}{\rho} \times\frac{\rho r_{2}}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)}\times\int_{0}^{1} dA_{1}(s) \\ &{}+\frac{\alpha+\beta}{\rho}\times\frac{\rho r_{2}}{3(\alpha +\beta)\int_{0}^{1} dA_{1}(s)}\times\int_{0}^{1} dA_{2}(s) \\ \leq&\frac{r_{2}}{3}+\frac{r_{2}}{3}+\frac{r_{2}}{3}=r_{2} \end{aligned}$$

and

$$\begin{aligned} \omega(A_{\phi}u) =&\max_{t\in[0,1]}\bigl\vert (A_{\phi}u)'(t)\bigr\vert \\ =&\max_{t\in[0,1]}\biggl\vert \int_{0}^{1} \frac{\partial G(t,s)}{\partial t}g\bigl(s,u(s),u_{s}(\cdot,\phi),u'(s) \bigr)\,ds \\ &{}+\frac{-\gamma}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s)\biggr\vert \\ \leq&\frac{l_{2}}{3\int_{0}^{1}\Lambda(s)\,ds}\times\int_{0}^{1}\Lambda (s)\,ds+\frac{\gamma}{\rho} \times\frac{\rho l_{2}}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)}\times\int _{0}^{1} dA_{1}(s) \\ &{}+\frac{\alpha}{\rho}\times\frac{\rho l_{2}}{3(\alpha+\beta )\int_{0}^{1} dA_{1}(s)}\times\int_{0}^{1} dA_{2}(s) \\ \leq&\frac{l_{2}}{3}+\frac{l_{2}}{3}+\frac{l_{2}}{3}=l_{2}. \end{aligned}$$

Thus, \(A_{\phi}u\in\bar{P}(\varpi,r_{2};\omega,l_{2})\) and \(A_{\phi} (\bar{P}(\varpi,r_{2};\omega,l_{2}) )\subset\bar {P}(\varpi,r_{2};\omega,l_{2})\). In addition, according to Lemma 4.3, we know that \(A_{\phi}:\bar{P}(\varpi,r_{2};\omega,l_{2})\rightarrow \bar{P}(\varpi,r_{2};\omega,l_{2})\) is completely continuous.

Secondly, we show that condition (B1) of Lemma 2.3 holds. We let \(u(t)=\frac{r}{\vartheta}\) for \(t\in[0,1]\). It is obvious that \(u(t)=\frac{r}{\vartheta}\in\bar{P}(\varpi,\frac{r}{\vartheta };\omega,l_{1})\) and \(\psi(u)=\frac{r}{\vartheta}>r\), and consequently

$$\biggl\{ u\in\bar{P}\biggl(\varpi,\frac{r}{\vartheta};\omega,l_{2};\psi ,r\biggr): \psi(u)>r \biggr\} \neq\emptyset. $$

For all \(u\in\bar{P}(\varpi,\frac{r}{\vartheta};\omega,l_{2};\psi,r)\), we have \(r\leq u(t)\leq\frac{r}{\vartheta}\), \(|u'(t)|\leq l_{2}\) for all \(t\in[\vartheta,1]\). Thus, by assumption (G6), we get

$$g\bigl(t,u(t),u_{t},u'(t)\bigr)>\frac{r}{\lambda\int_{0}^{1}G(s,s)\,ds}\quad \mbox{for } t\in[\vartheta,1]. $$

From the definition of the functional ψ and Lemma 3.2, we know that

$$\begin{aligned} \psi(A_{\phi}u) =&\min_{t\in[\vartheta,1]}\bigl\vert (A_{\phi}u) (t)\bigr\vert =\min_{t\in[\vartheta,1]}(A_{\phi}u) (t) \\ =& \min_{t\in[\vartheta,1]} \biggl\{ \int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi),u'(s)\bigr)\,ds \\ &{}+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr\} \\ \geq& \int_{0}^{1}\min_{0\leq t\leq 1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi),u'(s)\bigr)\,ds \\ \geq&\lambda \int_{0}^{1}G(s,s)g\bigl(s,u(s),u_{s}( \cdot,\phi),u'(s)\bigr)\,ds \\ >& \lambda\int_{0}^{1}G(s,s)\,ds\times \frac{r}{\lambda\int_{0}^{1}G(s,s)\,ds}=r. \end{aligned}$$

So, we obtain \(\psi(A_{\phi}u)>r\) for \(u\in\bar{P}(\varpi,\frac{r}{\vartheta};\omega,l_{2};\psi,r)\). Therefore, condition (B1) of Lemma 2.3 is satisfied.

Thirdly, we show that condition (B2) of Lemma 2.3 is satisfied. For all \(u\in\bar{P}(\varpi,r_{1};\omega,l_{1})\), we have \(0\leq u(t)\leq r_{1}\), \(-l_{1}\leq u'(t)\leq l_{1}\) for \(t\in[0,1]\). From assumption (G7) we obtain

$$g\bigl(t,u(t),u_{t},u'(t)\bigr)< \min \biggl\{ \frac{r_{1}}{3\int_{0}^{1}G(s,s)\,ds}, \frac{l_{1}}{3\int_{0}^{1}\Lambda(s)\,ds} \biggr\} \quad \mbox{for } t\in[0,1]. $$

Thus

$$\begin{aligned} \varpi(A_{\phi}u) =&\max_{t\in[0,1]}\bigl|(A_{\phi}u) (t)\bigr|=\max_{t\in [0,1]}(A_{\phi}u) (t) \\ =&\max_{t\in[0,1]} \biggl\{ \int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi ),u'(s)\bigr)\,ds \\ &{}+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr\} \\ < &\frac{r_{1}}{3\int_{0}^{1}G(s,s)\,ds}\times\int_{0}^{1}G(s,s) \,ds \\ &{}+\frac {\gamma+\delta}{\rho}\times\frac{\rho r_{1}}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)}\times\int_{0}^{1} dA_{1}(s) \\ &{}+\frac{\alpha+\beta}{\rho}\times\frac{\rho r_{1}}{3(\alpha +\beta)\int_{0}^{1} dA_{1}(s)}\times\int_{0}^{1} dA_{2}(s) \\ < &\frac{r_{2}}{3}+\frac{r_{2}}{3}+\frac{r_{2}}{3}=r_{1} \end{aligned}$$

and

$$\begin{aligned} \omega(A_{\phi}u) =&\max_{t\in[0,1]}\bigl|(A_{\phi}u)'(t)\bigr| \\ =&\max_{t\in[0,1]} \biggl|\int_{0}^{1} \frac{\partial G(t,s)}{\partial t}g\bigl(s,u(s),u_{s}(\cdot,\phi),u'(s) \bigr)\,ds \\ &{}+\frac{-\gamma}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr| \\ < &\frac{l_{1}}{3\int_{0}^{1}\Lambda(s)\,ds}\times\int_{0}^{1}\Lambda (s)\,ds+\frac{\gamma}{\rho} \times\frac{\rho l_{1}}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)}\times\int _{0}^{1} dA_{1}(s) \\ &{}+\frac{\alpha}{\rho}\times\frac{\rho l_{1}}{3(\alpha+\beta )\int_{0}^{1} dA_{1}(s)}\times\int_{0}^{1} dA_{2}(s) \\ < &\frac{l_{1}}{3}+\frac{l_{1}}{3}+\frac{l_{1}}{3}=l_{1}. \end{aligned}$$

We get \(A_{\phi}:\bar{P}(\varpi,r_{1};\omega,l_{1})\rightarrow \bar {P}(\varpi ,r_{1};\omega,l_{1})\), which means that (B2) in Lemma 2.3 is satisfied.

Finally, we show that condition (B3) of Lemma 2.3 holds. Indeed, according to Lemma 3.2, we have

$$\begin{aligned} \psi(A_{\phi}u) =& \min_{t\in[\vartheta,1]}\bigl\{ \bigl\vert (A_{\phi }u) (t)\bigr\vert \bigr\} =\min_{t\in[\vartheta,1]} \bigl\{ (A_{\phi}u) (t)\bigr\} \\ =& \min_{t\in[\vartheta,1]} \biggl\{ \int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi),u'(s)\bigr)\,ds \\ &{}+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr\} \\ \geq& \lambda\int_{0}^{1}G(s,s)g \bigl(s,u(s),u_{s}(\cdot,\phi),u'(s)\bigr)\,ds + \frac{\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s) \\ &{}+\frac{\alpha\vartheta +\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \\ \geq& \min \biggl\{ \lambda, \frac{\delta}{\rho}, \frac{\alpha\vartheta+\beta}{\rho} \biggr\} \biggl[\int_{0}^{1}G(s,s)g\bigl(s,u(s),u_{s}( \cdot,\phi),u'(s)\bigr)\,ds \\ &{}+\int_{0}^{1}h_{1}\bigl(s,u(s) \bigr)\,dA_{1}(s) +\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr] \end{aligned}$$
(4.8)

and

$$\begin{aligned} \varpi(A_{\phi}u) =&\max_{t\in[0,1]}\bigl\{ \bigl\vert (A_{\phi}u) (t)\bigr\vert \bigr\} =\max _{t\in[0,1]}\bigl\{ (A_{\phi}u) (t)\bigr\} \\ =&\max_{t\in[0,1]} \biggl\{ \int_{0}^{1}G(t,s)g \bigl(s,u(s),u_{s}(\cdot,\phi),u'(s)\bigr)\,ds \\ &{}+ \frac{\gamma(1-t)+\delta}{\rho}\int_{0}^{1}h_{1} \bigl(s,u(s)\bigr)\,dA_{1}(s)+\frac{\alpha t+\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr\} \\ \leq&\int_{0}^{1}G(s,s)g\bigl(s,u(s),u_{s}( \cdot,\phi),u'(s)\bigr)\,ds \\ &{}+\frac{\gamma+\delta}{\rho}\int _{0}^{1}h_{1}\bigl(s,u(s)\bigr) \,dA_{1}(s)+\frac{\alpha +\beta}{\rho}\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \\ \leq&\max \biggl\{ 1, \frac{\gamma+\delta}{\rho}, \frac{\alpha+\beta}{\rho} \biggr\} \biggl[ \int_{0}^{1}G(s,s)g\bigl(s,u(s),u_{s}( \cdot,\phi),u'(s)\bigr)\,ds \\ &{}+\int_{0}^{1}h_{1}\bigl(s,u(s) \bigr)\,dA_{1}(s) +\int_{0}^{1}h_{2} \bigl(s,u(s)\bigr)\,dA_{2}(s) \biggr]. \end{aligned}$$
(4.9)

For all \(u\in\bar{P}(\varpi,r_{2};\omega,l_{2};\psi,r)\) with \(\varpi(A_{\phi}u)>\frac{r}{\vartheta}\), in the light of (4.8), (4.9) and condition (G8), we have

$$\psi(A_{\phi}u)\geq\frac{\min \{\lambda, \frac{\delta }{\rho}, \frac{\alpha\vartheta+\beta}{\rho} \}}{ \max \{1, \frac{\gamma+\delta}{\rho}, \frac{\alpha+\beta }{\rho} \}}\times\varpi(A_{\phi}u)>r. $$

Therefore, condition (B3) of Lemma 2.3 is satisfied. So, all the conditions of Lemma 2.3 are satisfied. It follows from Lemma 2.3 and the assumption that \(g(t,0,0,0)\neq0\) on \([0,1]\) that \(A_{\phi}\) has at least three fixed points \(u_{1}\), \(u_{2}\) and \(u_{3}\) satisfying

$$\begin{aligned}& \max_{t\in[0,1]}\bigl\{ u_{1}(t)\bigr\} < r_{1}, \qquad \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{1}(t) \bigr\vert \bigr\} < l_{1},\qquad \\& r< \min_{t\in[\vartheta,1]}\bigl\{ u_{2}(t)\bigr\} \leq\max _{t\in[0,1]}\bigl\{ u_{2}(t)\bigr\} \leq r_{2}, \qquad \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{2}(t) \bigr\vert \bigr\} < l_{2}, \\& \min_{t\in[\vartheta,1]}\bigl\{ u_{3}(t)\bigr\} < r,\qquad r_{1}< \max_{t\in[0,1]}\bigl\{ u_{3}(t)\bigr\} < \frac{r}{\vartheta},\qquad l_{1}< \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{3}(t)\bigr\vert \bigr\} \leq l_{2}. \end{aligned}$$

The proof is complete. □

5 Some examples

In this section, we present some examples to illustrate our main results.

Example 5.1

Consider the boundary value problem of delayed nonlinear fractional differential equations as follows:

$$ \left \{ \textstyle\begin{array}{l} D^{\frac{5}{2}}u(t)+f(t,u,u_{t})=0, \quad t\in[0,1], \\ \frac{1}{2}u(0)-u'(0)=g_{1} (\int_{0}^{1}u(s)\,dA_{1}(s) ),\qquad u''(0)=0, \\ \frac{1}{2}u(1)+\frac{1}{2}u'(1)=g_{2} (\int_{0}^{1}u(s)\,dA_{2}(s) ), \\ u(s)=\phi(s), \quad s\in[-\tau,0], \end{array}\displaystyle \right . $$
(5.1)

where, \(q=\frac{5}{2}\), \(\alpha=\gamma=\delta=\frac{1}{2}\), \(\beta=1\), \(g_{1}(v)=g_{2}(v)=|\sin\frac{v}{4}|\), \(A_{1}(s)=A_{2}(s)=\frac{1}{3}s\), \(\tau>0\), \(\phi(s)\in C_{\tau}\) and

$$ f(t,u,v)= \left \{\textstyle\begin{array}{l@{\quad}l} \frac{\sin(\pi t)}{100}+3(u+v)^{3}, & t\in[0,1], u+v\leq 1, \\ \frac{\sin(\pi t)}{100}+\frac{u+v}{20}+10, & t\in[0,1], u+v>1. \end{array}\displaystyle \right . $$

It is easy to see that assumptions (H1)-(H4) hold and \(f(t,0,0)\neq0\) on \([0,1]\), \(g_{i}(0)=0\) (\(i=1,2\)), \(\rho=\alpha\gamma+\alpha\delta+\beta\gamma=1\), \(0< l_{1}=\frac{1}{4}<\frac{\rho}{3(\gamma+\delta)\int_{0}^{1} dA_{1}(s)}=1\) and \(0< l_{1}=\frac{1}{4}<\frac{\rho}{3(\alpha+\beta)\int_{0}^{1} dA_{2}(s)}=\frac{2}{3}\).

Take \(\theta=\frac{1}{4}\in(0,1)\), \(\xi=\frac{1}{2}\in(0,1)\), \(\sigma=2\in(1,+\infty)\), \(\kappa=4\in[3,+\infty)\), \(a=0.1\), \(b=1\) and \(c=50\) with \(0< a< b<\frac{c}{\mu^{2}}\). By a simple calculation, we obtain \(\lambda=\frac{4\alpha\gamma\delta[(q-2)\alpha+(q-1)\beta]\times \min\{1,\beta\}}{ [(q-1)\alpha\delta+\alpha\gamma-\beta\gamma]^{2}+4\alpha\beta \gamma[(q-1)\delta+\gamma]}=\frac{56}{81}<1\), \(\mu=\min \{\lambda,\frac{\gamma\theta+\delta}{\gamma +\delta},\frac{\alpha\theta+\beta}{\alpha+\beta} \} =\min \{\frac{56}{81},\frac{5}{8},\frac{3}{4} \}=\frac{5}{8}\), \(G(s,s) =\frac{(s+2)[2(1-s)^{\frac{3}{2}}+3(1-s)^{\frac{1}{2}}]}{8\Gamma (\frac{5}{2})}\), \(\int_{0}^{1}G(s,s)\,ds=\frac{116}{105\sqrt{\pi}}\approx0.623295\), \(\int_{\theta}^{1-\theta}G(\xi,s)\,ds=\int_{\frac{1}{4}}^{\frac {3}{4}}G(\frac{1}{2},s)\,ds=\frac{195\sqrt{3}-47}{480\sqrt{\pi }}\approx0.341746\) and

$$\begin{aligned}& f(t,u,v)= \frac{\sin(\pi t)}{100}+3(u+v)^{3}\leq0.034 \\& \hphantom{f(t,u,v)}< \frac{a}{\kappa\int_{0}^{1}G(s,s)\,ds}\approx0.040109\quad \mbox{for } (t,u,v) \in[0,1]\times[0,0.1]\times[0,0.1], \\& f(t,u,v)= \frac{\sin(\pi t)}{100}+\frac{u+v}{20}+10\geq10.107071 \\& \hphantom{f(t,u,v)}> \frac{b\sigma}{\mu\int_{\theta}^{1-\theta}G(\xi,s)\,ds}\approx 9.363679\quad \mbox{for } (t,u,v)\in \biggl[\frac{1}{4},\frac{3}{4} \biggr]\times \biggl[1, \frac{64}{25} \biggr]\times \biggl[1,\frac {64}{25} \biggr], \\& f(t,u,v)\leq \frac{\sin(\pi t)}{100}+\frac{u+v}{20}+10\leq15.01 \\& \hphantom{f(t,u,v)}< \frac{c}{\kappa\int_{0}^{1}G(s,s)\,ds}\approx20.054709\quad \mbox{for } (t,u,v) \in[0,1]\times[0,50]\times[0,50]. \end{aligned}$$

Thus, all the conditions of Theorem 3.1 are satisfied. According to Theorem 3.1, BVP (5.1) has at least three positive solutions \(u_{1}\), \(u_{2}\) and \(u_{3}\) such that

$$ \|u_{1}\|_{I}< 0.1< \|u_{3}\|_{I},\qquad \min_{\frac{1}{4}\leq t\leq\frac{3}{4}}u_{3}(t)< 1< \min_{\frac{1}{4}\leq t\leq\frac{3}{4}}u_{2}(t). $$

Example 5.2

Consider the following boundary value problem:

$$ \left \{ \textstyle\begin{array}{l} D^{\frac{5}{2}}u(t)+g(t,u,u_{t},u')=0, \quad t\in[0,1], \\ \frac{1}{2}u(0)-u'(0)=\int_{0}^{1}h_{1}(s,u(s))\,dA_{1}(s), \qquad u''(0)=0, \\ \frac{1}{2}u(1)+\frac{1}{2}u'(1)=\int_{0}^{1}h_{2}(s,u(s))\,dA_{2}(s), \\ u(s)=\phi(s),\quad s\in[-\tau,0], \end{array}\displaystyle \right . $$
(5.2)

where \(q=\frac{5}{2}\), \(\alpha=\gamma=\delta=\frac{1}{2}\), \(\beta=1\), \(h_{1}(t,u)=\frac{t^{2}u}{3}+ (\frac{u}{100} )^{3}t\), \(h_{2}(t,u)=\frac{t^{2}u}{5}\), \(A_{1}(s)=A_{2}(s)=\frac{1}{3}s^{2}\), \(\tau>0\), \(\phi(s)\in C_{\tau}\) and

$$ g(t,u,v,w)= \left \{\textstyle\begin{array}{l@{\quad}l} \frac{t}{50}+\frac{2}{3}(u^{3}+v^{3})+ (\frac{w}{100} )^{3}, & t\in[0,1], u\leq3, \\ \frac{t}{50}+ (\frac{w}{100} )^{3}+18, & t\in [0,1], u>3. \end{array}\displaystyle \right . $$

Clearly, assumptions (G1)-(G4) hold and \(g(t,0,0,0)\neq0\) on \([0,1]\).

Choose \(\vartheta=\frac{1}{5}\in(0,1)\), \(r=3\), \(r_{1}=\frac{1}{2}\), \(r_{2}=120\), \(l_{1}=\frac{3}{4}\) and \(l_{2}=90\). So \(0< r_{1} < r<\frac{r}{\vartheta}\) and \(0< l_{1}< l_{2}\). By calculating, we obtain \(\rho=\alpha\gamma+\alpha\delta+\beta\gamma=1\), \(\lambda=\frac{4\alpha\gamma\delta[(q-2)\alpha+(q-1)\beta]\times \min\{1,\beta\}}{ [(q-1)\alpha\delta+\alpha\gamma-\beta\gamma]^{2}+4\alpha\beta \gamma[(q-1)\delta+\gamma]}=\frac{56}{81}<1\), \(G(s,s)=\frac{(s+2)[2(1-s)^{\frac{3}{2}}+3(1-s)^{\frac {1}{2}}]}{8\Gamma(\frac{5}{2})}\), \(\int_{0}^{1}G(s,s)\,ds=\frac{116}{105\sqrt{\pi}}\approx0.623295\), \(\Lambda(s)=\frac{2(1-s)^{\frac{3}{2}}+15(1-s)^{\frac {1}{2}}}{8\Gamma(\frac{5}{2})}\), \(\int_{0}^{1}\Lambda(s)\,ds=\frac{9}{5\sqrt{\pi}}\approx1.015541\). Now, we show that conditions (G5)-(G8) are satisfied:

(G5):

\(g(t,u,v,w)\leq18.749<29.540904=\min\{ 64.175069,29.540904\}\approx\min \{\frac{r_{2}}{3\int_{0}^{1}G(s,s)\,ds}, \frac{l_{2}}{3\int_{0}^{1}\Lambda(s)\,ds} \}\) for \((t,u,v,w)\in[0,1]\times[0,120]\times[0,120]\times[-90,90]\); \(h_{1}(t,u)\leq41.728\leq90=\frac{\rho\min\{r_{2}, l_{2}\}}{3(\gamma +\delta)\int_{0}^{1} dA_{1}(s)}\), \(h_{2}(t,u)\leq24<60=\frac{\rho\min\{r_{2}, l_{2}\}}{3(\alpha+\beta )\int_{0}^{1} dA_{2}(s)}\) for \((t,u)\in[0,1]\times[0,120]\);

(G6):

\(g(t,u,v,w)\geq17.275>6.961849\approx\frac {r}{\lambda\int_{0}^{1}G(s,s)\,ds}\) for \((t,u,v,w)\in[\frac{1}{5},1] \times[3,15]\times[3,15]\times[-90,90]\);

(G7):

\(g(t,u,v,w)\leq0.186633<0.246174=\min\{ 0.267396,0246174\}\approx\min \{\frac{r_{1}}{3\int_{0}^{1}G(s,s)\,ds}, \frac{l_{1}}{3\int_{0}^{1}\Lambda(s)\,ds} \}\) for \((t,u,v,w)\in[0,1]\times[0,\frac{1}{2}]\times[0,\frac{1}{2}]\times [-\frac{3}{4},\frac{3}{4}]\); \(h_{1}(t,u)\leq0.016668<0.25=\frac{\rho\min\{r_{1}, l_{1}\}}{3(\gamma +\delta)\int_{0}^{1} dA_{1}(s)}\), \(h_{2}(t,u)\leq0.1<\frac{1}{3}=\frac{\rho\min\{r_{1}, l_{1}\}}{3(\alpha +\beta)\int_{0}^{1} dA_{2}(s)}\) for \((t,u)\in[0,1]\times[0,\frac{1}{2}]\);

(G8):

\(\frac{\min \{\lambda, \frac{\delta}{\rho }, \frac{\alpha\vartheta+\beta}{\rho} \}}{ \max \{1, \frac{\gamma+\delta}{\rho}, \frac{\alpha+\beta }{\rho} \}}=\frac{1}{3} >\frac{1}{5}=\vartheta\).

From the above, we see that all the conditions of Theorem 4.1 are satisfied. Hence, by Theorem 4.1, BVP (5.2) has at least three nonnegative solutions \(u_{1}\), \(u_{2}\) and \(u_{3}\) such that

$$\begin{aligned}& \max_{t\in[0,1]}\bigl\{ u_{1}(t)\bigr\} < \frac{1}{2}, \qquad \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{1}(t) \bigr\vert \bigr\} < \frac{3}{4}, \\& 3< \min_{t\in[\vartheta,1]}\bigl\{ u_{2}(t)\bigr\} \leq\max _{t\in[0,1]}\bigl\{ u_{2}(t)\bigr\} \leq 120, \qquad \max _{t\in[0,1]}\bigl\{ \bigl\vert u'_{2}(t)\bigr\vert \bigr\} < 90, \\& \min_{t\in[\vartheta,1]}\bigl\{ u_{3}(t)\bigr\} < 3,\qquad \frac{1}{2}< \max_{t\in[0,1]}\bigl\{ u_{3}(t)\bigr\} < 15,\qquad \frac{3}{4}< \max_{t\in[0,1]}\bigl\{ \bigl\vert u'_{3}(t)\bigr\vert \bigr\} \leq90. \end{aligned}$$

References

  1. Li, X, Chen, F, Li, X: Generalized anti-periodic boundary value problems of impulsive fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 18, 28-41 (2013)

    Article  MATH  MathSciNet  Google Scholar 

  2. Tariboon, J, Ntouyas, SK, Sudsutad, W: Positive solutions for fractional differential equations with three-point multi-term fractional integral boundary conditions. Adv. Differ. Equ. 2014, Article ID 28 (2014)

    Article  MathSciNet  Google Scholar 

  3. Zhang, X, Liu, L, Wu, Y: Multiple positive solutions of a singular fractional differential equation with negatively perturbed term. Math. Comput. Model. 55(3), 1263-1274 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  4. Chang, Y, Nieto, JJ: Some new existence results for fractional differential inclusions with boundary conditions. Math. Comput. Model. 49, 605-609 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  5. Nyamoradi, N, Bashiri, T: Multiple positive solutions for nonlinear fractional differential systems. Fract. Differ. Calc. 2(2), 119-128 (2012)

    Article  MathSciNet  Google Scholar 

  6. Zhao, Y, Sun, S, Han, Z, Li, Q: The existence of multiple positive solutions for boundary value problems of nonlinear fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 16(4), 2086-2097 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  7. Li, X, Liu, S, Jiang, W: Positive solutions for boundary value problem of nonlinear fractional functional differential equations. Appl. Math. Comput. 217(22), 9278-9285 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  8. Zhao, Y, Chen, H, Huang, L: Existence of positive solutions for nonlinear fractional functional differential equation. Comput. Math. Appl. 64(10), 3456-3467 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  9. Rehman, M, Khan, R, Asif, N: Three point boundary value problems for nonlinear fractional differential equations. Acta Math. Sci., Ser. B 31(4), 1337-1346 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  10. Zhou, Y, Jiao, F, Li, J: Existence and uniqueness for fractional neutral differential equations with infinite delay. Nonlinear Anal. TMA 71(7-8), 3249-3256 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  11. Ouyang, Z, Chen, YM, Zou, SL: Existence of positive solutions to a boundary value problem for a delayed nonlinear fractional differential system. Bound. Value Probl. 2011, Article ID 475126 (2011)

    Article  MathSciNet  Google Scholar 

  12. Zhao, KH, Gong, P: Positive solutions of m-point multi-term fractional integral BVP involving time-delay for fractional differential equations. Bound. Value Probl. 2015, Article ID 19 (2015)

    Article  MathSciNet  Google Scholar 

  13. Zhao, KH, Gong, P: Positive solutions for impulsive fractional differential equations with generalized periodic boundary value conditions. Adv. Differ. Equ. 2014, Article ID 255 (2014)

    Article  Google Scholar 

  14. Zhao, KH, Gong, P: Existence of positive solutions for a class of higher-order Caputo fractional differential equation. Qual. Theory Dyn. Syst. 14(1), 157-171 (2015)

    Article  MathSciNet  Google Scholar 

  15. Zhao, KH: Multiple positive solutions of integral BVPs for high-order nonlinear fractional differential equations with impulses and distributed delays. Dyn. Syst. 30(2), 208-223 (2015)

    Article  MathSciNet  Google Scholar 

  16. Zhao, KH: Impulsive boundary value problems for two classes of fractional differential equation with two different Caputo fractional derivatives. Mediterr. J. Math. (2015). doi:10.1007/s00009-015-0536-0

    Google Scholar 

  17. Webb, JRL: Positive solutions of some higher order nonlocal boundary value problems. Electron. J. Qual. Theory Differ. Equ. 2009, 29 (2009)

    Google Scholar 

  18. Webb, JRL, Infante, G: Nonlocal boundary value problems of arbitrary order. J. Lond. Math. Soc. 79(2), 238-258 (2009)

    MATH  MathSciNet  Google Scholar 

  19. Henderson, J, Luca, R: Existence and multiplicity of positive solutions for a system of fractional boundary value problems. Bound. Value Probl. 2014, Article ID 60 (2014)

    Article  MathSciNet  Google Scholar 

  20. Ahmad, B, Ntouyas, S, Alsaedi, A: New existence results for nonlinear fractional differential equations with three-point integral boundary conditions. Adv. Differ. Equ. 2011, Article ID 107384 (2011)

    Article  MathSciNet  Google Scholar 

  21. Feng, M, Ji, D, Ge, W: Positive solutions for a class of boundary-value problem with integral boundary conditions in Banach spaces. J. Comput. Appl. Math. 222, 351-363 (2008)

    Article  MATH  MathSciNet  Google Scholar 

  22. Jia, M, Liu, X: Three nonnegative solutions for fractional differential equations with integral boundary conditions. Comput. Math. Appl. 62, 1405-1412 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  23. Zhang, HE: Multiple positive solutions of nonlinear BVPs for differential systems involving integral conditions. Bound. Value Probl. 2014, Article ID 61 (2014)

    Article  Google Scholar 

  24. Chatthai, T, Jessada, T, Sotiris, KN: Impulsive fractional boundary-value problems with fractional integral jump conditions. Bound. Value Probl. 2014, Article ID 17 (2014)

    Article  Google Scholar 

  25. Zhao, KH, Gong, P: Positive solutions of Riemann-Stieltjes integral boundary problems for the nonlinear coupling system involving fractional-order differential. Adv. Differ. Equ. 2014, Article ID 254 (2014)

    Article  Google Scholar 

  26. Kilbas, A, Srivastava, H, Trujillo, J: Theory and Applications of Fractional Differential Equations. North-Holland Mathematics Studies, vol. 204. Elsevier, Amsterdam (2006)

    Book  MATH  Google Scholar 

  27. Podlubny, I: Fractional Differential Equations. Academic Press, New York (1993)

    Google Scholar 

  28. Legget, RW, Williams, RL: Multiple positive fixed points of nonlinear operators on ordered Banach spaces. Indiana Univ. Math. J. 28, 673-688 (1979)

    Article  MathSciNet  Google Scholar 

  29. Bai, ZB, Ge, WG: Existence of three positive solutions for some second-order boundary value problems. Comput. Math. Appl. 48, 699-707 (2004)

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

The author would like to thank the anonymous referees for their useful and valuable suggestions. This work is supported by the National Natural Sciences Foundation of Peoples Republic of China under grant (No. 11161025), Yunnan Province Natural Scientific Research Fund project (No. 2011FZ058).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kaihong Zhao.

Additional information

Competing interests

The author declares that they have no competing interests.

Author’s contributions

The author 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhao, K. Triple positive solutions for two classes of delayed nonlinear fractional FDEs with nonlinear integral boundary value conditions. Bound Value Probl 2015, 181 (2015). https://doi.org/10.1186/s13661-015-0445-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0445-y

Keywords