第一章-模态分析理论基础

IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China0/26姜节胜西北工业大学振动工程研究所第一章模态分析理论基础IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China模态分析理论基础是20世纪30年代机械阻抗与导纳的概念上发展起来。吸取了振动理论、信号分析、数据处理、数理统计、自动控制理论的有关营养，形成一套独特的理论。模态分析的最终目标是识别出系统的模态参数，为结构系统的振动分析、振动故障诊断和预报、结构动力特性的优化设计提供依据。IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China解析模态分析可用有限元计算实现，而试验模态分析则是对结构进行可测可控的动力学激励，由激振力和响应的信号求得系统的频响函数矩阵，再在频域或转到时域采用多种识别方法求出模态参数，得到结构固有的动态特性，这些特性包括固有频率、振型和阻尼比等。模态分析定义为：将线性时不变系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，坐标变换的变换矩阵为振型矩阵，其每列即为各阶振型。IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China有限元分析软件(如ANSYS、NASTRAN、SAP、MAC等)在结构设计中被普遍采用,但在设计中,由于计算模型和实际结构的误差,而且受到边界条件很难准确确定的影响,特别是结构的形状和动态特性很复杂时，有限元简化模型和计算的误差较大。通过对结构进行实验模态分析，可以正确确定其动态特性，并利用动态实验结果修改有限元模型，从而保证了在结构响应、寿命预计、可靠性分析、振动与噪声控制分析与预估以及优化设计时获得有效而正确的结果。IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,Chinaa.获得结构的固有频率,可避免共振现象的发生当外界激励力的频率等于振动系统的固有频率时，系统发生共振现象。此时系统最大限度地从外界吸收能量，导致结构过大有害振动。结构设计人员要设法使结构不工作在固有频率环境中。相反，共振现象并非总是有害的：振动筛、粉末碾磨机、打夯机和灭虫声发射装置等等就是共振现象的利用。结构设计人员此时要设法使这种器械工作在固有频率环境中，可以获得最大能量利用率。试验模态分析的典型应用IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,Chinab.为了应用模态叠加法求结构响应，确定动强度,和疲劳寿命分析告诉我们任何线性结构在已知外激励作用下他的响应是可以通过每个模态的响应迭加而成的。所以模态分析另一主要的应用是建立结构动态响应的预测模型，为结构的动强度设计及疲劳寿命的估计服务。c.载荷(外激励)识别由激励和模态参数预测响应的问题称为动力学正问题，反之由响应和模态参数求激励称为反问题。原则上只要全部的各阶模态参数都求得,由响应就可以求出外激励（称为载荷识别）。IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,Chinad.振动与噪声控制既然结构振动是各阶振型响应的迭加，只要设法控制相关频率附近的优势模态（改设计和加阻尼材料等或使用智能材料）就可以达到控制结构振动的目的。对汽车车厢内或室内辐射噪声的控制，道理也一样。车厢座舱或室内辐射噪声与其结构的振动特性（模态）关系密切，由于辐射噪声是由结构振动“辐射”出来的。控制了结构的振动，也就是实现了辐射噪声的控制。IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,Chinae.为结构动力学优化设计提供目标函数或约束条件动力学设计，即对主要承受动载荷而动特性又至关重要的结构，以动态特性指标作为设计准则，对结构进行优化设计。它既可在常规静力设计的结构上，运用优化技术，对结构的元件进行结构动力修改；也可从满足结构动态性能指标出发，综合考虑其它因素来确定结构的形状，乃至结构的拓扑（布局设计、开孔、增删元件）。动力学优化设计就是在结构总体设计阶段就应对结构的模态参数提出要求，避免事后修补影响全局。f.有限元模性修正与确认当今工程结构计算采用最广泛的计算模型就是有限元模型。再好的算法和软件都是建立在理想的结构物理参数和边界条件假设上的。结构有限元计算结果和试验往往存在不小差距。此时在模态试验可信的前提下，一般是以试验结果来对有限元模型进行修正和确认。经过修正和确认的有限元模型具有优化概念下的与试验结果最大的接近。可以进一步用于后继的响应、载荷和强度计算。IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China单自由度系统频响函数分析粘性阻尼系统•阻尼力（与振动速度成正比）：•强迫振动方程及其解•解的形式（s为复数）及拉氏变换：xcfdfkxxcxm...stXex)()()(2sfsxkcsmsIVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China0...kxxcxm02kcsms2002,11jsmk/0)2/(kmc阻尼比范围（0－1）内为欠阻尼无阻尼固有频率实部：衰减因子，反映系统阻尼虚部：有阻尼系统的固有频率自由振动IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,ChinaIVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China结构阻尼（滞后阻尼）系统•阻尼力：与位移成正比，相位比位移超前90度•结构阻尼系数•运动方程及拉氏变换jxfdgkg—为结构阻尼比或结构损耗因子fjxkxxm..)()(])1([2sfsxkjgmsIVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China•传递函数和频率响应函数kjgmssH)1(1)(2kjgmH)1(1)(2–用实部和虚部表示2222222)1()1(11)(ggjgkHg2与粘性阻尼系统相比频响函数形式相同和相互置换即可得各自表达式（1＋jg）k—复刚度IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China位移、速度和加速度传递函数•位移、速度和加速度频率响应函数•三者之间的关系•动刚度（位移阻抗）•动柔度（位移导纳）)()()(sfsxsHd)()()(sfsvsHv)()()(sfsasHa)()()(fxHd)()()(fvHv)()()(faHa)()()()()(22ddvaHHjHjHkcsmssZ2)(kcsmssH21)(IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China•质量阻抗、阻尼阻抗、刚度阻抗（位移、速度、加速度）•质量导纳、阻尼导纳、刚度导纳（位移、速度、加速度）左至右阻抗除，导纳乘jjIVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China•单自由度频响函数的特性曲线Bode图(幅频图和相频图)•幅频图:频响函数的幅值与频率的关系22)()()(IRHHH222211k•相频图:相位与频率的关系)1(2arctan)()(arctan2RIHHHIVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China阻尼愈大，在固有频率附近相位曲线的陡度越小时曲线始点约为1/k，为弹簧的导纳线；低频时外力主要由弹簧力来平衡；时，产生共振，幅值为此时惯性力与弹簧力平衡，激励力与阻尼力平衡时幅值下降，最后趋向于渐近线极值为0，高频时系统激励力主要由惯性力来平衡01,0k2/10m2/1IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China•实频图（结构阻尼）（粘滞阻尼）–两个极值点–半功率带宽2222)2()1(11)(kHR)1(41)(1212,12,1kHR2222)1(11)(gkHRkggkgHggR21})]1(1{[)1(1)(2111222,12,1212g半功率带宽反映阻尼大小阻尼越大，半功率带宽越大，反之亦然IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China•虚频图•（结构阻尼）•（粘性阻尼）•以结构阻尼为例：–系统共振时虚部达到最大值–系统共振时实部为零–半功率点处的值])1[()(222gkgHI])2()1[(2)(222kHI半功率的概念是针对功率（而非幅值）而言，在半功率点处，虚部正好为其最大值的一半，但幅值却为最大幅值的有效值。maxmax2,1)(707.0)(21)(HHHIVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,ChinaNyquist图－频响函数矢端轨迹图•结构阻尼系统Nyquist圆（导纳圆）》特点•起始点(频率为零)非原点，约在(1/k,-g/k)处，圆心坐标（0，-1/2kg）•初相角为arctan(-g)•圆的直径为虚部最大值1/(kg)•半径为实部最大值1/(2kg)•直径处对应半功率带宽两个频率点22222121))(()]([kgkgHHIR共振频率点max1ddsIVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China•粘滞阻尼系统–Nyquist图»特点»»(是变的，所以不是圆)22224141))(()]([kkHHIR桃子形，阻尼比越小轨迹圆越大在固有频率附近，曲线接近圆，仍可利用圆的特性IVEInstituteofVibrationEngineering,NorthwesternPolytechnicalUniversity,China速度与加速度频响函数特性曲线•关系回顾•幅频图)