振动理论及应用大作业
- 格式:docx
- 大小:174.43 KB
- 文档页数:5
振动理论在固体中的应用振动理论是物理学中一门重要的理论,它研究物体在受到外力作用之后的振动行为。
在固体中,振动理论被广泛应用于多个领域,包括材料科学、机械工程以及地震学等等。
本文将对振动理论在固体中的应用进行探讨,并且分析其对相关领域的影响。
首先,振动理论在材料科学领域中有着重要的应用。
例如,在材料强度研究中,振动理论提供了一种可行的方法来评估材料的强度和稳定性。
通过分析物体在振动过程中的应力分布情况,研究人员可以预测材料是否会出现疲劳、变形或者破裂等问题,从而为材料的设计和选用提供重要参考。
此外,振动理论还被广泛应用于声学领域。
例如,在音响系统设计中,振动理论被用来分析材料的声学性能,以提高音响效果和音质。
其次,振动理论在机械工程领域中也有着广泛的应用。
例如,在机械结构的设计中,振动理论可以帮助工程师预测和解决结构的振动问题。
通过对结构的分析,工程师可以合理安排结构的材料和几何形状,以避免共振和疲劳现象的发生。
此外,振动理论还可以用于机械故障诊断和预测。
通过监测机械设备的振动特征,工程师可以及时发现并解决潜在的故障,从而提高设备的可靠性和寿命。
此外,振动理论在地震学研究中也发挥着重要作用。
地震是地球内部能量的释放,会引起地壳的振动。
利用振动理论,地震学家可以研究地震波的传播规律和地震活动的机制。
通过分析地震波的振动特征,地震学家可以预测地震的发生时间、地点和强度,并提供重要的地震预警信息。
此外,振动理论还可以用于地震建筑和结构物的设计。
通过使用抗震材料和合理的结构设计,可以减缓地震对建筑物的破坏,提高人们的安全性。
综上所述,振动理论在固体中的应用非常广泛,涵盖了材料科学、机械工程以及地震学等领域。
通过对物体振动行为的研究和分析,研究人员可以预测材料的性能、解决结构的振动问题,甚至预测地震的发生。
随着科技的不断发展,振动理论在相关领域的应用还将不断深化和拓展,为科学研究和实际应用提供更多的可能性。
《振动力学》——习题第二章 单自由度系统的自由振动2-1 如图2-1 所示,重物1W 悬挂在刚度为k 的弹簧上并处于静止平衡位置,另一重物2W 从高度为h 处自由下落到1W 上且无弹跳。
试求2W 下降的最大距离和两物体碰撞后的运动规律。
解:222221v gW h W =,gh v 22=动量守恒:122122v gW W v g W +=,gh W W W v 221212+=平衡位置:11kx W =,kW x 11=1221kx W W =+,kW W x 2112+=故:kW x x x 21120=-= ()2121W W kgg W W k n +=+=ω故:tv t x txt x x n nn n nn ωωωωωωsin cos sin cos 12000+-=+-=xx 0x 1x 12平衡位置2-2 一均质等直杆,长为l ,重量为w ,用两根长h 的相同的铅垂线悬挂成水平位置,如图2-2所示。
试写出此杆绕通过重心的铅垂轴做微摆动的振动微分方程,并求出振动固有周期。
解:给杆一个微转角2a=h 2F =mg由动量矩定理:ah a mg a mg Fa M ml I M I 822cos sin 12122-=-≈⋅-====αθαθ其中12cossin ≈≈θααh l ga p ha mg ml n 22222304121==⋅+θθ g h a l ga h l p T n 3π23π2π222===2-3 一半圆薄壁筒,平均半径为R , 置于粗糙平面上做微幅摆动,如图2-3所示。
试求其摆动的固有频率。
图2-3 图2-42-4 如图2-4 所示,一质量m连接在一刚性杆上,杆的质量忽略不计,试求下列情况系统作垂直振动的固有频率:(1)振动过程中杆被约束保持水平位置;(2)杆可以在铅垂平面内微幅转动;(3)比较上述两种情况中哪种的固有频率较高,并说明理由。
图T 2-9 答案图T 2-9解:(1)保持水平位置:m kk n 21+=ω(2)微幅转动:mglllF2112+=mgl1l2xx2xx'mglll2121+=k2k1ml1l2()()()()()()()()()mgk k l l k l k l mgk k l l k l l k l l l k l mg k k l l k l k l l l l k l l mg l mgk l l l k l l l l l l k l l mg l l l l x x k F x x x 2122122212121221221121212221212211211121212122211211121221112111 ++=+-++=+-⋅+++=⎥⎦⎤⎢⎣⎡+-++++=+-+='+=故:()22212121221k l k l k k l l k e++=mk en =ω 2-5 试求图2-5所示系统中均质刚性杆AB 在A 点的等效质量。
振动理论及工程应用_天津大学中国大学mooc课后章节答案期末考试题库2023年1.振动问题属于动力学问题中的第二类问题,即已知主动力求()。
答案:运动2.振动是指物体在平衡位置附近所做的()。
答案:往复运动3.弹簧串联、等效刚度(),弹簧并联,等效刚度()。
答案:减小增加4.在建立单自由度弹簧—质量系统的运动微分方程时,当选择物块的静平衡位置为坐标原点,假设x轴正方向垂直向下,则物块的位移、速度和加速度的正方向如何确定()。
答案:都垂直向下5.质点或质点系的运动相互影响的现象叫做()。
答案:耦联6.激振力与受迫振动的位移相位差为()时,振动系统达到共振状态。
答案:90°7.小车重P在斜面自高度h处滑下与缓冲器相撞,斜面倾角为α,缓冲弹簧刚度系数为k。
如缓冲质量不计,斜面摩擦不计,小车碰撞后,系统的自由振动周期为()。
答案:8.在图示振动系统中,已知重为P的AB杆对O轴的回转半径为ρ,物块重为Q,两个弹簧的刚度系数均为k,当系统静止时,杆处于水平。
则此系统微振动的圆频率为:()答案:9.关于主振型的正交性,下列说法错误的是()答案:零固有圆频率对应的主振型不与系统的其他主振型关于质量矩阵和刚度矩阵正交10.关于主振型矩阵和正则振型矩阵的关系是()。
答案:将主振型矩阵的各列除以其对应主质量矩阵元素的平方根,得到的振型就是正则振型11.关于主振型矩阵和正则振型矩阵下列说法错误的是()。
答案:将主振型矩阵的各列除以其对应主刚度的平方根,得到的振型就是正则振型12.瑞利第一商用()方程求解,瑞利第二商用()方程求解。
答案:作用力位移13.瑞利法估算基频的结果是精确值的(),邓克莱法估算基频的结果是精确值的()答案:上限下限14.子空间迭代法是将()与()结合起来的计算方法,它对自由度数较大系统的前若干阶固有频率及主振型非常有效。
答案:里兹法矩阵迭代法15.一维单元应变位移关系矩阵B为:()答案:16.在杆的纵向振动中,要考虑的边界条件是()答案:位移和轴向力17.以下不属于梁横向振动的近似解法的是()答案:传递矩阵法18.下列哪些是主动控制的特点()。
(整理)《振动⼒学》课程作业.《振动⼒学》2015春节学期作业⼀、⽆阻尼⾃由振动1、如图所⽰,T型结构可绕⽔平轴O作微⼩摆动,已知摆动部分的质量为w,机构绕O轴的时(即机构处于平衡位置时),两弹簧⽆转动惯量为J,两弹簧的弹簧系数均为k,且当=0伸缩,试求该机构的摆动频率。
(答案:ω)2、如图所⽰,长度为L的刚性杆件,在O点铰⽀,⾃由端固定⼀质量为m的⼩球。
在距离铰⽀端a处,由两个刚度系数为k/2的弹簧将刚性杆件⽀持在铅垂⾯内。
求该系统的固有频率。
(忽略刚性杆件和弹簧的质量)(答案:ω)3、如图所⽰,悬臂梁长为L ,截⾯抗弯刚度为EI ,梁的⾃由端有质量为m 的质量块,弹簧刚度为k ,求系统的固有频率。
(答案:ω= )4、如图所⽰,半径为R 的均质半圆柱体,在⽔平⾯内只作滚动⽽不滑动的微摆动,求其固有⾓频率。
(答案:ω=)5、如图所⽰,抗弯刚度为623010(N m )EI =?? 的梁AB ,借弹簧⽀撑于A,B 两点处,弹簧系数均为300(/)k N m = 。
忽略梁的质量,试求位于B 点左边3m 处,重量为1000()W N = 的物块⾃由振动的周期。
(答案:T=0.533s )6、⼀个重W 的⽔箱,借助四根端点嵌固的竖置管柱⽀撑着。
每根柱⼦的长为L,抗弯刚度为EI 。
试求该⽔箱顺⽔平⽅向⾃由振动的周期。
(管柱的质量忽略不计)(答案:2T = )7、《结构动⼒学基础》,第2章课后习题,第1题、第2题、第8题⼆、有阻尼⾃由振动1、如图所⽰,库伦曾⽤下述⽅法测定液体的粘性系数'c :在弹簧上悬挂⼀薄板A ,先测出薄板在空⽓中的振动周期1T ,然后测出在待测粘性系数的液体中的振动周期2T 。
设液体对薄板的阻⼒等于2A 'c v ,其中2A 为薄板的表⾯⾯积,v 为薄板的速度。
如薄板重W ,试有测得的数据1T 和2T ,求出粘性系数'c 。
空⽓对薄板的阻⼒不计。
(答案:'c =)2、物体质量为2kg ,挂在弹簧下端。
机械振动学学习报告摘要:简述了机械振动学的发展历程,振动利用中的若干新工艺理论与技术,振动机械及其相关技术的应用与发展,介绍了振动在人类生活工作中起到了非常重要的作用。
通过对具体实例——单电机振动给料机的计算分析,得出机械振动对机器工作性能的影响。
并介绍了单自由度、多自由度的线性振动系统振动的基本理论和隔振的基本原理。
关键词:机械振动;振动给料机;线性振动系统Abstract:This paper describes the development course of study of mechanical vibration and the utilization of some new technology theory and technology. The vibration has played a very important role in human life and work. By analyzing the practical example-single motor , vibrating feeder calculation and analysis of mechanical vibration machine has influence on the performance. And introduced the single-degree-of-freedom, multi-freedom system vibration of the linear vibration of the basic theory and the basic principle of vibration isolation.Keywords:Mechanical vibration; Vibrates the feeding machine; Linear vibration system第一章绪论1.1振动振动学的发展振动振动学科是20世纪后半期逐渐形成和发展起来的一门新学科。
1. 如图所示,一直角均质细杆,水平部分杆长为l ,质量为 m ,竖直部分杆长为l 2,质量为m 2,细杆可绕直角顶点处的固定轴O 无摩擦地转动,水平杆的未端与刚度系数为k 的弹簧相连,平衡时水平杆处于水平位置.求杆作微小摆动时的周期.解:依题意,由能量法求系统固有频率.系统动能为222121ωωV H J J T +=,其中水平部分杆的转动惯量为231ml J H =,竖直部分杆的转动惯量为2238)2)(2(31ml l m J V ==.即22222222323)3831(21θωω ml ml ml ml T ==+⨯=以平衡位置为原点,计算系统的势能:竖直部分杆的重力势能为)cos 1(2)cos 1(222θθ-=-⋅=mgl lmg U V ; 弹簧与水平部分杆组成的系统是受重力影响的弹性系统.以平衡位置为零势能位置,重力势能与弹性力势能之和相当于由平衡位置处计算变形的弹簧的单独弹性力的势能,故这部分的势能为221kx U S H =+.当杆做微小摆动时有θl x ≈.因此2221θkl U S H =+. 所以2221)cos 1(2θθkl mgl U U U S H V +-=+=+. 由能量守恒定律知 0)(=+U T dt d ,即0)21)c o s 1(223(2222=+-+θθθkl mgl ml dt d 即0sin 2322=++θθθθθθ kl mgl ml 当杆作微小摆动时,0≠θ且θθ≈sin .上式整理得,032=++θθmlklmg系统固有频率ml kl mg n 32+=ω,系统微小摆动周期klmg mlT n n +==2322πωπ. .2.求如图所示的两种情况下的固有频率.解:⑴如左图所示,悬臂梁与弹簧的受力相同,故而可视为两弹性体的串联.当有集中力作用于悬臂梁的悬空端时,其刚度为33l EIk b =,因此整个系统的等效刚度为3333333klEI EIkk lEI k l EIk k k k k b b eq +=+⋅=+⋅= 所以,忽略悬臂梁的质量得,系统固有频率)3(33kl EI m EIkmk eq n +==ω⑵如右图所示,悬臂梁与弹簧的变形相同,故而可视为三弹性体的并联. 当有集中力作用于悬臂梁的悬空端时,其刚度为33lEIk b =,因此整个系统的等效刚度为332l EIk k k k k b eq +=++=. 所以,忽略悬臂梁的质量得,系统固有频率3323mlkl EI m k eqn +==ω3.均质杆AB ,质量为M ,长为l 3,B 端刚性连接一质量为m 的 物体,其大小可略去不计.AB 杆在O 处用铰链连接,并用弹簧刚度系数均为k 的两弹簧加以约束,如图所示.试求系统自由振动的频率.解:依题意,由能量法求系统固有频率. 系统动能为222121mv J T B +=ω. 其中均质杆的转动惯量为2222241432)3(121Ml Ml Ml l M l M J B =+=⎪⎭⎫⎝⎛⋅+=,均质杆的转动角速度为θω =; 集中质量m 的线速度为l l v ωω22=⋅= 即222222221221l m M l m l M T θωω ⎪⎭⎫ ⎝⎛+=+=以平衡位置为原点,计算系统的势能:弹簧与杆组成的系统是受重力影响的弹性系统.以平衡位置为零势能位置,重力势能与弹性力势能之和相当于由平衡位置处计算变形的弹簧的单独弹性力的势能,故势能为()2222212θθkl l k kx U =≈⨯= 由能量守恒定律知 0)(=+U T dt d ,即0)221(222222=++θθθkl ml Ml dt d 0≠θ ,整理得,042=++θθmM k 即系统固有频率mM kn 42+=ω所以系统自由振动的频率mM kf n 42212+==ππω.4.如图所示,质量为2m 的均质圆盘在水平面上可作无滑动的滚动,鼓轮绕轴的转动惯量为I ,忽略绳子的弹性、质量及各轴承间的摩擦力,求此系统的固有频率.解:依题意,由能量法求系统固有频率:系统的势能由两弹簧的弹性势能组成.右侧弹簧的弹性势能为22221x k U =;左侧弹簧的弹性势能为2222112211121)(21x R R k x R R k U ==.2m故系统的势能为221222121)(21x k k R R U U U +=+=系统的动能由小车平移的动能、圆盘平面运动的动能和鼓轮绕轴转动的动能组成. 其中,小车平移的动能为212112121xm v m T ==; 圆盘平面运动的动能为2222222222222224343))(21(21212121x m v m r v r m xm w J v m T ==+=+=鼓轮绕轴转动的动能为2222223321)(2121R x I R v I Iw T ===(3w 为鼓轮转动的角速度).故系统的动能为2222221321214321R x I x m x m T T T T ++=++= 设位移x 的变化规律为)sin(θω+=t A x n ,则有)cos(θω+=t Aw xn n . 因此系统最大势能为2212221max)(21A k k R R U +=; 系统最大动能为2222222221max214321n n n R A I A m A m T ωωω++= 由能量守恒定律知,max maxU T =.整理得,I R m R m R k R k n232222222212222112+++=ω 即系统的固有频率为IR m R m R k R k n 23222222221222211+++=ω5.在图示系统中以系统的平衡位置开始算起 , 盘的中央的位移当作广义坐标.假定盘很薄,并且做纯滚动.求系统的固有频率.解:依题意,由能量法求系统固有频率:系统的势能由两弹簧的弹性势能组成.(此系统是受重力影响的弹性系统.以平衡位置为零势能位置,重力势能与弹性力势能之和相当于由平衡位置处计算变形的弹簧的单独弹性力的势能),故系统的势能为22229)2(22121kx r r x k kx U p p =⋅⋅⋅+=. 系统的动能由圆盘平面运动的动能、鼓轮绕轴转动的动能和重物竖直运动的动能组成.其中圆盘平面运动的动能为2222214343212121xm mv r v mr mv T ==⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛+=; 鼓轮绕轴转动的动能为22222121pp pp r x I r v I T =⎪⎪⎭⎫⎝⎛⋅=; 重物竖直运动的动能为()2223342221x m r r v m v m T p p =⎪⎪⎭⎫ ⎝⎛⋅=⋅=. 故系统的动能为22232121419pp r x I x m T T T T +=++= 由能量守恒定律知0)(=+U T dt d ,即0)2921419(2222=++kx r xI x m dt d p p 整理得,021992=++x r I m kxpp 所以系统的固有频率为22199pp n r I m k+=ω.6.建立图示系统运动的微分方程.以θ作为广义坐标,并假定θ很小,试求系统的固有频率.解:依题意,θ与x 的关系有,L xL x 3443==θ,即43θL x = 对支点取矩有,)(43434343kx x c x M L L x c L kx L x M M O ++-=⋅-⋅-⋅-= 刚杆的转动惯量I 为:22163)23(121mL L m I =⋅= 由动量矩定理知,)(43kx x c x M L M I O++-== θ 将43θL x =代入上式整理得,03333=++++θθθMm k M m c 所以系统的运动微分方程为:03333=++++θθθMm k M m c .其无阻尼固有频率nω为:Mm kn 33+=ω.7.求图示系统微幅扭振的周期.两个摩擦轮可分别绕水平点1O 与2O 转动,互相吻合,不能相对滑动,在图示位置(半径A O 1与B O 2在同一水平线上),弹簧不受力,弹簧系数为1k 与2k , 摩擦轮可看为等厚均质圆盘,质量为1m 与2m .解:依题意,此系统为单自由度系统,取两摩擦轮的转角21,θθ为坐标.由能量法求系统固有频率.两摩擦轮互相吻合,不能滑动,所以2121,θθθθ B A B A r r r r ==.重力势能无变化.故系统势能U 为:21221222211)(21)(21)(21θθθA B A r k k r k r k U +=⋅+⋅=两摩擦轮的转动惯量21,J J 分别为:22221121,21B A r m J r m J ==.所以系统动能T 为:2122122222121222211)(41)21(21)21(212121θθθθθ A B A r m m r m r m J J T +=⋅⋅+⋅⋅=+=由能量守恒定律知 0)(=+U T dt d ,即0))(21)(41(2122121221=+++θθA A r k k r m m dt d 整理得,0)(2121211=+++θθm m k k所以系统固有频率n ω为:2121)(2m m k k n ++=ω.其微幅振动周期n T 为:)(2222121k k m m T nn ++==πωπ.8.轮子可绕水平轴转动,对转轴的转动惯量为0J ,轮缘绕有软绳,下端挂有重是P 的物体,绳与轮缘之间无滑动.在图示位置由水平弹簧维持平衡.半径R 与a 都是已知的.求微幅振动的周期.解:依题意,由能量法求固有频率.选取轮子转动的角速度θ为坐标.此系统是受重力影响的系统,故计算总势能可由平衡位置计算,不计重力势能.系统势能U 为:22221)(21θθka a k U =⋅=重物竖直运动的速度P v 为:θR v P =,所以系统的动能为: 202220221212121θθθ J g PR J v g P T P +⋅=+⋅⋅=由能量守恒定律知 0)(=+U T dt d ,即0)212121(222022=++⋅θθθka J g PR dt d 整理得,0022=++θθJ gPRka所以系统固有频率n ω为:022J gPRka n +=ω. 其微幅振动周期n T 为:20222kaJ gPR T nn +==πωπ9.求图示两个弹簧在点O 的等值弹簧系数,刚杆AB 可以在图示平面内绕点O偏转.解:依题意,此系统为二自由度系统.取刚杆AB 两端竖直向下的方向为广义坐标21,x x ,点O 处竖直向下为0x .则有)(1210x x ba ax x -++=. 假设在点O 处有竖直向下的集中载荷0P ,则对点A 取矩有,22220)(k P b a a x b a x k a P ⋅+=⇒+⋅=⋅; 对点B 取矩有,11110)(k P b a b x b a x k b P ⋅+=⇒+⋅=⋅ 所以221222101020121210)()()(1)(b a k k b k a k P k P b a b b k P b a a a b a b a bx ax x x b a a x x ++=⋅+⋅+⋅+⋅⋅+=++=-++=等值弹簧系数0k 为:2221221222102210000)()()(bk a k b a k k b k a k P b a k k P x P k ++=++==10.求图示系统的运动方程并求临界阻尼系数与有阻尼固有频率.解:依题意,选取刚杆转动的角度θ为运动坐标.阻尼c 处速度为θa ,弹簧k 处位移为θb . 系统的转动惯量为:2ma J =.由动量矩定理得,a a c b b k J ⋅⋅-⋅⋅-=)()(θθθ即θθθ222ca kb ma --= 所以系统的运动方程为:022=++θθθma kb m c 系统无阻尼固有频率为:22makb n =ω. 临界阻尼系数为:km a bmakb m m c n c 22222===ω 阻尼比ξ为:kmb acc c c 2==ξ 系统有阻尼固有频率为:222222222)2(24)2(11mc m ka b am c a km b m a kb kmb acn d -⎪⎭⎫ ⎝⎛=-=⋅-=⋅-=ωξω11.图示系统,设轮子无侧向摆动,且轮子与绳子间无滑动,不计绳子和弹簧的质量,轮子是均质的,半径为R ,质量为M ,重物质量m ,试列出系统微幅振动微分方程,求出其固有频率.解:依题意,由能量法求固有频率.位移坐标x 方向如图所示.轮子与绳间无滑动,其转动角速度w 为:Rxw =.轮子做平面复合运动,其动能M T 为 222243)()21(2121x M R x MR x M T M =⋅⋅+=.所以系统动能T 为222143x m x M T T T m M +=+= 此系统是受重力影响的系统,故计算总势能可由平衡位置计算,不计重力势能.系统势能U 为:222)2(21kx x k U =⋅=由能量守恒定律知0)(=+U T dt d ,即0)22143(222=++kx x m x M dt d 整理得,0238=++x mM k x. 此即所求微幅振动方程. 所以系统固有频率n ω为:mM kn 238+=ω.12.鼓轮:质量M ,对轮心回转半径ρ,在水平面上只滚不滑,大轮半径R ,小轮半径r ,弹簧刚度21,k k ,重物质量为m ,不计轮D 和弹簧质量,且绳索不可伸长.求系统微振动的固有频率.解:依题意,由能量法求固有频率.鼓轮质心所在的位移坐标x 方向如图所示.此系统是受重力影响的系统,故计算总势能可由平衡位置计算,不计重力势能.系统势能U 为:221)(21x k k U +=.(两弹簧并联,等效刚度eq k 为:21k k k eq +=) 鼓轮在水平面上只滚不滑,因而其速度瞬心为鼓轮与水平面的接触点.所以鼓轮滚动的角速度ω为:Rx =ω,绳子的线速度m v 为:x R rR r R v m +=+=ω)(,鼓轮的转动惯量2ρM J C =,鼓轮做平面复合运动.系统动能T 为:2222222)()(2121)(21212121RxM x M xR r R m J x M mv T C m ⋅+++=++=ρω 即222222)()(x RR M r R m T ρ+++= 由能量守恒定律知0)(=+U T dtd,即 0))(212)()((22122222=+++++x k k x RR M r R m dt d ρ 整理得,0)()()(222221=+++++x R M r R m R k k x ρ 所以系统固有频率n ω为:)()()(222221ρω++++=R M r R m R k k n .13.试用m 的坐标与m 2的坐标写出振系的运动微分方程.刚杆AB 的重量可以不计.解:依题意,整个系统的动能为:2221222121)2(2121x m x m x m x m T +=⋅+=整个系统的势能只计算弹性势能即可.故总势能为:21222121222121)2(21)2(21)(21)(21x x k x x k l l x x x k l l x x x k U -+-=⋅-++⋅--=因此系统的质量矩阵M 为:j i ij x x Tm m m ∂∂∂= ⎥⎦⎤⎢⎣⎡=2200M 系统的刚度矩阵K 为:j i ij x x U k k k k k ∂∂∂= ⎥⎦⎤⎢⎣⎡--=25445K 所以系统的运动微分方程为:⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡--+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡0054452002121x x k k k k x x m m即⎭⎬⎫=+-=-+0542045212211kx kx xm kx kx xm14.刚杆本身的质量可略去不计.再设三个质量都只能沿x方向运动。
振动与波动的基本理论与应用振动和波动是物理学中的重要概念,它们在自然界和工程应用中都起着重要的作用。
本文将介绍振动与波动的基本理论,并探讨它们在不同领域的应用。
一、振动的基本理论1.1 振动的定义与特性振动是物体在固定参考点附近的周期性运动。
振动可以由一个简谐运动表示,即一个物体在固定点附近以固定频率和幅度往复运动。
振动的三要素包括振动的频率、幅度和相位。
1.2 振动现象与振动的形成原因振动现象广泛存在于自然界和人工环境中,例如钟摆的摆动、音叉的发声等。
振动的形成原因包括外力的作用、内力的作用以及能量的转换等,不同的振动系统具有不同的形成原因。
1.3 振动的数学描述与分析方法振动可以通过数学方法进行描述和分析。
常见的方法包括使用振动方程、运动学方程和能量守恒等原理,以及应用傅里叶分析等数学工具对振动信号进行分析。
二、波动的基本理论2.1 波动的定义与特性波动是能量传播的过程,它沿介质中的传播路径传递能量而不传递物质。
波动的基本特性包括波长、振幅、频率和传播速度等。
2.2 波动的分类与性质根据波动的传播介质和振动方向,波动可以分为机械波和电磁波。
机械波需要介质进行传播,而电磁波可以在真空中传播。
波动还具有干涉、衍射和折射等特性。
2.3 波动的数学描述与分析方法波动可以使用波动方程进行数学描述和分析。
常见的波动方程有波动方程、亥姆霍兹方程和斯涅耳定律等,利用这些方程可以对波动信号进行定量分析。
三、振动与波动的应用3.1 振动与波动在物理学中的应用振动和波动在物理学中有广泛的应用,例如声波在空气中的传播、地震波的成因及其对地壳的影响等。
振动和波动理论也被广泛应用于材料研究、电磁学和量子力学等领域。
3.2 振动与波动在工程学中的应用振动和波动在工程学中有许多重要应用,例如结构动力学研究中的振动分析、声学设计中的声波传播模拟以及电磁波在通信系统中的传输等。
振动和波动理论的应用为工程领域提供了重要的工具和方法。
如图所示三自由度弹簧质量系统,试编写MA TLAB程序,求系统:
(1)求固有频率、主振型及正则振型;
(2)对初始条件的自由响应;
(3)对外激励的稳态响应。
参数选取:各个质量值;各段刚度值;初始条件;简谐激励
解:
(1)求固有频率、主振型及正则振型
取k1=k2=k3=k4=k,m1=m2=m,m3=m;初始条件为x1=1,x2=1,简谐激励为F1=sin t,F2=sin2t,F3=sin3t
程序:
%定义刚度矩阵和质量矩阵
%syms k m
k=1;m=1;
k1=k;k2=k;k3=k;k4=k;
m1=m;m2=m;m3=m;
k11=k1+k2;k12=-k2;k13=0;
k21=-k2;k22=k2+k3;k23=-k3;
k31=0;k32=-k3;k33=k3;
K=[k11,k12,k13;k21,k22,k23;k31,k32,k33];
M=[m1,0,0;0,m2,0;0,0,m3];
R=inv(M)*K;
%求特征值特征向量
D1=eig(R);
[p2,d]=sort(D1);
[V,Dm]=eig(R);
%求固有频率
fori=1:3
p(i,1)=sqrt(p2(i,1));
end
%求主振型和正则矩阵
fori=1:3
for j=1:3
A1(i,j)=V(i,d(j))/V(1,d(j));
end
end
Ap=A1;
Mp=Ap'*M*Ap;
Kp=Ap'*K*Ap;
fori=1:3
AN(:,i)=Ap(:,i)/sqrt(Mp(i,i));
end
MN=AN'*M*AN;
KN=AN'*K*AN;
p
Ap
AN
运行结果:
固有频率p,主振型矩阵Ap,正则振型矩阵AN如下:p =
0.4450
1.2470
1.8019
Ap =
1.0000 1.0000 1.0000
1.8019 0.4450 -1.2470
2.2470 -0.8019 0.5550
AN =
0.3280 0.7370 0.5910
0.5910 0.3280 -0.7370
0.7370 -0.5910 0.3280
(2)对初始条件的自由响应
程序:
t=[0,100];
y0=[1;0;0;0;1;0];
[T,X]=ode45(@f,t,y0);
figure;
subplot(3,1,1);
plot(T,X(:,1),'b')
title('自由响应1'),grid on
xlabel('时间t'),ylabel('位移x')
subplot(3,1,2);
plot(T,X(:,2),'r');
title('自由响应2'),grid off
xlabel('时间t'),ylabel('位移x')
subplot(3,1,3);
plot(T,X(:,3),'g');
title('自由响应3'),grid off
xlabel('时间t'),ylabel('位移x')
function dy=f(t,y)
k=1;m=1;
k1=k;k2=k;k3=k;k4=k;
m1=m;m2=m;m3=m;
k11=k1+k2;k12=-k2;k13=0;
k21=-k2;k22=k2+k3;k23=-k3;
k31=0;k32=-k3;k33=k3+k4;
dy=zeros(6,1);
dy(1)=y(4);
dy(2)=y(5);
dy(3)=y(6);
dy(4)=-1/m1*(k11*y(1)+k12*y(2)+k13*y(3)); dy(5)=-1/m2*(k21*y(1)+k22*y(2)+k23*y(3)); dy(6)=-1/m2*(k31*y(1)+k32*y(2)+k33*y(3));
结果:
(3)对外激励的稳态响应
%定义刚度矩阵和质量矩阵
%syms k m
k=1;m=1;w=1;
k1=k;k2=k;k3=k;k4=k;
m1=m;m2=m;m3=m;
k11=k1+k2;k12=-k2;k13=0;
k21=-k2;k22=k2+k3;k23=-k3;
k31=0;k32=-k3;k33=k3+k4;
K=[k11,k12,k13;k21,k22,k23;k31,k32,k33]; M=[m1,0,0;0,m2,0;0,0,m3];
%求特征值特征向量
R=inv(M)*K;
D1=eig(R);
[p2,d]=sort(D1);
[V,Dm]=eig(R);
%固有频率
for i=1:3
p(i,1)=sqrt(p2(i,1));
end
%求主振型和正则矩阵
for i=1:3
for j=1:3
A1(i,j)=V(i,d(j))/V(1,d(j));
end
end
Ap=A1;
Mp=Ap'*M*Ap;
Kp=Ap'*K*Ap;
for i=1:3
AN(:,i)=Ap(:,i)/sqrt(Mp(i,i));
end
MN=AN'*M*AN;
KN=AN'*K*AN;
F1=[1;0;0];
F2=[0;2;0];
F3=[0;0;3];
QN1=AN'*F1*sin(w*t);
QN2=AN'*F2*sin(2*w*t);
QN3=AN'*F3*sin(3*w*t);
b1=[1/(p2(1,1)-w^2),0,0;
0,1/(p2(2,1)-w^2),0;
0,0,1/(p2(3,1)-w^2)];
b2=[1/(p2(1,1)-(2*w)^2),0,0;
0,1/(p2(2,1)-(2*w)^2),0;
0,0,1/(p2(3,1)-(2*w)^2)];
b3=[1/(p2(1,1)-(3*w)^2),0,0;
0,1/(p2(2,1)-(3*w)^2),0;
0,0,1/(p2(3,1)-(3*w)^2)]; XN1=b1*QN1;
XN2=b2*QN2;
XN3=b3*QN3;
T=20;
t=1:0.1:T;
X=AN*XN1+AN*XN2+AN*XN3; figure;
subplot(3,1,1);
plot(t,X(1,:),'g');
title('稳态响应1'),grid on xlabel('时间t'),ylabel('位移x') subplot(3,1,2);
plot(t,X(2,:),'r:');
title('稳态响应2'),grid off xlabel('时间t'),ylabel('位移x') subplot(3,1,3);
plot(t,X(3,:),'b-.');
title('稳态响应3'),grid off xlabel('时间t'),ylabel('位移x')
结果:。