单自由度系统机械振动 1. 图示系统的轮和绳之间无相对滑动，只作纯 滚动，建立系统的运动微分方程，并求系统 的固有频率，圆盘转动惯量为J ，质量块的 质量为m ，弹簧刚度为K 。 2． 图所示，W=1000N ，k=2 104N/m ，图示位 置弹簧已承受初压力F 0=100N ，现将支承突 然撤去，重块落下后作自由振动时的振动位 移表达式？（取重力加速度g=10m/s 2） 3．如图所示为一台机器，其总质 量为M ，安装在一个弹簧和一 个阻尼器上，弹簧常数为k ，阻 尼系数为c 。机器工作时旋转中 心为O ，角速度为ω，不平衡 质量大小为m ，偏心距离为e 。 机器只能在垂直方向运动。求机器振动时传给地面的力的最大值。 W K
4．图示系统中，质量m 上受激励力为 F （t ）=sin ωt+10sin10ωt 时， 求质量m 的稳态响应 5. 图示系统的轮和绳之间无相对滑动，只作纯滚动，建立系统的运动微分方程，并求系统的固 有频率，圆盘转动惯量为J ，质量块的质量为m ， 弹簧刚度为K 6. 一重块与两弹簧相连，W=490N ，k=9800N/m ， 图示位置弹簧不受力，现将支承突然撤去，重块 落下后作自由振动时的振动位移表达式？ 7. 如图所示为一台机器，其总质量为m ，通过一个弹簧和一个阻尼器安装在基础上，弹 簧常数为k ，阻尼系数为c 。基础的运动为 y(t)=Ysin ωt ，机器只能在垂直方向运动。求 基础振动时传给机器的力的最大值。 W K K
8．图示系统中，质量m上受激励力为 F（t）=sinωt+10sin10ωt时， 求质量m的稳态响应。 9.一般振动问题，如图所示： 三类振动问题分别是： （1）振动分析，已知，求； （2）振动环境预测或载荷分析，已知，求； （3）系统识别，已知，求。 10. 振动问题的分类，根据自由度数分，有， 和。 11. 简谐振动x=Asin(ωt+φ)，其中的振动位移为，振幅 为， 振动频率为为，振动的初相位为 12. n个自由度振动系统有个固有频率，有个固有 振型， 其中的第i阶主振型有个节点。
目 录 1引言——机械振动的仿真原理 (1) 1.1 Matlab Simulink 功能简述 (1) 1.2机械振动的物理模型 (1) 1.2.1简谐振动的物理模型 (1) 1.2.2阻尼振动的物理模型 (2) 1.2.3受迫振动的物理模型 (2) 1.3 Matlab Simulink 仿真原理简述 (4) 2简谐振动方程的解及其模拟仿真 (5) 2.1简谐振动方程的求解 (5) 2.2简谐振动模型的仿真研究 (5) 2.2.1基本模型的建立 (5) 2.2.2 速度、加速度的监测 (7) 2.2.3 动能、势能、机械能监测 (8) 2.3简谐振动的图像分析 (9) 3阻尼振动方程的求解和仿真模拟 (11) 3.1弹簧振子做阻尼振动方程的求解 (11) 3.2弹簧振子做阻尼振动的模拟仿真研究 (11) 3.3阻尼振动的图像分析 (14) 4受迫振动的方程的求解和仿真模拟 (16) 4.1弹簧振子做受迫振动方程的求解 (16) 4.2弹簧振子做受迫振动的仿真模拟研究 (17) 4.2.1策动力频率0ωω<时弹簧振子的受迫振动仿真模拟 (17) 4.2.2策动力频率0ωω>时弹簧振子受迫振动的仿真模拟 (20) 4.2.3策动力频率0ωω=时弹簧振子的仿真模拟 (22) 4.3受迫振动的图像分析 (23) 5几点补充说明与仿真模拟中问题分析 (25) 5.1物理振动模型建立的补充说明 (25)
5.2方程求解中的补充说明 (25) 5.3仿真模拟中的问题分析 (25) 6结语 (27) 参考文献 (28) 附录 (29) 致 (30)
摘要 机械振动主要有简谐振动，阻尼振动，受迫振动三种。对三种振动建立模型，列出振动方程，再对三种振动给定初始条件，就可以利用Matlab Simulink功能对三种振动进行仿真模拟，得出振动的位移，速度，加速度，动能，势能，机械能随时间的变化关系图像。另外，我们对振动方程求解，得出振子位移关于时间的函数，再分别对其求一阶、二阶导数，就可以得出速度、加速度函数，再经过简单运算就可以得到动能、势能、机械能函数。我们再通过分析函数来分析其图像，再对比仿真模拟出的图像，就可以确定我们的仿真研究方法的可信度。 关键词：简谐振动；阻尼振动；受迫振动；共振
一、填空题 1、机械振动按不同情况进行分类大致可分成（线性振动）和非线性振动；确定性振动和（随机振动）；（自由振动）和强迫振动。 2、周期运动的最简单形式是（简谐运动），它是时间的单一（正弦）或（ 余弦）函数。 3、单自由度系统无阻尼自由振动的频率只与（质量）和（刚度）有关，与系统受到的激励无关。 4、简谐激励下单自由度系统的响应由（瞬态响应）和（稳态响应）组成。 5、工程上分析随机振动用（数学统计）方法，描述随机过程的最基本的数字特征包括均值、方差、（自相关函数）和（互相关函数）。 6、单位脉冲力激励下，系统的脉冲响应函数和系统的（频响函数）函数是一对傅里叶变换对，和系统的（传递函数）函数是一对拉普拉斯变换对。 2、在离散系统中，弹性元件储存( 势能)，惯性元件储存（动能），（阻尼）元件耗散能量。 4、叠加原理是分析（线性）系统的基础。 5、系统固有频率主要与系统的（刚度）和（质量）有关，与系统受到的激励无关。 6、系统的脉冲响应函数和（频响函数）函数是一对傅里叶变换对，和（传递函数）函数是一对拉普拉斯变换对。 7、机械振动是指机械或结构在平衡位置附近的（往复弹性）运动。 1．振动基本研究课题中的系统识别是指根据已知的激励和响应特性分析系统的性质，并可得到振动系统的全部参数。（本小题2分） 2．振动按激励情况可分为自由振动和强迫振动两类。（本小题2分）。 3．图（a）所示n个弹簧串联的等效刚度= k ∑ = n i i k1 1 1 ；图（b）所示n个粘性阻尼串联的等效粘 性阻尼系数= e C ∑ = n i i c1 1 1 。（本小题3分） （a）（b） 题一 3 题图 4．已知简谐振动的物体通过距离静平衡位置为cm x5 1 =和cm x10 2 =时的速度分别为s cm x20 1 = &和s cm x8 2 = &，则其振动周期= T；振幅= A10.69cm。（本小题4分） 5．如图（a）所示扭转振动系统，等效为如图（b）所示以转角 2 ?描述系统运动的单自由度 系统后，则系统的等效转动惯量= eq I 2 2 1 I i I+，等效扭转刚度= teq k 2 2 1t t k i k+。（本小题4分）
第2章多自由度系统的振动 基本要点： ①建立系统微分方程的几种方法； ②固有频率、固有振型的概念以及固有振型关于质量和刚度矩阵的加权正交性； ③多自由度系统运动的解耦—模态坐标变换及运用模态叠加法求解振动系统的响应。 引言 多自由度振动系统的几个工程实例；多自由度系统振动分析的特点；多自由度系统振动分析与单自由度系统的区别与联系。 §2.1多自由度系统的振动方程 ●方程的一般形式：质量矩阵、阻尼矩阵、刚度矩阵和激振力 §2.2建立系统微分方程的方法 ●影响系数：刚度影响系数、柔度影响系数 ●刚度矩阵法、柔度矩阵法及这两种方法的特点；Lagrange方程法 §2.3无阻尼系统的自由振动 ●二自由度系统的固有振动：固有频率、固有振型。 ●二自由度系统的自由振动 ●二自由度系统的运动耦合与解耦 弹性耦合，惯性耦合； 振动系统的耦合取决于坐标系的选择； ●多自由度系统的固有振动 固有振动的形式及条件：特征值、特征向量、模态质量、模态刚度； 固有振型的性质：关于质量矩阵和刚度矩阵的加权正交性； 刚体模态； ●运动的解耦：模态坐标变换（主坐标变换）。 ●多自由度系统的自由振动 §2.4无阻尼系统的受迫振动 ●频域分析：动刚度矩阵和频响函数矩阵，频响函数矩阵的振型展开式，系统反 共振问题。 ●时域分析：单位脉冲响应矩阵，任意激励下的响应，模态截断问题，模态加速 度法。 §2.5比例阻尼系统的振动 ●多自由度系统的阻尼：Rayleigh比例阻尼。 ●自由振动 ●受迫振动：频响函数矩阵，单位脉冲响应矩阵，任意激励下的响应。 §2.6一般粘性阻尼系统的振动
●自由振动：物理空间描述，状态空间描述。 ●受迫振动：脉冲响应矩阵，频响函数矩阵，任意激励下的响应。 思考题： ①刚度矩阵和柔度矩阵在什么条件下是互逆的两个矩阵？从物理上和数学两方面加以解 释？ ②为什么说模态质量、模态刚度的数值大小没有直接意义？ ③证明固有振型关于质量矩阵和刚度矩阵的加权正交性，并讨论其物理意义。 ④在实际的多自由度系统振动分析中，为什么要进行模态截断？ 参考书目 1.胡海岩，机械振动与冲击，航空工业出版社，2002 2.故海岩，机械振动基础，北京航空航天大学出版社，2005 3.季文美，机械振动，科学出版社，1985。（图书馆索引号：TH113.1/1010） 4.郑兆昌主编, 机械振动上册,机械工业出版社，1980。（图书馆索引号： TH113.1/1003-A） 5.Singiresu S R, Mechanical vibrations，Longman Prentice Hall, 2004(图书馆索引 号：TH113.1/WR32)
习 题 1-1一单层房屋结构可简化为题1-1图所示的模型，房顶质量为m ，视为一刚性杆；柱子高h ，视为无质量的弹性杆，其抗弯刚度为EJ 。求该房屋作水平方向振动时的固有频率。 解：由于两根杆都是弹性的，可以看作是两根相同的弹簧的并联。 等效弹簧系数为k 则 mg k δ= 其中δ为两根杆的静形变量，由材料力学易知 δ=3 24mgh EJ = 则 k = 3 24EJ h 设静平衡位置水平向右为正方向，则有 " m x kx =- 所以固有频率3 n 24mh EJ p = 1-2 一均质等直杆，长为 l ，重量为W ，用两根长h 的相同的铅垂线悬挂成水平位置，如题1-2图所示。试写出此杆绕通过重心的铅垂轴作微摆动的振动微分方程，并求出振动固有周期。 解：给杆一个微转角θ 2 a θ＝h α 2F cos α＝mg 由动量矩定理： a h a mg a mg Fa M ml I M I 822cos sin 12 1 2 2-=-≈?-=== =αθ αθ&& 题1-1图 题1-2图 θ F sin α 2 θα h mg θ
其中 12 cos sin ≈≈θ α α h l ga p h a mg ml n 2 2 2 2 2304121==?+θθ&& g h a l ga h l p T n 3π23π2π22 2=== 1-3求题1-3图中系统的固有频率，悬臂梁端点的刚度分别是k 1和k 3，悬臂梁的质量忽略不计。 解：悬臂梁可看成刚度分别为k 1和k 3的弹簧，因此，k 1与k 2串联，设总刚度为k 1ˊ。k 1ˊ与k 3并联，设总刚度为k 2ˊ。k 2ˊ与k 4串联，设总刚度为k 。即为 21211k k k k k += '，212132k k k k k k ++='，4 241213231421432421k k k k k k k k k k k k k k k k k k k k ++++++= ) (42412132314 214324212k k k k k k k k k k m k k k k k k k k k p ++++++= 1-4求题1-4图所示的阶梯轴一圆盘系统扭转振动的固有频率。其中J 1、J 2和J 3是三个轴段截面的极惯性矩，I 是圆盘的转动惯量，各个轴段的转动惯量不计，材料剪切弹性模量为G 。 解： 111/l GJ k = （1） 222/l GJ k = （2） 333/l GJ k = （3） )/(23323223l J l J J GJ k += （4） ) (/)()4)(3)(2(1/)(2332113221332122312l J l J Il l J J l J J l J J G P I k k P n n +++=+=知 ）由（ 题1-3图 题1-4图
习 题 2-1已知系统的弹簧刚度k =800 N/m ，作自由振动时的阻尼振动周期为1.8s ，相邻两振幅的比值 1 2 .41=+i i A A ，若质量块受激振力t t F 3cos 360)(=N 的作用，求系统的稳态响应。 解：由题意，可求出系统的运动微分方程为 t m x n x p x n 3cos 360 22 =++ 得到稳态解 )3cos(α-=t B x 其中 m k B B B 45.0360 4)1(02 2220 == +-= λζλ 222 122tg λζλ ωωα-=-= n p n 由 d nT i i A A e 2.41 === +η 489 .3π 2797 .0ln 8 .1ln ======d d d d d T p T n T nT η η 又 22n p p n d -= 有 579.32 22=+=n d n p n p p 45.51255.1298.0374 .0838 .01838.0223.02tg 103.1408 .045 .0838.0223.04)838.01(45 .0223.0579 .3797.0838.0579 .33 2 222===-??= == ??+-= === == =ααζω λB p n p n n 所以 x ＝1.103 cos(3t －51?27') 2-2一个无阻尼弹簧质量系统受简谐激振力作用，当激振频率ω1 =6rad/s 时，系统发生共振；给
质量块增加1 kg 的质量后重新试验，测得共振频率ω2 =5.86rad/s ，试求系统原来的质量及弹簧刚度。 解：设原系统的质量为m ，弹簧常数为k 由 m k p n = ，共振时m k p n ==1ω 所以 m k =6 ① 又由 当 86.51 2=+= =m k p n ω ② ①与②联立解出 m ＝20.69 kg ，k ＝744.84 N/m 2-3总质量为W 的电机装在弹性梁上，使梁产生静挠度st δ，转子重Q ，重心偏离轴线e ，梁重及阻尼可以不计，求转速为ω时电机在垂直方向上稳态强迫振动的振幅。 解：列出平衡方程可得： 222()sin sin()sin()st Q W W k x w e wt x g g W Q x kx w e wt g g kg Q x x w e wt W W ππ-σ+- =+=++=+ 所以：2n kg P W Q h w e W ==， 又因为st st W W k k =σ=σ即 22() st st B w e B W g w =σ-σ将结果代入Q = 即为所求的振幅 2-4如题2-4图所示，作用在质量块上的激振力t F t F ωsin )(0=，弹簧支承端有运动 t a x s ωco s =，写出系统的运动微分方程，并求稳态振动。 题2-4图
摘要 机械振动主要有简谐振动，阻尼振动，受迫振动三种。对三种振动建立模型，列出振动方程，再对三种振动给定初始条件，就可以利用Matlab Simulink功能对三种振动进行仿真模拟，得出振动的位移，速度，加速度，动能，势能，机械能随时间的变化关系图像。另外，我们对振动方程求解，得出振子位移关于时间的函数，再分别对其求一阶、二阶导数，就可以得出速度、加速度函数，再经过简单运算就可以得到动能、势能、机械能函数。我们再通过分析函数来分析其图像，再对比仿真模拟出的图像，就可以确定我们的仿真研究方法的可信度。 关键词：简谐振动；阻尼振动；受迫振动；共振
1引言——机械振动的仿真原理 1.1 Matlab Simulink功能简述 Simulink是基于Matlab的框图设计环境，可以用来对各种动态系统进行建模、分析和仿真，它的建模范围广泛，可以针对任何能用数学来描述的系统进行建模，例如航空航天动力学系统、卫星控制制导系统、通信系统、船舶及汽车等，其中包括了连续、离散，条件执行，事件驱动，单速率、多速率和混杂系统等。Simulink提供了利用鼠标拖放的方法来建立系统框图模型的图形界面，而且还提供了丰富的功能块以及不同的专业模块集合，利用Simulink几乎可以做到不书写一行代码即完成整个动态系统的建模工作。除此之外，Simulink还支持Stateflow，用来仿真事件驱动过程。 Simulink是从底层开发的一个完整的仿真环境和图形界面，是模块化了的编程工具，它把Matlab的许多功能都设计成一个个直观的功能模块，把需要的功能模块用连线连起来就可以实现需要的仿真功能了。也可以根据自己的需要设计自己的功能模块，Simulink功能强大，界面友好，是一种很不错的仿真工具[1]。 1.2机械振动的物理模型 物理学中的机械振动主要分为简谐振动、阻尼振动、受迫振动三种。下面我们根据这三种类型的振动建立物理模型来分别研究。 1.2.1简谐振动的物理模型
习 题 4-1 在题3-10中，设m 1=m 2=m ，l 1=l 2=l ，k 1=k 2=0，求系统的固有频率和主振型。 解：由题3-10的结果 22121111)(l g m l g m m k k +++ =，2 221l g m k -=，2212l g m k - =，2 2222l g m k k += 代入m m m ==21，021==k k ，l l l ==21 可求出刚度矩阵K 和质量矩阵M ??? ???=m m M 00；?? ?? ??????- - =l mg l mg l mg l mg K 3 由频率方程02=-M p K ，得 0322 =????? ??? ? ?-- - -=mp l mg l mg l mg mp l mg B 0242 2 2224 2 =+-∴l g m p l g m p m l g p ) 22(1-=∴ ，l g p )22(2+= 为求系统主振型，先求出adjB 的第一列 ???? ? ? ?? ??-=l mg mp l mg adjB 2 分别将频率值21p p 和代入，得系统的主振型矩阵为 ??????-=112) 1(A ?? ????+=112)2(A 题4-1图
4-2 题4-2图所示的均匀刚性杆质量为m 1，求系统的频率方程。 解：设杆的转角θ和物块位移x 为广义坐标。利用刚度影响系数法求刚度矩阵k 。 设0,1==x θ，画出受力图，并施加物体力偶与力 2111,k k ，由平衡条件得到， 222111a k b k k +=， a k k 221-= 设1,0==x θ，画出受力图，并施加物体力偶与力2212,k k ，由平衡条件得到， 12k a k 2-=， a k k 222= 得作用力方程为 ?? ? ???=??????????? ?--++????????????? ?00003122222 2122 1x a k a k a k a k b k x m a m θθ 由频率方程02=-M K p ，得 031 2 22222 212221=----+p m a k a k a k p a m a k b k 4-3 题4-3图所示的系统中，两根长度为l 的均匀刚性杆的质量为m 1及m 2，求系统的刚度矩阵和柔度矩阵，并求出当m 1=m 2=m 和 k 1=k 2=k 时系统的固有频率。 解：如图取21,θθ为广义坐标，分别画受力图。由动量矩定理得到， l l k l l k I 4 34343432 11111θθθ+-= 2 2434343432 2211122l l k l l k l l k I θθθθ--= 题4-3图 题4-2图
《机械振动基础》大作业 （2015年春季学期） 题目基于MATLAB求系统特性 姓名 学号 班级 专业机械设计制造及其自动化 报告提交日期
哈尔滨工业大学
报告要求 1.请根据课堂布置的2道大作业题，任选其一，拒绝雷同和抄袭； 2.报告最好包含自己的心得、体会或意见、建议等； 3.报告统一用该模板撰写，字数不少于3000字，上限不限； 4.正文格式：小四号字体，行距为倍行距； 5.用A4纸单面打印；左侧装订，1枚钉； 6.课程报告需同时提交打印稿和电子文档予以存档，电子文档由班 长收齐，统一发送至：。 7.此页不得删除。 评语： 成绩（15分）：教师签名： 年月日
求解多自由度矩阵的认识体会。 二、MATLAB程序图： >> m=[]; k1=[]; k=[]; c=[]; c1=[]; for i=1:9 a=input('输入质量矩阵m：'); m(i,i)=a; end ; for j=1:9 b=input('输入刚度系数k：'); k1(1,j)=b; end for l=1:8 k(l,l)=k1(l)+k1(l+1); k(9,9)=k1(9); k(l+1,l)=-k1(l+1); k(l,l+1)=-k1(l+1); k(9,8)=-k1(9);
k(8,9)=-k1(9); end ; syms w; B=k-w^2*m %系统的特征矩阵B Y=det(B); %展开行列式 W=solve(Y); %求解wh lW=length(W); [V,D]=eig(k,m); for I=1:9 for J=1:9 V(J,I)=V(J,I)/V(5,I); end end V W 三MATLAB结果输入输出： 程序输入内容： 输入质量矩阵m：1 输入质量矩阵m：2 输入质量矩阵m：3 输入质量矩阵m：4 输入质量矩阵m：5 输入质量矩阵m：6 输入质量矩阵m：7 输入质量矩阵m：8 输入质量矩阵m：9 输入刚度系数k：10 输入刚度系数k：11 输入刚度系数k：12 输入刚度系数k：13 输入刚度系数k：14 输入刚度系数k：15 输入刚度系数k：16 输入刚度系数k：17 输入刚度系数k：18
:单自由度系统自由衰减振动及固有频率、阻尼比的测定实验指导书 陈安远 （武汉大学力学实验教学中心） 1.实验目的 1、了解单自由度系统模型的自由衰减振动的有关概念； 2、学习用频谱分析信号的频率； 3、学习测试单自由度系统模型阻尼比的方法。 2.实验仪器及安装示意图 实验仪器:INV1601B型振动教学实验仪、INV1601T型振动教学实验台、加速度传感器、MSC-1力锤(橡胶头)、重块。 软件:INV1601型DASP软件。 图1实验系统示意图 3实验原理 单自由度系统的阻尼计算，在结构和测振仪器的分析中是很重要的。阻尼的计算常常通过衰减振动的过程曲线（波形）振幅的衰减比例来进行计算。衰减振动波形示于图2。用衰减波形求阻尼可以通过半个周期的相邻两个振幅绝对值之比，或经过一个周期的两个同方向
振幅之比,这两种基本方式进行计算。通常以一个周期的相邻两个振幅值之比为基准来计算的较多。两个相邻振幅绝对值之比,称为波形衰减系数。 图2衰减振动波形 1、对经过一个周期为基准的阻尼计算 每经过一个周期的振幅的比值为一常量： η=d nT i i e A A =+1 这个比例系数η表示阻尼振动的振幅(最大位移)按几何级数递减。衰减系数η常用来表示振幅的减小速率。叫做振幅减缩率或减幅系数。 如果用减幅系数η的自然对数来表示振幅的衰减则更加方便。 δ=ln （η）=ln d i i nT A A =+1=21ξπξ- δ称为振动的对数衰减率或对数减幅系数。可以利用δ来求得阻尼比ξ。 2、在小阻尼时,由于η很小；这样读数和计算误差较大,所以一般地取相隔若干个波峰序号的振幅比来计算对数衰减率和阻尼比。 4.实验步骤 1、仪器安装 参照仪器安装示意图安装好配重质量块，加速度传感器。 2、开机进入INV1601型DASP 软件的主界面, 进入单通道示波状态进行波形和频谱同时示波，见图2。 3400Hz 、采样点数为2K,标定值和工程单位等参数(按实际
第二章 单自由度系统的自由振动 本章以阻尼弹簧质量系统为模型,讨论单自由度系统的自由振动。 §2-1 无阻尼系统的自由振动 无阻尼单自由度系统的动力学模型如图所示。设质量为m ，单位是kg 。弹簧刚度为K ，单位是N ／m ，即弹簧单位变形所需的外力。弹簧在自由状态位置如图中虚线所示。当联接质量块后，弹簧受重力W=mg 作用而产生拉伸变形：，同时也产生弹簧恢复力K ，当其等于重力W 时，则处于静平衡位置，即 W=K 若系统受到外界某种初始干扰，使系统静平衡状态遭到破坏．则弹簧力不等于重力，这种不平衡的弹性恢复力，便使系统产生自由振动。首先建立座标，为简便起见，可选静平衡位置为座标原点，建立铅垂方向的座标x ，从原点算起，向下为正，向上为负，表示振动过程中质量块的位置。现设质量m 向下运动 到x ，此时弹簧恢复力为K(+x)，显然大于重力W ，由 于力不平衡，质量块在合力作用下，将产生加速度运动，故可按牛顿运动定律(作用于一个质点上所有力的合力，等于该质点的质量和沿合力方向的加速度的乘 积)，建立运动方程，取与x 正方向一致的力、加速度、速度为正，可列如下方程 改写为 0=+kx x m && (1-1-1 令 m k p = 2 (1-1-2) 单自由度无阻尼系统自由振动运动方程为 02=+x p x && (1-1-3) 设方程的特解为 st e x = 将上式代入(1-1-3)处特征方程及特征根为 ip s p s ±==+2,1220 则(1-1-3)的通解为 pt D pt C e C e C x ipt ipt sin cos 11+=+=- (1-1-4) C 、 D 为任意积分常数，由运动的初始条件确定，设t=0时 00,x x x x &&== （1-1-5） ()x m x k W F && =+?-= ∑量位静平衡位置 一自由度弹簧—质量系统 ? ==k mg W x &x )
第1章 单自由度系统的振动 1.1概述 机械振动是工程中常见的物理现象。悬挂在弹簧上的物体在外界干扰下所作的往复运动就是最简单直观的机械振动。广泛地说，各种机器设备及其零部件和基础，都可以看成是不同程度的弹性系统。例如桥梁在车辆通过时引起的振动，汽轮机、发电机由于转子不平衡引起的振动等。因此，机械振动就是在一定的条件下，振 动体在其平衡位置附近所作的往复性的机械运动。 实际中的振动系统是很复杂的。为了便于分析研究和运用数学工具进行计算，需要在满足工程要求的条件下，把实际的振动系统简化为力学模型。例如图示1.1-1 就是个最简单的单自由度质量(m )—弹簧(k )系统。 如果实际系统很复杂，要求的精度较高，简化的力学模型也就复杂。 振动系统中和参数的动态特性，可以用常系数线性微分方程来描述的，称为线性振动。但工程实际中也有很多振动系统是不能线性化的，如果勉强线性化，就会使系统的性质改变，所得的系统只能按非线性振动系统处理。 机械振动分析方法很多。对于简单的振动系统，可以直接求解其微分方程的通解。由于计算机进行数值计算非常方便，所以振动仿真是一种最直接的方法。 由于振动模型中尤其是多自由度振动很方便用矩阵微分方程来描述，所以MATLAB 语言在振动仿真中体现出十分优越的特性。 本章先介绍机械振动的单自由度、多自由度振动的基础，然后介绍仿真计算的各种计算公式，最后通过MATLAB 语言来实现。 1.2单自由度系统的振动 1.2.1 无阻尼自由振动 如图1.1-1所示的单自由度振动系统可以用如下微分方程描述 ： 0=+kx x m (1.2.1-1) 令 m k n = 2ω ，方程的通解为 t b t a x n n ωωcos sin += (1.2.1-2) 式(1.2.1-2)表示了图示(1.1-1)中质量m 的位置随时间而变化的函数关系，反映了振动的形式与特点，称为振动函数。 式(1.2.1-2)中，a 、b 为积分常数，它决定于振动的初始条件。如假定t =0时，质量块的位移 x=x 0，其速度 00V x x == ,则 00 ,x b V a n == ω 即 图 1.1-1
第3次作业题： 1、如图所示起重机小车，其质量为m 1=2220kg,在质心A 处用绳悬挂一重物B ，其质量为m 2=2040kg 。绳长l=14m,左侧弹簧是缓冲器，刚度系数k=852.6kN/m 。设绳和弹簧质量均忽略不计，当车连同重物B 以匀速v 0=1m/s 碰上缓冲器后，求小车和重物的运动。 2、两个质量块m 1和m 2用一弹簧k 相连，m 1的上端用绳子拴住，放在一个与水平面成а角的光滑斜面上，如习题下图所示。若t=0时突然割断绳子，两质量块将沿斜面下滑。试求瞬时t 两质量块的位置。 答案： α ωsin ]) (cos 2)([21222 221221g m m k t m t m m k m x +-++= αωsin ]) (cos 2)([21222 221222g m m k t m t m m k m x ++++= 3．如图，已知m 2＝2×m 1=m ，k 3=2k 1=2k 2=2k ，x 10=1.2,x 20=10x ＝20x ＝0，试求系统的固有频率，主振型以及相应。 答案：利用程序，易得 固有频率： ωn 1=3.162277rad/s ，ωn 2=5 rad/s 主振型： m 1 m 2 k 3 k 2 k 1
系统相应： t x 5cos 8.03.1622777t cos 4.01+= t x 5cos 4.03.1622777t cos 4.02-= 4．已知:?? ? ???=11009][m ,[c ]= ??????--11.01.01,][k =??????--905050110,)}({t f =? ?? ???21,激振力频率 ω=3rad/s,试求系统的稳态响应。 答案：利用给定程序，输入给定数据，即获得系统的稳态响应。 第四章 多自由度系统振动 §4-1 多自由度系统运动方程的建立 (引言：问题的提出。)工程中的机械振动问题，有一些可以简化成一个或两个自由度系统的振动问题，因此可以用前面几章中介绍的方法进行分析计算。但是也有很多问题不能采用这种过于简化的力学模型来进行分析。一般来说，各种机器及其零部件的质量和刚度都具有分布的性质，因此理论上都是无限多自由度系统，即为弹性体。但由于机器的结构比较复杂，若都按无限多自由度来处理，在数学上有很大的，甚至目前还无法解决的困难。因此，只好将系统的结构用一些离散的结构来理想化。这样就把弹性体变成数目有限个的离散单元组成的有限多自由度系统。 如前所述，振动系统有多少个自由度就有多少个固有频率和主振型，也就有多少阶主振动，因此弹性体就有无穷多阶主振动。但有意
第5章 两自由度系统的振动 应用单自由度系统的振动理论，可以解决机械振动中的一些问题。但是，工程中有很多实际问题必须简化成两个或两个以上自由度，即多自由度的系统，才能描述其机械振动的主要特征。多自由度系统的振动特性与单自由度系统的振动特性有较大的差别，例如，有多个固有频率、主振型、 主振动和多个共振频率等。本章主要介绍研究两自由度系统机械振动的基本方法。 如图5-1所示。平板代表车身，它的位置可以由质心C 偏离其平衡位置的铅直位移z 及平板的转角θ 来确定。这样，车辆在铅直面内的振动问题就被简化为一个两自由度的系统。 5.1 双质量弹簧系统的自由振动 5.1.1 运动微分方程 图5-2(a)表示两自由度的弹簧质量系统。略去摩擦力及其它阻尼，以它们各自的静平衡位置为坐标x 1、x 2的原点，物体离开其平衡位置的位移用x 1、x 2表示。两物体在水平方向的受力图如图5-2(b)所示，由牛顿第二定律得 ? ? ?=+-=-++00)(2212222212111x k x k x m x k x k k x m (5-1) 这就是两自由度系统的自由振动微分方程。习惯上写成下列形式 ???=+-=-+00212211dx cx x bx ax x (5-2) 显然此时 2 2 1 2 1 2 1,,m k d c m k b m k k a = == += 但对不同的系统， 式(5-2)中各系数的意义并不相同。 图5-2两自由度的弹簧质量系统
5.1.2 固有频率和主振型 根据微分方程的理论，设方程(5-2)的解，即两自由度无阻尼自由振动系统的解为 ?? ? ??+=+=)sin()sin(2211ααpt A x pt A x (5-3) 或写成以下的矩阵形式 )sin(2121α+??? ???????=??????????pt A A x x (5-4) 将式(5-4)代入式(5-2)，可得代数齐次方程组 ? ?? ???=????????????----002122 A A p d c b p a (5-5) 保证式(5-5)具有非零解的充分必要条件是式(5-5)的系数行列式等于零，即 0)(2 22 =----= ?p d c b p a p 展开后为 0)(24=-++-bc ad p d a p (5-6) 式(5-6)唯一确定了频率p 满足的条件，通常称为频率分程或特征方程。它是2p 的二次代数方程，它的两个特征根为 )(222 2 2 ,1bc ad d a d a p --??? ??++= bc d a d a +?? ? ??-+=2 22 (5-7) 由于式(5-7)确定的2p 的两个正实根仅取决于系统本身的物理性质，与运动的初始条件无关，因此p 称为系统的固有频率。较小的一个称为第一阶固有频率，较大的一个称为第二阶固有频率。 5.2.2 主振型 将固有频率p 1和p 2分别代入式(5-5)的任一式，可得到对应于它们的振幅比