基于改进小波包能量曲率的斜拉桥损伤识别方法及装置

1.本发明涉及桥梁技术领域，尤其涉及一种基于改进小波包能量曲率的斜拉桥损伤识别方法及装置。


背景技术：

2.我国作为世界上桥梁数量最多的国家，面临着巨大的桥梁运营管理任务，为了保障运营期间桥梁结构的安全性和良好的使用性能，需要采用桥梁健康监测系统及时准确的监测桥梁的健康状态。桥梁健康监测过程中，损伤识别是一个重要的环节。小波分析具有多分辨率分析的优点，而且在时频两域都具有表征信号特征的能力，这个能力对于损伤识别特别重要。
3.小波包分析能够为信号提供一种更加精细的分析方法，对低频和高频部分都能进一步分解。小波包分解具有很强的局部化分解能力，为提取和识别损伤类型提供了保障。因此，利用小波包分析进行斜拉桥的损伤识别具有重要意义。
4.现有技术中，基于小波包能量曲率和基于小波包能量变化率指标都能有效识别定位损伤，具有很强的损伤识别敏感性，但目前只是分别使用了两种指标，损伤识别精度有进一步提高的潜力。


技术实现要素：

5.本发明实施例提供了一种基于改进小波包能量曲率的斜拉桥损伤识别方法及装置，以解决当前损伤识别中损伤识别精度需要进一步提高的问题。
6.第一方面，本发明实施例提供了一种基于改进小波包能量曲率的斜拉桥损伤识别方法，包括：
7.获取斜拉桥的时频响应信号；
8.对所述时频响应信号进行小波包分解与重构，计算小波包能量曲率；
9.根据所述小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值；
10.根据所述损伤指标值诊断所述斜拉桥的损伤。
11.在一种可能的实现方式中，所述根据所述小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值，包括：
12.根据所述小波包能量曲率，计算小波包能量曲率平方和；
13.根据所述小波包能量曲率平方和，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值。
14.在一种可能的实现方式中，所述根据所述小波包能量曲率，计算小波包能量曲率平方和，包括：
15.根据计算所述小波包能量曲率平方和；
16.其中，wpeks表示所述小波包能量曲率平方和；ki表示所述时频响应信号经小波包分解成不同频段后，第i频段的小波包能量曲率；所述i表示频率参数，i＝0，1，
…
，2
j-1；所述j表示尺度参数，所述j为大于0的整数。
17.在一种可能的实现方式中，所述根据所述小波包能量曲率平方和，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值，包括：
18.根据计算小波包能量曲率平方和的相对变化率；
19.其中，wpeksr表示所述小波包能量曲率平方和的相对变化率；wpeksd表示斜拉桥损伤后的小波包能量曲率平方和；wpeksu表示斜拉桥损伤前的小波包能量曲率平方和；
20.将计算得到的wpeksr，确定为所述损伤指标值。
21.在一种可能的实现方式中，所述根据所述损伤指标值诊断所述斜拉桥的损伤，包括：
22.根据所述斜拉桥建立对应的模型；
23.在所述模型中，设置不同损伤程度工况，模拟出不同损伤程度工况下所述模型中各测点对应的损伤指标值；
24.根据所述不同损伤程度工况下对应的各测点的损伤指标值，确定损伤位置和所述损伤位置对应的损伤指标值；
25.根据所述损伤位置对应的损伤指标值和对应损伤程度工况下的损伤程度值，计算得到所述斜拉桥的损伤程度与损伤指标值的函数对应关系；
26.根据，基于所述斜拉桥的实测时频响应信号计算得到的损伤指标值，和所述函数对应关系，诊断所述斜拉桥的损伤。
27.在一种可能的实现方式中，所述对所述时频响应信号进行小波包分解与重构，计算小波包能量曲率，包括：
28.对所述时频响应信号进行小波包分解与重构，得到小波包能量；
29.根据所述小波包能量，计算所述小波包能量曲率。
30.在一种可能的实现方式中，，所述根据所述小波包能量，计算所述小波包能量曲率，包括：
31.根据计算所述小波包能量曲率；
32.其中，ki表示所述时频响应信号经小波包分解成不同频段后，第i频段的小波包能量曲率；表示所述时频响应信号经j层小波包分解后第i频段内的第k+1个测点的小波包能量；表示所述时频响应信号经j层小波包分解后第i频段内的第k个测点的小波包能量；表示所述时频响应信号经j层小波包分解后第i频段内的第k-1个测点的小波包能量；l表示测点间距；所述i表示频率参数，i＝0，1，
…
，2
j-1；所述j表示尺度参数，所述j为大于0的整数；所述k表示位置参数，所述k为大于0的整数。
33.第二方面，本发明实施例提供了一种基于改进小波包能量曲率的斜拉桥损伤识别装置，包括：
34.获取模块，用于获取斜拉桥的时频响应信号；
35.计算模块，用于对所述时频响应信号进行小波包分解与重构，计算小波包能量曲率；
36.所述计算模块，还用于根据所述小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值；
37.损伤识别模块，用于根据所述损伤指标值诊断所述斜拉桥的损伤。
38.第三方面，本发明实施例提供了一种终端，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上第一方面或第一方面的任一种可能的实现方式所述方法的步骤。
39.第四方面，本发明实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上第一方面或第一方面的任一种可能的实现方式所述方法的步骤。
40.本发明实施例提供一种基于改进小波包能量曲率的斜拉桥损伤识别方法及装置，通过获取斜拉桥的时频响应信号；对时频响应信号进行小波包分解与重构，计算小波包能量曲率；根据小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值；根据损伤指标值诊断斜拉桥的损伤。在对斜拉桥进行损伤识别时，将小波包能量曲率和小波包能量变化率两种指标相结合，提出了一种新的损伤指标—小波包能量曲率平方和的相对变化率。充分利用了小波包能量曲率和小波包能量变化率两种指标的敏感特性，能够有效提高损伤识别精度。
附图说明
41.为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
42.图1是本发明实施例提供的基于改进小波包能量曲率的斜拉桥损伤识别方法的实现流程图；
43.图2是本发明实施例中设置的损伤位置和测点的示意图；
44.图3是本发明实施例中设置的不同工况的损伤位置及损伤程度的汇总图；
45.图4是本发明实施例提供的不同工况下各测点损伤指标的数据汇总图；
46.图5(a)是本发明实施例提供的算例分析中，工况2对应的各测点损伤指标的点线图；
47.图5(b)是本发明实施例提供的算例分析中，工况3对应的各测点损伤指标的点线图；
48.图5(c)是本发明实施例提供的算例分析中，工况4对应的各测点损伤指标的点线图；
49.图5(d)是本发明实施例提供的算例分析中，工况5对应的各测点损伤指标的点线
图；
50.图5(e)是本发明实施例提供的算例分析中，工况6对应的各测点损伤指标的点线图；
51.图6是本发明实施例提供的不同工况下各种损伤指标的的相对差值对比图；
52.图7是本发明实施例提供的基于改进小波包能量曲率的斜拉桥损伤识别装置的结构示意图；
53.图8是本发明实施例提供的终端的示意图。
具体实施方式
54.以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本发明实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本发明的描述。
55.为使本发明的目的、技术方案和优点更加清楚，下面将结合附图通过具体实施例来进行说明。
56.图1为本发明实施例提供的基于改进小波包能量曲率的斜拉桥损伤识别方法的实现流程图，详述如下：
57.步骤101，获取斜拉桥的时频响应信号。
58.可选的，这里的时频响应信号可以是加速度响应信号，也可以是其他动力响应信号，此处不做具体限定。
59.步骤102，对时频响应信号进行小波包分解与重构，计算小波包能量曲率。
60.可选的，对时频响应信号进行小波包分解与重构，计算小波包能量曲率，包括：
61.步骤201，对时频响应信号进行小波包分解与重构，得到小波包能量。
62.小波包分解能对信号的振动频带进行多层次的划分，分解后的每一个测点均代表一个小波包系数。
63.小波包系数可以表示为其中，表示信号经j层小波包分解后，第i频段内第k个测点的小波包系数；表示小波包函数；s(t)表示时频响应信号；i表示频率参数，i＝0，1，
…
，2
j-1；j表示尺度(层数)参数，k表示位置参数，j、k均为大于0的整数。
64.需要说明的是，小波包函数具有正交性，也就是，对于经小波包分解后，同一层不同频段间的小波包函数满足：同频段间的小波包函数满足：其中，m、n分别表示不同频段。
65.对于任一时频响应信号s(t)，经过j层小波包分解后，可以根据(m≠n)计算小波包能量。
66.其中，ej表示时频响应信号s(t)经过j层小波包分解后的小波包能量；表示
时频响应信号s(t)经过j层小波包分解后，第m频段对应的小波包组分信号；表示时频响应信号s(t)经过j层小波包分解后，第n频段对应的小波包组分信号。
67.进一步地，根据计算各频段对应的小波包组分信号。其中，表示时频响应信号s(t)经过j层小波包分解后，第i频段对应的小波包组分信号。
68.考虑到小波包函数具有正交性，因此，上述计算小波包能量的公式可以简化为：由此，可以得到：其中，表示时频响应信号s(t)经过j层小波包分解后，第i个频段内的小波包能量。
69.步骤202，根据小波包能量，计算小波包能量曲率。
70.为计算任一频段的小波包能量曲率，首先需要计算当前频段内各个测点的小波包能量。
71.将位置参数k带入公式得到其中，表示时频响应信号s(t)经j层小波包分解后，第i频段内第k个测点的小波包能量，表示时频响应信号s(t)经j层小波包分解后，第i频段内第k个测点的信号分量。
72.进一步地，可以根据计算测点间距不等时的小波包能量曲率。其中，ki表示时频响应信号经小波包分解成不同频段后，第i频段的小波包能量曲率；表示时频响应信号经j层小波包分解后第i频段内的第k+1个测点的小波包能量；表示时频响应信号经j层小波包分解后第i频段内的第k个测点的小波包能量；表示时频响应信号经j层小波包分解后第i频段内的第k-1个测点的小波包能量；l
k+1
为第k+1个测点到第k个测点的测点间距；l
k-1
为第k-1个测点到第k个测点的测点间距；i表示频率参数，i＝0，1，
…
，2
j-1，j表示尺度参数，j为大于0的整数；k表示位置参数，k为大于0的整数。
73.在此基础上，设定测点间距相等，且测点间距均为l。则可以根据计算小波包能量曲率。
74.步骤103，根据小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定
为损伤指标值。
75.可选的，根据小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值，包括：
76.步骤301，根据小波包能量曲率，计算小波包能量曲率平方和。
77.进一步地，可以根据计算小波包能量曲率平方和；
78.其中，wpeks表示小波包能量曲率平方和；ki表示时频响应信号经小波包分解成不同频段后，第i频段的小波包能量曲率；i表示频率参数，i＝0，1，
…
，2
j-1，j表示尺度参数，j为大于0的整数。
79.步骤302，根据小波包能量曲率平方和，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值。
80.进一步地，可以根据计算小波包能量曲率平方和的相对变化率；将计算得到的wpeksr，确定为损伤指标值。
81.其中，wpeksr表示小波包能量曲率平方和的相对变化率，也就是本发明提出的新的损伤指标，wpeksd表示斜拉桥损伤后的小波包能量曲率平方和，wpeksu表示斜拉桥损伤前的小波包能量曲率平方和。这里的wpeksu可以从斜拉桥配备的数据库中获取。
82.步骤104，根据损伤指标值诊断斜拉桥的损伤。
83.可选的，根据损伤指标诊断斜拉桥的损伤，包括：
84.步骤401，根据斜拉桥建立对应的模型；
85.进一步地，可以根据斜拉桥建立基准有限元模型。
86.步骤402，在模型中，设置不同损伤程度工况，模拟出不同损伤程度工况下模型中各测点对应的损伤指标值；
87.在基准有限元模型中，设置不同的工况及测点，不同工况对应不同的损伤位置以及不同的损伤程度。
88.进一步地，针对不同工况，提取各个测点的时频响应信号，并根据上述计算损伤指标值的方法计算在不同工况下各个测点对应的损伤指标值。
89.步骤403，根据不同损伤程度工况下对应的各测点的损伤指标值，确定损伤位置和损伤位置对应的损伤指标值；
90.根据不同损伤程度工况下对应的各测点的损伤指标值可以确定损伤位置。具体地，在任一工况下，损伤指标值大于预设值的测点位置即确定为损伤位置，该测点对应的损伤指标值即为损伤位置对应的损伤指标值。这里的预设值，可以自行设定，在这里不限定具体数值。
91.进一步地，为更加直观准确的确定损伤位置，可以根据上述不同损伤程度工况下对应的各测点的损伤指标值建立测点位置和损伤指标值的点线图，损伤指标值出现明显峰值的测点位置即为损伤位置，该测点对应的损伤指标值即为损伤位置对应的损伤指标值。具体的点线图可以参考下面算例描述部分。
92.步骤404，根据损伤位置对应的损伤指标值和对应损伤程度工况下的损伤程度值，计算得到斜拉桥的损伤程度与损伤指标值的函数对应关系；
93.也就是说，根据上述得到的不同损伤程度工况下，损伤位置对应的损伤指标值，在matlab中进行拟合计算，得到斜拉桥中损伤位置的损伤程度与损伤指标值的函数对应关系。
94.步骤405，根据，基于斜拉桥的实测时频响应信号计算得到的损伤指标值，和函数对应关系，诊断斜拉桥的损伤。
95.实测斜拉桥的时频响应信号，并进一步计算各测点的损伤指标值，可以确定各测点位置是否为损伤位置。进一步地，将损伤位置的损伤指标值带入上述函数对应关系，对斜拉桥中损伤位置的实际损伤进行诊断。
96.下面结合一具体算例来说明该基于改进小波包能量曲率的斜拉桥损伤识别方法的可行性：
97.选取一座全长98.5米，总宽7.2米的斜拉桥建立基准有限元模型，斜拉桥主梁上共布置23个测点，主梁上设置两处80mm的损伤区域a和b，其中a中心在测点5处，b中心在测点14处，如图2所示。
98.如图3示，损伤程度设定10％、15％和20％三种情况，再结合两个损伤区域a和b，共设置6中损伤工况。
99.在不同的损伤工况下，分别提取各测点的加速度响应数据，并根据加速度响应数据计算各测点的损伤指标值。其中，不同工况下的各测点损伤指标值，请参考图4。
100.需要说明的是，本算例中计算损伤指标值时，经过循环试算，最终采用db20小波分解9层。但由于小波包分解层数对损伤指标值wpeksr识别损伤位置及损伤程度的影响较大，且对单损伤与多损伤的影响不同；小波包基对损伤指标值wpeksr识别损伤位置及损伤程度的影响不大，仅对多损伤工况下的小损伤识别程度有影响。具体应用时，需要先找到合适的小波包层数和小波包基。
101.参见图5(a)～图5(e)，根据图4中不同工况下的损伤指标值，建立点线图，图5(a)、图5(b)和图5(c)中可以直观看到5测点相比其他测点出现了明显峰值，说明损伤位置发生在5测点；图5(d)和图5(e)中可以直观看到5测点与14测点出现明显峰值，说明损伤发生在5测点和14测点，与预设的损伤位置一致。
102.由此，可以得出，根据损伤指标值wpeksr可以直观准确的识别出损伤位置。
103.进一步地，工况2中5测点的损伤指标值为0.09023，工况3中5测点的损伤指标值为0.14477，工况5中5测点与14测点的损伤指标值分别为0.09058、0.08863，工况6中5测点与14测点的损伤指标值分别为0.20824、0.20881，结合这四种工况的损伤程度与损伤位置的损伤指标值在matlab中拟合计算得出函数：y＝-153.07x3+66.7547x
2-8.3167x+0.4194，其中，y表示损伤程度，x表示损伤位置的损伤指标值wpeksr。利用上述函数关系求得工况4的损伤程度为0.2003817，与实际设定损伤程度的误差为0.19％。可见，上述函数关系能对斜拉桥结构的损伤程度进行准确识别。
104.进一步地，针对同一工况，采用小波包能量曲率差和小波包能量的相对变化率两种已有的损伤指标对基准有限元模型上的加速度响应信号进行分析，同样采用db20小波9层分解，并将两种指标的计算结果同本方法中的损伤指标wpeksr作比较。取同一损伤位置预设损伤程度比值为0.5的两工况与其对应损伤指标大小比值的相对差值来比较各损伤指标的识别精度。
105.其中，根据计算相对差值。其中，rer表示相对差值；q表示两工况对应的同一损伤位置的损伤指标比值。指标rer越接近于0，则说明其对应的损伤指标对损伤程度的识别效果越好，各损伤指标的rer值结果如表6所示。
106.参见图6可知，wpeksr的rer值小于其他损伤指标的rer值，说明wpeksr指标对于损伤识别的精度高于其他两项损伤指标。
107.本发明实施例通过获取斜拉桥的时频响应信号；对时频响应信号进行小波包分解与重构，计算小波包能量曲率；根据小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值；根据损伤指标值确定斜拉桥的损伤。在对斜拉桥进行损伤识别时，将小波包能量曲率和小波包能量变化率两种指标相结合，提出了一种新的损伤指标—小波包能量曲率平方和的相对变化率。充分利用了小波包能量曲率和小波包能量变化率两种指标的敏感特性，能够有效提高损伤识别精度。
108.应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。
109.以下为本发明的装置实施例，对于其中未详尽描述的细节，可以参考上述对应的方法实施例。
110.图7示出了本发明实施例提供的基于改进小波包能量曲率的斜拉桥损伤识别装置的结构示意图，为了便于说明，仅示出了与本发明实施例相关的部分，详述如下：
111.如图7所示，基于改进小波包能量曲率的斜拉桥损伤识别装置7包括：获取模块71、计算模块72和损伤识别模块73。
112.获取模块71，用于获取斜拉桥的时频响应信号。
113.这里的时频响应信号可以是加速度响应信号，也可以是其他动力响应信号，此处不做具体限定。
114.在一种可能的实现方式中，获取模块71，用于获取斜拉桥的加速度响应信号。
115.计算模块72，用于对时频响应信号进行小波包分解与重构，计算小波包能量曲率。
116.在一种可能的实现方式中，计算模块72，用于对时频响应信号进行小波包分解与重构，得到小波包能量。
117.计算模块72，用于对时频响应信号进行小波包分解与重构；
118.计算模块72，还用于根据计算时频响应信号s(t)经过j层小波包分解后，第i个频段内的小波包能量。其中，表示时频响应信号s(t)经过j层小波包分解后，第i个频段内的小波包能量；表示时频响应信号s(t)经过j层小波包分解后，第i频段对应的小波包组分信号；i表示频率参数，i＝0，1，
…
，2
j-1，j表示尺度(层数)参数，j为大于0的整数。
119.具体地，对于任一时频响应信号s(t)，经过j层小波包分解后，计算模块72，用于根据(m≠n)计算小波包能量。其中，ej表
示时频响应信号s(t)经过j层小波包分解后的小波包能量；表示时频响应信号s(t)经过j层小波包分解后，第m频段对应的小波包组分信号；表示时频响应信号s(t)经过j层小波包分解后，第n频段对应的小波包组分信号；m＝0，1，
…
，2
j-1，n＝0，1，
…
，2
j-1，m≠n；m、n分别表示不同频段。
120.考虑到小波包函数的正交特性，上述公式可以简化为由此得到，其中，表示时频响应信号s(t)经过j层小波包分解后，第i个频段内的小波包能量；表示时频响应信号s(t)经过j层小波包分解后，第i频段对应的小波包组分信号；i表示频率参数，i＝0，1，
…
，2
j-1，j表示尺度(层数)参数，j为大于0的整数。
121.计算模块72，还用于根据小波包能量，计算小波包能量曲率。
122.为计算任一频段的小波包能量曲率，计算模块72，还用于计算当前频段内各个测点的小波包能量。
123.将位置参数k带入上述公式得到其中，表示时频响应信号s(t)经j层小波包分解后，第i频段内第k个测点的小波包能量，表示时频响应信号s(t)经j层小波包分解后，第i频段内第k个测点的信号分量。
124.也就是说，计算模块72，用于根据计算任一频段内各个测点的小波包能量。
125.进一步地，计算模块72，还用于根据计算测点间距不等时的小波包能量曲率。其中，ki表示时频响应信号经小波包分解成不同频段后，第i频段的小波包能量曲率；表示时频响应信号经j层小波包分解后第i频段内的第k+1个测点的小波包能量；表示时频响应信号经j层小波包分解后第i频段内的第k个测点的小波包能量；表示时频响应信号经j层小波包分解后第i频段内的第k-1个测点的小波包能量；l
k+1
为第k+1个测点到第k个测点的测点间距；l
k-1
为第k-1个测点到第k个测点的测点间距；i表示频率参数，i＝0，1，
…
，2
j-1，j表示尺度参数，j为大于0的整数；k表示位置参数，k为大于0的整数。
126.在此基础上，设定测点间距相等，且测点间距均为l。则计算模块72，还用于根据
计算小波包能量曲率。
127.计算模块72，还用于根据小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值。
128.在一种可能的实现方式中，计算模块72，用于根据小波包能量曲率，计算小波包能量曲率平方和。
129.进一步地，计算模块72，用于根据计算小波包能量曲率平方和；
130.其中，wpeks表示小波包能量曲率平方和；ki表示时频响应信号经小波包分解成不同频段后，第i频段的小波包能量曲率；i表示频率参数，i＝0，1，
…
，2
j-1，j表示尺度参数，j为大于0的整数。
131.计算模块72，还用于根据小波包能量曲率平方和，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值。
132.进一步地，计算模块72，用于根据计算小波包能量曲率平方和的相对变化率；将计算得到的wpeksr，确定为损伤指标值；
133.其中，wpeksr表示小波包能量曲率平方和的相对变化率，也就是本发明提出的新的损伤指标，wpeksd表示斜拉桥损伤后的小波包能量曲率平方和，wpeksu表示斜拉桥损伤前的小波包能量曲率平方和，这里的wpeksu可以从斜拉桥配备的数据库中获取。
134.损伤识别模块73，用于根据损伤指标值诊断斜拉桥的损伤。
135.在一种可能的实现方式中，损伤识别模块73，用于根据斜拉桥建立对应的模型；
136.进一步地，损伤识别模块73，用于根据斜拉桥建立基准有限元模型。
137.损伤识别模块73，还用于在模型中，设置不同损伤程度工况，模拟出不同损伤程度工况下模型中各测点对应的损伤指标值；
138.进一步地，损伤识别模块73，用于在基准有限元模型中，设置不同的工况及测点，不同工况对应不同的损伤位置以及不同的损伤程度；同时，损伤识别模块73，还用于针对不同工况，提取各个测点的时频响应信号，并根据上述计算损伤指标值的方法计算在不同工况下各个测点对应的损伤指标值。
139.损伤识别模块73，还用于根据不同损伤程度工况下对应的各测点的损伤指标值，确定损伤位置和损伤位置对应的损伤指标值；
140.进一步地，损伤识别模块73，用于根据上述不同损伤程度工况下对应的各测点的损伤指标值建立测点位置和损伤指标值的点线图，损伤指标值出现明显峰值的测点位置即为损伤位置，该测点对应的损伤指标值即为损伤位置对应的损伤指标值。
141.损伤识别模块73，还用于根据损伤位置对应的损伤指标值和对应损伤程度工况下的损伤程度值，计算得到斜拉桥的损伤程度与损伤指标值的函数对应关系；
142.进一步地，损伤识别模块73，用于根据上述得到的不同损伤程度工况下，损伤位置对应的损伤指标值，在matlab中进行拟合计算，得到斜拉桥中损伤位置的损伤程度与损伤
指标值的函数对应关系。
143.损伤识别模块73，还用于根据，基于斜拉桥的实测时频响应信号计算得到的损伤指标值，和函数对应关系，诊断斜拉桥的损伤。
144.进一步地，损伤识别模块73，还用于实测斜拉桥的时频响应信号，并进一步计算各测点的损伤指标值，可以确定各测点位置是否为损伤位置。进一步地，将损伤位置的损伤指标值带入上述函数对应关系，对斜拉桥中损伤位置的实际损伤进行诊断。
145.本发明实施例通过获取模块，用于获取斜拉桥的时频响应信号；计算模块，用于对时频响应信号进行小波包分解与重构，计算小波包能量曲率；计算模块，还用于根据小波包能量曲率，计算小波包能量曲率平方和的相对变化率，确定为损伤指标值；损伤识别模块，用于根据损伤指标值确定斜拉桥的损伤。在对斜拉桥进行损伤识别时，将小波包能量曲率和小波包能量变化率两种指标相结合，提出了一种新的损伤指标—小波包能量曲率平方和的相对变化率。充分利用了小波包能量曲率和小波包能量变化率两种指标的敏感特性，能够有效提高损伤识别精度。
146.图8是本发明实施例提供的终端的示意图。如图8所示，该实施例的终端8包括：处理器80、存储器81以及存储在所述存储器81中并可在所述处理器80上运行的计算机程序82。所述处理器80执行所述计算机程序82时实现上述各个基于改进小波包能量曲率的斜拉桥损伤识别方法实施例中的步骤，例如图1所示的步骤101至步骤104。或者，所述处理器80执行所述计算机程序82时实现上述各装置实施例中各模块的功能，例如图7所示模块71至73的功能。
147.示例性的，所述计算机程序82可以被分割成一个或多个模块，所述一个或者多个模块被存储在所述存储器81中，并由所述处理器80执行，以完成本发明。所述一个或多个模块可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序82在所述终端8中的执行过程。例如，所述计算机程序82可以被分割成图7所示的模块71至73。
148.所述终端8可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述终端8可包括，但不仅限于，处理器80、存储器81。本领域技术人员可以理解，图8仅仅是终端8的示例，并不构成对终端8的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述终端还可以包括输入输出设备、网络接入设备、总线等。
149.所称处理器80可以是中央处理单元(central processing unit，cpu)，还可以是其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
150.所述存储器81可以是所述终端8的内部存储单元，例如终端8的硬盘或内存。所述存储器81也可以是所述终端8的外部存储设备，例如所述终端8上配备的插接式硬盘，智能存储卡(smart media card，smc)，安全数字(secure digital，sd)卡，闪存卡(flash card)等。进一步地，所述存储器81还可以既包括所述终端8的内部存储单元也包括外部存储设备。所述存储器81用于存储所述计算机程序以及所述终端所需的其他程序和数据。所述存
储器81还可以用于暂时地存储已经输出或者将要输出的数据。
151.所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本技术的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
152.在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
153.本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
154.在本发明所提供的实施例中，应该理解到，所揭露的装置/终端和方法，可以通过其它的方式实现。例如，以上所描述的装置/终端实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。
155.所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
156.另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
157.所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个基于改进小波包能量曲率的斜拉桥损伤识别方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、u盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖
区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括是电载波信号和电信信号。
158.以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。