基于统计过程控制的土木工程结构损伤检测方法与流程

本发明涉及一种土木工程结构损伤检测方法，特别是涉及一种基于统计过程控制的土木工程结构损伤检测方法。

背景技术：

一般来说，结构健康监测中的相关方法可分为5个层次：(1)确定结构损伤是否发生，可称为“损伤检测方法”，即本发明所涉及的方法；(2)确定结构损伤位置的方法；(3)确定结构损伤类型的方法；(4)确定结构损伤程度的方法；(5)确定结构剩余寿命的方法。现有结构健康监测领域的相关文献表明，从层次(1)至层次(5)，方法实施的难度依次增大。

对于土木工程领域的实际结构，即使是上述层次(1)的方法实施难度也较大。现有技术主要有两种，一是有限元模型方法，通过建立结构的物理模型，描述监测数据与结构参数之间的函数关系，从而根据监测数据的变化反演结构参数的变化，来确定结构损伤的发生；二是模态分析方法，将原始监测数据从时域转化为模态域，根据模态数据变化确定结构损伤的发生。

但是这两种现有方法都存在着许多不足，具体来讲，有限元模型方法旨在模拟实际结构中各微元体几何、力学、物理等方面的关系，因此随着结构复杂程度的增加，有限元模型建模难度极大提高，以致于通常难以建立合理的模型，模态分析方法需要从原始时域监测数据抽取出模态域特征量，这种处理方法必然会减少数据中的信息量。

技术实现要素：

发明目的：本发明要解决的技术问题是提供一种基于统计过程控制的土木工程结构损伤检测方法，其建模难度不会随结构复杂程度的增加而增加，并且它直接采用原始监测数据进行结构损伤检测，避免了数据信息量的损失，从而提高了方法的有效性。

技术方案：本发明所述的基于统计过程控制的土木工程结构损伤检测方法,包括以下步骤：

(1)利用数据监测设备采集土木工程的结构响应数据xk，共m×n个，1≤k≤m×n，并把xk分为m组，每组数据的个数都是n；

(2)对每组数据采用下式求均值与极差ri：

(3)均值极差ri的上下限分别按下式计算：

其中

其中，分别为的上限、中心线、下限，uclr、centerr、lclr分别为ri的上限、中心线、下限，a2、d3、d4是与n有关的统计参数；

(4)检查在之外的数据个数是否超过m·p1，检查在[uclr,lclr]之外ri的数据个数是否超过m·p2，其中p1是当结构处于正常状态下时在之外的概率，p2是当结构处于正常状态下时ri在[uclr,lclr]之外的概率，如果上述二者之一超过，即表明结构处于损伤状态，完成对结构损伤的检测。

所述的结构响应数据xk包括加速度、速度、位移和应变，且在相同时间间隔采集。

为了在结构损伤检测的有效性，uclr、lclr的设置原则为使所述的p1在0.27％至1％范围内，使所述的p2在0.27％至1％范围内。

进一步的，n＝4、a2＝0.729、d3＝0、d4＝2.282，或者n＝5、a2＝0.577、d3＝0、d4＝2.114、或者n＝6、a2＝0.483、d3＝0、d4＝2.004。

有益效果：本方法能够采用数据模型对土木工程结构进行损伤检测计算，模型复杂程度与实际结构没有直接关系，可较为方便地建立合理模型，并且本方法采用原始时域监测数据进行计算，避免了时域数据向频域数据转换过程中信息量的损失。使用本方法，能够准确计算得出国际结构控制协会与美国土木工程协会设计的基准结构全部工况下的损伤状态。

附图说明

图1是benchmark结构的构型图；

图2是本方法的流程图；

图3是无损状态下的检测结果数据图；

图4是有损状态1下的检测结果数据图。

具体实施方式

国际结构控制协会(internationalassociationforstructuralcontrol，简写为iasc)、美国土木工程协会(americansocietyofcivilengineers，简写为asce)提出了一个如图1所示的benchmark结构，作为比较各种结构损伤检测方法的试验研究平台，方便研究者们将不同方法应用于此结构，以评价所采用方法的优劣。

原始监测数据xk是在外部荷载的作用下，土木工程结构响应的采集数据。其中，外部荷载一般来源于地震、风等自然因素，若结构体量较小，也可使用人工荷载；结构响应包括加速度、速度、位移、应变等不同类型；数据采集一般需要利用数据监测设备，在相同时间间隔的时刻测得结构响应。本实施例中采用的原始监测数据均为matlab程序计算得到的模拟监测数据，原始监测数据均包含了噪信比为10％的计算误差，损伤情况采用结构中的斜支撑进行模拟，即：结构各层的8*4＝32根斜支撑全部存在的状态，对应于无损状态；去掉不同层的不同数量的斜支撑，对应于不同的损伤状态。

如图2所示，具体的检测方法是，第一步利用数据监测设备采集土木工程的结构响应数据xk，共m×n个，1≤k≤m×n，并把xk分为m组，每组数据的个数都是n；第二步对每组数据求出均值与极差ri，i＝1,2,…,m；第三步求出均值极差ri的上下限uclr1和lclr1；第四步对未知状态下的x’h按照第一步和第二步操作，1≤h≤m×n，得到和rj’；第五步将第三步所得上下限和rj’画在二维图上；第六步检查在之外的数据个数，如果超过m·0.27％则检测出结构处于损伤状态；第七步检查在[lclr1,uclr1]之外rj’的数据个数，如果超过m·0.27％则检测出结构处于损伤状态；如果第六步和第七步都没有超过m·0.27％则检测出结构处于无损状态。

首先在无损状态下进行检验。与uclr1、clr1、lclr1由无损状态下的数据计算得到，计算采用的监测数据为第一层顶板处的结构加速度响应，数据个数取为2000，把原始测试数据分为500组，每组数据的个数是4，a2＝0.729、d3＝0、d4＝2.282，p1＝p2＝0.27％。利用原始数据进行检测计算，均值极差ri、上限ucl、中心线center、下限lcl的数据图如图3所示，其中，上、下两条水平线分别为ucl、lcl，圆圈表示ri的超限数据。在正常状态下ri的最大超限数据个数可认为是500*0.27％＝1.35，在图3中ri的实际超限数据的个数分别为1和1，均小于500*0.27％＝1.35，即本方法正确指出了结构处于无损状态下。

第二在有损状态下进行检测。定义有损状态1为去掉第一层的所有斜支撑的状态，与无损状态下的情况相同，有损状态1下均值极差ri的上限ucl、中心线center、下限lcl采用与uclr1、clr1、lclr1表示，监测数据也采用第一层顶板处的结构加速度响应，采样频率为1000，数据个数取为2000，把原始测试数据分为500组，每组数据的个数是4，a2＝0.729、d3＝0、d4＝2.282，p1＝p2＝0.27％。利用原始数据进行检测计算，均值极差ri、上限ucl、中心线center、下限lcl的数据图如图4所示，其中，上、下两条水平线分别为ucl、lcl，圆圈表示ri的超限数据。图4中ri的实际超限数据的个数显然均大于正常状态下最大超限数据个数500*0.27％＝1.35，即本方法正确指出结构处于有损状态下。