CN103793579B

CN103793579B - 一种基于ei和逐步累积法的传感器优化布设方法

Info

Publication number: CN103793579B
Application number: CN201410072194.6A
Authority: CN
Inventors: 梁鹏; 贺敏; 王晓光; 杨国强; 鱼江英; 董向前; 孟修建; 桂晓华; 张立龙; 吴向男; 马旭明; 李斌
Original assignee: Changan University
Current assignee: Changan University
Priority date: 2014-02-28
Filing date: 2014-02-28
Publication date: 2017-08-04
Anticipated expiration: 2034-02-28
Also published as: CN103793579A

Abstract

本发明公开了一种基于EI和逐步累积法的传感器优化布设方法，该方法提供的公式可以直接用于计算不同结构类型、不同计算跨径桥梁的传感器优化布设位置。适用于桥梁结构计算分析技术领域。利用本发明的方法，MAC最大非对角元仅为EI法的五十分之一，逐步累积法的十分之一，且运行时间相对较短。因此，本发明解决了现有优化算法的不足。

Description

一种基于EI和逐步累积法的传感器优化布设方法

技术领域

本发明属于桥梁结构计算分析技术领域，涉及一种传感器优化布设方法，尤其是一种基于EI和逐步累积法的传感器优化布设方法。

背景技术

随着桥梁结构的跨径和规模越来越大，对桥梁结构在运营状态下进行健康监测是一个值得研究的热点问题。在桥梁健康监测系统中，传感器是整个系统的“眼睛”。通过传感器系统，桥梁健康监测系统可以实时地对结构进行监测，并传回数据，为后续工作打下基础。传感器优化布设最为设计结构健康监测系统的关键问题之一，是实现桥梁健康监测系统的第一步，对数据采集有着至关重要的影响。而传感器系统本身是一个“矛盾体”。一方面，为了对结构的健康状况做出合理的、科学的评价，需要使用尽量多的传感器，监测尽量多的醒目。另一方面，为了保证监测系统本身的可靠性和经济性，要求使用的监测设备不能过多。在已实施的桥梁健康监测系统中，有些系统规模过于庞大，系统本身的可靠性都无法保证；而有些结构的监测系统则想反，规模过于精简，以至于无法对结构的健康状态做出科学的评估。因此在应用中需要找到合适的方法寻求两者最佳的平衡点。

目前，传感器优化布置主流算法主要有2种：有效独立法（EI）和序列法（MAC）。有效独立法不计及结构质量分布的不均匀性，可能会出现选取测点能量低导致信息丢失的情况。而序列法在待选测点较多时，计算量较大。

发明内容

本发明的目的在于克服上述现有技术的缺点，提供一种基于EI和逐步累积法的传感器优化布设方法，其可以计算不同结构类型、不同计算跨径桥梁。

本发明的目的是通过以下技术方案来实现的：

这种基于EI和逐步累积法的传感器优化布设方法，包括以下步骤：

1）通过有限元计算，得到结构的模态向量矩阵，确定监测振型的数目，确定传感器的布置数目；

2）利用EI法对模态向量矩阵进行优化，选取初始测点，形成初始模态向量阵；并得到剩余测点矩阵,形成候选添加测点集合；

3）求解初始模态向量阵的置信度矩阵MAC，并求MAC矩阵的最大非对角元max；

4）从步骤2）确定的候选添加测点集合中选择一个点添加到初始测点中，将第k个自由度添加到测量自由度中，求其模态向量阵的〖(MAC]_tj)_k矩阵，并计算其最大非对角元max_k；其中k为自然数；

5）计算f（k）=max_k-max的值；

6）对所有候选测点重复步骤3）至步骤5），得到数组f（i），（i=1，…，m），并将|f（i）|（f（i）<0）最小值对应的测点加入到已有测点中；其中m为自然数；

7）以添加了测点后的测点集合作为初始测点，再重复步骤4）～步骤6），直到测点数目与传感器布置数目相同。

进一步的，利用EI法对模态向量矩阵进行优化，产生初始测点集合，然后通过逐步累积法选取测点，具体包括以下步骤：

1）利用设计图纸，建立桥梁有限元模型并进行分析计算，得到其模态向量矩阵Φ_s；

2）利用Φ_s计算得E矩阵；

将E矩阵的对角元素写成一个列向量：

E_D={E₁₁,E₂₂,…，E_ss}^T （2）

E是幂等矩阵，其对角线上第i个元素表示第i个自由度或测试点对矩阵Φ_s秩的贡献，也就是对矩阵A₀的贡献；如果对角元素E_ii=0，则表示在相应第i个测点上无法识别所关心的模态，如果对角元素E_ii=1，则表示相应测点是关键点，不能排除；因此，E代表候选传感器位置几何的有效独立分布，E对角线上的元素代表响应传感器候选测点对模态矩阵的线性无关的贡献；

3）获得矩阵E后，通过对角元素的大小来对各个候选点的优先顺序进行排序，用迭代算法每次排除相应对角元素最小的测点，再进行下一次迭代，直至达到初始测点数目，并得到与初始测点位置对应的模态矩阵；

4）将EI法选取测点的模态矩阵Φ（Φ∈R^n×m，一般有m<n，并且r(Φ)=m，即矩阵Φ列满秩）作为初始模态矩阵，将结构的剩余自由度组成的模态矩阵作为候选模态矩阵；记由Φ得到的MAC矩阵的最大非对角元为MAX=MAC_ij，增加的第k个自由度后MAC_ij变为(MAC_ij)_k，则有：

显然有：A>0，B>0，C>0；所以，要使MAC_ij降低，只需AD-BC<0；

5）分别增加剩余自由度中的每个自由度，计算增加后的MAC最大非对角元并计算增加前后MAC最大非对角元的差值；即计算过程4）中增加的第k个自由度前后MAC最大非对角元变化的差值，即f（k）＝max_k-max的值；

6）对所有的测点重复步骤4）至步骤5），得到数组f（i），（i=1，…，m），并将|f（i）|（f（i）<0）最小值对应的测点加入到已有测点中，形成新的初始测点集和新的初始模态向量矩阵；其中m为自然数；

本发明具有以下有益效果：

本发明得到的是次优解，由振动理论可知，如果待识别的模态个数为m，则m个传感器就可以满足参数识别的要求，其他的传感器事实上是出于振型可视化或者振型匹配的考虑，所以次优解是可以接受的。并且，利用本发明的方法，MAC最大非对角元仅为EI法的五十分之一，逐步累积法的十分之一，且运行时间相对较短。因此，本发明解决了现有优化算法的不足。

附图说明

图1为EI-逐步累积法的操作流程图；

图2为有限元模型主梁示意图；

图3为前4阶振型测点MAC值示意图；

图4为EI-逐步累积法位置示意图；

图5（a）为有效独立法（EI）位置示意图；

图5（b）为逐步累积法位置示意图；

图5（c）为EI-逐步累积法位置示意图。

具体实施方式

参见图1，本发明基于EI和逐步累积法的传感器优化布设方法，包括以下步骤：

2）利用EI法对模态向量矩阵进行优化，选取初始测点，形成初始模态向量阵；并得到剩余测点矩阵，形成候选添加测点集合；

5）计算f（k）＝max_k-max的值；

以上利用EI法对模态向量矩阵进行优化，产生初始测点集合，然后通过逐步累积法选取测点，具体包括以下步骤：

2）利用Φ_s计算得E矩阵；

将E矩阵的对角元素写成一个列向量：

E_D={E₁₁,E₂₂,…，E_ss}^T （2）

显然有：A>0，B>0，C>0；所以，要使MAC_ij降低，只需AD-BC<0；

下面结合附图对本发明做进一步详细描述：

如图2，以苏通长江大桥辅航道桥（140+268+140）m的连续刚构桥为例进行传感器优化布置。

1，建立全桥Midas/Civil有限元模型主梁158个单元，159个节点。如图2。

2，进行模态分析，提取前4阶竖向振型Ф。

3，前4阶振型测点MAC值见图3。

从图3可以发现，测点MAC矩阵中，最大非对角元值为0.0048676，相关振型为第1阶和第3阶，说明前4阶竖向振型具有良好的正交性。

4,确定传感器数量。

传感器的安装数量需要人为指定。设传感器数量为9个。

5，进行优化布置计算。

（1）对有限元得到的模态矩阵Ф。利用EI法选取初始测点，初始测点个数为模态监测阶数，同时得到剩余测点矩阵为候选测点集。

（2）计算初始测点模态矩阵的MAC值，并求MAC的最大非对角元max；

（3）从第（1）步确定的候选添加测点集合中选择一个点添加到初始测点中，例如第k个自由度添加到测量自由度中，求其模态向量的（MAC_ij）_k矩阵，并计算其最大非对角元max_k；

（4）计算f（k）＝max_k-max的值；

（5）对所有候选测点重复第3至第5步，得到数组f（i），（i=1，…，m），并将|f（i）|（f（i）<0）最小值对应的测点加入到已有测点中；

（6）以添加了测点后的测点集合作为初始测点，再重复（3）～（5），直到得到满意的测点数目，本例中满意测点数目为9个。

6，得出优化结果。

优化方案结果如下：

传感器安装节点为：15、16、70、73、86、88、90、144、145

传感器安装位置示意图：见图4

该安装方案下MAC最大非对角元：0.02189

计算过程所需时间为：1.7820s

7，结果对比。

利用现有方法对该桥进行传感器优化布设计算，并与本发明进行对比，结果如下：

各方法位置示意如图5。

由上表和图可知，EI法效果很差，即MAC最大非对角元很大，接近于1。且所得方案中5个测点相邻，造成了测量振型可视化程度下降。逐步累积法效果一般，但运行时间较长。而本发明效果最好，MAC最大非对角元仅为EI法的五十分之一，逐步累积法的十分之一，且运行时间相对较短。解决了现有优化算法的不足。

Claims

1.一种基于EI和逐步累积法的传感器优化布设方法，其特征在于，包括以下步骤：

1)通过有限元计算，得到结构的模态向量矩阵，确定监测振型的数目，确定传感器的布置数目，其中，传感器的安装数量需要人为指定；

2)利用EI法对模态向量矩阵进行优化，选取初始测点，形成初始模态向量阵，初始测点个数为模态监测阶数；并得到剩余测点矩阵,形成候选添加测点集合；

3)求解初始模态向量阵的置信度矩阵MAC，并求MAC矩阵的最大非对角元max；

4)从步骤2)确定的候选添加测点集合中选择一个点添加到初始测点中，将第k个自由度添加到测量自由度中，求其模态向量阵的(MAC_ij)_k矩阵，并计算其最大非对角元max_k；其中k为自然数；

5)计算f(k)＝max_k-max的值；

6)对所有候选测点重复步骤3)至步骤5)，得到数组f(i)，(i＝1，…，m)，并将|f(i)|(f(i)<0)最小值对应的测点加入到已有测点中；其中m为自然数；

7)以添加了测点后的测点集合作为初始测点，再重复步骤4)～步骤6)，直到测点数目与传感器布置数目相同；

利用EI法对模态向量矩阵进行优化，产生初始测点集合，然后通过逐步累积法选取测点，具体包括以下步骤：

a)利用设计图纸，建立桥梁有限元模型并进行分析计算，得到其模态向量矩阵Φ_s；

b)利用Φ_s计算得E矩阵；

将E矩阵的对角元素写成一个列向量：

E_D＝{E₁₁,E₂₂,…，E_ss}^T (2)

E是幂等矩阵，其对角线上第i个元素表示第i个自由度或测试点对矩阵Φ_s秩的贡献，也就是对矩阵A₀的贡献；如果对角元素E_ii＝0，则表示在相应第i个测点上无法识别所关心的模态，如果对角元素E_ii＝1，则表示相应测点是关键点，不能排除；因此，E代表候选传感器位置几何的有效独立分布，E对角线上的元素代表响应传感器候选测点对模态矩阵的线性无关的贡献；

c)获得矩阵E后，通过对角元素的大小来对各个候选点的优先顺序进行排序，用迭代算法每次排除相应对角元素最小的测点，再进行下一次迭代，直至达到初始测点数目，并得到与初始测点位置对应的模态矩阵；

d)将EI法选取测点的模态矩阵Φ(Φ∈R^n×m，一般有m<n，并且r(Φ)＝m，即矩阵Φ列满秩)作为初始模态矩阵，将结构的剩余自由度组成的模态矩阵作为候选模态矩阵；记由Φ得到的MAC矩阵的最大非对角元为MAX＝MAC_ij，增加的第k个自由度后MAC_ij变为(MAC_ij)_k，则有：

显然有：A>0，B>0，C>0；所以，要使MAC_ij降低，只需AD-BC<0；

e)分别增加剩余自由度中的每个自由度，计算增加后的MAC最大非对角元并计算增加前后MAC最大非对角元的差值；即计算过程d)中增加的第k个自由度前后MAC最大非对角元变化的差值，即f(k)＝max_k-max的值；

f)对所有的测点重复步骤d)至步骤e)，得到数组f(i)，(i＝1，…，m)，并将|f(i)|(f(i)<0)最小值对应的测点加入到已有测点中，形成新的初始测点集和新的初始模态向量矩阵；其中m为自然数；

g)以添加了测点后的测点集合作为初始测点，再重复步骤d)～步骤f)，直到测点数目与传感器布置数目相同。