论文有关弹簧振动.docVIP

论文题目： 弹簧振动方程及求解 作者： 王宝山（数学与应用数学·哈尔滨·150040） 摘要： 列出弹簧振动的方程，并逐步求出一般方程的解 关键词： 简谐运动、阻尼自由振动、阻尼强迫振动 英文摘要： Spring vibration equations listed, and gradually seek the general equation solution. 正文： 1.1一物体B系于挂在顶板上一弹簧的末端，（我们将假设弹簧的质量与这一物体的质量比较起来是可以笑的忽略不计的。）现在来求该物体在外力扰动时的运动微分方程式。当物体B不受外力扰动时，重力被作用与物体B上的弹簧的弹力所平衡而处于静止位置，把物体B静止的位置取位坐标轴x的原点o，向下方向去为正向。 若有一外力f（t）沿竖直方向作用与物体B上，那么物体B将离开静止位置O，记x=x(t)表示物体B在t时刻关于静止位置O的位移，于是, 分别表示物体B的速度和加速度。 由牛顿第二定律F=ma，m是体B的质量，a= 是体B的加速度，而F是作用于物体B上的合外力。 这时合外力F由如下几部分构成。 弹簧的恢复力f1,依胡克定律，弹簧恢复力f1与物体B的位移x成正比，即f1=-cx式中比例常数c（.0）叫作弹性系数，根据所取得坐标系，恢复力f1的方向与位移x的方向相反，所以上式右端添一负号。 空气阻力f2，当速度不太大时，空气阻力可取为与物体B的速度读成正比，亦即f2=-式中比例常数（.0）叫作阻尼系数，式中右边的负号是由于阻力f2的方向与物体B的速度dx/dt的方向相反。 外力f(t) 因此，我们得到 F=-cx-+ f(t) 从而我们得物体B在外力f(t)作用下的运动微分方程式 m ++cx=f(t) (4.1) 如果f(t)≡0，这时得到方程 m ++cx=0 （4.1）’ 而称弹簧的振动为阻尼自由振动。 如果f(t)≡0且=0，即假定没有外力且忽略阻力，这时得到方程 m +cx=0 （4.1）’’ 下面我们求解方程（4.1），（4.1）’，（4.1）’’，并阐明在各个情况下界的物理意义。 1.2.1简谐振动─无阻尼自由振动 令k2=，方程（4.1）’’变为 + k2=0 这是一个二阶常系数齐次方程。特征方程为2+k2=0，特征根是1,2=±ik，它的通解为 x=C1cos kt+C2sin kt （4.45） 其中C1，C2为任意常数。为了阐明上式物理意义，象三角学中常做的那样，我们把上式改写成如下形式： x=Asin(kt+α) (4.46) 由此可见，物体在平衡位置附近做简谐振动。量A称为振幅，幅角kt+α称为振动相位，位相在t=0时所取之值，即α，称为初位相，k=是固有振动频率，T==2为周期。易见，k仅与弹簧的刚度和物体的质量有关。将（4.46）对t微分，可以得物体的运动速度 ν==Akcos(kt+α) 1.2.2如果令= k2,=2n，则方程（4.1）’就变为 +2n+ k2=0 (4.47) 的形式。它是一个二阶常系数线性齐次方程。它的特征方程是2+2n+ k2 =0,特征根是1,2 =-n± (4.48) 现在分三种情况讨论： （1）n2- k20,这时对应于介质阻尼不太大的情形。如果令n2- k2=k12 则（4.48）为1,2 =-n±ik1的形式。这时，方程（4.47）的通解为 x=e-nt(C1cos k1t+C2sin k1t ) 用类似（4.46）的方法将它化为 x=A e-nt sin(k1t+α) (4.49) （2）n2- k2=0，这时通解为 x= e-nt(C1+C2t) (4.50) 此时运动不具振动性质，且当t→+时，x→0. (3) n2- k20，这时对应于介质阻尼相对较大的情形。令n2- k2=h2特征根为 1,2 =-n±h 因为hn，故这时两个特征根为负，通解为x=C1e-(n+h)t+C2e-(n-h)t 易见,此时运动不是周期运动，因而不具振动性质，且当t→+时，x→0. 1.2.3阻尼强迫振动 设作用与物体的外力为 f(t)=Qsin pt其中p，Q均为常量。这时，方程（4.1）具形式 +2n+ k2x=qsin pt (4.51) 其中2n=,k2=c/m,q=. 这是一个二阶常系数线