CN108363679B - 一种自动追踪结构模态参数的方法 - Google Patents
一种自动追踪结构模态参数的方法 Download PDFInfo
- Publication number
- CN108363679B CN108363679B CN201810211974.2A CN201810211974A CN108363679B CN 108363679 B CN108363679 B CN 108363679B CN 201810211974 A CN201810211974 A CN 201810211974A CN 108363679 B CN108363679 B CN 108363679B
- Authority
- CN
- China
- Prior art keywords
- ref
- mode
- modal
- modality
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000013598 vector Substances 0.000 claims abstract description 37
- 230000005284 excitation Effects 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000005314 correlation function Methods 0.000 claims description 8
- 238000013016 damping Methods 0.000 claims description 7
- 230000000295 complement effect Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005624 perturbation theories Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明属于结构健康监测技术领域,提供了一种自动追踪结构模态参数的方法。首先,利用自然激励技术处理随机响应,并采用特征系统实现算法结合稳定图识别在不同时段的模态参数;然后,根据指定被追踪时段内各模态观测向量与前一时段追踪后参考模态所构成的子空间之间的相关性,将被追踪时段内的模态划分为可追踪模态和不可追踪模态两类;最后,将可追踪模态按照模态观测向量相关性最大且频率偏差最小的原则,逐一归入具有相同结构特性的模态类;同时,将不可追踪模态与原参考模态的合集作为新的参考模态,用于后一时段的模态追踪。本发明算法无需人为设定阈值和参考模态,可实现模态参数的自动追踪。
Description
技术领域
本发明属于结构健康监测技术领域,涉及一种自动追踪结构模态参数的方法。
背景技术
结构的服役性能可通过结构模态参数随时间的变化来反映。因此,自动准确地识别出结构的模态参数十分必要。目前广泛采用的模态参数识别方法包括最小二乘复频域法、随机子空间法和特征系统实现算法等。为获取模态参数的实时变化,这些识别方法均将结构响应按时间分为许多子段,然后对每一子段的响应时程进行识别,进而获取各时间段下的模态参数。然而受激励水平、环境干扰以及算法稳定性的影响,各时段获取的模态数量未必相同且各阶模态并非一一对应。模态追踪技术的目的就是保证在不同时段内识别的结构模态能够保持正确的前后对应关系,不发生“模态交叉”现象。
现有的模态追踪方法主要分为三类:1)人工排序法:根据经验归类前后两个时段内识别出的模态参数,这种方法人为参与的工作量较大;2)容许限值法:依照经验设定频率偏差或模态置信度(Modal Assurance Criterion,简称MAC)容许限值来追踪,一般包括固定容限值和自适应调整容限值两种。这类方法若阈值设定不合理,会存在模态错误归类或丢失的问题;3)预测-校正法:基于摄动理论预测后一时段的模态参数,然后比较预测的模态参数与识别的模态参数。这类方法预测计算的效率较低,难以在实际大型工程中应用。因此,精确且无需人为分析的自动模态追踪技术具有重要的工程意义。
发明内容
本发明的目的是提供一种自动追踪结构模态参数的方法,解决在实际工程中实时提取模态参数时,由于经验阈值不合理造成模态追踪不准确或参考模态选择不合理造成模态缺失的问题。
本发明的技术方案是:提出一种自动追踪结构模态参数的方法,其特点是利用自然激励技术处理结构随机激励响应,并采用特征系统实现算法结合稳定图提取结构在不同时段的模态参数;以第一个时间段内识别的各阶模态参数为初始参考模态参数,对参考模态观测向量矩阵进行奇异值分解,获得参考模态子空间及其正交补空间,利用被追踪时段各模态观测向量与参考模态子空间的相关性,将其划分为可追踪模态和不可追踪模态类。对于可追踪模态,利用模态观测向量相关性最大以及频率偏差最小的原则进行追踪,同时,用不可追踪模态与原参考模态的合集来自适应更新参考模态,进行下一时段的模态追踪。
本发明的技术方案:
一种自动追踪结构模态参数的方法,步骤如下:
步骤一:获取不同时段下的模态参数
(1)针对待追踪结构,采集第h时段在第z测点的随机响应y(t)=[y1(t),y2(t),…,yz(t)]T,t=1,2,…,N,N为样本时程点数,利用自然激励方法获得各个时间延迟下的相关函数矩阵r(τ)如下:
式中:rij(τ)为测点i和测点j处加速度响应间的互相关函数;
根据得到的相关函数矩阵r(τ)构造Hankel矩阵Hms(k-1)和Hms(k):
令k=1,对矩阵Hms(k-1)利用特征系统实现算法求解各个模型阶次下的模态参数,包括系统特征值、频率、阻尼比、振型向量和模态观测向量;模型阶次从偶数δ开始,按照偶数倍δ递增到nuδ;
设定频率偏差限值ef,lim、阻尼比偏差限值eξ,lim和模态置信度指标限值eMAC,lim,将满足这三个限值的模态作为稳定模态;若相邻模型阶次下的两个稳定模态同时满足频率偏差小于ef,lim和模态置信度大于eMAC,lim,则将其归为同一类;对于模态数量大于限值ntol的类称为物理类,在各物理类内按照模态参数平均值选择物理模态代表,从而获得第h时段下各模态对应的模态参数,包括系统特征值λi,h、频率fi,h、阻尼比ξi,h、振型向量和模态观测向量上标T表示转置;
步骤二:追踪各时段的模态参数
当h=1,以此时段内获取的β阶模态作为初始参考模态,并将此时段的频率f1=[f1,1,f2,1,…,fβ,1]和模态观测向量矩阵W1=[w1,1,w2,1,…,wβ,1],分别标记为参考频率向量fref=[f1,ref,f2,ref,…,fβ,ref]和参考模态观测向量矩阵Wref=[w1,ref,w2,ref,…,wβ,ref];
对参考模态观测向量矩阵Wref进行奇异值分解,获得参考模态子空间U1及其正交补空间U2:
式中:H表示复共轭转置;由于同一时段下识别的各模态间线性无关,故矩阵Wref的秩与参考模态的阶数保持一致;
当h=2时段时,获取α阶模态,其中频率为f2=[f1,2,f2,2,…,fα,2],模态观测向量矩阵为W2=[w1,2,w2,2,…,wα,2];对于h=2时段下的模态j,分别计算其模态观测向量wj,2与参考模态子空间U1,以及与正交补空间U2的相关性wMOC:
若wMOC(U1,wj,2)≥wMOC(U2,wj,2),则将模态j定义为可追踪模态,否则定义为不可追踪模态;假定在识别的α阶模态中,共有η阶可追踪模态,η≤α,此时产生α-η阶不可追踪模态;将可追踪模态的频率重新表示为模态观测向量矩阵表示为若h=2时段下的第阶可追踪模态与第χ阶参考模态满足:
则追踪为同一类模态,其中MOC表示两个模态观测向量间的相关性;
此外,将第h=2时段下的α-η阶不可追踪模态与原参考模态一起,构造新的参考模态,用于下一时刻的模态追踪;此时,新的参考频率向量扩展为fref=[f1,ref,f2,ref,…,fβ,ref,fβ+1,ref,…,fβ+α-η,ref],参考模态观测向量矩阵扩展为Wref=[w1,ref,w2,ref,…,wβ,ref,wβ+1,ref,…,wβ+α-η,ref];对于h=3,4,…时刻,模态追踪步骤同h=2时的模态追踪步骤。
本发明的有益效果:无需选定初始参考模态,且在追踪过程中无需人工参与和设置阈值,可实现自动准确的追踪模态参数。
附图说明
图1是某大桥主梁竖向加速度传感器布置图。
图2是自动追踪主梁结构模态参数的结果。
具体实施方式
以下结合技术方案和附图,进一步阐明本发明的实施方式。
设在桥梁主梁上布置14个竖向加速度传感器,如图1所示。以100Hz的采样频率采集环境激励下的竖向加速度响应数据,采集时间为2016年9月1日到2016年9月7日,每次以一个小时的响应时程进行模态参数识别。
具体实施方案如下:
(1)将2016年9月1日0:00-1:00的结构随机响应所在时段定为h=1,首先对响应y(t)=[y1(t),y2(t),…,y14(t)]T,t=1,2,…,N,利用自然激励技术获得各个时间延迟下的相关函数矩阵,如公式(1)。
(2)令m=250,s=250;分别选取τ=1~499和τ=2~500处的相关函数矩阵构造Hankel矩阵Hms(0)和Hms(1),如公式(2)。
(3)令最小计算阶次为δ=4,并按照δ=4的倍数递增至280,共选定阶次nu=70;对Hankel矩阵Hms(0)和Hms(1),利用特征系统实现算法计算在各个计算阶次下的模态参数。
(4)设定频率误差限值ef,lim=5%、阻尼比误差限值eξ,lim=20%、模态置信度指标限值eMAC,lim=90%,将满足这三个限值的模态作为稳定模态。若相邻模型阶次下的两个稳定模态同时满足频率偏差小于ef,lim和模态置信度大于eMAC,lim,则归为同一类。对于模态数量大于限值ntol=0.5nu的类称为物理类,在各物理类内按照模态参数平均值选择物理模态代表。从而获得此时段下频率小于2.5Hz的β=15个模态,作为初始参考模态,其中参考频率:f1,ref=0.387Hz,f2,ref=0.648Hz,f3,ref=0.754Hz,f4,ref=0.932Hz,f5,ref=0.985Hz,f6,ref=1.060Hz,f7,ref=1.278Hz,f8,ref=1.321Hz,f9,ref=1.513Hz,f10,ref=1.605Hz,f11,ref=1.685Hz,f12,ref=1.954Hz,f13,ref=2.000Hz,f14,ref=2.038Hz,f15,ref=2.212Hz。
(5)对于h=2时段内识别出α=16阶模态,其中频率为f1,2=0.386Hz,f2,2=0.644Hz,f3,2=0.755Hz,f4,2=0.929Hz,f5,2=0.983Hz,f6,2=1.061Hz,f7,2=1.257Hz,f8,2=1.318Hz,f9,2=1.503Hz,f10,2=1.595Hz,f11,2=1.676Hz,f12,2=1.949Hz,f13,2=1.998Hz,f14,2=2.033Hz,f15,2=2.220Hz,f16,2=2.253Hz。
(6)对参考模态观测向量矩阵Wref进行奇异值分解,获得参考模态子空间U1及其正交补空间U2,并分别计算h=2时段内识别出的α=16阶模态与空间U1和U2的相关性。对于模态j=1,…,15,满足wMOC(U1,wj,2)≥wMOC(U2,wj,2),进而可通过公式(6)~(9)进行模态追踪;对于模态j=16,wMOC(U1,w16,2)=0.238,wMOC(U2,w16,2)=0.762,因wMOC(U1,wj,2)<wMOC(U2,wj,2),故为不可追踪模态。将模态j=16添加到原参考模态中,用于下一时段的追踪。此时参考模态阶数为16,其中,f16,ref=2.253Hz。对于h=3,4,…时段,重复步骤(1)~(4)进行模态识别,重复步骤(6)进行模态追踪,识别结果如图2所示。
Claims (1)
1.一种自动追踪结构模态参数的方法,其特征在于,步骤如下:
步骤一:获取不同时段下的模态参数
(1)针对待追踪结构,采集第h时段在第z测点的随机响应y(t)=[y1(t),y2(t),…,yz(t)]T,t=1,2,…,N,N为样本时程点数,利用自然激励方法获得各个时间延迟下的相关函数矩阵r(τ)如下:
式中:rpq(τ)为测点p和测点q处加速度响应间的互相关函数;
根据得到的相关函数矩阵r(τ)构造Hankel矩阵Hms(k-1)和Hms(k):
令k=1,对矩阵Hms(k-1)利用特征系统实现算法求解各个模型阶次下的模态参数,包括系统特征值、频率、阻尼比、振型向量和模态观测向量;模型阶次从偶数δ开始,按照偶数倍δ递增到nuδ,其中nu为选定的阶次数;
设定频率偏差限值ef,lim、阻尼比偏差限值eξ,lim和模态置信度指标限值eMAC,lim,将满足这三个限值的模态作为稳定模态;若相邻模型阶次下的两个稳定模态同时满足频率偏差小于ef,lim和模态置信度大于eMAC,lim,则将其归为同一类;对于模态数量大于限值ntol的类称为物理类,在各物理类内按照模态参数平均值选择物理模态代表,从而获得第h时段下各模态对应的模态参数,包括系统特征值λi,h、频率fi,h、阻尼比ξi,h、振型向量和模态观测向量上标T表示转置;
步骤二:追踪各时段的模态参数
当h=1,以此时段内获取的β阶模态作为初始参考模态,并将此时段的频率f1=[f1,1,f2,1,…,fβ,1]和模态观测向量矩阵W1=[w1,1,w2,1,…,wβ,1],分别标记为参考频率向量fref=[f1,ref,f2,ref,…,fβ,ref]和参考模态观测向量矩阵Wref=[w1,ref,w2,ref,…,wβ,ref];
对参考模态观测向量矩阵Wref进行奇异值分解,获得参考模态子空间U1及其正交补空间U2:
式中:H表示复共轭转置;由于同一时段下识别的各模态间线性无关,故矩阵Wref的秩与参考模态的阶数保持一致;
当h=2时段时,获取α阶模态,其中频率为f2=[f1,2,f2,2,…,fα,2],模态观测向量矩阵为W2=[w1,2,w2,2,…,wα,2];对于h=2时段下的模态j,分别计算其模态观测向量wj,2与参考模态子空间U1,以及与正交补空间U2的相关性wMOC:
若wMOC(U1,wj,2)≥wMOC(U2,wj,2),则将模态j定义为可追踪模态,否则定义为不可追踪模态;假定在识别的α阶模态中,共有η阶可追踪模态,η≤α,此时产生α-η阶不可追踪模态;将可追踪模态的频率重新表示为模态观测向量矩阵表示为若h=2时段下的第l阶可追踪模态与第χ阶参考模态满足:
则追踪为同一类模态,其中MOC表示两个模态观测向量间的相关性;
此外,将第h=2时段下的α-η阶不可追踪模态与原参考模态一起,构造新的参考模态,用于下一时刻的模态追踪;此时,新的参考频率向量扩展为fref=[f1,ref,f2,ref,…,fβ,ref,fβ+1,ref,…,fβ+α-η,ref],参考模态观测向量矩阵扩展为Wref=[w1,ref,w2,ref,…,wβ,ref,wβ+1,ref,…,wβ+α-η,ref];对于h=3,4,…时刻,模态追踪步骤同h=2时的模态追踪步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810211974.2A CN108363679B (zh) | 2018-03-12 | 2018-03-12 | 一种自动追踪结构模态参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810211974.2A CN108363679B (zh) | 2018-03-12 | 2018-03-12 | 一种自动追踪结构模态参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108363679A CN108363679A (zh) | 2018-08-03 |
CN108363679B true CN108363679B (zh) | 2021-04-20 |
Family
ID=63000284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810211974.2A Active CN108363679B (zh) | 2018-03-12 | 2018-03-12 | 一种自动追踪结构模态参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108363679B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111274630B (zh) * | 2020-01-15 | 2022-09-20 | 大连理工大学 | 一种用于工程结构柔度识别的物理模态提取方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101901209A (zh) * | 2010-02-04 | 2010-12-01 | 福州大学 | 基于改进emd和arma模型的结构响应分析方法 |
CN103344448A (zh) * | 2013-06-26 | 2013-10-09 | 中国路桥工程有限责任公司 | 一种桥梁结构损伤识别方法和系统 |
WO2015040483A1 (es) * | 2013-09-18 | 2015-03-26 | Comsa Emte, S.L. | Método y sistema de evaluación de la integridad estructural de torres de celosía o reticuladas |
-
2018
- 2018-03-12 CN CN201810211974.2A patent/CN108363679B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101901209A (zh) * | 2010-02-04 | 2010-12-01 | 福州大学 | 基于改进emd和arma模型的结构响应分析方法 |
CN103344448A (zh) * | 2013-06-26 | 2013-10-09 | 中国路桥工程有限责任公司 | 一种桥梁结构损伤识别方法和系统 |
WO2015040483A1 (es) * | 2013-09-18 | 2015-03-26 | Comsa Emte, S.L. | Método y sistema de evaluación de la integridad estructural de torres de celosía o reticuladas |
Non-Patent Citations (2)
Title |
---|
基于加速度二次协方差矩阵参数变化比法的环境振动下结构损伤识别;王立新 等;《振动与冲击》;20160831;第35卷(第8期);第143-158页 * |
基于小波变换的结构模态参数识别;伊廷华 等;《振动工程学报》;20060331;第19卷(第1期);第51-56页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108363679A (zh) | 2018-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110674752B (zh) | 一种基于隐马尔可夫模型的刀具磨损状态识别与预测方法 | |
Ding et al. | Structural damage detection using artificial bee colony algorithm with hybrid search strategy | |
Ferreira et al. | On the block maxima method in extreme value theory: PWM estimators | |
CN111222549B (zh) | 一种基于深度神经网络的无人机故障预测方法 | |
CN108415884B (zh) | 一种结构模态参数实时追踪方法 | |
CN105843073A (zh) | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 | |
WO2019173943A1 (zh) | 一种自动追踪结构模态参数的方法 | |
CN107084714A (zh) | 一种基于RoboCup3D的多机器人协作目标定位方法 | |
CN105117593A (zh) | 基于小波变换和粒子群改进灰色模型的短期风速预测方法 | |
CN106529185B (zh) | 一种古建筑位移的组合预测方法及系统 | |
CN110442911B (zh) | 一种基于统计机器学习的高维复杂系统不确定性分析方法 | |
CN110377942A (zh) | 一种基于有限高斯混合模型的多模型时空建模方法 | |
CN103577676A (zh) | 污水处理工艺综合评价的灰色加权法 | |
CN117313029A (zh) | 一种基于卡尔曼滤波参数提取与状态更新的多传感器数据融合方法 | |
CN109508788A (zh) | 一种基于arma模型的sdn流量预测方法 | |
CN106228026A (zh) | 一种基于最优退化特征量的剩余寿命预测算法 | |
CN115270239A (zh) | 基于动力特性和智能算法响应面法的桥梁可靠性预测方法 | |
CN107145937A (zh) | 基于弹性scad罚函数的回声状态网络时间序列预测方法 | |
CN108363679B (zh) | 一种自动追踪结构模态参数的方法 | |
Gokpinar et al. | Generalization of inclusion probabilities in ranked set sampling | |
CN110084324A (zh) | 基于极限学习机的卡尔曼滤波参数自适应更新方法 | |
Bedoui et al. | A new recursive algorithm for simultaneous identification of discrete time delay systems | |
CN109884893B (zh) | 一种多工序变量间动态时滞估计方法 | |
Yang et al. | An improved probabilistic neural network with ga optimization | |
CN112817282A (zh) | 一种基于数字孪生系统的高炉炉顶压力控制回路参数优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |