基于MATLAB的振动模态分析毕业设计论文文档格式.docx
《基于MATLAB的振动模态分析毕业设计论文文档格式.docx》由会员分享，可在线阅读，更多相关《基于MATLAB的振动模态分析毕业设计论文文档格式.docx（81页珍藏版）》请在冰豆网上搜索。
在国内也有越來越多的科学技术工作者参加到学习和倡导这种语言的行列中來。
应用MATLAB软件对选矿用振动筛的振动特性进行研究，可以充分发挥计算机技术的优势,为选矿用振动筛振动特性研究探索新的途径。
在工程振动中，确定系统固有频率与主振型是非常重要的。
固有频率是决定系统振动特性的重要物理量，它既是防止系统共振的依据，乂是多自由度系统解耦分析（模态分析）的前提，因此研究某系统振动时，首先要求出系统的固有频率。
主振型则为初步分析各构件的振动情况以及解耦分析奠定了基础。
对于多自由度振动系统，计算系统固有频率与主振型主要有2种方法〔9
（1）利用特征矩阵方程式与特征方程式求解；
（2）矩阵迭代法求解。
2种方法各有各的特色。
对于低自由度的振动系统，方法一容易、快捷。
但是在实际工程中，大多数振动
系统都是自由度较多，用特征矩阵方程式与特征方程式求解系统固有频率与主振型这种传缆的计算方法虽然从原则上可行，但当自由度增加时，惯性、刚度阵的阶数增高，计算量也急剧加大，这显然很不方便。
但采用矩阵迭代法，即使是自由度很大的振动系统，计算量也只不过是多进行矩阵迭代而己，而且假设的初始矩阵愈接近实际状况，迭代的次数愈少，相应的计算量也愈少。
1.2国内外研究现状
1.2.1机械振动理论的发展状况及应用现状
振动理论是力学的一个重要组成部分［2】，人类对振动现象的认识有悠久的历史。
振动力学的物理基础在17世纪己经奠定，到了18世纪，振动力学己从物理学中独立出来。
最主要的成就为线性振动理论的形成，它是与数学中的常微分方程和偏微分方程同步发展的。
目前，振动及系统按运动微分方程的形式分为以下两种。
线性振动：
描述其运动的方程为线性微分方程，相应的系统称为线性系统。
线性振动的一个重要特性是线性叠加原理成立。
非线性振动印描述其运动的方程为非绒性微分方程，相应的系统称为非线性系统。
对于非线性振动叠加原理不再成立。
在实际的振动机械或振动系统中，严格的讲，都是非线性的。
但是，建立振动系统的非线性力学模型难度大,求解困难，有些问题其至无解可求。
在实际的工程应用中，很多情况下在误差允许的范围之内用线性的方法解决复杂的近线性问题。
线性振动有确定的力学模型一一线性微分方程，可以求得准确的解，能够描述出振动系统的主要特征。
由于用线性振动的方法能够解决众多的工程实际问题，线性振动的理论一直倍受关注，并且在理论和实验方面己经得到很大的发展和成熟。
特别是多自由度系统的振动的理论，可以说既是振动力学的核心乂是应用得最广泛的振动理论。
线性振动在当今不仅是作为基础科学的力学的一个重要组成部分，而且正走上向工程科学发展的道路，它在航空、机械、船舶、车辆、建筑、水利等工业技术部门中占有愈來愈重要的地位。
线性振动的应用可分为两个方面：
一个方面是减少由于振动而造成的危害，目的在于减振其至于避免有害的振动；
另一个方面利用振动，如工业上常釆用的振动筛选、振动沉桩、振动输送以及按振动理论设计的测量传感器、地震仪等等就是这方面的典型例子。
选矿用振动筛是振动筛选设各中的一种，线性振动理论在选矿用振动筛的设计制造及生产运行中有着广泛的应用，有关这方面的内容将在下一节中详细介绍。
线性振动的理论在发展过程中产生了一个重要分支，那就是模态分析理论。
在对选矿用振动筛进行分析时，需要通过实验來验证理论的正确性，振动实验则需要用到模态分析技术。
模态分析技术从20世纪60年代后期发展至今己趋成熟"
I。
它和有限元分析技术一起，己成为结构动力学中两大支柱。
模态分析是结构动力学中的〜种“逆问题”分析方法，它与传统的“正问题”方法（主要是指有限元方法）不同，是建立在实验（或实测）的基础上，采用实验与理论相结合的方法来处理工程中的振动问题。
目前这一技术己发展成为解决工程中振动问题的重要手段，在机械、航空、航天、土木、建筑、造船、化工等
工程领域被广泛应用［5】。
近十年來，模态分析理论吸取了振动理论、信号处理、信号分析、数据处理、数理统计及自动控制理论中的有关“营养”，结合自身内容的发展，形成了一套独特的理论为模态分析及参数识别技术的发展奠定了理论基础。
模态分析的基础理论概念主要包括；
机械阻抗、导纳、传递函数（或频响函数）、实模态、复模态等。
模态测试技术主要釆用同时测量输入及输出的方法，对一个振动系统來说，可以表示成图1-1所示的
框图
输入（激励）
系
输出（响应）
统
图1.1模态分析框图
Fig.1.1ModalAnalysisDiagram
通过测量激励和响应，进行模念分析可以确定系统。
自从FFT问世以来，目前广泛釆用宽频带激振技术。
其中主要有脉冲、阶跃激励，快速正弦扫描等瞬态激励和纯随机、伪随机、周期随机、瞬态随机等激励方法。
此外，由于F弦慢扫描技术测试精度高，它仍不失为重要激励手段。
模态参数辨识的频域方法有：
分量分析法、导纳圆辨识方法、正交多项式曲线拟合、非线性优化辨识方法等。
模态参数辨识的时域方法与模态参数辨识的频域方法不同，它无需将所测得的响应与激励的时间历程信号变换到频域中去，而是直接在时域中进行参数辨识。
它与频域法相比，两者所采取的分析路线不同，如图1.2所示。
吋域
FFT
步页
传递函
参数
►
数估计
识另IJ
图1.2模态参数辨识分析路线框图
Fig.1.2Modalparameteridentificationoflinediagram
时域法比频域法发展较晚，但近几年來有长足的进展。
自70年代以來主要有:
Ibrahim时域法（简称LTD法）、最小二乘复指数法（LSCE法）、多参考点复指数法（PRCE法）、特征系统实现算法（ERA）。
模态分析技术在动态载荷识别、模型修正与结构动力修改中有广泛的应用，结构动态特征灵敏度分析是非常重要的方法之一。
模态综合技术主要有组合系统法和模态综合法。
随着电子技术与计算机技术的迅速发展，模态分析己成为解决复杂结构振动问题的主要工具，并与计算机辅助设计（CAD）.计算机辅助实验（CAT）相结合,进入产品设计阶段，作为计算机辅助工程中的重要环节，有着广泛的应用〔6】。
1.2.2MATLAB软件的发展状况及应用现状
MATLAB软件概述：
MATLAB的名称源自MatrixLaboiatoiy,是一门计算语言口⑺。
在工程计算领域，计算机技术的应用正逐步将科技人员从繁重的计算工作中解放出来。
在科学计算和工程应用的过程中，一些技术人员尝试用Basic,Foman以及C语言编制程序来减轻计算的工作量，但编制程序不仅要掌握所用语-的语法，还要对有关算法进行深入分析。
为了满足用户对工程数学计算的要求，MATLAB的功能、特点、应用范围：
MATLAB越來越广地被人们应用是源于它在求解方程、数值计算、程序编写上的优点，而它的这些优点是由它的功能和特点决定的。
MATLAB的主要功能丄1）数值计算功能，一条MATLAB语句相当于几十条C语言或Foman语言的语句。
（2）符号计算功能，利用MATLAB的符号计算功能可以清晰地获得解的表达式，对于避免出错和提高程序的可读性均有很大的帮助。
（3）数据分析和可视化功能，在科学计算和研究工作中，技术人员经常会遇到大量的原始数据，而对数据的分析往往难于入手。
MATLAB能将这些数据以图形的方式显示出來，不仅使数据间的关系清晰明了，而且对于揭示其内在本质往往有着非常重要的作用。
MATLAB提供了良好的用户界面，许多函数本身会自动绘制出图形，而且会自动选取坐标刻度，绘制出直角坐标、极坐标、对数坐标下的二维和三维图形，以及条形图、直方图、等高线图、饼形图、离散数据图和瀑布图等专用图形。
（4）文字处理功能。
MATLAB的主要特点：
（1）功能强大，MATLAB不但在数值计算和符号计算方面具有强大的功能，而且在计算结果的分析和数据可视化方面也有着其它类似软件难以匹敌的优势A】。
Notebook,Simulink功能以及各种专业工具箱将MATLAB的应用扩展到非常广的领域。
（2）界面友好、编程效率高，MATLAB的指令表达方式与标准教科书的数学表达式非常相近，用户不需要有较高的计算机编程基础，只要按照计算要求输入表达式，
MATLAB将为用户计算出结果。
同时使用MATLAB语言设计的程序，其编译和执行速度都超过了传统c和Foman语言设计的程序，在工程计算方面的编程效率也高于其它编程语言。
（3）扩展性强，MATLAB的最重要特点之一就是它的可扩展性。
这个特点使得用户能够自由地开发自己的应用程序。
这些年來，许多使用MATLAB的科学家、工程师和技术人员己经开发出相当多的不同领域的应用程序。
MATLAB的应用范围：
MATLAB由主包和各种工具箱组成。
主包是MATLAB的核心，工具箱是扩展的有专门功能的函数。
例如，控制系统工具箱应用于连续和离散系统设计、
频域和时域响应等控制领域；
信号处理工具箱应用于自适应去噪和压缩、谱分析和估计等信号处理领域；
通信工具箱应用于信号编码、调制解调等通信领域。
应用MATLAB的各种工具箱可以在很大程度上减小用户编程时的复杂度，因此MATLAB在很广的领域内得到了应用，其典型应用有；
自动控制、图像信号处理、生物医学工程、语音处理、雷达工程、信号分析、振动理论、时序分析与建模、化学、统计学、经济学等
1.3MATLAB语言的优点
MATLAB作为一个以矩阵和数组为核心计算的软件，对矩阵迭代法中的矩阵迭代计算尤其适合［1°
］。
就所查的资料看，以前的学者和研究人员迭代求解系统固有频率与主振型时,大部分都是用VisiualBasic或Fortran语言來编写程序〔⑴。
限于VisiualBasic或Fortian本身语句以及语法的局限性，用这种高级语言编写的程序涉及到选择合适的算法和编写冗长的语言代码以及键入和调试等一系列问题。
即使有现成的标准予程序可供调用，要在一些较复杂的、科研问题中编写一个完整的程序仍然是一个复杂的、技巧性很强的工作。
因此,用高级语言编写的程序一般代码段较长，需要调用的子程序较多，整个程序的通读性较差。
相反，MATLAB则有简洁、可读性强等优点。
1.4本文研究的内容
振动机械在国民经济中占有重要的位置，振动筛是振动机械中的重要一员。
一直以來有许多人对振动筛进行设计和研究，但是，振动筛的动态设计和计算机辅助设计近年來刚刚起步。
振动特性是振动筛非常重要的有别于非振动机械的一个本质特点，却往往被设计者和制造者简单化。
客观的说，一般的振动都是非线性的，但在许多情况下可以近似看作线性來处理。
线性振动理论不论从基础理论还是实验技术方面近年來都有很大的发展，特别是应用现代化振动测试仪器测量振动信号以及应用计算机软件來分析处理振动信号，为从事振动研究的科技人员带來了极大的方便。
把振动的理论应用到工程实际中去，切实解决工程中遇到的实际的振动问题是研究振动理论的根本目的。
需要对该力学模型进行深入的分析（借助MATLAB软件进行仿真分析）。
本文主要利用MATLAB对振动系统进行模拟分析对于虚拟抽象的理论图像化，处理单自由度振动的3个阻尼和强迫单自由度阻尼振动，多自由度系统振动矩阵迭代求解。
2单自由度系统的振动
2.1单自由度振动系统数学模型的建立问
建立和分析有粘性阻尼时的自由度振动微分方程。
以静平衡位置为原点建立如图坐标，由牛顿定律得运动方程为⑷：
//////
kx
-AA/WV-m
twc^cx+kx=O（2-1）
令
cCrk
2h=—®
=—
mm
其中〃称为衰减系数，单位为1/S；
©
是相应的无阻尼时的固有频率，式（2-1）可以写为：
（2-2）
x+2nx+a）^x=0
如果进一步令
（2-3）
其中无量纲的§
称为相对阻尼系数，则式（2・2）可写为:
（2-4）
x^lqcox+cd;
}x=0
为了求解，令
x=b（2-6）
代入（2-4）后得到特征方程：
52+2qG\s+=0（2-7）
他的两个特征根为：
力=_®
”土乌広-1（2-8）
根据相对阻尼系数g的不同大小，可以将阻尼分为三种状态：
g〉l时为过阻尼，§
=1时为临界阻尼，Ovgvl时为欠阻尼。
1）过阻尼状态
§
>
1,久与比是两个不等的负实根，
初始条件
=©
洁-1
（2-9）
x（Q）=xo,x（O）=xq
（2-10）
系统初始条件响应为
XO=e
严g湖+兀+呼"
asha/t）
CD
（2-11）
临界阻尼状态
g=l,s=-©
是二重根，
方程（2-4）
的通解为系统对式（2-10）的初始条件的响应为
x（t）=e~^[x0+（i0+a）nxQ）t]
（2-12）
欠阻尼状态
<
1,其中
初始条件响应
（2-13）
X0=£
~^nl（x0coscodt+m+Wosincodt）叫
（2-14）
2.2参数设定与求解
阻尼比g分别取；
应用Matlab对式（2-11）和（2-12）,（2-14）求解。
程序如下:
cleai;
fomiatcompact;
a=0.5;
t=0:
0.1:
18;
;
w0=l;
k=l;
xO=l;
wd=wO.*sqit（l-a*a）;
xl=wd
y=exp（-a*wO.*t）.*（xO.*cos（wd.*t）+（（xl+a*wd*xO）./wd）*sin（wd.*t））figure（l）,plot（t,y9lrf）;
lioldon
a=1.0;
w0=1;
wd=l;
xl=wd;
y=exp（-wd.*t）.*（xO十（x1+wd*xO）.*t）;
figure
（1）,plot（t,y/d,）;
holdon
a=2.0;
t=0:
wd=wO*sqrt（a*a-l）;
y=exp（-a*wO.*t）.*（x0.*cosh（wd.*t）+（x1+a*wO*xO）/wO.*sinli（t））;
figure（l）,plot（t,y/v,）;
结论：
图2-2为Matlab计算后给出的响应曲线，从中可以得到一些重要的结论冋：
在0vgvl的情况下，阶跃信号输入时，输出信号为衰减振荡，其振荡角频率（阻尼振荡角频率）为©
，幅值按指数衰减越大，阻尼越大，衰减越快。
1时，振荡系统等同于两个一阶系统串联。
此时虽然不产生振荡，但也需要经过较长时间才能达到稳态。
在一定的§
之下，欠阻尼系统能够更快地达到稳态值；
而过阻尼系统反应迟饨，动作缓慢，所以系统通常设计成欠阻尼系统，§
取值为2
图2-2
算例绘制无阻尼单自由度系统的固有频率和周期随静变形的变化曲线。
固有频率tv”和周期r”
取g=9.81加/“。
可以利用下列MATLAB程序画出％在0~0.5范围内©
和r”的变换曲线.
%Ex2_l7.111
g=9・81；
for1=1:
101t（i）=0.01+（0.5-0.01）*（i-l）/100;
w（i）=（g/t（i））A0.5;
tao（i）=2*pi*（t（i）/g）A0.5;
end
plot（t,w）;
gtext（rw_n,）;
holdon;
plot（t,tao）;
gtext（lT_nl）;
xlabel（'
delta_s_t\）;
title（rExample2.r）;
35
Example2.1
30
25
15
10
5
00.05
0.20.250.30.350.40.450.5
deltast
1
\
w
n
J■
*
f-
f—
■■■
Tn
20
2.3单自由度系统的强迫振动冏
简谐激励是激励形式中最简单的一种，虽然它在实际中存在的场合比较少但掌握系统对于简谐激励的响应的规律，是理解系统对周期激励或更一般形式激励的响应基础。
图所示的弹簧质量系统中，质量块上作用有简谐激振力
其中P。
为激振力幅，Q为激振频率。
以静平衡位置为坐标原点建立图示的坐标系。
从图的受力分析，得到运动微分方程为：
（2-16）
tnx+cx+kx=/?
0siiicot
由常微分方程理论知道，方程（3.2）的通解x由相应的齐次方程的通解勺和非齐次方程的任意特解©
两部分组成，即
（2-17）
当欠阻尼时，式中忑⑴为有阻尼自由振动，它的特点是振动频率为阻尼固有频率,振幅按指数规律衰减，称为瞬态振动或瞬态响应；
©
⑴是一种持续的等幅振动，它是由于简谐激励振力的持续作用而产生的，称之为稳态强迫振动或稳态振动，在间隔充分长时间考虑的振动就是这种稳态振动，而在刚受到外界激励时，系统的响应则是上述两种振动之和。
可见，系统受简谐激励后的响应可以分为两个阶段，一开始的过程称为过渡阶段，经充分长时间后，瞬态响应消失这时进入过渡阶段，经充分长时间后,瞬态响应消失，这是进入稳态阶段。
将方程（2-15）的两端同除以质量〃八并且令
（2-18）
其中§
为相对阻尼系数，©
为相应的无阻尼系统的固有频率，则方程（2J5）成
x+2cxotlx+a）^x=—snictX
m
（2-19）
上述方程特解可以通过x=3sin（曲-0）或者x=Acoscot+BsmcDt来求得，这里介
绍用复数方法求式（2・19）的特解。
先将式（2J9）写为下列的复数形式
x+2cconx+CD^x=—em
（2-20）
其中x是复数设复数形式的特解为
x=Be®
（2-21）
其中B称为复振幅，其意义是包含有相位的振幅。
将式（2・21）代入
（2-20）,解
记久为频
升级会员 