《一种提高随机损伤定位向量法监测结果准确性的方法.pdf》由会员分享，可在线阅读，更多相关《一种提高随机损伤定位向量法监测结果准确性的方法.pdf（9页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103884776 A (43)申请公布日 2014.06.25 CN 103884776 A (21)申请号 201410123708.6 (22)申请日 2014.03.28 G01N 29/04(2006.01) G01N 29/44(2006.01) (71)申请人 大连理工大学 地址 116024 辽宁省大连市甘井子区凌工路 2 号 (72)发明人 安永辉 欧进萍 (74)专利代理机构 大连理工大学专利中心 21200 代理人 赵连明 (54) 发明名称 一种提高随机损伤定位向量法监测结果准确 性的方法 (57) 摘要 本发明公开了一种提高随机损伤定位向量法。
2、 监测结果准确性的方法， 涉及到 SDLV 法定位效果 的改进， 属于结构健康监测技术领域。其特征是， 在SDLV法的步骤中增加一个选择观测矩阵C的步 骤从而可获得更准确的定位结果 ； 具体如下 ： 从 SDLV法的核心入手分析不同测点数比例时C矩阵 的选择与损伤定位结果之间的相关性 ； 进而确定 了适用于常规和稀疏测点布置方案下 SDLV 法的 C 矩阵选择依据， 即测点数百分比 ； 最后提出了两 种确定测点数百分比阈值的方法。本发明提高了 基于 SDLV 法的桁架结构损伤监测结果的准确性， 有利于桁架结构的实时安全监测在工程中的应 用。 (51)Int.Cl. 权利要求书 2 页 说明书 。
3、5 页 附图 1 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书5页 附图1页 (10)申请公布号 CN 103884776 A CN 103884776 A 1/2 页 2 1. 一种提高随机损伤定位向量法监测结果准确性的方法， 由被测试桁架、 加速度传感 器和计算机组成 ； 将双向加速度传感器安装在桁架结构的待监测子结构桁架节点上， 计算 机中的数据采集模块采集所有测点的加速度振动响应， 计算机中的数据分析模块进行数据 分析， 最终由计算机的损伤监测结果显示模块给出桁架结构的损伤位置 ； 具体随机损伤定 位向量法包括如下步骤 ： 步骤一、 数据采集模。
4、块采集桁架结构节点测点的纵向和竖向振动响应， 其中纵向为沿 着桁架结构跨径的方向 ； 步骤二、 数据分析模块中基于 NExT-ERA 或其他参数识别技术识别出前几阶振型 步骤三、 数据分析模块中建立基于模态参数的两种 C 矩阵计算模块 ： 第一种 C 矩阵计算模块 第二种 C 矩阵计算模块 其中是振型的共轭复数 ； 步骤四、 当系统输入全部未知时， 数据分析模块中构建如下 Q 矩阵 其中 A 是系统矩阵、 C 是输出矩阵， 均由测点处的振动数据识别得到 ； p 是常数， 其值为 0、 1、 2 时分别对应当输出数据为位移、 速度、 加速度时， 代表伪逆 ； 步骤五、 数据分析模块中确定 SDL。
5、V 向量的个数 ： 由于实测数据存在噪声并且用来计算 的模态数有限， 导致奇异值分解的 S 矩阵中出现了非零项 X， 即 SDLV 向量的个数 q 为满足上式的个数的一半， 向下取整 ； 步骤六、 数据分析模块中联合多 SDLV 向量进行损伤定位 ： 获得损伤定位向量后， 将其 作为荷载施加到无损结构有限元模型上， 求得各单元的应力即为单元特征应力 ， 实际应 用中将 WSI 符合下式的单元确定为损伤单元 其中 WSI 为权重应力指数， j为单元的特征应力， R 建议取 1， b 建议取不大于 0.1 的 数， 一般取 0.1 ； 其特征在于， 在步骤四之前 , 在数据分析模块中增加计算测点数。
6、百分比阈值 oa 的步 权 利 要 求 书 CN 103884776 A 2 2/2 页 3 骤， 并根据计算的测点数百分比阈值 oa 的大小从两个 C 矩阵中选择最优的 C 矩阵计算模 块， 从而使损伤定位结果更加准确 ： 当待测子结构的测点数百分比小于等于测点数百分比 阈值 oa 时 ( 即满足式 (7)， 选择第二种 C 矩阵计算模块 ； 否则， 待测子结构的测点数百分 比大于测点数百分比阈值 oa 时选择第一种 C 矩阵计算模块 ； nmeasured/Nall oa (7) 其中 nmeasured/Nall是测点数百分比， nmeasured是测点数， Nall是待测结构平面里所有。
7、节点 数， oa 是测点数百分比阈值。 2. 根据权利要求 1 所述的一种提高随机损伤定位向量法监测结果准确性的方法， 其特 征在于， 按如下方式确定测点数百分比阈值 oa ： 首先定义一个衡量损伤定位结果准确性的 标准 ； 然后在有限元模型中针对已知的某个损伤工况进行数值模拟并观察规律 ； 最终确定 测点数百分比阈值 oa， 工程应用中根据该测点数百分比阈值 oa 确定待监测子结构应选取 的 C 矩阵 ； 在数据分析模块中定义衡量损伤定位结果准确性的标准， 包括如下内容 ： 比较基于两 个 C 矩阵的损伤定位结果时， 符合下列条件的结果更准确 ： 1) 结果应成功检测到损伤单元 ； 2)结果。
8、应有较少的误报单元 ； 3)若前两个条件都已满足， 令threshold取值从0.1到1以0.1 的间隔取值， 然后定义两个数 Pu和 Pd， 并且 P=Pu+Pd： Pu是损伤定位结果中无损单元中应力 小于每一个阈值threshold的总个数 ； 当基于每个threshold进行判别时， d是损伤单元中应力 比另一个结果大的单元的个数， Pd=10d ； 当 P 值较小时， 认为该结果更准确 ； 在待测结构的有限元模型中针对已知的损伤工况进行数值模拟 ： 测点数逐渐增多， 选 择不同大小的子结构分别将不同的 C 矩阵模块代入 SDLV 法进行 100 次模拟计算， 针对不同 大小的子结构统计。
9、采用哪一个 C 矩阵模块时其损伤定位结果准确率更高 ； 即确定使用第一 个 C 矩阵可在多大的子结构中得到更准确的结果， 使用第二个 C 矩阵可以在多大的子结构 中得到更准确的结果， 最终确定测点数百分比阈值 oa。 3. 根据权利要求 1 所述的一种提高随机损伤定位向量法监测结果准确性的方法， 其特 征在于， 根据以下分析从损伤前后的柔度矩阵差 Q 矩阵的维数确定阈值 oa 的值 ； 当采用 第一种 C 矩阵计算模块计算矩阵 Q 时， 对 Q 进行奇异值分解后 V 矩阵的属性随着矩阵 Q 的维数改变而改变 ； 1)当方阵Q的维数小于某个维数阈值时， 奇异值分解之后， V矩阵虚部的所有列都为 。
10、0， 所以 V 是一个正交矩阵 ； 此时采用第二种 C 矩阵计算模块 ； 2) 然而， 随着矩阵 Q 维数变大， V 矩阵的最后几列元素的虚部不再为 0 ； 并且虚部不 为 0 的列数随着 Q 维数的增大而增加 ； 当这个维数增大到某个维数阈值时， V 的实部不再 是正交矩阵， 此时采用第一种 C 矩阵计算模块 ； 由于矩阵 Q 的维数与测点数相关， 根据这里提到的 Q 的维数阈值计算出测点百分 比阈值 oa。 4. 根据权利要求 1 所述的一种提高随机损伤定位向量法监测结果准确性的方法， 其特 征在于， 该方法用于桁架结构的常规测点布置方案或用于桁架结构的稀疏测点布置方案。 权 利 要 求 。
11、书 CN 103884776 A 3 1/5 页 4 一种提高随机损伤定位向量法监测结果准确性的方法 技术领域 0001 本发明属于结构健康监测技术领域， 涉及到一种桁架结构实时安全监测的重要方 法损伤定位向量 (SDLV) 法的进一步研究， 具体地说是从 SDLV 法的核心入手进一步研究从 而提出了一种可以提高 SDLV 法损伤定位结果准确性的方法。 背景技术 0002 结构健康监测技术是土木工程及机械和航空航天等领域的研究热点， 而损伤识别 是结构健康监测的核心技术 ； 基于振动信息的结构损伤识别通过对振动响应信号进行分析 去定位损伤， 一些损伤识别方法可实现在线实时监测。桁架结构 ( 网。
12、架结构 ) 是广泛用在 各类工业结构和工程结构中的常见结构类型， 但许多桁架结构如钢桁桥结构常在极端环境 中工作， 损伤常常发生并且难以避免。很多损伤识别方法都以各种结构动力指纹为基础构 建损伤指标， 其中结构模态频率和振型最基本也最常用。 0003 本发明的直接背景技术是随机损伤定位向量 (Stochastic Damage Locating Vectors) 法 ( 简称 SDLV 法 )， 该方法由美国东北大学的 Bernal 提出并验证， 较好地解决了 Bernal 之前提出的损伤定位向量 (DLV) 法在环境激励等未知输入下的应用问题， 为工程应 用提供了重要支持 ； Nagayam。
13、a 等、 Gao 等也对该法做了很多研究。SDLV 法主要用于桁架结 构、 网架结构等的损伤定位， 其损伤指标也是基于模态信息构建而来， 该法的主要过程是 ： 通过测得的振动响应求出结构 / 子结构的柔度矩阵， 根据损伤前后柔度矩阵的变化求出损 伤定位向量(一个荷载向量集合)并将其视为静力作用在无损结构的有限元模型上求出每 个单元的应力， 应力为零的单元即为损伤单元(理论上为零， 实际应用中应力小于0.1的单 元为损伤单元 )。 0004 基于 SDLV 法定位损伤有以下 6 个流程步骤 ： 0005 步骤一、 数据采集模块采集桁架结构节点测点的纵向和竖向振动响应 ( 如加速 度、 速度、 位。
14、移 )， 纵向为沿着桁架结构跨径的方向 ； 0006 步骤二、 数据分析模块中基于 NExT(James 等提出 )-ERA(Juang 和 Pappa 提出 ) 或 其他参数识别技术识别出前几阶振型 0007 步骤三、 数据分析模块中建立基于模态参数的两种 C 矩阵 (Nagayama 等提出 ) 计 算模块 ： 0008 第一种 C 矩阵计算模块 0009 0010 第二种 C 矩阵计算模块 0011 0012 其中是振型的共轭复数 ； 说 明 书 CN 103884776 A 4 2/5 页 5 0013 步骤四、 当系统输入全部未知时， 数据分析模块中构建如下 Q 矩阵 0014 00。
15、15 其中 A 是系统矩阵、 C 是输出矩阵， 均由测点处的振动数据识别得到 ； p 是常数， 其 值为 0、 1、 2 时分别对应当输出数据为位移、 速度、 加速度时， 代表伪逆。 0016 步骤五、 数据分析模块中确定 SDLV 向量的个数 ： 由于实测数据存在噪声并且用来 计算的模态数有限， 导致奇异值分解的 S 矩阵中出现了非零项 X， 即 0017 0018 SDLV 向量的个数 q 为满足上式的个数的一半， 向下取整。 0019 0020 步骤六、 数据分析模块中联合多 SDLV 向量进行损伤定位 ： 获得损伤定位向量后， 将其作为荷载施加到无损结构有限元模型上， 求得各单元的应力。
16、即为单元特征应力 ， 实 际应用中将 WSI 符合下式的单元确定为损伤单元 0021 0022 其中 WSI 为权重应力指数， j为单元的特征应力， R 建议取 1， b 建议取不大于 0.1 的数， 一般取 0.1。 0023 存在的问题是 ： 已有研究给出了两种构建 C 矩阵的方法， 但并没有对其做详细讨 论、 也没有给出其选择方法 ； 这造成基于 SDLV 法的一些损伤定位结果准确性不高， 如未能 识别出损伤单元 ( 漏报单元 )、 或者存在多识别的现象 ( 误报单元 ) 等 ； 同时同样的振动响 应数据， 若使用另一个 C 矩阵时却能得到更准确的损伤定位结果。那么如何根据具体的监 测对。
17、象选择合适的 C 矩阵才能得到更准确的损伤定位结果是 SDLV 法亟待解决的一个重要 问题。 发明内容 0024 为了降低基于 SDLV 法的损伤定位结果中一些单元被误判或漏判的概率， 发明了 一种提高随机损伤定位向量法监测结果准确性的方法， 本发明要解决的核心技术问题是提 供一种提高 SDLV 法损伤定位结果精确性的方法， 通过进一步研究观测矩阵 C， 发现损伤定 位结果的精确性与 C 矩阵的选择直接相关 ； 同时发现 C 的选择规则与测点数百分比阈值 oa 相关， 而测点数百分比阈值oa是直接影响到C矩阵计算模块选择的参数 ； 为此， 本发明提出 了两种可以确定测点数百分比阈值 oa 的方。
18、法。 0025 第一种是通过对待监测桁架结构中不同的子结构进行数值模拟来确定。在上述 SDLV法步骤四之前,在数据分析模块中增加计算测点数百分比阈值oa的步骤， 并根据计算 说 明 书 CN 103884776 A 5 3/5 页 6 的测点数百分比阈值 oa 的大小从两个 C 矩阵中选择最优的 C 矩阵计算模块， 从而使损伤定 位结果更加准确 ： 当待测子结构的测点数百分比小于等于测点数百分比阈值oa时(即满足 式 (7)， 选择第二种 C 矩阵计算模块 ； 否则， 待测子结构的测点数百分比大于测点数百分比 阈值 oa 时选择第一种 C 矩阵计算模块 ； 0026 nmeasured/Nal。
19、l oa (7) 0027 其中 nmeasured/Nall是测点数百分比， nmeasured是测点数， Nall是待测结构平面里所有节 点数， oa 是测点数百分比阈值。 0028 测点数百分比阈值 oa 的确定方式如下 ： 首先定义一个衡量损伤定位结果准确性 的标准 ； 然后在有限元模型中针对已知的某个损伤工况进行数值模拟并观察规律 ； 最终确 定测点数百分比阈值 oa， 工程应用中根据该测点数百分比阈值 oa 确定待监测子结构应选 取的 C 矩阵 ； 0029 在数据分析模块中定义衡量损伤定位结果准确性的标准， 包括如下内容 ： 比较基 于两个 C 矩阵的损伤定位结果时， 符合下列条。
20、件的结果更准确 ： 1) 结果应成功检测到损伤 单元 ； 2) 结果应有较少的误报单元 ； 3) 若前两个条件都已满足， 令 threshold取值从 0.1 到 1 以 0.1 的间隔取值， 然后定义两个数 Pu和 Pd， 并且 P=Pu+Pd： Pu是损伤定位结果中无损单元 中应力小于每一个阈值threshold的总个数 ； 当基于每个threshold进行判别时， d是损伤单元 中应力比另一个结果大的单元的个数， Pd=10d ； 当 P 值较小时， 认为该结果更准确 ； 0030 在待测结构的有限元模型中针对已知的损伤工况进行数值模拟 ： 测点数逐渐增 多， 选择不同大小的子结构分别将。
21、不同的 C 矩阵模块代入 SDLV 法进行 100 次模拟计算， 针 对不同大小的子结构统计采用哪一个 C 矩阵模块时其损伤定位结果准确率更高 ； 即确定使 用第一个 C 矩阵可在多大的子结构中得到更准确的结果， 使用第二个 C 矩阵可以在多大的 子结构中得到更准确的结果， 最终确定测点数百分比阈值 oa。 0031 第二种是根据以下分析从损伤前后的柔度矩阵差 Q 矩阵的维数确定测点数百 分比阈值 oa 的值 ； 当采用第一种 C 矩阵计算模块计算矩阵 Q 时， 对 Q 进行奇异值分解 后 V 矩阵的属性随着矩阵 Q 的维数改变而改变 ； 0032 1)当方阵Q的维数小于某个维数阈值时， 奇异。
22、值分解之后， V矩阵虚部的所有列 都为 0， 所以 V 是一个正交矩阵 ； 此时采用第二种 C 矩阵计算模块 ； 0033 2) 然而， 随着矩阵 Q 维数变大， V 矩阵的最后几列元素的虚部不再为 0 ； 并且虚 部不为 0 的列数随着 Q 维数的增大而增加 ； 当这个维数增大到某个维数阈值时， V 的实部 不再是正交矩阵， 此时采用第一种 C 矩阵计算模块 ； 0034 由于矩阵 Q 的维数与测点数相关， 根据这里提到的 Q 的维数阈值可以计算出 测点百分比阈值 oa。 0035 该方法既可用于桁架结构的常规测点布置方案， 也可用于桁架结构的稀疏测点布 置方案。 0036 本发明指出损伤识。
23、别结果的准确性与两个 C 矩阵计算模块相关， 并给出了 C 矩阵 计算模块的选择规则， 按此规则选择 C 矩阵可以获得更准确的损伤定位结果。本发明的有 益效果是提供了一种提高桁架结构损伤定位结果精确性的系统和方法， 并且实施简单， 具 有较好的实用价值。 说 明 书 CN 103884776 A 6 4/5 页 7 附图说明 0037 图 1 在桁架结构的两个不同大小的子结构中监测损伤单元 24 的示意图。 0038 图 2 噪声水平 5% 下基于第一种 C 矩阵计算模块的子结构 1 的损伤定位结果。 0039 图 3 噪声水平 5% 下基于第二种 C 矩阵计算模块的子结构 1 的损伤定位结果。
24、。 0040 图 4 噪声水平 5% 下基于第一种 C 矩阵计算模块的子结构 2 的损伤定位结果。 0041 图 5 噪声水平 5% 下基于第二种 C 矩阵计算模块的子结构 1 的损伤定位结果。 0042 图 2- 图 5 中的横坐标是桁架杆编号， 纵坐标是归一化累计特征应力， 若某杆的归 一化累计特征应力小于 0.1， 则其为损伤杆。 0043 图中 ： 子结构 1 是包含 7 个测点的实线矩形框中的子结构 ； 子结构 2 是包含 17 个 测点的虚线矩形框中的子结构。 具体实施方式 0044 下面结合附图和具体例子对本发明的技术方案进行详细说明 ： 0045 (1) 建立结构有限元模型 ：。
25、 以我国常见的在役贝雷梁式钢桁桥为例， 该模型由 108 个桁架节点和 312 根桁架杆单元连接而成， 下称该结构为桁架结构 A ； 如图 1 是该桁架结构 的一个竖向平面的示意图。 0046 (2) 针对具体结构确定其测点数百分比阈值 oa ： 从小 ( 三节点的三角形子结构 ) 到大选取不同大小的子结构按照 SDLV 法损伤定位的 6 个流程步骤 ( 见背景技术 ) 进行损 伤定位研究。 0047 测点数百分比阈值 oa 的第一种确定方法 ： 针对不同大小的待测子结构分别选择 不同的 C 矩阵模块分别代入 SDLV 法进行 100 次模拟计算并统计采用哪一个模块的结果准 确率更高 ； 即确。
26、定使用第一个 C 矩阵可在多大的子结构中得到更准确的结果， 使用第二个 C 矩阵可以在多大的子结构中得到更准确的结果。常规测点布置方案下 ( 即待测子结构所有 节点均为测点 ) 桁架结构 A 中不同大小的待测子结构 ( 即测点数百分比不同 ) 的损伤定位 结果呈现出的规律如表 1 所示 ； 表 1 中的规律适用于桁架结构 A 中任意的子结构。从表 1 中确定测点数百分比阈值 oa 为 25.00%。 0048 表 1 桁架结构 A 中不同大小子结构的损伤定位结果 0049 0050 测点数百分比阈值 oa 的第二种确定方法 ： 根据以下分析从损伤前后的柔度矩阵 差 Q 矩阵的维数来确定测点数百。
27、分比阈值 oa 的值。当基于第一种 C 矩阵计算模块计算 矩阵 Q 时， 对 Q 进行奇异值分解后 V 矩阵的属性随着矩阵 Q 的维数改变而改变。当 使用第一种C矩阵计算模块时， V矩阵有一个与复数矩阵Q的维数相关的属性 ： 1)当方阵 Q 的维数小于某个维数阈值时， 奇异值分解之后， V 矩阵虚部的所有列都为 0， 所以 V 是一 个正交矩阵。2) 随着矩阵 Q 维数变大， V 矩阵的最后几列元素的虚部不再为 0 ； 当这个维 说 明 书 CN 103884776 A 7 5/5 页 8 数增大到某个维数阈值时， V 的实部不再是正交矩阵。该维数阈值与测点数相关， 为 2 倍测 点数 *2 。
28、倍测点数， 进而可计算出测点百分比阈值 oa。 0051 (3) 观测矩阵 C 的选择规则 0052 当满足式 (7) 时， 第二种 C 矩阵计算模块能给出更精确的结果， 否则第一种 C 矩阵 计算模块的损伤定位结果更准确。 0053 nmeasured/Nall oa (7) 0054 其中 nmeasured是测点数， Nall是待测结构平面里所有的节点数， oa 是测点数百分比 阈值。 0055 (4) 损伤定位结果准确性的提高效果 0056 图1中实线、 虚线矩形框中分别为两个不同大小的子结构， 当单元24损伤时， 噪声 水平 5% 下采用两种 C 矩阵计算模块的损伤定位结果如图 2 。
29、5 所示 ： 当在子结构 1 中监测 损伤单元 24 时， 其测点百分比 7/5225.00%， 所以应使用第一种 C 矩阵计 算模块结果较好， 从图45中的结果也能看出图4中的结果准确识别出了损伤单元24(比 图 5 更准确 )。 0057 本发明不受桁架结构形式的约束， 同样适用于其他形式的桁架结构。 例如， 在某桁 架结构基于SDLV法的损伤定位中， 当研究某子结构时采用不同的振动数据的100次数值模 拟中， 使用第一种 C 矩阵计算模块的损伤定位结果约 90 次都比第二种 C 矩阵计算模块更准 确。这个结果说明了 C 矩阵的选择很重要， 同时说明所发明的损伤监测系统对基于 SDLV 法 的损伤定位结果准确性的提高有重要的积极作用。 说 明 书 CN 103884776 A 8 1/1 页 9 图 1 图 2图 3 图 4 图 5 说 明 书 附 图 CN 103884776 A 9 。