1、在国内也有越来越多的科学技术工作者参加到学习和倡导这种语言的行列中来。应用 MATLAB软件对选矿用振动筛的振动特性进行研究，可以充分发挥计算机技术的优势，为选矿用振 动筛振动特性研究探索新的途径。在工程振动中，确定系统固有频率与主振型是非常重要的。固有频率是决定系统振动特性的重 要物理量，它既是防止系统共振的依据，又是多自由度系统解耦分析（模态分析）的前提，因此研究 某系统振动时，首先要求出系统的固有频率。主振型则为初步分析各构件的振动情况以及解耦分析奠 定了基础。对于多自由度振动系统，计算系统固有频率与主振型主要有2种方法（1） : （1）利用特征 矩阵方程式与特征方程式求解；（2）矩阵迭
2、代法求解。2种方法各有各的特色。对于低自由度的振 动系统，方法一容易、快捷。但是在实际工程中，大多数振动系统都是自由度较多，用特征矩阵方程式与特征方程式求解系统固有频率与主振型这种传缆的 计算方法虽然从原则上可行，但当自由度增加时，惯性、刚度阵的阶数增高，计算量也急剧加大，这 显然很不方便。但采用矩阵迭代法，即使是自由度很大的振动系统，计算量也只不过是多进行矩阵迭 代而已，而且假设的初始矩阵愈接近实际状况，迭代的次数愈少，相应的计算量也愈少。1.2国内外研究现状1.2.1机械振动理论的发展状况及应用现状振动理论是力学的一个重要组成部分 ，人类对振动现象的认识有悠久的历史。振动力学的物 理基础在
3、17世纪已经奠定，到了 18世纪，振动力学已从物理学中独立出来。最 主要的成就为线性 振动理论的形成，它是与数学中的常微分方程和偏微分方程同步发展的。目前，振动及系统按运动微分方程的形式分为以下两种。线性振动：描述其运动的方程为线性微分方程，相应的系统称为线性系统。线性振动的一个重要 特，性是线性叠加原理成立。非线性振动:描述其运动的方程为非绒性微分方程，相应的系统称为非线 性系统。对于非线性振动叠加原理不再成立。在实际的振动机械 或振动系统中，严格的讲，都是非线 性的。但是，建立振动系统的非线性力学模型难度大，求解困难，有些问题甚至无解可求。在实际的 工程应用中，很多情况下在误差允许的范围
4、之内用线性的方法解决复杂的近线性问题。线性振动有确 定的力学模型一一线性微分方程，可以求得准确的解，能够描述出振动系统的主要特征O由于用线性 振动的方法能够解 决众多的工程实际问题，线性振动的理论一直倍受尖注，并且在理论和实验方面已 经得到很大的发展和成熟。特别是多自由度系统的振动的理论，可以说既是振动力学的核心又是应用 得最广泛的振动理论。线性振动在当今不仅是作为基础科学的力学的一个重要组成部 分，而且正走上 向工程科学发展的道路，它在航空、机械、船舶、车辆、建筑、水利等工 业技术部门中占有愈来愈重 要的地位。线性振动的应用可分为两个方面：一个方面是减少 由于振动而造成的危害，目的在于减振
5、甚至于避免有害的振动；另一个方面利用振动，如 工业上常采用的振动筛选、振动沉桩、振动输送以 及按振动理论设计的测量传感器、地震 仪等等就是这方面的典型例子。选矿用振动筛是振动筛选设各 中的 一种，线性振动理论在选矿用振动筛的设计制造及生产运行中有着广泛的应用，有尖这方面的内容将在下一节中详细介绍。 线性振动的理论在发展过程中产生了一个重要分支，那就是模态分析理论。在对选矿用振动筛进行分 析时，需要通过实验来验证理论的正确性，振动实验则需要用到模态分析技术。模态分析技术从20 世纪60年代后期发展至今已趋成熟。它和有限元分析技术一起，已成为结构动力学中两大支 柱。模态分析是结构动力学中的种 逆问
6、题”分析方法，它与传统的正问题”方法（主要是指有限元方法）不同，是建立在实验（或实测）的基础 上，采用实验与理论相结合的方法来处理工程中的振动问题。目前这一技术已发展成为解决工程中振 动问题的重要手段，在机械、航空、航天、土木、建筑、造船、化工等工程领域被广泛应用5。近十年来，模态分析理论吸取了振动理论、信号处理、信号分析、数据处 理、数理统计及自动控制理论中的有尖 营养”结合自身内容的发展，形成了一套独特的理论为模态分析及参数识别技术的发展奠定了理论基础。模态分析的基础理论概念主要包括； 机械阻抗、导纳、传递函数（或频响函数）、实模态、复模态等。模态测试技术主要采用同时测量输 入及输出的方法
7、，对一个振动系统来说，可以表示成图 1-1所示的框图 系 统图模态分析框图Fig. 1.1 Modal A nalysis Diagram通过测量激励和响应，进行模念分析可以确定系统。自从 FFT问世以来，目前广泛采用宽频带激振技术。其中主要有脉冲、阶跃激励，快速正弦扫描等瞬态激励和纯随机、伪随机、周期 随机、瞬态随机等激励方法。此外，由于 F弦慢扫描技术测试精度高，它仍不失为重要激励手段。模态参数辨识的频域方法有：分量分析法、导纳圆辨识方法、正交多项式曲线拟 合、非线性优化辨识方法等。模态参数辨识的时域方法与模态参数辨识的频域方法不同，它无需将所 测得的响应与激励的时间历程信号变换到频域中去
8、，而是直接在时域中进行参数辨识。它与频域法相 比，两者所采取的分析路线不同，如图 12所示。时域FFT频率传递函传递参数模态信号函数A数估计识另U图1.2模态参数辨识分析路线框图Fig.1.2 Modal parameter identification of line diagram时域法比频域法发展较晚，但近几年来有长足的进展。自 70年代以来主要有:Ibrahim时域法（简称LTD法）、最小二乘复指数法（LSCE法）、多参考点复指数法（PRCE法）、特征系统实现算法（ERA）。模态分析技术在动态载荷识别、模型修正与结构动力 修改中有广泛的应用，结构动态特征灵敏度分析是非常重要的方法之一。
9、模态综合技术主要有组合 系统法和模态综合法。随着电子技术与计算机技术的迅速发展，模态分析已成为解决复杂结构振动问 题的主要工具，并与计算机辅助设计（CAD） 计算机辅助实验（CAT）相结合，进入产品设计阶段，作为计算机辅助工程中的重要环节，有着广泛的应用 （6）122 MATLAB软件的发展状况及应用现状MATLAB软件概述：MATLAB的名称源自Matrix Laboratory，是一门计算语言口。在工程计算领域，计算机技术的应用正逐步将科技人员从繁重的计算工作中解放出来。在科学计算和 工程应用的过程中，一些技术人员尝试用 Basic，Fortran以及C语言编制程序来减轻计算的工作量，但编
10、制程序不仅要掌握所用语的语法，还要对有尖算法进行深入分析。为了 满足用户对工程数学计算的要求，MATLAB的功能、特点、应用范围：MATLAB越来越广地被人们 应用是源于它在求解方程、数值计算、程序编写上的优点，而它的这些优点是由它的功能和特点决定 的。MATLAB的主要功能：（1）数值计算功能，一条MATLAB语句相当于几十条C语言或Fortran 语言的语句。（2）符号计算功能，利用MATLAB的符号计算功能可以清晰地获得解的表达式，对于避 免出错和提高程序的可读性均有很大的帮助。（3）数据分析和可视化功能，在科学计算和研究工作 中，技术人员经常会遇到大量的原始数据，而对数据的分析往往难于
11、入手。MATLAB能将这些数 据以图形的方式显示出来，不仅使数据间的尖系清晰明了，而且对于揭示其内在本质往往有着非常重要的作用。MATLAB提供了良好的用户界面，许多函数本身会自动绘制出图形，而且会自动选取坐标 刻度，绘 制出直角坐标、极坐标、对数坐标下的二维和三维图形，以及条形图、直方图、等高线图、饼形图、 离散数据图和瀑布图等专用图形。 （4）文字处理功能。MATLAB的主要特点：（1）功能强大，MATLAB不但在数值计算和符号计算方面具有强大的功能，而且在计算 结果的分析和数据可视化方面也有着其它类似软件难以匹敌的优势 9。Notebook，Simulink功能以及各种专业工具箱将MAT
12、LAB的应用扩展到非常广的领域。（2）界面友好、编程 效率高，MATLAB的指令表达方式与标准教科书的数学表达式非常相近， 用户不需要有较高的计算机编程基础，只要按照计算要求输入表达式，MATLAB将为用户计算出结果。同时使用MATLAB语言设计的程序，其编译和执行速度都超 过了传统c和Fortran语言设计的程序，在工程计算方面的编程效率也高于其它编程语言。（3）扩 展性强，MATLAB的最重要特点之一就是它的可扩展性。这个特点使得用户能够自由地开发自己的 应用程序。这些年来，许多使用 MATLAB的科学家、工程师和技术人员已经开发出相当多的不同领域的应用程序。 MATLAB的应用范围：MA
13、TLAB由主包和各种工具箱组成。主包是MATLAB的核匚、，工具箱是扩展的有专门功能的函数。例如，控制 系统工具箱应用于连续和离散系统设计、频域和时域响应等控制领域；信号处理工具箱应用于自适应去噪和压缩、谱分析和估计等信号 处理领域；通信工具箱应用于信号编码、调制解调等通信领域。应用 MATLAB的各种工具箱可以在很大程度上减小用户编程时的复杂度，因此 MATLAB在很广的领域内得到了应用，其典型应用有；自动控制、图像信号处理、生物医学工程、语音处理、雷达工程、信 号分析、振动理论、时序分析与建模、化学、统计学、经济学等1.3 MATLAB语言的优点MATLAB作为一个以矩阵和数组为核心计算的
14、软件，对矩阵迭代法中的矩阵迭代计算尤其适合 10。就所查的资料看，以前的学者和研究人员迭代求解系统固有频率与主振型时，大部分都是用 Visiual Basic或Fortran语言来编写程序门。限于Visiual Basic或Fortran本 身语句以及语法的 局限性，用这种高级语言编写的程序涉及到选择合适的算法和编写冗长的语言代码以及键入和调试等 一系列问题。即使有现成的标准予程序可供调用，要在一些较复杂的、科研问题中编写一个完整的程 序仍然是一个复杂的、 技巧性很强的工作。因此，用高级语言编写的程序一般代码段较长，需要调用的子程序较多，整个程序的通读性较差。相反， MATLAB则有简洁、可读
15、性强等优点。1.4本文研究的内容振动机械在国民经济中占有重要的位置，振动筛是振动机械中的重要一员。一直以来有许多人 对振动筛进行设计和研究，但是，振动筛的动态设计和计算机辅助设计近年来刚刚起步。振动特f生是 振动筛非常重要的有别于非振动机械的一个本质特点，却往往被设计者和制造者简单化。客观的说， 一般的振动都是非线性的，但在许多情况下可以近似看 作线，性来处理。线t生振动理论不论从基础理论 还是实验技术方面近年来都有很大的发展，特别是应用现代化振动测试仪器测量振动信号以及应用计 算机软件来分析处理振动信号，为从事振动研究的科技人员带来了极大的方便。把振动的理论应用到工程实际中去，切实解决工程中
16、遇到的实际的振动问题是研究 振动理论的 根本目的。需要对该力学模型进行深入的分析（借助 MATLAB八件进行仿真分析）。本文主要利用MATLAB对振动系统进行模拟分析对于虚拟抽象的理论图像化，处理单自由度振动的3个阻尼和强迫单自由度阻尼振动，多自由度系统振动矩阵迭代求解。2单自由度系统的振动2.1单自由度振动系统数学模型的建立旧建立和分析有粘性阻尼时的自由度振动微分方程。以静平衡位置为原点建立如图坐定律得运动方程为购：2-2）st2-6）2-7）代入（Z4）后得到特征方程：,22S 2 乜.n S .n =0他的两个特征根为：2-8）临界阻尼5 o：:门时为欠阻尼1）过阻尼状态1，S1与s是两
17、个不等的负实根，令$ =n ：2 -1（2-9）初始条件x（0） = V（O） = xo（2-10）系统初始条件响应为x （t）二 e_nt （x） ch 4 a0 sh 4）（2-11）S - 一根据相对阻尼系数的不同大小，可以将阻尼分为三种状态：1时为过阻尼，=1时为临界阻尼状态=1 ,s二5是二重根，方程（2-4）的通解为系统对式 （0）的初始条件的响应为x（t）二 entxo（X o 一 nXo） t（2-12）欠阻尼状态1，其中 *n 7- 2（2-13）初始条件响应2.2参数设定与求解阻尼比 分别取；应用Matlab对式（211）和（212） , （2-14）求解。程序如下：cle
18、ar,format compact;a=0.5;t=0:0.1:18;w0=1;k=1 ;x0=1;wd=w0.*sqrt（1 a*a）;x 1=wdy=exp（-a*w0.*t）.*（x0.*cos（wd.*t）+（x1+a*wd*x0）./wd）*si n（wd.*t）figure（1）,plot（t,y/r,）;hold ona=1.0;w0=1 ;wd=1 ;x1=wd;y=exp（-wd.*t）.*（x0+（x1+wd*x0）.*t）;figure（1 ）,plot（t,y/d）;a=2.0;wd=w0*sqrt（a*a-1）;y=exp（-a*w0.*t）.*（x0.*cosh（w
19、d.*t）+（x1 +a*wO*xO）/wO.*si nh（t）;figure,plot（t,y,V）;结论：图2-2为Matlab计算后给出的响应曲线，从中可以得到一些重要的结论 叫在0 ： 1的情况下，阶跃信号输入时，输出信号为衰减振荡，其振荡角频率（阻尼 振荡角频率）为，幅值按指数衰减越大，阻尼越大，衰减越快。1时，振荡系统等同于两个一阶系统串联。此时虽然不产生振荡，但也需要经过较长时间 才能达到稳态。在一定的之下，欠阻尼系统能够更快地达到稳态值；而过阻尼系统反应迟饨，动作缓慢，所图22算例绘制无阻尼单自由度系统的固有频率和周期随静变形的变化曲线。固有频率-n和周期n取g =9.81 m
20、/s2。可以利用下列MATLAB?序画出r在00.5范围内二和n的变换曲线：%Ex2_17.mg=9.81;for i=1:101 t（i）=0.01 +（0.5-0.01 ）*（i-1 ）/100;w（i）=（g/t（i）F0.5;tao（i）=2*pi*（t（i）/gF0.5;end plot（t,w） ;gtext（w_nhold on ;plot（t,tao）;gtext（*T_ n1）;xlabel（,delta_s_t，）；title（,Example2.1,）；Example2.12.3单自由度系统的强迫振动佝简谐激励是激励形式中最简单的一种，虽然它在实际中存在的场合比较少但掌握
21、 系统对于简谐激励的响应的规律，是理解系统对周期激励或更一般形式激励的响应基础。图所示的弹簧质量系统中，质量块上作用有简谐激振力其中p。为激振力幅，为激振频率。以静平衡位置为坐标原点建立图示的坐标系。 从图的受力分析，得到运动微分方程为：（2-16）mx ex kx = p0 si n t由常微分方程理论知道，方程（3.2）的通解x由相应的齐次方程的通解Xh和非齐次方程的任意特解Xp两部分组成，即（2-17）x（t）=Xh（t） XP（t）当欠阻尼时，式中Xh（t）为有阻尼自由振动，它的特点是振动频率为阻尼固有频率，振幅按 指数规律衰减，称为瞬态振动或瞬态响应； Xp（t）是一种持续的等幅振动
22、，它是由于简谐激励振力的持续作用而产生的，称之为稳态强迫振动或稳态振动，在间隔充分长时间考 虑的振动就是这种稳态振动，而在刚受到外界激励时，系统的响应则是上述两种振动之和o可 见，系统受简谐激励后的响应可以分为两个阶段，一开始的过程称为过渡阶段，经充分长时间 后，瞬态响应消失这时进入过渡阶段，经充分长时间后，瞬态响应消失，这是进入稳态阶段。将方程（215）的两端同除以质量m，并且令-2 （2-18）m其中为相对阻尼系数为相应的无阻尼系统的固有频率，则方程（2-15）成x 2, r）x fx 二旦 sin t m（2-19）上述方程特解可以通过兀二Bsin（, t_）或者AcostBsint来求
23、得，这里 介绍用复数方法求式（2-19）的特解。先将式（2-19）与为下列的复数形式X +2SnX +%二山 6 也（2-20）其中x是复数设复数形式的特解为（2-21）其中B称为复振幅，其意义是包含有相位的振幅。将式（2-21 ）代入（2-20），解（2-22）记为频率比，它定义为则式（222）可以写成PoQ ik1 i2, ki（1 _2）2 （2.）2式中B=j1k（1 -2）2（2 ）22将式（2-24）代入（221 ），得到复数形式的特解为（2-23）（2-24）（2-25）（2-26）x=Bei（t）（2-27）比较方程（2-17）与（218），可知（2-19）中的位移x是（2-20）中复数x的虚部, 因此（225）的虚部就是方程（2-12）的特解，即有x= Bsin（ t - J（2-28）其中B为振幅，为相位差。由式（Z26）、2-23）及（2-24）得出稳态强迫振动有如下的基本特点：1）线性系统对简谐激励的稳态响应是频率等同于激励频率而相位滞后与激振力的 简谐振动；2）稳态响应的振幅及相位差只取决于系统本身的物理性质和激振的频率及力幅， 而与系统本身进入运动的方式无尖无阻尼系统对简谐激励的稳态响应可以从式（2-26）得出。
