CN105404740B - 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法 - Google Patents

一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法 Download PDF

Info

Publication number
CN105404740B
CN105404740B CN201510810161.1A CN201510810161A CN105404740B CN 105404740 B CN105404740 B CN 105404740B CN 201510810161 A CN201510810161 A CN 201510810161A CN 105404740 B CN105404740 B CN 105404740B
Authority
CN
China
Prior art keywords
matrix
sensitivity
static
joint index
damage
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.)
Expired - Fee Related
Application number
CN201510810161.1A
Other languages
English (en)
Other versions
CN105404740A (zh
Inventor
杨辰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Space Technology CAST
Original Assignee
China Academy of Space Technology CAST
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Academy of Space Technology CAST filed Critical China Academy of Space Technology CAST
Priority to CN201510810161.1A priority Critical patent/CN105404740B/zh
Publication of CN105404740A publication Critical patent/CN105404740A/zh
Application granted granted Critical
Publication of CN105404740B publication Critical patent/CN105404740B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,首先通过设置单元刚度系数与建模误差系数表征损伤与误差,其次基于静动力结果计算位移‑频率联合指标并修正,然后提取有限信息构建位移‑频率灵敏度矩阵并修正,接着应用修正的灵敏度矩阵与联合指标计算单元刚度改变量,最后改变模型误差系数计算不同误差的结果。本发明基于静动力联合指标,并考虑有限信息,通过对灵敏度方程进行修正,保证了灵敏度矩阵中的每一项都能反映静力位移采样自由度与固有频率采样阶次联合信息。本发明将会为大型复杂结构的静动信息联合损伤识别工作给出参考,并将针对众多自由度仅具备有限采样信息的损伤识别工作探索一种新的解决途径。

Description

一种基于联合指标与有限信息的修正灵敏度结构损伤识别 方法
技术领域
本发明涉及一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,尤其适用于具有众多自由度的大型复杂结构的损伤识别与健康监测工作。
背景技术
随着空间科学的进步、航天技术的发展及未来人类对空间领域的需求,现代航天器结构正在向着大型化,复杂化方向发展。包括载人飞船、深空探测卫星、空间太阳能电站与太阳帆等航天器结构,在复杂的空间服役环境中受到设计载荷作用以及各种突发性外在因素影响而面临结构的损伤积累问题,从而使结构的安全受到威胁。利用非破坏性方法来检测结构是否存在损伤,并对损伤进行定位与评估的损伤识别问题已受到国内外学术界、工程界广泛关注以及深入研究。在这其中,基于结构动力特性测试的损伤识别方法由于可以评价结构整体性、效益-造价比高、易于在线实现等特点,已经成为结构整体评价技术中最为重要的一种方法。
基于频率的损伤识别方法可以在一定程度上实现对结构的健康监测,特别是因为频率的测量方便准确,误差较低,因此可以用较为准确的频率输入参数进行损伤识别,但缺点是因为频率是全局量,其对局部损伤不敏感,或者不同的损伤会导致相同的频率变化,容易出现损伤识别“误诊”的情况。而针对静力位移的识别可以较大程度上的实现局部的损伤诊断工作,因为静力位移可以实现对结构的局部信息的精细刻画,但由于静力位移自由度的采样点众多,实现全部自由度信息的采样极其不现实。近些年,利用静力位移与固有频率的损伤识别方法逐渐引起了工程界的研究兴趣。如何利用有限静力位移响应联合少量固有频率数据就能实现较好的损伤识别结果始终是本发明所关心的重点。
因此,本发明结合静力位移与固有频率信息的联合指标,针对大型复杂结构采样的静力位移信息有限的约束,提出一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法。
发明内容
本发明的目的在于:克服现有技术的不足,提供一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,本发明能够对灵敏度方程进行修正,保证了灵敏度矩阵中的每一项都能反映静力位移采样自由度与固有频率采样阶次联合信息。本发明将会为大型复杂结构的静动信息联合损伤识别工作给出参考,并将针对众多自由度仅具备有限采样信息的损伤识别工作探索一种新的解决途径。
本发明采用的技术方案为:一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,步骤如下:
步骤(1)、计算单元刚度矩阵Ki,并通过设置单元刚度系数αi计算总体刚度矩阵K,
步骤(2)、设置实际模型Kc与参考模型Kr的建模误差系数β,
Kc=βKr (2)
步骤(3)、静力位移灵敏度识别方法中,具有n个自由度的静力学方程,
Ku=P (3)
通过设置的已知的静力载荷列向量P,计算是节点静力位移响应u,基于此构建静力位移的损伤识别灵敏度方程,
SDISP·△α=△u (4)
其中SDISP是静力位移灵敏度矩阵,具体表达式为,
[SDISP]i,k=(-K-1[K1u K2u…Kmu])i,k (5)
步骤(4)、在固有频率灵敏度识别方法中,具有n个自由度的无阻尼自由振动方程,
其中,与x(t)分别是加速度与位移向量,其特征方程为,
i=λii,i=1,2,...,n (7)
通过计算结构质量矩阵M,可以计算特征值λi(圆频率平方)与特征向量φi(模态振型),基于此构建固有频率的损伤识别灵敏度方程,
SFREQ·△α=△λ (8)
其中SFREQ是固有频率灵敏度矩阵,具体表达式为,
步骤(5)、通过步骤(3)的式(3)与步骤(4)的式(7)可以计算静力位移-固有频率联合指标矩阵△h,该矩阵的每一项△hi,j可以表达为,
△hi,j=△(uij),i=1,…,n j=1,…,n (10)
步骤(6)、提取步骤(5)联合指标矩阵△h中相应采样的静力位移自由度行与相应采样的固有频率列,
△hi,j=△(uij),i=1,…,p j=1,…,q (11)
步骤(7)、构建静力位移-固有频率灵敏度识别方法的特征方程,
S·△α=△h (12)
其中,S为静力位移-固有频率联合指标灵敏度矩阵;
步骤(8)、针对联合指标与有限信息的矩阵阶数,利用矩阵直积对步骤(7)灵敏度矩阵进行修正,
其中,所构造的算子为,
步骤(9)、考虑到联合指标的矩阵阶数,利用矩阵拉直对步骤(7)的联合指标矩阵进行修正,
步骤(10)、将步骤(8)的式(13)和(14),以及步骤(9)的式(15)代入步骤(7)的式(12),可以得到,
此时修正的灵敏度矩阵为pq×m维;
步骤(11)、由于Sij,k矩阵中的每一项可以表达为,
发现上式中存在两种独立的灵敏度矩阵中的相应元素,即,
即步骤(3)的式(5)与步骤(4)的式(9)计算结果代入(17),可计算Sij,k
步骤(12)、通过步骤(11)求得的静力位移-固有频率灵敏度矩阵S以及步骤(11)测得的相应矩阵△h,我们可以方便的计算结构刚度改变量△α,
式中为Moore-p广义逆,此时,每个单元的损伤度可计算得出。
步骤(13)、改变步骤(2)中的建模误差系数β,并重复步骤(3)-(12),计算得到在不同的建模误差下的结构刚度改变量△α。
其中,所述步骤(2)中建模误差系数β,是通过设置实际模型Kc与参考模型Kr的比值来设定的。
其中,所述步骤(8)中的灵敏度矩阵修正,是基于矩阵直积计算以及定义的相应联合指标与单元相对刚度系数灵敏度偏导数算子。
其中,所述步骤(9)中的联合指标修正,是基于矩阵拉直计算。
本发明与现有技术相比有益效果为:
(1)针对大型复杂结构信息贫乏,本发明采用静动信息联合的结构损伤识别方法。通过静力信息(位移、应变等)与动力信息(频率、振型等)对结构进行损伤识别的工作由来已久,静力识别方法方便但不精确,动力识别方法准确但较繁琐。如何能够综合二者信息进行损伤识别工作是最近国内外研究学者较为关注的。本发明正是基于此,将静动力信息进行联合,定义了一种新的指标矩阵,能够综合反映结构静动力学信息,对损伤识别的灵敏程度较静力、动力两种单工况识别方法有明显提高。
(2)在本发明定义的静动联合指标矩阵中,该矩阵的每一个元素能够充分反映每一个静力位移自由度信息与每一阶固有频率信息的综合特性,实现了对每一个静动力学信息的充分联合。
(3)针对大型复杂结构自由度众多的特点,采用有限的静力位移自由度以及有限的固有频率阶次就可以进行较为准确的损伤识别工作。对所有自由度全部进行采样不现实,本发明仅可应用有限自由度并结合有限频率便可进行识别工作,避免了大量采样带来的困难以及误差
(4)考虑联合指标矩阵的维度,本专利考虑了灵敏度矩阵维数与联合指标维数的修正,较好的实现了灵敏度识别方程的建立。对于联合指标,定义为矩阵,此时不满足灵敏度方程的对应计算维数,本发明对此进行了充分考虑,对联合指标以及灵敏度矩阵进行了修正,有效地完成了计算过程的实现性。
(5)通过设置建模误差系数来表征实际建立模型与参考模型的偏差,使从事结构损伤识别与健康监测工作的工程技术人员更能清楚的分析模型不确定性对识别工作的影响。本发明相较于在测量噪声或模态振型添加噪声的模拟,更能从建模初始层面体现误差对结果的影响。
附图说明
图1为本发明流程图;
图2为3跨15杆桁架示意图;
图3为β=1.10时识别损伤单元损伤度柱状图;
图4为不同β时识别损伤单元与未损伤单元损伤度柱状图。
具体实施方式
下面结合附图及实施例对本发明做详细说明。
本发明一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,如图1所示,步骤如下:
步骤(1)、如图2所示的3跨-15杆平面桁架结构,其中所有横纵向桁架单元长度l=0.5m,杆的横截面积A=4×10-4m2;材料密度ρ=7670kg/m3,弹性模量E=2×1011Pa,单元数量m=15;设置第6号杆单元具有15%的损伤度。基于上述参数计算单元刚度矩阵Ki,并通过设置单元刚度系数αi计算总体刚度矩阵K,
步骤(2)、设置实际模型Kc与参考模型Kr的建模误差系数β=1.10,
Kc=βKr (2)
步骤(3)、在静力位移灵敏度识别方法中,具有n=12个自由度的静力学方程,
Ku=P (3)
在第6节点处施加水平与竖直载荷各1N,建立静力载荷列向量P,计算是节点静力位移响应u,基于此构建静力位移的损伤识别灵敏度方程,
SDISP·△α=△u (4)
其中SDISP是静力位移灵敏度矩阵,具体表达式为,
[SDISP]i,k=(-K-1[K1u K2u…Kmu])i,k (5)
步骤(4)、在固有频率灵敏度识别方法中,具有n=12个自由度的无阻尼自由振动方程,
其中,与x(t)分别是加速度与位移向量,其特征方程为,
i=λii,i=1,2,…,12 (7)
通过计算结构质量矩阵M,可以计算特征值λi(圆频率平方)与特征向量φi(模态振型),基于此构建固有频率的损伤识别灵敏度方程,
SFREQ·△α=△λ (8)
其中SFREQ是固有频率灵敏度矩阵,具体表达式为,
步骤(5)、通过步骤(3)的式(3)与步骤(4)的式(7)可以计算静力位移-固有频率联合指标矩阵△h,该矩阵的每一项△hi,j可以表达为,
△hi,j=△(uij),i=1,2,…,12 j=1,2,…,12 (10)
步骤(6)、提取步骤(5)联合指标矩阵△h中相应采样的静力位移自由度行(选择3y、5y、6x共3个自由度的静力响应)与相应采样的固有频率列(前3阶固有频率),
△hi,j=△(uij),i=1,…,3 j=1,…,3 (11)
步骤(7)、构建静力位移-固有频率灵敏度识别方法的特征方程,
S·△α=△h (12)其中,S为静力位移-固有频率联合指标灵敏度矩阵。
步骤(8)、针对联合指标与有限信息的矩阵阶数,利用矩阵直积对步骤(7)灵敏度矩阵进行修正,
其中,所构造的算子为,
步骤(9)、针对联合指标的矩阵阶数,利用矩阵拉直对步骤(7)的联合指标矩阵进行修正,
步骤(10)、将步骤(8)的式(13)和(14),以及步骤(9)的式(15)代入步骤(7)的式(12),可以得到,
此时修正的灵敏度矩阵为pq×m=9×15维。
步骤(11)、由于Sij,k矩阵中的每一项可以表达为,
发现上式中存在两种独立的灵敏度矩阵中的相应元素,即,
即步骤(3)的式(5)与步骤(4)的式(9)计算结果代入(17),可计算Sij,k
步骤(12)、通过步骤(11)求得的静力位移-固有频率灵敏度矩阵S以及步骤(9)测得的相应矩阵△h,计算结构刚度改变量△α,
式中为Moore-p广义逆,此时,每个单元的损伤度可计算得出,如图3所示。
步骤(13)、改变步骤(2)中的建模误差系数β∈[1.00,1.25],并重复步骤(3)-(12),计算得到在不同的建模误差下的结构刚度改变量△α,如图4所示。
本发明未详细说明部分属于本领域技术人员公知常识。

Claims (4)

1.一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,其特征在于,该方法步骤如下:
步骤(1)、计算单元刚度矩阵Ki,并通过设置单元刚度系数αi计算总体刚度矩阵K,
步骤(2)、设置实际模型Kc与参考模型Kr的建模误差系数β,
Kc=βKr (2)
步骤(3)、在静力位移灵敏度识别方法中,具有n个自由度的静力学方程,
Ku=P (3)
通过设置的已知的静力载荷列向量P,计算节点静力位移响应u,基于此构建静力位移的损伤识别灵敏度方程,
SDISP·Δα=Δu (4)
Δα为结构刚度改变量;
其中SDISP是静力位移灵敏度矩阵,具体表达式为,
[SDISP]i,k=(-K-1[K1u K2u … Kmu])i,k (5)
步骤(4)、在固有频率灵敏度识别方法中,具有n个自由度的无阻尼自由振动方程,
其中,与x(t)分别是加速度与位移向量,其特征方程为,
i=λii,i=1,2,...,n (7)
通过计算结构质量矩阵M,可以计算特征值λi与特征向量φi,其中,特征值λi代表圆频率平方,特征向量φi代表模态振型,基于此构建固有频率的损伤识别灵敏度方程,
SFREQ·Δα=Δλ (8)
其中SFREQ是固有频率灵敏度矩阵,具体表达式为,
[SFREQ]j,k=φj TKkφj (9)
步骤(5)、通过步骤(3)的式(3)与步骤(4)的式(7)可以计算静力位移-固有频率联合指标矩阵Δh,该矩阵的每一项Δhi,j可以表达为,
Δhi,j=Δ(uij),i=1,…,n;j=1,…,n (10)
步骤(6)、提取步骤(5)联合指标矩阵Δh中相应采样的静力位移自由度行与相应采样的固有频率列,
Δhi,j=Δ(uij),i=1,…,p;j=1,…,q (11)
步骤(7)、构建静力位移-固有频率灵敏度识别方法的特征方程,
S·Δα=Δh (12)
其中,S为静力位移-固有频率联合指标灵敏度矩阵;
步骤(8)、针对联合指标与有限信息的矩阵阶数,利用矩阵直积对步骤(7)灵敏度矩阵进行修正,
其中,所构造的算子为,
步骤(9)、针对联合指标的矩阵阶数,利用矩阵拉直对步骤(7)的联合指标矩阵进行修正,
步骤(10)、将步骤(8)的式(13)和(14),以及步骤(9)的式(15)代入步骤(7)的式(12),可以得到,
此时修正的灵敏度矩阵为pq×m维;
步骤(11)、由于Sij,k矩阵中的每一项可以表达为,
发现上式中存在两种独立的灵敏度矩阵中的相应元素,即,
即步骤(3)的式(5)与步骤(4)的式(9)计算结果代入(17),可计算Sij,k
步骤(12)、通过步骤(11)求得的静力位移-固有频率灵敏度矩阵S以及步骤(9)测得的相应矩阵Δh,计算结构刚度改变量Δα,
式中为Moore-p广义逆,此时,每个单元的损伤度可计算得出;
步骤(13)、改变步骤(2)中的建模误差系数β,并重复步骤(3)-(12),计算得到在不同的建模误差下的结构刚度改变量Δα。
2.根据权利要求1所述的一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,其特征在于:所述步骤(2)中建模误差系数β,是通过设置实际模型Kc与参考模型Kr的比值来设定的。
3.根据权利要求1所述的一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,其特征在于:所述步骤(8)中的灵敏度矩阵修正,是基于矩阵直积计算以及定义的相应联合指标与单元相对刚度系数灵敏度偏导数算子。
4.根据权利要求1所述的一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法,其特征在于:所述步骤(9)中的联合指标修正,是基于矩阵拉直计算。
CN201510810161.1A 2015-11-20 2015-11-20 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法 Expired - Fee Related CN105404740B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510810161.1A CN105404740B (zh) 2015-11-20 2015-11-20 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510810161.1A CN105404740B (zh) 2015-11-20 2015-11-20 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法

Publications (2)

Publication Number Publication Date
CN105404740A CN105404740A (zh) 2016-03-16
CN105404740B true CN105404740B (zh) 2018-07-03

Family

ID=55470227

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510810161.1A Expired - Fee Related CN105404740B (zh) 2015-11-20 2015-11-20 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法

Country Status (1)

Country Link
CN (1) CN105404740B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106257529B (zh) * 2016-08-03 2019-06-18 中国空间技术研究院 基于区间有效独立法及其可能度计算的传感器配置方法
CN106289947B (zh) * 2016-08-31 2021-12-24 河北新途科技有限公司 一种轻质高强梁结构损伤识别方法
CN107133389A (zh) * 2017-04-14 2017-09-05 中国空间技术研究院 面向结构健康监测基于区间鲁棒优化的传感器配置方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104484502A (zh) * 2014-11-21 2015-04-01 华中科技大学 一种基于正向子结构的有限元模型修正方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104484502A (zh) * 2014-11-21 2015-04-01 华中科技大学 一种基于正向子结构的有限元模型修正方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《Damage identification by the eigenparameter decomposition of structural flexibility change》;Q. W. Yang,J. K. Liu;《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING》;20081110;第444-459页 *
《On an iterative general-order perturbation method for multiple structural damage detection》;C.N. Wong, W.D. Zhu, G.Y. Xu;《Journal of Sound and Vibration》;20041231;第363-386页 *
《基于动静测试的融合灵敏度结构损伤识别方法》;杨辰,王晓军,邱志平;《应用力学学报》;20131030;第30卷(第5期);第669-673页 *
《结构损伤识别的联合灵敏度方法研究》;邵晓蓉,杨秋伟,赵卫;《机械强度》;20121231;第317-321页 *

Also Published As

Publication number Publication date
CN105404740A (zh) 2016-03-16

Similar Documents

Publication Publication Date Title
WO2020192623A2 (zh) 一种旋转叶片非接触式动应变场测量方法及其系统
CN108195535B (zh) 基于非线性激振特征的螺栓结合部松动检测方法及系统
Esfandiari et al. Finite element model updating using frequency response function of incomplete strain data
Pak Wing shape sensing from measured strain
Guan et al. Improved damage detection method based on element modal strain damage index using sparse measurement
CN104344993B (zh) 一种构件承载力和材料性能参数的试验测量方法
CN106289947B (zh) 一种轻质高强梁结构损伤识别方法
CN105404740B (zh) 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法
CN107862170B (zh) 一种基于动态缩聚的有限元模型修正方法
CN110108430B (zh) 一种纤维增强复合材料梁的高精度损伤定位方法
CN103148971B (zh) 一种测试超高压管式反应器热套端部结构局部应力场的方法
CN102323058B (zh) 基于转轴应变信号的汽轮发电机组轴承载荷识别装置和方法
CN110068406A (zh) 基于静载识别的四边简支薄板结构光纤应变场重构方法
CN104502010A (zh) 未知边界条件下索力测试方法及装置
CN102252792B (zh) 一种杆件绝对轴力测试方法
CN112326787A (zh) 一种基于专属测试车多点快速静态采集的梁式桥识别方法
CN110133103A (zh) 基于平面脉冲声波激振的纤维增强复合材料参数辨识方法
CN118095016B (zh) 基于逆元法与有限元刚度矩阵构建的载荷识别方法
Hu et al. Multi-type sensor placement and response reconstruction for building structures: Experimental investigations
CN106383003A (zh) 基于柔度识别的索结构索力的测量方法及测量系统
Wu et al. Crack diagnosis method for a cantilevered beam structure based on modal parameters
CN105548259B (zh) 一种卫星结构热稳定性试验方法
CN107101800B (zh) 输电铁塔模型的风力测量设备及其测量方法
CN113639941A (zh) 一种测试桥梁应变柔度矩阵的方法
Tsushima et al. Monitoring multi-axial vibrations of flexible rockets using sensor-instrumented reference strain structures

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180703

Termination date: 20181120