振动力学课程设计.docVIP

实用标准文案 PAGE 0 精彩文档 振动力学 课程设计 3k3m4m7mk 3k 3m 4m 7m k 5k x1 x2 x3 4k 年级：工程力学09级02班 姓名：陶昶 学号：20091210220 1.在图示振动系统中，建立系统的振动微分方程，要求写出详细的过程。 m m2 m1 m3 k2 k1 k4 k3 2.求系统的振动固有频率。 3.计算系统的振动模态，绘制主振型的示意图。 4.计算系统的主质量、主刚度和简正振型矩阵。 5.初始条件为：,位移单位为m,速度单位为m/s。求系统自由振动的响应。 6.在质量为1的物体上作用简谐力，求系统强迫振动的响应。 7.在质量为3的物体上作用非周期激励力，为单位阶跃函数，求系统强迫振动的响应。 8.在固定端和第1个物体之间安装一个阻尼系数为的阻尼器，在第1个和第2个物体之间安装一个阻尼系数为的阻尼器，在第2个和第3个物体之间安装一个阻尼系数为的阻尼器，在第3个物体和固定端之间安装一个阻尼系数为的阻尼器。已知：。建立系统的有阻尼振动微分方程，计算系统的阻尼矩阵、模态阻尼矩阵。 9.用瑞利法估算系统的基频。 10.用传递矩阵法计算系统的固有频率。 3k3m 3k 3m 4m 7m k 5k x1 x2 x3 4k 1. 分别以两物体的平衡位置为坐标原点，取 x1 、x2和x3 为广义坐标，由牛顿第二定律得 自由振动微分方程为 写成矩阵形式为 上式即为系统的振动微分方程。 2. 令 则振动微分方程可以写作 令 即 展开得系统的本征方程 用MATLAB求解该方程得 3.广义本征值方程为 1模态矩阵为 1 主振型示意图如下图所示 0.88390.4613 0.8839 0.4613 1 1 -1.6009-0.9703 -1.6009 -0.9703 48.6668 48.6668 1 1 -21.0338 -21.0338 4.用MATLAB求主质量和主刚度 上式即为简正振型矩阵。 5.用MATLAB求得模态矩阵的逆矩阵为 主质量和主刚度为 主坐标为 得到用主坐标表示的动力学方程 将原坐标的初始条件化为主坐标的初始条件 则主坐标表示的系统自由振动规律为 转换为实际坐标表示的系统自由振动规律为 6.系统的动力学方程为 其中 令 代入上述方程后得到 计算与主坐标对应的激励力幅值，得到 列出解耦的主坐标受迫振动方程 其中 解出 其中 转换为原坐标的受迫振动规律 7.作用于系统上的激励力为 变换为主坐标的激励力为 列出主坐标动力学方程 写出各主坐标的响应函数 脉冲响应矩阵为 其中 同理可求得 8.分别以两物体的平衡位置为坐标原点，取 x1 、x2和x3 为广义坐标，由牛顿第二定律得 则自由振动微分方程为 写成矩阵形式为 上式即为有阻尼振动微分方程。阻尼矩阵为 用MATLAB计算模态阻尼矩阵得 用MATLAB计算系统的柔度矩阵和动力矩阵 近似取，带入用动力矩阵表示的瑞利商公式，得到 ，与基频精确值相比，相对误差约为0.01%。 10.系统各个子传递矩阵为 且 且 为解题方便，使状态向量无量纲化，即令 则各个子传递矩阵改写为 其中。 考虑边界条件有，始端，末端。 所以有。 导出。 解该三次方程得 附件（主要MATLAB程序） 1.求频率 clear %清除变量 clc syms v %