CN103646177B - 一种考虑温度影响的土木工程结构损伤预警方法 - Google Patents
一种考虑温度影响的土木工程结构损伤预警方法 Download PDFInfo
- Publication number
- CN103646177B CN103646177B CN201310690423.6A CN201310690423A CN103646177B CN 103646177 B CN103646177 B CN 103646177B CN 201310690423 A CN201310690423 A CN 201310690423A CN 103646177 B CN103646177 B CN 103646177B
- Authority
- CN
- China
- Prior art keywords
- matrix
- sample
- model coefficient
- standard deviation
- model
- 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
- 230000006378 damage Effects 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000000513 principal component analysis Methods 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims description 47
- 239000000470 constituent Substances 0.000 claims description 11
- 230000009467 reduction Effects 0.000 claims description 7
- 230000004044 response Effects 0.000 abstract description 7
- 230000001133 acceleration Effects 0.000 abstract description 6
- 230000008859 change Effects 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000036541 health Effects 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 4
- 208000027418 Wounds and injury Diseases 0.000 description 3
- 208000014674 injury Diseases 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 238000000700 time series analysis Methods 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明涉及土木工程结构损伤识别领域,特别涉及一种考虑温度影响的土木工程结构损伤预警方法,利用时间序列分析中的AR模型和多元统计分析中的主成分分析并结合标准差控制图进行结构的损伤预警研究。首先利用AR模型对结构损伤前后的加速度响应数据进行拟合,并提取模型系数,其次利用主成分分析去除温度对AR模型系数的影响,最后利用标准差控制图进行损伤预警。本发明理论严谨,方法新颖,方案合理,可操作性强。该方法的优点是直接利用加速度响应,不需要有限元模型和模态参数,属于数据驱动的方法,适于进行实时的结构健康监测,同时该方法在温度变化的条件下,具有的抗噪声能力,成功地进行了损伤预警。
Description
技术领域
本发明涉及土木工程结构损伤识别领域,特别涉及一种考虑温度影响的土木工程结构损伤预警方法。
背景技术
结构损伤会引起动力响应特征参数的变化,而温度的改变同样也会导致动力响应特征参数发生变化,有时甚至会淹没损伤引起的动力响应特征参数的变化,从而导致基于振动的结构损伤识别方法失效。
吴森等,基于AR模型和主成分分析的损伤识别方法,振动、测试与诊断,2012.10,32(5),847-845.提出了一种基于结构加速度时间序列的损伤识别方法。首先提取结构损伤前后的AR模型系数,构成多个原始数据矩阵;其次,对多个原始数据矩阵分别进行主成分分析得到前两阶主成分,并建立相应的椭圆控制图,以前两阶主成分在控制椭圆中的分布情况来判别结构是否存在损伤,但是,该方法未考虑温度对土木工程结构损伤预警的影响。
发明内容
本发明针对现有技术中的不足,提供一种考虑温度影响的土木工程结构损伤预警方法。
为消除温度的影响,本专利利用时间序列分析中的AR模型和多元统计分析中的主成分分析并结合标准差控制图进行结构的损伤预警研究。首先利用AR模型对结构损伤前后的加速度响应数据进行拟合,并提取模型系数,其次利用主成分分析去除温度对AR模型系数的影响,最后利用标准差控制图进行损伤预警。四层框架数值模拟结果表明,在温度变化的条件下,该方法成功地进行了损伤预警。
本发明采取的技术方案是:
一种考虑温度影响的土木工程结构损伤预警方法,具体步骤如下:
(1)观测值{Xt(t=1,2,…,p)}为一平稳的、零均值的随机过程,若Xt的取值与前p步的各个取值Xt-1、Xt-2、…Xt-p有关,则可以用差分方程来表示它们的内在联系:
式中,p代表模型阶次,为自回归系数(AR模型系数),表示t-i时刻的“历史值”Xt-i对t时刻值Xt的贡献(或者影响),{ɑt}为均值为0、方差为σɑ 2白噪声序列。该模型所描述的{xi}是对其自身过去数值进行回归,故称为自回归模型,简称AR模型。
(2)通过对AR模型系数进行主成分分析,求主成分分析结果残差的方法,可减少或剔除温度对AR模型系数的影响。
利用主成分分析减少或剔除温度影响的过程简述如下:
a)假设由AR模型系数构建的样本矩阵Y∈Rn×N,对样本矩阵Y的协方差矩阵进行奇异值分解:
YYT=USUT (2)
上式中U为正交矩阵,即UUT=I,矩阵U中第i列向量定义Y的第i个主成分;矩阵S为奇异值矩阵。
b)前m阶主成分的取值通过指标I来确定:
式中,I应大于一定门槛值e(e可取值为90%),门槛值说明前m个主成分能反映原变量Y信息量的百分比,也就是贡献率达到e。
c)对样本矩阵Y进行降维(从n维降到m维),方程可以表示为:
X=TY (4)
其中,矩阵X∈Rm×N为得分矩阵,T∈Rm×n为荷载矩阵,矩阵T由矩阵U中的前m列组成。
d)将降维得到的矩阵X再投影至原向量空间,计算得到样本矩阵Y的残差E。
残差E即为剔除了温度影响的AR模型系数,残差E的计算过程即为对样本矩阵Y剔除温度影响的过程。
(3)利用结构损伤前后剔除温度影响的AR模型系数构造欧氏距离作为样本,计算样本的标准差S,将标准差S作为统计量,绘制标准差控制图进行结构损伤预警;
定义无损状态下的AR模型系数与损伤状态下的AR模型系数的欧氏距离如下:
式中,k=(1,2,…,p),Ek u为无损状态下剔除温度影响的p个AR模型系数所组成的向量,Ek d为损伤状态下剔除温度影响的p个AR模型系数所组成的向量;
建立控制图需要计算控制图的中心限和上下控制界限。设欧氏距离d为正态总体样本X,X~N(μ,σ),x为取自总体样本X的样本,标准差S作为样本统计量。由《均值及标准偏差控制图及其应用》查得标准差图控制限计算公式如下:
其中,n为样本容量;c4为与样本大小有关的常数,可以由《均值及标准偏差控制图及其应用》附表查得;S为样本的标准差。
在控制图中可由点出界的个数来进行损伤判别。利用3σ准则取显著性水平为0.9973来控制置信水平的上下限。在结构正常状况下,连续100个点中有3个或多于3个点超出控制界限的概率为0.0026,是小概率事件,此概率值接近0.0027(1-0.9973=0.0027),与3σ相近,因此采用的损伤判断准则可以近似为:连续100个点中有3个或者多于3个点超出控制界限就可以判定结构发生损伤。
本发明的有益效果是:
本发明理论严谨,方法新颖,方案合理,可操作性强。该方法的优点是直接利用加速度响应,不需要有限元模型和模态参数,属于数据驱动的方法,适于进行实时的结构健康监测,同时该方法在温度变化的条件下,具有的抗噪声能力,成功地进行了损伤预警。
附图说明
图1为四层框架的三维有限元模型示意图;
1为柱1;2为柱2;3为柱3;4为柱4;
5为X向面二;6为Y向面二;7为X向面一;8为Y向面一;
9为损失杆件;10为采集集点;
图2为钢的弹性模量和温度的关系的示意图;
图3为无噪声时剔除温度影响前后的S控制图,
其中:(a)为去除温度影响前的S控制图;
(b)为去除温度影响后的S控制图;
图4为噪声程度为3%时剔除温度影响前后的S控制图;
其中:(a)为去除温度影响前的S控制图;
(b)为去除温度影响后的S控制图;
图5为本发明方法流程图。
具体实施方式
本发明的具体实施方式如下:
实施例1:
一种温度因素影响下的土木工程结构损伤预警方法,具体步骤如下:
(1)建立AR模型:
式中,p代表模型阶次,为自回归系数(AR模型系数),表示t-i时刻的“历史值”Xt-i对t时刻值Xt的贡献或影响,{ɑt}为均值为0、方差为σɑ 2白噪声序列;
(2)通过对上述AR模型系数进行主成分分析,求主成分分析结果残差的方法,减少环境温度对AR模型系数的影响;
a)假设由AR模型系数构建的样本矩阵Y∈Rn×N,对样本矩阵Y的协方差矩阵进行奇异值分解:
YYT=USUT (2)
上式中U为正交矩阵,即UUT=I,矩阵U中第i列向量定义Y的第i个主成分;矩阵S为奇异值矩阵。
b)前m阶主成分的取值通过指标I来确定:
式中,I大于门槛值e,e取值为90%,m取值为1;
c)对样本矩阵Y进行降维,方程可以表示为:
X=TY (4)
其中,矩阵X∈Rm×N为得分矩阵,T∈Rm×n为荷载矩阵,矩阵T由矩阵U中的前m列组成;
d)将降维得到的矩阵X再投影至原向量空间,计算得到样本矩阵Y的残差E:
残差E即为剔除了环境温度影响的AR模型系数;
(3)利用结构损伤前后剔除温度影响的AR模型系数构造欧氏距离作为样本,计算样本的标准差S,将标准差S作为统计量,绘制标准差控制图进行结构损伤预警;
定义无损状态下的AR模型系数与损伤状态下的AR模型系数的欧氏距离如下:
式中,k=(1,2,…,p),Ek u为无损状态下剔除温度影响的p个AR模型系数所组成的向量,Ek d为损伤状态下剔除温度影响的p个AR模型系数所组成的向量;
建立控制图需要计算控制图的中心限和上下控制界限,设欧氏距离d为正态总体样本X,X~N(μ,σ),x为取自总体样本X的样本,标准差S作为样本统计量,由《常规控制图》查得标准差图控制限计算公式如下:
其中,n为样本容量;c4为与样本大小有关的常数,可以由《常规控制图》附表查得;S为样本的标准差,
在控制图中由点出界的个数来进行损伤判别,连续100个点中有3个或者多于3个点超出控制界限就可以判定结构发生损伤。
实施例2:
采用ANSYS建立了一个四层框架的三维有限元模型,如图1所示,基本参数:密度为7800Kg/m3,泊松比μ=0.3。柱、横梁以及支撑采用BEAM4单元,顶板采SHELL63单元,顶板上堆积质量采用MASS21单元进行模拟。
该模型共32个结点,192个自由度,50个BEAM4单元,1个SHELL3单元,4个MASS21单元,7种单元截面类型,模型与基础固接。
本实施例以Y向面Ⅱ(6)的斜撑分别发生不同程度损伤作为研究对象,损伤是通过弹性模量的折减来实现的,温度变化范围为-20℃~40℃,钢材的弹性模量和温度关系如图2所示。
考虑温度影响下的数值模拟通过设定结构的弹性模量来实现,温度变化范围为-20℃~40℃,每隔5℃进行一次环境设定并采样,分别为-20℃,-15℃,-10℃,-5℃,0℃,5℃,10℃,15℃,20℃,25℃,30℃,35℃,40℃,共计13个温度采样。
以损失杆件(9)分别损伤0、0.3、0.4、0.5、0.7为研究对象。
对每种损伤程度,在-20℃~40℃范围内每隔5℃分别采集节点(10)的前20s的10000个加速度响应数据,这样每种工况下共得到13组不同温度下的加速度响应数据,将每组10000个数据均分为5个子块,利用Matlab编制的程序获得每个子块的前4阶AR模型系数,然后利用主成分分析去除温度对各工况下AR模型系数的影响,因为利用主成分分析处理各工况每个子块的AR模型系数的方法相同,仅以损失杆件(9)损伤0.3时第一个子块的AR模型系数为例,列出利用主成分分析方法剔除温度影响前后的AR模型系数,分别见表1和表2。
按式(7)计算每种损伤工况在不同温度下的欧式距离,将各温度下欧式距离作为样本,计算欧氏距离的标准差S;
将标准差带入公式(8,9)计算控制界限,从而得到剔除温度影响前后的S控制图,如图3-图4所示。图3为无噪声时剔除温度影响前后的S控制图,图4为噪声程度为3%时剔除温度影响前后的S控制图。其中,图3(a)、图4(a)控制上限值为0.13,控制下限值为0;图3(b)、图4(b)的控制上限值为0.03,控制下限值为0。
表1损伤杆件损伤30%时未剔除温度影响前AR模型参数
表2损伤杆件损伤30%时剔除温度影响后AR模型参数
通过图3(a)和图3(b),图4(a)和图4(b)对比可得,无论有无噪声时,未去除温度影响时,结构无损时有部分点散落在控制上限之外,而结构发生损伤时有部分点散落在控制上限以内,温度的影响使S控制图不能准确反应结构的损伤状况。而利用主成分分析去除温度影响后,结构无损状态下的统计量点均分布在控制上限以内,发生损伤时的统计量点则分布在控制上限之外,能够准确地判断结构的损伤状况,从而发出损伤预警。
Claims (1)
1.一种考虑温度影响的土木工程结构损伤预警方法,具体步骤如下:
(1)建立AR模型:
式中,p代表模型阶次,为自回归系数,AR模型系数,表示t-i时刻的“历史值”Xt-i对t时刻值Xt的贡献或影响,{ɑt}为均值为0、方差为σɑ 2白噪声序列;
(2)通过对上述AR模型系数进行主成分分析,求主成分分析结果残差,减少环境温度对AR模型系数的影响;
a)假设由AR模型系数构建的样本矩阵Y∈Rn×N,对样本矩阵Y的协方差矩阵进行奇异值分解:
YYT=USUT (2)
上式中U为正交矩阵,即UUT=I,矩阵U中第i列向量定义Y的第i个主成分;矩阵S为奇异值矩阵;
b)前m阶主成分的取值通过指标I来确定:
式中,I大于门槛值e,e取值为90%;
c)对样本矩阵Y进行降维,方程可以表示为:
X=TY (4)
其中,矩阵X∈Rm×N为得分矩阵,T∈Rm×n为荷载矩阵,矩阵T由矩阵U中的前m列组成;
d)将降维得到的矩阵X再投影至原向量空间,计算得到样本矩阵Y的残差E:
残差E即为剔除了温度影响的AR模型系数;
(3)利用结构损伤前后剔除温度影响的AR模型系数构造欧氏距离作为样本,计算样本的标准差S,将标准差S作为统计量,绘制标准差控制图进行结构损伤预警;
定义无损状态下的AR模型系数与损伤状态下的AR模型系数的欧氏距离如下:
式中,k=(1,2,…,p),Ek u为无损状态下剔除温度影响的p个AR模型系数所组成的向量,Ek d为损伤状态下剔除温度影响的p个AR模型系数所组成的向量;
建立控制图需要计算控制图的中心限和上下控制界限,设欧氏距离d为正态总体样本X,X~N(μ,σ),x为取自总体样本X的样本,标准差S作为样本统计量,由《均值及标准偏差控制图及其应用》查得标准差图控制限计算公式如下:
其中,n为样本容量;c4为与样本大小有关的常数,可以由《均值及标准偏差控制图及其应用》附表查得;S为样本的标准差,
在控制图中由点出界的个数来进行损伤判别,连续100个点中有3个或者多于3个点超出控制界限就可以判定结构发生损伤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310690423.6A CN103646177B (zh) | 2013-12-16 | 2013-12-16 | 一种考虑温度影响的土木工程结构损伤预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310690423.6A CN103646177B (zh) | 2013-12-16 | 2013-12-16 | 一种考虑温度影响的土木工程结构损伤预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103646177A CN103646177A (zh) | 2014-03-19 |
CN103646177B true CN103646177B (zh) | 2016-08-17 |
Family
ID=50251390
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310690423.6A Active CN103646177B (zh) | 2013-12-16 | 2013-12-16 | 一种考虑温度影响的土木工程结构损伤预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103646177B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297004A (zh) * | 2014-09-18 | 2015-01-21 | 天津大学 | 基于ar-arx模型的桥梁实时损伤预警方法 |
CN106338372A (zh) * | 2016-09-19 | 2017-01-18 | 中国海洋大学 | 基于残差应变能的海洋平台损伤定位方法及系统 |
CN108197603A (zh) * | 2018-01-30 | 2018-06-22 | 青岛理工大学 | 考虑环境因素影响的基于奇异谱分析的结构损伤预警方法 |
CN110544311B (zh) * | 2018-05-29 | 2023-04-25 | 百度在线网络技术(北京)有限公司 | 安全警示方法、装置及存储介质 |
CN109840386B (zh) * | 2019-03-05 | 2022-09-23 | 中国人民解放军海军勤务学院 | 基于因子分析的损伤识别方法 |
CN109993211B (zh) * | 2019-03-05 | 2022-09-23 | 中国人民解放军海军勤务学院 | 基于结构加速度ar模型系数的损伤识别方法 |
CN110795693A (zh) * | 2019-09-24 | 2020-02-14 | 集美大学 | 一种噪声监测数据自动校正方法、装置及计算机可读存储介质 |
CN114664058B (zh) * | 2022-01-29 | 2023-08-18 | 上海至冕伟业科技有限公司 | 一种消防水系统的整体故障预警系统及方法 |
CN117216844B (zh) * | 2023-09-12 | 2024-03-26 | 汕头大学 | 一种桥梁结构损伤检测方法、系统和存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101221104A (zh) * | 2007-10-16 | 2008-07-16 | 吴智深 | 基于分布式应变动态测试的结构健康监测方法 |
CN101561379A (zh) * | 2009-05-13 | 2009-10-21 | 清华大学 | 一种用于结构损伤检测的敲击扫描方法 |
US7853339B2 (en) * | 2006-09-29 | 2010-12-14 | Fisher-Rosemount Systems, Inc. | Statistical signatures used with multivariate analysis for steady-state detection in a process |
CN101923027A (zh) * | 2010-05-28 | 2010-12-22 | 清华大学 | 一种结构损伤检测系统、设备以及结构损伤检测方法 |
CN102183524A (zh) * | 2011-01-10 | 2011-09-14 | 哈尔滨工业大学 | 用于土木工程结构表观缺陷评定的双ccd检测方法及系统 |
-
2013
- 2013-12-16 CN CN201310690423.6A patent/CN103646177B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7853339B2 (en) * | 2006-09-29 | 2010-12-14 | Fisher-Rosemount Systems, Inc. | Statistical signatures used with multivariate analysis for steady-state detection in a process |
CN101221104A (zh) * | 2007-10-16 | 2008-07-16 | 吴智深 | 基于分布式应变动态测试的结构健康监测方法 |
CN101561379A (zh) * | 2009-05-13 | 2009-10-21 | 清华大学 | 一种用于结构损伤检测的敲击扫描方法 |
CN101923027A (zh) * | 2010-05-28 | 2010-12-22 | 清华大学 | 一种结构损伤检测系统、设备以及结构损伤检测方法 |
CN102183524A (zh) * | 2011-01-10 | 2011-09-14 | 哈尔滨工业大学 | 用于土木工程结构表观缺陷评定的双ccd检测方法及系统 |
Non-Patent Citations (2)
Title |
---|
PCA-based detection of damage in time-varying systems;A.Bellino等;《Mechanical Systems and Signal Processing》;20101031;第24卷(第7期);2250-2260 * |
基于控制图法的结构损伤识别方法研究;纪志远;《中国优秀博硕学位论文全文数据库-工程科技Ⅱ辑》;20060115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103646177A (zh) | 2014-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103646177B (zh) | 一种考虑温度影响的土木工程结构损伤预警方法 | |
CN110175434B (zh) | 一种基于卷积神经网络的铁路扣件系统损伤检测方法 | |
CN104712542B (zh) | 一种基于物联网的往复压缩机敏感特征提取与故障诊断方法 | |
CN102663264B (zh) | 桥梁结构健康监测静态参数的半监督协同评估方法 | |
CN106555788A (zh) | 基于模糊处理的深度学习在液压装备故障诊断中的应用 | |
CN112149316A (zh) | 基于改进的cnn模型的航空发动机剩余寿命预测方法 | |
CN105629958B (zh) | 一种基于子时段mpca‑svm的间歇过程故障诊断方法 | |
CN109558635B (zh) | 一种基于单元模态应变能灵敏度的结构区间不确定性损伤识别方法 | |
CN109115879B (zh) | 一种基于模态振型和卷积神经网络的结构损伤识别方法 | |
CN108629369A (zh) | 一种基于Trimmed SSD的尿沉渣有形成分自动识别方法 | |
CN102829967A (zh) | 一种基于回归模型系数变化的时域故障识别方法 | |
CN103884776B (zh) | 一种提高随机损伤定位向量法监测结果准确性的方法 | |
CN111595541A (zh) | 一种海量振动传递率数据卷积神经网络处理的多维结构损伤识别方法 | |
CN116738868B (zh) | 一种滚动轴承剩余寿命预测方法 | |
CN108717548A (zh) | 一种面向传感器动态增加的行为识别模型更新方法及系统 | |
CN116168720A (zh) | 一种电机声音异常故障诊断方法、系统、及可存储介质 | |
CN105203327B (zh) | 一种应用于发动机气路分析的气路测量参数选择方法 | |
CN102620980B (zh) | 一种利用神经网络预报板材成形性能的方法 | |
CN112418334B (zh) | 一种铁路动车裙板格栅变形故障识别方法 | |
CN116957361B (zh) | 一种基于虚实结合的船舶任务系统健康状态检测方法 | |
CN104318079B (zh) | 一种基于故障演化分析的故障预测特征选择方法 | |
Khan et al. | Continuous modal parameter identification of cable-stayed bridges based on a novel improved ensemble empirical mode decomposition | |
CN116757085A (zh) | 基于残差连接卷积循环网络的冲击荷载作用位置反演方法 | |
CN114118864B (zh) | 一种大型群众体育赛事综合评估系统 | |
CN114580239B (zh) | 一种考虑不确定性的桥梁损伤识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |