模态参数辨识方法.ppt

*,/26,IVE,Institute of Vibration Engineering,Northwestern Polytechnical University,China,模态参数辨识的频域方法,张永强 高级工程师,靖江泰斯特电子有限公司,西北工业大学 振动工程研究所,分析分量法,导纳园识别方法,正交多项式曲线拟合,分量分析法,将频响函数分成实部风量和虚部分量进行分析。,基本公式和主模态概念,N,自由度结构系统结构，,p,点激励,l,点响应的实模态频响函数,等效刚度,与测点与激励点有关,主,模态附近频响函数,若模态密度不很大，各阶模态比较远,则其余模态频响函数在该主模态附近较小，且几乎不随频率变化,因此在第,r,阶模态附近可用剩余模态表示频响函数,剩余模态的实部,与虚部,主模态：当趋于某阶模态的固有频率时，该模态起主要作用此时起主要作用的模态成为主导模态，或叫主模态。,实频图与虚频图,剩余模态与频率无关,在实频图和虚频图上相,当于将横坐标平移一距离,此平行线又名剩余柔度线,模态参数的确定,固有频率的确定,实频线与剩余柔度线交点确定,虚频线的峰值确定,峰值较尖，确定容易,剩余柔度尺寸无影响,因此用虚频峰值确定更好,阻尼比 或 的确定,用半功率带宽来确定,结构阻尼系统阻尼比系数,粘性阻尼系统阻尼比系数,模态振型的确定,对 主模态（不含剩余柔度）,测出,L,个测点的 值（,l,1,，,2,，,3,，,，,L,）,单点激振时 为常数，所以上式即为模态振型。,对激励点归一化的振型,模态刚度的确定,取原点频响函数且对原点归一化,模态刚度,模态质量的确定,分量分析法的特点,简单方便，许多信号分析仪有实虚频图分析能力；,当模态密度不高时，有一定的精度；,峰值有误差时，直接影响辨识精度；,模态密集时，用半功率带宽确定模态阻尼，误差较大；,模态密集时剩余模态不能用复常数表示，辨识精度受影响；,图解法受图解精度的影响。,基本方法,在共振峰值附近取,6,10,个频率点；,用它们拟合,Nyquist,导纳圆图；,拟合方法是理论圆上的数值与实测值之间的误差最小；,根据导纳圆特点辨识模态参数；,相邻模态耦合较松忽略剩余模态时多自由度系统第,r,阶模态附近可按单自由度处理，近似公式如下：,导纳圆辨识方法,固有频率,圆上使 具有最大值的点,模态振型,p,点激励,l,点响应的位移导纳虚部的大小为导纳圆直径；,由各测点导纳圆半径组成的归一化向量组成。,模态阻尼,根据导纳圆几何关系即可求出,模态刚度,取原点频响函数,对原点归一化：,求出原点导纳圆半径后由,此时有,模态质量,多自由度多模态情况,导纳圆有多个，不一定是整圆；,根据各段圆弧拟合各导纳圆；,如图有三个模态的固有频率。,导纳圆方法的优缺点,不仅利用频响函数峰值信息；,同时利用附近很多点的信息；,没有峰值信息也可求出固有频率；,此法任建立在主模态基础之上，模态密集时误差较大；,精度受图解精度限制；,许多仪器都能显示导纳圆图。,基本思想,计算机方法,构造多项式传递函数,测得,L,个频率点的频响函数,构造理论值和测量频率点值的误差,使误差值最小求得多项式的系数,得到拟合的频响函数,n,2N,m=n-2,s=,j,w,正交多项式曲线拟合,极点：令分母多项式为零，此时的,s,为极点,由极点可求的固有频率和阻尼比,留数（振型）：,正交多项式曲线拟合,参数辨识的关键间题是如何使频响函数的实测值与理论值之间的总方误差最小，即进行搜索，寻找一个使目标函数最小的最优值。若第,j,个频率点,j,处的频响函数值以,Hj,表示，而实测到的值用 表示，则误差,i,可表示如下：,由于,Hj,的表达式中，其分母含有待辨识的固有频率与模态阻尼比，因而,i,按其本质来说，对待辨识的参数是非线性的。假若我们构造一个误差函数,ej,，使,ej,与待识别参数成线性关系，进而进行最佳搜索，则为线性优化方法，本节所述即为线性优化方法之一。,一、频响函数的有理分式多项式由式,由式 ，传递函数矩阵,H(s,),可表示如下：,阻抗矩阵,Z,（,s,）的阶次为,NxN,，其中第,i,行，第,j,列的元素,Zij,为,可见阻抗矩阵,Z,（,s,）的行式列的展开式中,s,的最高次,n=2N,，可表达如下,式中的,adjZ,（,s,）为,Z,（,s,）的伴随矩阵，即,Z,（,s,）的代数余子式矩阵的转置矩阵。其元素,s,的最高阶次,m,将比,n,少,2.,因此传递函数矩阵中第,l,行，第,p,列的元素为,将上式分子，分母各除以,bn,，且令,则,其中,dn,=1,。,我们用,FFT,分析仪进行频响函数测量，则,s=,jw,。上式可进一步表示为,二、以正交多项式分别对分子、分母进行曲线拟合,若被拟合的测髦频率点数为,L,，引入负频率概念，使,为,-L,、,-1,，,1,，,L,，共,2L,个频率点，且令,对频响函数的实验值以 表示，则,为书写方便，省去角标,lp,，引人理论模型值与测量值之间误差,为使误差对系数,ck,，,dk,线性化，对上式进行变换得,定义总方差为,式中,E=(e-L,、,e-1,，,e1,，,eL)T,，角标,H,表示共轭转置。,可见，误差向量,E,可表示为,式中,为求得,J,的极小值，可使,J,对矢量,C,和,D,求偏导为,0,由式 可得，,式中：,若,P,（,j,）和,q,（,j,）满足下列条件：,则式 中的矩阵,Y,和,Z,将变为单位矩阵,展开上式，分别求得,D,和,C,满足式的条件，实际上就是使,pk,(j,）和,qk,（,j,）,(K=0,1,，,）为正交多项式。从而使式的系数矩阵的对角线为单位矩阵。列向量,C,及,D,能单独求解，且方程的阶次几乎是原式的一半，从而降低了系数矩阵的条件数，使病态问题有所改善。这就是采用正交多项式进行曲线拟合的好处之一，且,C,、,D,分别求解，系数矩阵阶次较原有的低一半左右，因而运算速度明显加快。,三、极点和留数在已知分子、分母多项式系数后，令分母多项式为零，此时的,s,值，即为极点，由极点值可进一步求得固有频率和阻尼比。,P,点激励，,l,点响应的传递函数,Hlp(s,),可表示为：,式中,Alpr,为,Hlp,（,s,）在极点,sr,的留数，在已知极点的情况下，把上式两边各乘以（,s,sr,），并令,s,sr,，则,式中，,D,为,D,（,s,）对,s,的一阶导数，同理，可求得,四、剩余模态的影响上述公式只适用于待识别的模态数与系统的自由度数相一致的情况。实际结构往往并非如此，除了模态稀疏的情况下，仍可应用上述计算步骤外，一般均需考虑待辨识频段以外的模态影响，即所谓剩余模态的影响。,Richardson,曾用增加分子多项式阶数的办法来解决这个问题。若在待辨识的频段内，包含两个模态峰值，分母多项式阶数,n,取为,2N=4,，而分子多项式阶数取为,m=2N+3=7,，当然也可增加分母多项式的阶次，同样可以取得较好的效果。,五、关于总体曲线拟合,在单点频响函数测试时，往往是一点激励，再测量原点及其余各点响应。大多数曲线拟合是在局部”意义上进行的，也就是说在一组测量中，每个测量单独地进行曲线拟合，对每一模态进行模态频率、阻尼和留数的估计，因此，每个模态要估计四个参数（计算复留数作为两个参数），若有五个模态的参数需被辨识，那么，总共有,20,个未知参数需被估计。一般而论，精确地估计阻尼和留数比精确地估计频率更难，而留数义总是和阻尼估计紧紧地联系在一起的。若阻尼有较大误差，势必导致留数估计有较大误差。为了减少误差，可以把曲线拟合过程分为两个步骤：第一步，估计频率和阻尼，不同测量点的同一模态的固有频率和阻尼应是相同的。故借助所有测量数据，先求得固有频率和阻尼；第二步，用已知的频率和阻尼值，进一步分别估计振型或留数。这种估计过程，我们称为总体曲线拟合。,总体曲线拟合的优点是借助所有测量数据这比从单个测量的曲线拟合能得到更精确的频率和阻尼的估计,。再因阻尼是作为第一步拟合的结果，是已知的，所以，通常在第二步估计中能更精确地估计留数。,总体频率和阻尼的直接估计方法可如下进行：因为同一模态，不同测点的频响函数，其表达式的分母均一样，因而可以把所有测点的频响函数进行平均，然后以此平均的频响函数求得总体频率和阻尼。,【,例,】,有一四自由度系统，以正交多项式对原点位移导纳进行曲线拟合，拟合频段选在,19,一,23Hz,之间的密集模态频段内。为说明该方法对具有一定噪声污染的测量数据的效性，为此分别对频响函数的实部及虚部加以,5,、,10,的随机噪声。,表列出固有频率与阻尼比的拟合结果。图为具有,5,的随机噪声的频响函数、经拟合后的理论曲线（如虚线所示）与实际曲线的比较。,由上所述，正交多项式曲线拟合方法即使在相邻模态藕合只差,0.8Hz,，噪声污染达,10,的频响函数情况下进行曲线拟合，也能得到较好的效果，从所给出的数据中可以看出，固有频率与阻尼比几乎相等，而振型稍有差别。,