机械振动学实验指导书.docVIP

实验名称 单自由度系统数值模拟 实验目的、要求 1．熟悉单自由度系统强迫振动特性和求解方法； 2．掌握强迫振动系统的计算机模拟仿真方法。 实验设备及仪器 1. 计算机 2. Matlab软件 3. c语言 实验步骤 1．利用如右图所示的受力分析，得出单自由度系统强迫振动的运动方程。 物体沿水平方向振动，取物体无扰力下的静平衡位置为坐标原点，水平向右为x轴正向，建立如图所示的坐标系。受力情况如图，其激励力为：，其中， 称为激励力的力幅，为常值。 为激励频率，为常值。 根据牛顿第二定律，得到单自由度系统强迫振动的运动方程： 2．对方程进行求解。 令，，， 则原方程可以变形为： 这是一个非齐次二阶常系数微分方程，根据微分方程理论，它的解由两部分组成 其中，代表齐次微分方程的解，简称齐次解，当时，由前面的单自由度阻尼自由振动可得： 其中：，称为衰减振动的固有频率。 ， 由于激励为简谐的，根据微分方程的理论，上述微分方程有如下形式的特解： 其中的可以如下得到： 将 代入到微分方程中， 要想等式成立，等式两边对应的和系数应该相等（和正交） 可以采用如下写法： 从而： 则系统的稳态解： 微分方程的通解为： 其中：， 右端第一项是齐次解，代表衰减的自由振动；第二项是特解，代表与激励力同频率的简谐运动。自由振动，在运动开始后很短的时间内迅速消失，通常可以不加考虑。强迫振动却不因阻尼而衰减，它的振幅与相角也与运动的初始条件无关，对于一定的振动系统，与是激励力的幅值和激励频率的函数，只要和保持不变，与是常值。强迫振动是稳态振动，通常称为稳态响应，特解也称为稳态解。在自由振动消失后，代表物体的全部运动。 令，称为系统的放大因子 3. 利用计算机程序MATLAB或者C语言编写程序对单自由度系统强迫振动的特性进行模拟仿真。 4．取时，绘制单自由度系统放大因子M与频率比关系的响应曲线。 四、实验结果与分析： 1.程序如下： clear r=0:0.1:4; for x=0.1:0.1:1; a=sqrt((1-r.^2).^2+(2*x*r).^2); M=a.\1; plot(r,M),hold on; end title(单自由度系统放大因子M与频率比r的关系图0701014230 赵建鑫); xlabel(r轴); ylabel(M轴); grid on axis([0 4 0 5]); text(1.03,4.5,\leftarrow0.1); text(1,2.5,\leftarrow0.2); text(1.083,1.5,\leftarrow0.3); text(1,1,\leftarrow0.5); text(1,0.732,\leftarrow0.7); text(1,0.5,\leftarrow1.0); 2.图形如下： 五、实验心得 通过本次上机实验我熟悉了单自由度系统强迫振动特性和求解方法；掌握了强迫振动系统的计算机模拟仿真方法。 实验名称 两自由度系统数值模拟 实验目的、要求 1．熟悉两自由度系统的分析和求解方法； 2．掌握两自由度系统计算机求解和表达方法。 实验设备及仪器 1. 计算机 2. Matlab软件 3. c语言 实验步骤 利用图中所示的受力分析，得出双自由度系统强迫振动的运动方程。 原理：如图所示，由质量m1和弹簧k1组成的系统是主系统。为了在相当宽的工作速度范围内，使主系统的振动能够减小到要求的强度，设计了由质量m2和弹簧k2组成的系统，这两个系统组成两自由度系统。其运动方程为： + + = 2. 求解方程。 解此方程得： 式中 = 因而有： 式中 引入下列符号： 如果，可以得到无阻尼吸振器的主系统振幅的表达式： 3. 利用计算机程序MATLAB或者C语言编写程序对双自由度系统强迫振动的特性进行模拟仿真。 4．取，=1时，时；以及时，绘制两自由度系统X1/X2的绝对值与频率比关系的响应曲线。 四、实验结果与分析 clear r=0.6:0.01:1.3; X11=(1-r.^2).^2; X01=((1-r.^2).^2-0.05*(r.^2)).^2; X12=(1-r.^2).^2+0.04*(r.^2); X02=((1-r.^2).^2-0.05*(r.^2)).^2+0.04*(r.^2).*((1-1.05*(r.^2)).^2); X13=(1-r.^2).^2+0.4096*(r.^2); X03=((1-r.^2).^2-0.05*(r.^2)).^2+0.4096*(r.^2).*((1-1.05*(r.^2)).^2); X14=1; X04=(1-1.05*(r.^2)).^2;