1-《机械振动基础》大作业,基于matlab的多自由度振动讲义.docVIP

《机械振动基础》大作业 （2016年春季学期） 题目 多自由度振动系统的固有频率和固有阵型 姓 名 学 号 班 级 专 业 报告提交日期 哈尔滨工业大学 报告要求 请根据课堂布置的2道大作业题，任选其一，拒绝雷同和抄袭； 报告最好包含自己的心得、体会或意见、建议等； 报告统一用该模板撰写，字数不少于3000字，上限不限； 正文格式：小四号字体，行距为1.25倍行距； 用A4纸单面打印；左侧装订，1枚钉； 课程报告需同时提交打印稿和电子文档予以存档，电子文档由班长收齐，统一发送至：liuyingxiang868@hit.edu.cn。 此页不得删除。 评语： 成绩（15分）： 教师签名： 年 月 日 基于MATLAB的对多自由度振动系统的固有频率和固有阵型的分析 一、利用MATLAB编程求七自由度系统的固有频率和固有振型 模型如下图所示 系统中各质量为 （kg）：m1=2;m2=4;m3=5;m4=6;m5=5;m6=4;m7=2; 各处弹性系数为 （N/m）k1=5;k2=5;k3=5;k4=5;k5=5;k6=5;k7=5; 二、实验程序如下： clear close m1=2; m2=4; m3=5; m4=6; m5=5; m6=4; m7=2; k1=5; k2=5; k3=5; k4=5; k5=5; k6=5; k7=5; m=[m1,0,0,0,0,0,0; 0,m2,0,0,0,0,0; 0,0,m3,0,0,0,0; 0,0,0,m4,0,0,0; 0,0,0,0,m5,0,0; 0,0,0,0,0,m6,0; 0,0,0,0,0,0,m7]; k=[k1+k2,-k2,0,0,0,0,0; -k2,k2+k3,-k3,0,0,0,0; 0,-k3,k3+k4,-k4,0,0,0; 0,0,-k4,k4+k5,-k5,0,0; 0,0,0,-k5,k5+k6,-k6,0; 0,0,0,0,-k6,k6+k7,-k7; 0,0,0,0,0,-k7,k7]; [V,D]=eig(k,m) for j=1:1:7 w(j)=sqrt(D(j,j)); for i=1:1:7 absV(i,j)=abs(V(i,j)); end end max=(absV); for j=1:1:7 for i=1:1:7 V(i,j)=V(i,j)/max(j); end end w figure x=1:7; for a=1:7 subplot(3,4,a),plot(x,V(x,a)); hold on; grid on; title(振型图); end 将程序在MATLAB上运行，得到结果如下 图（1）多自由度系统的固有频率 得到振型图如下 图（2） 多自由度系统的固有阵型 图（2）中分别为本实验条件下，模拟出来小车的振动的阵型图。 四、心得、体会 1、对多自由度系统的认识和学习 多自由度系统，是指必须通过两个以上的独立广义坐标才能够描述系统运动特性的系统，或者说是自由度数目多于一个，但又不属于连续弹性体（自由度数目为无穷多个）的系统。在实际工程问题中，大多都是多自由振动的系统，由此可见，我们的大作业非常有实际的应用价值。 在学习本章节的时候，就有一个问题出现，多自由度的系统离我们很遥远，至少实在数学计算中，一个复杂的两自由度系统就需要用公式推导一段时间，才能得出相应的参数（运动的微分方程和相关的固有频率），期间需要列出功和能的方程、求导以及化简等步骤。到后来学习到多自由度系统时，就会发现，开式我们对多自由度系统的学习也只是没有找到方法而已，在引入新的方法后，替代了原有的直接法列方程，对我们有了很大的帮助，就拿刚度、柔度计算发来说就它们的规律而言，就要比直接法快的很多。所谓直接法，就是直接应用动力学的基本定律或定理（例如牛顿第二定律或达朗伯原理）建立系统运动微分方程的方法。当我们使用刚度系数法，解决简单的多自由度系统就很简单了，刚度系数法又称单位位移法，是把一个动力学系统当作一个静力系统来看待，用静力学方法确定出系统所有的刚度系数（刚度矩阵中元素），借助于这些系数建立系统的运动方程。对于本题来说，分别给每一个质量块一个单位的位移，之后再去分析，用隔离体分析法建立的运动方程，就能将矩阵中的参数求解出来。推广来说的话，用刚度系数法同样也可以建