基于Matlab-Simulink的机械振动仿真研究

1、目 录1引言机械振动的仿真原理11.1 Matlab Simulink功能简述11.2机械振动的物理模型1简谐振动的物理模型1阻尼振动的物理模型2受迫振动的物理模型21.3 Matlab Simulink仿真原理简述42简谐振动方程的解与其模拟仿真52.1简谐振动方程的求解52.2简谐振动模型的仿真研究5基本模型的建立52.2.2 速度、加速度的监测72.2.3 动能、势能、机械能监测82.3简谐振动的图像分析93阻尼振动方程的求解和仿真模拟113.1弹簧振子做阻尼振动方程的求解113.2弹簧振子做阻尼振动的模拟仿真研究113.3阻尼振动的图像分析144受迫振动的方程的求解和仿真模拟164.1弹簧振子做受迫振动方程的求解164.2弹簧振子做受迫振动的仿真模拟研究17策动力频率时弹簧振子的受迫振动仿真模拟17策动力频率时弹簧振子受迫振动的仿真模拟20策动力频率时弹簧振子的仿真模拟224.3受迫振动的图像分析235几点补充说明与仿真模拟中问题分析255.1物理振动模型建立的补充说明255.2 方程求解中的补充说明255.3 仿真模拟中的问题分析256结语27参考文献28附录29致谢30摘

2、要机械振动主要有简谐振动,阻尼振动,受迫振动三种.对三种振动建立模型,列出振动方程,再对三种振动给定初始条件,就可以利用Matlab Simulink功能对三种振动进行仿真模拟,得出振动的位移,速度,加速度,动能,势能,机械能随时间的变化关系图像.另外,我们对振动方程求解,得出振子位移关于时间的函数,再分别对其求一阶、二阶导数,就可以得出速度、加速度函数,再经过简单运算就可以得到动能、势能、机械能函数.我们再通过分析函数来分析其图像,再对比仿真模拟出的图像,就可以确定我们的仿真研究方法的可信度.关键词：简谐振动；阻尼振动；受迫振动；共振AbstractThe mechanical vibration mainly includes three kinds of vibrations: the simple harmonic oscillator, the pure damping vibration and the forced damping vibration. According to the three vibrations, we can build models for t

3、hem, list the vibration equations. After giving the initial conditions, we can carry out analog simulations to test the three vibrations by using Matlab Simulink functions, and obtain the images of vibration displacement, velocity, acceleration, kinetic energy, potential energy, mechanical energy variation over time. In addition, we e to the function of oscillator displacement over time after we seek the solution of equations. And then we can gain the derivatives of their first order and second

4、order respectively, so we can draw functions of speed and acceleration. Furthermore, we can get the functions of kinetic energy, potential energy, mechanical energy by a simple operation. Let us analyze the images by analyzing the functions, and then we make a parison between the images of analog simulations and them. Finally, we can determine the credibility of the simulation methodology. Keywords: Simple harmonic oscillator; Damping vibration; Forced vibration; Resonance / 1引言机械振动的仿真原理1.1 Matl

5、ab Simulink功能简述Simulink是基于Matlab的框图设计环境,可以用来对各种动态系统进行建模、分析和仿真,它的建模范围广泛,可以针对任何能用数学来描述的系统进行建模,例如航空航天动力学系统、卫星控制制导系统、通信系统、船舶与汽车等,其中包括了连续、离散,条件执行,事件驱动,单速率、多速率和混杂系统等.Simulink提供了利用鼠标拖放的方法来建立系统框图模型的图形界面,而且还提供了丰富的功能块以与不同的专业模块集合,利用Simulink几乎可以做到不书写一行代码即完成整个动态系统的建模工作.除此之外,Simulink还支持Stateflow,用来仿真事件驱动过程.Simulink是从底层开发的一个完整的仿真环境和图形界面,是模块化了的编程工具,它把Matlab的许多功能都设计成一个个直观的功能模块,把需要的功能模块用连线连起来就可以实现需要的仿真功能了.也可以根据自己的需要设计自己的功能模块,Simulink功能强大,界面友好,是一种很不错的仿真工具1.1.2机械振动的物理模型物理学中的机械振动主要分为简谐振动、阻尼振动、受迫振动三种.下面我们根据这三种类型的振动建立

6、物理模型来分别研究.1.2.1简谐振动的物理模型图1弹簧振子做简谐振动物理实验模型如上图所示,弹簧振子在附近做简谐振动.已知弹簧振子质量为,所受合力为,弹簧劲度系数为,则有：.又由牛顿第二定律有：于是可以得到： 令,则可得： 方程的解即为弹簧振子在时刻时的振动位移,一阶导数即为弹簧振子在时刻时振动速度,其二阶导数即为弹簧振子在时刻时的加速度.1.2.2阻尼振动的物理模型如图1,若弹簧振子在x轴上受到粘滞阻尼的作用力,则弹簧振子做阻尼的振动.设弹簧振子受到的阻尼力为：式中g为阻尼系数,与物体的形状以与周围性质有关.弹簧振子受到的弹力为,则对弹簧振子,有牛顿第二定律有：整理后得： 令,则有：这就是阻尼振动的振动方程.其解即为弹簧振子在时刻t时的振动位移,一阶导数即为弹簧振子在时刻t时振动速度,其二阶导数即为弹簧振子在时刻t时的加速度.1.2.3受迫振动的物理模型如图2,弹簧振子在附近做阻尼振动.已知弹簧振子质量为,弹簧劲度系数为.平行于轴的平面对弹簧振子有阻尼力的作用.对弹簧振子施加一外加激励力,设,则称为谐激励力,其中为外施激励频率,t是持续时间. 对弹簧振子受力分析,其所受弹力为：.由

7、于阻尼振动是振幅随时间不断减少的振动.能量减少的原因是有粘滞阻尼和辐射阻尼.为方便,均视为粘滞阻尼.则弹簧图2弹簧振子在外加激励力作用下做阻尼受迫振动振子所受阻尼力为： 式中g为阻尼系数,与物体的形状以与周围性质有关.则对弹簧振子,由牛顿第二定律有：对式变形可得:令,为固有频率,为阻尼因数,则可变为：方程的解就是时刻时弹簧振子的位移,其一阶导数即为弹簧振子在时刻时振动速度,其二阶导数即为弹簧振子在时刻时的加速度2.我们记为相对阻尼系数或阻尼比.根据阻尼对系统振动的影响,振动响应分为弱阻尼、强阻尼和临界阻尼三种情况,这里仅讨论弱阻尼的情况.1.3 Matlab Simulink仿真原理简述在得到弹簧振子的简谐振动、阻尼振动和受迫振动方程后,通过这三个方程,我们可以用高等数学的方法求出这三个方程的通解.同时,我们可以用Matlab的计算功能求出它们的通解.这三个方程的通解表示振子位移随时间的变化情况.我们得到的这三个方程,前两个为二阶常系数线性齐次微分方程,第三个为二阶常系数非齐次微分方程.根据这三个方程,我们可以通过MatlabSimulink中的各种模块模拟弹簧振子的位移、速度、加速度

8、,再添加一个平方模块,设置好系数,就可以模拟振子动能、势能、机械能,用线连接各模块,这样流程图就做好了.设置好各模块的参数后,再设置好系统环境变量,点击运行,通过示波器模块就可以模拟出相应的图像曲线3.图像的横坐标均表示时间,纵坐标相应为位移、速度、加速度、动能、势能、机械能.图像表示这些物理量随时间变化关系.通过这三种情况方程的通解,我们可以分析振子位移随时间变化情况,再和模拟出的图像对比分析.对方程通解求一阶导,就可以得到振子速度随时间变化关系,分析出速度随时间变化情况,再和模拟出的图像对比分析.同样我们可以求出方程通解的二阶导数,这就是振子加速度随时间变化关系,分析函数特征,再和模拟出的函数图像对比分析,就可以分析模拟出的图像是否正确,与其与理论符合情况.2简谐振动方程的解与其模拟仿真2.1简谐振动方程的求解这里,我们设系统初始条件为s时,.通过高等数学方法解这个齐次微分方程可得：式中.则速度表达式为:,将初始条件代入式,可得：这就是满足初始条件的简谐振动方程的解.由式我们可以得出弹簧振子位移随时间的变化情况.振子周期为.时,振子位移正向最大位移出,即图1中的A位置,此时振子速度

9、为0,加速度最大；经,振子向负方向运动到平衡位置,此时振子速度最大,加速度为0；再经,振子继续向负方向运动到负的最大位移处,此时速度为0,加速度最大；再经过,振子向正向运动到平衡位置,此时速度最大,加速度为0；最后经过,振子回到初始位置,即正的最大位移处,完成一个周期的振动.通过matlab7.0符号运算,可以得出该微分方程的解,相关程序见附录程序1.2.2简谐振动模型的仿真研究2.2.1基本模型的建立我们设系统的固有频率,则.于是式变为：.打开SimulinkLibraryBrowser,选择新建按钮,根据所需要模拟的运动方程选取模块,其中包括Subtract、Intergrator、Gain以与Scope模块,需要注意的是将Subtract模块中的Listofsigns改为-,以便让前面的符号为负,为了使前的系数为0.5,将Gainl中的值设为0.5,速度项系数Gain设为0.设置位移模块的初始值为4,速度模块的初始值设为0,加速度模块的初始值设为0.这样,几个关键模块的属性就根据方程的需要设置好了.运用Gain1将和0.5相乘.运用Subtract使前面的符号为负.运用Intergrator将积分为,将积分为.Scope为示波器输出模块.最后,将各个模块按照方程的需要逐一连接,如图3所示4. 图3简谐振动位移仿真模拟流程图 图4简谐振动位移仿真模拟图像单击右键,选择Configuration Parameters设置系统的运行环境,初始运

《基于Matlab-Simulink的机械振动仿真研究》由会员m****分享，可在线阅读，更多相关《基于Matlab-Simulink的机械振动仿真研究》请在金锄头文库上搜索。