受迫振动终结报告

  • 格式:docx
  • 大小:300.14 KB
  • 文档页数:8

阻尼振动与受迫振动 实验报告

生物医学工程系 生医8班 鲍小凡 2008013215 一、 实验目的 1. 观测阻尼振动,学习测量振动系统基本参数的方法; 2. 研究受迫振动的幅频特性和相频特性,观察共振现象; 3. 观测不同阻尼对受迫振动的影响。

二、 实验原理 1. 有粘滞阻尼的阻尼振动 弹簧和摆轮组成一振动系统,设摆轮转动惯量为J,粘滞阻尼的阻尼力矩大小定义为角速度dθ/dt与阻尼力矩系数γ的乘积,弹簧劲度系数为k,弹簧的反抗力矩为-kθ。忽略弹簧的等效转动惯量,可得转角θ的运动方程为

220ddJkdtdt

记ω0为无阻尼时自由振动的固有角频率,其值为ω0=k/J ,定义阻尼系数β=γ/(2J),则上式可以化为:

2220ddkdtdt

小阻尼即2200时,阻尼振动运动方程的解为 22

0exp()cosiittt (*)

由上式可知,阻尼振动角频率为220d,阻尼振动周期为2ddT 2. 周期外力矩作用下受迫振动的解 在周期外力矩Mcosωt激励下的运动方程和方程的通解分别为

22cosddJkMtdtdt

22

0expcoscosiimtttt

这可以看作是状态(*)式的阻尼振动和频率同激励源频率的简谐振动的叠加。一般t>>τ后,就有稳态解

cosmtt

稳态解的振幅和相位差分别为

2

2222

0

/4mMJ



 220

2arctan

其中,φ的取值范围为(0,π),反映摆轮振动总是滞后于激励源支座的振动。 3. 电机运动时的受迫振动运动方程和解 弹簧支座的偏转角的一阶近似式可以写成

cosmtt

式中αm是摇杆摆幅。由于弹簧的支座在运动,运动支座是激励源。弹簧总转角为cosmtt。于是在固定坐标系中摆轮转角θ的运动方程为

22cos0mddJktdtdt

也可以写成 22cosmddJkktdtdt

于是得到 202222204mm







由θm的极大值条件0m可知,当外激励角频率2202时,系统发生共振,θm有极大值220m。 引入参数02kJ,称为阻尼比。 于是,我们得到

222

0012mm

020

2arctan1

4. 描述阻尼振动的常用参量

方程中的参量 转动惯量 劲度系数 阻尼力矩系数 固有角频率 外激励角频率 符号 J k γ /okJ ω

下表给出了文献中描述阻尼振动的常用参量及其计算公式。 序号 名称 符号及公式 1 阻尼系数 /(2)J 2 阻尼比 /2/okJ 3 阻尼振动周期 222/doT

4 时间常数 2/1/J 5 品质因数 1/(2)Q 6 对数减缩率 2/2/1dT

三、 实验任务和实验步骤 1. 调整仪器使波耳共振仪处于工作状态。 2. 测量最小阻尼时的阻尼比ζ和固有角频率ω0。 3. 测量其他2种或3种阻尼状态的振幅,并求ζ。 4. 测定受迫振动的幅频特性和相频特性曲线。

四、 实验仪器和实验步骤。 波尔共振仪

1. 打开电源开关,关断电机和闪光灯开关,阻尼开关置于“0”档,光电门H、I可以手动微调,避免和摆轮或者相位差盘接触。手动调整电机偏心轮使有机玻璃转盘F上的0位标志线指示0度,亦即通过连杆E和摇杆M使摆轮处于平衡位置。然后拨动摆轮使偏离平衡位置150至200度,松开手后,检查摆轮的自由摆动情况。正常情况下,震动衰减应该很慢。 2. 开关置于“摆轮”,拨动摆轮使偏离平衡位置150至200度后摆动,由大到小依次读取显示窗中的振幅

值θj;周期选择置于“10”位置,按复位钮启动周期测量,体制时读取数据10dT。并立即再次启动周

期测量,记录每次过程中的10dT的值。 (1)逐差法计算阻尼比ζ; (2)用阻尼比和振动周期Td计算固有角频率ω0。 3. 依照上法测量阻尼(2、3、4)三种阻尼状态的振幅。求出ζ、τ、Q。 4. 开启电机开关,置于“强迫力”,周期选择置于“1”,调节强迫激励周期旋钮以改变电机运动角频率ω,选择2个或3个不同阻尼比(和任务3中一致),测定幅频和相频特性曲线,注意阻尼比较小(“0”和“1”档)时,共振点附近不要测量,以免振幅过大损伤弹簧;每次调节电机状态后,摆轮要经过多次摆动后振幅和周期才能稳定,这时再记录数据。要求每条曲线至少有12个数据点,其中要包括共振点, 即φ=π/2的点。 五、 实验结果记录及作图 1、 测量最小阻尼时的阻尼比ζ和固有角频率ω0。由大到小依次读取的振幅值θ1,θ2,…,θj,…θ50如下数据表所示:

序号 θj lnθj 周期 序号 θj lnθj 周期 Dj=lnθj+25-lnθj 1 152 5.024 15.509 26 124 4.820 15.478 -0.204 2 151 5.017 27 123 4.812 -0.205 3 150 5.011 28 122 4.804 -0.207 4 148 4.997 29 121 4.796 -0.201 5 147 4.990 30 120 4.787 -0.203 6 146 4.984 31 119 4.779 15.461 -0.205 7 145 4.977 32 118 4.771 -0.206 8 144 4.970 33 117 4.762 -0.208 9 143 4.963 34 117 4.762 -0.201 10 141 4.949 35 115 4.745 -0.204 11 140 4.942 15.495 36 115 4.745 -0.197 12 139 4.934 37 114 4.736 -0.198 13 138 4.927 38 113 4.727 -0.200 14 137 4.920 39 112 4.718 -0.202 15 136 4.913 40 111 4.710 -0.203 16 135 4.905 41 110 4.700 15.447 -0.205 17 133 4.890 42 109 4.691 -0.199 18 133 4.890 43 108 4.682 -0.208 19 131 4.875 44 107 4.673 -0.202 20 130 4.868 45 107 4.673 -0.195 21 129 4.860 15.478 46 105 4.654 -0.206 22 128 4.852 26 105 4.654 -0.198 23 127 4.844 27 104 4.644 -0.200 24 126 4.836 28 103 4.635 -0.201 25 125 4.828 29 102 4.625 -0.203

计算: 1225ln()…=123.166,262750ln()…=118.105

𝑏=ln⁡(θ26θ27…θ50)−ln⁡(θ1θ2…θ25)25×25≈−8.098×10−3

∆b≈Sb=1I√∑(Dj−D̅)2I−1≈0.007 由公式20.52(1)b计算出ζ=0.001289,2322223/24/[]()(4/1)bbb0.001114 dT=1.548s,5100.0010.001ddTTss 22/(1)odT=4.059

0222222220023/2222()()()()(1)1ddTTddd

TTT

=0.0026rad/s。

2、 测量其他2种或3种阻尼状态的振幅,求出ζ、τ、Q。 阻尼2 θj lnθj 周期 阻尼2 θj lnθj 周期 Dj=lnθj+5-lnθj 1 139 4.934 1.551 6 84 4.431 1.543 -0.503 2 125 4.828 1.549 7 76 4.331 1.542 -0.497 3 113 4.727 1.548 8 68 4.220 1.541 -0.507 4 101 4.615 1.546 9 62 4.127 1.540 -0.488 5 93 4.533 1.544 10 56 4.025 1.539 -0.508 计算: ln(θ1θ2…θ5)=23.637,ln(θ6θ7…θ10)=21.134

𝑏=ln(θ6θ7…θ10)−ln(θ1θ2…θ5)5×5≈−0.1001

∆b≈Sb=1I√∑(Dj−D̅)2I−1≈0.003 由公式20.52(1)b计算出ζ=0.01594,2322223/24/[]()(4/1)bbb0.000264 dT=1.544s,5100.0010.001ddTTss 22/(1)odT=4.069

0222222220023/2222()()()()(1)1ddTTddd

TTT

=0.0026rad/s。

阻尼3 θj lnθj 周期 阻尼3 θj lnθj 周期 Dj=lnθj+5-lnθj 1 160 5.075 1.553 6 72 4.277 1.542 -0.798 2 136 4.913 1.551 7 58 4.060 1.540 -0.853 3 120 4.787 1.548 8 45 3.807 1.539 -0.980 4 101 4.615 1.546 9 30 3.401 1.539 -1.214 5 87 4.466 1.544 10 15 2.708 1.544 -1.758 计算: ln(θ1θ2…θ5)=23.856,ln(θ6θ7…θ10)=18.253

𝑏=ln(θ6θ7…θ10)−ln(θ1θ2…θ5)5×5≈−0.2241