Processing math: 0%

数学物理学报, 2023, 43(5): 1427-1439

一类双质子弱耦合碰撞系统的对称周期解

王梓欢,, 王超,*

盐城师范学院数学与统计学院 江苏盐城 224001

Symmetric and Periodic Solutions for a Class of Weakly Coupled Systems Composed of Two Particles with Obstacles

Wang Zihuan,, Wang Chao,*

School of Mathematics and Statistics, Yancheng Teachers University, Jiangsu Yancheng 224001

通讯作者: * 王超,Email: wangchaosudamath@163.com

收稿日期: 2022-01-29   修回日期: 2023-02-15  

基金资助: 国家自然科学基金(12071410)

Received: 2022-01-29   Revised: 2023-02-15  

Fund supported: NSFC(12071410)

作者简介 About authors

王梓欢,Email:2539778968@qq.com

摘要

考虑一类具有两个自由度的弱耦合对称碰撞方程的对称碰撞周期解的存在性、重性问题. 在一类关于时间映射的超线性条件下证明了方程无穷多个对称碰撞调和解和对称碰撞次调和解的存在性. 同时, 还给出了一个适合两个自由度的对称碰撞方程的对称碰撞周期解存在的充分条件.

关键词: 对称碰撞方程; 权函数; 对称碰撞周期解; 时间映射; Poincareˊ映射

Abstract

The problems of the existence and multiplicity of symmetric periodic solutions with impact for a class of weakly coupled systems of two degrees of freedom with obstacles are concerned. Under some superlinear assumption on time-mapping, the existence of infinite symmetric harmonic solutions and symmetric subharmonic solutions with impacts of the equation are proved. Furthermore, a sufficient condition for the existence of even and periodic bouncing solution is given for the coupled symmetric impact equations of two degrees of freedom.

Keywords: Symmetric equations with obstacles; Weight functions; Symmetric periodic solutions with impacts; Time-mapping; Poincaré mapping

PDF (686KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

王梓欢, 王超. 一类双质子弱耦合碰撞系统的对称周期解[J]. 数学物理学报, 2023, 43(5): 1427-1439

Wang Zihuan, Wang Chao. Symmetric and Periodic Solutions for a Class of Weakly Coupled Systems Composed of Two Particles with Obstacles[J]. Acta Mathematica Scientia, 2023, 43(5): 1427-1439

1 引言

本文考虑具有完全弹性碰撞的弱耦合对称碰撞系统

\begin{equation} \left\{\begin{array}{lll} & x_i''+q_i(t)g_i(x_i)=p_i(t,x), \hbox{当}\ x_i(t)>0,\\ &x_i(t)\geq0,\\ &\hbox{若}\ x_i(t_0)=0, \hbox{则}\ x_i'(t_0+)=-x_i'(t_0-), \end{array}\right. i=1,2, \end{equation}
(1.1)

对称碰撞周期解的存在性和多解性问题, 其中 x=(x_1,x_2)\in(\Bbb{R}^+)^2,q_i(t):[2\pi]\rightarrow \Bbb R^+(i=1,2) 为连续的 2\pi-周期有界变差偶函数, 且满足

(q_0) (1) \forall t\in[\pi], q_i(t)\geq0 \{t\in[{\rm \pi}]:q_i(t)>0\}\neq\emptyset;

(2) 集合 \{t\in[\pi]:q_i(t)\equiv0\} 只有有限个连通分支;

(3) 若在 (c,d)\subset[\pi]q_i>0q_i(c)=0 (或 q_i(d)=0), 则 q_ic 的一个右邻域内单调不减 (或相应地 q_id 的一个左邻域内单调不增).

p_i(t,x):\Bbb R\times{(\Bbb {R}^+)^2}\rightarrow \Bbb R\ (i=1,2) 是连续的且关于 t2\pi-周期的偶函数; 当 x_i=0 时有 p_i(t,x)=0.

g_i(x_i):\Bbb {R}^+\rightarrow \Bbb R\ (i=1,2) 是局部 Lipschitz 连续函数, 满足

(g_0)\lim\limits_{x_i\to+\infty}g_i(x_i)=+\infty

(\tau_0^+)\lim\limits_{e\to+\infty}\tau_{g_i}^+(e)=0,

其中 \Bbb{R}=(-\infty,+\infty), \Bbb{R}^+=[0,+\infty),\tau_{g_i}^+(e)=\sqrt{2}\int_{0}^{G_{i,+}^{-1}(e)} \frac{1}{\sqrt{e-G_i(s)}}{\rm d}s,G_i(x_i)=\int_0^{x_i} g_i(s){\rm d}s 为方程

x_i''+g_i(x_i)=0

的位势函数.

定义1.1 连续函数 x:\Bbb R\rightarrow ({\Bbb R^+})^2 称为方程 (1.1) 的碰撞解, 如果对 i=1,2 都有

|x_i(t)|+|x'_i(t)|>0, \quad \forall t\in\Bbb{R},

且存在一组双向无穷序列 \{t^j_i\}_j\subset\Bbb R 使得对所有的 j\in \Bbb Z, 以下性质成立

(1) x_i(t^j_i)=0;

(2) x_i'(t^j_i+)=-x_i'(t^j_i-);

(3) x_i\in C^2((t^j_i,t^{j+1}_i),\Bbb R^+) 是方程 (1.1) 通常意义下的解.

我们称 x(t):=(x_1(t),x_2(t)) 是方程 (1.1) 的一个偶的周期碰撞解, 如果 x(t) 是方程 (1.1) 的一个周期碰撞解且

对任意的 t\in\Bbb{R}

x_i(-t)=x_i(t), \quad i=1,2.

由于碰撞方程具有广泛的应用背景, 有关碰撞方程的实验和数值的研究得到了大量的成果. 在数学理论研究方面, 近年来针对二阶碰撞方程也得到了许多成果. 例如, 针对碰撞周期解问题, 主要的方法有: 通过将辅助方程的周期解取极限的方法来研究(见文献 [1,2,3,4,5,6]); 通过把周期解的存在性转化为碰撞相平面上的后继映射或 \rm Poincar\acute{e} 映射的不动点来进行研究(见文献 [7,8,9]); 运用变分法来证明碰撞周期解的存在性(见文献 [10,11,12,13])等. 此外, 还有拟周期解和有界解(见文献 [14,15,16,17,18])、混沌动力行为(见文献 [19,20])等一些结果.

在上述关于碰撞问题的数学理论研究中, 研究对象主要是二阶碰撞方程, 而对于高维碰撞系统的研究没有看到发表的成果.

常微分方程

x''+f(t,x)=0,\ x\in\Bbb{R}^n

对称周期解问题由于其在天体力学和物理上的应用背景而受人关注. 当 n=1 时, 在超线性条件

\lim\limits_{|x|\to+\infty}\frac{g(x)}{x}=+\infty

下, Nakajima[21] 中证明了方程

x''+g(x)=e(t)

无穷多个偶的调和解、任意阶的偶次调和解的存在性和偶次调和解的稠密性分布结果. 在文献 [22] 中, 针对一个自由度的碰撞方程

\left\{\begin{array}{lll} & x''+q(t)g(x)=0, \hbox{当}\ x(t)>0,\\ &x(t)\geq0,\\ &\hbox{若}\ x(t_0)=0, \hbox{则}\ x'(t_0+)=-x'(t_0-),\end{array}\right.

在条件 (q_0)(\tau_0^+) 下得到了与文献 [21] 中类似的结果. 当权函数 q(t) 变号时, 文献 [19] 在超线性条件 (\tau_0^+) 下得到了无穷多个对称碰撞调和解和次调和解的存在性. 一个有趣的问题是, 以往关于二阶方程的对称周期解的相关结论能不能推广到高维的系统中去, 特别地, 对于即使只有两个自由度的弱耦合碰撞方程 (1.1), 有没有对称的碰撞周期解存在. 事实上, 由于没有现成的研究框架, 研究高维碰撞系统碰撞周期解的存在性和多解性不是一件容易的工作. 例如, 由于碰撞周期解空间的不完备性使得拓扑度方法无法直接使用; 非光滑临界点理论往往只能获得碰撞周期解的存在性而无法获得多解性结果, 而且变分法不能够了解系统中的每个质子的具体碰撞信息.

本文主要研究在 (q_0) 和超线性条件 (\tau_0^+) 下方程 (1.1) 对称碰撞周期解的存在性和多解性问题. 在下文中我们是这样安排的: 在第 2 节, 我们首先引进适当的坐标变换, 把碰撞相平面转化为全平面, 通过研究等价系统对应周期解的存在性来得到原系统的碰撞周期解. 其次, 给出对称周期解存在的充分条件. 在第 3 节, 运用相平面分析的方法, 验证了等价系统的 Poincaré 映射在相空间的每一个分平面内的适当大的圆环上具有扭转性. 从而, 在第 4 节, 通过运用一个二维平面上的拉伸扭转零点定理证明了对称周期解存在的充分条件是适用的, 即证明了方程 (1.1) 无穷多个偶的调和解和次调和解的存在性.

x_i'=y_i, i=1,2, 则方程 (1.1) 等价于碰撞系统

\begin{equation} \left\{\begin{array}{lll} & x_i'=y_i, y_i'=-q_i(t)g_i(x_i)+p_i(t,x), \hbox{当}\ x_i(t)>0,\\ &x_i(t)\geq0,\\ &\hbox{若}\ x_i(t_0)=0, \hbox{则}\ x_i'(t_0+)=-x_i'(t_0-),\end{array}\right. i=1,2, \end{equation}
(1.2)

易见, 碰撞系统 (1.2) 的解定义在碰撞相空间 \{(x,y): (x,y)\in \Bbb{R}^{4},\ x_i\geq 0, i=1,2\} 上.

定义1.2 称一个满足 x'(0):=(x_1'(0),x_2'(0))=(0,0) 的偶的碰撞周期解 x(t) 的初值点 (x(0),x'(0)) 为一个碰撞 \varepsilon-点, 其中 x(0):=(x_1(0),x_2(0)). 特别地, 若一个初值点为碰撞 \varepsilon-点的偶的碰撞周期解 x(t) 的最小正周期为 2m\pi(m\in \Bbb{Z}^+), 称此 \varepsilon-点为 m 阶的.

注1.1 在下文中, 为了简洁起见, 通常把 (x_1'(t),x_2'(t))=(0,0) 记为 x'(t)=0, 以及 (x(t),x'(t)):=(x_1(t),x_1'(t),x_2(t),x_2'(t)), t\in\Bbb{R}.

本文中, 总假设 \Bbb{R}^+=[0,+\infty), \Bbb{R}_0^+=(0,+\infty),\Bbb{Z}_0^+:=\{m\in\Bbb{Z}:m>0\}.

2 坐标变换和 \varepsilon-点存在的充分条件

对任意的 k\in \Bbb{Z},\mathcal{I}_k^1:=(4k\pi-\frac{3}{2}\pi, 4k\pi+\frac{\pi}{2}),\mathcal{I}_k^2:=(4k\pi+\frac{\pi}{2}, 4k\pi+\frac{5}{2}\pi).

i=1,2,\theta_i\in \mathcal{I}_k^1r_i>0 时, 定义 \phi_{i,1}: (r_i,\theta_i)\rightarrow (x_i=x_{i,1}(r_i,\theta_i),y_i=y_{i,1}(r_i,\theta_i)), 其中

\begin{equation} \left\{ \begin{array}{lll} x_i=r_i\cos(\frac{\theta_i}{2}+\frac{\pi}{4}),\\ y_i=r_i\sin(\frac{\theta_i}{2}+\frac{\pi}{4}); \end{array} \right. \end{equation}
(2.1)

\theta_i\in \mathcal{I}_k^2r_i>0 时, 定义 \phi_{i,2}: (r_i,\theta_i)\rightarrow (x_i=x_{i,2}(r_i,\theta_i),y_i=y_{i,2}(r_i,\theta_i)), 其中

\begin{equation} \left\{ \begin{array}{lll} x_i=r_i\cos(\frac{\theta_i}{ 2}-\frac{ 3}{4}\pi),\\ y_i=r_i\sin(\frac{\theta_i}{ 2}-\frac{ 3}{ 4}\pi). \end{array} \right. \end{equation}
(2.2)

易见, 对每一个固定的 k\in \Bbb{Z}, 映射 \phi_{i,1}\phi_{i,2} 分别是从 \{(r_i,\theta_i): r_i>0, 4k\pi-\frac{3}{2}\pi<\theta_i<4k\pi+\frac{\pi}{2}\}\{(r_i,\theta_i): r_i>0, 4k\pi+\frac{\pi}{2}<\theta_i<4k\pi+\frac{5}{2}\pi\} 到右半平面 \{(x_i,y_i):(x_i,y_i)\in \Bbb{R}^2, x_i>0\} 上连续的一一映射.

易证, 对 i=1,2 及任意的点列 \{(r_i^{(n)},\theta_i^{(n)})\}_n (r_i^{(n)}>0\hbox{且}\theta_i\in \mathcal{I}_k^1 (或 \theta_i\in \mathcal{I}_k^2)), 若当 n\rightarrow +\infty 时有 r_i^{(n)}\rightarrow r_i^{0}>0\theta_i^{(n)}\rightarrow (4k\pi-\frac{3}{2}\pi)^+(\theta_i^{(n)}\rightarrow (4k\pi+\frac{\pi}{2})^+), 则有

\lim\limits_{n\rightarrow +\infty}x_i^{(n)}=0,\quad \lim\limits_{n\rightarrow +\infty}y_i^{(n)}=-r_i^0.

同理, 若当 n\rightarrow +\infty 时有 r_i^{(n)}\rightarrow r_i^0>0\theta_i^{(n)}\rightarrow (4k\pi-\frac{3}{2}\pi)^-(\theta_i^{(n)}\rightarrow (4k\pi+\frac{\pi}{2})^-), 则有

\lim\limits_{n\rightarrow +\infty}x_i^{(n)}=0,\quad \lim\limits_{n\rightarrow +\infty}y_i^{(n)}=r_i^0.

(0,-r_i^0)(0,r_i^0) 可以用来描述系统(1.2)的碰撞解在碰撞时刻的左右极限.

引理2.1 设系统 (1.2) 的解 (x_i(t),y_i(t))(i=1,2)t=t_0 时发生碰撞, 则碰撞对应到 \theta_i(t_0)=2\tilde{k}\pi-\frac{3}{2}\pi, 即发生一次碰撞相当于 \theta_i(t) 取到一次 2\tilde{k}\pi-\frac{3\pi}{2}(\tilde{k}\in {\Bbb Z}).

i=1,2, 根据碰撞解的定义和变换 (2.1)-(2.2) 式, 若系统 (1.2) 的解 (x_i(t),y_i(t))t=t_0 时发生碰撞, 则 \theta_i(t_0)=4k\pi-\frac{3}{2}\pi,\ r(t_0)=\bar {r_i}>0 或者 \theta_i(t_0)=4k\pi+\frac{\pi}{2},\ r(t_0)=\bar {r_i}>0, 无论哪种情况, 都有 \theta_i(t_0)=2\tilde{k}\pi-\frac{3\pi}{2}(\tilde{k}\in {\Bbb Z}). 证毕.

因此, 在坐标变换 (2.1)-(2.2) 式下, 碰撞系统可写成极坐标下的形式

\begin{equation} \left\{ \begin{array}{l} r_i'=r_i\sin(\frac{\theta_i}{2}+\frac{\pi}{4})\cos(\frac{\theta_i}{2}+\frac{\pi}{4}) [1-\frac{ q_i(t)g_i(x_i)}{ x_i}]+p_i(t,x)\sin(\frac{\theta_i}{2}+\frac{\pi}{4}),\\ \theta_i'=-2\sin^2(\frac{\theta_i}{2}+\frac{\pi}{4}) -\frac{ 2[q_i(t)g_i(x_i)-p_i(t,x)]}{ \ r_i}\cos(\frac{\theta_i}{2}+\frac{\pi}{4}), \\ \end{array}\right. \ \ \theta_i\in \mathcal{I}_k^1;\\ \end{equation}
(2.3)
\begin{equation} \left\{ \begin{array}{l} r_i'=r_i\sin(\frac{\theta_i}{2}+\frac{5\pi}{4})\cos(\frac{\theta_i}{2}+\frac{5\pi}{4}) [1-\frac{ q_i(t)g_i(x_i)}{ x_i}]+p_i(t,x)\sin(\frac{\theta_i}{2}+\frac{5\pi}{4}),\\ \theta_i'=-2\sin^2(\frac{\theta_i}{2}+\frac{5\pi}{4}) -\frac{ 2[q_i(t)g_i(x_i)-p_i(t,x)]}{ r_i}\cos(\frac{\theta_i}{2}+\frac{5\pi}{4}), \\ \end{array}\right. \ \ \theta_i\in \mathcal{I}_k^2,\\ \end{equation}
(2.4)

\begin{equation} \left\{ \begin{array}{l} \theta_i(t_0)=2k_0\pi-\frac{3\pi}{2}\\ r_i(t_0)= \lim_{t\rightarrow t_0}r_i(t)\\ \end{array}\right., \ \ \ \mbox{当}\ \lim_{t\rightarrow t_0}\theta_i(t)=2k_0\pi-\frac{3\pi}{2}, k_0\in \Bbb{Z},\\ \end{equation}
(2.5)
i=1,2.

注意到, (2.3) 和 (2.4) 式描述了碰撞系统的经典解的性质, 而 (2.5) 式则描述了碰撞解在碰撞时刻的性质.

定义2.1 定义在区间 I\subset \Bbb{R} 上的连续函数 (r(t),\theta(t))\in\Bbb{R}^{4} 称为系统 (2.3)-(2.4) 和 (2.5) 式的一个连续解, 如果对每一个 i=1,2

1.r_i(t)>0, 对任意的 t\in I;

2. 对任意的 t\in I \ \hbox{且}\ \theta_i\in \mathcal{I}_k^1, k\in\Bbb{Z}, 都有 (r_i(t),\theta_i(t)) 满足 (2.3) 式;

3. 对任意的 t\in I\ \hbox{且}\ \theta_i\in \mathcal{I}_k^2, k\in\Bbb{Z}, 都有 (r(t),\theta(t)) 满足 (2.4) 式;

4. 对所有的 t_0\in I\lim\limits_{t\rightarrow t_0}\theta_i(t)=4k\pi-\frac{3}{2}\pi\lim\limits_{t\rightarrow t_0}\theta_i(t)=4k\pi+\frac{\pi}{2}, (2.5) 式成立.

注2.1 假设 (r(t),\theta(t)) 为系统 (2.3)-(2.4) 和 (2.5) 式的定义在区间 I\subset \Bbb{R} 上的一个连续解, 则在坐标变换 (2.1)-(2.2) 式下, (x(t),y(t)) 为系统 (1.2) 的定义在区间 I\subset \Bbb{R} 上的一个解. 特别地, 当 I=\Bbb{R} 时, (x(t),y(t)) 为碰撞系统 (1.2) 的满足定义 1.1 的一个碰撞解.

i=1,2, 假设 t_0 为解 (x_i(t),y_i(t)) 的碰撞时刻, 则当 t\neq t_0 时, (2.3)\mbox{、}(2.4) 式右端关于 r_i\mbox{、}\theta_i 是连续的. 当 t\rightarrow t_0 时, \theta'_-(t_0)=\theta'_+(t_0)=-2, 由导数存在性定理知, \theta'(t_0)=-2, 因此振子不会在零点停留.

注2.2 碰撞系统 (1.2) 定义在区间 I\subset\Bbb{R} 上的解 (x_1(t),x_2(t),y_1(t),y_2(t)) 对应到系统 (2.3)-(2.5) 定义在区间 I 上的一个连续解 (r_1(t),r_2(t),\theta_1(t),\theta_2(t)),(r_i(t),\theta_i(t)) 是由 (2.3) 式在区间 I\cap \mathcal{I}_k^1 上的连续解和 (2.4) 式在区间 I\cap \mathcal{I}_k^2 上的连续解交替相接而成的 (i=1,2).\lim\limits_{t\rightarrow t_1}\theta_i(t)=4k\pi+\frac{\pi}{2} 时有相同的结论.

引理2.2 在坐标变换 (2.1)-(2.2) 式下, x_i(r_i,\theta_i),y_i(r_i,\theta_i) 关于 \theta_i2\pi-周期的, 即 x_i(r_i,\theta_i+2\pi)=x_i(r_i,\theta_i), y_i(r_i,\theta_i+2\pi)=y_i(r_i,\theta_i).

事实上, 不妨设 \theta_i\in \mathcal{I}_k^1(k\in\Bbb{Z}), 则有 x_i(r_i,\theta_i)=r_i\cos(\frac{\theta_i}{2}+\frac{\pi}{4}). 由于 \theta_i+2\pi\in \mathcal{I}_k^2, 从而 x_i(r_i,\theta_i+2\pi)=r_i\cos(\frac{\theta_i+2\pi}{2}-\frac{3}{4}\pi)= r_i\cos(\frac{\theta_i}{2}+\frac{\pi}{4})=x_i(r_i,\theta_i).\theta_i\in \mathcal{I}_k^2 时, 则有 \theta_i+2\pi\in \mathcal{I}_{k+1}^1, 易证 x_i(r_i,\theta_i+2\pi)=x_i(r_i,\theta_i). 同理可证 y_i 关于 \theta 也是 2\pi-周期的. 证毕.

注2.3 引理 2.2 说明, 系统 (2.3)-(2.5) 满足初始条件

(r(0),\theta(0))=(r_{10},\theta_{10},r_{20},\theta_{20})

的解 (r(t;r_{10},\theta_{10},r_{20},\theta_{20}),\theta(t;r_{10},\theta_{10},r_{20},\theta_{20})) 与满足初始条件

(r(0),\theta(0))=(r_{10},\theta_{10}+2k_1\pi,r_{20},\theta_{20}+2k_2\pi)

的解 (r(t;r_{10},\theta_{10}+2k_1\pi,r_{20},\theta_{20}+2k_2\pi), \theta(t;r_{10},\theta_{10}+2k_1\pi,r_{20},\theta_{20}+2k_2\pi)),k_1、k_2\in\Bbb{Z}, 对应到碰撞系统 (1.2) 的同一个解 (x(t;x_{10},y_{10},x_{20},y_{20}),y(t;x_{10},y_{10},x_{20},y_{20})), 其中, 当 \theta_{i0}\in \mathcal{I}_k^1 时, (x_{i0},y_{i0}) 由 (2.1) 式取得, 当 \theta_{i0}\in \mathcal{I}_k^2 时, (x_{i0},y_{i0}) 由 (2.2) 式取得 (i=1,2).

在下文中, 对任一个 t_0\in \Bbb{R}q=(a_1,b_1,a_2,b_2)\in (\Bbb{R}^+\times \Bbb{R})^2,x(t;t_0,a_1,b_1,a_2,b_2) 是方程 (1.1) 的满足初始条件

x(t_0)=a,\quad x'(t_0)=b

的解, 其中 a=(a_1,a_2)\in(\Bbb{R}^+)^2, b=(b_1,b_2)\in\Bbb{R}^2, 满足 |a_i|+|b_i|>0, i=1,2. 特别地, 若 t_0=0, b=(0,0),x(t,a):=x(t;0,(a_1,0,a_2,0)).

引理2.3m 是一个正整数, a=(a_1,a_2)\in \Bbb{R}^+_0\times\Bbb{R}^+_0,x(t,a) 是系统 (1.1) 的在 [-4m\pi, 4m\pi] 上有定义的解, 满足

|(x_i(t),x_i'(t))|>0, \forall t\in [-4m\pi, 4m\pi] (i=1,2).

\begin{equation} x_i'(m\pi,a)=0, i=1,2, \end{equation}
(2.6)

x(t,a) 是系统 (1.1) 的一个偶的 2m\pi-周期碰撞解. 特别地, 碰撞 \varepsilon-(a_1,0,a_2,0)m 阶的当且仅当对每一个 k\in\{1,\cdots, m-1\}, 都存在 j=j(k)\in\{1,2\} 使得 x_j'(k\pi)\neq0.

在证明主要结论之前, 先证明一个命题: 若 u(t) 是方程 (1.1) 的在 [-4m\pi, 4m\pi] 上有定义的解且满足

u_i'(k\pi)=0,\quad i=1,2,

其中 k\in\Bbb{Z}0\leq k\leq m, 则对于每一个 t\in[-2m\pi, 2m\pi]u(2k\pi-t)=u(t). 事实上, 首先, 当 t\in[-2m\pi, 2m\pi] 时, 2k\pi-t\in[-2m\pi,4m\pi], 从而 u(2k\pi-t)[-2m\pi, 2m\pi] 上有定义. 其次, 若令 v(t)=u(2k\pi-t),

v_i(t)=u_i(2k\pi-t),\quad i=1,2,

则对每一个 i\in\{1,2\}, 在满足 \forall t\in I, v_i(t)>0 的区间 I\subset[-2m\pi, 2m\pi] 内, 有

\begin{eqnarray*} v_i''(t)+q_i(t)g_i(v_i(t))&=&u_i''(2k\pi-t)+q_i(2k\pi-t)g_i(u_i(2k\pi-t))\\ &=&p_i(2k\pi-t,u(2k\pi-t))\\ &=&p_i(t,v(t)). \end{eqnarray*}

同时, 对那些使得 v_i(t_0)=0t_0\in I 而言 (此时有 u_i(2k\pi-t_0)=0), 易见 v'_{i+}(t_0)=\lim\limits_{t\rightarrow t_0^+}v_i'(t)=-\lim\limits_{t\rightarrow t_0^+}u_i'(2k\pi-t)=\lim\limits_{t\rightarrow t_0^-}u_i'(2k\pi-t)=-\lim\limits_{t\rightarrow t_0^-}v_i'(t)=-v'_{i-}(t_0).

所以, v(t) 也是方程 (1.1) 定义在 [-2m\pi, 2m\pi] 上的一个解. 又对每一个 i=1,2,t=k\pi 时, 有 v_i(k\pi)=u_i(k\pi),u_i'(k\pi)=0u_i(k\pi)\neq 0. 从而 v_i'(k\pi)=-u_i'(k\pi)=0, 故由解的存在唯一性定理知 v(t)=u(t),u(2k\pi-t)=u(t).

由上述命题, 对 i=1,2, 若取 k=0 时有 x'(0,a)=(0,0), 从而对 \forall t\in [-2m\pi, 2m\pi], x_i(-t,a)=x_i(t,a). 同时, 又因为 x'(m\pi)=(0,0),k=m 且令 t=-t, 则有

x_i(2m\pi-t,a)=x_i(t,a)=x_i(-t,a)=x_i(2m\pi+t,a), \quad \forall t\in[-2m\pi,2m\pi].

从而有 x_i'(-2m\pi)=-x_i'(2m\pi)=x_i'(0)=0, 从而 x_i'(2m\pi)=x_i'(0)=x_i'(-2m\pi). 注意到 x_i(2m\pi)=x_i(0)=x_i(-2m\pi), 由解的存在唯一性及系统的周期性知结论成立. 证毕.

由注 2.3 和引理 2.3, 我们有

引理2.4m 是一个正整数, (r(t;r_{10},\theta_{10},r_{20},\theta_{20}),\theta(t;r_{10},\theta_{10},r_{20},\theta_{20})) 是系统 (2.3)-(2.5) 在 [-4m\pi, 4m\pi] 上有定义的连续解, 其中 \theta_{i0}=-\frac{\pi}{2}(\rm mod\, 2\pi). 如果

\begin{equation} \theta_i(m\pi;r_0,\theta_0)=-\frac{\pi}{2}\ (\rm mod\, 2\pi), \quad i=1,2, \end{equation}
(2.7)

(r_{10},0,r_{20},0) 是方程 (1.2) 的一个碰撞 \varepsilon-点. 特别地, (r_{10},0,r_{20},0)m 阶的碰撞 \varepsilon-点当且仅当对每一个 k\in\{1,\cdots, m-1\}, 都存在 j=j(k)\in\{1,2\} 使得 \theta_j(k\pi;r_0,\theta_0)\neq-\frac{\pi}{2}(\rm mod\, 2\pi).

假设 \mathcal{P} 为碰撞系统 (1.2) 的时间 \pi\rm Poincar\acute{e} 映射, 即

\mathcal{P}(x_{10},y_{10},x_{20},y_{20})=(x(\pi;x_{10},y_{10},x_{20},y_{20}),y(\pi;x_{10},y_{10},x_{20},y_{20})),

其中 (x(t;x_{10},y_{10},x_{20},y_{20}),y(t;x_{10},y_{10},x_{20},y_{20})) 是方程 (1.2) 的定义在 [0,\tilde{b})(\tilde{b}>\pi) 上的连续解, 满足 (x(0;x_{10},y_{10},x_{20},y_{20}),y(0;x_{10},y_{10},x_{20},y_{20}))=(x_{10},y_{10},x_{20},y_{20})\in(\Bbb{R}^+\times\Bbb{R})^2. 相应地, 记

\mathcal{\tilde{P}}: (r_{10},\theta_{10},r_{20},\theta_{20})\rightarrow (r(\pi;r_{10},\theta_{10},r_{20},\theta_{20}),\theta(\pi;r_{10},\theta_{10},r_{20},\theta_{20}))

为系统 (2.3)-(2.5) 的时间 \pi\rm Poincar\acute{e} 映射.

3 相平面分析

类似于文献 [23] 中的证明方法, 并做适当的修改, 可以证明

引理3.1 设条件 (g_0)(q_0) 成立, 则对任意给定的区间 [\hat{a},\hat{b}], 都存在正数 R_i^*=R_i^*(\hat{a},\hat{b}) (i=1,2), 使得对任意给定的 t_0\in[\hat{a},\hat{b}]r_{i0}> R_i^*\ (i=1,2), 系统 (2.3)-(2.5) 以 (r_0,\theta_0)\in (\Bbb{ R}^+\times \Bbb{ R})^2 为初值的连续解 (r_1(t;t_0,r_0,\theta_0),\theta_1(t;t_0,r_0,\theta_0), r_2(t;t_0,r_0,\theta_0),\theta_2(t;t_0,r_0,\theta_0))[\hat{a},\hat{b}] 上可延拓, 其中 (r_0,\theta_0):=(r_{10},\theta_{10},r_{20},\theta_{20}).

注3.1 初值问题解的可延拓性既要求解在有限时间内不会出现 blow-up 现象, 同时还要求不会趋向于原点. 条件 (q_0) 可以确保前一个要求满足, R_i^* 的定义可以确保后一个要求满足.

推论3.1(g_0)(q_0) 成立, 则对任一个 T>0 都存在常数 \alpha_i>R_i^*(-T,T)\beta_i(\alpha_i,T)>\alpha_i(i=1,2) 使得下列结论成立:

(1) \forall t_0\in\Bbb{R},(r(t),\theta(t)) 为方程 (2.3)-(2.5) 的满足 R_i^*(-T,T)<r_i(t_0)\leq \alpha_i 的解, 则对所有的 |t-t_0|\leq Tr_i(t)\leq \beta_i(\alpha_i,T);

(2) 对 \forall t_0\in\Bbb{R},(r(t),\theta(t)) 为方程 (2.3)-(2.5) 的满足 r_i(t_0)\geq \beta_i(\alpha_i,T) 的解, 则对所有的 |t-t_0|\leq Tr_i(t)\geq \alpha_i.

注3.2 在引理 3.1 的条件中, r_{i0}> R_i^* 可以保证在区间 [\hat{a},\hat{b}]r_i(t)>0\ (i=1,2).

注3.3 由推论 3.1, 对任意的 m\in\Bbb{Z}_0^+\alpha_i>R_i^*(-2m\pi,2m\pi)a_i>\beta_i(\alpha_i,2m\pi),(r(t,a),\theta(t,a)) 是方程 (2.3)-(2.5) 的在区间 [-2m\pi,2m\pi] 上有定义且满足 r_i(0)=a_i\theta_i(0)=-\frac{\pi}{2}(\rm mod\,\pi) 的连续解, 则对所有 |t|\leq 2m\pi 都有 r_i(t,a)>\alpha_i>0.

注3.4 可以验证, 对 i=1,2I\subset[\pi], 如果 \forall t\in I 都有 q_i(t)>0, 则当 r_i(t) 充分大时有 \theta_i'(t)<0, 从而解 (r_i(t),\theta_i(t)) 在区间 I 内在 u_i-v_i 分平面上是顺时针旋转的.

下面, 对 r_i>0,\ \theta_i\in\Bbb{R}\ (i=1,2) 定义变量变换

\begin{equation} u_i=r_i\cos \theta_i,\ v_i=r_i\sin \theta_i. \end{equation}
(3.1)

(r(t),\theta(t)) 为 (2.3)-(2.5) 式的一个以 (r(\hat{a}),\theta(\hat{a}))=(r_0,\theta_0)=(r_{10},\theta_{10},r_{20},\theta_{20}) 为初值的连续解.

由引理 3.1 知, 当 r_i(\hat{a})>R_i^*(-2\pi,2\pi) 时, 存在一个正数 M_i^0,\forall t\in[\hat{a},\hat{b}]\subset[\pi],r_i(t)>M_i^0.e=(e_1, e_2), 其中 e_i:=(r_{i,0},\theta_{i,0})\in {\Bbb R}^+\times {\Bbb R}\ (i=1,2), 从而可以定义

{\rm rot}_{[\hat{a},\hat{b}]}e_i:=\int_{\hat{a}}^{\hat{b}} -\theta_i' {\rm d}t =\theta_i(\hat{a})-\theta_i(\hat{b}),\quad i=1,2.

易见 \theta_i(\hat{b})-\theta_i(\hat{a}) 是向量 (u_i(t),v_i(t)) 随时间 t 在区间 [\hat{a},\hat{b}] 上增加时在 u_i-v_i 分平面上所转的角度, 它与 \theta_i(t) 的选择无关.

下面, 我们将估计方程 (2.3)-(2.5) 的一个轨道 (r(t),\theta(t)) 在每一个分平面 u_i-v_i 上绕原点一整圈所需时间. 对 i=1,2, 由于 \theta_i(t) 在碰撞时刻连续, \theta_i'(t) 在非碰撞时刻连续且在碰撞时刻极限存在, 所以由导数极限定理 \theta_i'(t) 在碰撞时刻存在. 为此, 下面我们在涉及到 \theta_i(t) 的估计时就不再区分碰撞时刻和非碰撞时刻.

0<\theta_0<\frac{\pi}{2}, u_i o v_i 相平面分成三个区域: \mathbf{I}:=\{(r_i,\theta_i): \frac{\pi}{2}-\theta_0<\theta_i<\frac{\pi}{2}\},\mathbf{II}:=\{(r_i,\theta_i): -\frac{3}{2}\pi+\theta_0\leq\theta_i\leq\frac{\pi}{2}-\theta_0\}\mathbf{III}:=\{(r,\theta): -\frac{3}{2}\pi\leq\theta_i<-\frac{3}{2}\pi+\theta_0\}. 下面我们计算轨道通过每个区域所需的时间. 在下列引理中, 记

\theta_*=\left\{\begin{array}{ll} \frac{\pi}{4},\ -\frac{3}{2}\pi<\theta(\rm{mod}\ 4\pi)<\frac{\pi}{2},\\ \frac{5}{4}\pi,\ \frac{\pi}{2}<\theta(\rm{mod}\ 4\pi)<\frac{5}{2}\pi.\end{array}\right.

引理3.2 假设条件 (g_0) 成立. 设 \theta_0>0\delta>0 为任意的两个常数且 \theta_0<\delta,(r(t),\theta(t)) 是方程 (2.3)-(2.5) 定义在 [\pi] 上的解. 则对 i=1,2 都存在正数 r_i^0>R_i^*:=R_i^*(0,\pi), 使得对每一个区间 [t_1,t_2]\subset[\pi],t_0\in[t_1,t_2], 只要 r_i(t_0)\geq r_i^0 且对任意的 t\in[t_1,t_2](r_i(t),\theta_i(t))\in\mathbf{I},

t_2-t_1\leq \delta.

\delta>0 充分小, 取 0<\theta_0<\delta. 在区域 \mathbf{I}\sin^2(\frac{\pi}{2}-\frac{\theta_0}{2})\leq \sin^2(\frac{\theta_i}{2}+\theta_*)\leq1, 若取 \theta_0 充分小可使 \sin^2(\frac{\theta_i}{2}+\theta_*) >\frac{3}{4}.

由条件 (g_0) 可知存在 d_i>0, 使得 x_i\geq d_i 时有 g_i(x_i)\geq0. 此时有

-\frac{2q_i(t)g_i(x_i)}{r_i}\cos(\frac{\theta_i}{2}+\theta_*)\leq0.

x_i<d_i 时, 由于 q_i(t)\mbox{、}p_i(t,x) 有界, 取 r_0>R_i^* 适当大, 可使当 t\in[\pi]时有

|-\frac{q_i(t)g_i(x_i)}{r_i}\cos(\frac{\theta_i}{2}+\theta_*)|<\frac{3}{32} \quad \hbox{和}\quad |-\frac{p_i(t,x)}{r_i}\cos(\frac{\theta_i}{2}+\theta_*)|<\frac{3}{32}.

\theta_i'=-2\sin^2(\frac{\theta_i}{2}+\theta_*) -\frac{ 2[q_i(t)g_i(x_i)-p_i(t,x)]}{ \ r_i}\cos(\frac{\theta_i}{2}+\theta_*) \hbox{知},

则有 \theta_i'<-\frac{ 9}{ 8}, 从而

t_2-t_1\leq|\frac{\theta_0} {\frac{5}{4}}|\leq\frac{8}{9}\theta_0<\delta.

证毕.

类似可证 \mathbf{III} 的情况.

类似于文献 [22,引理 3.4] 的证明方法, 可以证明 \mathbf{II} 的情况.

引理3.3g_i:\Bbb{R}^+\rightarrow\Bbb{R} 为局部 Lipschitz 连续函数满足条件 (g_0)(\tau_0^+),q_i(t):[\pi]\rightarrow\Bbb R^+ 为连续有界变差函数且满足 (q_0). 则对 \forall \delta>0\forall \theta_0>0 都存在 r_{i,0}=r_{i,0}(\delta,\theta_0)>R_i^*, 使对每一个 [t_1,t_2]\subset I 和对方程 (2.3)-(2.5) 的每一个解 (r(t),\theta(t)), 满足对所有 t\in[t_1,t_2],r_i(t)>r_{i,0}(r_i(t),\theta_i(t))\in\mathbf{II}, 可得 t_2-t_1\leq\delta.

命题3.1g_i,q_i,p_i\ (i=1,2) 满足引理 3.3 中的条件. 则对每个 i\in\{1,2\} 而言, 任取正数 H_i>0, 都存在一个正数 R_i=R_i(H_i)>R_i^* 使得如果 (r(t),\theta(t)) 为方程 (2.3)-(2.5) 的定义在 [\pi] 上的一个解且满足 r_i(0)\geq R_i, 则有 \theta_i(0)-\theta_i(\pi)\geq 2\pi H_i.

任取一个区间 [a',b']\subset[\pi], 则由向量场 (2.3)-(2.5) 的性质易见

\theta_i(a')-\theta_i(b')>-2\pi,

从而

\theta_i(0)-\theta_i(\pi)=((\theta_i(0)-\theta_i(c))+(\theta_i(c)-\theta_i(d)) +(\theta_i(d)-\theta_i(\pi))> (\theta_i(c)-\theta_i(d))-4\pi.

为了估计 \theta_i(c)-\theta_i(d), 我们考虑 c 后轨线与 v_i-正半轴第一次相交的时刻 t^*, 则我们将计算时间 [t^*,d] 内轨道的圈数, 这只要看轨道与 v_i-正半轴相交的次数就可以了.

而且引理 3.2、3.3 将有助于计算轨道转一圈所需的时间. 事实上, 假设我们固定一个 \delta>0\mbox{、}\theta_{i,0}>0, 由引理 3.3, 存在一个 R_0=R_0(\delta,\theta_{i,0}), 使方程 (2.3)-(2.5) 的每一个解, 只要 (r_i(t),\theta_i(t))\in\mathbf{II}r_i(t)\geq R_0 则解穿过 \mathbf{II} 的时间少于 \delta; 另一方面, 只要 0<\theta_{i,0}<\delta, 则解穿过 \mathbf{I}\mathbf{III} 的时间分别少于 \delta. 由弹性性质, 取 N>0使得只要 r_i(t_0)\geq {\rm max} \{N,R_i^*\}=R_i(H_i), 就有对 \forall t\in[\pi],r_i(t)\geq R_0. 从而转一圈所需时间少于 3\delta. 另一方面轨道从出发到到达 v_i-正半轴之间的一段时间 t^*-c 少于 3\delta, 因为轨道至多穿过区域 \mathbf{I},\mathbf{II}, \mathbf{III} 各一次. 从而, 考虑到方程 (1.2), 如果轨道顺时针方向旋转了至少 H_i+2 圈则一定有

\theta_i(0)-\theta_i(\pi)\geq2\pi H.

所以当

3\delta+3\delta(H+2)\leq d-c,\ \hbox{即}\ \delta\leq\frac{d-c}{3H+9}\ \hbox{时},

解可以在 [c,d] 内旋转至少 H_i+2 圈. 得证.

4 无穷个偶周期碰撞解的存在性

为了证明周期碰撞解的存在性, 我们需要下面的零点定理.

引理4.1 [24,定理1.2] 设 \widehat{ABCD}\subset \Bbb{R}^2 是一个矩形, E\widehat{ABCD} 围成的闭区域, \varphi_1, \varphi_2:E\rightarrow\Bbb{R} 是两个连续函数, 满足

\varphi_1(z_1)\cdot\varphi_1(z_2)<0,\quad \forall z_1\in\widehat{AB}, z_2\in\widehat{CD}

\varphi_2(z_3)\cdot\varphi_2(z_4)<0,\quad \forall z_1\in\widehat{BC}, z_2\in\widehat{DA},

则在 E 中至少存在一个点 p 使得 \varphi_1(p)=\varphi_2(p)=0.

定理4.1 假设 (g_0)、(\tau_0^+)(q_0) 成立, 则方程 (1.1) 有无穷个偶的 2\pi-周期碰撞解 \{x(t,a_1^k,a_2^k)\}_{k=1}^{\infty}, 满足

\lim\limits_{k\rightarrow+\infty}a_i^k=+\infty,\quad i=1,2.

由引理 3.1 和推论 3.1, 对任意的 R_i^0>R_i^*(-4\pi,4\pi)\ (i=1,2), 系统 (2.3)-(2.5) 以 (R_1^0,\theta_{10},R_2^0,\theta_{20})\in (\Bbb{ R}^+\times \Bbb{ R})^2 为初值的连续解 (r(t;0,R_0,\theta_0),\theta(t;0,R_0,\theta_0))[-4\pi,4\pi]上可延拓且满足

r_i(t;0,R_0,\theta_0)>0,\quad \forall t\in [-4\pi,4\pi],

其中 R_0=(R_1^0,R_2^0),\theta_0=(\theta_{10},\theta_{20})\in \Bbb{R}^2. 由命题 3.1, 可取适当大的 R_i^*(-4\pi,4\pi)\ (i=1,2), 使得

\theta_i(\pi;R_0,\theta_0)-\theta_i(0;R_0,\theta_0)<0, \quad i=1,2.

任取 c_1^1>R_1^*(-4\pi,4\pi), c_2^1>R_2^*(-4\pi,4\pi), 则易知, 存在正数 R_i'>c_i^1>R_i''>0\ (i=1,2), 使得

0<R_1''\leq r_1(t;c_1^1,\theta_{10},r_2,\theta_{20})\leq R_1',\quad \forall t\in[\pi],

其中 r_2\geq c_2^1, (\theta_{10},\theta_{20})\in \Bbb{R}^2, 以及

0<R_2''\leq r_2(t;r_1,\theta_{10},c_2^1,\theta_{20})\leq R_2',\quad \forall t\in[\pi],

其中 r_1\geq c_1^1, (\theta_{10},\theta_{20})\in \Bbb{R}^2.

易见, 存在两个正数 M_1、M_2>0, 使得

|\theta_1'(t;c_1^1,\theta_{10},r_2,\theta_{20})|\leq M_1,\quad |\theta_2'(t;r_1,\theta_{10},c_2^1,\theta_{20})|\leq M_2,

其中 (\theta_{10},\theta_{20})\in\Bbb{ R}^2. 从而 \int_0^{\pi}\theta_i'{\rm d}t 为有界值, 这样, 一定存在 j_1^1、j_2^1\in \Bbb{Z}^+, 使得

\begin{equation} -j_1^1\pi<\theta_1(\pi;c_1^1,\theta_{10},r_2,\theta_{20})-\theta_1(0;c_1^1,\theta_{10},r_2,\theta_{20})<0 \end{equation}
(4.1)

对任意的 r_2\geq c_2^1(\theta_{10},\theta_{20})\in \Bbb{R}^2 成立,

\begin{equation} -j_2^1\pi<\theta_2(\pi;r_1,\theta_{10},c_2^1,\theta_{20})-\theta_2(0;r_1,\theta_{10},c_2^1,\theta_{20})<0 \end{equation}
(4.2)

对任意的 r_1\geq c_1^1, (\theta_{10},\theta_{20})\in \Bbb{R}^2 成立.

再由命题 3.1, 可取 d_1^1>R_1'd_2^1>R_2', 使得

\begin{equation} \theta_1(\pi;d_1^1,\theta_{10},r_2,\theta_{20})-\theta_1(0;d_1^1,\theta_{10},r_2,\theta_{20})<-j_1^1\pi \end{equation}
(4.3)

对任意的 r_2\geq c_2^1(\theta_{10},\theta_{20})\in \Bbb{R}^2 成立,

\begin{equation} \theta_2(\pi;r_1,\theta_{10},d_2^1,\theta_{20})-\theta_2(0;r_1,\theta_{10},d_2^1,\theta_{20})<-j_2^1\pi \end{equation}
(4.4)

对任意的 r_1\geq c_1^1, (\theta_{10},\theta_{20})\in \Bbb{R}^2 成立.

现在, 定义连续映射

\begin{array}{lll} f_1: [c_1^1,d_1^1]\times[c_2^1,d_2^1]\rightarrow \Bbb{R}\\ (r_1,r_2)\mapsto \theta_1(\pi;r_1,-\frac{\pi}{2},r_2,-\frac{\pi}{2})-\theta_1(0;r_1,-\frac{\pi}{2},r_2,-\frac{\pi}{2})+j_1^1\pi, \end{array}

\begin{array}{lll} f_2: [c_1^1,d_1^1]\times[c_2^1,d_2^1]\rightarrow \Bbb{R}\\ (r_1,r_2)\mapsto \theta_2(\pi;r_1,-\frac{\pi}{2},r_2,-\frac{\pi}{2})-\theta_2(0;r_1,-\frac{\pi}{2},r_2,-\frac{\pi}{2})+j_2^1\pi. \end{array}

易见, f_1f_2 都是矩形 E:=[c_1^1,d_1^1]\times[c_2^1,d_2^1]\subset\Bbb{R}^2\Bbb{R} 上的连续函数. 由 (4.1) 和 (4.3) 式知

f_1(c_1^1,r_2)>0, \ c_2^1\leq r_2\leq d_2^1,
f_1(d_1^1,r_2)<0, \ c_2^1\leq r_2\leq d_2^1;

以及由 (4.2) 和 (4.4) 式知

f_2(r_1,c_2^1)>0, \ c_1^1\leq r_1\leq d_1^1,
f_2(r_1,d_2^1)<0, \ c_1^1\leq r_1\leq d_1^1.

由引理 4.1 知, 至少存在一个点 (a_1^1,a_2^1)\in E 使得 f_1(a_1^1,a_2^1)=0, f_2(a_1^1,a_2^1)=0,

\theta_1(\pi;a_1^1,-\frac{\pi}{2},a_2^1,-\frac{\pi}{2})-\theta_1(0;a_1^1,-\frac{\pi}{2},a_2^1,-\frac{\pi}{2})=-j_1^1\pi,
\theta_2(\pi;a_1^1,-\frac{\pi}{2},a_2^1,-\frac{\pi}{2})-\theta_2(0;a_1^1,-\frac{\pi}{2},a_2^1,-\frac{\pi}{2})=-j_2^1\pi.

由引理 2.4, x(t;a_1^1,0,a_2^1,0) 是方程 (1.1) 的 2\pi-碰撞周期解. 并且, 由对称性知

\theta_i(2\pi;a_1^1,a_2^1)-\theta_i(0;a_1^1,a_2^1)=-2j_i^1 \pi,\quad i=1,2,

即在一个周期内 x_1(t;a_1^1,0,a_2^1,0) 恰好有 j_1^1 次完全弹性碰撞, x_2(t;a_1^1,0,a_2^1,0) 恰好有 j_2^1 次完全弹性碰撞.

接下来, 选取 c_1^2> d_1^1 以及 c_2^2> d_2^1, 类似地, 存在 j_1^2、j_2^2\in\Bbb{Z}^+(j_1^2> j_1^1、j_2^2> j_2^1), 使得

\begin{equation} -j_1^2\pi<\theta_1(\pi;c_1^2,\theta_{10},r_2,\theta_{20})-\theta_1(0;c_1^2,\theta_{10},r_2,\theta_{20})<0 \end{equation}
(4.5)

对任意的 r_2\geq c_2^2(\theta_{10},\theta_{20})\in \Bbb{R}^2 成立,

\begin{equation} -j_2^2\pi<\theta_2(\pi;r_1,\theta_{10},c_2^2,\theta_{20})-\theta_2(0;r_1,\theta_{10},c_2^2,\theta_{20})<0 \end{equation}
(4.6)

对任意的 r_1\geq c_1^2, (\theta_{10},\theta_{20})\in \Bbb{R}^2 成立.

同理, 存在 d_1^2>c_1^2d_2^2>c_2^2, 使得

\begin{equation} \theta_1(\pi;d_1^2,\theta_{10},r_2,\theta_{20})-\theta_1(0;d_1^2,\theta_{10},r_2,\theta_{20})<-j_1^2\pi \end{equation}
(4.7)

对任意的 r_2\geq c_2^2(\theta_{10},\theta_{20})\in \Bbb{R}^2 成立,

\begin{equation} \theta_2(\pi;r_1,\theta_{10},d_2^2,\theta_{20})-\theta_2(0;r_1,\theta_{10},d_2^2,\theta_{20})<-j_2^2\pi \end{equation}
(4.8)

对任意的 r_1\geq c_1^2, (\theta_{10},\theta_{20})\in \Bbb{R}^2 成立. 从而类似可证, 至少存在一个点 (a_1^2,a_2^2)\in[c_1^2,d_1^2]\times [c_2^2,d_2^2] 使得

\theta_1(\pi;a_1^2,-\frac{\pi}{2},a_2^2,-\frac{\pi}{2})-\theta_1(0;a_1^2,-\frac{\pi}{2},a_2^2,-\frac{\pi}{2})=-j_1^2\pi,
\theta_2(\pi;a_1^2,-\frac{\pi}{2},a_2^2,-\frac{\pi}{2})-\theta_2(0;a_1^2,-\frac{\pi}{2},a_2^2,-\frac{\pi}{2})=-j_2^2\pi.

由引理 2.4, x(t;a_1^2,0,a_2^2,0) 是方程 (1.1) 的 2\pi-碰撞周期解. 并且, 在一个周期内 x_i(t;a_1^2,0,a_2^2,0) 恰好有 j_i^2(i=1,2) 次完全弹性碰撞.

类似讨论, 可以证明存在两个正数列

d_i^2<c_i^3<d_i^3<c_i^4<d_i^4<\cdots (i=1,2)

和向量列 \{(a_1^k,a_2^k)\}_{k=3}^\infty(i=1,2) 以及两个正整数列 j_i^3<j_i^4<\cdots (i=1,2) 满足

\theta_1(\pi;a_1^k,-\frac{\pi}{2},a_2^k,-\frac{\pi}{2})-\theta_1(0;a_1^k, -\frac{\pi}{2},a_2^k,-\frac{\pi}{2})=-j_1^k\pi,
\theta_2(\pi;a_1^k,-\frac{\pi}{2},a_2^k,-\frac{\pi}{2})-\theta_2(0;a_1^k, -\frac{\pi}{2},a_2^k,-\frac{\pi}{2})=-j_2^k\pi,

其中 c_i^k\leq a_i^k\leq d_i^k, k=3,4,\cdots ; i=1,2.

注意到 \lim\limits_{k\rightarrow+\infty}j_i^k=+\infty\ (i=1,2), 由推论 3.1 和命题 3.1 知

\lim\limits_{k\rightarrow+\infty}a_i^k=+\infty,\quad i=1,2.

证毕.

由定理 4.1 的证明过程易得

推论4.1 假设条件 (g_0)、(\tau_0^+)(q_0) 成立, 则存在一个适当大的正整数 j_0 使得对于任意的正整数 j>j_0, 方程 (1.1) 至少存在一个偶的 2\pi-周期碰撞解 (x_1(t),x_2(t)) 满足 x_1(t)x_2(t)[2\pi] 内都恰好具有 j次碰撞. 同时, 存在一列正整数

j_0<j_1<j_2<\cdots <j_k<\cdots,\quad k\in\Bbb{Z}^+,

使得方程 (1.1) 至少存在一个偶的 2\pi-周期碰撞解 (x_1^k(t),x_2^k(t)) 满足 x_1^k(t)x_2^k(t)[2\pi] 内都恰好具有 j_k次碰撞且

\lim\limits_{k\rightarrow+\infty}|x_i^k(t)|+|x_i'^k(t)|=+\infty,\quad i=1,2.

定理4.2 假设 (g_0)、(\tau_0^+)(q_0), 则对任意的 m\in \Bbb{Z}_0^+, 方程 (1.1) 有无穷个偶的 2m\pi-周期碰撞解 \{x(t;a_1^k,a_2^k)\}_{k=1}^{\infty} 满足

\lim\limits_{k\rightarrow+\infty}a_i^k=+\infty\ (i=1,2),

并且, (a_1^k,0,a_2^k,0) 都是 m 阶的碰撞 \varepsilon-点.

类似于定理 4.1 证明中的构造方法, 可以找到四个正数列 \{c_i^k\}_k、\{d_i^k\}_k\ (i=1,2) 和两个素数列 q_i^1<q_i^2<\cdots <q_i^k<\cdots\ (i=1,2), 满足

c_1^{(i)}<d_1^{(i)}<c_2^{(i)}<d_2^{(i)}<\cdots <c_k^{(i)}<d_k^{(i)}<\cdots

以及

\lim\limits_{k\rightarrow+\infty}c_i^k=+\infty\ (i=1,2),

使得

\theta_i(m\pi;c_1^k,-\frac{\pi}{2},r_2,-\frac{\pi}{2})-\theta_i(0;c_1^k,-\frac{\pi}{2},r_2,-\frac{\pi}{2})>-q_1^k\pi,\quad \forall r_2\in[c_2^k,d_2^k],
\theta_i(m\pi;d_1^k,-\frac{\pi}{2},r_2,-\frac{\pi}{2})-\theta_i(0;d_1^k,-\frac{\pi}{2},r_2,-\frac{\pi}{2})<-q_1^k \pi,\quad \forall r_2\in[c_2^k,d_2^k];
\theta_i(m\pi;r_1,-\frac{\pi}{2},c_2^k,-\frac{\pi}{2})-\theta_i(0;r_1,-\frac{\pi}{2},c_2^k,-\frac{\pi}{2})>-q_2^k \pi,\quad \forall r_1\in[c_1^k,d_1^k],
\theta_i(m\pi;r_1,-\frac{\pi}{2},d_2^k,-\frac{\pi}{2})-\theta_i(0;r_1,-\frac{\pi}{2},d_2^k,-\frac{\pi}{2})<-q_2^k \pi,\quad \forall r_1\in[c_1^k,d_1^k].

从而, 由引理 4.1 知, 存在数列 \{a_i^k\}_k,\ c_i^k<a_i^k<d_i^k(k=1,2,\cdots )\ (i=1,2), 使得

\quad \theta_i(m\pi;a_1^k,-\frac{\pi}{2},a_2^k,-\frac{\pi}{2})+\frac{\pi}{2}=-q_i^k \pi,

从而

\theta_i(m\pi,a_1^k,-\frac{\pi}{2},a_2^k,-\frac{\pi}{2})=-\frac{\pi}{2}\ ({\rm mod}\,\pi)\ (i=1,2).

由引理 2.4 知, x(t;a_k^{(1)},a_k^{(2)})(k=1,2,\cdots ) 是方程 (1.1) 的 2m\pi-偶周期碰撞解. 由对称性易证 \theta_i(2m\pi;a_k^{(1)},a_k^{(2)})-\theta_i(0;a_k^{(1)},a_k^{(2)})=-2q_i^k \pi,x_i(t;a_1^k,a_2^k) 在一个周期内恰好发生 q_i^k 次碰撞. 易见, (a_1^k,0,a_2^k,0)m 阶碰撞 \varepsilon-点. 证毕.

参考文献

Lazer A, Mckenna P.

Periodic bounding for a forced linear spring with obstacle

Diff Inte Equa, 1992, 5: 165-172

[本文引用: 1]

Bonheune D, Fabry C.

Periodic motions in impact oscillators with perfectly elastic bouncing

Nonlinearity, 2002, 15: 1281-1298

DOI:10.1088/0951-7715/15/4/314      URL     [本文引用: 1]

Qian D.

Large amplitude periodic bouncing for impact oscillators with damping

Proc Amer Math Soc, 2005, 133: 1797-1804

DOI:10.1090/proc/2005-133-06      URL     [本文引用: 1]

Sun X, Qian D.

Periodic bouncing solutions for attractive singular second-order equations

Nonlinear Anal, 2009, 71: 4751-4757

DOI:10.1016/j.na.2009.03.049      URL     [本文引用: 1]

Liu Q, Wang Z.

Periodic impact behavior of a class of Hamiltonian oscillators with obstacles

J Math Appl Anal, 2010, 65: 67-74

[本文引用: 1]

Fonda A, Sfecci A.

Periodic bouncing solutions for nonlinear impact oscillators

Adv Nonlinear Stu, 2013, 13: 179-189

[本文引用: 1]

Qian D, Torres P.

Bouncing solutions of an equation with attractive singularity

Proc Roy Soc Edingburgh Sect A, 2004, 134: 201-213

[本文引用: 1]

Qian D, Torres P.

Periodic motions of linear impact oscillators via successor map

SIAM J Math Anal, 2005, 36: 1707-1725

DOI:10.1137/S003614100343771X      URL     [本文引用: 1]

丁卫, 钱定边.

碰撞 Hamiltonian 系统的无穷小周期解

中国科学: 数学, 2010, 40: 563-574

DOI:10.1360/za2010-40-6-563      URL     [本文引用: 1]

Ding W, Qian D B.

Infinitesimal periodic solutions of impact Hamiltonian systems

Scientia Sinica Mathematica, 2010, 40: 563-574

DOI:10.1360/za2010-40-6-563      URL     [本文引用: 1]

Jiang M.

Periodic solutions of second order differential equations with an obstacle

Nonlinearity, 2006, 19: 1165-1183

DOI:10.1088/0951-7715/19/5/007      URL     [本文引用: 1]

Ding W.

Subharmonic solutions of sublinear second order systems with impacts

J Math Anal Appl, 2011, 379: 538-548

DOI:10.1016/j.jmaa.2011.01.036      URL     [本文引用: 1]

Ding W, Qian D, Wang C, Wang Z.

Existence of Periodic Solutions of sub-linear Hamiltonian Systems

Acta Math Sin Eng Ser, 2016, 32: 621-632

[本文引用: 1]

Wu X, Li X.

Existence and multiplicity of solutions for a class of forced vibration problems with obstacles

Nonlinear Anal, 2009, 71: 3563-3570

DOI:10.1016/j.na.2009.02.021      URL     [本文引用: 1]

Zharnitsky V.

Invariant tori in Hamiltonian systems with impacts

Comm Math Phys, 2000, 211: 289-302

DOI:10.1007/s002200050813      URL     [本文引用: 1]

Ortega R. Dynamics of a forced oscillator having an obstacle//Benci V, Cerami G, Degiovanni M, Fortunato D. Variational and Topological Methods in the Study of Nonlinear Phenomena. Progre Nonlinear Differential Equations Appl 49. Boston: Birkhauser Press, 2001: 75-87

[本文引用: 1]

Qian D, Sun X.

Invariant tori for asymptotically linear impact oscillators

Science in China Series A, 2006, 49: 669-687

[本文引用: 1]

Wang Z, Liu Q, Qian D.

Existence of quasi-periodic solutions and Littlewood's boundedness problem of sub-linear impact oscillators

Nonlinear Anal: TMA, 2011, 74: 5606-5617

DOI:10.1016/j.na.2011.05.046      URL     [本文引用: 1]

Piao P, Sun X.

Boundedness of solutions for a class of impact oscillators with time-denpendent polynomial potentials

Comm Pure Appl Anal, 2013, 13: 645-655

DOI:10.3934/cpaa      URL     [本文引用: 1]

Wang C, Qian D, Liu Q.

Impact oscillators of Hill's type with indefinite weight: Periodic and chaotic dynamics

Discrete and Continuous Dynamical Systems, 2016, 36: 2305-2328

DOI:10.3934/dcdsa      URL     [本文引用: 2]

Herrera A, Torres P.

Periodic solutions and chaotic dynamics in forced impact oscillators

SIAM J Appl Dyn Syst, 2013, 12: 383-414

DOI:10.1137/120880902      URL     [本文引用: 1]

Nakajima F.

Even and periodic solution of the equation \ddot{u}+g(u)=e(t)

J Differential Equation, 1990, 83: 277-299

DOI:10.1016/0022-0396(90)90059-X      URL     [本文引用: 2]

王超.

一类超线性 Hill 型对称碰撞方程的周期运动

中国科学: 数学, 2014, 44: 235-248

DOI:10.1360/N012013-00093      URL     [本文引用: 2]

Wang C.

The periodic motions of a class of symmetric superlinear Hill's impact equations

Scientia Sinica Mathematica, 2014, 44: 235-248

DOI:10.1360/N012013-00093      URL     [本文引用: 2]

Papini D. Boundary value problems for second order differential equations with super-linear terms: a topological approch. Udline: Udline University, 2000

[本文引用: 1]

王超, 刘期怀, 钱定边, 王志国.

拓扑定理及其在超线性脉冲方程中的应用

中国科学: 数学, 2014, 44: 957-968

DOI:10.1360/012014-42      URL     [本文引用: 1]

Wang C, Liu Q H, Qian D B, Wang Z G.

A topological theorem and its application to the superlinear equations including impulses for the existence of periodic solutions

Scientia Sinica Mathematica, 2014, 44: 957-968

DOI:10.1360/012014-42      URL     [本文引用: 1]

/