Skip to main content

Positive solutions of m-point multi-term fractional integral BVP involving time-delay for fractional differential equations

Abstract

In the paper, we establish sufficient conditions for the existence and multiplicity of positive solutions to a class of higher-order delayed nonlinear fractional differential equations with m-point multi-term fractional integral boundary conditions. The results are established by converting the problem into an equivalent integral equation and applying fixed point theorems of the cone expansion and compression of norm type. Our study improves the previous results in the literature. As an application, an example is also provided to illustrate our main results.

1 Introduction

The fractional differential equation has a significant role to play in many fields such as physics, chemistry, aerodynamics, electrodynamics of a 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. The fractional differential equation also serves as an excellent tool for the description of hereditary properties of various materials and processes. In consequence, the subject of fractional differential equations is gaining much importance and attention. There are a large number of papers dealing with the existence or multiplicity of solutions or positive solutions of initial or boundary value problems for some differential equations. For details and examples, see [1–23] and the references therein. In [8, 9, 16, 18–23], the authors have discussed the existence of multiple positive solutions for boundary value problems of integer or fractional differential equations.

Integral boundary conditions have various applications in applied fields such as blood flow problems, chemical engineering, thermoelasticity, underground water flow, population dynamics, and so on. For more details of nonlocal and integral boundary conditions, see [5, 10–12, 20] and the references therein. In addition, it is well known that delay arises naturally in practical systems due to the transmission of signal or the mechanical transmission. Though the theory of ordinary differential equations with delays is mature, not much has been done for fractional differential equations with delays; see [6, 10, 13, 14].

In [11], Tariboon et al. are concerned with the existence of at least one, two or three positive solutions for the boundary value problem with three-point multi-term fractional integral boundary conditions,

$$ \left \{\begin{array}{l} D^{q}_{0+} u(t) +f(t,u(t))=0,\quad 1< q\leq2, t\in[0,1], \\ u(0)=0, \qquad u(1)=\sum_{i=1}^{m}\alpha_{i}(I^{p_{i}}_{0+}u)(\eta ),\quad 0<\eta<1, \end{array} \right . $$
(1.1)

where \(f:[0,1]\times[0,\infty)\rightarrow [0,\infty)\). \(D^{q}_{0+}\) is the standard Riemann-Liouville derivative of order q. \(I^{p_{i}}_{0+}\) is the Riemann-Liouville fractional integral of order \(p_{i}>0\). \(\alpha_{i}\geq0\) (\(i=1,2,\ldots, m-2\)) are real constants.

To the best of our knowledge, no one has studied the existence of multiple positive solutions with delayed nonlinear fractional differential equations. In this article, motivated by the above-mentioned papers, we study the existence of multiple positive solutions for the following higher-order delayed nonlinear fractional differential equation with m-point multi-term fractional integral boundary conditions:

$$ \left \{\begin{array}{l} D_{0+}^{q} u(t) +f(t,u(t-\tau_{1}),u(t+\tau_{2}))=0,\quad t\in[0,1], \\ u^{(j)}(0)=0, \qquad u(1)=\sum_{i=1}^{m-2}\alpha _{i}(I^{p_{i}}_{0+}u)(\zeta_{i}), \quad 0\leq j\leq n-2, 0<\zeta_{i}<1, \end{array} \right . $$
(1.2)

subject to the following initial conditions:

$$ \left \{\begin{array}{l} au(t)-bu'(t) =\eta(t), \quad t\in[-\tau_{1},0], \\ cu(t)+du'(t) =\xi(t),\quad t\in[1,1+\tau_{2}], \end{array} \right . $$
(1.3)

where \(0\leq\tau_{1}\), \(\tau_{2}<\theta\in(0,1/2)\) are suitably small. \(D_{0^{+}}^{q}\) is the standard Riemann-Liouville derivative of order \(n-1< q\leq n\), \(n\geq3\). \(a, b, c, d\geq0\), \(\eta(t)\in C([-\tau_{1},0], [0,+\infty))\), \(\xi(t)\in C([1,1+\tau_{2}], [0,+\infty))\), \(f\in C([0,1]\times[0,+\infty)\times[0,+\infty),[0,+\infty))\). \(I^{p_{i}}_{0+}\) is the Riemann-Liouville fractional integral of order \(p_{i}\). \(\alpha_{i}\geq0\) (\(i=1,2,\ldots,m-2\)) are real constants such that \(\sum_{i=1}^{m-2}\frac{\alpha_{i}\zeta_{i}^{p_{i}+q-1}\Gamma(q)}{\Gamma (p_{i}+q)}<1\).

Remark 1.1

By (1.2) and (1.3), we have \(\eta(0)=0\) and \(\xi(1)=cu(1)+du'(1)\).

Remark 1.2

When \(1< q\leq2\), \(\tau_{1}=\tau_{2}=0\), \(\varsigma_{i}\equiv\eta\), \(j=0\), then (1.2) and (1.3) degenerate into (1.1). So our models extend (1.1).

2 Preliminaries

For the convenience of the reader, we present here the necessary definitions from fractional calculus theory. These definitions and properties can be found in the literature.

Definition 2.1

(see [24, 25])

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

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

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

Definition 2.2

(see [24, 25])

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

$$D_{0+}^{\alpha}u(t)=\frac{1}{\Gamma(n-\alpha)}\frac {d^{n}}{dt^{n}}\int _{0}^{t}(t-s)^{n-\alpha-1}u(s)\, ds, $$

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

Lemma 2.1

(see [17])

Assume that \(u\in C(0,1)\cap L(0,1)\) with a fractional derivative of order \(\alpha>0\) that belongs to \(u\in C(0,1)\cap L(0,1)\). Then

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

for some \(C_{i} \in\mathbb{R}\), \(i=1,2,\ldots,n\), where n is the smallest integer greater than or equal to α.

Lemma 2.2

(see [26])

Let E be a Banach space, \(P\subseteq E\) is a cone, and \(\Omega_{1}\), \(\Omega_{2}\) are two bounded open balls of E centered at the origin with \(0\in\Omega_{1}\) and \(\overline{\Omega}_{1}\subset\Omega_{2}\). Suppose that \(A: P\cap(\overline{\Omega}_{2}\setminus\Omega _{1})\rightarrow P\) is a completely continuous operator such that either

  1. (i)

    \(\|Au\| \leq\|u\|\), \(u\in P\cap\partial \Omega_{1}\), and \(\|Au\| \geq\|u\|\), \(u\in P\cap\partial\Omega_{2}\), or

  2. (ii)

    \(\|Au\| \geq\|u\|\), \(u\in P\cap\partial \Omega_{1}\), and \(\|Au\| \leq\|u\|\), \(u\in P\cap\partial\Omega_{2}\)

holds. Then A has at least one fixed point in \(P\cap(\overline{\Omega}_{2}\setminus\Omega_{1})\).

For simplicity, we put

$$ \Omega=1-\sum_{i=1}^{m-2} \frac{\alpha_{i}\zeta_{i}^{p_{i}+q-1}\Gamma (q)}{\Gamma(p_{i}+q)}. $$
(2.1)

Now we present the Green’s function for BVP (1.2)-(1.3).

Lemma 2.3

Let \(\sum_{i=1}^{m-2}\frac{\alpha_{i}\zeta_{i}^{p_{i}+q-1}\Gamma(q)}{\Gamma (p_{i}+q)}<1\), \(\alpha_{i}\geq0\), \(p_{i}>0\) (\(i=1,2,\ldots,m-2\)), and \(h\in C([0,1],\mathbb{R})\), \(1\leq n-1< q \leq n\). The unique solution of

$$ \left \{\begin{array}{l} D_{0+}^{q} u(t)+h(t)=0,\quad t\in[0,1], \\ u^{(j)}(0)=0,\qquad u(1)=\sum_{i=1}^{m-2}\alpha _{i}(I^{p_{i}}_{0+}u)(\zeta_{i}), \quad j=0,1,\ldots,n-2, 0<\zeta_{i}<1, \end{array} \right . $$
(2.2)

is

$$u(t)=\int_{0}^{1}G(t,s)h(s)\, ds, $$

where \(G(t,s)\) is the Green’s function given by

$$ G(t,s)=g(t,s)+\sum_{i=1}^{m-2} \frac{\alpha_{i}t^{q-1}}{\Omega \Gamma(p_{i}+q)}g_{i}(\zeta_{i},s), $$
(2.3)

where

$$ g(t,s)= \frac{1}{\Gamma(q)}\left \{\begin{array}{l@{\quad}l} t^{q-1}(1-s)^{q-1}-(t-s)^{q-1}, & 0\leq s\leq t\leq1, \\ t^{q-1}(1-s)^{q-1}, & 0\leq t\leq s\leq1 \end{array} \right . $$
(2.4)

and

$$ g_{i}(\zeta_{i},s)= \left \{\begin{array}{l@{\quad}l} \zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}-(\zeta_{i}-s)^{p_{i}+q-1}, & 0\leq s\leq\zeta_{i}<1, \\ \zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}, & 0<\zeta_{i}\leq s\leq1. \end{array} \right . $$
(2.5)

Proof

Applying Lemma 2.1, BVP (2.2) can be expressed as an equivalent integral equation

$$ u(t)=-\frac{1}{\Gamma(q)}\int_{0}^{t}(t-s)^{q-1}h(s) \, ds+C_{1}t^{q-1}+C_{2}t^{q-2}+ \cdots+C_{n}t^{q-n}, $$
(2.6)

for \(C_{i} \in\mathbb{R}\) (\(i=1,2,\ldots,n\)) \(\in\mathbb{R}\). \(u^{(j)}(0)=0\) (\(j=0,1,2,\ldots,n-2\)) implies that \(C_{2}=C_{3}=\cdots=C_{n}=0\). Taking the Riemann-Liouville fractional integral of order \(p_{i}>0\) for (2.6), we get

$$\begin{aligned} \bigl(I^{p_{i}}_{0+}u\bigr) (t)&=\int_{0}^{t} \frac{(t-s)^{p_{i}-1}}{\Gamma(p_{i})} \biggl(C_{1}s^{q-1}-\int_{0}^{s} \frac{(s-r)^{q-1}}{\Gamma(q)}h(r)\,dr \biggr)\,ds \\ &=C_{1}\int_{0}^{t}\frac{(t-s)^{p_{i}-1}s^{q-1}}{\Gamma(p_{i})} \,ds-\int_{0}^{t}\frac{(t-s)^{p_{i}-1}}{\Gamma(p_{i})}\int _{0}^{s}\frac {(s-r)^{q-1}}{\Gamma(q)}h(r)\,ds\,dr \\ &=C_{1}\frac{t^{p_{i}+q-1}\Gamma(q)}{\Gamma(p_{i}+q)}-\frac{1}{\Gamma (p_{i}+q)}\int_{0}^{t}(t-s)^{p_{i}+q-1}h(s) \,ds. \end{aligned}$$

\(u(1)=\sum_{i=1}^{m-2}\alpha_{i}(I^{p_{i}}_{0+}u)(\zeta_{i})\) yields

$$\begin{aligned}& C_{1}-\int_{0}^{1}\frac{(1-s)^{q-1}}{\Gamma(q)}h(s) \,ds \\& \quad =C_{1}\sum_{i=1}^{m-2} \frac{\alpha_{i}\zeta_{i}^{p_{i}+q-1}\Gamma(q)}{\Gamma(p_{i}+q)} -\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Gamma(p_{i}+q)}\int_{0}^{\zeta _{i}}( \zeta_{i}-s)^{p_{i}+q-1}h(s)\,ds. \end{aligned}$$

Then we have

$$ C_{1}=\frac{1}{\Omega} \Biggl[\int_{0}^{1} \frac{(1-s)^{q-1}}{\Gamma(q)}h(s)\,ds -\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Gamma(p_{i}+q)}\int_{0}^{\zeta _{i}}( \zeta_{i}-s)^{p_{i}+q-1}h(s)\,ds \Biggr]. $$

Therefore, the unique solution of BVP (2.2) is written as

$$\begin{aligned} u(t) =&-\frac{1}{\Gamma(q)}\int_{0}^{t}(t-s)^{q-1}h(s) \,ds+\frac{1}{\Omega \Gamma(q)}\int_{0}^{1}(1-s)^{q-1}t^{q-1}h(s) \,ds \\ &{}-\frac{1}{\Omega}\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Gamma (p_{i}+q)}\int_{0}^{\zeta_{i}}( \zeta_{i}-s)^{p_{i}+q-1}t^{q-1}h(s)\,ds. \end{aligned}$$

Hence, by taking into account (2.1), we have

$$\begin{aligned} u(t) =&-\frac{1}{\Gamma(q)}\int_{0}^{t}(t-s)^{q-1}h(s) \,ds+\frac{1}{\Omega \Gamma(q)}\int_{0}^{1}(1-s)^{q-1}t^{q-1}h(s) \,ds \\ &{}-\frac{1}{\Omega}\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Gamma (p_{i}+q)}\int_{0}^{\zeta_{i}}( \zeta_{i}-s)^{p_{i}+q-1}t^{q-1}h(s)\,ds \\ &{}+\frac{1}{\Gamma(q)}\int_{0}^{1}(1-s)^{q-1}t^{q-1}h(s) \,ds-\frac {1}{\Gamma(q)}\int_{0}^{1}(1-s)^{q-1}t^{q-1}h(s) \,ds \\ =&\frac{1}{\Gamma(q)}\int_{0}^{1}(1-s)^{q-1}t^{q-1}h(s) \,ds-\frac {1}{\Gamma(q)}\int_{0}^{t}(t-s)^{q-1}h(s) \,ds \\ &{}+\sum_{i=1}^{m-2}\frac{\alpha_{i}t^{q-1}}{\Omega\Gamma (p_{i}+q)} \biggl(\int_{0}^{1}\zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}h(s) \,ds -\int_{0}^{\zeta_{i}}(\zeta_{i}-s)^{p_{i}+q-1}h(s) \,ds \biggr) \\ =&\int_{0}^{1}g(t,s)h(s)\,ds+\int _{0}^{1}\sum_{i=1}^{m-2} \frac {\alpha_{i}t^{q-1}}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s)h(s)\,ds \\ =&\int_{0}^{1}G(t,s)h(s)\,ds. \end{aligned}$$

The proof is complete. □

Lemma 2.4

The Green’s function \(G(t,s)\) defined by (2.3)-(2.5) has the following properties:

  1. (A1)

    \(G(t,s)\in C([0,1]\times[0,1])\) and \(G(t,s)\geq0\), for all \(t, s \in[0,1]\).

  2. (A2)

    \(G(t,s)\leq\max_{0\leq t\leq1}G(t,s)\leq G(s)\) and \(G(t,s)\geq\min_{0\leq t\leq1}G(t,s)\geq\sigma(t)G(s)\), where

    $$ G(s)=\frac{s(1-s)^{q-1}}{\Gamma(q-1)}+\sum_{i=1}^{m-2} \frac {\alpha_{i}}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s),\qquad \sigma(t)= \frac{t^{q-1}(1-t)}{q-1}. $$
  3. (A3)

    If \(\theta\in(0,1/2)\), then \(\min_{t\in[\theta,1-\theta]}G(t,s)\geq\sigma(\theta)G(s)\).

Proof

From the expression of \(G(t,s)\), it is obvious that \(G(t,s)\in C([0,1]\times[0,1])\). To prove \(G(t,s)\geq0\), we will show that \(g(t,s), g_{i}(\zeta_{i},s)\geq0\), \(i=1,2,\ldots,m-2\), for all \(0\leq s, t \leq1\).

Let \(k_{1}(t,s)=t^{q-1}(1-s)^{q-1}-(t-s)^{q-1}\) for \(0\leq s\leq t\leq1\), then we have

$$k_{1}(t,s)=(t-ts)^{q-1}-(t-s)^{q-1} \geq(t-s)^{q-1}-(t-s)^{q-1}=0. $$

For \(0\leq t\leq s\leq1\), \(k_{2}(t,s)=(t-ts)^{q-1}\geq0\). Therefore, \(g(t,s)\geq0\) for all \(0\leq s ,t \leq1\). Now, let \(k_{3}^{i}(\zeta_{i},s)=\zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}-(\zeta_{i}-s)^{p_{i}+q-1}\) for \(0\leq s\leq\zeta_{i}<1\), then we have

$$\begin{aligned} k_{3}^{i}(\zeta_{i},s)&=\zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}- \zeta _{i}^{p_{i}+q-1} \biggl(1-\frac{s}{\zeta_{i}} \biggr)^{p_{i}+q-1} \\ &\geq\zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}- \zeta_{i}^{p_{i}+q-1}(1-s)^{p_{i}+q-1} \\ &=\zeta_{i}^{p_{i}+q-1} \bigl((1-s)^{q-1}-(1-s)^{p_{i}+q-1} \bigr)\geq0. \end{aligned}$$

For \(0< \zeta_{i}\leq s\leq1\), \(k_{4}^{i}(\zeta_{i},s)=\zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}\geq0\). Therefore \(g_{i}(\zeta_{i},s)\geq0\), \(i=1,2,\ldots,m-2\), for all \(0\leq s\leq1\).

In the following, we prove (A2). When \(s\leq t\), we have \(1-s\geq1-t\), then

$$\begin{aligned} k_{1}(t,s)&=t^{q-1}(1-s)^{q-1}-(t-s)^{q-1}=(q-1) \int_{t-s}^{t-ts}x^{q-2}\, dx \\ &\leq(q-1) (t-ts)^{q-2}\bigl((t-ts)-(t-s)\bigr) \\ &=(q-1)t^{q-2}(1-s)^{q-2}(1-t)s \\ &\leq(q-1)s(1-s)^{q-1} \end{aligned}$$

and

$$\begin{aligned} k_{1}(t,s)&=t^{q-1}(1-s)^{q-1}-(t-s)^{q-1} \\ &=(t-ts)^{q-2}(t-ts)-(t-s)^{q-2}(t-s) \\ &\geq(t-ts)^{q-2}(t-ts)-(t-ts)^{q-2}(t-s) \\ &=t^{q-2}(1-s)^{q-2}(1-t)s \\ &\geq t^{q-1}(1-t)s(1-s)^{q-1}. \end{aligned}$$

When \(t\leq s\), we derive from \(q>2\) that

$$\begin{aligned} k_{2}(t,s)&=t^{q-1}(1-s)^{q-1}\leq(q-1)t^{q-2}t(1-s)^{q-1} \\ &\leq(q-1)t^{q-2}s(1-s)^{q-1} \\ &\leq(q-1)s(1-s)^{q-1} \end{aligned}$$

and

$$ k_{2}(t,s)=t^{q-1}(1-s)^{q-1}\geq t^{q-1}(1-t)s(1-s)^{q-1}. $$

Thus,

$$ \frac{t^{q-1}(1-t)s(1-s)^{q-1}}{\Gamma(q)}\leq g(t,s)\leq\frac{(q-1)s(1-s)^{q-1}}{\Gamma(q)}=\frac {s(1-s)^{q-1}}{\Gamma(q-1)}. $$

From the above analysis, we have for \(0\leq s\leq1\),

$$\begin{aligned} G(t,s)&\leq\max_{0\leq t\leq1}G(t,s)=\max_{0\leq t\leq1} \Biggl(g(t,s)+\sum_{i=1}^{m-2} \frac{\alpha _{i}t^{q-1}}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s) \Biggr) \\ &\leq\frac{s(1-s)^{q-1}}{\Gamma(q-1)}+\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s)=G(s) \end{aligned}$$

and

$$\begin{aligned} G(t,s)&\geq\min_{0\leq t\leq1}G(t,s)=\min_{0\leq t\leq1} \Biggl(g(t,s)+\sum_{i=1}^{m-2} \frac{\alpha _{i}t^{q-1}}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s) \Biggr) \\ &\geq\frac{t^{q-1}(1-t)s(1-s)^{q-1}}{\Gamma(q)}+\sum_{i=1}^{m-2} \frac{\alpha_{i}t^{q-1}}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s) \\ &\geq\min_{0\leq t\leq1} \biggl\{ \frac {t^{q-1}(1-t)}{q-1},t^{q-1} \biggr\} \Biggl(\frac {(q-1)s(1-s)^{q-1}}{\Gamma(q)}+\sum_{i=1}^{m-2} \frac{\alpha _{i}}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s) \Biggr) \\ &=\frac{t^{q-1}(1-t)}{q-1}G(s)=\sigma(t)G(s). \end{aligned}$$

Now we prove (A3). By (A1) and (A2), we have

$$ \min_{\theta\leq t\leq 1-\theta}G(t,s)\geq\min_{\theta\leq t\leq 1-\theta}\sigma(t)G(s) $$

and

$$ \frac{\partial\sigma(t)}{\partial t}=\frac{(q-1)t^{q-2}(1-t)-t^{q-1}}{q-1} =\frac{(q-1)t^{q-2}-qt^{q-1}}{q-1}=\frac{t^{q-2}(q-1-qt)}{q-1}. $$

It is clear that \(\sigma(t)\) is increasing in \(t\in[0,\frac{q-1}{q}]\) and decreasing in \(t\in[\frac{q-1}{q},1]\), respectively. By \(q>2\), we have

$$ \frac{q-1}{q}=1-\frac{1}{q}>1-\frac{1}{2}=\frac{1}{2}> \theta, $$

thus, \(\sigma(t)\) is increasing in \(t\in[\theta,\frac{q-1}{q}]\), for \(\frac{q-1}{q}\leq1-\theta\) or \(\frac{q-1}{q}\geq1-\theta\), and we have \(\min_{\theta\leq t\leq 1-\theta}\sigma(t)=\min\{\sigma(\theta),\sigma(1-\theta)\}\). For \(0<\theta<1-\theta<1\), \(q>2\), we get

$$ \sigma(\theta)-\sigma(1-\theta)=\frac{\theta^{q-1}(1-\theta )}{q-1}-\frac{(1-\theta)^{q-1}\theta}{q-1} = \frac{\theta^{q-1}(1-\theta)}{q-1} \biggl(1- \biggl(\frac{1-\theta }{\theta} \biggr)^{q-2} \biggr)<0. $$

Therefore

$$ \min_{\theta\leq t\leq 1-\theta}G(t,s)\geq\min_{\theta\leq t\leq 1-\theta}\sigma(t)G(s)= \sigma(\theta)G(s). $$

This proof is complete. □

3 Existence of multiple positive solutions

In this section, we will consider the existence of multiple positive solutions for the BVP (1.2)-(1.3).

Let \(E= \{u(t):u(t)\in C[-\tau_{1},1+\tau_{2}]\}\) denote a real Banach space with the norm \(\|\cdot\|\) defined by \(\|u\|= \max_{-\tau_{1}\leq t\leq1+\tau_{2}} |u(t) |\). Define the cone \(P\subset E\) by \(P = \{ u\in E:u(t)\geq0 \}\). Let

$$\begin{aligned}& K = \bigl\{ u\in P: u(t)\geq\sigma(t)\|u\|, \forall t\in[0,1] \bigr\} , \end{aligned}$$
(3.1)
$$\begin{aligned}& K_{r} =\bigl\{ u\in K:\|u\|< r\bigr\} , \qquad \partial K_{r}=\bigl\{ u\in K:\|u\| =r\bigr\} . \end{aligned}$$
(3.2)

Suppose that \(u(t)\) is a solution of (1.2)-(1.3); according to Lemma 2.3 and Remark 1.1, it can be written as

$$ u(t)=\left \{\begin{array}{l@{\quad}l} u(\tau_{1}, t), & -\tau_{1}\leq t\leq0, \\ \int_{0}^{1}G(t,s)f(s,u(s-\tau_{1}),u(s+\tau_{2}))\,ds, & 0\leq t\leq 1, \\ u(\tau_{2}, t), & 1\leq t\leq1+\tau_{2}, \end{array} \right . $$

where

$$ u(\tau_{1}, t)=\left \{\begin{array}{l@{\quad}l} e^{\frac{a}{b}t} (\frac{1}{b}\int_{t}^{0}e^{-\frac{a}{b}s}\eta (s)\,ds+u(0) ), & -\tau_{1}\leq t\leq0, b\neq0, \\ \frac{\eta(t)}{a}, & -\tau_{1}\leq t\leq0, b=0 \end{array} \right . $$

and

$$ u(\tau_{2}, t)=\left \{\begin{array}{l@{\quad}l} e^{-\frac{c}{d}t} (\frac{1}{d}\int_{1}^{t}e^{\frac{c}{d}s}\xi (s)\,ds+e^{\frac{c}{d}}u(1) ), & 1\leq t\leq1+\tau_{2}, d\neq0, \\ \frac{\xi(t)}{c}, & 1\leq t\leq1+\tau_{2}, d=0. \end{array} \right . $$

Throughout this paper, we assume that \(v_{0}(t)\) is the solution of (1.2)-(1.3) with \(f\equiv0\). Clearly, \(v_{0}(t)\) can be expressed as follows:

$$ v_{0}(t)=\left \{\begin{array}{l@{\quad}l} v_{0}(\tau_{1}, t), & -\tau_{1}\leq t\leq0, \\ 0, & 0\leq t\leq1, \\ v_{0}(\tau_{2}, t), & 1\leq t\leq1+\tau_{2}, \end{array} \right . $$

where

$$ v_{0}(\tau_{1}, t)=\left \{\begin{array}{l@{\quad}l} \frac{e^{\frac{a}{b}t}}{b}\int_{t}^{0}e^{-\frac{a}{b}s}\eta(s)\,ds, & -\tau_{1}\leq t\leq0, b\neq0, \\ \frac{\eta(t)}{a}, & -\tau_{1}\leq t\leq0, b=0 \end{array} \right . $$

and

$$ v_{0}(\tau_{2}, t)=\left \{\begin{array}{l@{\quad}l} \frac{e^{-\frac{c}{d}t}}{d}\int_{1}^{t}e^{\frac{c}{d}s}\xi(s)\,ds, & 1\leq t\leq1+\tau_{2}, d\neq0, \\ \frac{\xi(t)}{c}, & 1\leq t\leq1+\tau_{2}, d=0. \end{array} \right . $$

Obviously, \(v_{0}(t)\geq0\) for each \(t\in[-\tau_{1},1+\tau_{2}]\).

Let \(u(t)\) be a solution of (1.2)-(1.3) and \(v(t)=u(t)-v_{0}(t)\). Noting that \(v(t)\equiv u(t)\) for \(0\leq t\leq1\), we have

$$ v(t)=\left \{\begin{array}{l@{\quad}l} v(\tau_{1}, t), & -\tau_{1}\leq t\leq0, \\ \int_{0}^{1}G(t,s)f(s,(v+v_{0})(s-\tau_{1}),(v+v_{0})(s+\tau_{2}))\,ds, & 0\leq t\leq1, \\ v(\tau_{2}, t), & 1\leq t\leq1+\tau_{2}, \end{array} \right . $$

where

$$\begin{aligned}& (v+v_{0}) (s-\tau_{1})=v(s-\tau_{1})+v_{0}(s- \tau_{1}), \\& (v+v_{0}) (s+\tau_{2})=v(s+\tau_{2})+v_{0}(s+ \tau_{2}), \\& v(\tau_{1}, t)=\left \{\begin{array}{l@{\quad}l} e^{\frac{a}{b}t}v(0), & -\tau_{1}\leq t\leq0, b\neq0, \\ 0, & -\tau_{1}\leq t\leq0, b=0 \end{array} \right . \end{aligned}$$

and

$$ v(\tau_{2}, t)=\left \{\begin{array}{l@{\quad}l} e^{-\frac{c}{d}(t-1)}v(1), & 1\leq t\leq1+\tau_{2}, d\neq0, \\ 0, & 1\leq t\leq1+\tau_{2}, d=0. \end{array} \right . $$

Define an operator \(A:E\rightarrow E\) as follows:

$$ (Av) (t)=\left \{\begin{array}{l@{\quad}l} (B_{1}v)(t), & -\tau_{1}\leq t\leq0, \\ \int_{0}^{1}G(t,s)f(s,(v+v_{0})(s-\tau_{1}),(v+v_{0})(s+\tau_{2}))\,ds, & 0\leq t\leq1, \\ (B_{2}v)(t), & 1\leq t\leq1+\tau_{2}, \end{array} \right . $$
(3.3)

where

$$ (B_{1}v) (t)=\left \{\begin{array}{l@{\quad}l} e^{\frac{a}{b}t}\int_{0}^{1}G(0,s)f(s,(v+v_{0})(s-\tau_{1}),(v+v_{0})(s+\tau _{2}))\,ds, & -\tau_{1}\leq t\leq0, b\neq0, \\ 0, & -\tau_{1}\leq t\leq0, b=0 \end{array} \right . $$

and

$$ (B_{2}v) (t)=\left \{\begin{array}{l@{\quad}l} e^{-\frac{c}{d}}\int_{0}^{1}f(s,(v+v_{0})(s-\tau_{1}),(v+v_{0})(s+\tau_{2})) \\ \quad {}\times G(1,s)\,ds, & 1\leq t\leq1+\tau_{2}, d\neq0, \\ 0, & 1\leq t\leq1+\tau_{2}, d=0. \end{array} \right . $$

It is easy to derive that u is a positive solution of BVP (1.2)-(1.3) if \(v=u-v_{0}\) is a nontrivial fixed point of \(A: K \rightarrow K\), where \(v_{0}\) is defined as before.

Lemma 3.1

\(A:K \rightarrow K\) defined by (3.3) is completely continuous.

Proof

For \(v\in K\), we find from Lemma 2.4 and the definition of A that \(0\leq(Av)(t)\leq(Av)(0)\), for \(t\in[-\tau_{1},0]\), and \(0\leq(Av)(t)\leq(Av)(1)\) for \(t\in[1,1+\tau_{2}]\). Thus, \(\|Av\|=\|Av\|_{[0,1]}=\max_{0\leq t\leq1}|(Av)(t)|\). It follows from Lemma 2.4 that

$$\begin{aligned} (Av) (t)&\geq\sigma(t)\int_{0}^{1}G(s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\geq\sigma(t)\|Av\|_{[0,1]}=\sigma(t)\|Av\|. \end{aligned}$$

Thus, \(A(K)\subset K\). In addition, since f is continuous, it follows that A is continuous.

Let \(Q_{1} \in K\) be bounded, that is, there exists a positive constant \(M_{1}>0\) such that \(\|v\|_{[-\tau_{1},1+\tau_{2}]}\leq M_{1}\) for all \(v \in Q_{1}\). Then \(\|v+v_{0}\|_{[-\tau_{1},1+\tau_{2}]}\leq M_{1}+M_{0}\triangleq M_{2}\) for \(v\in Q_{1}\), where \(v_{0}\) is defined as before. Define a set \(Q_{2}\subset E\) as follows:

$$ Q_{2}=\bigl\{ \psi\in E:\Vert \psi \Vert _{[-\tau_{1},1+\tau_{2}]}\leq M_{2}\bigr\} . $$

Hence, \(\max_{-\tau_{1}\leq t\leq1+\tau_{2}}|(v+v_{0})(t)|\leq \|v+v_{0}\|_{[-\tau_{1},1+\tau_{2}]}\leq \|v\|_{[-\tau_{1},1+\tau_{2}]}+\|v_{0}\|_{[-\tau_{1},1+\tau_{2}]}\leq M_{2}\). Noting that f is continuous on \([0,1]\times[0,M_{2}]\times[0,M_{2}]\), there exists a constant \(M_{3}>0\) such that on \([0,1]\times[0,M_{2}]\times[0,M_{2}]\),

$$ \bigl\vert f\bigl(t,(v+v_{0}) (t-\tau_{1}),(v+v_{0}) (t+\tau_{2})\bigr)\bigr\vert \leq M_{3}. $$

Therefore,

$$ \|Av\|_{[-\tau_{1},1+\tau_{2}]}=\|Av\|_{[0,1]}\leq M_{3}\int _{0}^{1}G(s)\,ds. $$

Hence, \(A (Q_{1})\) is bounded.

Finally, we show the operator A is equicontinuous. For \(v\in Q_{1}\), we have

$$\begin{aligned}& (Av)'(t)=\int_{0}^{1} \frac{\partial}{\partial t}G(t,s)f\bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+\tau_{2})\bigr)\,ds, \quad t\in [0,1], \\& (Av)'(t)=\frac{a}{b}e^{\frac{a}{b}t}\int _{0}^{1}G(0,s)f\bigl(s,(v+v_{0}) (s- \tau_{1}),(v+v_{0}) (s+\tau_{2})\bigr)\,ds,\quad t \in[-\tau_{1},0] \end{aligned}$$

and

$$ (Av)'(t)=-\frac{c}{d}e^{-\frac{c}{d}(t-1)}\int _{0}^{1}G(1,s)f\bigl(s,(v+v_{0}) (s- \tau_{1}),(v+v_{0}) (s+\tau_{2})\bigr)\,ds,\quad t \in[1,1+\tau_{2}]. $$

In the light of \(f\leq M_{3}\) and

$$\begin{aligned} \frac{\partial}{\partial t}G(t,s)&\leq\frac{(q-1)t^{q-2}(1-s)^{q-1}+(q-1)(t-s)^{q-2}}{\Gamma (q)}+\sum _{i=1}^{m-2}\frac{\alpha_{i}(q-1)t^{q-2}}{\Omega \Gamma(p_{i}+q)}g_{i}( \zeta_{i},s) \\ &\leq\frac{1}{\Gamma(q-1)}+\sum_{i=1}^{m-2} \frac{\alpha _{i}(q-1)}{\Omega\Gamma(p_{i}+q)}g_{i}(\zeta_{i},s)\triangleq M', \end{aligned}$$

we have \(\|(Av)'\|\leq M\) for some positive constant M. Thus, for \(-\tau_{1}\leq t_{1}< t_{2}\leq1+\tau_{2}\), we have

$$ \bigl\Vert (Av) (t_{2})-(Av) (t_{1})\bigr\Vert \leq\int _{t_{1}}^{t_{2}}\bigl\vert (Av)'(s)\bigr\vert \,ds\leq M(t_{2}-t_{1})\rightarrow 0, \quad \mbox{as }t_{1}\rightarrow t_{2}. $$

Therefore, for any \(\epsilon>0\), there exists \(\delta=\delta(\epsilon)>0\) which is independent of \(t_{1}\), \(t_{2}\), and v such that \(\|(Av)(t_{2})-(Av)(t_{1})\|\leq\epsilon\), whenever \(|t_{2}-t_{1}|\leq\delta\). Thus, \(A(Q_{1})\) is equicontinuous. In view of the Ascoli-Arzela theorem, we can easily see that \(A:K \rightarrow K\) is a completely continuous operator. The proof is complete. □

Further we make the following assumptions for \(f(t,x,y)\):

(H1):

There exists a constant \(r_{1}>0\) such that \(0\leq x\leq r_{1}+\|v_{0}\|_{[-\tau_{1},0]}\), \(0\leq y\leq r_{1}+\|v_{0}\|_{[1,1+\tau_{2}]}\), and \(0\leq t\leq1\) implies \(f(t,x,y)<\rho_{1} r_{1}\), where \(\rho_{1}=\frac{1}{\int_{0}^{1}G(s)\,ds}\).

(H2):

There exists a constant \(r_{2}>0\) such that \(\sigma _{1}r_{2}\leq x\leq r_{2}\), \(\sigma_{2}r_{2}\leq y\leq r_{2}\), and \(\theta\leq t\leq1-\theta\) implies \(f(t,x,y)>\rho_{2} r_{2}\), where \(\rho_{2}=\frac{1}{\sigma(\theta)\int_{\theta}^{1-\theta}G(s)\,ds}\), \(\sigma_{1}=\min_{\theta-\tau_{1}\leq t\leq1-\theta-\tau_{1} }|\sigma(t)|\), \(\sigma_{2}=\min_{\theta+\tau_{2}\leq t\leq1-\theta+\tau_{2} }|\sigma(t)|\).

(H3):

\(f_{\infty}=\liminf_{x+y\rightarrow +\infty}\min_{t\in[0,1]}\frac{f(t,x,y)}{x+y}=\infty\).

(H4):

\(f^{\infty}=\limsup_{x+y \rightarrow +\infty}\max_{t\in[0,1]}\frac{f(t,x,y)}{x+y}=0\).

(H5):

\(f_{0}=\liminf_{x+y\rightarrow 0}\min_{t\in[0,1]}\frac{f(t,x,y)}{x+y}=\infty\).

(H6):

\(f^{0}=\limsup_{x+y \rightarrow 0}\max_{t\in[0,1]}\frac{f(t,x,y)}{x+y}=0\).

Theorem 3.1

Assume that (H1), (H2), and (H3) are satisfied. If \(r_{1}>r_{2}>0\), then BVP (1.2)-(1.3) has at least two positive solutions \(u_{1}\) and \(u_{2}\) such that

$$ 0< r_{2}<\|u_{1}\|_{[0,1]}<r_{1}< \|u_{2}\|_{[0,1]}. $$

Proof

Let \(A:K\rightarrow K\) be the cone preserving completely continuous that is defined by (3.3). Let \(\Phi_{r_{1}}=\{v\in E:\|v\|< r_{1}\}\), then for any \(v\in K\cap\partial\Phi_{r_{1}}\), we get

$$\begin{aligned}& 0\leq(v+v_{0}) (t-\tau_{1})\leq \|v\|+\|v_{0} \|_{[-\tau_{1},1]}=r_{1}+\|v_{0}\|_{[-\tau_{1},0]},\quad 0\leq t \leq1, \\& 0\leq(v+v_{0}) (t+\tau_{2})\leq \|v\|+\|v_{0} \|_{[0,1+\tau_{2}]}=r_{1}+\|v_{0}\|_{[1,1+\tau_{2}]},\quad 0\leq t \leq1. \end{aligned}$$

Thus, from (A3) and (H1), we have

$$\begin{aligned} \|Av\| &=\|Av\|_{[0,1]}=\max_{t\in[0,1]}(Av) (t) \\ & =\max_{t\in[0,1]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\leq\int_{0}^{1}G(s)f\bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+\tau_{2})\bigr)\,ds \\ &< \rho_{1}r_{1}\int_{0}^{1}G(s) \,ds=r_{1}=\|v\|. \end{aligned}$$

Therefore,

$$ \|Av\|< \|v\|, \quad v\in K\cap\partial\Phi_{r_{1}}. $$
(3.4)

Let \(\Phi_{r_{2}}=\{v\in E:\|v\|< r_{2}\}\), then for any \(v\in K\cap\partial\Phi_{r_{2}}\), we have

$$\begin{aligned}& r_{2}=\|v\|\geq v(t-\tau_{1})\geq \sigma(t- \tau_{1})\|v\|\geq\sigma_{1}r_{2},\quad \theta\leq t \leq 1-\theta, \\& r_{2}=\|v\|\geq v(t+\tau_{2})\geq \sigma(t+ \tau_{2})\|v\|\geq\sigma_{2}r_{2},\quad \theta\leq t \leq 1-\theta. \end{aligned}$$

Thus, from (H2) and (A3) of Lemma 2.4, we get

$$\begin{aligned} \|Av\| &=\|Av\|_{[0,1]}=\max_{t\in[0,1]}(Av) (t) \\ & =\max_{t\in[0,1]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ & \geq\min_{t\in[\theta,1-\theta]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\geq\sigma(\theta)\int_{0}^{1}G(s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\geq\sigma(\theta)\int_{\theta}^{1-\theta }G(s)f \bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &=\sigma(\theta)\int_{\theta}^{1-\theta}G(s)f\bigl(s,v(s- \tau _{1}),v(s+\tau_{2})\bigr)\,ds \\ &> \rho_{2}r_{2}\sigma(\theta)\int_{\theta}^{1-\theta}G(s) \,ds=r_{2}=\|v\|. \end{aligned}$$

So

$$ \|Av\|>\|v\|, \quad v\in K\cap\partial\Phi_{r_{2}}. $$
(3.5)

Choose \(L>0\) such that

$$ (\sigma_{1}+\sigma_{2})L\sigma(\theta)\int _{\theta}^{1-\theta}G(s)\,ds>1. $$
(3.6)

From (H3), there exists \(R_{1}>0\) such that

$$ f(t,x,y)\geq L(x+y), \quad x, y\geq0, x+y\geq R_{1}, 0 \leq t\leq1. $$
(3.7)

Choose

$$ R_{0}>\max \biggl\{ r_{1},\frac{R_{1}}{2\min\{\sigma_{1},\sigma_{2}\}} \biggr\} . $$

Let \(\Phi_{R}=\{v\in E:\|v\|< R,R\geq R_{0}\}\), then for any \(v\in K\cap\partial\Phi_{R}\), we have

$$\begin{aligned}& v(t-\tau_{1})\geq\sigma(t-\tau_{1})\|v\|\geq \sigma_{1}R\geq\min\{ \sigma_{1},\sigma_{2} \}R_{0}>\frac{1}{2}R_{1},\quad \theta\leq t\leq1- \theta, \\& v(t+\tau_{2})\geq\sigma(t+\tau_{2})\|v\|\geq \sigma_{2}R\geq\min\{ \sigma_{1},\sigma_{2} \}R_{0}>\frac{1}{2}R_{1},\quad \theta\leq t\leq1- \theta. \end{aligned}$$

Then, from (3.6) and (3.7) we have

$$\begin{aligned} \|Av\| &=\|Av\|_{[0,1]}=\max_{t\in[0,1]}(Av) (t) \\ & =\max_{t\in[0,1]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ & \geq\min_{t\in[\theta,1-\theta]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\geq\sigma(\theta)\int_{0}^{1}G(s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\geq\sigma(\theta)\int_{\theta}^{1-\theta }G(s)f \bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &=\sigma(\theta)\int_{\theta}^{1-\theta}G(s)f\bigl(s,v(s- \tau _{1}),v(s+\tau_{2})\bigr)\,ds \\ &\geq L\sigma(\theta)\int_{\theta}^{1-\theta}G(s)\bigl[v(s- \tau _{1})+v(s+\tau_{2})\bigr]\,ds \\ &\geq(\sigma_{1}+\sigma_{2})RL\sigma(\theta)\int _{\theta }^{1-\theta}G(s)\,ds>R=\|v\|. \end{aligned}$$

Therefore,

$$ \|Av\|> \|v\|, \quad v\in K\cap\partial\Phi_{R}. $$
(3.8)

Applying Lemma 2.2 to (3.4) and (3.5) yields the result that A has a fixed point \(v_{1}\in K\cap(\overline{\Phi}_{r_{1}}\setminus\Phi_{r_{2}})\) with \(v_{1}(t)\geq \sigma_{1} \|u\|>0\), \(t\in[0,1]\). Similarly, Lemma 2.2 associated with (3.4) and (3.8) shows that A has another fixed point \(v_{2}\in K\cap(\overline{\Phi}_{R}\setminus \Phi_{r_{1}})\) with \(v_{2}(t)\geq\sigma_{2} \|u\|>0\), \(t\in[0,1]\), which means that \(u_{1}(t)=v_{1}(t)+v_{0}(t)\) and \(u_{2}(t)=v_{2}(t)+v_{0}(t)\) are two positive solutions of BVP (1.2)-(1.3). Since

$$ \|u_{i}\|_{[0,1]}=\|v_{i}+v_{0} \|_{[0,1]}=\|v_{i}\|_{[0,1]}=\|Av_{i} \|_{[0,1]}=\| Av_{i}\|=\|v_{i}\|,\quad i=1,2, $$

it follows that \(u_{1}(t)\) and \(u_{2}(t)\) satisfy

$$ 0< r_{2}<\|u_{1}\|_{[0,1]}=\|v_{1} \|<r_{1}<\|v_{2}\|=\|u_{2}\|_{[0,1]}. $$

The proof is complete. □

Theorem 3.2

Assume that (H1), (H3), and (H5) are satisfied. There exist constants \(R>r_{1}>r>0\), then BVP (1.2)-(1.3) has at least two positive solutions \(u_{1}\) and \(u_{2}\) such that

$$ 0< r<\|u_{1}\|_{[0,1]}<r_{1}<\|u_{2} \|_{[0,1]}. $$

Proof

Choose \(0< r<r_{1}<R\), let \(\Phi_{r}=\{v\in E:\|v\|< r,r<r_{1}\}\). For the same \(L>0\) satisfying (3.6), (H5) implies that

$$ f(t,x,y)\geq L(x+y), \quad x, y\geq0, x+y\leq2r, 0\leq t\leq1. $$
(3.9)

Then for any \(v\in K\cap\partial\Phi_{r}\), we have

$$\begin{aligned}& r=\|v\|\geq v(t-\tau_{1})\geq\sigma(t-\tau_{1})\|v\|\geq \sigma _{1}r,\quad \theta\leq t\leq1-\theta, \\& r=\|v\|\geq v(t+\tau_{2})\geq\sigma(t+\tau_{2})\|v\|\geq \sigma _{2}r,\quad \theta\leq t\leq1-\theta. \end{aligned}$$

Then, from (3.6) and (3.9) we have

$$\begin{aligned} \|Av\| &=\|Av\|_{[0,1]}=\max_{t\in[0,1]}(Av) (t) \\ & =\max_{t\in[0,1]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\geq\sigma(\theta)\int_{0}^{1}G(s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\geq\sigma(\theta)\int_{\theta}^{1-\theta }G(s)f \bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &= \sigma(\theta)\int_{\theta}^{1-\theta}G(s)f\bigl(s,v(s- \tau _{1}),v(s+\tau_{2})\bigr)\,ds \\ &\geq L\sigma(\theta)\int_{\theta}^{1-\theta}G(s)\bigl[v(s- \tau _{1})+v(s+\tau_{2})\bigr]\,ds \\ &\geq(\sigma_{1}+\sigma_{2})rL\sigma(\theta)\int _{\theta }^{1-\theta}G(s)\,ds>r=\|v\|. \end{aligned}$$

Therefore,

$$ \|Av\|> \|v\|, \quad v\in K\cap\partial\Phi_{r}. $$
(3.10)

Applying Lemma 2.2 to (3.4) and (3.10) yields that A has a fixed point \(v_{1}\in K\cap(\overline{\Phi}_{r_{1}}\setminus\Phi_{r})\) with \(v_{1}(t)\geq \sigma_{1} \|u\|>0\), \(t\in[0,1]\). Similarly, Lemma 2.2 associated with (3.4) and (3.8) yields the result that A has another fixed point \(v_{2}\in K\cap(\overline{\Phi}_{R}\setminus \Phi_{r_{1}})\) with \(v_{2}(t)\geq\sigma_{2} \|u\|>0\), \(t\in[0,1]\). This means that \(u_{1}(t)=v_{1}(t)+v_{0}(t)\) and \(u_{2}(t)=v_{2}(t)+v_{0}(t)\) are two positive solutions of BVP (1.2)-(1.3). Since

$$ \|u_{i}\|_{[0,1]}=\|v_{i}+v_{0} \|_{[0,1]}=\|v_{i}\|_{[0,1]}=\|Av_{i} \|_{[0,1]}=\| Av_{i}\|=\|v_{i}\|,\quad i=1,2, $$

it follows that \(u_{1}(t)\) and \(u_{2}(t)\) satisfy

$$ 0< r<\|u_{1}\|_{[0,1]}=\|v_{1}\|<r_{1}< \|v_{2}\|=\|u_{2}\|_{[0,1]}. $$

The proof is complete. □

Theorem 3.3

Assume that (H2), (H4), and (H6) are satisfied. There exist constants \(R>r_{2}>r>0\), then BVP (1.2)-(1.3) has at least two positive solutions \(u_{1}\) and \(u_{2}\) such that

$$ 0< r<\|u_{1}\|_{[0,1]}<r_{2}<\|u_{2} \|_{[0,1]}. $$

Proof

Choose \(0< r<r_{2}<R\). By (H4), for any \(0<\varepsilon<\frac{1}{2\int_{0}^{1}G(s)\,ds}\), there exists \(R'>0\) such that

$$ f(t,x,y)\leq\varepsilon(x+y),\quad x, y\geq0, x+y\geq R', 0\leq t \leq1. $$

Putting

$$ C\triangleq\max_{0\leq t\leq1}\max_{0\leq x,y,x+y\leq R'}\bigl\vert f(t,x,y)\bigr\vert +1, $$

then

$$ f(t,x,y)\leq\varepsilon(x+y)+C,\quad x, y\geq0, 0\leq t\leq1. $$

Choose

$$ R_{0}>\bigl(C+2\varepsilon\|v_{0}\|\bigr)\frac{\int_{0}^{1}G(s)\,ds}{1-2\varepsilon\int_{0}^{1}G(s)\,ds}. $$

Let \(\Phi_{R}=\{v\in E: \|v\|< R, R\geq\max\{r_{2},R_{0}\}\}\). Then, for any \(v\in K\cap\partial\Phi_{R}\), we have

$$\begin{aligned} \|Av\| &=\|Av\|_{[0,1]}=\max_{t\in[0,1]}(Av) (t) \\ & =\max_{t\in[0,1]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\leq\int_{0}^{1}G(s)f\bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+\tau_{2})\bigr)\,ds \\ &\leq\int_{0}^{1}G(s) \bigl[\varepsilon \bigl((v+v_{0}) (s-\tau _{1})+(v+v_{0}) (s+ \tau_{2}) \bigr)+C \bigr]\,ds \\ &\leq\int_{0}^{1}G(s) \bigl[2\varepsilon \bigl(\|v\|+ \|v_{0}\| \bigr)+C \bigr]\,ds \\ &=2\varepsilon R\int_{0}^{1}G(s)\,ds+\bigl(C+2 \varepsilon\|v_{0}\|\bigr)\int_{0}^{1}G(s)\,ds < R=\|v\|. \end{aligned}$$

Therefore,

$$ \|Av\|< \|v\|, \quad v\in K\cap\partial\Phi_{R}. $$
(3.11)

Let \(\Phi_{r}=\{v\in E:\|v\|< r,r<r_{2}\}\). By (H6), for any \(0<\varepsilon<\frac{1}{2\int_{0}^{1}G(s)\,ds}\), there exists \(0< r'<r\) such that

$$ f(t,x,y)\leq\varepsilon(x+y),\quad x, y\geq0, x+y< r', 0\leq t\leq1. $$

Putting

$$ C'\triangleq\max_{0\leq t\leq1}\max_{0\leq x,y,r'\leq x+y\leq 2r} \bigl\vert f(t,x,y)\bigr\vert +1, $$

then

$$ f(t,x,y)\leq\varepsilon(x+y)+C', \quad x, y\in \partial \Phi_{r}, 0\leq t\leq1. $$

Choose

$$ r>\bigl(C'+2\varepsilon \Vert v_{0}\Vert \bigr) \frac{\int_{0}^{1}G(s)\,ds}{1-2\varepsilon\int_{0}^{1}G(s)\,ds}. $$

Then, for any \(v\in K\cap\partial\Phi_{r}\), we have

$$\begin{aligned} \|Av\| &=\|Av\|_{[0,1]}=\max_{t\in[0,1]}(Av) (t) \\ & =\max_{t\in[0,1]}\int_{0}^{1}G(t,s)f \bigl(s,(v+v_{0}) (s-\tau _{1}),(v+v_{0}) (s+ \tau_{2})\bigr)\,ds \\ &\leq\int_{0}^{1}G(s)f\bigl(s,(v+v_{0}) (s-\tau_{1}),(v+v_{0}) (s+\tau_{2})\bigr)\,ds \\ &\leq\int_{0}^{1}G(s) \bigl[\varepsilon \bigl((v+v_{0}) (s-\tau _{1})+(v+v_{0}) (s+ \tau_{2}) \bigr)+C' \bigr]\,ds \\ &\leq\int_{0}^{1}G(s) \bigl[2\varepsilon \bigl(\|v\|+ \|v_{0}\| \bigr)+C' \bigr]\,ds \\ &=2\varepsilon r\int_{0}^{1}G(s)\,ds+ \bigl(C'+2\varepsilon\|v_{0}\|\bigr)\int _{0}^{1}G(s)\,ds < r=\|v\|. \end{aligned}$$

So,

$$ \|Av\|< \|v\|, \quad v\in K\cap\partial\Phi_{r}. $$
(3.12)

Applying Lemma 2.2 to (3.5) and (3.12) yields the result that A has a fixed point \(v_{1}\in K\cap(\overline{\Phi}_{r_{1}}\setminus\Phi_{r})\) with \(v_{1}(t)\geq \sigma_{1} \|u\|>0\), \(t\in[0,1]\). Similarly, from Lemma 2.2 associated with (3.5) and (3.11) one derives that A has another fixed point \(v_{2}\in K\cap(\overline{\Phi}_{R}\setminus \Phi_{r_{1}})\) with \(v_{2}(t)\geq\sigma_{2} \|u\|>0\), \(t\in[0,1]\). This means that \(u_{1}(t)=v_{1}(t)+v_{0}(t)\) and \(u_{2}(t)=v_{2}(t)+v_{0}(t)\) are two positive solutions of BVP (1.2)-(1.3). Since

$$ \|u_{i}\|_{[0,1]}=\|v_{i}+v_{0} \|_{[0,1]}=\|v_{i}\|_{[0,1]}=\|Av_{i} \|_{[0,1]}=\| Av_{i}\|=\|v_{i}\|,\quad i=1,2, $$

it follows that \(u_{1}(t)\) and \(u_{2}(t)\) satisfy

$$ 0< r<\|u_{1}\|_{[0,1]}=\|v_{1}\|<r_{2}< \|v_{2}\|=\|u_{2}\|_{[0,1]}. $$

This completes the proof. □

We account for the control functions

$$\begin{aligned}& \varphi(r)=\max\bigl\{ f(t,x,y):(t,x,y)\in[0,1]\times[0,r+r_{0}]\times \bigl[0,r+r'_{0}\bigr]\bigr\} , \\& \psi(r)=\min\bigl\{ f(t,x,y):(t,x,y)\in[\theta,1-\theta]\times[\sigma _{1}r,r]\times[\sigma_{2}r,r]\bigr\} , \end{aligned}$$

where \(r_{0}=\|v_{0}\|_{[-\tau_{1},0]}\), \(r'_{0}=\|v_{0}\|_{[1,1+\tau_{2}]}\).

Theorem 3.4

Suppose that there exist two positive numbers \(\xi_{2}<\xi_{1}\) such that one of the following conditions is satisfied:

(B1):

\(\varphi(\xi_{1})< \rho_{1}\xi_{1}\), \(\psi(\xi _{2})>\rho_{2}\xi_{2}\).

(B2):

\(\psi(\xi_{1})>\rho_{2}\xi_{1}\), \(\varphi(\xi_{2})<\rho _{1}\xi_{2}\).

Then BVP (1.2)-(1.3) has at least one positive solution \(u\in K\) such that

$$ \xi_{2}<\|u\|_{[0,1]}<\xi_{1}. $$

Proof

Because of the similarity of the proof, we prove only this theorem under condition (B1). By assumption (B1), we have

$$ f(t,x,y)\leq\varphi(\xi_{1})< \rho_{1}\xi_{1}, \qquad f(t,x,y)\geq\psi(\xi_{2})>\rho_{2}\xi_{2}, $$

which are the assumptions (H1) and (H2). By Theorem 3.1, we find that A has a fixed point \(v\in K\cap(\overline{\Phi}_{\xi_{1}}\setminus\Phi_{\xi_{2}})\), which means that (1.2)-(1.3) has at least one positive solution u and \(\xi_{2}<\|u\|_{[0,1]}<\xi_{1}\). This completes the proof. □

Similarly, we can obtain the existence of multiple positive solutions for BVP (1.2)-(1.3).

Theorem 3.5

Suppose that there exist three positive numbers \(\xi_{3}<\xi_{2}<\xi_{1}\) such that one of the following conditions is satisfied:

(B3):

\(\varphi(\xi_{1})<\rho_{1}\xi_{1}\), \(\psi(\xi_{2})>\rho _{2}\xi_{2}\), \(\varphi(\xi_{3})<\rho_{1}\xi_{3} \).

(B4):

\(\psi(\xi_{1})>\rho_{2}\xi_{1}\), \(\varphi(\xi_{2})<\rho _{1}\xi_{2}\), \(\psi(\xi_{3})>\rho_{2}\xi_{3} \).

Then BVP (1.2)-(1.3) has at least two positive solutions \(u_{1}, u_{2}\in K\) such that

$$ \xi_{3}<\|u_{1}\|_{[0,1]}<\xi_{2}< \|u_{2}\|_{[0,1]}<\xi_{1}. $$

Theorem 3.6

Suppose that there exist four positive numbers \(\xi_{4}<\xi_{3}<\xi_{2}<\xi_{1}\) such that one of the following conditions is satisfied:

(B5):

\(\varphi(\xi_{1})<\rho_{1}\xi_{1}\), \(\psi(\xi_{2})>\rho _{2}\xi_{2}\), \(\varphi(\xi_{3})<\rho_{1}\xi_{3}\), \(\psi(\xi_{4})>\rho_{2}\xi _{4} \).

(B6):

\(\psi(\xi_{1})>\rho_{2}\xi_{1}\), \(\varphi(\xi_{2})<\rho _{1}\xi_{2}\), \(\psi(\xi_{3})>\rho_{2}\xi_{3}\), \(\varphi(\xi_{4})<\rho _{1}\xi_{4}\).

Then BVP (1.2)-(1.3) has at least three positive solutions \(u_{1}, u_{2}, u_{3}\in K\) such that

$$ \xi_{4}<\|u_{1}\|_{[0,1]}<\xi_{3}< \|u_{2}\|_{[0,1]}<\xi_{2}<\|u_{3}\| _{[0,1]}<\xi_{1}. $$

Theorem 3.7

Suppose that there exist \(n+1\) positive numbers \(\xi_{n+1}<\xi_{n}<\cdots<\xi_{2}<\xi_{1}\) such that one of the following conditions is satisfied:

(B7):

\(\varphi(\xi_{2k-1})< \rho_{1}\xi_{2k-1}\), \(\psi (\xi_{2k})>\rho_{2}\xi_{2k}\), \(k=1,2,\ldots, [\frac {n+2}{2} ]\).

(B8):

\(\psi(\xi_{2k-1})>\rho_{2}\xi_{2k-1}\), \(\varphi (\xi_{2k})<\rho_{1}\xi_{2k}\), \(k=1,2,\ldots, [\frac {n+2}{2} ] \).

Then BVP (1.2)-(1.3) has at least n positive solutions \(u_{i}\in K\) (\(i=1,2,\ldots,n\)) such that

$$ \xi_{i+1}<\|u_{i}\|_{[0,1]}<\xi_{i}. $$

4 Example

Consider the following four-point BVP of delayed nonlinear fractional differential equations:

$$\begin{aligned}& \left \{\begin{array}{l} D_{0+}^{\frac{5}{2}}u(t)+f (t,u (t-\frac{1}{6} ),u (t+\frac{1}{9} ) ), \quad t\in[0,1], \\ u(0)=u'(0)=0, \qquad u(1)=8(I_{0+}^{\frac{1}{2}}u) (\frac {1}{4} )+5(I_{0+}^{\frac{3}{2}}u) (\frac{1}{2} ), \end{array} \right . \end{aligned}$$
(4.1)

where \(f(t,x,y)=\frac{4(t^{2}+t+1)}{3}+\frac{x^{2}+y^{2}}{8}\), \(0\leq t\leq1\), \(x, y\geq0\), \(q=5/2\), \(\tau_{1}=\frac{1}{6}\), \(\tau_{2}=\frac{1}{9}\), \(\alpha_{1}=8\), \(\alpha_{2}=5\), \(p_{1}=\frac{1}{2}\), \(p_{2}=\frac{3}{2}\), \(\zeta_{1}=\frac{1}{4}\), \(\zeta_{2}=\frac{1}{2}\), \(m=4\). Choosing \(\theta=\frac{1}{4}\), \(r_{1}=2\), \(r_{2}=\frac{1}{50}\), \(r_{0}=r'_{0}=\frac{1}{5}\).

By a simple calculation, we get

$$\begin{aligned}& \sigma_{1}=\min_{\frac{1}{12}\leq t\leq \frac{7}{12}}\sigma(t)=\sigma \biggl( \frac{1}{12} \biggr)=\frac {11\sqrt{3}}{1\text{,}286}, \\& \sigma_{2}=\min_{\frac{13}{36}\leq t\leq \frac{31}{36}}\sigma(t)=\sigma \biggl( \frac{31}{36} \biggr)=\frac {155\sqrt{31}}{11\text{,}664}, \\& \Omega=1-\sum_{i=1}^{2}\frac{\alpha_{i}\zeta_{i}^{p_{i}+q-1}\Gamma (q)}{\Gamma(p_{i}+q)} \approx0.5827572. \end{aligned}$$

By Lemmas 2.3 and 2.4 and the aid of a computer, we obtain

$$\begin{aligned} \rho_{1} &= \biggl(\int_{0}^{1}G(s)\,ds \biggr)^{-1} = \Biggl(\int_{0}^{1} \frac{s(1-s)^{q-1}}{\Gamma(q-1)} \,ds+\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Omega\Gamma(p_{i}+q)}\int_{0}^{1}g_{i}( \zeta_{i},s)\,ds \Biggr)^{-1} \\ &= \Biggl(\frac{q-1}{\Gamma(q+2)}+\sum_{i=1}^{m-2} \frac{\alpha _{i}}{\Omega\Gamma(p_{i}+q)} \biggl(\int_{0}^{1}\zeta _{i}^{p_{i}+q-1}(1-s)^{q-1}\,ds-\int_{0}^{\zeta_{i}}( \zeta _{i}-s)^{p_{i}+q-1}\,ds \biggr) \Biggr)^{-1} \\ &= \Biggl(\frac{q-1}{\Gamma(q+2)}+\sum_{i=1}^{m-2} \frac{\alpha _{i}\zeta_{i}^{p_{i}+q-1}(p_{i}+q(1-\zeta_{i}))}{\Omega\Gamma (p_{i}+q+1)q} \Biggr)^{-1} \\ &= \biggl(\frac{8}{35\sqrt{\pi}}+\frac{19}{240\Omega}+\frac{11}{384\Omega} \biggr)^{-1}\approx3.025561 \end{aligned}$$

and

$$\begin{aligned} \rho_{2} =& \biggl(\sigma(\theta)\int_{\theta}^{1-\theta}G(s) \,ds \biggr)^{-1} \\ =& \biggl(\sigma \biggl(\frac{1}{4} \biggr) \biggr)^{-1} \Biggl(\int_{\frac{1}{4}}^{\frac{3}{4}}\frac{s(1-s)^{q-1}}{\Gamma(q-1)} \,ds+\sum _{i=1}^{m-2}\frac{\alpha_{i}}{\Omega\Gamma(p_{i}+q)}\int _{\frac{1}{4}}^{\frac{3}{4}}g_{i}(\zeta_{i},s) \,ds \Biggr)^{-1} \\ =&16 \Biggl(\frac{117\sqrt{3}-23}{1\text{,}120\sqrt{\pi}}+\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Omega\Gamma(p_{i}+q)} \\ &{}\times \biggl(\int_{\frac{1}{4}}^{\frac{3}{4}} \zeta_{i}^{p_{i}+q-1}(1-s)^{q-1}\,ds-\int _{\frac{1}{4}}^{\zeta_{i}}(\zeta_{i}-s)^{p_{i}+q-1} \,ds \biggr) \Biggr)^{-1} \\ =&16 \Biggl(\frac{117\sqrt{3}-23}{1\text{,}120\sqrt{\pi}}+\sum_{i=1}^{m-2} \frac{\alpha_{i}}{\Omega\Gamma(p_{i}+q)} \biggl[\frac {\zeta_{i}^{p_{i}+q-1}(3^{q}-1)}{q4^{q}}-\frac{(\zeta_{i}-\frac {1}{4})^{p_{i}+q}}{p_{i}+q} \biggr] \Biggr)^{-1} \\ =&16 \biggl(\frac{117\sqrt{3}-23}{1\text{,}120\sqrt{\pi}}+\frac{9\sqrt{3}-1}{320\Omega}+\frac{72\sqrt{3}-13}{6\text{,}144\Omega } \biggr)^{-1}\approx80.029634. \end{aligned}$$

Therefore, for \((t,x,y)\in[\frac{1}{4}, \frac{3}{4}]\times[\frac{11\sqrt{3}}{64\text{,}800}, \frac{1}{50}]\times[\frac{155\sqrt{31}}{583\text{,}200}, \frac {1}{50}]\), we have

$$\begin{aligned} f(t,x,y)&=\frac{4(t^{2}+t+1)}{3}+\frac{x^{2}+y^{2}}{8}\geq f \biggl(\frac{1}{4}, \frac{11\sqrt{3}}{64\text{,}800}, \frac{155\sqrt {31}}{583\text{,}200} \biggr) \\ &>f \biggl(\frac{1}{4}, 0, 0 \biggr)=1.75>\rho_{2}r_{2} \approx1.600593. \end{aligned}$$

For \((t,x,y)\in[0,1]\times[0,\frac{11}{5}]\times[0,\frac{11}{5}]\), we obtain

$$ f(t,x,y)=\frac{4(t^{2}+t+1)}{3}+\frac{x^{2}+y^{2}}{8}\leq f \biggl(1, \frac{11}{5}, \frac{11}{5} \biggr)=5.21<\rho _{1}r_{1} \approx6.051122 $$

and

$$\begin{aligned} f_{\infty}&=\liminf_{x+y\rightarrow +\infty}\min_{t\in[0,1]} \frac{f(t,x,y)}{x+y} \\ &=\liminf_{x+y\rightarrow +\infty}\min_{t\in[0,1]}\frac{1}{x+y} \biggl(\frac {4(t^{2}+t+1)}{3}+\frac{x^{2}+y^{2}}{8} \biggr) \\ &=\liminf_{x+y\rightarrow +\infty}\frac{1}{x+y} \biggl(\frac{4}{3}+ \frac{x^{2}+y^{2}}{8} \biggr) \\ &\geq\liminf_{x+y\rightarrow +\infty}\frac{1}{x+y} \biggl( \frac{4}{3}+\frac{(x+y)^{2}}{16} \biggr)=\infty. \end{aligned}$$

Thus (H1)-(H3) hold. With the use of Theorem 3.1, BVP (4.1) has at least two positive solutions \(u_{1}\) and \(u_{2}\) such that \(0<\frac{1}{50}<\|u_{1}\|_{[0,1]}<2<\|u_{2}\|_{[0,1]}\).

References

  1. Babakhani, A, Gejji, V: Existence of positive solutions of nonlinear fractional differential equations. J. Math. Anal. Appl. 278, 434-442 (2003)

    Article  MATH  MathSciNet  Google Scholar 

  2. Lakshmikantham, V, Leela, S: Nagumo-type uniqueness result for fractional differential equations. Nonlinear Anal. 71, 2886-2889 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  3. Feng, M, Ge, W: Existence results for a class of nth order m-point boundary value problems in Banach spaces. Appl. Math. Lett. 22, 1303-1308 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  4. Chang, Y-K, 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. Ahmad, B, Nieto, JJ: Existence results for nonlinear boundary value problems of fractional integrodifferential equations with integral boundary conditions. Bound. Value Probl. 2009, Article ID 708576 (2009)

    MathSciNet  MATH  Google Scholar 

  6. Abbas, S: Existence of solutions to fractional order ordinary and delay differential equations and applications. Electron. J. Differ. Equ. 2011, 9 (2011)

    Article  MathSciNet  Google Scholar 

  7. Bai, Z: On positive solutions of a nonlocal fractional boundary value problem. Nonlinear Anal. 72, 916-924 (2010)

    Article  MATH  MathSciNet  Google Scholar 

  8. Bai, C, Ma, J: Eigenvalue criteria for existence of multiple positive solutions to boundary value problems of second-order delay differential equations. J. Math. Anal. Appl. 301(2), 457-476 (2005)

    Article  MATH  MathSciNet  Google Scholar 

  9. Hong, CH, Yeh, CC, Lee, CF: Existence of positive solutions for functional differential equations. Comput. Math. Appl. 40(6), 783-792 (2000)

    Article  MATH  MathSciNet  Google Scholar 

  10. Anguraja, A, Latha Maheswari, M: Existence of solutions for fractional impulsive neutral functional infinite delay integro-differential equations with nonlocal conditions. J. Nonlinear Sci. Appl. 5, 271-280 (2012)

    MathSciNet  Google Scholar 

  11. 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, 28 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  12. Zhang, XQ, Wang, L, Sun, Q: Existence of positive solutions for a class of nonlinear fractional differential equations with integral boundary conditions and a parameter. Appl. Math. Comput. 226, 708-718 (2014)

    Article  MathSciNet  Google Scholar 

  13. Benchohra, M, Henderson, J, Ntouyas, SK, Ouahab, A: Existence results for fractional order functional differential equations with infinite delay. J. Math. Anal. Appl. 338, 1340-1350 (2008)

    Article  MATH  MathSciNet  Google Scholar 

  14. Liao, C, Ye, H: Existence of positive solutions of nonlinear fractional delay differential equations. Positivity 13, 601-609 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  15. Goodrich, C: Existence of a positive solution to a class of fractional differential equations. Appl. Math. Lett. 23, 1050-1055 (2010)

    Article  MATH  MathSciNet  Google Scholar 

  16. 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 

  17. Bai, Z, Lv, H: Positive solutions for boundary value problem of nonlinear fractional differential equation. J. Math. Anal. Appl. 311, 495-505 (2005)

    Article  MATH  MathSciNet  Google Scholar 

  18. Zhao, D, Liu, Y: Multiple positive solutions for nonlinear fractional boundary value problems. Sci. World J. 2013, Article ID 473828 (2013)

    Google Scholar 

  19. Xu, X, Jiang, D, Yuan, C: Multiple positive solutions for the boundary value problem of a nonlinear fractional differential equation. Nonlinear Anal., Theory Methods Appl. 71(10), 4676-4688 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  20. Jiang, J, Liu, L, Wu, Y: Multiple positive solutions of singular fractional differential system involving Stieltjes integral conditions. Electron. J. Qual. Theory Differ. Equ. 2012, 43 (2012)

    Article  MathSciNet  MATH  Google Scholar 

  21. Zhao, KH, Gong, P: Existence of positive solutions for a class of higher-order Caputo fractional differential equation. In: Qualitative Theory of Dynamical Systems (2014). doi:10.1007/s12346-014-0121-0

    Google Scholar 

  22. 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, 254 (2014)

    Article  MathSciNet  Google Scholar 

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

    Article  MathSciNet  Google Scholar 

  24. Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of Fractional Differential Equations. North-Holland Mathematics Studies, vol. 204. Elsevier, Amsterdam (2006)

    Book  MATH  Google Scholar 

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

    MATH  Google Scholar 

  26. Guo, D, Lakshmikantham, V, Liu, X: Nonlinear Integral Equations in Abstract Spaces. Mathematics and Its Applications, vol. 373. Kluwer Academic, Dordrecht (1996)

    Book  MATH  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the anonymous referees for their useful and valuable suggestions. This work was 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 authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally and significantly in writing this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhao, K., Gong, P. Positive solutions of m-point multi-term fractional integral BVP involving time-delay for fractional differential equations. Bound Value Probl 2015, 19 (2015). https://doi.org/10.1186/s13661-014-0280-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-014-0280-6

MSC

Keywords