模态参数辨识方法docx.docx
《模态参数辨识方法docx.docx》由会员分享,可在线阅读,更多相关《模态参数辨识方法docx.docx(17页珍藏版)》请在冰豆网上搜索。
模态参数辨识方法docx
模态参数辨识的频域方法
张永强高级工程师
靖江泰斯特电子有限公司
西北工业大学振动工程研究所
•分析分量法
•导纳园识别方法
•正交多项式曲线拟合
•分量分析法
•将频响函数分成实部风量和虚部分量进行分析。
-基本公式和主模态概念
・N自由度结构系统结构,p点激励1点响应的实模态频响函数
•主模态:
当趋于某阶模态的固有频率时,该模态起主要作用此时起主要作用的模态成为主导模态,或叫主模态。
•主模态附近频响函数
-若模态密度不很大,各阶模态比较远
-则其余模态频响函数在该主模态附近较小,且几乎不随频率变化
-因此在第r阶模态附近可用剩余模态表示频响函数
1研
©2)2+g;7(l—©2)2+g;
+(用+用)
©(莎)=
H阿二
心丄(1-研r+g;
L—g「
剩余模态的实部
与虚部
实频图与虚频图
・剩余模态与频率无关
・在实频图和虚频图上相
・当于将横坐标平移一距离
•此平行线又名剩余柔度线
二模态参数的确定2
・固有频率的确定
-实频线与剩余柔度线交点确定
-虚频线的峰值确定
-峰值较尖,确定容易
-剩余柔度尺寸无影响
-因此用虚频峰值确定更好
•阻尼比©或S的确定
・用半功率带宽来确定
A0.二◎—©
・结构阻尼系统阻尼比系数
gr=a或gr=
©
•粘性阻尼系统阻尼比系数
~T
或'二啦旦
2a)『
-模态振型的确定
•对®=1主模态(不含剩余柔度)
K&.
•测出L个测点的H佔*值(1=1,2,3,L)
•单点激振时―羞为常数,所以上式即为模态振型。
•对激励点归一化的振型
侃(丽1)}广-亡{呱訂
-模态刚度的确定
・取原点频响函数且对原点归一化
H1(co=1)=-
"K&Krgr
1
•模态刚度
K「1
-模态质量的确定M严电
-分量分析法的特点
・简单方便,许多信号分析仪有实虚频图分析能力;
・当模态密度不高时,有一定的精度;
:
・峰值有误差时,直接影响辨识精度;:
・模态密集时,用半功率带宽确定模态阻尼,误差较大;
•模态密集时剩余模态不能用复常数表示,辨识精度受影响;
・图解法受图解精度的影响。
导纳圆辨识方法
•基本方法
•在共振峰值附近取6—10个频率点;
•用它们拟合Nyquist导纳圆图;
・拟合方法是理论圆上的数值与实测值之间的误差最小;
・根据导纳圆特点辨识模态参数;
・相邻模态耦合较松忽略剩余模态时多自由度系统第『阶模态附近可按单自由度处理,近似公式如下:
-固有频率
圆上使攵仏⑵具有最大值的点
co-1时,—max
dco
-模态振型
p点激励1点响应的位移导纳虚部的大小为导纳圆直径;由各测点导纳圆半径组成的归一化向量组成。
K。
{%}广2恢讣
rSr
lP
心1,2,3,…,厶
-模态阻尼
・根据导纳圆儿何关系即可求
co^-cox2
£UL
r®tan乞+如乞
22
-模态刚度
・取原点频响函数
・对原点归一化:
伟广1
•求岀原点导纳圆半径后由
Ker二一
导纳圆求模态阻尼
2(j
此时有
1/K=(pcpIK=l/K
erfpr^prr、
K二K
rer
-模态质量Mr=Krl^.
-多自由度多模态情况
•导纳圆有多个,不一定是整圆;
•根据各段圆弧拟合各导纳圆;
•如图有三个模态的固有频率。
-导纳圆方法的优缺点
•不仅利用频响函数峰值信息;
•同时利用附近很多点的信息;
•没有峰值信息也可求出固有频率;
•此法任建立在主模态基础之上,模态密集时误差较大;
•精度受图解精度限制;
•许多仪器都能显示导纳圆图。
正交多项式曲线拟合
-基本思想
•计算机方法
•构造多项式传递函数
•测得L个频率点的频响函数
•构造理论值和测量频率点值的误差
•使误差值最小求得多项式的系数
•得到拟合的频响函数
*=N(s)
D©口:
极点:
令分母多项式为零,此时的S为极点由极点可求的固有频率和阻尼比
N兀($)二》r=l
留数(振型):
A-N⑸lprD(s)
工帀曲扌以
参数辨识的关键间题是如何使频响函数的实测值与理论值之间的总方误差最小,即进行搜索,寻找一个使目标函数最小的最优值。
若第j个频率点创处的频响函数值以Hj表示,而实测到的值用康示,则滾差疥可表示如下:
由于Hj的表达式中,其分母含有待辨识的固有频率与模态阻尼比,因而密按其本质来说,对待辨识的参数是非线性的。
假若我们构造一个误差函数“使ej与待识别参数成线性关系,进而进行最佳搜索,则为线社优化方法,本节所述即为线性优化方法之一。
一、频响函数的有理分式多项式由式
由式H(s)-Z-[(s)=(52M十虑十KP1递函数矩阵H⑸可表示如下:
Hd)=[孔)]-〜謝翳|①
阻抗矩阵Z(s)的阶次为NxN,其中第i行,第j列的元素Zij为
Ztj(s)=171屛+CtJS+kij
可见阻抗矩阵Z(s)的行式列的展开式中s的最高次n=2N,可表达如下
det[Z(5)J=十6js十妇"十…+b/71二D(s)
①式中的adjZ(s)为Z(s)的伴随矩阵,即Z(s)的代数余子式矩阵的转置矩阵。
其元素s的最高阶次m将比n少2•因此传递函数矩阵中第1行,第p列的元素为
N($)_uq++a2$十刃
D($)&Q+6]54-力2百十…+b„sn
将上式分子,分母各除以bn,且令
加(、)=I■〃I(”)=W几(Q-昇;q(心)=1,生(£)二♦…,q“(C-宀
则
匕^'kPk(、)
Hlp(s)=丄
其中dn=1o
我们用FFT分析仪进行频响函数测量,则s=jwo上式可进一步表示为
\CK加(je)
二△•
二、以正交多项式分别对分子、分母进行曲线拟合
若被拟合的测髦频率点数为L,引入负频率概念,使CD为CD-L....
CD-I,col,・・・coL,共2L—个频率点,且令
3T=—3八
H%(j3J=He(・T)=砒(仙〉
1=1*2.…
对频响函数的实验值以表示,则
舗p(T)胡八-沁)=用;(込)
i二1■2•…丄
为书写方便,省去角标lp,引人理论模型值与测量值之间误差
e,-H(T)-H(j®j)=
Wg(js)
导_片(冋)
•为使误差对系数Ck,dk线性化,对上式进行变换得
H
K・山
>:
“你(3〉LH(jtt)•)[L〃KqK(je,)+sdT)」
K0KD
•定义总方差为
J=二EhE
i"-t
•式中E=(e-L>...e-1,el,...eL)T,角标H表示共轨转置。
•可见,误差向量E可表示为
E=PC-QD-W
式中
A(j®i)Qo(jsJ
L)
•
…灯(j3・・L〉
…pw(jaj_|)
/>0(j3[)
■
•
P\(js)
…几(T)
•
■
■
•••PmSi.)
r)H(j^
…-£
L)宵(j3
-]>7iOi)
…H(ja_|
p-、
…片(jwj偽
月(问匸)他(血丄)
)知[(闷・jL
…虫血』务・](瓯)「
+=[宵(闷“)%(j®l),…,片(j®-i)q”(j9--i)
H(ja>[)q兀(jaJ•…月(T.)S(jeL>]P
•为求得J的极小值,可使J对矢量C和D求偏导为0昙"1=0.
X=-畑pHQ};
・式中:
Y=y(PHP十FTP*);z=y(ef,ae'e*);G=Re(PllW);
F=Re(QHB0o
•若P(jco)和q(jco)满足下列条件:
II,s=t
b・$Hf,(a.z=0,1.…
X|H(j^)|,g;(j©心(皿〉=工丨片(j仙)|Zg/jQg;(问)二r=wt
•则式雷:
监H;|中的矩阵Y和Z将变为单位矩阵
展开上式,分别求得D和C
(I-XtX)D=F-XtG,
C=G-XDo
满足③式的条件,实际上就是使pk(j(D)和qk(jco)(K=0,1,…)为正交多项式。
从而使式②的系数矩阵的对角线为单位矩阵。
列向量C及D能单独求解,且方程的阶次几乎是原式②的一半,从而降低了系数矩阵的条件数,使病态问题有所改善。
这就是采用正交多项式进行曲线拟合的好处之一,且C、D分别求解,系数矩阵阶次较原有的低一半左右,因而运算速度明显加快。
三、极点和留数
在已知毎子、分母多项式系数后,令分母多项式为零,此时的s值,即为极点,由极点值可进一步求得固有频率和阻尼比。
P点激励,1点响应的传递函数Hip⑸可表示为:
•式中Alpr为Hip(s)在极点sr的留数,在已知极点的情况下,把上式两边各乘以(s—sr),并令s—sr,贝U
A伽=limH/z>(5)(5-sr)--|…=
‘r.・・r
s=sr
N($)
1]($-》)($-5;〉
F-I
2(%)
(•»~£;)]](»一$丿)(»—$/)丿二IJ^r
X
N(5)
D£)|—
ULI垂l|伞槿杰的影响
上孤公』只'适用于待识别的模态数与系统的自由度数相一致的情况。
实际结构往往并非如此,除了模态稀疏的情况下,仍可应用上述计算步骤外,一般均需考虑待辨识频段以外的模态影响,即所谓剩余模态的影响。
Richardson曾用增加分了多项式阶数的办法來解决这个问题。
若在待辨识的频段内,包含两个模态峰值,分母多项式阶数n取为2N=4,而分子多项式阶数取为m=2N+3=7,当然也可增加分母多项式的阶次,同样可以取得较好的效果。
五、关于总体曲线拟合
在单点频响函数测试时,往往是一点激励,再测量原点及其余各点响应。
大多数曲线拟合是在'“局部”意义上进行的,也就是说.在一组测量中,每个测量单独地进行曲线拟合,对每一模态进彳亍模态频率、阻尼和留数的估计,因此,每个模态要估计四人参数(计賢复留数作为两个参数),若有五个模态的参数需被辨识,那么,总共有20个未知参数需被估计。
一般而论,精确地估计阻尼和留数比精确地估计频率更难,而留数义总是和阻尼估计紧紧地联系在一起的。
若阻尼有较大误差,势必导致留数估计有较大误差。
为了减少误差,可以把曲线拟合过程分为两个步骤:
第一步,估计频率和阻尼,不同测量点的同一模态的固有频率和阻尼应是相同的。
故借助所有测量数据,先求得固有频率和阻尼;第二步,用已知的频率和阻尼值,进一步分别估计振型或留数。
这种估计过程,我们称为总体曲线拟合。
总体曲线拟合的优点是借助所有测量数据.这比从单个测量的曲线拟合能得到更精确的频率和阻尼的估计・。
再因阻尼是作为第一步拟合的结果,是已知的,所以,通常在第二步估计中能更精确地估计留数。
总体频率和阻尼的直接估计方法可如下进行:
因为同一模态,不同豔嶠翊数疇離翱觸篇蠢麴腸飜聽点的频响
【例】
有一四自由度系统,以正交多项式对原点位移导纳进行曲线拟合,拟合频段选在19-23Hz之间的密集模态频段内。
为说明该方法对具有一定噪声污染的测量数据的效性,为此分别对频响函数的实部及虚部加以5%、10%的随机噪声。
表列出固有频率与阻尼比的拟合结果。
图为具有5%的随机噪声的频响函数、经拟合后的理论曲线(如虚线所示)与实际曲线的比较。
0曙噪声
5%嗓简
10%噪祢
固右報率(HJ
2069
20.68
20.67
阻尼比
0.0!
0.01
0.01
固有频率(Hz)
2148
21.48
2L49
宠・「泯心
阻尼比
0.0!
0.009
0.009
由上所述,正交多项式曲线拟合方法即使在相邻模态藕合只差0.8Hz,噪声污染达10%的频响函数情况下进行曲线拟合,也能得到较好的效果,从所给出的数据中可以看出,固有频率与阻尼比几乎相等,而振型稍有差别。
原点位移异纳曲线拟合比较