单自由度振动分析

结构动力学三级项目班级：冶金五班小组成员：邱林凯李海洋张富张富增指导老师：王健2017年4月18日目录摘要 (2)单自由度系统的振动 (3)单自由度振动系统数学模型的建立 (3)参数设定与求解 (5)单自由度系统的强迫振动 (8)本章小结 (17)总结与心得 (17)摘要振动系统问题是个比较虚拟的问题，比较抽象的理论分析，对于问题的分析可以实体化建立数学模型，通过MATLAB可以转化成为图像。

单自由度频率、阻尼、振型的分析，我们可以建立数学模型，最后通过利用MATLAB编程实现数据图形；多自由度主要研究矩阵的迭代求解，我们在分析抽象的理论的同时根据MATLAB编程实现数据的迭代最后可以得到所要的数据，使我们的计算更加简便。

关键词：振动系统；单自由度；MATLAB；多自由度前言振动系统是研究机械振动的运动学和动力学，研究单自由系统的振动有着实际意义，因为工程上有许多问题通过简化，用单自由度系统的振动理论就能得到满意的结果。

模态是振动系统的一种固有振动特性，模态一般包含频率、振型、阻尼。

利用MATLAB编程并验证程序的正确性。

通过程序的运行，能快速获得多自由度振动系统的固有频率以及主振型，为设计人员提供了防止系统共振的理论依据，也为初步分析各构件的振动情况以及解耦分析系统响应奠定了基础。

在结构动力学中，单自由度系统的振动是最简单的运动，但这部分又十分重要。

因为从中可得到有关振动理论的一些基本的概念和解决问题的方法，同时它也适用于更为复杂的振动问题，是分析多自由度体系振动问题的基础。

因此，搞清楚了单自由度系统的振动，将有助于我们提高分析和解决其他各种振动问题的能力。

另外在实际工程中，确实有许多振动问题，可简化为单自由度问题，或近似地用单自由度理论去分析解决。

单自由度系统的振动单自由度振动系统数学模型的建立建立和分析有粘性阻尼时的自由度振动微分方程。

以静平衡位置为原点建立如图坐标，由牛顿定律得运动方程为[13]:0=++kx x c xm (2-1) 令mkm c n n ==2,2ω 其中n 称为衰减系数，单位为s 1；n ω是相应的无阻尼时的固有频率，式(2-1)可以写为：022=++x x n x n ω (2-2)如果进一步令xxc kxxmmcmnnως=(2-3)其中无量纲的ς称为相对阻尼系数，则式(2-2)可写为：022=++x x x n n ωςω (2-4)为了求解，令st e x = (2-6)代入(2-4)后得到特征方程：0222=++n n s s ωςω (2-7)他的两个特征根为：122,1-±-=ςωςωn n s (2-8)根据相对阻尼系数ς的不同大小，可以将阻尼分为三种状态：1>ς时为过阻尼，1=ς时为临界阻尼，10<<ς时为欠阻尼。

过阻尼状态1>ς，1s 与2s 是两个不等的负实根，令12*-=ςωωn （2-9）初始条件00)0(,)0(x xx x == （2-10） 系统初始条件响应为)()(**00*0t sh x xt ch x e t x n t n ωωςωωςω++=- （2-11）临界阻尼状态n s ως-==,1是二重根，方程（2-4）的通解为系统对式（2-10）的初始条件的响应为])([)(000t x xx e t x n t n ωω++=- （2-12） 欠阻尼状态1<ς，其中21ςωω-=n d （2-13）初始条件响应)sin cos ()(000t x xt x e t x d dn d t n ωωςωωςω++=- （2-14）参数设定与求解阻尼比ς分别取；应用Matlab 对式(2-11)和(2-12),（2-14）求解。

程序如下：clear,format compact;a=0.5;t=0:0.1:18;;w0=1;k=1;x0=1;wd=w0.*sqrt(1-a*a);x1=wdy=exp(-a*w0.*t).*(x0.*cos(wd.*t)+((x1+a*wd*x0)./wd)*sin(wd.*t)) figure(1),plot(t,y,'r');hold on a=1.0;t=0:0.1:18; w0=1;wd=1;x1=wd;y=exp(-wd.*t).*(x0+(x1+wd*x0).*t); figure(1),plot(t,y,'d');hold ona=2.0;t=0:0.1:18;w0=1;wd=w0*sqrt(a*a-1);y=exp(-a*w0.*t).*(x0.*cosh(wd.*t)+(x1+a*w0*x0)/w0.*sinh(t));figure(1),plot(t,y,'v');hold on 结论：图2-2为Matlab 计算后给出的响应曲线，从中可以得到一些重要的结论：在10<<ς的情况下，阶跃信号输入时，输出信号为衰减振荡，其振荡角频率（阻尼振荡角频率）为d ω，幅值按指数衰减越大，阻尼越大，衰减越快。

1>ς时，振荡系统等同于两个一阶系统串联。

此时虽然不产生振荡，但也需要经过较长时间才能达到稳态。

在一定的ς之下，欠阻尼系统能够更快地达到稳态值；而过阻尼系统反应迟饨，动作缓慢，所以系统通常设计成欠阻尼系统，ς取值为2024681012141618-0.20.20.40.60.811.21.4t(s)x (t )a=1.0a=2.0a=0.5图2-2算例绘制无阻尼单自由度系统的固有频率和周期随静变形的变化曲线。

固有频率n ω和周期n τstn gδω=，gstn δπτ2=取2/81.9s m g =。

可以利用下列MATLAB 程序画出st δ在0~0.5范围内nω和n τ的变换曲线：%Ex2_17.m g=9.81;for i=1:101 t(i)=0.01+(0.5-0.01)*(i-1)/100;w(i)=(g/t(i))^0.5; tao(i)=2*pi*(t(i)/g)^0.5; endplot(t,w);gtext('w_n'); hold on;plot(t,tao);gtext('T_n'); xlabel('delta_s_t'); title('Example2.1');0.050.10.150.20.250.30.350.40.450.505101520253035w nT ndelta stExample2.1单自由度系统的强迫振动简谐激励是激励形式中最简单的一种，虽然它在实际中存在的场合比较少但掌握系统对于简谐激励的响应的规律，是理解系统对周期激励或更一般形式激励的响应基础。

图所示的弹簧质量系统中，质量块上作用有简谐激振力t P t P ωsin )(0= (2-15)其中0P 为激振力幅，ω为激振频率。

以静平衡位置为坐标原点建立图示的坐标系。

从图的受力分析，得到运动微分方程为：t p kx x c xm ωsin 0=++ (2-16) 由常微分方程理论知道，方程(3.2)的通解x 由相应的齐次方程的通解h x 和非齐次方程的任意特解p x 两部分组成，即)()()(t x t x t x p h += (2-17)当欠阻尼时，式中)(t x h 为有阻尼自由振动，它的特点是振动频率为阻尼固有频率，振幅按指数规律衰减，称为瞬态振动或瞬态响应；)(t x p 是一种持续的等幅振动，它是由于简谐激励振力的持续作用而产生的，称之为稳态强迫振动或稳态振动，在间隔充分长时间考虑的振动就是这种稳态振动，而在刚受到外界激励时，系统的响应则是上述两种振动之和。

可见，系统受简谐激励后的响应可以分为两个阶段，一开始的过程称为过渡阶段，经充分长时间后，瞬态响应消失这时进xm P(tkxc ckx mP(t入过渡阶段，经充分长时间后，瞬态响应消失，这是进入稳态阶段。

将方程（2-15）的两端同除以质量m ，并且令22n mc ςω= （2-18） 其中ς为相对阻尼系数，n ω为相应的无阻尼系统的固有频率，则方程（2-15）成为t mp x x xn n ωωςωsin 202=++ （2-19） 上述方程特解可以通过)sin(ϕω-=t B x 或者t B t A x ωωsin cos +=来求得，这里介绍用复数方法求式（2-19）的特解。

先将式（2-19）写为下列的复数形式ti n n e mp x x xωωςω022=++ （2-20） 其中x 是复数设复数形式的特解为t i Be x ω= （2-21）其中B 称为复振幅，其意义是包含有相位的振幅。

将式（2-21）代入（2-20），解得ωςωωωn n i m P B 21220+-=（2-22）记λ为频率比，它定义为nωωλ=（2-23） 则式（2-22）可以写成ϕςλλςλλi i Be e kp i k P B --=+-=+-=222020)2()1(1211（2-24）式中2220)2()1(1ςλλ+-=kp B （2-25）2112λςλϕ-=-tg （2-26） 将式（2-24）代入（2-21），得到复数形式的特解为)(ϕω-=t i Be x （2-27）比较方程（2-17）与（2-18），可知（2-19）中的位移x 是（2-20）中复数x 的虚部，因此（2-25）的虚部就是方程（2-12）的特解，即有)sin(ϕω-=t B x （2-28）其中B 为振幅，ϕ为相位差。

由式（2-26）、2-23）及（2-24）得出稳态强迫振动有如下的基本特点：1)线性系统对简谐激励的稳态响应是频率等同于激励频率而相位滞后与激振力的简谐振动；2)稳态响应的振幅及相位差只取决于系统本身的物理性质和激振的频率及力幅，而与系统本身进入运动的方式无关。

无阻尼系统对简谐激励的稳态响应可以从式（2-26）得出。

当n ωω<时，得到1<λ，0=ϕ，这时t k P x ωλsin 1120-=（2-29） 当n ωω>时，得到1>λ，πϕ=，这时)sin(1120πωλ--=t k P x （2-23）式（2-21）也可以写成（2-22）的形式，这时相位差反映在振幅2011λ-k P 的符号中。

上述结果也可以由直接设t B x ωsin =并代入下列方程而得到：t P kx xm ωsin 0=+ （2-24） 为了具体讨论影响稳定响应的振幅和相位差的各种因素，记kP B 00=（2-25） 0B 实际是质量块在激振力幅静作用下的最大位移。

再引入无量纲的振幅放大因子β，它定义为2220)2()1(1ςλλβ+-==B B （2-26） 由式（2-26）和（2-19）可以分别画出以相对阻尼系数ς为参数的曲线——λβ-曲线与λϕ-曲线，前者称为幅频响应曲线，后者称为相频响应曲线如图所示 程序如下for kesai=[0.05,0.10,0.15,0.25,0.375,0.50,1.0] lamda=0:0.01:5.0;beta=1./(sqrt((1-lamda.^2).^2+(2*kesai*lamda).^2)); plot(lamda,beta) hold on endaxis([0 5 0 3]);00.51 1.522.533.544.550.511.522.53频率比振幅放大因子1.00.50.3750.050.150.25偏心质量引起的强迫振动振幅放大因子2222)2()1(ςλλλ+-=me MB (2-27)程序如下：for kesai=[0.05,0.10,0.15,0.25,0.375,0.50] lamda=0:0.01:5.0;beta=lamda./(sqrt((1-lamda.^2).^2+(2*kesai*lamda).^2)); plot(lamda,beta) hold on endaxis([0 5 0 3]);00.51 1.522.533.544.550.511.522.53频率比M B /m e0.50.100.150.250.3750.501.0支撑运动引起的强迫振动振幅放大因子2222)2()1()2(1ςλλςλβ+-+==a B (2-28) 程序如下：for kesai=[0.05,0.10,0.15,0.25,0.375,0.50,1.0]lamda=0:0.01:5.0;beta=sqrt((1+(2*kesai*lamda).^2)./((1-lamda.^2).^2+(2*kesai*lamda).^2));plot(lamda,beta)hold on endaxis([0 5 0 3]);00.51 1.522.533.544.550.511.522.53频率比振幅放大因子0.050.100.150.250.3750.501.01.141算例利用MATLAB ，绘制弹簧-质量系统在简谐力作用下的响应曲线。