基于柔度变化率与模式匹配的震后结构损伤快速识别方法与流程

技术特征：
1.一种基于柔度变化率与模式匹配的震后结构损伤快速识别方法，其特征在于，所述震后结构损伤快速识别方法包括下列步骤:s1、利用有限元软件，对识别对象结构建立结构的有限元模型；由损伤位置和损伤程度构成的排列组合的损伤情况形成损伤模式库中的所有损伤工况，并对所有损伤工况进行模态分析，得到各种损伤工况下识别对象结构的模态振型、固有频率；s2、计算各种损伤工况的柔度变化率，建立带损伤位置和损伤程度标签的以柔度变化率为参数的损伤模式库，柔度变化率的计算过程如下：式中，f为柔度矩阵，f
ij
的物理意义是在位置j测量点施加单位集中力时，在位置i测量点产生的位移，r为模态阶数,ω
r
为识别对象结构第r阶的固有频率，φ
i(r)
、φ
j(r)
分别表示位置i测量点的第r阶模态振型值与位置j测量点的第r阶模态振型值的乘积，i＝1,2,
…
,m，j＝1,2,
…
,n；计算无损和损伤工况的柔度矩阵变化量δf，计算公式如下：δf＝f
d
‑
f
u
式中，f
u
、f
d
分别为无损和损伤工况的柔度矩阵；通过柔度矩阵对角线元素计算柔度矩阵变化率δf，计算公式如下：其中，diag(
·
)为柔度矩阵的对角线元素，f
iiu
、f
iid
分别表示无损和损伤工况柔度矩阵的对角元素，[]
′
表示行向量的转置；损伤模式库中各种损伤的匹配向量如下：δf
p
＝[δf
1 δf
2 δf3…
δf
p
…
δf
q
]其中，δf
p
表示损伤模式库中的第p种损伤工况的柔度变化率，p＝1,2,
…
,q；s3、在待进行识别对象结构上安装加速度传感器，测量识别对象结构的加速度响应；s4、对所测识别对象结构的动力响应进行模态识别，得到识别对象结构的第r阶实测固有频率ω
tr
和振型φ
t(r)
＝[φ
t1(r) φ
t2(r)
…
φ
tj(r)
…
φ
tn(r)
]'，其中，φ
tj(r)
为位置j测量第r阶实测振型值，j＝1,2,
…
,n；s5、将实测振型和固有频率按照步骤s2计算得到识别对象结构实测柔度矩阵t：
其中，ω
tr
为识别对象结构第r阶的固有频率，φ
ti(r)
、φ
tj(r)
表示识别对象结构位置i测量点的第r阶模态振型值与位置j测量点的第r阶模态振型值的乘积，计算柔度变化率向量形成测试匹配向量集δt
k
：δt
k
＝[δt
1 δt
2 δt3…
δt
k
…
δt
p
]其中，δt
k
表示测试集中的第k种待测工况，k＝1,2,
…
,p；s6、将测试匹配向量集中的匹配向量δt
k
依次与步骤s2中建立的损伤模式库中各种损伤的匹配向量δf
p
进行相似度度量；s7、确定相似度度量的匹配结果最小值对应的损伤模式库中编号下的损伤工况即为该实测工况的损伤。2.根据权利要求1所述的基于柔度变化率与模式匹配的震后结构损伤快速识别方法，其特征在于，所述步骤s6中采用欧氏距离法进行相似度度量，计算公式给如下：其中，δf
p
(a)为损伤模式库中第a个损伤工况的柔度变化率向量；δt
k
(b)为测试集中第b个待测工况的柔度变化率向量。3.根据权利要求1所述的基于柔度变化率与模式匹配的震后结构损伤快速识别方法，其特征在于，所述步骤s7中，当出现两个较为接近的最小值时，则认为损伤模式库中对应的这两种损伤均为实测工况的损伤。4.根据权利要求1所述的基于柔度变化率与模式匹配的震后结构损伤快速识别方法，其特征在于，所述步骤s4中，利用复模态指示函数法cmif对所测识别对象结构的动力响应进行模态识别。

技术总结
本发明公开了一种基于柔度变化率与模式匹配的震后结构损伤快速识别方法，步骤如下：建立结构的有限元模型，并进行各种损伤工况的模态分析；计算损伤的柔度变化率形成损伤模式库；对结构动力响应进行模态识别，并计算其匹配向量形成测试集；将测试集中的匹配向量依次与模式库进行相似度度量，根据匹配结果判断损伤位置和损伤程度。本发明的方法计算量小、计算速度快；可剔除激励的影响，使得建立的损伤模式库不受外界环境因素的影响；能高效快速地识别出损伤位置和损伤程度，弥补了传统基于柔度矩阵识别结构损伤方法不能量化损伤程度的不足。因此，本发明的方法适用于结构震后损伤与安全评估，对震后的紧急救援和应急决策具有重大的意义。重大的意义。重大的意义。


技术研发人员：聂振华 李钻锋 马宏伟
受保护的技术使用者：暨南大学
技术研发日：2021.07.12
技术公布日：2021/11/9