数学建模单摆的运动分析.doc

对大幅度单摆运动周期公式的研究
\
摘要
单摆作为经典的力学模型,已被众多学者加以研究。多数关于单摆的研究都以摆角小于5度为限,将单摆运动周期近似拟合为在重力与绳拉力下合力的简谐运动的周期。这阻碍了我们对单摆运动周期的进一步探索。此次我们先用力学原理对单摆运动做一般性分析,再通过数学手法简化运动公式,从而得出一般情况下单摆的周期公式,并使用插值法求解出最终结果。
关键字 大幅度 单摆运动 周期 公式
一问题重述
通常对于小幅度(θ≤5。)的单摆运动周期可以近似拟用简谐运动周期公式求解。现在我们试图探究如何求解大幅度单摆运动的周期,并推导出近似公式。
二问题分析
单摆的摆球在重力,摆线拉力的联合作用下做大幅度摆动(θ≥5。)。其运动轨迹可以通过力学分析得到基本运动公式,并以此推导出周期公式。最后通过数学手段简化得出数学解析式。
三基本假设
1空气对单摆运动的阻力与浮力就是如此之小,以至于可以忽略且并不对问题的研究产生交大影响。
2摆线就是一根柔软且无弹性的轻线。
四符号说明
意义
符号
单位
小球质量
m
kg
摆线长度
L
m
小球半径
r
m
重力加速度
g
m/s2
摆角
θ
rad
时间
t
s
角速度
ω
rad/s
五模型建立与求解

一个质量为 ｍ 的小球由一根轻质的长度为 Ｌ 的刚性细绳悬挂在一个固定的支架上(小球半球远远小于细绳长度),小球在重力的作用下可在垂直平面内来回摆动( 不考虑空气阻力),单摆的受力分析图如下:
(即求解大幅摆角单摆运动周期的解析式)
由牛顿第二定律:
d2θdt2+glsinθ=0 (１)
式(１)就是关于θ(角位移)、ｇ(重力加速度)、l(摆长)的一般普遍公式。若给定初始条件,式(１)的任意精度的数值解就是可以求出来的、当θ≤ 5。式可由sinθ≈θ近似求解。但就是当θ≥5。由于误差增大,不能再由上述近似条件求解。通过数值模拟求解的方法可得。
当单摆的摆动角度＞５°,由于系统的机械能守恒,从能量的观点出发也可以求解单摆周期的精确解,这样就不需要详细讨论式(１)非线性微分方程。这时小球运动的速度需要用θ来表示,选择质点运动的轨迹的最低点为势能零点,初始条件为:
θ|t=0＝+θ0
dθdt|t=0＝０
由此可得:
mgl(1-cosθ0 )=12ml2dθdt2+mgl(1-cosθ) (2)
对(2)式数学分析求解,由(2)式得:
dθdt=∓2gl(cosθ-cosθ0) (3)
上式中 ＋(－)表示质点逆时针(顺时针)摆动。
分离变量得摆动周期:
T=22 lg 0θ0dθ(cosθ -cosθ0 ) (4)
又 cosθ=1-2sin(θ/2)2
设sinφ=sin(θ02)sin(θ 2), k=sin(θ02) 则(4)式改写为:
T=22 lg 0θ0dφ(1-k2sinφ2) (5)
其中|θ|<π,k<1。若给定一个摆动幅度(振幅),利用计算机编