对称非线性弹簧振子的周期特性.doc

对称非线性弹簧振子的周期特性 倪亚贤 ,董慎行 ()苏州大学 物理系 ,江苏 苏州 215006 ( ) 摘要 :通过计算机编程 Quick Basic描绘对称非线性弹簧振子振动的特性曲线 ,使难懂的物理过程变得直观 、 形象 . 对称非线性弹簧振子的振动是一种周期性振动 ,但不是严格的简谐振动 ,其振动周期随非线性系数 、振幅的 变化而偏离简谐振动的周期 . 关键词 :非线性 ;弹簧振子 ;简谐振动 ;振动周期 () 文章编号 :100020712 20030420022203 中图分类号 :O 322 文献标识码 :A 1 引言 (任何一个物理量 物体的位置 、电流强度 、电压 、电 ) 场强度 、磁场强度等在某一定值附近的反复变化都称 为振动 . 简谐振动是最简单 、最基本的振动 , 而弹簧振 子是最常见的简谐振动 . 弹簧振子的恢复力与位移成 正比 ,方向总是指向平衡位置 ,如下式所示 : ( )1 f = - k x 实际上只有当位移不很大时 , 这种线性关系才成 立 . 当位移较大时 , 尽管在弹性形变范围内 , 恢复力与 [ 1 ] 位移之间将呈现非线性关系: 2 3 ( )f = - kx - kx - kx + ? 2 1 2 3 图 1 对称非线性弹簧振子 恢复力与位移的关系 当 k= 0 , k?0 时 , 恢复力在松弛点两边还是对 2 3 称的 . 当 k 0 , 恢复力比线性关系所预期的值小 , 称为 3 , 具有简谐振动的特征 . 而当振幅和非线性项 大致相同 非线性软弹簧 . 当 k 0 , 恢复力比线性关系所预期的 3 3 x 的系数| k′| 逐渐增大 , 恢复力逐渐偏离线性 , 振动将 值大 , 称为非线性硬弹簧 , 在这样的恢复力作用下 , 振 不再是严格的简谐振动 . k′ 0 , 保守力具有渐硬特性 子的运动仍是周期性的 . 下面讨论 k?0 , k= 0 , k? 1 2 3 非线性 . k′ 0 , 保守力具有渐软特性非线性 . 0 的情况 , 非线性弹簧振子的恢复力 3f = - k x - kx′ ( )3 2 对称非线性弹簧振子的振动曲线 它与位移的关系可以通过 Quick Basic 编程 , 演示其动 [ 1 ] 态的绘图过程 , 如图 1 所示. 由 于 非 线 性 微 分 方 程 除 少 数 可 以 求 出 精 确 的 解 从图 1 可以看出 , 对称非线性弹簧振子的恢复力 , 目前还没有适用于各类方程的通用的精确解析法 , 外 确实仍具有对称性 , 振子运动对称地分布在平衡位置 因此大部分只能借助计算机按所需精度求数值解 , 或 [ 1 ] 两侧 . 在振幅不大的情况下 , 与线性弹簧振子的恢复力 用近 似 方 法 求 解 , 如 微 扰 法、渐 近 法 、相 平 面 法 收稿日期 :2001 - 09 - 17 ;修回日期 :2002 - 06 - 10 ( ) 基金项目 :苏州大学青年教师基金资助项目 Q3108221 ( ) 作者简介 :倪亚贤 1978 —,女 ,江苏南通人 ,苏州大学物理系硕士生. [ 2 ] , 比线性弹簧振子的抛物线陡 . 当 k′ 振子能量有所增大 ?. 本文主要采用计算机进行数值计算描点来呈 0 , 软弹簧势能曲线也对称 , 但能量较线性时来得小 . 非线性振动的特性 . 当振幅不大时 , 三条势能曲线差别不大 , 近似重叠 , 因 对称非线性弹簧振子的振动曲线方程为 3 x? + k x + kx′ = 0 ( )m4 此 ,小振幅 时 , 弹 簧 振 子 的 振 动 可 以 近 似 看 作 简 谐 振 是一个二阶非线性微分方程 , 难于求得解析解 . 我们 动 . 以采用欧拉迭代法进行数值求解 , 通过 Quick Basic 4 对称非线性弹簧振子的振动周期与非线性系数和 序绘出振动曲线 . 为简化计算 , 设 m = 1 , k = 1 , 并取 振幅的关系=′ 0 和| k′| = 0 . 000 008 , 对应的振动曲线如图 2 所 . 由非线性振子的振动曲线可以发现 , 非线性系数 k对振动周′期有影响 . 我们可以对其振动方程进行分 ( ) 析 . 由振动方程 4可知 , 令 m = 1 , 方程化为 3 ( )x? + k x + kx′ = 0 7 在这里| k′| ν k , 即非线性项与线性项相比很小时 , 可 以通过近似方法求解 , 尽管所得解不甚精确 , 但也能反 映非线性振动的主要特征 . 求解的方法很 多 , 常 用 的 [ 1 ] 一种是微扰法, 也称逐次近似法 , 这种方法是将恢复 图 2 对称非线性弹簧振子的振动曲线 力中的非线性成