matlab仿真弹簧波,基于Matlab/Simulink的三弹簧谐振子微振动的仿真实验

第 29 卷 第 6 期 大 学 物 理 实 验 Vol． 29 No． 6 2016 年 12 月 PHYSICAL EXPEＲIMENT OF COLLEGE Dec． 2016 收稿日期: 2016-07-14 基金项目: 江苏省现代教育技术研究课题(2013-Ｒ-25643) 文章编号:1007-2934(2016)06-0103-05 基于 Matlab/Simulink 的三弹簧谐振子微振动的仿真实验 张 林 (南京林业大学，江苏 南京 210037) 摘 要: 建立了三弹簧振子耦合系统微振动的实验模型，利用 Matlab/Simulink 仿真软件对该实验进行仿真建模，讨论系统在不同实验参数下从周期振动到混沌的各种动力学状态。仿真实验结果表明， 参数的选取对系统微振动的动力学行为有着很大的影响。 关 键 词: 谐振子;微振动;Simulink 仿真 中图分类号: O 4-39 文献标志码: A DOI:10． 14139/j． cnki． cn22-1228． 2016． 006． 028 在物理实验教学中，一维弹簧谐振子的微振动是学生认识客观世界复杂振动的起点和基础。在无外力作用时，单个弹簧振子在平衡位置附近做简单的周期性无阻尼振动，而对于多个弹簧振子组成的耦合系统，其微振动在不同参数下可表现为周期振动、准周期振动以及无规则振动(混沌)，因此多个弹簧振子耦合系统的微振动引起了人们广泛的研究兴趣。文献［ 1］首先使用了Matlab 软件对三自由度的微振动进行了计算和讨论，并模拟了系统的动力学过程。随后文献［2-6］分别使用数值计算和虚拟仿真等方法对多个弹簧振子耦合系统的复杂振动进行了较深入的研究。 为了让学生更好地认识到多弹簧耦合振动的复杂性和可调制性，本文建立了一个三弹簧振子耦合系统的实验模型，并利用 Matlab/Simulink 仿真软件讨论了不同实验参数对系统微振动状态的影响。 1 一维弹簧振子系统微振动实验的simulink 建模 建立如图 1 所示的一维三弹簧振子的实验系统:三个弹簧振子的质量均为 m ，弹簧的劲度系数分别为 k1 、k2 和 k3 ，三个谐振子在微振动过程中相对自身的平衡位置的位移分别为 x1 、x2 和 x3 ，取 x1 ， x2 和 x3 为广义坐标，则系统的动力学方程可表示为: 图 1 一维三弹簧谐振子系统微振动实验的示意图 m d2x1 dt2 + k1x1 + k2 x1 － x ( ) 2 = 0 m d2x2 dt2 + k2 x2 － x ( ) 1 + k3 x2 － x ( ) 3 = 0 m d2x3 dt2 + k3 x3 － x ( ) 2 =        0 (1) t = 0 时刻，系统的初位置和初速度分别定义 为: x1 = x10 x2 = x20 x3 = x{ 30 和 v1 = v10 v2 = v20 v3 = v{ 30 (2) 要得到系统在不同实验条件下的动力学方程，就必须求解上述的二阶常系数微分方程组，我们是通过建立 Matlab/Simulink 仿真模型来求解系统的微振动状态。 为了简化模型、减少参数，我们将方程组(1) 式改写为: d2x1 dt2 = － k1 m x1 + k2 m x2 － x ( ) 1 d2x2 dt2 = － k2 m x2 － x ( ) 1 + k3 m x3 － x ( ) 2 d2x3 dt2 = k3 m x2 － x ( )          3 (3) 打开 Matlab/S