Skip to main content

Periodic solutions of Rayleigh equations with singularities

Abstract

In this paper, we study the existence of periodic solutions of Rayleigh equations with singularities \(x''+f(t, x')+g(x)=p(t)\). By using the limit properties of the time map, we prove that the given equation has at least one 2π periodic solution.

1 Introduction

In this paper, we are concerned with the existence of periodic solutions of singular Rayleigh equations

$$ x''+f\bigl(t, x'\bigr)+g(x)=p(t), $$
(1.1)

where \(g: (0, +\infty)\to\mathbf{R}\) is continuous and has a singularity at the origin, \(f: \mathbf{R^{2}}\to\mathbf{R^{2}}\) is continuous and 2π periodic with respect to the first variable t, \(p: \mathbf{R}\to\mathbf{R}\) is continuous and 2π periodic.

Equation (1.1) can be used to model the oscillations of a clarinet reed [1]. The dynamic behaviors of (1.1) have been widely investigated due to their applications in many fields such as physics, mechanics, and the engineering technique fields (see [28] and the references therein). Recently, the periodic problem of equations with singularities has been studied widely because of their background in applied sciences (see [913] and the references therein).

When \(f\equiv0\), (1.1) is a conservation system

$$ x''+g(x)=p(t). $$
(1.2)

Assume that g satisfies

$$(\mathrm{h}_{1}) \quad\lim_{x\to0^{+}}g(x)=-\infty, $$

and

$$\frac{n^{2}}{4}< \liminf_{x\to+\infty}\frac{g(x)}{x}\leq\limsup _{x\to +\infty}\frac{g(x)}{x}< \frac{(n+1)^{2}}{4}; $$

moreover, the primitive function G of g satisfies

$$(\mathrm{h}_{2})\quad \lim_{x\to0^{+}}G(x)=+\infty, $$

where \(G(x)=\int_{1}^{x}g(u)\,du\). It was proved in [10] that (1.2) has at least one 2π periodic solution.

It is well known that time maps play an important role in studying the existence and multiplicity of periodic solutions of (1.2) [14, 15]. Assume that g satisfies

$$(\mathrm{h}_{3}) \quad\lim_{x\to+\infty}g(x)=+\infty. $$

Condition (h3) implies that there exists a constant \(d>0\) such that

$$ g(x)>\sup\bigl\{ \bigl|f(t, 0)\bigr|+\bigl|p(t)\bigr|:t\in[0, 2\pi]\bigr\} , \quad\mbox{for } x\geq d. $$
(1.3)

Let us consider the autonomous system

$$x''+g(x)=0, $$

or its equivalent system

$$ x'=y,\qquad y'=-g(x). $$
(1.4)

The first integral of (1.4) is the curve

$$\Gamma_{c}: \frac{1}{2}y^{2}+G(x)=G(c), $$

where c is an arbitrary constant. From conditions (h i ) (\(i=1,2,3\)) we know that, for \(c>0\) sufficiently large, \(\Gamma_{c}\) is a closed curve. Let \((x(t), y(t))\) be any solution of (1.4) whose orbit is \(\Gamma_{c}\). Clearly, this solution is periodic. Let \(T(c)\) denote the least positive period of this solution. It is not hard to calculate

$$T(c)=\sqrt{2}\int_{d(c)}^{c}\frac{dx}{\sqrt{G(c)-G(x)}}, $$

where \(0< d(c)< c\), \(G(d(c))=G(c)\), \(\lim_{c\to+\infty}d(c)=0\). From [10] we know that, if conditions (h i ) (\(i=1,2,3\)) hold, then

$$\lim_{c\to+\infty}\int_{d(c)}^{1} \frac{dx}{\sqrt{G(c)-G(x)}}=0. $$

Now, let us set

$$ \tau(c)=\sqrt{2}\int_{1}^{c}\frac{dx}{\sqrt{G(c)-G(x)}}. $$
(1.5)

In this paper, we deal with the existence of periodic solutions of (1.1) by using the asymptotic properties of the time map τ. Assume that the limit

$$(\mathrm{h}_{4})\quad\lim_{|x|\to+\infty}\frac{f(t,x)}{x}=0 $$

holds uniformly for \(t\in[0, 2\pi]\). We obtain the following result.

Theorem 1.1

Assume that conditions (h i ) (\(i=1,2,3,4\)) hold. Then (1.1) possesses at least one 2π periodic solution provided that the inequality

$$\limsup_{c\to+\infty}\tau(c)>2\pi $$

holds.

Using Theorem 1.1, we can obtain the following corollary.

Corollary 1.2

Assume that conditions (h i ) (\(i=1,2,3,4\)) hold. Then (1.1) possesses at least one 2π periodic solution provided that the inequality

$$\liminf_{x\to+\infty}\frac{2G(x)}{x^{2}}< \frac{1}{4} $$

holds.

Throughout this paper, we always use the notations:

$$\begin{aligned}& \|x\|_{\infty}=\max\bigl\{ \bigl|x(t)\bigr|:t\in[0, 2\pi]\bigr\} ,\qquad \|x \|_{1}=\int_{0}^{2\pi}\bigl|x(t)\bigr|\,dt,\\& \|x \|_{2}=\biggl(\int_{0}^{2\pi}x^{2}(t)\,dt \biggr)^{\frac{1}{2}} \end{aligned}$$

for any continuous 2π periodic function \(x(t)\). For a function \(I(c, \cdot)\), the notation \(I=o(1)\) means that, for \(c\to+\infty\), \(I\to0\) holds uniformly with respect to the other variables.

2 A continuation lemma

It is well known that the continuation theorem plays a key role in studying the existence of periodic solutions of ordinary differential equations. Now we shall introduce a continuation lemma for (1.1). To this end, we consider the equivalent system of (1.1),

$$ x' = y,\qquad y' =-\bigl(g(x) + f (t, y)-p(t)\bigr). $$
(2.1)

Now, we embed system (2.1) into a family of equations with one parameter \(\lambda\in[0, 1]\),

$$ x' =\lambda y,\qquad y' =-\lambda\bigl(g(x) + f (t,\lambda y)-p(t)\bigr). $$
(2.2)

Lemma 2.1

Assume that conditions (h i ) (\(i=1,2,3,4\)) hold. Suppose that there exists a constant \(\zeta\geq d\) (d is given in (1.3)) such that, if \((x(t), y(t))\) is a 2π-periodic solution of system (2.2) for some \(\lambda\in(0, 1)\), then

$$\max\bigl\{ x(t):t\in[0, 2\pi]\bigr\} \neq\zeta,\quad t\in[0, 2\pi]. $$

Then system (2.1) has at least one 2π-periodic solution.

We shall use a classical consequence of Mawhin’s continuation theorem [16], Theorem 7.2 to prove Lemma 2.1. For the reader’s convenience, we restate it here.

Lemma 2.2

Let \(\Psi=\Psi(t, z; \lambda):[0, 2\pi]\times\mathbf{R}^{m}\times[0, 1]\to\mathbf{R}^{m}\) be a continuous function and let \(\Omega\subset\mathbf{R}^{m}\) be a (non-empty) open bounded set (with boundary Ω and closure Ω̄). Assume the following conditions:

  1. (1)

    for any 2π-periodic solution \(z(t)\) of \(z'=\lambda\Psi(t, z; \lambda)\) with \(\lambda\in(0, 1)\), such that \(z(t)\in\bar{\Omega}\), for all \(t\in[0, 2\pi]\), it follows that \(z(t)\in\Omega\), for all \(t\in[0, 2\pi]\);

  2. (2)

    \(\Psi_{0}(z)\neq0\), for each \(z\in\partial\Omega\) and \(d_{B}(\Psi_{0}, \Omega, 0)\neq0\), where

    $$\Psi_{0}(z)=\frac{1}{2\pi}\int_{0}^{2\pi} \Psi(t, z; 0)\,dt,\quad \textit{for } z\in\mathbf{R}^{m}. $$

Then the equation \(z'=\Psi(t, z; 1)\) has at least one 2π-periodic solution and \(z(t)\in\bar{\Omega}\), for all \(t\in[0, 2\pi]\).

Proof of Lemma 2.1

We shall use Lemma 2.2 to prove this continuation lemma. Set

$$\bar{f}=\frac{1}{2\pi}\int_{0}^{2\pi}f(t,0)\,dt,\qquad \bar {p}=\frac{1}{2\pi}\int_{0}^{2\pi}p(t)\,dt. $$

Then there exists \(\tilde{t}\in[0, 2\pi]\) such that

$$\bar{f}-\bar{p}=\frac{1}{2\pi}\int_{0}^{2\pi} \bigl(f(t, 0)-p(t)\bigr)\,dt=f(\tilde{t}, 0)-p(\tilde{t}). $$

From condition (h1) we know that there exists a constant \(0< d_{0}< d\) such that

$$g(x)< -\sup\bigl\{ \bigl|f(t, 0)\bigr|+\bigl|p(t)\bigr|:t\in\mathbf{R}\bigr\} , \quad 0< x\leq d_{0}. $$

Therefore, we have

$$g(x)< -\bigl(\bigl|f(\tilde{t}, 0)\bigr|+\bigl|p(\tilde{t})\bigr|\bigr) \leq-\bigl|f(\tilde{t}, 0)-p(\tilde{t})\bigr|=-|\bar{f}-\bar{p}|,\quad 0< x\leq d_{0}. $$

Meanwhile, we have

$$g(x)>\bigl|f(\tilde{t}, 0)\bigr|+\bigl|p(\tilde{t})\bigr|\geq \bigl|f(\tilde{t}, 0)-p(\tilde{t})\bigr|=|\bar{f}-\bar{p}|,\quad x\geq d. $$

We claim that there exist constants \(0<\varepsilon<d_{0}\) and \(c>0\) such that, if \((x(t), y(t))\) is a 2π periodic solution of (2.2) with \(x(t)\leq\zeta\), \(t\in[0, 2\pi]\), then

$$\varepsilon< x(t)< \zeta, \qquad -c\leq y(t)\leq c,\quad t\in[0, 2\pi]. $$

Integrating the second equality of (2.2) on \([0, 2\pi]\) and applying the first equality of (2.2), we get

$$\int_{0}^{2\pi}g\bigl(x(t)\bigr)\,dt =-\int _{0}^{2\pi}f\bigl(t, x'(t)\bigr)\,dt+\int _{0}^{2\pi}p(t)\,dt. $$

Then we obtain

$$-\int_{I_{1}}g\bigl(x(t)\bigr)\,dt=\int_{I_{2}}g \bigl(x(t)\bigr)\,dt+\int_{0}^{2\pi}f\bigl(t, x'(t)\bigr)\,dt-\int_{0}^{2\pi}p(t)\,dt, $$

where \(I_{1}=\{t\in[0, 2\pi]: 0< x(t)< d_{0}\}\), \(I_{2}=\{t\in[0, 2\pi]: d_{0}\leq x(t)\leq\zeta\}\). Hence, we have

$$\begin{aligned} \int_{0}^{2\pi}\bigl|g\bigl(x(t)\bigr)\bigr|\,dt=&- \int_{I_{1}}g\bigl(x(t)\bigr)\,dt+ \int_{I_{2}}\bigl|g\bigl(x(t)\bigr)\bigr|\,dt \\ \leq& 2 \int_{I_{2}}\bigl|g\bigl(x(t)\bigr)\bigr|\,dt+ \int_{0}^{2\pi}\bigl|f\bigl(t, x'(t)\bigr)\bigr|\,dt+ \int_{0}^{2\pi}\bigl|p(t)\bigr|\,dt \\ \leq& 2M_{1}+ \int_{0}^{2\pi}\bigl|f\bigl(t, x'(t)\bigr)\bigr|\,dt+\|p\|_{1}, \end{aligned}$$
(2.3)

where \(M_{1}=2\pi\cdot\max\{|g(x)|:d_{0}\leq x\leq\zeta\}\).

Let us take a fixed constant δ satisfying \(0<(1+\frac{\pi}{\sqrt{3}})\delta<1\). From (h4) we see that there exists \(R_{\delta}>0\) such that, for any \(|s|\geq R_{\delta}\) and \(t\in[0, 2\pi]\),

$$\bigl|f(t,s)\bigr|\leq\delta|s|. $$

Set

$$M_{2}=\max\bigl\{ \bigl|f(t,s)\bigr|: t\in[0, 2\pi], |s|\leq R_{\delta} \bigr\} . $$

Then we see that, for any \((t, s)\in\mathbf{R}^{2}\),

$$ \bigl|f(t, s)\bigr|\leq\delta|s|+M_{2}. $$
(2.4)

From (2.3) and (2.4) we get

$$\int_{0}^{2\pi}\bigl|g\bigl(x(t)\bigr)\bigr|\,dt \leq2M_{1}+\delta\int_{0}^{2\pi}\bigl|x'(t)\bigr|\,dt+2 \pi M_{2}+\|p\|_{1}. $$

Set \(M=2M_{1}+2\pi M_{2}+\|p\|_{1}\). Then we obtain

$$ \int_{0}^{2\pi}\bigl|g\bigl(x(t)\bigr)\bigr|\,dt\leq\delta \sqrt{2\pi}\bigl\| x'\bigr\| _{2}+M. $$
(2.5)

From (2.2) we know that \(x(t)\) satisfies the equation as follows:

$$ x''(t)+\lambda^{2}\bigl(f\bigl(t, x'(t)\bigr)+g\bigl(x(t)\bigr)-p(t)\bigr)=0. $$
(2.6)

Multiplying (2.6) by \(x(t)-\bar{x}\) with \(\bar{x} =\frac{1}{2\pi}\int_{0}^{2\pi}x(s)\,ds\), and integrating the equality on \([0, 2\pi]\), we get from (2.5) and (2.6)

$$\begin{aligned} \int_{0}^{2\pi}\bigl|x'(t)\bigr|^{2}\,dt={}& \lambda^{2}\int_{0}^{2\pi}f\bigl(t, x'(t)\bigr) \bigl(x(t)-\bar{x}\bigr)\,dt+\lambda^{2}\int _{0}^{2\pi}g\bigl(x(t)\bigr) \bigl(x(t)-\bar {x} \bigr)\,dt\\ &{}-\lambda^{2}\int_{0}^{2\pi}p(t) \bigl(x(t)-\bar{x}\bigr)\,dt \\ \leq{}& \int_{0}^{2\pi}\bigl|f\bigl(t, x'(t) \bigr) \bigl(x(t)-\bar{x}\bigr)\bigr|\,dt+\int_{0}^{2\pi}\bigl|g \bigl(x(t)\bigr) \bigl(x(t)-\bar{x}\bigr)\bigr|\,dt\\ &{}+\int_{0}^{2\pi}\bigl|p(t) \bigl(x(t)-\bar{x}\bigr)\bigr|\,dt \\ \leq{}& \delta\int_{0}^{2\pi}\bigl|x'(t)\bigr|\bigl|x(t)- \bar{x}\bigr|\,dt+M_{2}\int_{0}^{2\pi}\bigl|x(t)- \bar{x}\bigr|\,dt \\ &{}+\|x-\bar{x}\|_{\infty}\int_{0}^{2\pi} \bigl(\bigl|g\bigl(x(t)\bigr)\bigr|+\bigl|p(t)\bigr|\bigr)\,dt \\ \leq{}& \bigl(\delta \bigl\| x'\bigr\| _{2}+M_{2}\sqrt{2\pi} \bigr)\|x-\bar{x}\|_{2}+\bigl(\delta\sqrt{2\pi }\bigl\| x' \bigr\| _{2}+M+\|p\|_{1}\bigr)\|x-\bar{x}\|_{\infty}. \end{aligned}$$

Using the Wirtinger inequality and the Sobolev inequality, we have

$$\|x-\bar{x}\|_{2}\leq\bigl\| x'\bigr\| _{2}, \qquad\|x-\bar{x} \|_{\infty}\leq \sqrt{\frac{\pi}{6}}\bigl\| x'\bigr\| _{2}. $$

Then we get

$$\begin{aligned}\int_{0}^{2\pi}\bigl|x'(t)\bigr|^{2}\,dt \leq \delta\biggl(1+\frac{\pi}{\sqrt{3}}\biggr)\bigl\| x' \bigr\| _{2}^{2}+ \biggl(M_{2}\sqrt{2\pi}+\sqrt { \frac{\pi}{6}}\bigl(M+\|p\|_{1}\bigr) \biggr)\bigl\| x' \bigr\| _{2},\end{aligned}$$

which means that

$$\bigl\| x'\bigr\| _{2}\leq\gamma\bigl\| x' \bigr\| _{2}+c_{0}, $$

where \(\gamma=\delta(1+\frac{\pi}{\sqrt{3}})\), \(c_{0}=M_{2}\sqrt{2\pi}+\sqrt{\frac{\pi}{6}}(M+\|p\|_{1})\). Since \(0<\gamma<1\), we have

$$ \bigl\| x'\bigr\| _{2}\leq\frac{c_{0}}{1-\gamma}:=c_{1}. $$
(2.7)

Integrating the first equation of (2.2) on \([0, 2\pi]\) and noticing \(\lambda\in(0, 1]\), we get

$$\int_{0}^{2\pi}y(t)\,dt=0, $$

which implies that there exists \(t_{0}\in[0, 2\pi]\) such that \(y(t_{0}) = 0\). Then we get from (2.4), (2.5), and (2.7)

$$\begin{aligned} \bigl|y(t)\bigr|&\leq\bigl|y(t_{0})\bigr|+\int_{0}^{2\pi}\bigl|y'(t)\bigr|\,dt \leq \int_{0}^{2\pi}\bigl|g\bigl(x(t)\bigr)\bigr|\,dt+\int _{0}^{2\pi}\bigl|f\bigl(t, x'(t)\bigr)\bigr|\,dt+\int _{0}^{2\pi}\bigl|p(t)\bigr|\,dt \\ &\leq 2\sqrt{2\pi}\delta c_{1}+M+2\pi M_{2}+\|p \|_{1}:=c. \end{aligned}$$

Therefore,

$$ \|y\|_{\infty}\leq c. $$
(2.8)

Let \(x(t_{*})\) (\(t_{*}\in[0, 2\pi]\)) be the minimum of \(x(t)\). Then we have \(x'(t_{*})=0\) and \(x''(t_{*})\geq0\). Since \(x(t_{*})\) satisfies

$$x''(t_{*})+\lambda^{2}\bigl(f(t_{*}, 0)+g \bigl(x(t_{*})\bigr)-p(t_{*})\bigr)=0, $$

we have

$$f(t_{*}, 0)+g\bigl(x(t_{*})\bigr)\leq p(t_{*}). $$

Hence,

$$g\bigl(x(t_{*})\bigr)\leq-f(t_{*}, 0)+p(t_{*})\leq\bigl|f(t_{*}, 0)\bigr|+\bigl|p(t_{*})\bigr| \leq \sup \bigl\{ \bigl|f(t, 0)\bigr|+\bigl|p(t)\bigr|:t\in\mathbf{R}\bigr\} , $$

which implies

$$ x(t_{*})< d. $$
(2.9)

Let \(x(t^{*})\) (\(t^{*}\in[0, 2\pi]\)) be the maximum of \(x(t)\). Then we have \(x'(t^{*})=0\) and \(x''(t^{*})\leq0\). Similarly, we can obtain

$$ x\bigl(t^{*}\bigr)>d_{0}. $$
(2.10)

From (2.9) and (2.10) we see that there exists \(\bar{t}\in[0, 2\pi]\) such that

$$ d_{0}\leq x(\bar{t})\leq d. $$
(2.11)

In what follows, we shall prove that there exists \(0<\varepsilon<d_{0}\) such that, for any 2π periodic solution \((x(t), y(t))\) of (2.2) with \(x(t)\leq\zeta\), \(t\in[0, 2\pi]\),

$$\varepsilon< x(t)< \zeta,\quad t\in[0, 2\pi]. $$

The right inequality \(x(t)<\zeta\) (\(t\in[0, 2\pi]\)) follows directly from the condition \(\max\{x(t):t\in[0, 2\pi]\}\neq\zeta\) and \(x(t)\leq\zeta\), \(t\in[0, 2\pi]\). Next, we prove the left inequality. Otherwise, there exist a sequence \(\{\lambda_{n}\}\) with \(\lambda_{n}\in(0, 1]\) and a sequence of 2π periodic solutions of (2.2) \(\{(x_{n}(t), y_{n}(t))\}\) (with \(\lambda=\lambda_{n}\) in (2.2)), satisfying \(x_{n}(t)\leq\zeta\), \(t\in[0, 2\pi]\), and

$$\min_{t\in[0, 2\pi]}x_{n}(t)\to0,\quad n\to\infty. $$

Without loss of generality, we assume that, for every n,

$$ \min_{t\in[0, 2\pi]}x_{n}(t)< d_{0}. $$
(2.12)

Set \(\varepsilon_{n}=x_{n}(t_{n})=\min_{t\in[0, 2\pi]}x_{n}(t)\), \(t_{n}\in [0, 2\pi]\). From (2.11) and (2.12) we see that there exists \(\alpha_{n}\in(t_{n}, t_{n}+2\pi)\) such that

$$x(\alpha_{n})=d_{0},\qquad \varepsilon_{n}< x_{n}(t)< d_{0},\quad t\in(t_{n}, \alpha_{n}). $$

Since \((x_{n}(t), y_{n}(t))\) satisfies the equation

$$y_{n}'(t) =-\lambda_{n} \bigl(g \bigl(x_{n}(t)\bigr) + f \bigl(t,\lambda_{n} y_{n}(t)\bigr)-p(t)\bigr), $$

we have

$$y_{n}(t)y_{n}'(t) =-\lambda_{n} y_{n}(t) \bigl(g\bigl(x_{n}(t)\bigr) + f \bigl(t, \lambda_{n} y_{n}(t)\bigr)-p(t)\bigr). $$

Recalling \(x_{n}'(t)=\lambda_{n} y_{n}(t)\), we get

$$ y_{n}(t)y_{n}'(t) =-\bigl(g \bigl(x_{n}(t)\bigr)+f \bigl(t,x_{n}'(t) \bigr)-p(t)\bigr)x_{n}'(t). $$
(2.13)

Integrating both sides of (2.13) over the interval \([t_{n}, \alpha_{n}]\) and using the fact \(x_{n}'(t_{n})=\lambda_{n} y_{n}(t_{n})=0\), we obtain

$$\frac{1}{2}y_{n}^{2}(\alpha_{n})=-\int _{\varepsilon_{n}}^{d_{1}}g(s)\,ds-\int_{t_{n}}^{\alpha_{n}}f \bigl(t, x_{n}'(t)\bigr)x_{n}'(t)\,dt+ \int_{t_{n}}^{\alpha_{n}}p(t)x_{n}'(t)\,dt. $$

Therefore, we get

$$\begin{aligned} -\int_{\varepsilon_{n}}^{d_{1}}g(s)\,ds={}&\biggl\vert \int _{\varepsilon_{n}}^{d_{1}}g(s)\,ds\biggr\vert \leq \frac{1}{2}y_{n}^{2}(\alpha_{n})+\int _{0}^{2\pi}\bigl\vert f\bigl(t, x_{n}'(t) \bigr)x_{n}'(t)\bigr\vert \,dt \\ &{}+\int_{0}^{2\pi} \bigl\vert p(t)x_{n}'(t)\bigr\vert \,dt. \end{aligned}$$
(2.14)

From (h2) we have

$$ \lim_{n\to\infty}\int_{\varepsilon_{n}}^{d_{1}}g(s)\,ds=- \infty. $$
(2.15)

Next, we shall estimate the right hand side of (2.14). First, it follows from (2.8) that we have

$$y_{n}^{2}(\alpha_{n})\leq c^{2}. $$

Meanwhile, according to (2.4) and (2.7), we get

$$\begin{aligned} \int_{0}^{2\pi}\bigl\vert f\bigl(t, x_{n}'(t)\bigr)x_{n}'(t)\bigr\vert \,dt&\leq \delta\int_{0}^{2\pi}x_{n}^{\prime2}(t)\,dt+M_{2} \int_{0}^{2\pi}\bigl|x_{n}'(t)\bigr|\,dt \\ &\leq \delta\bigl\| x_{n}'\bigr\| _{2}^{2}+ \sqrt{2\pi}M_{2}\bigl\| x_{n}'\bigr\| _{2} \\ &\leq \delta c_{1}^{2}+\sqrt{2\pi}M_{2}c_{1}. \end{aligned}$$

Obviously, we have

$$\int_{0}^{2\pi}\bigl\vert p(t)x_{n}'(t) \bigr\vert \,dt\leq \sqrt{2\pi}\|p\|_{\infty}\bigl\| x_{n}' \bigr\| _{2}\leq \sqrt{2\pi}\|p\|_{\infty}c_{1}. $$

Hence, the right hand side of (2.14) is bounded. This conclusion contradicts (2.15).

To use Lemma 2.2, we define an open bounded set \(\Omega=\{(x, y): \varepsilon< x<\zeta, -c-1<y<c+1\}\), and a map \(S: (0, +\infty)\times \mathbf{R}\to\mathbf{R}^{2}\), \(S(x, y)=(y, -g(x)-\bar{f}+\bar{p})\). Then, for any 2π-periodic solution \((x(t), y(t))\) of system (2.2), such that \((x(t), y(t))\in\bar{\Omega}\), for all \(t\in[0, 2\pi]\), we have \((x(t), y(t))\in\Omega\), for all \(t\in[0, 2\pi]\). Therefore, the first condition of Lemma 2.2 is satisfied. Obviously, S does not vanish outside the rectangle Ω. Furthermore, the Brouwer degree of S, \(d_{B}(S, \Omega, 0)\), is defined and \(d_{B}(S, \Omega, 0)=d_{B}(g, (\varepsilon, \zeta), \bar{p}-\bar{f})=1\) because g is continuous and \(g(\varepsilon)<\bar{p}-\bar{f}\), \(g(\zeta)>\bar{p}-\bar{f}\). According to Lemma 2.2, system (2.1) has at least one 2π periodic solution. □

Lemma 2.3

[14]

Assume that \(g: \mathbf{R}\to \mathbf{R}\) is continuous and \(\lim_{|x|\to+\infty}\operatorname{sgn}(x)g(x)=+\infty\). Then, for any constant \(\nu\in\mathbf{R}\),

$$\lim_{c\to\pm\infty}\frac{\tau_{g}(\nu, c)}{\tau_{g}(c)}=1, $$

where

$$\tau_{g}(c)=2\biggl\vert \int_{0}^{c} \frac{dx}{\sqrt{2(\tilde{G}(c)-\tilde{G}(x))}}\biggr\vert ,\qquad \tau_{g}(\nu, c)=2\biggl\vert \int _{0}^{c}\frac{dx}{\sqrt{2(\tilde{G}(c)-\tilde{G}(x)+\nu(c-x))}}\biggr\vert $$

with \(\tilde{G}(x)=\int_{0}^{x}g(s)\,ds\).

Remark 2.4

When \(g: [0, +\infty)\to\mathbf{R}\) is continuous and satisfies \(\lim_{x\to+\infty}g(x)=+\infty\), we can also define \(\tau_{g}(c)\) and \(\tau_{g}(\nu, c)\) for \(c>0\) large enough. In this case, we know from Lemma 2.3 that, for any constant ν,

$$\lim_{c\to+\infty}\frac{\tau_{g}(\nu, c)}{\tau_{g}(c)}=1. $$

When \(g: (0, +\infty)\to\mathbf{R}\) is continuous and \(\lim_{x\to+\infty}g(x)=+\infty\), we can get a similar estimate. Under this condition, it is noted that g may have a singularity at the origin, \(x=0\), namely, \(\lim_{x\to0^{+}}g(x)=-\infty\). For any constant \(\nu\in\mathbf{R}\) and sufficiently large \(c\geq1\), let us set

$$\tau_{g}^{+}(\nu, c)=2\int_{1}^{c} \frac{dx}{\sqrt{2(G(c)-G(x)+\nu(c-x))}}, $$

where \(G(x)=\int_{1}^{x}g(s)\,ds\). Then we have

$$ \lim_{c\to+\infty}\frac{\tau_{g}^{+}(\nu, c)}{\tau(c)}=1, $$
(2.16)

where τ is defined by (1.5).

In fact, let us consider a function \(g_{0}: [0, +\infty)\to\mathbf{R}\), \(g_{0}(x)=g(x+1)\), \(x\geq0\). Obviously, \(g_{0}\) is continuous on the interval \([0, +\infty)\) and satisfies \(\lim_{x\to+\infty}g_{0}(x)=+\infty\). Then we have, for \(x\geq0\),

$$\tilde{G}_{0}(x)=\int_{0}^{x} g_{0}(s)\,ds=\int_{1}^{x+1}g(s)\,ds=G(x+1). $$

According to Lemma 2.3, we get

$$\lim_{c\to+\infty}\frac{\tau_{g_{0}}(\nu, c)}{\tau_{g_{0}}(c)}=1. $$

When \(c>0\) is large enough, we have

$$\begin{aligned} \tau_{g_{0}}(\nu, c)&= 2\int_{0}^{c} \frac{dx}{\sqrt{2(\tilde{G}_{0}(c)-\tilde{G}_{0}(x)+\nu (c-x))}} \\ &= 2\int_{0}^{c}\frac{dx}{\sqrt{2(G(c+1)-G(x+1)+\nu(c-x))}} \\ &= 2\int_{1}^{c+1}\frac{dx}{\sqrt{2(G(c+1)-G(x)+\nu(c+1-x))}}= \tau_{g}^{+}(\nu, c+1). \end{aligned}$$

Similarly, we have

$$\tau_{g_{0}}(c)=\tau(c+1). $$

Consequently, we get

$$\lim_{c\to+\infty}\frac{\tau_{g}^{+}(\nu, c+1)}{\tau(c+1)}=1. $$

Therefore, the conclusion (2.16) holds.

3 Proof of Theorem 1.1

In this section, we shall use the continuation Lemma 2.1 given in Section 2 to prove Theorem 1.1.

Proof of Theorem 1.1

Let us set

$$\tau=\limsup_{c\to+\infty}\tau(c)>2\pi. $$

Then there exist \(0<\varepsilon_{0}<\frac{1}{3}(\tau-2\pi)\) and a sequence \(\{c_{n}\}\) with \(\lim_{n\to\infty}c_{n}=+\infty\) such that, for every n,

$$\tau(c_{n})>\tau-\varepsilon_{0}>2\pi+2 \varepsilon_{0}. $$

We shall prove that the condition of Lemma 2.1 is satisfied for \(\zeta=c_{n}\) with n sufficiently large.

Let \((x(t), y(t))\) be any 2π periodic solution of (2.2) for some \(\lambda\in(0, 1]\) and suppose that, for n large enough,

$$x\bigl(t^{*}\bigr)=\max_{t\in[0, 2\pi]}x(t) =c_{n}>d, $$

where d is given in (1.3). Assume that \(x(t_{*})\) (\(t_{*}\in[0, 2\pi]\)) is a local minimum of \(x(t)\). From the proof of Lemma 2.1

$$x(t_{*})< d. $$

Then there exists an interval \([\alpha, \beta]\subset[0, 2\pi]\) containing \(t^{*}\), with \(\alpha=\alpha(x, \lambda)\), \(\beta=\beta(x, \lambda)\) such that

$$x(\alpha)=x(\beta)=d;\qquad x(t)>d, \quad t\in(\alpha, \beta) $$

and

$$y\bigl(t^{*}\bigr)=0;\qquad y(t)>0,\quad t\in\bigl[\alpha, t^{*}\bigr), \qquad y(t)< 0,\quad t\in \bigl(t^{*}, \beta \bigr]. $$

From (2.2) we have

$$ y(t)y'(t)+\lambda\bigl(f\bigl(t, \lambda y(t)\bigr)+g\bigl(x(t) \bigr)-p(t)\bigr)y(t)=0. $$
(3.1)

Integrating both sides of (3.1) on the interval \([t, t^{*}]\) with \(\alpha\leq t\leq t^{*}\), we have

$$ y^{2}(t)=2\bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(t) \bigr)\bigr)+2\lambda\int_{t}^{t^{*}}f\bigl(\tau, \lambda y(\tau)\bigr)y(\tau)\,d\tau-2\lambda\int_{t}^{t^{*}}p( \tau)y(\tau)\,d\tau. $$
(3.2)

From (h4) we know that, for any sufficiently small \(\varepsilon>0\), there is a constant \(M_{\varepsilon}>0\) such that, for any \((t, y)\in R^{2}\),

$$ \bigl|f(t, y)\bigr|\leq\varepsilon|y|+M_{\varepsilon}. $$
(3.3)

Since \(y(t)>0\), \(t\in[\alpha, t^{*}]\), it follows from (3.2) and (3.3) that, for \(t\in[\alpha, t^{*}]\),

$$\begin{aligned} y^{2}(t)\leq& 2\bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G \bigl(x(t)\bigr)\bigr)+2\int_{t}^{t^{*}}\bigl|f\bigl(\tau, \lambda y(\tau)\bigr)\bigr|\bigl|\lambda y(\tau)\bigr|\,d\tau+2\int_{t}^{t^{*}}\bigl|p( \tau)\bigr|\bigl|\lambda y(\tau)\bigr|\,d\tau \\ &\leq 2\bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(t)\bigr)\bigr)+2\varepsilon \int_{t}^{t^{*}}y^{2}(\tau)\,d\tau +2\int _{t}^{t^{*}}\bigl(\bigl|p(\tau)\bigr|+M_{\varepsilon}\bigr)\bigl| \lambda y(\tau)\bigr|\,d\tau \\ &\leq 2\bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(t)\bigr)\bigr)+2\varepsilon \int_{t}^{t^{*}}y^{2}(\tau)\,d\tau +M_{\varepsilon}'\int_{t}^{t^{*}}x'( \tau)\,d\tau \\ &= 2\bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(t)\bigr)\bigr)+2\varepsilon\int _{t}^{t^{*}}y^{2}(\tau)\,d\tau +M_{\varepsilon}'\bigl(x\bigl(t^{*}\bigr)-x(t)\bigr), \end{aligned}$$

where \(M_{\varepsilon}'=M_{\varepsilon}+\|p\|_{\infty}\). Let us set

$$\phi(t)=\int_{t}^{t^{*}}y^{2}(\tau)\,d\tau. $$

Then we have

$$\phi'(t)=-y^{2}(t). $$

Hence,

$$ -\phi'(t)-2\varepsilon\phi(t)\leq 2\bigl(G\bigl(x\bigl(t^{*}\bigr) \bigr)-G\bigl(x(t)\bigr)\bigr)+M_{\varepsilon}'\bigl(x\bigl(t^{*} \bigr)-x(t)\bigr). $$
(3.4)

Multiplying both sides of (3.4) by \(e^{2\varepsilon t}\) and integrating over the interval \([t, t^{*}]\) yields

$$-\int_{t}^{t^{*}}\bigl[\phi(\tau)e^{2\varepsilon\tau} \bigr]'\,d\tau\leq \int_{t}^{t^{*}}\bigl[2 \bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(\tau)\bigr)\bigr)+M_{\varepsilon}' \bigl(x\bigl(t^{*}\bigr)-x(\tau )\bigr)\bigr]e^{2\varepsilon \tau}\,d\tau. $$

Since \(\phi(t^{*})=0\), we have

$$\phi(t)e^{2\varepsilon t}\leq\int_{t}^{t^{*}}\bigl[2 \bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(\tau)\bigr)\bigr)+M_{\varepsilon }' \bigl(x\bigl(t^{*}\bigr)-x(\tau)\bigr)\bigr]e^{2\varepsilon \tau}\,d\tau. $$

From \(x'(t)=\lambda y(t)\geq0\), \(t\in[\alpha, t^{*}]\) we know that \(x(t)\) is increasing on the interval \([\alpha, t^{*}]\). Therefore, we get, for \(t\in[\alpha, t^{*}]\),

$$\phi(t)e^{2\varepsilon t}\leq e^{2\varepsilon t^{*}}\int_{t}^{t^{*}} \bigl[2\bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(\tau)\bigr) \bigr)+M_{\varepsilon}'\bigl(x\bigl(t^{*}\bigr)-x(\tau )\bigr)\bigr]\,d\tau. $$

Furthermore,

$$\phi(t)\leq2\pi e^{4\pi \varepsilon}\bigl[2\bigl(G\bigl(x\bigl(t^{*}\bigr)\bigr)-G \bigl(x(t)\bigr)\bigr)+M_{\varepsilon}'\bigl(x\bigl(t^{*}\bigr)-x(t) \bigr)\bigr]. $$

Consequently, we can get, for \(t\in[\alpha, t^{*}]\),

$$y^{2}(t)\leq \bigl(1+\kappa(\varepsilon)\bigr)\bigl[2\bigl(G\bigl(x \bigl(t^{*}\bigr)\bigr)-G\bigl(x(t)\bigr)\bigr)+M_{\varepsilon}' \bigl(x\bigl(t^{*}\bigr)-x(t)\bigr)\bigr], $$

where \(\kappa(\varepsilon)=4\pi\varepsilon e^{4\pi\varepsilon}\). Recalling \(x'(t)=\lambda y(t)\) and \(y(t)>0\) for \(t\in[\alpha, t^{*}]\), we have

$$x'(t)\leq \sqrt{1+\kappa(\varepsilon)}\sqrt{2\bigl(G \bigl(x\bigl(t^{*}\bigr)\bigr)-G\bigl(x(t)\bigr)\bigr)+M_{\varepsilon }' \bigl(x\bigl(t^{*}\bigr)-x(t)\bigr)}. $$

Hence,

$$ \frac{x'(t)}{\sqrt{1+\kappa(\varepsilon)}\sqrt {2(G(x(t^{*}))-G(x(t)))+M_{\varepsilon}'(x(t^{*})-x(t))}}\leq 1. $$
(3.5)

Integrating both sides of (3.5) over interval \([\alpha, t^{*}]\) yields

$$\frac{1}{\sqrt{1+\kappa(\varepsilon)}}\int_{d_{2}}^{c_{n}} \frac{dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}\leq t^{*}-\alpha. $$

Similarly, we can get

$$\frac{1}{\sqrt{1+\kappa(\varepsilon)}}\int_{d_{2}}^{c_{n}} \frac{dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}\leq \beta-t^{*}. $$

Therefore, we obtain

$$\frac{2}{\sqrt{1+\kappa(\varepsilon)}}\int_{d_{2}}^{c_{n}} \frac{dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}\leq \beta-\alpha. $$

Using (h3) we can easily derive that, for \(n\to\infty\),

$$\int_{1}^{d_{2}} \frac{dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}=o(1). $$

Then we have

$$\frac{2}{\sqrt{1+\kappa(\varepsilon)}}\int_{1}^{c_{n}} \frac{dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}+o(1)\leq \beta-\alpha. $$

It follows from Remark 2.4 that

$$\lim_{n\to\infty}\frac{1}{\tau(c_{n})}\int_{1}^{c_{n}} \frac{2 dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}=1. $$

Consequently, we have

$$\int_{1}^{c_{n}} \frac{2 dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}= \tau(c_{n}) \bigl(1+o(1)\bigr). $$

Furthermore,

$$\int_{1}^{c_{n}} \frac{2dx}{\sqrt{2(G(c_{n})-G(x))+M_{\varepsilon}'(c_{n}-x)}}\geq 2\pi+2 \varepsilon_{0}+o(1). $$

Since \(\lim_{\varepsilon\to0^{+}}\sqrt{1+\kappa(\varepsilon)}=1\), there exist a sufficiently small \(\varepsilon>0\) and a sufficiently large n such that, if \(\max_{[0, 2\pi]}x(t)=c_{n}\), then

$$\beta-\alpha>2\pi+\varepsilon_{0}, $$

which contradicts with the inequality \(\beta-\alpha<2\pi\). Then we find \(\zeta=c_{n}\) for n sufficiently large. Consequently, from the continuation Lemma 2.1, we know that (2.1) has at least one 2π periodic solution. □

Proof of Corollary 1.2

Let us denote \(\rho=\liminf_{x\to+\infty}\frac{2G(x)}{x^{2}}<\frac{1}{4}\). Then there exists \(\varepsilon>0\) such that \(\rho_{\varepsilon}=\rho+\varepsilon\in(\rho, \frac{1}{4})\). Define

$$\psi(x)=\rho_{\varepsilon} x^{2}-2G(x),\quad x\geq1. $$

Therefore, we have

$$\limsup_{x\to+\infty}\psi(x)=\limsup_{x\to+\infty}x^{2} \biggl(\rho_{\varepsilon }-\frac{2G(x)}{x^{2}}\biggr)=+\infty. $$

It follows that there exists a sequence \(\{c_{n}\}\) with \(\lim_{n\to+\infty}c_{n}=+\infty\) such that

$$\psi(x)\leq\psi(c_{n}),\quad x\in(1, c_{n}). $$

Consequently,

$$2\bigl(G(c_{n})-G(x)\bigr)\leq\rho_{\varepsilon} \bigl(c_{n}^{2}-x^{2}\bigr),\quad x\in(1, c_{n}). $$

Hence, we have

$$\tau(c_{n})=2\int_{1}^{c_{n}} \frac{dx}{\sqrt{2(G(c_{n})-G(x))}}\geq 2\int_{1}^{c_{n}} \frac{dx}{\sqrt{\rho_{\varepsilon}(c_{n}^{2}-x^{2})}} = \frac{2}{\sqrt{\rho_{\varepsilon}}}\biggl(\frac{\pi}{2}-\arcsin \frac{1}{c_{n}}\biggr). $$

As a result, we get

$$\limsup_{n\to+\infty}\tau(c_{n})\geq \frac{\pi}{\sqrt{\rho_{\varepsilon}}}>2 \pi, $$

which implies that \(\limsup_{c\to+\infty}\tau(c)>2\pi\). According to Theorem 1.1, (1.1) has at least one 2π periodic solution. □

Remark 3.1

In [12], the existence of periodic solutions of the Hamiltonian systems of the type

$$ x'=g_{1}(t, y),\qquad y'=-g_{2}(t, x) $$
(3.6)

was studied. A similar result was obtained (see [12], Corollary 3.13) for system (3.6). However, this corollary cannot be applied directly to obtain the main results of this paper because the asymptotic behavior of the primitive G of the nonlinearity g is treated in present paper.

References

  1. Radhakrishnan, S: Exact solutions of Rayleigh’s equation and sufficient conditions for inviscid instability of parallel, bounded shear flows. Z. Angew. Math. Phys. 45, 615-637 (1994)

    Article  MATH  MathSciNet  Google Scholar 

  2. Smith, HL: On the small oscillations of the periodic Rayleigh equation. Q. Appl. Math. 44, 223-247 (1986)

    MATH  Google Scholar 

  3. Smith, RA: Period bounds for generalized Rayleigh equation. Int. J. Non-Linear Mech. 6, 271-277 (1977)

    Article  Google Scholar 

  4. Omari, P, Villari, G: On a continuation lemma for the study of a certain planar system with applications to Liénard and Rayleigh equations. Results Math. 14, 156-173 (1998)

    Article  MathSciNet  Google Scholar 

  5. Ma, T, Wang, Z: A continuation lemma and its applications to periodic solutions of Rayleigh differential equations with subquadratic potential conditions. J. Math. Anal. Appl. 385, 1107-1118 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  6. Ma, T: Periodic solutions of Rayleigh equations via time-maps. Nonlinear Anal. 75, 4137-4144 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  7. Habets, P, Torres, PJ: Some multiplicity results for periodic solutions of a Rayleigh differential equation. Dyn. Contin. Discrete Impuls. Syst., Ser. A Math. Anal. 8, 335-351 (2001)

    MATH  MathSciNet  Google Scholar 

  8. Si, J: Invariant tori and unboundedness for second order differential equations with asymmetric nonlinearities depending on the derivatives. Nonlinear Anal. 67, 3098-3115 (2007)

    Article  MATH  MathSciNet  Google Scholar 

  9. Habets, P, Sanchez, L: Periodic solutions of some Liénard equations with singularities. Proc. Am. Math. Soc. 109, 1035-1044 (1990)

    MATH  MathSciNet  Google Scholar 

  10. Wang, Z: Periodic solutions of the second order differential equations with singularities. Nonlinear Anal. 58, 319-331 (2004)

    Article  MATH  MathSciNet  Google Scholar 

  11. Sfecci, A: Nonresonance conditions for radial solutions of nonlinear Neumann elliptic problems on annuli. Rend. Ist. Mat. Univ. Trieste 46, 255-270 (2014)

    MathSciNet  Google Scholar 

  12. Fonda, A, Sfecci, A: A general method for the existence of periodic solutions of differential systems in the plane. J. Differ. Equ. 252, 1369-1391 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  13. Fonda, A, Toader, R: Periodic orbits of radially symmetric Keplerian-like systems: a topological degree approach. J. Differ. Equ. 244, 3235-3264 (2008)

    Article  MATH  MathSciNet  Google Scholar 

  14. Opial, Z: Sur les périodes des solutions de l’équation différentielle \(x''+g(x)=0\). Ann. Pol. Math. 10, 49-72 (1961)

    MATH  MathSciNet  Google Scholar 

  15. Fonda, A, Zanolin, F: On the use of time-maps for the solvability of nonlinear boundary value problems. Arch. Math. 59, 245-259 (1992)

    Article  MATH  MathSciNet  Google Scholar 

  16. Mawhin, J: Equivalence theorems for nonlinear operator equations and coincidence degree theory for some mappings in locally convex topological vector spaces. J. Differ. Equ. 12, 610-636 (1972)

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors are grateful to the referees for many valuable suggestions to make the paper more readable. Research supported by National Nature Science Foundation of China, No. 11501381 and the Grant of Beijing Education Committee Key Project, No. KZ201310028031.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zaihong Wang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ZW proved a continuation lemma for Rayleigh equations. TM participated in obtaining a prior estimate and helped to draft the manuscript. All 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

Wang, Z., Ma, T. Periodic solutions of Rayleigh equations with singularities. Bound Value Probl 2015, 154 (2015). https://doi.org/10.1186/s13661-015-0427-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0427-0

MSC

Keywords