\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2014 (2014), No. 136, pp. 1--14.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2014 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2014/136\hfil Second-order differential inclusions]
{Periodic solutions for second-order differential inclusions
 with nonsmooth potentials under weak AR-conditions}

\author[L. Chen, Q. Zhang, G. Li \hfil EJDE-2014/136\hfilneg]
{Lizhen Chen, Qinghua Zhang, Gang Li}  % in alphabetical order

\address{Lizhen Chen \newline
Department of Applied Mathematics,
Shanxi University of Finance and Economics,
Taiyuan, Shanxi, 030006, China}
\email{chenlz409@126.com}

\address{Qinghua Zhang (corresponding author)\newline
Department of Mathematics,
Nantong University,
Nantong, Jiangsu 226007, China}
\email{zhangqh1971@126.com}

\address{Gang Li \newline
Department of Mathematics,
Yangzhou University,
Yangzhou, Jiangsu 225002, China}
\email{yzgangli@163.com}

\thanks{Submitted  April 3, 2014. Published June 11, 2014.}
\subjclass[2000]{34K37, 46E35, 47H10}
\keywords{Sobolev space; periodic solution; locally Lipschitz potential;
\hfill\break\indent AR-condition; nonsmooth C-condition; 
the least action principle; mountain pass lemma}

\begin{abstract}
 In this article, we study a periodic second-order differential
 inclusions with locally Lipschitz potentials. By means of the least
 action principle and the minimax principle of nonsmooth type, we
 prove the existence of two  nontrivial periodic solutions
 under the weak AR-conditions. The method developed in this paper
 can be applied for studying second-order differential inclusions
 of periodic type, and for elliptic equations with Neumann boundary
 condition.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{remark}[theorem]{Remark}
\allowdisplaybreaks


\section{Introduction and preliminaries}

The AR-condition, which was introduced  by
Abrosstti-Rabinowitz (see Rabinowitz \cite{ra}), is used to
investigate the PS-condition (or C-condition) for an energy
functional, which probably has critical points (refer to
Papageorgiou-Papageorgiou \cite{pp}, Hu-Papageorgiou \cite{hup}, and
Nikolaos \cite{ni}). A nonsmooth function
$f:T\times\mathbb{R}^{N}\to\mathbb{R}$ of Caratheodory type,
which growth is constrained by the $r-$th power of $|x|$ ($r>1$), is
said to be satisfying the AR-condition, if
\begin{itemize}
\item[(H0)]  there exist $\mu>p>1$ and $M_0>0$ such that for a.e.
$t\in T$ and all $x\in\mathbb{R}^{N}$ with $|x|\geq M_0$, the
inequality $0<\mu f(t,x)\leq(u,x)$, or equivalently $\mu
f(t,x)\leq-f^{0}(t,x,-x)$ (see \cite{pp}) holds for all
$u\in\partial f(t,x)$.
\end{itemize}
It is well known that (cf. \cite{kop},  \cite[p. 9]{ra}, or
\cite[p. 93]{mw}), any function $f(t,x)$ satisfying the AR-condition is
$p$-superlinear definitely; i.e., there are $\beta_{i}>0$, $i=1,2$,
such that $f(t,x)\geq\beta_{1}|x|^{\mu}-\beta_{2} $ for almost all $t\in
T$ and all $x\in\mathbb{R}^{N}$.

If $r=p$, then the AR-condition can be weakened as (please refer to
Kourogenis-Papageorgiou \cite{kop})
\begin{itemize}
\item[(H1)]
\[
\liminf_{|x|\to\infty}\frac{(u,x)-pf(t,x)}{|x|^{\alpha}}>0,
\quad \alpha\in(0,p)
\]
uniformly for all $u\in\partial f(t,x)$,  and a.e. $t\in T$,
\end{itemize}
or even as
\begin{itemize}
\item[(H2)]
\[
\liminf_{|x|\to\infty}((u,x)-pf(t,x))=+\infty,
\]
uniformly for all $u\in\partial f(t,x)$,  and a.e. $t\in T$ (cf
Li-Zhou \cite{lz}).
\end{itemize}

In this article, we study the periodic problem
\begin{equation} \label{e1.1}
\begin{gathered}
-(J_{p}(x'(t)))'+a(t)J_{p}(x(t))\in\partial f(t,x(t))\quad
\text{a.e. on }T=[0,b],\\
x(0)=x(b),\quad x'(0)=x'(b),\quad 1<p<+\infty
\end{gathered}
\end{equation}
where $J_{p}$ denotes the $p$-Laplacian defined by
\[
J_{p}(x)=\begin{cases}
|x|^{p-2}x,& \text{if }x\neq0,\\
0,&\text{if } x=0,
\end{cases}
\]
and the nonsmooth potential $f(t,x)$ satisfying the AR-like
condition
\begin{itemize}
\item[(H3)] there is a $d_{*}>0$, such that for a.e. $t\in T$, inequality
\begin{equation}
(u,x)-pf(t,x)\geq d_{*}
\end{equation}
 for all $u\in\partial f(t,x)$ with $|x|\geq M_0$.
\end{itemize}
Since condition (H3) is weaker than all those mentioned above,
it is not used here to verify the C-condition for the energy functional anymore. 
Instead, an extra condition (see $H(f)_{1}$ $(vii)$) which describe 
the asymptotic property of the potentials is taken into account. 
In this situation, the weak AR-condition is used only to investigate the asymptotic
property of the energy functional in the infinite, and the former
role of it to verify the nonsmooth C-condition has been replaced by
the additional one.

For the convenience of the reader, we firstly overview briefly the
theory of nonsmooth analysis, which is established by Clarke \cite{cl}.

Let $X$ be a Banach space, on which, we define a real function $f$,
which is said to be locally Lipschitz at $x$ (or equivalently,
Lipschitz near $x$) with the rank $L$, if there is a neighborhood
$U$ of $x$, such that
$$
|f(x_{1})-f(x_{2})|\leq L\|x_{1}-x_{2}\|
$$
for all $x_{i}\in U$, $i=1,2$.

Attached to any locally Lipschitz function $f$, we can define the
Clarke derivative at $x$ along the direction $v\in X$:
\[
f^{\circ}(x,v)=\limsup_{y\to x,\;
\lambda\downarrow0}\frac{f(y+\lambda v)-f(y)}{\lambda}.
\]

Let $v$ run over the space $X$, we get a function
$f^{\circ}(x,\cdot)$, which is subadditive, positively homogeneous
and uniformly Lipschitz with the rank $L$. So there exists a
$w^{*}-$compact and convex subset $\partial f(x)$ of $X^{*}$, for
which, $f^{\circ}(x,\cdot)$ is the corresponding support function;
i.e.,
\[
\partial f(x)=\{\xi\in X^{*}:\langle\xi,v\rangle\leq f^{\circ}(x,v),\;
 \forall v\in X\},
\]
or conversely,
\[
f^{\circ}(x,v)=\sup\{\langle\xi,v\rangle:\xi\in\partial f(x)\}
\]
with the supermum being reached for each $v\in X$. In this
definition, $\partial f(x)$ is called the Clarke subdifferential of
$f$ at the point $x$.

Here we list two properties of $\partial f(x)$:
\begin{itemize}
\item[(1)] $\partial(\lambda f)(x)=\lambda\partial f(x)$,
for all $\lambda\in\mathbb{R}$, $x\in X$,
\item[(2)] if $g$ is another function locally Lipschitz on $X$, then
for all $x\in X$, the following inclusion holds
\[
\partial(f+g)(x)\subseteq\partial f(x)+\partial g(x),
\]
and equality holds if one the functions $f,g$ is strictly
differentiable.
\end{itemize}
For other properties, please refer to \cite{cl}.

In what follows, we need the nonsmooth C-condition (introduced by
Cerami originally) for a locally Lipschitz function $f$; that is,
any sequence $\{x_{n}\}\in X$ with
\begin{itemize}
\item[(1)]  $\{\varphi(x_{n})\}$ bounded, and

\item[(2)]  $(1+\|x_{n}\|)m(x_{n})\to0$, as $n\to\infty$
\end{itemize}
contains a strongly convergent subsequence.

In this definition, the notation $m(x)$ denotes
$\inf_{\xi\in\partial\varphi(x)}\|\xi\|_{*}$, which can be attained
by some $\xi\in\partial\varphi(x)$ because $\partial\varphi(x)$ is
$w^{*}-$compact and the norm $\|\cdot\|_{*}$ is $w^{*}-$lsc (cf.
\cite{kop}).

This article is organized as follows. 
In Sections 2,3, we pay attention
to the periodic differential system of scalar type, where the
potential function $f$ satisfying the AR-like condition (H1) or
(H3). By truncating the potential function, investigating the
asymptotic properties of the corresponding energy functional in the
infinite area, and using the critical point theory of nonsmooth type
developed by Chang \cite{cha} and Kourogenis-Papageorgiou
\cite{kop}, we find two distinct positive solutions for this
problem. And in Section 4, we turn to deal with the periodic
differential inclusion in $\mathbb{R}^{N}$ ($N\geq1$) under the
weaker AR-conditions, and deduce the multiple existence of
solutions. Our work space in this paper is
$$
W_{\rm per}^{1,p}(T,\mathbb{R}^{N})=\{x\in W^{1,p}
(T,\mathbb{R}^{N}):x(0)=x(b)\}.
$$
For other treatment of positive solutions for the periodic
differential systems, please refer to \cite{pp2} and \cite{cs}.


\section{Properties of the energy functional}

Firstly, we give some hypotheses on $f$ and $a$, which will be used
in the following paragraphs.
\begin{itemize}
\item[(HF)] $f: T\times\mathbb{R}\to\mathbb{R}$ is a real
function satisfying $f(\cdot,0)\in L^{\infty}(T,\mathbb{R}_{+})$,
and $0\in\partial f(t,0)$ a.e. on $T$, and
\begin{itemize}
\item[(i)]  for each $x\in\mathbb{R}$, the function $t\to f(t,x)$
is measurable,

\item[(ii)] as the function of $x$, $f(t,x)$ is locally Lipschitz for
a.e. $t\in T$, and

\item[(iii)] for almost all $t\in T$, and all $x\in\mathbb{R}$, 
\begin{equation} \label{e2.1}
|u|\leq a_0(t)+d_0|x|^{p-1}
\end{equation}
for all $u\in\partial f(t,x)$, where $a_0\in
L^{\infty}(T,\mathbb{R}_{+})$, and $d_0>0$,

\item[(iv)] there are constants $M_0,d>0$ (without loss of generality,
we may assume $M_0\geq1$ here) and $0<\alpha<p$, such that
\begin{equation}
ux-pf(t,x)\geq dx^{\alpha}
\end{equation}
for a.e. $t\in T$, and all $u\in\partial f(t,x)$ with $x\geq M_0$,
or weakly

\item[(iv')] for a.e. $t\in T$, 
\begin{equation}
ux-pf(t,x)\geq d
\end{equation}
for all $u\in\partial f(t,x)$ with $x\geq M_0$,

\item[(v)] $\limsup_{x\to0^{+}}x^{-p}f(t,x)\leq0$, uniformly for
a.e $t\in T$,

\item[(vi)] there exists $x_0\geq M_0$ satisfying
\begin{equation} \label{e2.4}
\frac{7}{6}\int_0^{b}f(t,x_0)dt>\frac{1}{p}\|a\|_{1}x_0^{p},
\end{equation}

\item[(vii)]
\begin{equation} \label{e2.5}
x_0^{p}\sup_{x\geq x_0}\frac{f(t,x)}{x^{p}}\leq\frac{3}{4}(\frac{d}{p+2}+\frac{1}{b}\int_0^{b}f(t,x_0)dt)
\end{equation}
uniformly for a.e. $t\in T$.
\end{itemize}

\item[(HA)] $a\in L^{1}(T,\mathbb{R}_{+})$ and
$\int_0^{b}a(t)dt>d_0$.
\end{itemize}

\begin{remark} \label{rem2.1} \rm
Combining (iii) with the mean value theorem for locally
Lipschitz function (please refer to \cite[Theorem 2.3.7]{cl}),
we can deduce the following estimate for $f(t,x)$:
\begin{equation} \label{e2.6}
|f(t,x)|\leq a_0(t)|x|+\frac{d_0}{p}|x|^{p},
\end{equation}
and then for each $\varepsilon>0$, there is a constant $C_{\varepsilon}>0$, 
such that
\begin{equation} \label{e2.7}
|f(t,x)|\leq C_{\varepsilon}+\frac{d_0+\varepsilon}{p}|x|^{p},
\end{equation}
for a.e. $t\in T$, and all $x\in\mathbb{R}$.
\end{remark}

\begin{remark} \label{rem2.2} \rm
Some properties of $f$ can be derived from  (HF) (v) and  (vii).
The first one is $f(t,0)\leq0$, and hence $f(t,0)=0$ a.e. on $T$,
since $f(t,0)$ is nonnegative. And for every $\varepsilon>0$ and $\mu>p$, there
is a $C_{\varepsilon,\mu}>0$ such that
\begin{equation}
f(t,x)\leq\varepsilon x^{\mu}+C_{\varepsilon,\mu}
\end{equation}
for a.e. $t\in T$ and all $x\in\mathbb{R}$.
\end{remark}

\begin{remark} \label{rem2.3} \rm
Due to hypothesis (HA), we can deduce by contradiction that, for
each nonnegative number $\varepsilon$ smaller than
$\int_0^{b}a(t)dt$, there is $c_{\varepsilon}>0$ (depends on
$\varepsilon$) for which, the following inequality
\begin{equation}
\frac{1}{p}\|x'\|_{p}^{p}+\frac{1}{p}\int_0^{b}(a(t)-\varepsilon)|x(t)|^{p}dt\geq
c_{\varepsilon}\|x\|_{1,p}^{p},
\end{equation}
holds on $W_{\rm per}^{1,p}(T)$.
\end{remark}

Let us define the truncated function $f_{1}$ of $f$ as follows:
\begin{equation}
f_{1}(t,x)=\begin{cases}
f(t,x),& \text{if }x\geq0,\\
0,&\text{if } x<0.
\end{cases}
\end{equation}
Evidently, $f_{1}(\cdot,x)$ is measurable for each $x\in\mathbb{R}$,
and $f_{1}(t,\cdot)$ is locally Lipschitz for a.e $t\in T$. In view
of the chain rules (cf \cite[Theorem 2.3.9]{cl}), we have
\begin{equation}
\partial f_{1}(t,x)=\begin{cases}
\partial f(t,x),& \text{if }x>0,\\
co([0,1]\cdot\partial f(t,x)),& \text{if }x=0,\\
0,&\text{if }x<0
\end{cases}
\end{equation}
a.e.on $T$, where $[0,1]\cdot\partial f(t,x)=\{\lambda
y:\lambda\in[0,1], y\in\partial f(t,0)\}$. Thus $f_{1}$ has all the
properties as $f$ has under hypotheses (HF).

If we define $\phi(x)=\int_0^{b}f_{1}(t,x(t))dt$,
$ x\in L^{p}(T)$, then we get a locally Lipschitz functional with its
subdifferential $\partial\phi(x)\subseteq\int_0^{b}\partial
f_{1}(t,x(t))dt$, which means for every $\xi\in\partial\phi(x)$,
there exists a $L^{q}(T)$ selection $u$ of 
$\partial f_{1}(\cdot,x(\cdot))$ ($\frac{1}{p}+\frac{1}{q}=1$), such that
\begin{equation}
(\xi,y)_{p,q}=\int_0^{b}u(t)y(t)dt
\end{equation}
for all $y\in L^{p}(T)$ (see \cite[Theorem 2.7.5]{cl}).

Since $W_{\rm per}^{1,p}(T)$ is a dense subspace of $L^{p}(T)$ with the
inclusion mapping continuous, if we set
$\phi_0=\phi|_{W_{\rm per}^{1,p}(T)}$, then we obtain another locally
Lipschitz functional $\phi_0$ defined on $W_{\rm per}^{1,p}(T)$ with
$\partial\phi_0(x)=\partial\phi(x)=S_{\partial
f_{1}(\cdot,x(\cdot))}^{q}$ holding for all $x\in W_{\rm per}^{1,p}(T)$
(see the corollary of  \cite[Theorem 2.3.10]{cl}).

Now, we can consider the energy functional on $W_{\rm per}^{1,p}(T)$
associated with the periodic system \eqref{e1.1}:
\begin{eqnarray}
\Phi(x)=\frac{1}{p}\int_0^{b}|x'(t)|^{p}dt
+\frac{1}{p}\int_0^{b}a(t)|x(t)|^{p}dt-\int_0^{b}f_{1}(t,x(t))dt.
\end{eqnarray}
Based on the discussion in \cite{zh} and the finite sum rule (see
Section 1), we know that $\Phi$ is locally Lipschitz and its Clarke
subdifferential can be represented by
\begin{equation}
\partial\Phi(x)=\mathcal{J}_{p}x+a(\cdot)J_{p}(x(\cdot))-\partial\phi_0(x),
\end{equation}
where
$\mathcal{J}_{p}:W_{\rm per}^{1,p}(T)\to(W_{\rm per}^{1,p}(T))^{*}$
is a demicontinuous and $(S)_{+}$ operator defined by
\[
\langle\mathcal{J}_{p}x,y\rangle
=\int_0^{b}J_{p}(x'(t))y'(t)dt,\quad \forall\;x,y\in W_{\rm per}^{1,p}(T).
\]

\begin{lemma} \label{lem2.4} 
Under hypotheses {\rm (HF) (i)--(iv)} and {\rm (HA)}, the energy functional
$\Phi$ satisfies the nonsmooth $C$-condition.
\end{lemma}

\begin{proof}
 Suppose $\{\Phi(x_{n})\}$ is bounded and
$(1+\|x_{n}\|_{1,p})m(x_{n})\to0$ as $n\to\infty$
for some sequence $\{x_{n}\}$ of $W_{\rm per}^{1,p}(T)$, we will show
that $\{x_{n}\}$ has a convergent subsequence. For this purpose, we
need to introduce two truncated functions for any function $x$.
\begin{equation}
x^{+}=\max\{0,x\},\quad x^{-}=\max\{0,-x\}.
\end{equation}

As we know that, both $x^{+}$ and $x^{-}$ are all in
$W_{\rm per}^{1,p}(T)$ whenever $x\in W_{\rm per}^{1,p}(T)$ (please refer to
Marcus-Mizel \cite{mm} or Evans-Gariepy \cite{eg}, p.130). Moreover,
for the derivatives of $x^{+}$ and $x^{-}$, we have
\begin{gather*}
(x^{+})'(t)=\begin{cases} 
x'(t),&\text{if }x>0,\\
0,&\text{if }x\leq0,
\end{cases} \\
(x^{-})'(t)=\begin{cases}
0,&\text{if }x\geq0,\\
-x'(t),&\text{if }x<0.
\end{cases}
\end{gather*}
Consequently, the following relations between $x$ and $x^{+},x^{-}$
hold:
\begin{gather}
\max\{\|x^{+}\|_{1,p},\|x^{-}\|_{1,p}\}\leq\|x\|_{1,p},\quad
\|x\|_{1,p}\leq\|x^{+}\|_{1,p}+\|x^{-}\|_{1,p}, \label{e2.19}\\
\int_0^{b}J_{p}(x'(t))(x^{+})'(t)dt
 =\int_{\{x>0\}}|x'(t)|^{p}dt=\int_0^{b}|(x^{+})'(t)|^{p}dt, \label{e2.20}\\
\int_0^{b}J_{p}(x'(t))(x^{-})'(t)dt=-\int_{\{x<0\}}|x'(t)|^{p}dt
=-\int_0^{b}|(x^{-})'(t)|^{p}dt. \label{e2.21}
\end{gather}

Since the subdifferential $\partial\Phi(x)$ is a $w^{*}$-compact
subset of $(W_{\rm per}^{1,p}(T))^{*}$, its minimal norm can be reached
by some elements in it. Hence for each $n\in\mathbb{N^{*}}$, there
is $\xi_{n}\in\partial\Phi(x_{n})$, such that
\begin{gather}
\|\xi_{n}\|_{*}=m(x_{n})\to0, \label{e2.22}\\
\max\{|\langle\xi_{n},x_{n}\rangle|,|\langle\xi_{n},x_{n}^{-}\rangle|\}
\leq\|x_{n}\|_{1,p}m(x_{n})\to 0 \label{e2.23}
\end{gather}
as $n\to\infty$. Therefore we can find constant $M>0$,
sequences of positive numbers $\{\varepsilon_{n}\}$ with
$\varepsilon_{n}\downarrow0$, and functions  $\{u_{n}\}$ with
$u_{n}\in S_{\partial f_{1}(\cdot,x_{n}(\cdot))}^{q}$, satisfying
\begin{gather}
-Mp\leq\int_0^{b}|x_{n}'(t)|^{p}dt+\int_0^{b}a(t)|x_{n}(t)|^{p}dt
 -\int_0^{b}pf_{1}(t,x_{n}(t))dt\leq Mp, \label{e2.24}\\
-\varepsilon_{n}\leq-\int_0^{b}|x_{n}'(t)|^{p}dt
 -\int_0^{b}a(t)|x_{n}(t)|^{p}dt+\int_0^{b}u_{n}(t)x_{n}(t)dt\leq
 \varepsilon_{n}, \label{e2.25}\\
-\varepsilon_{n}\leq\int_0^{b}|(x_{n}^{-})'(t)|^{p}dt
 +\int_0^{b}a(t)|x_{n}^{-}(t)|^{p}dt\leq\varepsilon_{n}. \label{e2.26}
\end{gather}
The above inequality tells us that $\{x_{n}^{-}\}$ is bounded in
$W_{\rm per}^{1,p}(T)$ by some constant $M_{1}>0$, and inequalities
\eqref{e2.22}, \eqref{e2.23} jointly give
\begin{equation}
-Mp-\varepsilon_{n}\leq\int_0^{b}(u_{n}(t)x_{n}(t)-pf_{1}(t,x_{n}(t)))dt
\leq Mp+\varepsilon_{n}.
\end{equation}
For the intermediate terms in \eqref{e2.25}, we have
\begin{equation} \label{e2.28}
\begin{aligned}
&\int_0^{b}(u_{n}(t)x_{n}(t)-pf_{1}(t,x_{n}(t)))dt\\
&= \int_{\{x_{n}>0\}}(u_{n}(t)x_{n}(t)-p f(t,x_{n}(t)))dt \\
&\leq \Big(\int_{\{0<x_{n}<M_0\}}+\int_{\{x_{n}\geq
M_0\}}\Big)(u_{n}(t)x_{n}(t)-p f(t,x_{n}(t)))dt,
\end{aligned}
\end{equation}
and
\begin{equation} \label{e2.29}
\int_{\{0<x_{n}< M_0\}}(u_{n}(t)x_{n}(t)-p f(t,x_{n}(t)))dt\leq C
\end{equation}
for some $C>0$. Thus, by the AR-like condition (HF) (iv), we have
\begin{equation} \label{e2.30}
d\int_{\{x_{n}\geq M_0\}}x_{n}^{\alpha}(t)dt\leq
Mp+\varepsilon_{n}+C,
\end{equation}
which yields the boundedness of  $\{x_{n}^{+}\}$ in $L^{\alpha}(T)$.

Next, select $p<\mu<p+\alpha$ in \eqref{e2.4}, then there exists $0<\theta<1$,
such that
\begin{equation}  \label{e2.31}
\frac{1}{\mu}=\frac{1-\theta}{\alpha}+\frac{\theta}{p+\alpha}.
\end{equation}
Thus, using the interpolation inequality, we have
\begin{equation}  \label{e2.32}
\|x_{n}^{+}\|_{\mu}\leq\|x_{n}^{+}\|_{\alpha}^{1-\theta}\|x_{n}^{+}\|_{p+\alpha}^{\theta}\leq
C_{1}\|x_{n}^{+}\|_{p+\alpha}^{\theta}.
\end{equation}
This inequality, together with the boundedness of $\Phi\{x_{n}\}$
and $\{x_{n}^{-}\}$ produce the following estimates
\begin{equation} \label{e2.33}
\begin{aligned}
&\|(x_{n}^{+})'\|_{p}^{p}+\int_0^{b}a(t)|x_{n}^{+}(t)|^{p}dt \\
&\leq Mp+C_{2}M_{1}^{p}+\int_{\{x_{n}>0\}}(a_{1}(t)+c_{1}|x_{n}'(t)|^{\mu})dt \\
&\leq Mp+C_{2}M_{1}^{p}+b\|a_{1}\|_{\infty}+c_{1}\|x_{n}^{+}\|_{\mu}^{\mu}\\
&\leq C_{3}+C_{4}\|x_{n}^{+}\|_{p+\alpha}^{\theta\mu} \\
&\leq C_{3}+C_{5}\|x_{n}^{+}\|_{\infty}^{\theta\mu}\leq
C_{3}+C_{6}\|x_{n}^{+}\|_{1,p}^{\theta\mu}.
\end{aligned}
\end{equation}
Finally, using \eqref{e2.7} with $\varepsilon=0$, we have
\begin{equation}
\|x_{n}^{+}\|_{1,p}^{p}\leq
C_{7}+C_{8}\|x_{n}^{+}\|_{1,p}^{\theta\mu}.
\end{equation}
All the constants $C_{i}>0$, $i=1,2,\cdots,8$ in estimates \eqref{e2.30},
\eqref{e2.31} and \eqref{e2.32}  are independent of $n$.

It is easily to show that $\theta\mu<p$ (please see \cite{kop} for a
reference), therefore the consequence  $\{x_{n}^{+}\}$, and hence
$\{x_{n}\}$ is bounded in $W_{\rm per}^{1,p}(T)$.

On account of the compact embedding of $W_{\rm per}^{1,p}(T)$ into
$C(T)$, we deduce that $\{x_{n}\}$ has a subsequence, which is still
denoted by $\{x_{n}\}$, converging to $x$ weakly in
$W_{\rm per}^{1,p}(T)$, and strongly in $C(T)$. Consequently,
\begin{gather*}
\int_0^{b}a(t)(J_{p}(x_{n}(t)),x_{n}(t)-x(t))dt\to0, \\
\int_0^{b}u_{n}(t)(x_{n}(t)-x(t))dt\to0
\end{gather*}
as $n\to\infty$.

Using $x_{n}-x$ as a test function in \eqref{e2.21}, we have
\begin{equation}
\limsup_{n\to\infty}\langle\mathcal{J}_{p}x_{n},x_{n}-x\rangle\leq0,
\end{equation}
which yields $x_{n}\to x$ in $W_{\rm per}^{1,p}(T)$, since
$\mathcal {J}_{p}$ is of class $(S)_{+}$. 
Thus the proof is complete. 
\end{proof}

If we replaced (HF) (iv) by (vii), then we can verify the
nonsmooth C-condition of $\Phi$ more briefly, i.e.

\begin{lemma} \label{lem2.5}
 Under Hypotheses {\rm (HF) (i)--(iii), (v), (vii)} together with 
{\rm (HA)}, the energy functional
$\Phi$ also satisfies the nonsmooth $C$-condition.
\end{lemma}

Proof of this lemma is omitted here, and a similar proof can be find
in Section 4.

\section{Existence of positive solutions in the scalar case}

To derive the existence results for positive solutions of
periodic system \eqref{e1.1}, we need the following three lemmas.

\begin{lemma} \label{lem3.1}
If hypotheses {\rm (HF) (i)--(iii), (v), (HA)} are satisfied,
then there exist $\rho,\beta>0$, such that
\begin{equation}
\inf_{\|x\|_{1,p}=\rho}\Phi(x)\geq\beta.
\end{equation}
\end{lemma}

\begin{proof} 
Based on the inequality \eqref{e2.6}, we know that
\begin{equation}
\Phi(x)\geq\frac{1}{p}\|x'\|_{p}^{p}+\frac{1}{p}\int_0^{b}a(t)|x(t)|^{p}dt
-\int_0^{b}[\frac{\varepsilon}{p}|x(t)|^{p}+C_{\varepsilon,\mu}|x(t)|^{\mu}]dt,
\end{equation}
where $\varepsilon$ can be selected smaller than
$\int_0^{b}a(t)dt$ so that we can invoking \eqref{e2.7} to derive the
lower estimates for functional $\Phi$; that is,
\begin{equation} \label{e3.3}
\begin{aligned}
\Phi(x)
&\geq c_{\varepsilon}\|x\|_{1,p}^{p}-C_{\varepsilon,\mu}\|x\|_{\mu}^{\mu}\\
&\geq c_{\varepsilon}\|x\|_{1,p}^{p}-c_{2}\|x\|_0^{\mu}\\
&\geq c_{\varepsilon}\|x\|_{1,p}^{p}-c_{3}\|x\|_{1,p}^{\mu},
\end{aligned}
\end{equation}
for some constants $c_{i}>0$, $i=2,3$ independent of 
$x\in W_{\rm per}^{1,p}(T)$. Consider that $\mu>p$, so if $\rho>0$ is taken
small enough, we then have $\inf_{\|x\|_{1,p}=\rho}\Phi\geq\beta$
for some $\beta>0$. 
\end{proof}

The proof of the next lemma is much similar to that of Lemma \ref{lem4.5}, and
here it is omitted.

\begin{lemma} \label{lem3.2}
If {\rm (HF)$_1$ (i)--(iii), (iv'), (vi), (HA)} hold,
then for the fixed number $x_0>0$ taken in $(vi)$, there is
an $r>1$ making the value of $\Phi$ at $rx_0$ negative.
\end{lemma}

\begin{remark} \label{rem3.3} \rm
Lemma \ref{lem3.2} remains true if (HF) (iv') is replaced by (iv),
since the latter is weaker than the former.
\end{remark}

\begin{lemma} \label{lem3.4}
Under conditions {\rm (HF) (i)--(iii),  (HA)}, every critical point
$x$ of $\Phi$, which satisfies $0\in\partial\Phi(x)$, lies in the
set
$$
C_{\rm per}^{1}(T)_{+}=\{x\in C^{1}(T):x(0)=x(b),\;x'(0)
=x'(b),\text{ and }x(t)\geq0,\;\forall t\in T\}
$$ 
and solves problem \eqref{e1.1}.
\end{lemma}

\begin{proof} 
A simple computation shows that every critical point $x$ of
$\Phi$ lies in $C_{\rm per}^{1}(T)$ with the equality
\begin{equation}
\mathcal {J}_{p}x+a(\cdot)J_{p}(x(\cdot))-u=0
\end{equation}
holding for some $u\in S_{\partial
f_{1}(\cdot,x(\cdot))}^{q}$ (cf \cite{pp} or \cite{pp2}), which
means
\begin{equation} \label{e3.5}
\langle\mathcal {J}_{p}x,y\rangle+\int_0^{b}
a(t)J_{p}(x(t))y(t)dt-\int_0^{b}u(t)y(t)dt=0
\end{equation}
for all $y\in W_{\rm per}^{1,p}(T)$.

Take $y=x^{-}$ as the test function in \eqref{e3.5}, we have
\begin{equation}
\int_0^{b}|(x^{-})'(t)|^{p}dt+\int_0^{b}a(t)|x^{-}(t)|^{p}dt=0,
\end{equation}
hence $\|(x^{-})'\|_{p}=0$ and
$\int_0^{b}a(t)|x^{-}(t)|^{p}dt=0$, which produce $x^{-}=0$ in
$W_{\rm per}^{1,p}(T)$ since $\int_0^{b}a(t)dt>0$. Therefore,
$x(t)=x^{+}(t)\geq0$ for all $t\in T$.

By reviewing \eqref{e3.5} again, we can also find $u(t)=0$ a.e. on
$\{x=0\}$, this fact combining with $0\in\partial f(t,0)$ a.e. on
$T$, leads to the desired results $u\in S_{\partial
f(\cdot,x(\cdot))}^{q}$ and then $x$ is a solution of problem \eqref{e1.1}.
\end{proof}

\begin{theorem} \label{thm3.5}
Under hypotheses {\rm (HF) (i)--(vi),  (HA)}, the periodic problem
\eqref{e1.1} has a nontrivial and nonnegative solution in
$C_{\rm per}^{1}(T)_{+}$.
\end{theorem}

\begin{proof} 
Lemmas \ref{lem2.4}, \ref{lem3.1}, \ref{lem3.2} together with the ``Mountain Pass
Lemma" of nonsmooth type \cite{kop,cha} yield a critical point
$x$ of the energy functional $\Phi$ with $\Phi(x)=c$, where
\begin{gather*}
c=\inf_{\gamma\in\Gamma}\max_{t\in[0,1]}\Phi(\gamma(t)), \\
\Gamma=\{\gamma\in C([0,1],W_{\rm per}^{1,p}(T)):\gamma(0)=0,\gamma(1)=rx_0\}.
\end{gather*}
Obviously, $c\geq\beta>0$, which infers $x\neq0$. Finally, invoking
Lemma \ref{lem3.4}, we have $x\in C_{\rm per}^{1}(T)_{+}$ and solves problem
\eqref{e1.1}.
\end{proof}

\begin{theorem} \label{thm3.6}
Putting hypotheses {\rm (HF) (i)--(iii), (iv'), (vii), (HA)} together, 
the periodic problem \eqref{e1.1} has a nontrivial and
nonnegative solution which minimizes the energy functional on the
whole place $W_{\rm per}^{1,p}(T)$.
\end{theorem}

\begin{proof} 
For the minimizer, we appeal to the lower estimates of
the energy functional $\Phi$ in view of \eqref{e2.5} and \eqref{e2.7}; i.e,
\begin{align*}
\Phi(x)&\geq p^{-1}\|x'\|_{p}^{p}+p^{-1}\int_0^{b}a(t)|x(t)|^{p}dt
-\varepsilon p^{-1}\|x\|_{p}^{p}-C_{\varepsilon,\mu}b\\
&\geq c_{\varepsilon}\|x\|_{1,p}^{p}-C_{\varepsilon,\mu}b
\end{align*}
if $0<\varepsilon<\int_0^{b}a(t)dt$. Considering that
$W_{\rm per}^{1,p}(T)$ is reflexive, and $\Phi$ is weakly
semicontinuous, we can use the least action principle of nonsmooth
type, to find another critical point of $\Phi$, which realizes the
whole minimum of $\Phi$ on $W_{\rm per}^{1,p}(T)$, and solves
 problem \eqref{e1.1} simultaneously.
 \end{proof}

\begin{theorem} \label{thm3.7}
If hypotheses {\rm (HF) (i)--(iii), (iv'), (v)-(vii), (HA)}
are all satisfied, the periodic problem \eqref{e1.1} has at least two
nontrivial and nonnegative solutions in $C_{\rm per}^{1}(T)_{+}$.
\end{theorem}

\begin{proof} 
Using Theorems \ref{thm3.5} and \ref{thm3.6}, we can find two solutions of \eqref{e1.1} 
denoted by $x_{1}$ and $x_{2}$ respectively. By reviewing the proofs 
of the two theorems, we can also find $\Phi(x_{1})=c$ and 
$\Phi(x_{2})=\inf_{x\in W_{\rm per}^{1,p}(T)}$. 
On the other hand, Lemmata \ref{lem3.1} and \ref{lem3.2} tell us that, 
$\inf_{x\in W_{\rm per}^{1,p}(T)}\Phi(x)<0<c$, so the two 
solutions $x_{1}$ and $x_{2}$ are different and nontrivial,
 which completes the theorem. 
\end{proof}

\begin{remark} \label{rem3.8} \rm
According to the strong maximum principle (cf. Vazquez \cite{va}),
if condition \eqref{e2.1} in $H(f)(iii)$ is replaced by
\begin{equation}
|u|\leq\hat{c}_0|x|^{p-1}\ \ (\hat{c}_0>0)\quad
\forall u\in\partial f(t,x),
\end{equation}
then every solution of problem \eqref{e1.1} can be positive everywhere (see
also \cite{hup}).
\end{remark}

\begin{remark} \label{rem3.9} \rm
If $a(t)\equiv0$, and {\rm H(f)(vii)} is replaced by
\begin{itemize}
\item[(vii')] there is a function $\theta\in L^{\infty}(T)$ with
$\theta(t)\leq0$ and $\int_0^{b}\theta(t)dt<0$, satisfying
\begin{equation}
\limsup_{x\to+\infty}\frac{pf(t,x)}{x^{p}}\leq\theta(t)
\end{equation}
uniformly for a.e. $t\in T$,
\end{itemize}
then the same conclusions as in Theorem \ref{thm3.6} and \ref{thm3.7} can be reached.
\end{remark}

\section{Existence results in the common case}

In this section we assume that  $N\geq1$, and the function
$f:T\times\mathbb{R}^{N}\to\mathbb{R}$ satisfies
\begin{itemize}
\item[(Hf1)] \begin{itemize}
\item[(i)]  for each $x\in\mathbb{R}^{N}$, $f(\cdot,x)$
is measurable, and $f(\cdot,0)\in L^{\infty}(T,\mathbb{R}_{+})$,

\item[(ii)] for a.e. $t\in T$, $f(t,\cdot)$ is locally Lipschitz, and
$0\in\partial f(t,0)$,

\item[(iii)] for almost all $t\in T$, all $x\in\mathbb{R}^{N}$ and all
$u\in\partial f(t,x)$, the following inequality
\begin{equation}
|u|\leq a_0(t)+d_0|x|^{p-1}
\end{equation}
holds, where $a_0\in L^{\infty}(T,\mathbb{R}_{+})$, and $d_0>0$,

\item[(iv)] there are constants $M_0\geq1$, $d>0$, such that
\begin{equation}
(u,x)-pf(t,x)\geq d|x|^{\alpha}
\end{equation}
for a.e. $t\in T$, and all $u\in\partial f(t,x)$ with $|x|\geq
M_0$, or weakly

\item[(iv')] for a.e. $t\in T$, 
\begin{equation}
(u,x)-pf(t,x)\geq d
\end{equation}
for all $u\in\partial f(t,x)$ with $|x|\geq M_0$,

\item[(v)] $\limsup_{|x|\to0}|x|^{-p}f(t,x)\leq0$, uniformly for
a.e $t\in T$, and


\item[(vi)] there exists $x_0\in\mathbb{R}^{N}$ with $|x_0|>M_0$ satisfying
\begin{equation} \label{e4.4}
\frac{7}{6}\int_0^{b}f(t,x_0)dt>\frac{1}{p}\|a\|_{1}|x_0|^{p},
\end{equation}

\item[(vii)]
\begin{equation}  \label{e4.5}
|x_0|^{p}\sup_{|x|\geq|x|_0}\frac{f(t,x)}{|x|^{p}}\leq\frac{3}{4}(\frac{d}{p+2}+\frac{1}{b}\int_0^{b}f(t,x_0)dt)
\end{equation}
uniformly for a.e. $t\in T$.
\end{itemize}
\end{itemize}

\begin{remark} \label{rem4.1} \rm
Similar to Section 2, we can derive some inequalities from
(Hf1) (iii) (v) and (vii):
There exists $a_{1}\in L^{\infty}(T,\mathbb{R}_{+})$, and $c_{1}>0$, such that
\begin{equation}
|f(t,x)|\leq a_{1}(t)+c_{1}|x|^{\mu}
\end{equation}
for a.e. $t\in T$, and all $x\in\mathbb{R}^{N}$, and
\begin{gather}
f(t,x)\leq\frac{\varepsilon}{p}|x|^{p}+C_{\varepsilon,\mu}, \\
f(t,x)\leq\frac{\varepsilon}{p}|x|^{p}+C_{\varepsilon,\mu}|x|^{\mu}
\end{gather}
for a.e. $t\in T$ and all $x\in\mathbb{R}^{N}$, where $\mu\geq p$ and
$C_{\varepsilon,\mu}>0$ depend on $\varepsilon$ and $\mu$.
\end{remark}

Associated with the periodic system \eqref{e1.1} with $N\geq1$, the energy
functional on $W_{\rm per}^{1,p}(T)$ is defined by
\begin{equation} \label{e4.9}
\Phi(x)=\frac{1}{p}\int_0^{b}|x'(t)|^{p}dt+\frac{1}{p}
\int_0^{b}a(t)|x(t)|^{p}dt-\int_0^{b}f(t,x(t))dt.
\end{equation}
Its Clarke subdifferential can be represented by
\begin{equation}
\partial\Phi(x)=\mathcal
{J}_{p}x+a(\cdot)J_{p}(x(\cdot))-\partial\phi_0(x),
\end{equation}
where $\phi_0(x)=\int_0^{b}f(t,x(t))dt$ for all
 $x\in W_{\rm per}^{1,p}(T,\mathbb{R}^{N})$.
Following the steps as in Sections 2  and 3, we have the following results.

\begin{lemma} \label{lem4.2}
 Under conditions {\rm (Hf1) (i)--(iv), (HA)}, 
the energy functional $\Phi$ satisfies the nonsmooth $C$-condition.
\end{lemma}

\begin{lemma} \label{lem4.3} 
Under conditions {\rm (Hf1) (i)--(iii), (v), (vii), (HA)}, 
the energy functional $\Phi$ satisfies the nonsmooth $C$-condition.
\end{lemma}

\begin{proof} 
Suppose $\{\Phi(x_{n})\}$ is bounded and
$(1+\|x_{n}\|_{1,p})m(x_{n})\to0$ as $n\to\infty$,
we need only to show the boundedness of $\{x_{n}\}$ in
$W_{\rm per}^{1,p}(T,\mathbb{R}^{N})$. Similar to Lemma \ref{lem2.4}, 
there is a constant $M>0$, and functions $u_{n}\in S_{\partial
f(\cdot,x_{n}(\cdot))}^{q}$, $\xi_{n}=\mathcal
{J}_{p}x_{n}+a(\cdot)J_{p}(x_{n}(\cdot))-u_{n}$ 
with $\|\xi_{n}\|_{*}=m(x_{n})$, such that
\begin{gather}
-M\leq p^{-1}\int_0^{b}|x_{n}'(t)|^{p}dt
+p^{-1}\int_0^{b}a(t)|x_{n}(t)|^{p}dt-\int_0^{b}f(t,x_{n}(t))dt
\leq M,  \label{e4.11}\\
|\langle\xi_{n},x_{n}\rangle|
\leq(1+\|x_{n}\|_{1,p})m(x_{n})\to 0 \label{e4.12}
\end{gather}
Inequality \eqref{e4.9} together with \eqref{e4.5} and \eqref{e2.7}, lead to the
 estimate
\begin{align*}
c_0\|x_{n}\|_{1,p}^{p}
&\leq p^{-1}\|x_{n}'\|_{p}^{p}+p^{-1}\int_0^{b}a(t)|x_{n}(t)|^{p}dt\\
&\leq M+\int_0^{b}f(t,x_{n}(t))dt\\
&\leq M+C_{\varepsilon,\mu}b+\varepsilon p^{-1}\|x_{n}\|_{p}^{p}.
\end{align*}
Select $\varepsilon>0$ smaller than $c_0p$ in \eqref{e4.11}, we can
deduce the boundedness of $\{x_{n}\}$ in
$W_{\rm per}^{1,p}(T,\mathbb{R}^{N})$.
 \end{proof}

\begin{lemma} \label{lem4.4}
There exist $\rho,\beta>0$, such that
\begin{equation}
\inf_{\|x\|_{1,p}=\rho}\Phi(x)\geq\beta.
\end{equation}
provided {\rm (Hf1) (i)--(iii), (vi), (HA)} hold.
\end{lemma}

\begin{lemma} \label{lem4.5}
If {\rm (Hf1) (i)--(iii), (iv'), (vi), (HA)} hold,
then for the fixed point $x_0$ taken in $(vi)$, there is an $r\geq1$ 
making the value of $\Phi$ at $rx_0$ negative.
\end{lemma}

\begin{proof} 
Take any $\nu\geq p+2$, and consider the function
\begin{equation}
\lambda(r)=r^{-\nu}f(t,rx_0),\quad r\geq1.
\end{equation}
It is easy to check that $\lambda$ is locally Lipschitz on
$[1,+\infty)$, and according to the chain rule of multiple
functions, we have
\begin{equation}
\lambda'(r)=-\nu r^{-\nu-1}f(t,rx_0)+r^{-\nu}(u,x_0)
\end{equation}
for some $u\in\partial f(t,rx_0)$ and a.e. $r>1$, thus invoking (iv'), 
we obtain
\begin{equation}
\lambda'(r)\geq-(\nu-p) r^{-\nu-1}f(t,rx_0)+dr^{-\nu-1}
=-(\nu-p) r^{-1}\lambda(r)+dr^{-\nu-1}
\end{equation}
a.e. on $(1,+\infty)$.
Notice that
\begin{equation} \label{e4.19}
r^{-1}\lambda(r)= r^{-2} r^{-\nu+1}f(t,rx_0)\leq d_{1}r^{-2},
\end{equation}
where
\[
d_{1}=\frac{3}{4}(\frac{d}{p+2}+\frac{1}{b}\int_0^{b}f(t,x_0)dt)
\geq\sup_{r\geq1}r^{-\nu+1}f(t,rx_0)
\]
is derived from \eqref{e4.5}. Thus we obtain
\begin{equation} \label{e4.20}
\lambda'(r)\geq-d_{1}(\nu-p)r^{-2}+dr^{-\nu-1}.
\end{equation}
For each $r>1$, integrating both sides of the inequality \eqref{e4.20} 
on $[1,r]$, we have
\begin{equation}
\lambda(r)-\lambda(1)\geq-d_{1}(\nu-p)(1-r^{-1})+\frac{d}{\nu}(1-r^{-\nu})
\end{equation}
which lead to the estimate
\begin{equation}
f(t,rx_0)\geq r^{\nu}f(t,x_0)-d_{1}(\nu-p)(r^{\nu}-r^{\nu-1})
+\frac{d}{\nu}(r^{\nu}-1).
\end{equation}
Consequently, for the energy functional $\Phi$, we have
\begin{equation} \label{e4.23}
\begin{aligned}
\Phi(rx_0)&= \frac{|x_0|^{p}r^{p}}{p}\int_0^{b}a(t)dt-\int_0^{b}f(t,rx_0)dt\\
&\leq \frac{\|a\|_{1}|x_0|^{p}}{p}r^{p}-r^{\nu}\int_0^{b}f(t,x_0)dt
 +bd_{1}(\nu-p)(r^{\nu}-r^{\nu-1})-\frac{bd}{\nu}(r^{\nu}-1)\\
&= (bd_{1}(\nu-p)-B-A)r^{\nu-2}[r^{2}-\frac{bd_{1}(\nu-p)}{bd_{1}(\nu-p)-B-A}r\\
&\quad +\frac{K}{bd_{1}(\nu-p)-B-A}r^{-\nu+p+2}+\frac{B}{bd_{1}(\nu-p)-B-A}],
\end{aligned}
\end{equation}
where $A=\int_0^{b}f(t,x_0)dt$, $B=bd/\nu$ and $K=p^{-1}\|a\|_{1}|x_0|^{p}$.

If $K<A$, i.e. $p^{-1}\|a\|_{1}|x_0|^{p}<\int_0^{b}f(t,x_0)dt$, 
then $\Phi(x_0)<0$, and the desired result follows.

If $A\leq K<7A/6$, then take $\nu=p+2$, $r=1+s\delta$, $k=K/(A+B)$ and 
$\delta=(B+A)/(bd_{1}(\nu-p)-B-A)$, and the part in the square brackets 
in the last line of \eqref{e4.23} can be represented by:
\begin{align*}
g(s)
&= (1+s\delta)^{2}-(1+\delta)(1+s\delta)
+k\delta+\frac{B\delta}{A+B}\\
&= (k-(1-s)(1+s\delta)+\frac{B}{A+B})\delta,\quad s\in[0,1].
\end{align*}
Thanks to \eqref{e4.19}, we have $\delta=2$, and 
$g(s)=2(k-(1-s)(1+2s)+B/(A+B))$ which reaches its minimum at $1/4$. 
On the other hand, from \eqref{e4.5}, we can deduce that $B\geq A/3$, 
thus, on account of \eqref{e4.4}, we finally have
\begin{align*}
g(\frac{1}{4})
&= 2(k-\frac{1}{8}-\frac{A}{A+B})\\
&= 2(K-A-\frac{A+B}{8})/(A+B)\\
&\leq 2(K-\frac{7}{6}A)/(A+B)<0,
\end{align*}
which leads to $\Phi((3/2)x_0)<0$ and the proof has been completed. 
\end{proof}

\begin{remark} \label{rem4.6} \rm
In this proof, the desired inequality $f(t,rx_0)\geq r^{\nu} f(t,x_0)- Cr^{p}$ 
for some $\nu>p$ and $C>0$ could not be obtained, since it contradicts 
hypothesis (HF) (vii). In this sense, hypotheses (HF) (vii), (H0) and 
$\int_0^{b}f(t,x_0)dt>0$ are not compatible.
\end{remark}

\begin{lemma} \label{lem4.7}
Under conditions {\rm (Hf1) (i)--(iii),  (HA)}, every critical
point $x$ of $\Phi$ lies in the set
$$
C_{\rm per}^{1}(T,\mathbb{R}^{N})=\{x\in C^{1}(T,\mathbb{R}^{N}):x(0)
=x(b),\; x'(0)=x'(b)\}
$$ 
and solves problem \eqref{e1.1}.
\end{lemma}

Combining all the results, we have the following statement.

\begin{theorem} \label{thm4.8}
If hypotheses {\rm (Hf1) (i)--(iv),  (HA)} are satisfied, then
the periodic problem \eqref{e1.1} has a nontrivial solution in
$C_{\rm per}^{1}(T,\mathbb{R}^{N})$.
\end{theorem}

\begin{remark} \label{rem4.9} \rm
This theorem is a revised result of \cite[Theorem 1]{pp}.
\end{remark}

\begin{theorem} \label{thm4.10}
Under hypotheses {\rm (Hf1) (i)--(iii), (iv'), (vii), (HA)}, 
the energy functional has a minimizer in
$W_{\rm per}^{1,p}(T,\mathbb{R}^{N})$ which solves problem \eqref{e1.1}.
\end{theorem}

\begin{theorem} \label{thm4.11}
Putting all the hypotheses {\rm (Hf1) (i)--(iii), (iv'), (v)--(vii), (HA)} together, 
the periodic problem \eqref{e1.1} has at
least two nontrivial solutions in $C_{\rm per}^{1}(T,\mathbb{R}^{N})$.
\end{theorem}

\begin{corollary} \label{coro4.12}
If hypotheses {\rm (Hf1) (i)--(iii),
(v), (vii), (HA)} hold, and if in addition, we assume $A>K$ 
(see the proof of Lemma \ref{lem4.5}), then the periodic problem \eqref{e1.1} has at
least two nontrivial solutions in $C_{\rm per}^{1}(T,\mathbb{R}^{N})$.
\end{corollary}

\begin{corollary} \label{coro4.13}
If hypotheses {\rm (Hf1) (i)--(iii), (iv'),
(v), (vii), (HA)} hold,  then there exists a $0<\varepsilon_0<7A/6K$, 
such that for each $\varepsilon\in(0,\varepsilon_0)$, the periodic problem
\begin{gather*}
-(J_{p}(x'(t)))'+\varepsilon a(t)J_{p}(x(t))\in\partial f(t,x(t))
\quad \text{a.e. on }T,\\
x(0)=x(b),\quad x'(0)=x'(b),\quad 1<p<+\infty
\end{gather*}
has at least two nontrivial solutions in $C_{\rm per}^{1}(T,\mathbb{R}^{N})$.
\end{corollary}

There are some interesting topics in the study of differential inclusions
with nonsmooth potentials by means of nonsmooth analysis and critical point theory,
such as multiple solutions of elliptic inclusions with p-Laplacian, 
and infinite many  solutions for  the hemivariational inequalities etc. 
For the related discussions, please
refer to \cite{mmm} and \cite{dm}  with references therein.

\subsection*{Acknowledgments}
This research was supported  by the the National Natural Science Foundation
of China (11271316).

\begin{thebibliography}{99}

\bibitem{cl} F. Clarke;
Optimization and Nonsmooth Analysis,
Wiley- Interscience, New York, 1983.

\bibitem{cha} K.-C. Chang;
 Variational methods for non-differential
functionals and their applications to partial differential
equations, \emph{J. Math. Anal. Appl.}, \textbf{80} (1981), 102-129.

\bibitem{cs} J. Cheng, Y. Shao;
 The positive solutions of boundary value problems for a lass of one-dimensional 
$p$-Laplacian, \emph{Nonlinear Anal.}, \textbf{68} (2008), 883-891.

\bibitem{dm} G. D'Agu\`{\i}, G. Molica Bisci;
 Infinitely many solutions for perturbed hemivariational inequalities,
\emph{Bound. Value Probl.}, 2010, Art. ID 363518, 1-15.

\bibitem{eg} L. C. Evans, R. H. Gariepy;
 Measure theory and Fine Properties of Functions, CRC press, Boca
Raton, 1992.

\bibitem{hup} S. Hu, N. S. Papageorgiou;
Positive periodic and homoclinic solutions for nonlinear differential equations with
nonsmooth potential, \emph{Positivity}, \textbf{10} (2006),
343-363.

\bibitem{kop} N. C. Kourogenis, N. S. Papageorgiou;
 Nonsmooth critical point theory and nonlinear elliptic equations at resonance, 
\emph{Kodai Math. J}, \textbf{23} (2000), 108-135.

\bibitem{lz} G. Li, H. S. Zhou;
Asymptotically linear Dirichlet problem for the $p$-Laplacian, 
\emph{Nonlinear Anal.}, {\textbf{43}} (2001), 1043-1055.

\bibitem{mw} J. Mawhin,  M. Willem;
Critical point theory and Hamiltonian systems, Springer-Verlag, New York, 1989.

\bibitem{mm} M. Marcus, V. J. Mizel;
Every suerposition operator mapping one Sobolev space into another is continuous, 
\emph{J. Functional Analysis.}, \textbf{33} (1979), 217-229.

\bibitem{mmm} S. A. Marano, G. Molica Bisci, D. Motreanu;
 Multiple solutions for a class of elliptic hemivariational inequalities,
\emph{J. Math. Anal. Appl.}, \textbf{337} (2008), 85-97.

\bibitem{ni} H. Nikolaos;
 On Neumann hemivariational inequalities, \emph{Abstract and Applied Anal.},
\textbf{7} (2002), 103-112.

\bibitem{pp} E. H. Papageorgiou, N. S. Papageorgiou;
 Non-linear second-order periodic systems with nonsmooth potential, 
\emph{Proc. indian Acad. Sci.}, \textbf{114} (2003), 269-297.

\bibitem{pp2} N. S. Papageorgiou, F. Papalini;
 Pairs of positive solutions for the periodic scalar $p$-Laplacian, 
\emph{J. Fixed Point Theorey Appl.}, \textbf{5} (2009), 157¨C184.

\bibitem{ra} P. Rabinowitz;
 Minimax methods in critical point theory with applications to differential equations,
 \emph{CBMS Regional Conf. Ser. in Math., Amer. Math. Soc., Providence, R. L.,}
 \textbf{65} (1986).

\bibitem{va} S. Varga;
 Existence and infinitely many solutions for an abstract class of hemivariational 
inequalities,  \emph{J. Inequalities Appl.}, \textbf{2} (2005), 89-105.

\bibitem{zh} Q. Zhang;
 Second order differential inclusions with a series of
convex functions, \emph{Nonlinear Anal.}, \textbf{72} (2010), 4053-4062.

\end{thebibliography}

\end{document}
