Skip to main content

On the periodic boundary value problem for Duffing type fractional differential equation with p-Laplacian operator

Abstract

By using the continuation theorem of coincidence degree theory, we study the existence of solutions of Duffing type fractional differential equations with a p-Laplacian operator. Under certain nonlinear growth conditions of the nonlinearity, we obtain a new result on the existence of solutions.

1 Introduction

Fractional calculus is a generalization of ordinary differentiation and integration on an arbitrary order that may be noninteger. Fractional differential equations have been proved to be valuable tools in the modeling of many phenomena in various fields of science and engineering. We can see numerous applications in viscoelasticity, neurons, electrochemistry, control, etc. (see [1–6]). Recently, with the intensive development of the theory of fractional calculus itself and its applications, there have many important results of fractional differential equations on initial value problems, and boundary value problems at nonresonance and resonance (see [7–12]).

In the study of the turbulent flow in a porous medium, Leibenson (see [13]) introduced the p-Laplacian equation as follows:

$$ \bigl(\phi_{p} \bigl(x'(t)\bigr) \bigr)'=f\bigl(t,x(t),x'(t)\bigr), $$
(1.1)

where \(\phi_{p}(s)=|s|^{p-2}s\), \(p>1\). Obviously, \(\phi_{p}\) is invertible and its inverse operator is \(\phi_{q}\), where \(q>1\) is a constant such that \(1/p+1/q=1\). In the past few decades, many important results relative to (1.1) with certain boundary value conditions have been obtained. See the papers [14–20] and the references therein. However, as far as we know, work on the existence of solutions for periodic boundary value problems (PBVPs for short) of fractional differential equations was discussed less.

The aim of this paper is to concentrate on the following periodic boundary value problem for Duffing type fractional differential equations with p-Laplacian operator:

$$ \textstyle\begin{cases} D_{0^{+}}^{\beta}(\phi_{p}(D_{0^{+}}^{\alpha}x(t)) )+g(t,x(t))=e(t), \quad t\in[0,T],\\ x(0)=x(T), \quad\quad D_{0^{+}}^{\alpha}x(0)=D_{0^{+}}^{\alpha}x(T), \end{cases} $$
(1.2)

where \(0<\alpha,\beta\leq1\), \(D_{0^{+}}^{\alpha},D_{0^{+}}^{\beta}\) are Caputo fractional derivatives, \(T>0\) is a given constant, and \(g:[0,T]\times \mathbb{R}\rightarrow\mathbb{R}\), \(e:[0,T]\rightarrow\mathbb{R}\) are continuous. Throughout this paper, we assume that

$$ \int_{0}^{T} (T-s)^{\beta-1}e(s)\,ds=0. $$

The choice of periodic boundary conditions is motivated by the difficulty in the study of the PBVP

$$ \textstyle\begin{cases} D_{0^{+}}^{\beta}(\phi_{p}(D_{0^{+}}^{\alpha}x(t)) )=h(t),\quad t\in [0,T],\\ x(0)=x(T),\quad\quad D_{0^{+}}^{\alpha}x(0)=D_{0^{+}}^{\alpha}x(T). \end{cases} $$
(1.3)

As we know, PBVP (1.3) is not solvable for each \(h\in C([0,T],\mathbb{R})\), and, when solvable, has no unique solution because \(x(t)+c\), \(\forall c\in\mathbb{R}\) is a solution together with \(x(t)\). A trivial necessary condition for the solvability of PBVP (1.3) is that \(\int_{0}^{T}(T-s)^{\beta-1}h(s)\,ds=0\).

Notice that \(D_{0^{+}}^{\beta}(\phi_{p}(D_{0^{+}}^{\alpha}) )\) is a nonlinear operator, so the coincidence degree theory for linear differential operators is invalid in the direct application to it.

The rest of this paper is organized as follows. In Section 2, we describe the fractional differential operator and some lemmas. In Section 3, some sufficient conditions for the existence of solutions for PBVP (1.2) are established, and a new result on the existence of solutions is obtained. Finally, in Section 4, an example is given to illustrate the main result.

2 Preliminaries

Some definitions of the fractional derivative have emerged over the years (see [21, 22]), and in this paper we restrict our attention to the use of the Caputo fractional derivative. In this section, we introduce some basic definitions and lemmas which will be used in what follows. For details, we refer the reader to [21–25].

Definition 2.1

The Riemann-Liouville fractional integral operator 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 side integral is pointwise defined on \((0,+\infty)\), where \(\Gamma(\cdot) >0\) is the Gamma function.

Definition 2.2

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

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

where n is the smallest integer greater than or equal to α, provided that the right side integral is pointwise defined on \((0,+\infty)\).

Lemma 2.1

(see [23])

Let \(\alpha>0\). Assume that u, \(D_{0^{+}}^{\alpha}u\in L([0,T],\mathbb{R})\). Then the following equality holds:

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

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

Lemma 2.2

(see [24])

For any \(u,v\geq0\),

$$ \phi_{p}(u+v)\leq \textstyle\begin{cases} \phi_{p}(u)+\phi_{p}(v),& \textit{if } p< 2;\\ 2^{p-2} (\phi_{p}(u)+\phi_{p}(v) ),& \textit{if } p\geq2. \end{cases} $$

Now we briefly recall some notations and an abstract existence result, which can be found in [25].

Let X, Y be real Banach spaces, \(L: \operatorname{dom}L\subset X\rightarrow Y\) be a Fredholm operator with index zero, and \(P: X\rightarrow X\), \(Q:Y\rightarrow Y \) be projectors such that

$$\begin{aligned}& \operatorname{Im}P=\operatorname{Ker}L,\quad\quad \operatorname{Ker}Q=\operatorname{Im}L, \\& X=\operatorname{Ker}L\oplus\operatorname{Ker}P,\quad\quad Y=\operatorname{Im}L\oplus\operatorname{Im}Q. \end{aligned}$$

It follows that

$$ L|_{\operatorname{dom}L\cap\operatorname{Ker}P}: \operatorname{dom}L\cap\operatorname{Ker}P\rightarrow \operatorname{Im}L $$

is invertible. We denote the inverse by \(K_{P}\).

If Ω is an open bounded subset of X such that \(\operatorname{dom}L\cap \bar{\Omega}\neq\varnothing\), then the map \(N:X\rightarrow Y\) will be called L-compact on \(\bar{\Omega}\) if \(QN(\bar{\Omega})\) is bounded and \(K_{P}(I-Q)N:\bar{\Omega}\rightarrow X\) is compact.

Lemma 2.3

(see [25])

Let X and Y be two Banach spaces, \(L:\operatorname{dom}L\subset X\rightarrow Y\) be a Fredholm operator with index zero, \(\Omega\subset X\) be an open bounded set, and \(N:\bar {\Omega}\rightarrow Y\) be L-compact on \(\bar{\Omega}\). Suppose that all of the following conditions hold:

  1. (1)

    \(Lx\neq\lambda Nx\), \(\forall x\in\partial\Omega\cap\operatorname{dom}L\), \(\lambda\in(0,1)\);

  2. (2)

    \(QNx\neq0\), \(\forall x\in\partial\Omega\cap\operatorname{Ker}L\);

  3. (3)

    \(\operatorname{deg}(JQN,\Omega\cap\operatorname{Ker}L,0)\neq0\), where \(J: \operatorname{Im}Q\rightarrow\operatorname{Ker}L\) is an isomorphism map.

Then the equation \(Lx=Nx\) has at least one solution on \(\bar{\Omega}\cap\operatorname{dom}L\).

3 Existence result

For making use of the continuation theorem to study the existence of solutions for PBVP (1.2), we consider a system as follows:

$$ \textstyle\begin{cases} D_{0^{+}}^{\alpha}x_{1}(t)=\phi_{q} (x_{2}(t)),\\ D_{0^{+}}^{\beta}x_{2}(t)=e(t)-g(t,x_{1}(t)),\\ x_{1}(0)=x_{1}(T), \quad\quad x_{2}(0)=x_{2}(T). \end{cases} $$
(3.1)

Clearly, if \(x(\cdot)=(x_{1}(\cdot),x_{2}(\cdot))^{\mathrm{T}}\) is a solution of PBVP (3.1), then \(x_{1}(\cdot)\) must be a solution of PBVP (1.2). So, to prove PBVP (1.2) has solutions, we only need to show that PBVP (3.1) has solutions.

In this paper, we take \(X=\{x=(x_{1},x_{2})^{\mathrm{T}}\mid x_{1},x_{2}\in C([0,T],\mathbb{R})\}\) with the norm \(\|x\|=\max\{\|x_{1}\|_{0},\|x_{2}\|_{0}\} \), where \(\|x_{i}\|_{0}=\max_{t\in[0,T]}|x_{i}(t)|\) (\(i\in\{1,2\}\)). By means of the linear functional analysis theory, we can prove X is a Banach space.

Define the operator \(L:\operatorname{dom}L\subset X\rightarrow X\) by

$$ Lx= \begin{pmatrix} D_{0^{+}}^{\alpha}x_{1}\\ D_{0^{+}}^{\beta}x_{2} \end{pmatrix} , $$
(3.2)

where

$$\begin{aligned} \operatorname{dom}L =&\bigl\{ x\in X\mid D_{0^{+}}^{\alpha}x_{1},D_{0^{+}}^{\beta}x_{2}\in C\bigl([0,T],\mathbb{R}\bigr), \\ &x_{1}(0)=x_{1}(T),x_{2}(0)=x_{2}(T) \bigr\} . \end{aligned}$$

Let \(N:X\rightarrow X\) be defined by

$$ Nx(t)= \begin{pmatrix} \phi_{q} (x_{2}(t))\\ e(t)-g(t,x_{1}(t)) \end{pmatrix} , \quad\forall t\in[0,T]. $$
(3.3)

It is easy to see that PBVP (3.1) can be converted to the operator equation

$$ Lx=Nx,\quad x\in\operatorname{dom}L. $$

Now let us introduce some lemmas.

Lemma 3.1

Let L be defined by (3.2), then

$$\begin{aligned}& \operatorname{Ker}L=\bigl\{ x\in X\mid x(t)=c, \forall t\in[0,T],c\in \mathbb{R}^{2}\bigr\} , \end{aligned}$$
(3.4)
$$\begin{aligned}& \operatorname{Im}L= \Biggl\{ y\in X \Biggm| \begin{pmatrix} \int_{0}^{T} (T-s)^{\alpha-1} y_{1}(s)\,ds\\ \int_{0}^{T} (T-s)^{\beta-1} y_{2}(s)\,ds \end{pmatrix}=0 \Biggr\} . \end{aligned}$$
(3.5)

Proof

Obviously, from Lemma 2.1, we can see that (3.4) holds.

If \(y\in\operatorname{Im}L\), then there exists \(x\in\operatorname{dom}L\) such that \(y=Lx\). That is, \(y_{1}(t)=D_{0^{+}}^{\alpha}x_{1}(t)\), \(y_{2}(t)=D_{0^{+}}^{\beta}x_{2}(t)\). By using Lemma 2.1, we have

$$\begin{aligned}& x_{1}(t)=c_{1}+\frac{1}{\Gamma(\alpha)}\int _{0}^{t}(t-s)^{\alpha-1}y_{1}(s)\,ds,\quad c_{1}\in\mathbb{R}, \\& x_{2}(t)=c_{2}+\frac{1}{\Gamma(\beta)}\int_{0}^{t}(t-s)^{\beta-1}y_{2}(s)\,ds, \quad c_{2}\in\mathbb{R}. \end{aligned}$$

From conditions \(x_{1}(0)=x_{1}(T)\), \(x_{2}(0)=x_{2}(T)\), we obtain

$$ \begin{pmatrix} \int_{0}^{T} (T-s)^{\alpha-1} y_{1}(s)\,ds\\ \int_{0}^{T} (T-s)^{\beta-1} y_{2}(s)\,ds \end{pmatrix}=0. $$
(3.6)

On the other hand, suppose \(y\in X\) and satisfies (3.6). Let \(x_{1}(t)=I_{0^{+}}^{\alpha}y_{1}(t)\), \(x_{2}(t)=I_{0^{+}}^{\beta}y_{2}(t)\). Obviously \(x_{1}(0)=x_{1}(T)\), \(x_{2}(0)=x_{2}(T)\). Hence \(x=(x_{1},x_{2})^{\mathrm{T}}\in\operatorname {dom}L\) and \(Lx=y\). So \(y\in\operatorname{Im}L\). The proof is complete. □

Lemma 3.2

Let L be defined by (3.2), then L is a Fredholm operator of index zero. The projectors \(P:X \rightarrow X\) and \(Q:X \rightarrow X\) can be defined as

$$\begin{aligned}& Px(t)=x(0), \quad\forall t\in[0,T], \\& Qy(t)=\begin{pmatrix} \frac{\alpha}{T^{\alpha}}\int_{0}^{T} (T-s)^{\alpha-1}y_{1}(s)\,ds\\ \frac{\beta}{T^{\beta}}\int_{0}^{T} (T-s)^{\beta-1}y_{2}(s)\,ds \end{pmatrix}:=\begin{pmatrix} (Qy)_{1}\\ (Qy)_{2} \end{pmatrix}, \quad\forall t\in[0,T]. \end{aligned}$$

Furthermore, the operator \(K_{P}: \operatorname{Im}L\rightarrow{\operatorname{dom}L\cap \operatorname{Ker}P}\) can be written as

$$ K_{P}y=\begin{pmatrix} I_{0^{+}}^{\alpha}y_{1}\\ I_{0^{+}}^{\beta}y_{2} \end{pmatrix}, $$

which is \((L|_{\operatorname{dom}L\cap\operatorname{Ker}P})^{-1}\).

Proof

For any \(y \in X\), we have

$$ Q^{2}y=Q\begin{pmatrix} (Qy)_{1}\\ (Qy)_{2} \end{pmatrix} =\begin{pmatrix} (Qy)_{1}\frac{\alpha}{T^{\alpha}}\int_{0}^{T}(T-s)^{\alpha-1}\,ds\\ (Qy)_{2}\frac{\beta}{T^{\beta}}\int_{0}^{T}(T-s)^{\beta-1}\,ds \end{pmatrix} =Qy. $$

Let \(y^{*}=y-Qy=\bigl( {\scriptsize\begin{matrix}{}y_{1}^{*}\cr y_{2}^{*} \end{matrix}} \bigr)\), then we get

$$\begin{aligned} \int_{0}^{T}(T-s)^{\alpha-1}y_{1}^{*}(s)\,ds =& \int_{0}^{T}(T-s)^{\alpha-1}y_{1}(s)\,ds- \int_{0}^{T}(T-s)^{\alpha-1}(Qy)_{1}\,ds \\ = &\frac{T^{\alpha}}{\alpha}\bigl((Qy)_{1}-\bigl(Q^{2}y \bigr)_{1}\bigr)=0. \end{aligned}$$

Similarly, we have \(\int_{0}^{T}(T-s)^{\beta-1}y_{2}^{*}(s)\,ds=0\). So \(y^{*}\in \operatorname{Im}L\). Hence \(X=\operatorname{Im}L+\operatorname{Im}Q\). Since \(\operatorname{Im}L\cap \operatorname{Im}Q=\{0\}\), we have \(X=\operatorname{Im}L\oplus\operatorname{Im}Q\). Thus

$$\begin{aligned}& \operatorname{dim} \operatorname{Ker}L=\operatorname{dim}\operatorname{ Im}Q= \operatorname{codim}\operatorname{Im}L=2. \end{aligned}$$

This means that L is a Fredholm operator of index zero.

From the definition of \(K_{P}\), for \(y\in\operatorname{Im}L\), we have

$$ LK_{P}y=\begin{pmatrix} D_{0^{+}}^{\alpha}I_{0^{+}}^{\alpha}y_{1}\\ D_{0^{+}}^{\beta}I_{0^{+}}^{\beta}y_{2} \end{pmatrix}=y. $$

On the other hand, for \(x\in\operatorname{dom}L\cap\operatorname{Ker}P\), we have \(x_{1}(0)=x_{2}(0)=0\). By Lemma 2.1, we get

$$ K_{P}Lx=\begin{pmatrix} x_{1}-x_{1}(0)\\ x_{2}-x_{2}(0) \end{pmatrix}=x. $$

So, we know that \(K_{P}=(L_{\operatorname{dom}L\cap\operatorname{Ker}P})^{-1}\). The proof is complete. □

Lemma 3.3

Let N be defined by (3.3). Assume \(\Omega\subset X\) is an open bounded subset such that \(\operatorname{dom}L\cap \bar{\Omega}\neq\varnothing\), then N is L-compact on \(\bar{\Omega}\).

Proof

By the continuity of \(\phi_{q}\), e, g, we find that \(QN(\bar{\Omega})\) and \(K_{P} (I-Q)N(\bar{\Omega})\) are bounded. Moreover, there exists a constant \(M>0\) such that \(\|(I-Q)Nx\|\leq M\), \(\forall x\in\bar{\Omega}\), \(t\in[0,T]\). Thus, in view of the Arzelà-Ascoli theorem, we only need prove that \(K_{P}(I-Q)N(\bar{\Omega})\subset X\) is equicontinuous.

For \(0\leq t_{1}< t_{2}\leq T\), \(x\in\bar{\Omega}\), we have

$$\begin{aligned}& K_{P}(I-Q)Nx(t_{2})-K_{P}(I-Q)Nx(t_{1}) \\& \quad= \begin{pmatrix} I_{0^{+}}^{\alpha}((I-Q)Nx)_{1}(t_{2})-I_{0^{+}}^{\alpha}((I-Q)Nx)_{1}(t_{1})\\ I_{0^{+}}^{\beta}((I-Q)Nx)_{2}(t_{2})-I_{0^{+}}^{\beta}((I-Q)Nx)_{2}(t_{1}) \end{pmatrix}. \end{aligned}$$

From \(\|(I-Q)Nx\|\leq M\), \(\forall x\in\bar{\Omega}\), \(t\in[0,T]\), we can see that

$$\begin{aligned}& \biggl\vert \frac{1}{\Gamma(\alpha)} \biggl(\int_{0}^{t_{2}}(t_{2}-s)^{\alpha -1} \bigl((I-Q)Nx\bigr)_{1}(s)\,ds -\int_{0}^{t_{1}}(t_{1}-s)^{\alpha-1} \bigl((I-Q)Nx\bigr)_{1}(s)\,ds \biggr)\biggr\vert \\& \quad\leq \frac{M}{\Gamma(\alpha)} \biggl(\int_{0}^{t_{1}} \bigl((t_{1}-s)^{\alpha -1}-(t_{2}-s)^{\alpha-1} \bigr)\,ds+\int_{t_{1}}^{t_{2}}(t_{2}-s)^{\alpha-1}\,ds \biggr) \\& \quad= \frac{M}{\Gamma(\alpha+1)}\bigl(t_{1}^{\alpha}-t_{2}^{\alpha}+2(t_{2}-t_{1})^{\alpha}\bigr) \\& \quad\leq \frac{M}{\Gamma(\alpha+1)}\bigl(t_{2}^{\alpha}-t_{1}^{\alpha}+2(t_{2}-t_{1})^{\alpha}\bigr). \end{aligned}$$

Since \(t^{\alpha}\) is uniformly continuous on \([0,T]\), we find that \((K_{P}(I-Q)N(\bar{\Omega}))_{1}\subset C([0,T],\mathbb{R})\) is equicontinuous. A similar proof can show that \((K_{P}(I-Q)N(\bar{\Omega}))_{2}\subset C([0,T],\mathbb{R})\) is equicontinuous. So we find that \(K_{P}(I-Q)N:\bar{\Omega}\rightarrow X\) is compact. The proof is complete. □

Now we give the main result as regards the existence of solutions for PBVP (1.2).

Theorem 3.1

Assume that:

(H1):

there exists a constant \(d_{1}>0\) such that

$$ (-1)^{i}xg(t,x)>0 \quad\bigl(i\in\{1,2\}\bigr), \forall t\in[0,T], |x|>d_{1}; $$
(H2):

there exist a constant \(d_{2}>0\) and nonnegative functions \(a,b\in C([0,T],\mathbb{R})\) such that

$$ \bigl\vert g(t,x)\bigr\vert \leq a(t)|x|^{p-1}+b(t),\quad \forall t \in[0,T],|x|>d_{2}. $$

Then PBVP (1.2) has at least one solution, provided that

$$ \begin{aligned} &\gamma_{1}:=\frac{2^{p} T^{\beta+{\alpha p-\alpha}} \|a\|_{0}}{\Gamma (\beta+1)(\Gamma(\alpha+1))^{p-1}}< 1, \quad\textit{if } p< 2; \\ & \gamma_{2}:=\frac{2^{2p-2} T^{\beta+{\alpha p-\alpha}} \|a\|_{0}}{\Gamma (\beta+1)(\Gamma(\alpha+1))^{p-1}}< 1, \quad\textit{if } p\geq2. \end{aligned} $$
(3.7)

Proof

Set

$$ \Omega'=\bigl\{ x\in\operatorname{dom}L\mid Lx=\lambda Nx, \lambda\in(0,1)\bigr\} . $$

For \(x\in\Omega'\), we get \(Nx\in\operatorname{Im}L\). So by (3.5), we have

$$\begin{aligned}& \int_{0}^{T}(T-s)^{\alpha-1} \phi_{q}\bigl(x_{2}(s)\bigr)\,ds=0, \\& \int_{0}^{T}(T-s)^{\beta-1}\bigl(e(s)-g \bigl(s,x_{1}(s)\bigr)\bigr)\,ds=0. \end{aligned}$$

From the integral mean value theorem and \(\int_{0}^{T}(T-s)^{\beta -1}e(s)\,ds=0\), there exist constants \(\zeta,\eta\in(0,T)\) such that \(x_{2}(\zeta)=0\), \(g(\eta,x_{1}(\eta))=0\). Together with the condition (H1), we have \(|x_{1}(\eta)|\leq d_{1}\). By Lemma 2.1, we have

$$ x_{1}(t)=x_{1}(\eta)-I_{0^{+}}^{\alpha}D_{0^{+}}^{\alpha}x_{1}(\eta )+I_{0^{+}}^{\alpha}D_{0^{+}}^{\alpha}x_{1}(t), $$

which, together with

$$\begin{aligned} \bigl\vert I_{0^{+}}^{\alpha}D_{0^{+}}^{\alpha}x_{1}(t)\bigr\vert =&\frac{1}{\Gamma(\alpha)}\biggl\vert \int _{0}^{t}(t-s)^{\alpha-1}D_{0^{+}}^{\alpha}x_{1}(s)\,ds\biggr\vert \\ \leq&\frac{1}{\Gamma(\alpha)}\bigl\Vert D_{0^{+}}^{\alpha}x_{1}\bigr\Vert _{0}\cdot\frac {1}{\alpha}t^{\alpha}\\ \leq&\frac{T^{\alpha}}{\Gamma(\alpha+1)}\bigl\Vert D_{0^{+}}^{\alpha}x_{1}\bigr\Vert _{0}, \quad\forall t\in[0,T], \end{aligned}$$

and \(|x_{1}(\eta)|\leq d_{1}\), yields

$$ \|x_{1}\|_{0}\leq d_{1}+ \frac{2T^{\alpha}}{\Gamma(\alpha+1)}\bigl\Vert D_{0^{+}}^{\alpha}x_{1} \bigr\Vert _{0}. $$
(3.8)

On the other hand, if \(x\in\Omega'\), we have

$$ \textstyle\begin{cases} D_{0^{+}}^{\alpha}x_{1}(t)=\lambda\phi_{q} (x_{2}(t)),\\ D_{0^{+}}^{\beta}x_{2}(t)=\lambda(e(t)-g(t,x_{1}(t))). \end{cases} $$
(3.9)

From the first equation of (3.9), we get \(x_{2}(t)=\phi_{p}(\lambda ^{-1} D_{0^{+}}^{\alpha}x_{1}(t))\). By substituting it into the second equation of (3.9), we get

$$ D_{0^{+}}^{\beta}\bigl(\phi_{p}\bigl(D_{0^{+}}^{\alpha}x_{1}(t)\bigr) \bigr)=\lambda^{p}e(t)- \lambda^{p}g(t,x_{1}):= \lambda^{p}N_{g}x_{1}(t). $$

Thus, by Lemma 2.1, we obtain

$$ \phi_{p}\bigl(D_{0^{+}}^{\alpha}x_{1}(t)\bigr)=c_{0}+\lambda^{p}I_{0^{+}}^{\beta}N_{g}x_{1}(t). $$
(3.10)

Then we have

$$ x_{1}(t)=c_{1}+I_{0^{+}}^{\alpha}\phi_{q}\bigl(c_{0}+\lambda^{p}I_{0^{+}}^{\beta}N_{g}x_{1}\bigr) (t). $$

By the boundary condition \(x_{1}(0)=x_{1}(T)\), we get

$$ \frac{1}{\Gamma(\alpha)}\int_{0}^{T}(T-s)^{\alpha-1} \phi_{q}\bigl(c_{0}+\lambda ^{p}I_{0^{+}}^{\beta}N_{g}x_{1}(s)\bigr)\,ds=0. $$

Obviously, there exists a constant \(\xi\in(0,T)\) such that \(\phi _{q}(c_{0}+\lambda^{p}I_{0^{+}}^{\beta}N_{g}x_{1}(\xi))=0\), which implies that \(c_{0}=-\lambda^{p}I_{0^{+}}^{\beta}N_{g}x_{1}(\xi)\). By substituting it into (3.10), we have

$$ \phi_{p}\bigl(D_{0^{+}}^{\alpha}x_{1}(t)\bigr)=-\lambda^{p}I_{0^{+}}^{\beta}N_{g}x_{1}(\xi )+\lambda^{p}I_{0^{+}}^{\beta}N_{g}x_{1}(t). $$
(3.11)

From the hypothesis (H2), we get

$$\begin{aligned} \bigl\vert I_{0^{+}}^{\beta}N_{g}x_{1}(t) \bigr\vert =&\frac{1}{\Gamma(\beta)}\biggl\vert \int_{0}^{t}(t-s)^{\beta-1}e(s)\,ds -\int_{0}^{t}(t-s)^{\beta-1}g \bigl(s,x_{1}(s)\bigr)\,ds\biggr\vert \\ \leq&\frac{T^{\beta}}{\Gamma{(\beta+1)}}\|e\|_{0} +\frac{1}{\Gamma(\beta)}\int _{0}^{t}(t-s)^{\beta-1}\bigl\vert g \bigl(s,x_{1}(s)\bigr)\bigr\vert \,ds \\ \leq&\frac{T^{\beta}}{\Gamma{(\beta+1)}}\bigl(\Vert e\Vert _{0}+\|a \|_{0}\|x_{1}\| _{0}^{p-1}+\|b \|_{0}+G_{d_{2}}\bigr), \quad\forall t\in[0,T], \end{aligned}$$

where \(G_{d_{2}}=\max \{|g(t,x)|\mid t\in[0,T],|x|\leq d_{2}\} \). Together with (3.8), (3.11), and

$$\bigl|\phi_{p}\bigl(D_{0^{+}}^{\alpha}x_{1}(t)\bigr)\bigr|=\bigl|D_{0^{+}}^{\alpha}x_{1}(t)\bigr|^{p-1}, $$

we have

$$\begin{aligned} \bigl\Vert D_{0^{+}}^{\alpha}x_{1}\bigr\Vert _{0}^{p-1} \leq&\frac{2T^{\beta}}{\Gamma(\beta+1)}\bigl(\Vert e\Vert _{0}+\|a\|_{0}\|x_{1}\|_{0}^{p-1} +\|b\|_{0}+G_{d_{2}}\bigr) \\ \leq& \frac{2T^{\beta}}{\Gamma(\beta+1)} \biggl[\Vert e\Vert _{0}+\|b\| _{0}+G_{d_{2}} +\|a\|_{0} \biggl(d_{1}+\frac{2T^{\alpha}}{\Gamma(\alpha+1)}\bigl\Vert D_{0^{+}}^{\alpha}x_{1}\bigr\Vert _{0} \biggr)^{p-1} \biggr]. \end{aligned}$$

If \(p<2\), by using Lemma 2.2, we get

$$\begin{aligned} \bigl\Vert D_{0^{+}}^{\alpha}x_{1}\bigr\Vert _{0}^{p-1} \leq&\frac{2T^{\beta}}{\Gamma(\beta+1)}\bigl(\Vert e\Vert _{0}+\|b\|_{0}+G_{d_{2}}\bigr) \\ &{}+\frac{2T^{\beta} \Vert a\Vert _{0}}{\Gamma(\beta+1)} \biggl(d_{1}^{p-1} + \frac{(2T^{\alpha})^{p-1}}{(\Gamma(\alpha+1))^{p-1}}\bigl\Vert D_{0^{+}}^{\alpha}x_{1} \bigr\Vert _{0}^{p-1} \biggr) \\ =&A_{1}+\frac{2^{p} T^{\beta+{\alpha p-\alpha}} \Vert a\Vert _{0}}{\Gamma(\beta +1)(\Gamma(\alpha+1))^{p-1}}\bigl\Vert D_{0^{+}}^{\alpha}x_{1}\bigr\Vert _{0}^{p-1}, \end{aligned}$$

where \(A_{1}=\frac{2T^{\beta}}{\Gamma(\beta+1)}(\|e\|_{0}+\|b\| _{0}+G_{d_{2}})+\frac{2T^{\beta}\|a\|_{0}}{\Gamma(\beta+1)}d_{1}^{p-1}\). Then, from (3.7), we have

$$ \bigl\Vert D_{0^{+}}^{\alpha}x_{1}\bigr\Vert _{0}\leq \biggl(\frac{A_{1}}{1-\gamma_{1}} \biggr)^{q-1}:=M_{1}. $$

Thus, from (3.8), we get

$$ \Vert x_{1}\Vert _{0}\leq d_{1}+\frac{2T^{\alpha}}{\Gamma(\alpha+1)}M_{1}. $$
(3.12)

If \(p\geq2\), similar to the above argument, let \(A_{2}=\frac{2T^{\beta}}{\Gamma(\beta+1)}(\|e\|_{0}+\|b\|_{0}+G_{d_{2}})+\frac{2^{p-1}T^{\beta}\|a\| _{0}}{\Gamma(\beta+1)}d_{1}^{p-1}\), we obtain

$$ \|x_{1}\|_{0}\leq d_{1}+ \frac{2T^{\alpha}}{\Gamma(\alpha+1)}M_{2}, $$
(3.13)

where \(M_{2}=(\frac{A_{2}}{1-\gamma_{2}})^{q-1}\). So, combining (3.12) with (3.13), we get

$$ \|x_{1}\|_{0}\leq\max \biggl\{ d_{1}+\frac{2T^{\alpha}}{\Gamma(\alpha +1)}M_{1},d_{1}+ \frac{2T^{\alpha}}{\Gamma(\alpha+1)}M_{2} \biggr\} :=M. $$
(3.14)

From the second equation of (3.9) and Lemma 2.1, we have

$$ x_{2}(t)=c+\lambda I_{0^{+}}^{\beta}N_{g}x_{1}(t), $$

which together with \(x_{2}(\zeta)=0\) yields

$$ x_{2}(t)=-\lambda I_{0^{+}}^{\beta}N_{g}x_{1}( \zeta)+\lambda I_{0^{+}}^{\beta}N_{g}x_{1}(t). $$

Then we have

$$ \|x_{2}\|_{0}\leq\frac{2T^{\beta}}{\Gamma(\beta+1)}\bigl(\Vert e\Vert _{0}+G_{M}\bigr):=\overline{M}, $$

where \(G_{M}=\max\{|g(t,x)| \mid t\in[0,T],|x|\leq M\}\). Together with (3.14), we obtain

$$ \|x\|=\max\bigl\{ \Vert x_{1}\Vert _{0}, \|x_{2}\|_{0}\bigr\} \leq\max\{M,\overline{M} \}:=M_{0}. $$

Let \(\Omega=\{x\in X\mid\|x\|< M_{0}+1\}\). From the above argument, we know that the equation

$$ Lx=\lambda Nx,\quad \forall\lambda\in(0,1) $$

has no solution on \(\partial\Omega\cap\operatorname{dom}L \). So the condition (1) of Lemma 2.3 is satisfied. Next the other two conditions of Lemma 2.3 are to be verified.

For \(x\in\operatorname{Ker}L\), we have \(x_{1}(t)=c_{1}\), \(x_{2}(t)=c_{2}\), \(\forall t\in [0,T]\), \(c_{1},c_{2}\in\mathbb{R}\). If \(QNx=0\), we obtain

$$\begin{aligned}& \int_{0}^{T}(T-s)^{\alpha-1} \phi_{q}(c_{2})\,ds=0, \\& \int_{0}^{T}(T-s)^{\beta-1}g(s,c_{1})\,ds=0. \end{aligned}$$

From the first equality, we get \(c_{2}=0\). From the second equality and (H1), we have \(|c_{1}|\leq d_{1}\). So \(\|x\|=\max\{|c_{1}|,|c_{2}|\} \leq d_{1}< M_{0}+1\). Then the condition (2) of Lemma 2.3 is satisfied.

Define the operators \(J:\operatorname{Im}Q\rightarrow\operatorname{Ker}L\) by

$$ J(x_{1},x_{2})^{\mathrm{T}}=\bigl((-1)^{i+1}x_{2},x_{1} \bigr)^{\mathrm{T}}, $$

and \(F: [0,1]\times\bar{\Omega}\rightarrow X\) by

$$ F(\mu,x) = \mu x+(1-\mu)JQNx = \begin{pmatrix} \mu x_{1}+(-1)^{i}(1-\mu)\frac{\beta}{T^{\beta}}\int_{0}^{T}(T-s)^{\beta -1}g(s,x_{1}(s))\,ds\\ \mu x_{2}+(1-\mu)\frac{\alpha}{T^{\alpha}}\int_{0}^{T}(T-s)^{\alpha-1} \phi_{q}(x_{2}(s))\,ds \end{pmatrix}, $$

where \(i\in\{1,2\}\). Let \(x\in\operatorname{Ker}L\) satisfying \(F(\mu,x)=0\), we get \(x_{1}(t)=c_{1}\), \(x_{2}(t)=c_{2}\), \(\forall t\in[0,T]\), \(c_{1},c_{2}\in\mathbb{R}\), and

$$\begin{aligned}& \mu c_{1}+(-1)^{i}(1-\mu) \frac{\beta}{T^{\beta}}\int_{0}^{T}(T-s)^{\beta -1}g(s,c_{1})\,ds=0, \end{aligned}$$
(3.15)
$$\begin{aligned}& \mu c_{2}+(1-\mu)\phi_{q}(c_{2})=0. \end{aligned}$$
(3.16)

From (3.16), we get \(c_{2}=0\) because \(c_{2}\) and \(\phi_{q}(c_{2})\) have the same sign. From (3.15), if \(\mu=0\), we get \(|c_{1}|\leq d_{1}\) because of (H1). If \(\mu\in(0,1]\), we also get \(|c_{1}|\leq d_{1}\). In fact, if \(|c_{1}|>d_{1}\), in view of (H1), one has

$$ \mu c_{1}^{2}+(1-\mu)\frac{\beta}{T^{\beta}}\int _{0}^{T}(T-s)^{\beta -1}(-1)^{i}c_{1}g(s,c_{1})\,ds>0, $$

which contradicts (3.15). From the argument above, we obtain \(\|x\| < M_{0}+1\). Thus

$$ F(\mu,x)\neq0, \quad\forall(\mu,x)\in[0,1]\times(\partial\Omega\cap \operatorname{Ker}L). $$

Hence, by the homotopy property of the degree, we have

$$\begin{aligned} \operatorname{deg}(JQN,\Omega\cap\operatorname{Ker}L,0) =&\operatorname{deg}\bigl(F(0,\cdot),\Omega\cap \operatorname{Ker}L,0\bigr) \\ =&\operatorname{deg}\bigl(F(1,\cdot),\Omega\cap\operatorname{Ker}L,0\bigr) = \operatorname{deg}(I,\Omega\cap\operatorname{Ker}L,0)\neq0. \end{aligned}$$

So the condition (3) of Lemma 2.3 is satisfied.

Consequently, by using Lemma 2.3, the operator equation \(Lx=Nx\) has at least one solution \(x(\cdot)=(x_{1}(\cdot),x_{2}(\cdot))^{\mathrm{T}} \) on \(\bar{\Omega}\cap\operatorname{dom}L\). Namely, PBVP (1.2) has at least one solution \(x_{1}(\cdot)\). The proof is complete. □

4 An example

In this section, we will give an example to illustrate our main result.

Example 4.1

Consider the following PBVP for a fractional p-Laplacian equation:

$$ \textstyle\begin{cases} D_{0^{+}}^{\frac{3}{4}} (\phi_{4} (D_{0^{+}}^{\frac{1}{2}}x(t) ) )-\frac{1}{120}x^{3}(t)+\frac{1}{2}=(1-t)^{\frac{1}{4}}\sin2\pi t,\quad t\in[0,1],\\ x(0)=x(1), \quad\quad D_{0^{+}}^{\frac{1}{2}} x(0)=D_{0^{+}}^{\frac{1}{2}} x(1). \end{cases} $$
(4.1)

Corresponding to PBVP (1.2), we get \(p=4\), \(\alpha=1/2\), \(\beta =3/4\), \(T=1\), \(e(t)=(1-t)^{\frac{1}{4}}\sin2\pi t\), and

$$ g(t,x)=-\frac{1}{120}x^{3}+\frac{1}{2}. $$

Choose \(a(t)=\frac{1}{120}\), \(b(t)=1\). By a simple calculation, we obtain

$$\begin{aligned}& xg(t,x)=-\frac{x}{120}\bigl(x^{3}-60\bigr)< 0,\quad \forall t \in[0,1], |x|>4, \\& \gamma_{2}=\frac{2^{3}/120}{\Gamma(\frac{3}{4}+1)} \biggl[\frac{2}{\Gamma(\frac{1}{2}+1)} \biggr]^{3}< 1. \end{aligned}$$

Obviously, PBVP (4.1) satisfies all assumptions of Theorem 3.1. Hence, PBVP (4.1) has at least one solution.

References

  1. Diethelm, K, Freed, AD: On the solution of nonlinear fractional order differential equations used in the modeling of viscoelasticity. In: Keil, F, Mackens, W, Voss, H, Werther, J (eds.) Scientific Computing in Chemical Engineering II-Computational Fluid Dynamics, Reaction Engineering and Molecular Properties, pp. 217-224. Springer, Heidelberg (1999)

    Google Scholar 

  2. Lundstrom, BN, Higgs, MH, Spain, WJ, Fairhall, AL: Fractional differentiation by neocortical pyramidal neurons. Nat. Neurosci. 11, 1335-1342 (2008)

    Article  Google Scholar 

  3. Glockle, WG, Nonnenmacher, TF: A fractional calculus approach of self-similar protein dynamics. Biophys. J. 68, 46-53 (1995)

    Article  Google Scholar 

  4. Hilfer, R: Applications of Fractional Calculus in Physics. World Scientific, Singapore (2000)

    Book  MATH  Google Scholar 

  5. Mainardi, F: Fractional calculus: some basic problems in continuum and statistical mechanics. In: Carpinteri, A, Mainardi, F (eds.) Fractals and Fractional Calculus in Continuum Mechanics, pp. 291-348. Springer, Wien (1997)

    Chapter  Google Scholar 

  6. Kirchner, JW, Feng, X, Neal, C: Fractal stream chemistry and its implications for contaminant transport in catchments. Nature 403, 524-526 (2000)

    Article  Google Scholar 

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

    Article  MathSciNet  Google Scholar 

  8. Delbosco, D, Rodino, L: Existence and uniqueness for a nonlinear fractional differential equation. J. Math. Anal. Appl. 204, 609-625 (1996)

    Article  MathSciNet  MATH  Google Scholar 

  9. Kilbas, AA, Trujillo, JJ: Differential equations of fractional order: methods, results and problems-I. Appl. Anal. 78, 153-192 (2001)

    Article  MathSciNet  MATH  Google Scholar 

  10. Kilbas, AA, Trujillo, JJ: Differential equations of fractional order: methods, results and problems-II. Appl. Anal. 81, 435-493 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  11. Agarwal, RP, O’Regan, D, Stanek, S: Positive solutions for Dirichlet problems of singular nonlinear fractional differential equations. J. Math. Anal. Appl. 371, 57-68 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  12. Bai, Z: On solutions of some fractional m-point boundary value problems at resonance. Electron. J. Qual. Theory Differ. Equ. 2010, Article ID 37 (2010)

    Google Scholar 

  13. Leibenson, LS: General problem of the movement of a compressible fluid in a porous medium. Izv. Akad. Nauk Kirg. SSSR 9, 7-10 (1983) (in Russian)

    MathSciNet  Google Scholar 

  14. Chen, T, Liu, W, Yang, C: Antiperiodic solutions for Liénard-type differential equation with p-Laplacian operator. Bound. Value Probl. 2010, Article ID 194824 (2010)

    Article  MathSciNet  Google Scholar 

  15. Jiang, D, Gao, W: Upper and lower solution method and a singular boundary value problem for the one-dimensional p-Laplacian. J. Math. Anal. Appl. 252, 631-648 (2000)

    Article  MathSciNet  MATH  Google Scholar 

  16. Lian, L, Ge, W: The existence of solutions of m-point p-Laplacian boundary value problems at resonance. Acta Math. Appl. Sin. 28, 288-295 (2005)

    MathSciNet  Google Scholar 

  17. Liu, B, Yu, J: On the existence of solutions for the periodic boundary value problems with p-Laplacian operator. J. Syst. Sci. Math. Sci. 23, 76-85 (2003)

    MATH  Google Scholar 

  18. Pang, H, Ge, W, Tian, M: Solvability of nonlocal boundary value problems for ordinary differential equation of higher order with a p-Laplacian. Comput. Math. Appl. 56, 127-142 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  19. Su, H, Wang, B, Wei, Z, Zhang, X: Positive solutions of four-point boundary value problems for higher-order p-Laplacian operator. J. Math. Anal. Appl. 330, 836-851 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  20. Zhang, J, Liu, W, Ni, J, Chen, T: Multiple periodic solutions of p-Laplacian equation with one-side Nagumo condition. J. Korean Math. Soc. 45, 1549-1559 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  21. Podlubny, I: Fractional Differential Equations. Academic Press, San Diego (1999)

    MATH  Google Scholar 

  22. Samko, SG, Kilbas, AA, Marichev, OI: Fractional Integrals and Derivatives: Theory and Applications. Gordon and Breach, Switzerland (1993)

    Google Scholar 

  23. Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam (2006)

    Google Scholar 

  24. Lian, H: Boundary value problems for nonlinear ordinary differential equations on infinite intervals, Doctoral thesis (2007)

  25. Gaines, R, Mawhin, J: Coincidence Degree and Nonlinear Differential Equations. Springer, Berlin (1977)

    MATH  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the anonymous referee for his/her valuable comments, which have improved the correctness and presentation of the manuscript. This work is supported by the National Natural Science Foundation of China (11271364).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hua Jin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All the authors contributed equally and significantly in writing this article. All the authors 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

Jin, H., Liu, W. On the periodic boundary value problem for Duffing type fractional differential equation with p-Laplacian operator. Bound Value Probl 2015, 144 (2015). https://doi.org/10.1186/s13661-015-0408-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0408-3

MSC

Keywords