力学与实践, 2021, 43(1): 108-111 DOI: 10.6052/1000-0879-20-151

教育研究

滑块碰撞动力学与圆周率的关联1)

李开玮,2)

广东理工学院工业自动化系,广东肇庆 526100

THE RELATION BETWEEN COLLISION DYNAMICS OF SLIDER AND PI 1)

LI Kaiwei,2)

Guangdong Polytechnic College, Zhaoqing 526100, Guangdong, China

通讯作者: 2)李开玮,讲师,研究方向为大学物理教学。E-mail:likaiwei_2006@163.com

责任编辑: 胡漫

收稿日期: 2020-04-20   网络出版日期: 2021-02-08

基金资助: 1)广东理工学院科技项目资助.  2020GKJZD003

Received: 2020-04-20   Online: 2021-02-08

作者简介 About authors

摘要

滑块碰撞与圆周率$\pi$看似不相关的两个事物却具有某种联系,本文研究了滑块碰撞动力学,利用计算工具MATLAB,发现在理想情况下,滑块的碰撞次数与$\pi$有一定的关联,并在理论上证明了计算结果。

关键词: 滑块碰撞 ; 圆周率 ; 碰撞次数

Abstract

The collision of slider seems not relevant to Pi, but there is a certain connection between the two things really. In this paper, the collision dynamics of slider is studied. By using MATLAB, the connection between the number of collision of slider and Pi is established. The calculation results are proved theoretically.

Keywords: collision of slider ; Pi ; number of collision

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

本文引用格式

李开玮. 滑块碰撞动力学与圆周率的关联1). 力学与实践[J], 2021, 43(1): 108-111 DOI:10.6052/1000-0879-20-151

LI Kaiwei. THE RELATION BETWEEN COLLISION DYNAMICS OF SLIDER AND PI 1). MECHANICS IN ENGINEERING[J], 2021, 43(1): 108-111 DOI:10.6052/1000-0879-20-151

滑块碰撞是牛顿力学中常见的问题,一般通过动量守恒和能量变化的关系去研究碰撞过程;而圆周率$\pi$是几何数学的问题,为圆的周长与直径之比,历史上我国数学家刘徽通过割圆术算出$\pi$的小数点后第5位3.14159,之后数学家们计算出的小数点后位数越来越多[1-2]。$\pi$的数字为:3.14159265358979323846264338327950288419716939937510$\cdots$,滑块碰撞与 这两者看起来毫无联系,但有趣的是,作者在研究碰撞问题时,发现滑块的碰撞次数与$\pi$极其相似,并分析了它们的内在关联。

1 问题来源

图1所示,水平光滑的地面上放置小木块$m$,大木块$M$,左端是固定的墙壁,初始时刻$m$静止,$M$以初速度$v_{0}$向左运动,将与$m$发生碰撞,之后$m$获得速度向左运动,将与墙壁发生碰撞反弹,所有碰撞均没有动能损失,求碰撞次数。

图1

图1   问题示意图


解析: 首先考虑最简单情况,若两滑块质量相等$m=M$,则$M$向左运动第一次与$m$碰撞,根据动量守恒和能量守恒,将发生速度传递,接下来$m$以$v_{0}$向左运动与墙壁第二次碰撞,速度反向向右,大小不变,之后向右运动,与$M$第三次碰撞,$m$停止,$M$以$v_{0}$向右运动。总碰撞次数为3。

接下来讨论滑块质量不相等的情况,设$M/m=k$,以水平向左为正方向,为了描述的需要,称两滑块之间的碰撞为碰撞,称$m$与墙壁的碰撞为反射,设第$n$次碰撞(包括反射)后,$m$与$M$速度分别为$u_{n}$和$v_{n} $,对第一次碰撞,根据动量守恒和能量守恒有

$\begin{eqnarray} \left. {{\begin{array}{l} {mu_{1} +Mv_{1} =Mv_{0} }\\ {\dfrac{1}{2}mu_{1}^{2} +\dfrac{1}{2}Mv_{1}^{2} =\dfrac{1}{2}Mv_{0}^{2} }\\ \end{array} }} \right\} \end{eqnarray} $

解式(1)得

$\begin{eqnarray} \left. {{\begin{array}{l} {u_{1} =\dfrac{2k}{k+1}v_{0} } \\[2mm] {v_{1} =\dfrac{k-1}{k+1}v_{0} } \\ \end{array} }} \right\} \end{eqnarray} $

接下来$m$与墙壁发生第一次反射,速度变为

$\begin{eqnarray} \left. {{\begin{array}{l} {u_{2} =-u_{1} } \\ {v_{2} =v_{1} } \\ \end{array} }} \right\} \end{eqnarray} $

第一次反射后,$m$将向右运动与$M$发生第二次碰撞,同样根据动量守恒和能量守恒

$\begin{eqnarray} \left. {{\begin{array}{*{20}c} {mu_{3} +Mv_{3} =mu_{2} +Mv_{2} } \\ {\dfrac{1}{2}mu_{3}^{2} +\dfrac{1}{2}Mv_{3}^{2} =\dfrac{1}{2}mu_{2}^{2} +\dfrac{1}{2}Mv_{2}^{2} } \\ \end{array} }} \right\} \end{eqnarray} $

解式(4)得

$\begin{eqnarray} \left. {{\begin{array}{*{20}c} {u_{3} =\dfrac{(1-k)u_{2} +2kv_{2} }{k+1}} \\[2mm] {v_{3} =\dfrac{2u_{2} +(k-1)v_{2} }{k+1}} \\ \end{array} }} \right\} \end{eqnarray} $

将$u_{n}$,$v_{n} $写成向量$\left( {u_{n} \;v_{n} }\right)$,式(3)和式(5)可变为矩阵运算形式

$\begin{eqnarray} \left\{ {{\begin{array}{*{20}c} {u_{2} } \\ {v_{2} } \\ \end{array} }} \right\}=\left[ {{\begin{array}{*{20}c} {-1} & 0 \\ 0 & 1 \\ \end{array} }} \right]\left\{ {{\begin{array}{*{20}c} {u_{1} } \\ {v_{1} } \\ \end{array} }} \right\} \end{eqnarray} $
$\begin{eqnarray} \left\{ {{\begin{array}{*{20}c} {u_{3} } \\ {v_{3} } \\ \end{array} }} \right\}=\frac{1}{k+1}\left[ {{\begin{array}{*{20}c} {1-k} & {2k} \\ 2 & {k-1} \\ \end{array} }} \right]\left\{ {{\begin{array}{*{20}c} {u_{2} } \\ {v_{2} } \\ \end{array} }} \right\} \end{eqnarray} $

对比式(2)与式(7),发现式(2)也可以写为类似式(7)的样子

$\begin{eqnarray} \left\{ {{\begin{array}{*{20}c} {u_{1} } \\ {v_{1} } \\ \end{array} }} \right\}=\frac{1}{k+1}\left[ {{\begin{array}{*{20}c} {1-k} & {2k} \\ 2 & {k-1} \\ \end{array} }} \right]\left\{ {{\begin{array}{*{20}c} 0 \\ {v_{0} } \\ \end{array} }} \right\} \end{eqnarray} $

接下来第二次反射,第三次碰撞分别与式(3)和式(4)类似,只需将式(3)和式(4)中速度下标加1,因此结论也与式(6)和式(7)相似,因此每次碰撞和反射后,两滑块速度向量变换矩阵分别为

$\begin{eqnarray} \frac{1}{k+1}\left[ {{\begin{array}{*{20}c} {1-k} & {2k} \\ 2 & {k-1} \\ \end{array} }} \right],\ \ \left[ {{\begin{array}{*{20}c} {-1} & 0 \\ 0 & 1 \\ \end{array} }} \right] \end{eqnarray} $

当两滑块碰撞后速度满足$v_{n} <0$ 且$\left| {u_{n} } \right|<\left| {v_{n} }\right|$时,碰撞将不再发生。

根据式(1)$\sim$式(9),作者编写了MATLAB函数文件,如图2所示,利用该函数文件,计算了不同$k$值下的碰撞次数N_collide,如表1所示。根据表1数据,发现当$k=100^{N}$,N_collide与圆周率$\pi $的小数点后$N$位数字一样。

图2

图2   MATLAB函数文件


表1   碰撞次数随着$k$的变化

新窗口打开| 下载CSV


2 理论证明

表1可以猜测,N_collide $=[(10^N)\cdot Pi]$,$[x]$为小于$x$的最大整数,接下来证明这个猜测。

以水平面与墙壁交点为原点,水平向右为正方向,设初始时刻,$m$与$M$位置为$x_{0}$,$y_{0}$,$t$时刻$m$与$M$位置为$x(t)$,$y(t)$,根据图1

$\begin{eqnarray} 0\leqslant x(t)\leqslant y(t) \end{eqnarray} $

当$x(t)=0$时,$m$与墙壁碰撞;当$x(t)=y(t)$时,$m$与$M$碰撞。如图3所示,点$(x(t),y(t))$即可描述$t$时刻两滑块的位置,随着时间的推移,点$(x(t), y(t))$的移动形成位置相空间的轨迹图。第一次碰撞前,$m$位置不变,$M$向左运动,因此轨迹为竖直向下的直线,直至碰到直线$y=x$,发生碰撞,之后两滑块均向左运动,故$x(t)$,$y(t)$均减小,轨迹变为斜向下,直至$x(t)=0$,$m$与墙壁碰撞反射,反射后$m$速度反向,大小不变,因此在图像中有

$\begin{eqnarray} \varphi =\psi \end{eqnarray} $

图3

图3   两滑块相空间轨迹图


之后,将发生第2次滑块间碰撞,反射$\cdots\cdots$。当第$n$次碰撞后,轨迹直线方向为斜向右上方,且接近与$y=x$平行时,将不再与$y=x$相交,即碰撞结束,轨迹与$y$轴、直线$y=x$交点个数之和为碰撞发生次数。由图3可知

$\begin{eqnarray} \theta_{1} =45^{\circ},\ \ \tan \varphi =u_{1} /v_{1},\ \ \theta_{2}=45^{\circ}+\varphi, \cdots \end{eqnarray} $

根据式(12)可以计算出碰撞次数,但计算有点繁琐,主要是因为式(12)中第二式涉及两滑块速度,而每次碰撞后,滑块速度均发生变化,为了计算的简单,接下来我们对相空间$(x(t),y(t))$ 作一个变形,使式(12)中第二式消失。

将$(x(t),y(t))$变形为$(X(t),Y(t))$,变换关系为

$\begin{eqnarray} \left( {X\ \ Y} \right)=(x\ \ y)\left[ {{\begin{array}{*{20}c} {\sqrt m } & 0 \\ 0 & {\sqrt M } \\ \end{array} }} \right] \end{eqnarray} $

故$X(t)=\sqrt m x(t)$,$Y(t)=\sqrt M y(t)$,$X-Y$ 相空间轨迹如图4所示,$x-y$ 相空间中的直线$y=x$变为$X-Y$中的$Y=\sqrt {{M}/{m}} X$,它与$Y$轴夹角为

$\begin{eqnarray} \tan \alpha =\sqrt {\frac{m}{M}} \end{eqnarray} $

图4

图4   $X-Y$相空间轨迹图


当$m$与墙壁碰撞后,$m$速度反向,故图4中$\beta_{1} =\beta_{2} $,当两滑块碰撞时,如图4(a),$r$为直线$Y=\sqrt {{M}/{m}}x$方向向量$\left( {\sqrt m ,\sqrt M } \right)$,$w$为碰前入射向量$(\sqrt m u,\sqrt M v)$,${w}'$为碰后反射向量$(\sqrt m {u}',\sqrt M {v}')$,根据动量守恒,能量守恒式(4)可以改写为

$\begin{eqnarray} \left. {{\begin{array}{l} r\cdot w=r\cdot {w}'={\rm cons1}\\ |w|=|w'| ={\rm cons2} \\ \end{array} }} \right\} \end{eqnarray} $

由式(15)可得入射向量$w$与$r$夹角同反射向量${w}'$与$r$夹角相等,即

$\begin{eqnarray} \varphi_{1}=\varphi_{2} \end{eqnarray} $

由以上分析可以得到图4中各个角度关系

$\begin{eqnarray} \left.\begin{array}{l} \alpha_{1} =\alpha_{2} =\alpha\\ \beta_{1}=\alpha +\alpha_{2} =2\alpha\\ \alpha_{3}=\alpha +\beta_{2} =3\alpha\\ \qquad \quad \vdots \end{array}\right\} \end{eqnarray} $

由式(17)可知第$n$次碰撞前,轨迹方向与直线$Y=\sqrt {{M}/{m}}x$夹角为$n\alpha $,应满足条件

$\begin{eqnarray} n\alpha <\pi \end{eqnarray} $

因此碰撞次数应为

$\begin{eqnarray} n=\left[ {\frac{\pi }{\alpha }} \right] \end{eqnarray} $

由式(14)和式(19)可得

$\begin{eqnarray} n=\left[ {\frac{\pi }{\arctan \sqrt {{m}/{M}} }} \right]=\left[ {\frac{\pi }{\arctan \sqrt {{1}/{k}} }} \right] \end{eqnarray} $

将$k=100^{N}$ 代入式(20)得

$\begin{eqnarray} n=\left[ {\frac{\pi }{\arctan (10^{-N})}} \right] \end{eqnarray} $

将$\arctan x$在$x=0$ 附近作泰勒展开可得

$\begin{eqnarray} \arctan x=x-\frac{x^{3}}{3}+\frac{x^{5}}{5}-\frac{x^{7}}{7}+\cdots \end{eqnarray} $

$\begin{eqnarray} &&\dfrac{1}{\arctan x}-\dfrac{1}{x}=\dfrac{x-\arctan x}{\arctan x\cdot x}=\\&&\qquad \dfrac{\dfrac{1}{3}x^{3}-\dfrac{1}{5}x^{5}+\dfrac{1}{7}x^{7}-\cdots}{x\left( {x-\dfrac{1}{3}x^{3}+\dfrac{1}{5}x^{5}-\cdots} \right)}= \\&&\qquad x\dfrac{\dfrac{1}{3}-\dfrac{1}{5}x^{2}+\dfrac{1}{7}x^{4}-\cdots}{1\mbox{-}x^{2}\left( {\dfrac{1}{3}-\dfrac{1}{5}x^{2}+\cdots} \right)} \end{eqnarray} $

由式(23)可得

$\begin{eqnarray*} \lim\limits_{x\to 0} \left( {\frac{1}{\arctan x}-\frac{1}{x}} \right)=0,\ \ 0<\frac{1}{\arctan x}-\frac{1}{x}<x \end{eqnarray*}$

因此式(21)变为$n=\left[ {\pi \cdot 10^{N}} \right]$ 。

3 结语

本文通过MATLAB计算了滑块碰撞次数,发现当$k=100^{N}$,碰撞次数为$[\pi \times 10^{N}]$,即碰撞次数与圆周率小数点后数字一样,并利用相空间轨迹法巧妙地证明了该联系。将力学过程用几何图像表示出来有助于清晰理解整个运动过程,使计算变得简化。

参考文献

鞠实儿, 张一杰. 刘徽和祖冲之曾计算圆周率的近似值吗? 中国科技史杂志, 2019(4):389-401

[本文引用: 1]

Ju Shi'er, Zhang Yijie.

Did Liu Hui and Zu Chongzhi compute the approximation of $\pi$

The Chinese Journal for the History of Science and Technology, 2019(4):389-401 (in Chinese)

[本文引用: 1]

Mikami Y. The Development of Mathematics in China and Japan. New York: Chelsea Publishing Company, 1913

[本文引用: 1]

/