东南大学物理实验报告受迫振动.docx
物理实验报告标题:受迫振动的研究实验摘要:振动是自然界中最常见的运动形式之一,由受迫振动引发的共振现象在日常生活和工程技术中极为普遍。共振现象在许多领域有着广泛的应用,例如,众多电声器件需要利用共振原理设计制作。它既有实用价值,也有破坏作用。本实验采用玻耳共振仪定量测定了阻尼振动的振幅比值,绘制了受迫振动的幅频特性和相频特性曲线,并分析了阻尼对振动的影响以及受迫振动的幅频特性和相频特性。另外,实验中利用了频闪法来测定动态的相位差。目录1引言32.实验方法32.1实验原理32.1.1受迫振动32.1.2共振42.1.3阻尼系数的测量52.2实验仪器63实验内容、结果与讨论73.1测定电磁阻尼为0情况下摆轮的振幅与振动周期的对应关系73.2研究摆轮的阻尼振动83.3测定摆轮受迫振动的幅频与相频特性曲线,并求阻尼系数93.4比较不同阻尼的幅频与相频特性曲线134.总结155.参考文献151引言振动是自然界中最常见的运动形式之一,由受迫振动引发的共振现象在日常生活和工程技术中极为普遍。共振现象在许多领域有着广泛的应用,例如为研究物质的微观结构,常采用核共振方法。但是共振现象也有极大的破坏性,减震和防震是工程技术和科学研究的一项重要任务。表征受迫振动性质的是受迫振动的振幅频率特性和相位频率特性(简称幅频和相频特性)。本实验采用玻耳共振仪定量测定了阻尼振动的振幅比值,绘制了受迫振动的幅频特性和相频特性曲线,并分析了阻尼对振动的影响以及受迫振动的幅频特性和相频特性。2.实验方法2.1实验原理2.1.1受迫振动本实验中采用的是玻耳共振仪,其构造如图1所示:图一铜质圆形摆轮系统作受迫振动时它受到三种力的作用:蜗卷弹簧B提供的弹性力矩-k,轴承、空气和电磁阻尼力矩-bddt,电动机偏心系统经卷簧的外夹持端提供的驱动力矩M=M0cost。根据转动定理,有Jd2dt2=-k-bddt+M0cost (1)式中,J为摆轮的转动惯量,M0为驱动力矩的幅值,为驱动力矩的角频率,令02=kJ,2=bJ,m=M0J则式(1)可写为d2dt2+2ddt+02=mcost (2)式中为阻尼系数,0为摆轮系统的固有频率。在小阻尼(2-2)条件下,方程(2)的通解为:=ae-tcos0t+a+bcos(t+)此解为两项之和,可见摆轮的受迫振动分为两个分运动。随着时间推移,阻尼振动项可以衰减到忽略不计。另一项表示与驱动力同频率且振幅为b的振动。可见,虽然刚开始振动比较复杂,但是在不长的时间之后,受迫振动会到达一种稳定的状态,成为一种简谐振动,可以表示成:=bcost+ (3)振幅b和初相位(为受迫振动的角位移与驱动力矩之间的相位差)既与振动系统的性质与阻尼情况有关,也与驱动力的频率和力矩的幅度M0有关,而与振动的初始条件无关(初始条件只是影响达到稳定状态所用的时间)。b与由下述两项决定:b=m(02-2)2+422 (4)=arctan-202-2 (5)2.1.2共振由极值条件b=0可以得出,当驱动力的角频率为=02-22时,受迫振动的振幅达到最大值,产生共振。共振时的共振的角频率r、振幅r和相位差r分别为:r=02-22 r=m202-2 (6)r=arctan(-02-22)由上式可以看出,阻尼系数越小,共振的角频率r越接近于系统的固有频率0,共振振幅r也越大,振动的角位移的相位滞后于驱动力矩的相位越接近于/2。图2和图3给出了不同阻尼系数的条件下受迫振动系统的振幅的频率相应(幅频特性)曲线和相位差的频率响应(相频特性)曲线。图3图2 2.1.3阻尼系数的测量(1)由振动系统作阻尼振动时的振幅比值求阻尼系数摆轮A如果只受到蜗卷弹簧B提供的弹性力矩-k,轴承、空气和电磁阻尼力矩-bddt,阻尼较小(2<02)时,振动系统作阻尼振动,对应的振动方程和方程的解为:d2dt2+2ddt+02=0=ae-tcosat+aa=02-2由于阻尼振动的振幅随时间按指数律衰减,对相隔n个周期的两振幅之比取对数,则有:ln0n=lnae-tae-(t+nT)=nT实际的测量之中,可以以此来算出值。其中,n为阻尼振动的周期数,0为计时开始时振动振幅,n为的n次振动时振幅,T为阻尼振动时周期。(2)由受迫振动系统的幅频特性曲线求阻尼系数(只适合于202时的情况)由幅频特性可以看出,弱阻尼202情况下,共振峰附近01,+020,由(4)和(6)可得:br=202-2(02-2)2+422(-0)2+2当b=r2时,由上式可得:-0±。在幅频特性曲线上可以直接读出b=r2 处对应的两个横坐标+0和-0,从而可得: +-=2 (8)2.2实验仪器本实验使用的玻耳共振仪由共振仪和控制仪两部分组成,并用电缆互联。玻耳共振仪的示意图见前文图1。振动系统由铜质圆形摆轮A与弹簧B构成,弹簧的一端固定在机架支柱上,另一端与摆轮轴相联,在弹簧弹性力作用下,摆轮可绕轴自由往复振动。外激励是由转速十分稳定的可调电机的偏心轴通过连杆 和摆杆 加到振动系统上。当电机匀速转动时,可看作是一种简谐激励。若改变电机转速,就相当于改变激励的周期。电磁阻尼由阻尼线圈 产生,调节线圈电流可以改变电磁铁气隙中磁场,以达到改变阻尼力矩的作用。角度读数盘上方处也装有光电门,与控制电路相连接,可以用来测量强迫力矩的周期。图4共振仪部分的结构如图4所示。左边是振幅显示窗,显示三位数字的摆轮振幅;右边时间显示窗,显示5位数字振动周期,精度为10-3s。“摆轮、强迫力”和“周期选择”开关,分别用来测量摆轮强迫力矩的1次或者10次周期所需的时间。电机转速调节旋钮用来改变强迫力周期,它是通过精确改变电机转速来达到,其精度仅供参考。阻尼选择开关用来改变阻尼线圈直流电位的大小,。电机开关用来控制电机转动,当测量阻尼系数和摆轮固有频率与振幅关系时,电机开关处于断状态。3实验内容、结果与讨论3.1测定电磁阻尼为0情况下摆轮的振幅与振动周期的对应关系实验数据如下:振幅周期T0/s1501.570 1401.5711311.5721191.5731101.5741011.575891.576811.577671.578541.578501.578 对这些数据进行做点画图,得到:-T0的关系曲线作图日期:2014.9.3作者:岳晨涛周期T0振幅 发现振幅为54和50时数据偏离较大,去掉这两个点后经行线性拟合如下:作图日期:2014.9.3作者:岳晨涛周期T0振幅-T0的线性拟合曲线 能够知道振幅和周期的关系式为:T0=-9.80179*10-5*+1.584763.2研究摆轮的阻尼振动测得数据如下:振幅16014613412311310495877972周期T*10/s15.74515.74515.74515.74515.74515.74515.74515.74515.74515.745根据公式:ln0n=lnae-tae-(t+nT)=nT得: 15.745*=0.7985所以: =0.05071 3.3测定摆轮受迫振动的幅频与相频特性曲线,并求阻尼系数选择“阻尼1”时,实验数据:摆轮振幅振动周期*10/s振动周期*1/s相位差原始周期T0T0/T =0/5816.3171.631723.01.57910.96776716.2171.621727.51.57820.97327816.1301.613033.01.57710.97789716.0101.601042.01.57530.983911115.9401.594049.01.57390.987412615.8701.587059.01.57240.990813515.8211.582165.51.57150.993314715.7091.570990.01.57040.999714815.6841.568497.01.57031.001214415.6551.5655107.01.57061.003312615.6101.5610121.01.57241.007312415.5921.5592127.01.57261.008611015.5841.5584132.01.57451.01037815.4461.5446149.01.57711.02116215.3511.5351156.01.57871.02845315.2641.5264160.01.57961.03484315.1461.5146163.51.58051.0435做出这种情况下的幅频特性曲线,并进行拟合:0作图日期:2014.9.3作者:岳晨涛振幅幅频特性曲线拟合图拟合的具体参数:Equationy=y0 + (A/(w*sqrt(PI/2)*exp(-2*(x-xc)/w)2)Adj. R-Square0.98783ValueStandard ErrorAy049.968982.72179Axc0.999233.49071E-4Aw0.026670.00127AA3.167080.20432Asigma0.01334AFWHM0.03141AHeight94.73532能够得到拟合出的振幅与0的关系式为:=49.9690+3.16710.026672e-2(0-0.9992)20.02672所以,拟合曲线的峰值为:r=144.7189则: b=r2=102.3317用b带入关系式中能够解出:+0=1.01374,-0=0.98466用b带入在3.1中得到的关系式,能够解出0=3.9900就能够解出+=4.0448,-=3.9288根据式(8):+-=2能够得出阻尼系数:=0.05800能够看出,用这种方法求得的阻尼系数与在3.2中求得的阻尼系数是基本吻合的,但是仍然存在一些误差。做出这种情况下的相频特性曲线:相位差0作图日期:2014.9.3作者:岳晨涛相频特性曲线图 选择“阻尼2”时,实验数据:摆轮振幅振动周期*10/s振动周期*1/s相位差原始周期T0T0/T =0/5816.3131.6313251.57910.96806716.2081.6208301.57820.97378416.0791.6079381.57650.980510515.9551.5955501.57450.986812515.8391.5839661.57250.992813615.7811.5781751.57140.995813715.7421.5742861.57130.998213915.7181.5718901.57110.999613815.7011.5701961.57121.000712915.6361.56361151.57211.00549615.5361.55361361.57541.01407815.4591.54591461.57711.02026515.3811.53811531.57841.02625415.2811.52811571.57951.03364715.2021.52021611.58021.0394做出这种情况下的幅频特性曲线,并进行拟合:作图日期:2014.9.3作者:岳晨涛0振幅幅频特性曲线拟合图拟合的具体参数:Equationy=y0 + (A/(w*sqrt(PI/2)*exp(-2*(x-xc)/w)2)Adj. R-Square0.99539ValueStandard ErrorAy049.955321.82357Axc0.999372.68947E-4Aw0.027328.51631E-4AA2.98980.14043Asigma0.01366AFWHM0.03216AHeight87.32559能够得到拟合出的振幅与0的关系式为:=49.9553+2.98980.027322e-2(0-0.9994)20.027322能计算出拟合曲线的峰值为:r=137.2728则: b=r2=97.0665用b带入关系式中能够解出:+0 与 -0用计算“阻尼1”情况下的阻尼系数的方法,同样能够得到“阻尼2”情况下的阻尼系数为:=0.06051。能够看出“阻尼2”与“阻尼1”相比,阻尼系数有所增大这种情况下的相频特性曲线如下:相位差0作图日期:2014.9.3作者:岳晨涛相频特性曲线图3.4比较不同阻尼的幅频与相频特性曲线将“阻尼1”和“阻尼2”的幅频特性曲线拟合图画在同一坐标轴下:振幅作图日期:2014.9.3作者:岳晨涛0 比较之后能够发现,“阻尼1” 的幅频特性曲线拟合图和“阻尼2”的相比峰值更高,即在共振点附近的振幅更大。另外,“阻尼1” 的幅频特性曲线拟合图在共振点两侧的振幅变化更大。但是,在远离共振点的区域内,两者的图像时基本重合的。将“阻尼1”和“阻尼2”的相频特性曲线图画在同一坐标轴下:作图日期:2014.9.3作者:岳晨涛0相位差比较后能够发现“阻尼1” 的相频特性曲线图在共振点附近变化较快,更加陡峭,而“阻尼2”更加的平缓。能够得到:阻尼越大,图像越平缓。4.总结本次实验的结果基本是合理的,但其中“阻尼1”阻尼下两种方法计算出的阻尼系数有较大的误差,现分析如下:1. 幅频特性曲线的作图需要借助查阅“测定电磁阻尼为0情况下摆轮的振幅与振动周期的对应关系”这一项实验内容中的-T0的线性拟合曲线。但是在本次实验中的此项实验数据中明显有两个点和其他点不呈线性关系,有较大的偏离,故将这两个点去掉后进行了拟合。可能这样的拟合方式使得拟合结果不准确,对后面的幅频特性曲线的高斯拟合及阻尼系数的计算产生了较大的影响。2. 可能在“测定摆轮受迫振动的幅频与相频特性曲线,并求阻尼系数”此项内容实验的数据测量操作中,因为操作的失误或者实验室中其他因素的影响导致所测数据欠准确。而本次实验的其他内容方面,如不同阻尼的幅频和相频特性曲线比较等,及本都是符合客观事实的。5.参考文献1钱锋、潘人培,大学物理实验(修订版),北京:高等教育出版社,2005年