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

实验名称 单自由度系统数值模拟 实验目的、要求 1．熟悉单自由度系统强迫振动特性和求解方法； 2．掌握强迫振动系统的计算机模拟仿真方法。 实验设备及仪器 1. 计算机 2. Matlab软件 3. c语言 实验步骤 1．利用如右图所示的受力分析，得出单自由度系统强迫振动的运动方程。 物体沿水平方向振动，取物体无扰力下的静平衡位置为坐标原点，水平向右为x轴正向，建立如图所示的坐标系。受力情况如图，其激励力为：，其中， 称为激励力的力幅，为常值。 为激励频率，为常值。 根据牛顿第二定律，得到单自由度系统强迫振动的运动方程： 2．对方程进行求解。 令，，， 则原方程可以变形为： 这是一个非齐次二阶常系数微分方程，根据微分方程理论，它的解由两部分组成 其中，代表齐次微分方程的解，简称齐次解，当时，由前面的单自由度阻尼自由振动可得： 其中：，称为衰减振动的固有频率。 ， 由于激励为简谐的，根据微分方程的理论，上述微分方程有如下形式的特解： 其中的可以如下得到： 将 代入到微分方程中， 要想等式成立，等式两边对应的和系数应该相等（和正交） 可以采用如下写法： 从而： 则系统的稳态解： 微分方程的通解为： 其中：， 右端第一项是齐次解，代表衰减的自由振动；第二项是特解，代表与激励力同频率的简谐运动。自由振动，在运动开始后很短的时间内迅速消失，通常可以不加考虑。强迫振动却不因