\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2015 (2015), No. 204, pp. 1--17.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2015 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2015/204\hfil Impulsive Hamiltonian systems]
{Nontrivial periodic solutions to second-order impulsive Hamiltonian systems}

\author[J. R. Graef, S. Heidarkhani, L. Kong \hfil EJDE-2015/204\hfilneg]
{John R. Graef, Shapour Heidarkhani, Lingju Kong}

\address{John R. Graef \newline
Department of Mathematics,
University of Tennessee at Chattanooga,
Chattanooga, TN 37403, USA}
\email{John-Graef@utc.edu}

\address{Shapour Heidarkhani \newline
Department of Mathematics,
Faculty of Sciences, Razi University, Kermanshah 67149, Iran}
\email{s.heidarkhani@razi.ac.ir}

\address{Lingju Kong \newline
Department of Mathematics,
University of Tennessee at Chattanooga,
Chattanooga, TN 37403, USA}
 \email{Lingju-Kong@utc.edu}

\thanks{Submitted June 26, 2015. Published August 10, 2015.}
\subjclass[2010]{34B15, 47J10}
\keywords{Nontrivial periodic solution; impulsive Hamiltonian system;
\hfill\break\indent critical point theory; mountain pass theorem; 
 variational methods}

\begin{abstract}
 Based on variational methods and critical point theory,
 we study the existence of nontrivial periodic
 solutions to a class of second-order impulsive Hamiltonian
 systems. A unique feature of the approach used here is that we use
 a combination of techniques to obtain the existence of multiple solutions.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{remark}[theorem]{Remark}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{example}[theorem]{Example}
\allowdisplaybreaks

\section{Introduction}

We wish to give sufficient conditions for the existence of nontrivial 
periodic solutions to the second-order impulsive Hamiltonian system
\begin{equation}  \label{1}
\begin{gathered}
-\ddot{u}(t)+A(t)u(t)= \lambda \nabla
F(t,u(t))+\nabla H(u(t)), \quad\text{a.e. } t\in [0,T],
\\
\Delta (\dot{u}_i(t_j))=I_{ij}(u_i(t_j)),\quad  i=1,2,\ldots,N,\;
  j=1,2,\ldots,p,\\
u(0)- u(T)= \dot{u}(0)- \dot{u}(T)=0,
\end{gathered}
\end{equation}
where $N\geq 1$, $p \geq 2$, $u=(u_1,\ldots,u_N)$, $T>0$, $\lambda>0$ 
is a parameter,  $A:[0,T]\to \mathbb{R}^{N\times N}$ is a
continuous map from the interval $[0,T]$ to the set of $N\times N$
symmetric matrices, $t_j$, $j=1,2,\ldots,p$, are the
instants at which the impulses occur,
$0=t_0<t_{1}<\ldots<t_{p}<t_{p+1}=T$, and $\Delta(\dot{u}_i(t_j))=
\dot{u}_i(t_j^+)-\dot{u}_i(t_j^-)=\lim_{t\to t_j^+}\dot{u}_i(t)
-\lim_{t\to t_j^-}\dot{u}_i(t)$.
Without further mention, the following conditions are assumed to hold throughout 
the remainder of this article.
The functions $I_{ij}:\mathbb{R}\to\mathbb{R}$ are Lipschitz continuous with the
Lipschitz constants $L_{ij}> 0$, i.e.,
\begin{equation}  \label{new-1}
|I_{ij}(s_1)-I_{ij}(s_2)| \leq L_{ij}|s_1-s_2|
\end{equation}
for every $s_1$, $s_2\in\mathbb{R}$, and $I_{ij}(0)=0$ for
$i=1,2,\ldots,N$, $j=1,2,\ldots,p$. 
In addition, $F:[0,T]\times{\mathbb{R}}^N\to\mathbb{R}$ is measurable
 with respect to $t$ for all $u\in{\mathbb{R}}^N$,
continuously differentiable in $u$ for almost every $t\in[0,T]$,
$F(t,0,\ldots,0)=0$ for $t\in[0,T]$, and satisfies the
standard summability condition
\begin{equation}  \label{2}
\sup_{|\xi|\leq a}\max\{|F(\cdot,\xi)|,\ |\nabla F(\cdot,\xi)|\}\in L^1([0,T])
\end{equation}
for any $a>0$. Also, the function $H:{\mathbb{R}}^N\to\mathbb{R}$ is continuously
differentiable, $\nabla H$ is Lipschitz continuous with the Lipschitz 
constant $L> 0$, i.e.,
\begin{equation}  \label{new400}
|\nabla H(\xi_1)-\nabla H(\xi_2)| \leq L|\xi_1-\xi_2|
\end{equation}
for every $\xi_1$, $\xi_2\in\mathbb{R}^N$,
\begin{equation}  \label{new401}
H(0,\ldots,0)=0, \quad \text{and} \quad \nabla H(0,\ldots,0)=0.
\end{equation}
Assuming that $\nabla F:[0,T]\times{\mathbb{R}}^N\to\mathbb{R}$ is
continuous implies that condition \eqref{2} is satisfied.

The study of multiplicity of solutions of Hamiltonian systems, 
as a special case of dynamical systems,
is important mathematically
as well as being interesting from a practical view point since
these systems form a natural framework for
mathematical models of many natural phenomena in fluid
mechanics, gas dynamics, nuclear physics, relativistic mechanics,
etc. For background, theory, and applications of
Hamiltonian systems, we refer the reader to
\cite{CES,Ra1,T}. Inspired by the monographs \cite{MW,Ra2}, the
existence and multiplicity of periodic solutions for Hamiltonian
systems have been investigated using variational methods by many authors; 
for example, see
\cite{C,CR,D1,F,GA,IJ,New5,New2,WZ,ZL2,ZT,ZoL,ZZ} and the references
therein.

In recent years, investigating the existence of solutions to impulsive boundary 
value problems has become increasingly important due to their role in models 
of such things as spacecraft control, impact mechanics, physics, chemistry, 
chemical engineering, population
dynamics, biotechnology, economics, and inspection process in
operations research. We refer the reader to  \cite{BS,BHN,GHO,LBS,SP} for a 
general discussion
of impulsive differential equations and their applications. 
There have been many approaches used to study
existence of solutions of impulsive differential equations,
such as fixed point theory, topological degree, continuation methods, 
coincidence degree theory,
upper and lower solution methods, and the monotone
iterative method; see, for example, \cite{AO,FN,LJ} and
references contained therein. Recently, in \cite{BD1,NO,TG,ZD}, the authors used
critical point theory to study the existence and multiplicity of solutions
of impulsive problems.

Very recently, a great deal of work has been done on the existence
of multiple solutions to second-order impulsive Hamiltonian
systems. In \cite{CH,SCNO}, based on variational methods and
critical point theory, the existence of multiple solutions to
second-order impulsive Hamiltonian systems was established. We
also refer the interested reader to \cite{LS,SCN,ZhL} in which
second order Hamiltonian systems with impulsive effects have been
examined. In \cite{GHK}, the present authors used variational methods and critical
point theory to study the existence of infinitely many periodic
solutions to a class of perturbed second-order impulsive Hamiltonian systems.
For a discussion of multiple solutions to boundary value problems via variational 
methods and critical point theory, we refer the
reader to \cite{D2,HFK}.

Our results here are motivated by the recent papers \cite{CH,D2,SCNO}. 
We begin by obtaining the existence
of a nontrivial periodic solution by combining algebraic conditions on $F$ and
$H$. Another result, Theorem \ref{t7} below,
is concerned with the existence of three periodic solutions. We obtain it by 
combining the use of algebraic conditions
on the functions $F$ and $H$ to obtain the existence of two distinct solutions 
and then applying the mountain
pass theorem of Pucci and Serrin to obtain the third solution. 
This approach of combining techniques to obtain multiple
solutions of boundary value problems is somewhat unique.

\section{Preliminaries}

For a given non-empty set $X$ and two functionals
$\Phi$, $\Psi:X\to\mathbb{R}$, we define the functions
\begin{gather*}
  \vartheta(r_1,r_2)=\inf_{v\in \Phi^{-1}(r_1,r_2)}
\frac{\sup_{u\in \Phi^{-1}(r_1,r_2)}\Psi(u)-\Psi(v)}
{r_2-\Phi(v)}, 
\\
\rho_{1}(r_1,r_2)=\sup_{v\in \Phi^{-1}(r_1,r_2)}\frac{\Psi(v)-
\sup_{u\in \Phi^{-1}(-\infty,r_1)}\Psi(u)}{\Phi(v)-r_1} 
\end{gather*}
for all $r_1$, $r_2\in\mathbb{R}$ with $r_1<r_2$, and
$$  
\rho_{2}(r)=\sup_{v\in \Phi^{-1}(r,\infty)}\frac{\Psi(v)-
\sup_{u\in \Phi^{-1}(-\infty,r)}\Psi(u)}{\Phi(v)-r}  
$$
for all $r\in\mathbb{R}$. We also define the functional 
$I_{\lambda}: X\to\mathbb{R}$ by
$$  
I_{\lambda} : \Phi - \lambda \Psi.  
$$
The following two results are due to Bonanno; they will be used 
in the proofs of our main results.

\begin{theorem}[{\cite[Theorem 5.1]{B}}]  \label{t1}
 Let $X$ be a real Banach space, $\Phi:X\to\mathbb{R}$ be a
sequentially weakly lower semicontinuous, coercive, and
continuously G\^ateaux differentiable functional whose G\^ateaux
derivative admits a continuous inverse on $X^\ast$, and let
$\Psi:X\to\mathbb{R}$ be a continuously G\^ateaux differentiable
functional whose G\^ateaux derivative is compact. Assume there
are constants $r_1$, $r_2\in\mathbb{R}$, $r_1<r_2$, such that
$$  \vartheta(r_1,r_2)<\rho_{1}(r_1,r_2).  $$
Then, for each
$\lambda\in(\frac{1}{\rho_{1}(r_1,r_2)},\frac{1}{\vartheta(r_1,r_2)})$
there exists $u_{0,\lambda}\in\Phi^{-1}(r_1,r_2)$ such that
$I_{\lambda}(u_{0,\lambda})\leq I_\lambda(u)$ for all $u\in
\Phi^{-1}(r_1,r_2)$ and $I'_\lambda(u_{0,\lambda})=0$.
\end{theorem}

\begin{theorem}[{\cite[Theorem 5.3]{B}}]  \label{t2}
Let $X$ be a real Banach space, $\Phi:X\to\mathbb{R}$ be a
continuously G\^ateaux differentiable
functional whose G\^ateaux derivative admits a continuous inverse on
$X^\ast$, and let $\Psi:X\to\mathbb{R}$ be a continuously G\^ateaux
differentiable functional whose G\^ateaux derivative is compact. Choose
$r$ so that $\inf_{X}\Phi<r<\sup_X\Phi$,
$\rho_{2}(r)>0$, and for each $\lambda>\frac{1}{\rho_{2}(r)}$, the
functional $I_{\lambda}:=\Phi-\lambda \Psi$ is coercive. Then, for
each $\lambda\in(\frac{1}{\rho_{2}(r)},+\infty)$, there exists
$u_{0,\lambda}\in\Phi^{-1}(r,+\infty)$ such that
$I_{\lambda}(u_{0,\lambda})\leq I_\lambda(u)$ for all  $u\in
\Phi^{-1}(r,+\infty)$ and $I'_\lambda(u_{0,\lambda})=0$.
\end{theorem}

We assume throughout that the matrix $A$ satisfies the following conditions:
\begin{itemize}
\item[(M1)] $A(t)=(a_{kl}(t))$, $k=1,\ldots,N$, $l=1,\ldots,N$, is a symmetric matrix with $a_{kl}\in
L^\infty{[0,T]}$ for any $t\in[0,T]$;

\item[(M2)] There exists $\delta>0$ such that
$(A(t)\xi,\xi)\geq\delta|\xi|^2$ for any $\xi\in\mathbb{R}^N$ and
a.e. $t\in[0,T]$, where $(\cdot,\cdot)$ denotes the inner product
in $\mathbb{R}^N$.
\end{itemize}

Next, we recall some basic concepts that will be used in what follows. Set
\begin{align*}   
E= \Big\{&u : [0,T]\to {\mathbb{R}}^N: u
\text{ is absolutely continuous}, \\
& u(0)= u(T),\; \dot{ u} \in L^2([0,T],{\mathbb{R}}^N)\Big\} 
\end{align*}
with the inner product
$$  
\prec u,v\succ_E=\int_0^{T} [(\dot{u}(t),\dot{v}(t))+
(u(t),v(t))]dt. 
 $$
The corresponding norm is defined by
$$ 
\|u\|_E=\int_0^{T} (|\dot{u}(t)|^2+ |u(t)|^2)dt \quad \text{for all }  u\in E.
  $$
For every $u$, $v\in E$, we define
$$ 
\prec u,v\succ=\int_0^{T}[(\dot{u}(t),\dot{v}(t))+ (A(t)u(t),v(t))]dt,
$$
and we observe that conditions (M1) and (M2) ensure that this defines an inner
 product in $E$.
Then $E$ is a separable and reflexive Banach space with the norm
$$  
\|u\|= \prec u,u\succ^{1/2} \quad \text{for all }  u\in E. 
$$
Clearly, $E$ is an uniformly convex Banach space.

A simple computation shows that
\[
(A(t)\xi,\xi)=\sum_{k,l=1}^Na_{kl}(t)\xi_k\xi_l
\leq\sum_{k,l=1}^N\|a_{kl}\|_{L^{\infty}}|\xi|^2
\]
for every $t\in[0,T]$ and $\xi\in \mathbb{R}^N$. Along
 with condition (A2), this implies
\begin{equation}  \label{3}
\sqrt{m}\|u\|_E\leq\|u\|\leq\sqrt{M}\|u\|_E,
\end{equation}
where
$m=\min\{1,\delta\}$ and $M=\max\{1,\sum_{k,l=1}^{N}\|a_{kl}\|_\infty\}$, 
which means the norm $\|\cdot\|$ is equivalent to the norm $\|\cdot\|_E$.
Since $(E,\|\cdot\|)$ is compactly embedded in $C([0,T],\mathbb{R}^N)$ 
(see \cite{MW}), there
exists a positive constant $c$ such that
\begin{equation}  \label{4}
\|u\|_\infty \leq c\| u\|,
\end{equation}
where $\|u\|_\infty =\max_{t\in[0,T]}\mid u(t)\mid$ and
$c=\sqrt{\frac{2}{m}}\max\{\frac{1}{\sqrt{T}},\sqrt{T}\}$ (see
\cite{CH}).

For $u\in E$, $\Delta \dot{u}(t)=\dot{u}(t^{+})-\dot{u}(t^-)=0$ does not
necessarily hold for every $t\in(0,T)$, and the derivative $\dot{u}$
may possess discontinuities. This leads to the impulsive effects.

Next, we define what we mean by a solution of \eqref{1}.

\begin{definition} \rm
A function $u\in\{u\in E: \dot{u}\in (W^{1,2}(t_j,t_{j+1}))^N, \;
j=0,1,2,\ldots,p\}$ is said to be a classical solution of the
problem \eqref{1} if $u$ satisfies the differential equation, the impulse relations,
and the boundary conditions given in problem \eqref{1}.
\end{definition}

\begin{definition} \rm
By a weak solution of the problem \eqref{1}, we mean any $u\in E$ such that
\begin{equation}  \label{N1}
\begin{aligned}
&\int_0^{T} \Big[(\dot{u}(t),\dot{v}(t))+(A(t)u(t),v(t))
 -(\nabla H(u(t)),v(t))\Big]dt \\
&+\sum_{j=1}^p\sum_{i=1}^{N}I_{ij}(u_i(t_j))v_i(t_j)
 -\lambda\int_0^T(\nabla F(t,u(t)),v(t))dt= 0
\end{aligned}
\end{equation}
for every $v\in E$.
\end{definition}

An important relationship between a weak solution and a classical 
solution of \eqref{1} is given in the next lemma.

\begin{lemma}[{\cite[Lemma 2.2]{GHK}}]  \label{new lemma}
If $u\in E$ is a weak solution of \eqref{1}, then $u$ is a
classical solution of \eqref{1}.
\end{lemma}

In what follows, unless stated otherwise, by a solution of \eqref{1} we
will always mean a classical solution.
Without further mention, we will assume throughout that
$$  
K:=c^2(2LT+\sum_{j=1}^{p}\sum_{i=1}^{N}L_{ij})<1.  
$$
The following proposition is needed in the proofs of our main
results.

\begin{proposition}   \label{p1}
Let $J:E\to E^{*}$ be the operator defined by
\begin{align*}
J(u)v&=\int_0^{T} \big[(\dot{u}(t),\dot{v}(t))+(A(t)u(t),v(t))
 -(\nabla H(u(t)),v(t))\big]dt \\
&\quad +\sum_{j=1}^p\sum_{i=1}^{N}I_{ij}(u_i(t_j))v_i(t_j)
\end{align*}
for every $u$, $v\in E$. Then $J$ admits a continuous inverse on
$E^{*}$.
\end{proposition}

\begin{proof} 
Since $-L|\xi|^{2}\leq(\nabla H(\xi),\xi) \leq L|\xi|^{2}$ for every
$\xi\in\mathbb{R}^N$, and $ -L_{ij}|s|^{2}\leq I_{ij}(s)s \leq
L_{ij}|s|^{2} $ for every $s\in\mathbb{R}$ and all
$i=1,2,\ldots,N$, $j=1,2,\ldots,p$, in view of \eqref{4}, we have
\begin{align*}
J(u)u &= \int_0^{T} \Big[(\dot{u}(t),\dot{u}(t))
 +(A(t)u(t),u(t))-(\nabla H(u(t)),u(t))\Big]dt\\
&\quad +\sum_{j=1}^p\sum_{i=1}^{N}I_{ij}(u_i(t_j))u_i(t_j)\\
&\geq \Big(1-c^2LT-c^2\sum_{j=1}^{p}\sum_{i=1}^{N}L_{ij}\Big)\|u\|^2\\
&> \Big(1-K\Big)\|u\|^2.
\end{align*}
Since $K<1$, $J$ is coercive. Now for any $u$, $v\in E$,
\begin{align*}
\langle J(u)-J(v), u- v\rangle
&= \int_0^{T}(\dot{u}(t)-\dot{v}(t),\dot{u}(t)- \dot{v}(t))dt\\
&\quad + \sum_{j=1}^p\sum_{i=1}^{N}(I_{ij}(u_i(t_j))-I_{ij}(v_i(t_j)))(u_i(t_j)-v_i(t_j))\\
&\quad -\int_0^{T}
(\nabla H(u(t))-\nabla H(v(t)),u(t)- v(t))dt\\
&\geq \Big(1-c^2LT-c^2\sum_{j=1}^{p}
\sum_{i=1}^{N}L_{ij}\Big)\|u- v\|^{2}\\
&>(1-K)\|u- v\|^{2},
\end{align*}
so $J$ is uniformly
monotone. By \cite[Theorem 26.A (d)]{Z}, $J^{-1}$ exists and is
continuous on $E^*$.
\end{proof}

\section{Main results}

Our first existence result is contained in the following theorem.
For a given function $w\in E$ and a given nonnegative constant $r$
with
$$  
r\neq\frac{1}{2}(1+K)\|w\|^2,  
$$
we set
$$  
a_{w}(r):=\frac{ \int_0^{T}
\max_{|\xi|\leq c(\frac{2r}{1-K})^{1/2}}F(t,\xi)dt -
\int_0^{T}F(t,w(t))dt}{r-\frac{1}{2}(1+K)\|w\|^2}.
 $$

\begin{theorem}  \label{t3}
Assume that there exist constants $r_1 \geq 0$ and
$r_2 >0$, and a function $w\in E$ such that
\begin{itemize}
\item[(A1)] \ $\big(\frac{2r_1}{1-K}\big)^{1/2}
<\|w\|<\big(\frac{2r_2}{1+K}\big)^{1/2}$,

\item[(A2)]  $a_w(r_2)<a_w(r_1)$.
\end{itemize}
Then, for each
$\lambda\in\big(\frac{1}{a_w(r_1)},\ \frac{1}{a_w(r_2)} \big)$,
problem \eqref{1} has a non-trivial periodic solution $u^{\ast} \in E$ such that
$$  
r_{1}<\frac{1}{2}\|u^{\ast}\|^{2} +\sum_{j=1}^{p}
\sum_{i=1}^N\int_0^{u^{\ast}_i(t_j)}I_{ij}(s)ds
-\int_0^{T}H(u^{\ast}(t))dt <r_{2}.
$$
\end{theorem}

\begin{remark} \rm
In the above theorem, and in the results below, by $u^{*}$ we mean 
the vector $(u^{*}_1, u^{*}_2, \dots, u^{*}_N)$.
\end{remark}

\begin{proof}
Choose $\lambda$ as in the conclusion of the theorem. To apply Theorem
\ref{t1} to our problem, we take $X=E$ and define the
functionals $\Phi$, $\Psi$, $I_{\lambda}:X\to \mathbb{R}$ by
\begin{gather*} 
\Phi(u)= \frac{1}{2}\|u\|^{2}
+\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{u_i(t_j)}I_{ij}(s)ds-\int_0^{T}H(u(t))dt, \\
\Psi(u)= \int_0^{T}F(t,u(t))dt, \\
 I_{\lambda}(u) = \Phi(u) - \lambda \Psi(u)
\end{gather*}
for every $u\in X$. It is well known that $\Psi$ is a G\^{a}teaux differentiable
functional whose G\^{a}teaux derivative at the point $u\in X$ is
the functional $\Psi'(u)\in X^{*}$ given by
\begin{equation}  \label{new500}
 \Psi'(u)v= \int_0^T(\nabla F(t,u(t)),v(t))dt
\end{equation}
for every $v\in X$, and that
$\Psi':X \to X^{*}$ is a compact operator. Moreover, $\Phi$ is a
G\^{a}teaux differentiable functional whose G\^{a}teaux derivative
at the point $u\in X$ is the functional $\Phi'(u)\in X^{*}$ given by
\begin{equation}  \label{new501}
\begin{aligned}
\Phi'(u)v
&=\int_0^{T} \Big[(\dot{u}(t),\dot{v}(t))+(A(t)u(t),v(t))
-(\nabla H(u(t)),v(t))\Big]dt \\
&\quad +\sum_{j=1}^p\sum_{i=1}^{N}I_{ij}(u_i(t_j))v_i(t_j)
\end{aligned}
\end{equation}
for every $v = (v_1,v_2, \dots, v_N) \in X$. 
Furthermore, $\Phi$ is sequentially weakly lower
semicontinuous (see \cite{GHK}). From \eqref{new400} and \eqref{new401}, 
we have $|H(\xi)|\leq L|\xi|^2$ for all $\xi\in\mathbb{R}^N$.
From \eqref{new-1}, \eqref{4}, and the fact that $I_{ij}(0) = 0$, we have
\begin{equation}\label{5}
\frac{1}{2}(1-K)\|u\|^2
\leq\Phi(u)\leq\frac{1}{2}(1+K)\|u\|^2
\end{equation}
for $u\in X$.
Condition (A1) together with \eqref{5} implies
$$   
r_1<\Phi(w)<r_2.  
$$
From \eqref{4} and \eqref{5}, for each $u\in X$,
\begin{align*}
\Phi^{-1}(-\infty,r_{2})
&=\{u\in X :  \Phi(u)< r_{2}\} \\
&\subseteq \big\{u\in X :  \frac{1}{2}(1-K)\|u\|^2< r_{2}\big\}\\
&\subseteq \big\{ u\in X : |u(t)|\leq c\big(\frac{2r_{2}}{1-K}\big)^{1/2}
\text{ for each }  t\in [0,T]\big\},
\end{align*}
and it follows that
$$  
\sup_{u\in\Phi^{-1}(-\infty,r_{2})}\Psi(u)
=\sup_{u\in\Phi^{-1}(-\infty,r_{2})}\int_0^{T}F(t,u(t))dt
\leq \int_0^{T} \max_{|\xi|\leq
c(\frac{2r_{2}}{1-K})^{1/2}}F(t,\xi)dt. 
 $$
Therefore,
\begin{align*}
\vartheta(r_1,r_2)
&\leq \frac{\sup_{u\in \Phi^{-1}(-\infty,r_2)}\Psi(u)-\Psi(w)}{r_2-\Phi(w)}\\
&\leq \frac{\int_0^{T} \max_{|\xi|\leq
c\big(\frac{2r_{2}}{1-K}\big)^{1/2}}F(t,\xi)dt-\Psi(w)}{r_2-\Phi(w)}\\
&\leq \frac{ \int_0^{T} \max_{|\xi|\leq
c\big(\frac{2r_{2}}{1-K}\big)^{1/2}}F(t,\xi)dt
-\int_0^{T}F(t,w(t))dt}{r_{2}-\frac{1}{2}(1+K)\|w\|^2}\\
&= a_w(r_2).
\end{align*}
On the other hand, arguing as before,
\begin{align*}
\rho(r_1,r_2)
&\geq \frac{\Psi(w)-\sup_{u\in \Phi^{-1}(-\infty,r_1)}\Psi(u)}{\Phi(w)-r_1}\\
&\geq \frac{\Psi(w)-\int_0^{T} \max_{|\xi|\leq
c\big(\frac{2r_{1}}{1-K}\big)^{1/2}}F(t,\xi)dt}{\Phi(w)-r_1}\\
&\geq \frac{ \int_0^{T}F(t,w(t))dt-
\int_0^{T} \max_{|\xi|\leq
c\big(\frac{2r_{1}}{1-K}\big)^{1/2}}F(t,\xi)dt}
{\frac{1}{2}(1+K)\|w\|^2-r_{1}}\\
&= a_w(r_1).
\end{align*}
Hence, from condition (A1), we have
$\vartheta(r_1,r_2)<\rho(r_1,r_2).$ Therefore, by Theorem \ref{t1}, 
for each $\lambda\in\big(\frac{1}{a_w(r_1)},\frac{1}{a_w(r_2)} \big)$, 
the functional $I_{\lambda}$ admits at least one critical point $u^{\ast}\in X$ 
such that $r_1<\Phi(u^{\ast})<r_2$, that is, $u^{\ast}$ is a
nontrivial local minimum for $I_{\lambda}$ in $X$.

Since weak solutions of problem \eqref{1} are precisely the
solutions of the equation $I_{\lambda}'(u)=0$ (see \eqref{N1}, \eqref{new500}, 
and \eqref{new501}), $u^{\ast}$ is a weak solution of problem \eqref{1}.
 In view of Lemma \ref{new lemma}, this completes the proof of the theorem.
\end{proof}

The following corollary provides a sufficient condition for applying
Theorem \ref{t3} that does not require knowledge of two constants
$r_1$, $r_2$ and a test function $w$ satisfying (A1) and (a$_2$).

Let
$$  
D=\frac{(T-t_p)^2}{t_1t_p^{2}}+\frac{t_1}{3t_p^2}(t_p^2+t_pT+T^2)+(t_p-t_1)
+\frac{T-t_p}{t_p^2}+\frac{1}{3t_p^2}(T^3-t_p^3)>0,  
$$
and for a given nonnegative constant $\theta$ and a positive constant
$\eta$, with
$$  
(1-K)\theta^{2} \neq c^{2}(1+K)DM\eta^2,  
$$
let
$$
  b_{\eta}(\theta):=\frac{\int_0^{T}
\max_{|\xi|\leq \theta}F(t,\xi)dt -
\int_{t_{1}}^{t_{p}}F(t,\eta\mathcal{\varepsilon})dt}{\frac{1}{2}(1-K)\theta^{2}
-\frac{1}{2}c^{2}(1+K)DM\eta^2}  
$$
where $\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$.

\begin{corollary}  \label{c1}
Assume that there exist constants $\theta_1 \geq 0$, $\theta_2 >0$, and 
$\eta >0$ with
$\frac{\theta_{1}}{c\sqrt{Dm}}< \eta 
< \frac{\theta_{2}}{c}\sqrt{\frac{1-K}{DM(1+K)}}$ such
that
\begin{itemize}
\item[(A3)]  $F(t,\xi)\geq 0$ for each $t\in[0,t_1]\cup[t_p,T]$
and $|\xi|\leq\frac{\eta T}{t_p}$,

\item[(A4)]  $b_\eta(\theta_2)<b_\eta(\theta_1)$.
\end{itemize}
Then, for each $\lambda\in\big(\frac{1}{c^{2}}\frac{1}{b_\eta(\theta_1)},
\frac{1}{c^{2}}\frac{1}{b_\eta(\theta_2)} \big)$, the problem
\eqref{1} has a non-trivial periodic
solution $u^{\ast} \in E$ such that
\begin{align*}   
\frac{1}{2}(1-K)\big(\frac{\theta_{1}}{c}\big)^{2}
&<\frac{1}{2}\|u^{\ast}\|^{2}
+\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{u^{\ast}_i(t_j)}I_{ij}(s)ds \\
&\quad -\int_0^{T}H(u^{\ast}(t))dt
<\frac{1}{2}(1-K)\left(\frac{\theta_{2}}{c}\right)^{2}.  
\end{align*}
\end{corollary}

\begin{proof} 
Choose $r_1=\frac{1}{2}(1-K)(\frac{\theta_{1}}{c})^{2}$,
$r_{2}=\frac{1}{2}(1-K)(\frac{\theta_{2}}{c})^{2}$,
and
\begin{equation}  \label{8}
w(t)= \begin{cases}
(T+\frac{t_p-T}{t_1}t)\frac{\eta\mathcal{\varepsilon}}{t_p},& t\in[0,t_1),\\
\eta\mathcal{\varepsilon} , & t\in [t_1,t_p],\\
\frac{\eta\mathcal{\varepsilon}}{t_p}t, &t\in(t_p,T].
\end{cases}
\end{equation}
Then $w\in E$ and $\|w\|_E^2=D\eta^2$. By \eqref{3},
\begin{equation}  \label{9}
Dm\eta^2\leq\|w\|^2\leq DM\eta^2,
\end{equation}
and this together with the condition on $\eta$ implies
(A1) is satisfied. Moreover, since 
$0\leq w(t)\leq \frac{\eta T}{t_p}$ for each $t\in [0,T]$ and (A3) holds, we have
\begin{equation}  \label{10}
\int_0^{t_{1}} F\Big(t,\Big(T+\frac{t_1-T}{t_p}t\Big)
 \frac{\eta\mathcal{\varepsilon}}{t_p}\Big)dt
+ \int_{t_{p}}^{T}F\Big(t,\frac{\eta\mathcal{\varepsilon}}{t_p}t\Big)dt\geq 0.
\end{equation}
Therefore, from \eqref{9} and \eqref{10} it follows that
\begin{gather*}
a_{w}(r_2)= \frac{ \int_0^{T} \max_{|\xi|\leq
c(\frac{2r_{2}}{1-K})^{1/2}}F(t,\xi)dt
-\int_0^{T}F(t,w(t))dt}{r_{2}-\frac{1}{2}(1+K)\|w\|^2}\leq
c^2b_\eta(\theta_2), 
\\
c^2b_\eta(\theta_1)\leq \frac{
\int_0^{T}F(t,w(t))dt- \int_0^{T} \max_{|\xi|\leq
c(\frac{2r_{1}}{1-K})^{1/2}}F(t,\xi)dt}
{\frac{1}{2}(1+K)\|w\|^2-r_{1}}=a_{w}(r_1).
\end{gather*}
Therefore, (A4) implies that (A2) is satisfied.
Hence, by Theorem \ref{t3},the conclusion of the corollary follows.
\end{proof}

An easy consequence of Corollary \ref{c1} is the following existence result.

\begin{corollary}  \label{c2}
In addition to {\rm (A3)}, assume there exist $\theta >0$ and $\eta >0$
with $\eta< \frac{\theta}{c}\sqrt{\frac{1-K}{DM(1+K)}}$ such
that
\begin{itemize}
\item[(A5)] 
\[
\frac{ \int_0^{T} \max_{|\xi| \leq \theta}F(t,\xi)dt}{\theta^{2}} < {\frac{1-K}
{c^{2}(1+K)DM}} \frac{
\int_{t_{1}}^{t_{p}}F(t,\eta\mathcal{\varepsilon})dt}{ \eta^2},
\]
where
$\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$.
\end{itemize}
Then, for each 
\[
\lambda\in\Big(\frac{ (1+K)DM\eta^2}{
2\int_{t_{1}}^{t_{p}}F(t,\eta\mathcal{\varepsilon})dt},\
\frac{(1-K)\theta^{2}}{ 2c^{2}\int_0^{T}
\max_{|\xi|\leq \theta}F(t,\xi)dt} \Big),
\]
problem \eqref{1} has a non-trivial periodic solution
$u^{\ast} \in E$ such that
$$   
0<\frac{1}{2}\|u^{\ast}\|^{2}
+\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{u^{\ast}_i(t_j)}I_{ij}(s)ds
-\int_0^{T}H(u^{\ast}(t))dt <\frac{1}{2}(1-K)\big(\frac{\theta}{c}\big)^{2}.
$$
\end{corollary}

\begin{proof}
Choosing $\theta_1=0$ and $\theta_2=\theta$, we have
\begin{align*}
b_{\eta}(\theta)
&< \frac{\big(1-\frac{c^{2}(1+K)DM\eta^2}
{(1-K)\theta^{2}}\big)  \int_0^{T}
\max_{|\xi|\leq \theta}F(t,\xi)dt}{\frac{1}{2}(1-K)\theta^{2}
-\frac{1}{2}c^{2}(1+K)DM\eta^2}\\
&=\frac{ \int_0^{T} \max_{|\xi|\leq
\theta}F(t,\xi)dt}{\frac{1}{2}(1-K)\theta^{2}}\\
&< \frac{1} {\frac{1}{2}c^{2}(1+K)DM}\frac{
\int_{t_{1}}^{t_{p}}F(t,\eta\mathcal{\varepsilon})dt}{\eta^2}
=b_{\eta}(0).
\end{align*}
In particular,
$$  
b_{\eta}(\theta)<\frac{ \int_0^{T} \max_{|\xi|\leq
\theta}F(t,\xi)dt}{\frac{1}{2}(1-K)\theta^{2}}.  
$$
The conclusion then follows from Corollary \ref{c1}.
\end{proof}

Next, we present an application of Theorem \ref{t2} that will be
used to obtain multiple solutions to problem \eqref{1}.

\begin{theorem}  \label{t4}
Assume there exist a constant $\bar{r} >0$ and a function $\bar{w}$ with
$\frac{2\bar{r}}{1-K} <\|\bar{w}\|^2$ such that
\begin{itemize}
\item[(B1)]   $ \int_0^{T} \max_{|\xi|\leq
c(\frac{2\bar{r}}{1-K})^{1/2}}F(t,\xi)dt<
\int_0^{T}F(t,\bar{w}(t))dt$;

\item[(B2)]   $\limsup_{|\xi|\to +\infty}\frac{F(t,\xi)}{|\xi|^{2}} \leq 0$ 
uniformly for $t\in [0,T]$.
\end{itemize}
Then, for each
$\lambda\in\left(\bar{\lambda},+\infty\right)$, where
\[
\bar{\lambda}:=\frac{ \frac{1}{2}(1+K)\|\bar{w}\|^2-\bar{r}}{
\int_0^{T}F(t,\bar{w}(t))dt- \int_0^{T}
\max_{|\xi|\leq c(\frac{2\bar{r}}{1-K})^{1/2}}F(t,\xi)dt},
\]
problem \eqref{1} has at least one non-trivial
periodic solution $\bar{u}\in E$ such that
$$  
\frac{1}{2}\|\bar{u}\|^{2} +\sum_{j=1}^{p}
\sum_{i=1}^N\int_0^{\bar{u}_i(t_j)}I_{ij}(s)ds
-\int_0^{T}H(\bar{u}(t))dt >\bar{r}.
$$
\end{theorem}

\begin{proof} 
Choose $\lambda$ as in the conclusion of the theorem. Taking $X$ and the
functionals $\Phi$ and $\Psi$ as in the proof of Theorem \ref{t3},
we see that all the regularity assumptions required in Theorem
\ref{t2} are satisfied. By (B2), there is a
constant $\epsilon$ and a function $h_{\epsilon}\in L^1([0,T])$
with $0<\epsilon<\frac{1-K}{2\lambda c^2}$ such that
\begin{equation}  \label{11}
F(t,\xi)\leq \epsilon |\xi|^{2}+h_{\epsilon}(t)\quad \text{for all }
t\in [0,T],\; \xi\in \mathbb{R}^N.
\end{equation}
From the definitions of $\Phi$ and $\Psi$, \eqref{4}, \eqref{5} and \eqref{11},
 we obtain
\begin{align*}
I_{\lambda}(u) 
&\geq \frac{1}{2}(1-K)\|u\|^2-\lambda\epsilon
\int_0^{T}|u(t)|^2dt-\lambda \int_0^{T}h_{\epsilon}(t)dt\\
&\geq \frac{1}{2}\Big(1-K-\lambda\epsilon c^{2}\Big)\|u\|^2-\lambda
\|h_{\epsilon}\|_{L^1([0,T])}.
\end{align*}
Since $1-K-\lambda\epsilon c^{2}>0$, the functional $I_{\lambda}$ is coercive.
Arguing as in the proof of Theorem \ref{t3} shows that
$$   
\rho_{2}(\bar{r})\geq\frac{
\int_0^{T}F(t,\bar{w}(t))dt- \int_0^{T}
\max_{|\xi|\leq c(\frac{2\bar{r}}{1-K})^{1/2}}F(t,\xi)dt}
{\frac{1}{2}(1+K)\|\bar{w}\|^2-\bar{r}} > 0  
$$
by (B1) and (B2). By Theorem \ref{t2}, the functional
$I_{\lambda}$ admits at least one local minimum $\bar{u}\in
X$ such that $\frac{1}{2}\|\bar{u}\|^{2}
+\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{\bar{u}_i(t_j)}I_{ij}(s)ds-\int_0^{T}H(\bar{u}(t))dt
>\bar{r}$, and the conclusion follows.
\end{proof}

The following corollary provides a sufficient condition for
applying Theorem \ref{t4} that does not require knowledge of a
constant $\overline{r}$ and a test function $\overline{w}$
satisfying (B1) and (B2).

\begin{corollary}  \label{c3}
Assume that {\rm (A3)}, and {\rm (B2)} hold and there exist positive constants
$\bar{\theta}$ and $\bar{\eta}$ with 
$\frac{\bar{\theta}}{c\sqrt{Dm}}<\bar{\eta}$ such that
\begin{itemize}
\item[(B3)]  
\[ \int_0^{T} \max_{|\xi|\leq
\bar{\theta}}F(t,\xi)dt <
\int_{t_{1}}^{t_{p}}F(t,\bar{\eta}\mathcal{\varepsilon})dt
\]
where $\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$.
\end{itemize}
Then, for each $\lambda\in(\bar{\lambda'},+\infty)$,
where 
\[
\bar{\lambda'}:=\frac{ \frac{1}{2}(1+K)DM\bar{\eta}^2-
\frac{1}{2}(1-K) \big(\frac{\bar{\theta}}{c}\big)^{2}}
{\int_{t_{1}}^{t_{p}}F(t,\bar{\eta}\mathcal{\varepsilon})dt-
\int_0^{T} \max_{|\xi|\leq \bar{\theta}}F(t,\xi)dt},
\]
problem \eqref{1} has at least one non-trivial periodic solution $\bar{u}\in E$ 
such that
$$
\frac{1}{2}\|\bar{u}\|^{2}
+\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{\bar{u}_i(t_j)}I_{ij}(s)ds
-\int_0^{T}H(\bar{u}(t))dt
>\frac{1}{2}(1-K)\big(\frac{\bar{\theta}}{c}\big)^{2}.
$$
\end{corollary}

\begin{proof} 
Choose $\bar{r}=\frac{1}{2}(1-K)(\frac{\bar{\theta}}{c})^{2}$ and let
$\bar{w}$ be as in \eqref{8} with $\eta$ replaced by $\bar{\eta}$.
The conclusion follows from an application of Theorem \ref{t4}.
\end{proof}

Next, we point out some results for which the function $F$ is in factored form.
To be precise, consider the problem
\begin{equation}  \label{12}
\begin{gathered} 
-\ddot{u}(t)+A(t)u(t)= \lambda b(t)\nabla
G(u(t))+\nabla H(u(t)), \quad  a.e. t\in [0,T],
\\
\Delta (\dot{u}_i(t_j))=I_{ij}(u_i(t_j)),\quad  i=1,2,\ldots,N,\; j=1,2,\ldots,p,\\
u(0)- u(T)= \dot{u}(0)- \dot{u}(T)=0,
\end{gathered}
\end{equation}
where $b\in L^1([0,T])$, $b(t)\geq 0$ a.e. $t\in[0,T]$,
$b\not\equiv 0$, $G\in C^1(\mathbb{R}^{N},\mathbb{R})$,
$G(0,\ldots,0)=0$, and each component of the vector 
$\nabla G:\mathbb{R}^{N}\to\mathbb{R}^N$ is a
non-negative continuous function.

\begin{remark}  \label{r1} \rm
Since the first term on the right hand side of the equation in \eqref{12} 
is nonnegative, any weak solution of \eqref{12} is nonnegative.
To see this, assume that the set
$\mathcal{A}=\big\{t \in [0,T] : u_0(t)<0\big\}$ is non-empty and of positive
 measure. Let $\bar v(t)=\min\{0, u_0(t)\}$
for all $t \in [0,T]$. Clearly, $\bar v \in E$. Using the fact that $u_0$ 
is a weak solution of \eqref{12}, we have
\begin{align*}
&\int_0^{T} \Big((\dot{u}_0(t),\dot{\bar{v}}(t))
 +(A(t)u_0(t),\bar{v}(t))-(\nabla H(u_0(t)),\bar{v}(t))\Big)dt \\
&+\sum_{j=1}^p\sum_{i=1}^{N}I_{ij}(u_{0i}(t_j))\bar{v}_i(t_j)\\
&=\lambda\int_0^T(b(t)\nabla G(u_0(t)),\bar{v}(t))dt.
\end{align*}
Thus,
\begin{align*}
0& \leq \Big(1-c^2LT-c^2\sum_{j=1}^{p}\sum_{i=1}^{N}L_{ij}\Big)\int_{\mathcal{A}}[
(\dot{u}_0(t),\dot{u}_0(t))+ (A(t)u_0(t),u_0(t))]dt \\
&\leq \int_\mathcal{A} \Big((\dot{u}_0(t),\dot{u}_0(t))+(A(t)u_0(t),u_0(t))-(\nabla
H(u_0(t)),u_0(t))\Big)dt\\
&\quad +\sum_{j=1}^p\sum_{i=1}^{N}I_{ij}(u_{0i}(t_j))u_{0i}(t_j) \leq 0.
\end{align*}
Now $c^2LT+c^2\sum_{j=1}^{p}\sum_{i=1}^{N}L_{ij}<1$, so $u_0=0$ 
is in $\mathcal{A}$, which is a contradiction.
\end{remark}

We will now present some existence results that are consequences of 
Corollaries \ref{c1}, \ref{c2}, and \ref{c3}, respectively.
For a given nonnegative constant $\theta$ and a positive constant $\eta$, with
$$  
(1-K)\theta^{2} \neq c^{2}(1+K)DM\eta^2,  $$
define
$$   
c_{\eta}(\theta):=\frac{ \max_{|\xi|\leq \theta}G(\xi) \int_0^{T}
b(t)dt -G(\eta\mathcal{\varepsilon})
\int_{t_{1}}^{t_{p}}b(t)dt}{\frac{1}{2}(1-K)\theta^{2}
-\frac{1}{2}c^{2}(1+K)DM\eta^2},   
$$
where $\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$.

\begin{corollary}  \label{c4}
Assume that there exist constants $\theta_1 \geq 0$,
$\theta_2 >0$, and $\eta >0$ with
$\frac{\theta_{1}}{c\sqrt{Dm}}<\eta< \frac{\theta_{2}}{c}\sqrt{\frac{1-K}{DM(1+K)}}$ 
such that
\begin{itemize}
\item[(A6)]  $c_\eta(\theta_2)<c_\eta(\theta_1)$.
\end{itemize}
Then, for each
$\lambda\in\big(\frac{1}{c^{2}}\frac{1}{c_\eta(\theta_1)},\
\frac{1}{c^{2}}\frac{1}{c_\eta(\theta_2)} \big)$, problem
\eqref{12} has a positive periodic
solution $u^{\ast} \in E$ such that
\begin{align*} 
\frac{1}{2}(1-K)\big(\frac{\theta_{1}}{c}\big)^{2}
&<\frac{1}{2}\|u^{\ast}\|^{2}
 +\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{u^{\ast}_i(t_j)}I_{ij}(s)ds
 -\int_0^{T}H(u^{\ast}(t))dt \\
&<\frac{1}{2}(1-K)\big(\frac{\theta_{2}}{c}\big)^{2}.  
\end{align*}
\end{corollary}

\begin{corollary}   \label{c5}
Assume that there exist constants $\theta >0$ and $\eta >0$
with $\eta< \frac{\theta}{c}\sqrt{\frac{1-K}{DM(1+K)}}$ such that
\begin{itemize}
\item[(A7)]  
\[
\frac{ \max_{|\xi|\leq \theta}G(\xi) \int_0^{T}
 b(t)dt }{\theta^{2}} < {\frac{1-K}{c^{2}(1+K)DM}} 
\frac{ G(\eta\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt}{ \eta^2},
\]
where  $\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$.
\end{itemize}
Then, for each 
\[
\lambda\in\Big(\frac{ (1+K)DM\eta^2}{
2G(\eta\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt},\
\frac{(1-K)\theta^{2}}{ 2c^{2} \max_{|\xi|\leq
\theta}G(\xi)\int_0^{T}b(t)dt} \Big),
\]
 problem \eqref{12} has a positive periodic solution
$u^{\ast} \in E$ such that
$$  
\frac{1}{2}\|u^{\ast}\|^{2} +\sum_{j=1}^{p}
\sum_{i=1}^N\int_0^{u^{\ast}_i(t_j)}I_{ij}(s)ds
-\int_0^{T}H(u^{\ast}(t))dt
<\frac{1}{2}(1-K)\big(\frac{\theta}{c}\big)^{2}.  
$$
\end{corollary}

\begin{corollary}  \label{c6}
Assume  there exist constants $\bar{\theta} >0$ and $\bar{\eta} >0$ with
$\frac{\bar{\theta}}{c\sqrt{Dm}}<\bar{\eta}$ such that
\begin{itemize}
\item[(B4)]
\[\max_{|\xi|\leq \bar{\theta}}G(\xi)\int_0^{T}b(t)dt<
G(\bar{\eta}\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt,
\]
where $\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$;

\item[(B5)]  $\limsup_{|\xi|\to +\infty
}\frac{G(\xi)}{|\xi|^{2}}\leq 0$.
\end{itemize}
Then, for each
$\lambda\in(\hat{\lambda},+\infty)$, where
\begin{equation} \label{new-201}
 \hat{\lambda} :=\frac{\frac{1}{2}(1+K)DM\bar{\eta}^2-\frac{1}{2}(1-K)
(\frac{\bar{\theta}}{c})^{2}}{
G(\bar{\eta}\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt-
\max_{|\xi|\leq \bar{\theta}}G(\xi)\int_0^{T}b(t)dt},
\end{equation}
problem \eqref{11} has at least one
positive periodic solution $\bar{u}\in E$ such that
$$ 
 \frac{1}{2}\|\bar{u}\|^{2} +\sum_{j=1}^{p}\sum_{i=1}^N
\int_0^{\bar{u}_i(t_j)}I_{ij}(s)ds-\int_0^{T}H(\bar{u}(t))dt
>\frac{1}{2}(1-K)\big(\frac{\bar{\theta}}{c}\big)^{2}. 
 $$
\end{corollary}

One consequence of Corollary \ref{c5} is the following existence
result.

\begin{theorem} \label{t5}
Assume that
\begin{equation}  \label{13}
\lim_{x\to 0^{+}}\frac{ \max_{|\xi|\leq x}G(\xi)}{|x|^{2}}=+\infty.
\end{equation}
Then, for each $\lambda\in(0, \lambda^*)$, where
\[
\lambda^*:=\frac{ 1-K}{ 2c^{2}\int_0^{T}
b(t)dt}\sup_{\theta>0}\frac{\theta^{2}}{ \max_{|\xi|\leq
\theta}G(\xi)},
\]
 problem \eqref{12} has a positive periodic solution.
\end{theorem}

\begin{proof} For fixed $\lambda\in(0, \lambda^*)$, there exists a
 positive constant $\theta$ such that
$$  
\lambda<\frac{1-K}{2
c^{2}\int_0^{T} b(t)dt}\frac{\theta^{2}}{ \max_{|\xi|\leq
\theta}G(\xi)}.  
$$
Moreover, by \eqref{13}, we can choose $\eta >0$ satisfying 
$\eta< \frac{\theta}{c}\sqrt{\frac{1-K}{DM(1+K)}}$ such that
$$  
\frac{(1+K)DM} {2\lambda\int_{t_{1}}^{t_{p}}b(t)dt} 
< \frac{G(\eta\mathcal{\varepsilon})}{\eta^{2}}, 
$$
where $\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$. 
The conclusion then follows from Corollary \ref{c5}.
\end{proof}

The following examples illustrates some of our results.

\begin{example}  \rm
Take $N=1$ and consider the problem
\begin{equation} \label{14}
\begin{gathered} 
-u''(t)+u(t)= \lambda b(t)g(u(t))+h(u(t)), \quad \text{a.e. } t\in [0,3],
\\
\Delta (u'(t_j))=I_j(u(t_j)),\quad  j=1,2,\\
u(0)- u(3)= u'(0)- u'(3)=0,
\end{gathered}
\end{equation}
where $b(t)=e^t$ for every $t\in[0,3]$, $t_1=1$, $t_2=2$,
$g(x)=1+3x|x|e^{x^4}+4x|x|^5e^{x^4}$, $h(x)=\frac{1}{12}x^+$,
$x^+=\max\{x,0\}$, and $I_j(x)=\frac{1}{8}x$ for $j=1,2$ for every
$x\in\mathbb{R}$. It is easy to see that
$G(x)=x+|x|^3e^{x^4}$, and
$$  
\lim_{x\to 0^{+}}\frac{ \max_{|\xi|\leq x}G(\xi)}{x^{2}}=+\infty.  
$$
Moreover, since $c=\sqrt{6}$, $L=\frac{1}{72}$, and
$L_{1j}=\frac{1}{48}$ for $j=1,2$, we see that $K=\frac{3}{4}<1$.
Hence, applying Theorem \ref{t5}, for each 
$\lambda\in\Big(0,\frac{1}{48(e^3-1)(1+e)}\Big)$,
problem \eqref{14} has a positive periodic solution.
\end{example}

\begin{example} \rm
Let $N=2$, $p=2$, $T=3$, $t_1=1$, and $t_2=2$. Let $A:[0,3]\to
\mathbb{R}^{2\times 2}$ be the identity matrix, let
$G(\xi_1,\xi_2)=\xi_1+\xi_2+\frac{1}{4}\xi_1^4+\frac{1}{4}\xi_2^4$
for all $(\xi_1,\xi_2)\in \mathbb{R}^2$, $b\in L^1([0,3])$ be a
positive function, $I_{ij}(s)=\frac{1}{96}s(1+e^{-s})$ for all
$s\in\mathbb{R}$, for $i=1,2$ and $j=1,2$, and
$H(\xi_1,\xi_2)=\frac{1}{72\sqrt{2}}
(\frac{1}{2}\xi_1^2+\xi_1\xi_2+\frac{1}{2}\xi_2^2)$
for all $(\xi_1,\xi_2)\in \mathbb{R}^2$. It is clear that
$$  
\lim_{x\to 0^{+}}\frac{ \max_{|\xi|\leq x}G(\xi)}{x^{2}}=+\infty.  
$$
Moreover, since $c=\sqrt{6}$,
$L=\frac{1}{72\sqrt{2}}$ and $L_{ij}=\frac{1}{96}$ for 
$i=1,2$, $j=1,2$, we have $K=\frac{2+\sqrt{2}}{4\sqrt{2}}<1$. Hence,
applying Theorem \ref{t5}, for each
$$ 
\lambda\in\Big(0,
\frac{ 1-\frac{2+\sqrt{2}}{4\sqrt{2}}}{ 12\int_0^{T}
b(t)dt}\sup_{\theta>0}\frac{\theta^{2}}{ \max_{|\xi|\leq
\theta}\big(\xi_1+\xi_2+\frac{1}{4}\xi_1^4+\frac{1}{4}\xi_2^4\big)}\Big),  
$$
problem \eqref{12} has a positive periodic solution.
\end{example}

Our next theorem is for the existence of three positive
periodic solutions to problem \eqref{12}. It is based on Corollaries \ref{c5} 
and \ref{c6}.
We use a combination of algebraic conditions on the functions $G$ and $H$ 
that give two local minimums for the functional $I_{\lambda}$, and then we 
apply the Pucci-Serrin mountain pass lemma to obtain
the third solution.

\begin{theorem}  \label{t7}
Let {\rm (B5)}  hold and assume there
exist constants $\theta >0$, $\eta >0$, $\bar{\theta} >0$, and
$\bar{\eta} >0$ with
$$ 
 c\sqrt{\frac{DM(1+K)}{1-K}}\eta
<\theta\leq\overline{\theta}<c\sqrt{Dm}\overline{\eta}  
$$
such that {\rm (A7)} and {\rm (B4)} hold.
If
\begin{equation}  \label{15}
\frac{ \max_{|\xi|\leq \theta}G(\xi) \int_0^{T} b(t)dt
}{\theta^{2}}<\frac{1-K}{2{c}^{2}}\frac{
\max_{|\xi|\leq \bar{\theta}}G(\xi) \int_0^{T} b(t)dt-
G(\bar{\eta}\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt}{
\frac{1}{2}(1-K)
(\frac{\bar{\theta}}{c})^{2}-\frac{1}{2}(1+K)DM\bar{\eta}^2},
\end{equation}
where
$\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$,
then, for each
$$  
\lambda\in\Lambda:=\Big(\max\Big\{\hat{\lambda},
\frac{(1+K)DM\eta^2}{2
G(\eta\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt}\Big\},\,
\frac{(1-K)(\frac{\theta}{c})^{2}}{2 \max_{|\xi|\leq
\theta}G(\xi)\int_0^{T} b(t)dt }\Big),
$$
where $\hat{\lambda}$ is given in \eqref{new-201}, problem \eqref{12} 
has at least three positive periodic solutions.
\end{theorem}

\begin{proof}
First we observe that \eqref{15} implies $\Lambda\neq\emptyset$.
Fix $\lambda\in\Lambda$. Using Corollary \ref{c5}, we obtain the
first positive periodic solution $u^{\ast}$ as a local minimum of the 
functional $I_\lambda$ with
$$   
\frac{1}{2}\|u^{\ast}\|^{2}
+\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{u^{\ast}_i(t_j)}I_{ij}(s)ds
-\int_0^{T}H(u^{\ast}(t))dt
<\frac{1}{2}(1-K)\big(\frac{\theta}{c}\big)^{2}. 
$$
Corollary \ref{c6} guarantees a second positive periodic solution $\bar{u}$ with
$$  
\frac{1}{2}\|\bar{u}\|^{2}
+\sum_{j=1}^{p}\sum_{i=1}^N\int_0^{\bar{u}_i(t_j)}I_{ij}(s)ds
-\int_0^{T}H(\bar{u}(t))dt
>\frac{1}{2}(1-K)\big(\frac{\bar{\theta}}{c}\big)^{2}.  
$$
The mountain pass theorem of Pucci and Serrin (\cite{PS}) then ensures 
the existence of a third positive periodic solution.
\end{proof}

As a consequence of Theorem \ref{t7} we have the following result.

\begin{theorem}  \label{t8}
Assume that
\begin{gather}  \label{16}
\limsup_{|x|\to 0^+}\frac{\max_{|\xi|\leq x}G(\xi)}{|x|^2}=+\infty, \\
 \label{17}
\limsup_{|\xi|\to +\infty }\frac{G(\xi)}{|\xi|^{2}}=0,
\end{gather}
and there are constants $\bar{\theta} >0$ and $\bar{\eta} >0$
with $\frac{\bar{\theta}}{c\sqrt{Dm}}<\bar{\eta}$ such that
\begin{equation}  \label{18}
\frac{ \max_{|\xi|\leq \bar{\theta}}G(\xi) \int_0^{T}
b(t)dt }{\bar{\theta}^{2}}<\frac{1-K}
{c^{2}(1+K)DM}\frac{
G(\bar{\eta}\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt}{
\bar{\eta}^2}
\end{equation}
where $\mathcal{\varepsilon}=(1,0,\dots,0)\in\mathbb{R}^N$.
Then, for each
$$   
\lambda \in \Big(\frac{(1+K)DM\bar{\eta}^2}{
2G(\bar{\eta}\mathcal{\varepsilon})\int_{t_{1}}^{t_{p}}b(t)dt},\
\frac{ (1-K)\bar{\theta}^{2}}{ 2c^{2} \max_{|\xi|\leq
\bar{\theta}}G(\xi)\int_0^{T}b(t)dt}\Big),
$$
\eqref{12} has at least three positive periodic solutions.
\end{theorem}

\begin{proof} 
We can easily observe from \eqref{17} that (B5) is satisfied.
Moreover, by choosing $\eta$ small
enough and $\theta=\bar{\theta}$, we see that \eqref{16} implies condition 
(A7) holds, and \eqref{18} implies
(B4) and \eqref{15} hold. We thus have the conclusion of the theorem.
\end{proof}

In conclusion, we would like to mention that we believe that this approach 
of combining techniques to obtain multiple
solutions of boundary value problems, with or without impulses, 
will prove to be a valuable strategy.

\subsection*{Future research}
 One direction for future work would be to extend the results here to the case 
where the right hand side of the equation in \eqref{1} is not continuous. 
In this regard, we refer the reader to the paper 
of Molica Bisci and Repov\v s \cite{New3}. Another possibility
is to consider the situation where the equation in \eqref{1} is 
replaced by an inclusion.

\subsection*{Acknowledgments}
This article was written while the second author was visiting The University
of Tennessee at Chattanooga.


\begin{thebibliography}{99}

\bibitem{AO} R. P. Agarwal, D. O'Regan;
\emph{Multiple nonnegative solutions for second order impulsive differential 
equations,}  Appl. Math. Comput. 114 (2000), 51--59.

\bibitem{BD1} L. Bai, B. Dai;
\emph{Application of variational method to a class of Dirichlet boundary 
value problems with  impulsive effects,} J. Franklin Inst. 348 (2011), 2607--2624.

\bibitem{BS} D. Bainov, P. Simeonov;
\emph{Systems with Impulse Effect, Ellis
Horwood Series: Mathematics and Its Applications,} Ellis Horwood, Chichester, 1989.

\bibitem{BHN} M. Benchohra, J. Henderson, S. Ntouyas;
\emph{Theory of Impulsive Differential Equations,}
 Contemporary Mathematics and Its Applications, Vol 2, Hindawi, New York, 2006.

\bibitem{B} G. Bonanno;
 \emph{A critical point theorem via the Ekeland variational principle,} 
Nonlinear Anal. 75 (2012), 2992--3007.

\bibitem{CH}H. Chen, Z. He;
 \emph{New results for perturbed Hamiltonian systems with impulses,}
Appl. Math. Comput. 218 (2012), 9489--9497.

\bibitem{C} G. Cordaro;
 \emph{Three periodic solutions to an eigenvalue problem for a class of second order
Hamiltonian systems,} Abstr. Appl. Anal. 18 (2003), 1037--1045.

\bibitem{CR} G. Cordaro, G. Rao;
\emph{Three periodic solutions for perturbed second
order Hamiltonian systems,} J. Math. Anal. Appl. 359 (2009), 780--785.

\bibitem{CES} V. Coti-Zelati, I. Ekeland, E. Sere;
\emph{A variational approach to homoclinic orbits
in Hamiltonian systems,} Math. Ann. 288 (1990), 133--160.

\bibitem{D1} G. D'Agu\`{\i}, G. Molica Bisci;
\emph{Three non-zero solutions for elliptic Neumann problems,}
 Anal. Appl. 9 (2011), 383--394.

\bibitem{D2} G. D'Agu\`{\i};
\emph{Multiplicity results for nonlinear mixed boundary value problem},
 Bound. Value Prob. 2012 (2012), No. 134.

\bibitem{F} F. Faraci;
 \emph{Multiple periodic solutions for second order systems with changing sign 
potential}, J. Math. Anal. Appl. 319 (2006), 567--578.

\bibitem{FN} D. Franco, J. J. Nieto;
\emph{Maximum principle for periodic impulsive first order problems,} 
J. Comput. Appl. Math. 88 (1998) 149--159.

\bibitem{GHK} J. R. Graef, S. Heidarkhani, L. Kong;
\emph{Infinitely many periodic solutions to a class of perturbed second-order
impulsive Hamiltonian systems,} preprint.

\bibitem{GHO} J. R. Graef, J. Henderson, A. Ouahab;
\emph{Impulsive Differential Inclusions, A Fixed Point Approach},
De Gruyter Series in Nonlinear Analysis and Applications Vol. 20, 
De Gruyter, Berlin, 2013.

\bibitem{GA} H. Gu, T. An;
\emph{Existence of infinitely many periodic solutions for second-order Hamiltonian
systems,} Electron. J. Differential Equations, vol. 2013 (2013), No. 251, pp. 1--10.

\bibitem{HFK} S. Heidarkhani, M. Ferrara, S. Khademloo;
\emph{Nontrivial solutions for one-dimensional fourth-order
Kirchhoff-type equations,} Mediterr. J. Math., to appear.

\bibitem{IJ} M. Izydorek, J. Janczewska;
\emph{Homoclinic solutions for a class of second order
Hamiltonian systems,} J. Differential Equations 219 (2005), 375--389.

\bibitem{LBS} V. Lakshmikantham, D. D. Bainov, P. S. Simeonov;
\emph{Theory of Impulsive Differential Equations,} World Scientific, Singapore, 1989.

\bibitem{LS} F.-f. Liao, J. Sun;
\emph{Variational approach to impulsive problems: A
survey of recent results,} Abstr. Appl. Anal. Volume 2014, Article ID 382970,
 11 pages.

\bibitem{LJ} X. N. Lin, D. Q. Jiang;
\emph{Multiple positive solutions of Dirichlet boundary value problems 
for second order impulsive differential equations,} 
J. Math. Anal. Appl. 321 (2006), 501--514.

\bibitem{MW} J. Mawhin, M. Willem;
\emph{Critical Point Theory and Hamiltonian Systems,} Springer, New York, 1989.

\bibitem{New5} G. Molica Bisci;
\emph{Fractional equations with bounded primitive,}  
Appl. Math. Letters 27 (2014), 53--58.

\bibitem{New2} G. Molica Bisci, D. Repov\v s;
\emph{Multiple solutions for elliptic equations involving a general 
operator in divergence form,} Ann. Acad. Fenn. Math. 39 (2014), 259--273.

\bibitem{New3} G. Molica Bisci, D. Repov\v s;
\emph{Nonlinear algebraic systems with discontinuous terms,} 
J. Math. Anal. Appl. 398 (2013), 846--856.

\bibitem{NO} J. J. Nieto, D. O'Regan;
\emph{Variational approach to impulsive differential equations,}
 Nonlinear Anal. Real World Appl. 10 (2009), 680--690.

\bibitem{PS} P. Pucci, J. Serrin;
\emph{A mountain pass theorem,} J. Differential Equations 60 (1985) 142-149.

\bibitem{Ra1} P. H. Rabinowitz;
\emph{Homoclinic orbits for a class of Hamiltonian systems,} 
Proc. Roy. Soc. Edinb. 114 (1990), 33--38.

\bibitem{Ra} P. H. Rabinowitz;
\emph{Minimax Methods in Critical Point Theory
with Applications to Differential Equations,} 
CBMS Reg. Conf. Ser. Math. 65, Amer. Mat. Soc., Providence, 1986.

\bibitem{Ra2} P. H. Rabinowitz;
\emph{Variational methods for Hamiltonian systems,} in: Handbook of
Dynamical Systems, vol. 1, North-Holland, 2002, Part 1, 
Chapter 14, pp. 1091--1127.

\bibitem{R1} B. Ricceri;
\emph{A general variational principle and some of its applications,} 
J. Comput. Appl. Math. 113 (2000), 401--410.

\bibitem{SP} A. M. Samoilenko, N. A. Perestyuk;
\emph{Impulsive Differential Equations,} World Scientific, Singapore, 1995.

\bibitem{SCN} J. Sun, H. Chen, J. J. Nieto;
 \emph{Infinitely many solutions for second-order Hamiltonian system with 
impulsive effects,} Math. Comput. Model. 54 (2011), 544--555.

\bibitem{SCNO} J. Sun, H. Chen, J. J. Nieto, M. Otero-Novoa;
\emph{The multiplicity of solutions for perturbed
second-order Hamiltonian systems with impulsive effects,} 
Nonlinear Anal. 72 (2010), 4575--4586.

\bibitem{T} C. Tang;
\emph{Periodic solutions for nonautonomous second order systems with
sublinear nonlinearity,} Proc. Amer. Math. Soc. 126 (1998), 3263--3270.

\bibitem{TG} Y. Tian, W. Ge;
\emph{Applications of variational methods to boundary-value problem
for impulsive differential equations,} Proc. Edinb. Math. Soc. 51 (2008), 509--527.

\bibitem{WZ} Z. Wang, J. Zhang;
\emph{Periodic solutions of a class of second order non-autonomous
Hamiltonian systems,} Nonlinear Anal. 72 (2010), 4480--4487.

\bibitem{Z} E. Zeidler;
\emph{Nonlinear Functional Analysis and its Applications,} 
Vol. II, Springer, Berlin-Heidelberg-New York, 1985.

\bibitem{ZD} D. Zhang, B. Dai;
\emph{Existence of solutions for nonlinear impulsive differential equations 
with Dirichlet boundary conditions,} Math. Comput. Modelling 53 (2011), 1154--1161.

\bibitem{ZL2} Q. Zhang, C. Liu;
\emph{Infinitely many periodic solutions for second order
Hamiltonian Systems,} J. Differential. Eqs. 251 (2011), 816--833.

\bibitem{ZT} Q. Zhang, X. Tang;
\emph{New existence of periodic solutions for second order non-autonomous 
Hamiltonian systems,} J. Math. Anal. Appl. 369 (2010), 357--367.

\bibitem{ZhL} J. Zhou, Y. Li;
\emph{Existence of solutions for a class
of second order Hamiltonian systems with impulsive effects,}
Nonlinear Anal. 72 (2010), 1594--1603.

\bibitem{ZoL} W. Zou, S. Li;
\emph{Infinitely many solutions for Hamiltonian systems,} 
J. Differential Equations 186 (2002), 141--164.

\bibitem{ZZ} X. Zhang, Y. Zhou;
\emph{Periodic solutions of non-autonomous second order Hamiltonian systems,} 
J. Math. Anal. Appl. 345 (2008), 929--933.

\end{thebibliography}

\end{document}
