中国振动工程学会模态分析高级研修班讲课资料(第一章).pptx

Institute of Vibration Engineering 振动工程研究所振动工程研究所 1姜节胜西北工业大学 振动(zhndng)工程研究所模态分析理论(lln)基础第一页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 2 模态分析理论基础是20世纪30年代机械阻抗与导纳的概念上发展起来(q li)。吸取了振动理论、信号分析、数据处理、数理统计、自动控制理论的有关营养，形成一套独特的理论。模态分析的最终目标是识别出系统(xtng)的模态参数，为结构系统(xtng)的振动分析、振动故障诊断和预报、结构动力特性的优化设计提供依据。第二页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 3 解析模态分析可用有限元计算实现，而试验模态分析则是对结构进行可测可控的动力学激励，由激振力和响应的信号求得系统的频响函数矩阵，再在频域或转到时域采用多种识别方法求出模态参数，得到结构固有的动态特性(txng)，这些特性(txng)包括固有频率、振型和阻尼比等。模态分析定义为：将线性时不变系统振动微分方程组中的物理模态分析定义为：将线性时不变系统振动微分方程组中的物理(wl)(wl)坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，坐标变换的变换矩阵为振型矩阵，其及模态参数描述的独立方程，坐标变换的变换矩阵为振型矩阵，其每列即为各阶振型。每列即为各阶振型。第三页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 4 有限元分析软件(如ANSYS、NASTRAN、SAP、MAC等)在结构设计中被普遍采用,但在设计中,由于计算模型和实际结构的误差,而且受到边界条件很难准确确定的影响,特别是结构的形状和动态特性很复杂时，有限元简化模型和计算的误差较大。通过对结构进行实验模态分析，可以正确确定其动态特性，并利用动态实验结果(ji gu)修改有限元模型，从而保证了在结构响应、寿命预计、可靠性分析、振动与噪声控制分析与预估以及优化设计时获得有效而正确的结果(ji gu)。第四页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 5a.a.获得结构的固有频率获得结构的固有频率,可避免共振现象的发生可避免共振现象的发生 当外界激励力的频率等于振动系统的固有频率时，系当外界激励力的频率等于振动系统的固有频率时，系统发生共振现象。此时系统最大限度地从外界吸收能量，统发生共振现象。此时系统最大限度地从外界吸收能量，导致结构过大有害振动。结构设计人员要设法使结构不工导致结构过大有害振动。结构设计人员要设法使结构不工作在固有频率环境中。作在固有频率环境中。相反，共振现象并非总是有害的：振动筛、粉末碾磨相反，共振现象并非总是有害的：振动筛、粉末碾磨机、打夯机和灭虫声发射装置等等就是机、打夯机和灭虫声发射装置等等就是(jish)(jish)共振现象共振现象的利用。结构设计人员此时要设法使这种器械工作在固有的利用。结构设计人员此时要设法使这种器械工作在固有频率环境中，可以获得最大能量利用率。频率环境中，可以获得最大能量利用率。试验试验(shyn)模态分析的典型应模态分析的典型应用用第五页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 6b.b.为了应用模态叠加法求结构响应，确定动强度为了应用模态叠加法求结构响应，确定动强度(qingd),(qingd),和疲劳寿命和疲劳寿命 分析告诉我们任何线性结构在已知外激励作用下他的响分析告诉我们任何线性结构在已知外激励作用下他的响应是可以通过每个模态的响应迭加而成的。所以模态分析另应是可以通过每个模态的响应迭加而成的。所以模态分析另一主要的应用是建立结构动态响应的预测模型，为结构的动一主要的应用是建立结构动态响应的预测模型，为结构的动强度强度(qingd)(qingd)设计及疲劳寿命的估计服务。设计及疲劳寿命的估计服务。c.c.载荷载荷(zi h)(zi h)(外激励外激励)识别识别 由激励和模态参数预测响应的问题称为动力学正问由激励和模态参数预测响应的问题称为动力学正问题，反之由响应和模态参数求激励称为反问题。原则上题，反之由响应和模态参数求激励称为反问题。原则上只要全部的各阶模态参数都求得只要全部的各阶模态参数都求得,由响应就可以求出外由响应就可以求出外激励（称为载荷激励（称为载荷(zi h)(zi h)识别）。识别）。第六页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 7d.d.振动与噪声控制振动与噪声控制 既然结构振动是各阶振型响应的迭加，只要设法控制相关既然结构振动是各阶振型响应的迭加，只要设法控制相关频率附近的优势模态（改设计和加阻尼材料等或使用智能材料）频率附近的优势模态（改设计和加阻尼材料等或使用智能材料）就可以达到控制结构振动的目的。就可以达到控制结构振动的目的。对汽车对汽车(qch)(qch)车厢内或室内辐射噪声的控制，道理也一车厢内或室内辐射噪声的控制，道理也一样。车厢座舱或室内辐射噪声与其结构的振动特性（模态）关样。车厢座舱或室内辐射噪声与其结构的振动特性（模态）关系密切，由于辐射噪声是由结构振动系密切，由于辐射噪声是由结构振动“辐射辐射”出来的。控制了出来的。控制了结构的振动，也就是实现了辐射噪声的控制。结构的振动，也就是实现了辐射噪声的控制。第七页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 8e.e.为结构为结构(jigu)(jigu)动力学优化设计提供目标函数或约束动力学优化设计提供目标函数或约束条件条件 动力学设计，即对主要承受动载荷而动特性又至关重动力学设计，即对主要承受动载荷而动特性又至关重要的结构要的结构(jigu)(jigu)，以动态特性指标作为设计准则，对结，以动态特性指标作为设计准则，对结构构(jigu)(jigu)进行优化设计。它既可在常规静力设计的结构进行优化设计。它既可在常规静力设计的结构(jigu)(jigu)上，运用优化技术，对结构上，运用优化技术，对结构(jigu)(jigu)的元件进行的元件进行结构结构(jigu)(jigu)动力修改；也可从满足结构动力修改；也可从满足结构(jigu)(jigu)动态性动态性能指标出发，综合考虑其它因素来确定结构能指标出发，综合考虑其它因素来确定结构(jigu)(jigu)的形的形状，乃至结构状，乃至结构(jigu)(jigu)的拓扑（布局设计、开孔、增删元的拓扑（布局设计、开孔、增删元件）。动力学优化设计就是在结构件）。动力学优化设计就是在结构(jigu)(jigu)总体设计阶段总体设计阶段就应对结构就应对结构(jigu)(jigu)的模态参数提出要求，避免事后修补的模态参数提出要求，避免事后修补影响全局。影响全局。f.f.有限元模性修正与确认有限元模性修正与确认 当今工程结构计算采用最广泛的计算模型就是有限元模型。当今工程结构计算采用最广泛的计算模型就是有限元模型。再好的算法和软件都是建立在理想的结构物理参数和边界条件再好的算法和软件都是建立在理想的结构物理参数和边界条件假设上的。结构有限元计算结果和试验往往存在假设上的。结构有限元计算结果和试验往往存在(cnzi)(cnzi)不小不小差距。此时在模态试验可信的前提下，一般是以试验结果来对差距。此时在模态试验可信的前提下，一般是以试验结果来对有限元模型进行修正和确认。经过修正和确认的有限元模型具有限元模型进行修正和确认。经过修正和确认的有限元模型具有优化概念下的与试验结果最大的接近。可以进一步用于后继有优化概念下的与试验结果最大的接近。可以进一步用于后继的响应、载荷和强度计算。的响应、载荷和强度计算。第八页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 9单自由度系统(xtng)频响函数分析粘性阻尼系统阻尼力（与振动速度成正比）：强迫振动方程(fngchng)及其解解的形式（s为复数）及拉氏变换：第九页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 10自由(zyu)振动阻尼比阻尼比范围范围(fnwi)(fnwi)（0 01 1）内为欠阻尼内为欠阻尼无阻尼无阻尼(zn)(zn)固有频率固有频率实部：衰减因子，反映系统阻尼实部：衰减因子，反映系统阻尼虚部：有阻尼系统的固有频率虚部：有阻尼系统的固有频率 第十页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 11第十一页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 12结构阻尼(zn)（滞后阻尼(zn)）系统阻尼(zn)力：与位移成正比，相位比位移超前90度结构阻尼(zn)系数运动方程及拉氏变换g g 为结构为结构(jigu)(jigu)阻尼比或结构阻尼比或结构(jigu)(jigu)损耗因子损耗因子第十二页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 13传递函数和频率响应(pn l xin yn)函数用实部和虚部表示(biosh)与粘性(zhn xn)阻尼系统相比频响函数形式相同和 相互置换即可得各自表达式（1jg）k 复刚度第十三页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 14位移、速度和加速度传递函数(hnsh)位移、速度和加速度频率响应函数(hnsh)三者之间的关系动刚度（位移阻抗）动柔度（位移导纳）第十四页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 15质量(zhling)阻抗、阻尼阻抗、刚度阻抗（位移、速度、加速度）质量(zhling)导纳、阻尼导纳、刚度导纳（位移、速度、加速度）左至右 阻抗(zkng)除 ，导纳乘第十五页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 16单自由度频响函数的特性(txng)曲线Bode图(幅频图和相频图)幅频图:频响函数(hnsh)的幅值与频率的关系相频图:相位(xingwi)与频率的关系第十六页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 17阻尼愈大，在固有频率附近(fjn)相位曲线的陡度越小 时曲线始点约为1/k，为弹簧的导纳线；低频(dpn)时外力主要由弹簧力来平衡；时，产生共振，幅值为此时惯性力与弹簧力平衡，激励力与阻尼力平衡 时幅值下降，最后趋向于渐近线 极值为0，高频时系统激励力主要由惯性力来平衡 第十七页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 18实频图（结构阻尼和粘滞阻尼）两个(lin)极值点半功率带宽半功率(gngl)带宽反映阻尼大小阻尼越大，半功率(gngl)带宽越大，反之亦然第十八页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 19虚频图以结构阻尼为例：系统(xtng)共振时虚部达到最大值系统(xtng)共振时实部为零半功率点处的值半功率的概念是针对功率（而非幅值）而言，在半功率点处，虚部正好(zhngho)为其最大值的一半，但幅值却为最大幅值的有效值。第十九页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 20Nyquist图频响函数(hnsh)矢端轨迹图结构阻尼系统Nyquist圆（导纳圆）特点起始点(频率为零)非原点，约在(1/k,-g/k)处，圆心坐标（0，-1/2kg）初相角为arctan(-g)圆的直径为虚部最大值1/(kg)半径为实部最大值1/(2kg)直径处对应(duyng)半功率带宽两个频率点共振频率点第二十页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 21粘滞阻尼(zn)系统Nyquist图特点桃子桃子(to zi)形，阻尼比越小形，阻尼比越小轨迹圆越大轨迹圆越大在固有频率附近，曲线在固有频率附近，曲线接近圆，仍可利用接近圆，仍可利用(lyng)圆圆的特性的特性第二十一页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 22速度与加速度频响函数特性曲线关系(gun x)回顾幅频图第二十二页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 23实频图与虚频图第二十三页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 24Nyquist图第二十四页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 25不同激励下频响函数的表达式要点频响函数反映系统输入输出之间的关系(gun x)表示系统的固有特性线性范围内它与激励的型式与大小无关在不同类型激励力的作用下其表达形式常不相同简谐激励激励力响应位移频响函数第二十五页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 26 周期激励非正弦(zhngxin)周期力，如方波、锯齿波，周期为T响应的傅氏展开频响函数（定义为各频率点上的值）均包含幅值与相位(xingwi)两个量第二十六页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 27瞬态激励(jl)一般瞬态输入傅氏变换相应输出傅氏变换相应频响函数单位脉冲激励(jl)频响函数第二十七页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 28随机(su j)激励输入自相关函数输入自功率谱密度输入输出互相关函数互功率谱密度函数频响函数第二十八页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 29多自由度系统的频响函数分析两类系统约束系统自由系统约束系统2自由度运动方程(fngchng)（无阻尼）傅氏变换第二十九页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 30频响函数矩阵原点频响函数第i点的响应与第i点的激励(jl)之间的频响函数跨点频响函数第i点的响应与第j点的激励(jl)之间的频响函数原点频响函数特性原点频响函数第三十页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 31曲线及特性曲线及特性两个共振频率点（对应于分母为零）两个共振频率点（对应于分母为零）一个反共振点（分子为零）一个反共振点（分子为零）反共振是局部现象（仅仅反共振是局部现象（仅仅 振幅为零，因为此时振幅为零，因为此时(c sh)频响函数的其他频响函数的其他项均不为零）。项均不为零）。机架线第三十一页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 32一般多自由度约束(yush)系统N自由度约束(yush)系统有N个共振频率,（N1）个反共振频率对原点函数共振反共振交替出现对跨点频响函数无此规律一般两个距离远的跨点出现反共振的机会比较近的跨点少机架线第三十二页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 33自由(zyu)系统两自由(zyu)度系统运动方程（无阻尼）频响函数矩阵第三十三页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 34曲线及特性 时系统(xtng)产生刚体运动零频为刚体模态反共振点一个共振点高频时以高阶质量线为渐进线，趋向于零零阶等效零阶等效(dn xio)质量质量机架线第三十四页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 35一般多自由度系统(xtng)频响函数曲线一般总结共振于反共振频率满足以下关系（如果有零频则算第一阶）机架线第三十五页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 36多自由度系统(xtng)模态分析与模态参数 （基本理论及方法）比例阻尼线定常系统(xtng)物理坐标下的运动方程M M、C C、K K均为均为NNNN矩阵矩阵方程包含物理方程包含物理(wl)(wl)坐标坐标为耦合方程为耦合方程第三十六页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 37 传递函数(hnsh)和频响函数(hnsh)矩阵 拉氏变换第三十七页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 38 模态坐标下的运动(yndng)方程-任意(rny)l点的响应为各阶模态响应的线性组合振型矩阵振型矩阵(j zhn)（模态矩阵（模态矩阵(j zhn)）第第r阶振型阶振型（模态向量（模态向量)模态坐标模态坐标-模态坐标下的运动方程第三十八页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 39无阻尼(zn)自由振动特征方程全部模态第r阶模态模态正交性主模态：各阶模态主空间：各阶模态向量所组成的空间主坐标：相应的模态坐标第r阶模态的惯性力对第s阶模态位移所做的功为零；或第r阶模态的弹性力对第s阶模态位移(wiy)所做的功为零第三十九页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 40 模态质量(zhling)和模态刚度-模态刚度(n d)特定归一化情况(qngkung)（模态质量归一)它们的具体值没有太大的意义，取决于振型归一化，这是因为振型只是振动形态，没有振幅的意思。这三个振型（模态向量）是等价的-模态质量第四十页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 410000 解偶后的运动(yndng)方程第四十一页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 42-比例阻尼(zn)系统运动(yndng)方程比例(bl)阻尼模态阻尼M、K对称，所以C也对称，也具有正交性第四十二页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 43解偶运动方程(fngchng)（模态坐标下）对第r阶模态模态频率、模态向量、模态质量模态频率、模态向量、模态质量模态刚度、模态阻尼模态刚度、模态阻尼(zn)总称模态参数总称模态参数第四十三页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 44多自由度系统实模态分析实模态条件各点振动相位差为零，或为180度与无阻尼和比例阻尼系统等价(dngji)实模态下响应模态坐标物理坐标测点l的响应第四十四页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 45单点激励频响函数(hnsh)单点p激励l点响应测量l点与激励点p之间的频响函数(hnsh)频响函数(hnsh)与激励力大小无关第四十五页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 46几个概念(ginin)等效刚度等效质量等效质量与等效刚度的关系等效(dn xio)刚度与测点与激励点有关第四十六页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 47计及刚体位移下的频响函数(hnsh)刚体运动的频响函数(hnsh)考虑刚体位移下的频响函数(hnsh)第四十七页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 48剩余柔度认为是与频率(pnl)无关的常数也可认为是频率(pnl)的线性函数模态截断频响函数的合成频响函数为单个模态之叠加模态截断只关心(gunxn)前几阶和十几阶模态忽略高阶模态的影响所截模态数一般大于被分析模态数的两倍频响函数第四十八页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 49多自由度系统复模态分析特点(tdin)各点相位差不一定是0度或180度（与实模态不同）振型系数为复数结构阻尼系统结构阻尼材料内部阻尼滑移阻尼（接头、螺钉、铆钉、衬垫等）运动方程及拉氏变换R为结构(jigu)阻尼矩阵第四十九页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 50传递函数和频响函数矩阵一般表达(biod)特征解的正交性G为结构为结构(jigu)损耗因子矩阵损耗因子矩阵 GKR (IjG)K为复刚度矩阵为复刚度矩阵 模态矩阵(j zhn)（振型矩阵(j zhn)）模态质量矩阵(j zhn)、模态刚度矩阵(j zhn)都是复数第五十页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 51振型向量正则(zhn z)化频响函数矩阵具体表达第五十一页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 52一般粘性阻尼系统运动方程状态向量和状态方程阻尼矩阵不能在N维主空间(kngjin)解偶，需采用状态空间(kngjin)法引入状态向量状态方程扩展(kuzhn)为2N空间第五十二页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 53自由振动(zhndng)特征方程方程特征值(2N个）特征向量2N维空间系统的复模态频率(pnl)和复振型向量，共轭成对第五十三页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 54模态正交性矩阵(j zhn)表示第五十四页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 55正交性矩阵(j zhn)表达第五十五页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 56模态坐标(zubio)下的解利用正交性解偶后的方程振型叠加解模态坐标(zubio)t0时的模态坐标(zubio)向量第五十六页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 57l点的瞬时(shn sh)位移 第五十七页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 58复模态特性(txng)复共轭特性(txng)特征值与特征向量均为复数，共轭成对，共2N个复模态的正交性复特征向量在2N维空间中正交；而实模态在N维空间中正交复模态解偶性系统运动方程在2N维状态空间解偶，而实模态在N维空间解偶；复模态运动特征系统各点有无规律的相位差，而实模态则为0或180度各点不同时通过平衡点，而实模态则同时通过平衡位置各点的振动频率和周期仍相同，由 决定，对一定模态它是常数系统振动无一定振型，节点也不是固定的，而作周期性移动，这与实模态截然不同自由振动时衰减振动，各点衰减率相同，由 决定，这点与实模态相同第五十八页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 59复模态传递函数表达式模态坐标下运动方程(fngchng)传递函数矩阵频率响应矩阵第五十九页，共60页。Institute of Vibration Engineering 振动工程研究所振动工程研究所 60第六十页，共60页。