
Citation: | LI Yongle, PAN Junzhi, TI Zilong, RAO Gang. Inversion Method of Vortex-Induced Vibration Amplitude for Long-Span Bridges with Partially Installed Noise Barrier[J]. Journal of Southwest Jiaotong University, 2023, 58(1): 183-190. doi: 10.3969/j.issn.0258-2724.20210172 |
The sectional model test in wind tunnels is often used to measure the vortex-induced vibration (VIV) of long-span bridges. Since the sectional model test is based on two-dimensional theory, when the bridge has different aerodynamic configurations along the span due to the partial installation of noise barriers, it is difficult to measure the VIV response directly through the sectional model test. Based on the empirical linear VIV model, an assessment method of VIV between the sectional model and prototype bridge that considers the effects of multiple aerodynamic configurations is proposed. Firstly, the sectional model test is performed on the models with and without barriers respectively. Then, the prototype response of the noise barriers installed and not installed along the span is investigated by ANSYS harmonic analysis, and the corresponding amplitude of the vortex-induced force is obtained. Finally, according to the actual installation position of the noise barrier along the span, the vortex-induced force is imposed on the bridge and the prototype response with the partially installed noise barrier is obtained. In addition, based on the method in this paper, various noise barrier installation schemes are numerically simulated. The results indicate that fully enclosed noise barrier will significantly reduce the aerodynamic performance of the main girder and the overall VIV will be affected by partial installation of barrier to a large degree. The method in this paper can estimate the prototype response of multi-aerodynamic configurations bridges through the results of sectional model tests. The installation of the noise barrier should be arranged on the side span as far as possible under the conditions of noise reduction. If the arrangement length exceeds the position of the bridge tower, it should be shortened as much as possible to reduce the vortex-induced response.
随着中速磁悬浮列车的广泛商业化应用[1],作为其关键核心技术的电磁悬浮控制技术已经成熟稳定[2]. 商业化运营过程中随着规模和复杂性的迅速增加,对中速磁悬浮列车悬浮系统的安全性和可靠性的需求不断增长. 在过去几十年中,对磁浮列车悬浮系统的性能监测尤其是异常检测的研究引起了极大的兴趣[3-4].
磁浮列车悬浮系统的异常指在列车实际运行中,悬浮系统出现工作状态与期望状态不匹配,但没有导致故障或失效的情况. 悬浮系统的异常具有出现时刻未知,持续时间短,通过控制器调节可以恢复稳定的特点. 对悬浮系统异常的精确检测可以提前获取系统的运行情况,有利于提前安排预防性维护计划,降低工程事故率[5].
针对磁悬浮列车的异常检测研究通常可以分为针对轨道的检测和针对悬浮系统的检测. 针对磁浮列车轨道的异常检测问题,Zhang等[6]提出基于小波与分形相结合的高速磁悬浮列车长轨道磁异常检测方法,以发现高速磁悬浮列车长定子铁芯的局部短路故障;Deng等[7]对高温超导磁悬浮列车运行时永磁导轨的不规则性进行测量和表征,并设计了一种不规则性检测设备;罗茹丹[8]提出了“2+1”的检测方法实现对长定子行波主漏磁场异常情况的动态检测;林国斌等[9]提出一种基于车轨状态监测的悬浮冗余控制系统,用于在控制列车稳定悬浮的同时对车轨状态进行检测,并通过对悬浮间隙和电磁铁振动情况的分类和学习确定状态类型,具有实时性好、容错性高等优点;杨杰等[10]提出一种基于悬挂式磁悬浮轨道交通系统的轨道维护设备,在对轨道状态进行取样和检测后对轨道异常实现诊断. 文献[11]收集了振动条件下高温超导钉扎系统的水平和垂直加速度数据.
针对磁浮列车悬浮系统的异常检测问题,朱跃欧等[12]提出一种磁浮列车悬浮控制器的异常预警方法与检测系统,通过分析悬浮间隙值所在区间,实现了车辆异常状况的在线检测与预警;王平等[13]提出了一种基于改进典型相关分析的中低速悬浮系统异常检测方法;通过运营线数据验证了该方法能获得更高的检测率;王平等[14]提出了一种基于超球体高斯分布的悬浮系统异常检测方法,该方法对系统异常检测的结果较为精确. 以新型悬挂式永磁磁浮列车研究对象,马政[15]设计并实现了一套基于物联网的永磁磁浮列车远程监控系统,实现了对列车悬浮状态的异常监控;Wang等[16]引入了支持向量机方法,利用悬浮传感器采集的实时信号来检测磁悬浮列车悬浮系统的异常状态;Ma等[17]将磁悬浮列车悬浮控制系统收集的数据用于训练机器学习模型,并根据数据确定磁悬浮列车的实际运行状态,以达到实时监控的目的.
然而,上述方法存在2个问题:1) 大部分方法是基于单变量的异常检测,没有充分利用系统的运行数据. 磁浮列车悬浮系统的异常检测大多都是只基于悬浮间隙,而实际运行中电流、加速度数据同样包含系统的异常信息. 因此,如何应用列车悬浮系统运行时产生的多变量数据进行异常检测是一个关键问题. 2) 需要有异常的先验信息,而实际系统的异常先验信息是不完全可知的,这无疑会降低异常检测的准确率;此外,面向实际工程系统的异常检测模块或系统在设计时需考虑异常漏检率和异常误检率的折中,这需要从数理的角度进行分析.
基于上述分析,本文将基于参数化残差的异常检测方法引入悬浮系统的异常检测问题,以最低异常漏检率为设计目标,提出了磁悬浮列车悬浮系统基于最低漏检率的残差生成器设计方法及异常检测方法. 结合磁悬浮列车悬浮系统的运行数据,验证了该残差生成器能够在满足期望误报率的情况下,实现异常漏检率的最小化.
系统出现异常时的状态空间模型[18]为
{x(k+1)=Ax(k)+Bu(k)+Bff(k)+w(k),y(k)=Cx(k)+Du(k)+Dff(k)+v(k), | (1) |
式中:x(k)为时刻k的系统状态;u(k)为时刻k的系统输入;y(k)为时刻k的系统输出;A、B、C、D为一定维数的系统参数矩阵;Bf和Df为异常参数矩阵;f(k)为异常向量;w(k)为系统的过程噪声;v(k)为系统的输出噪声.
则状态检测量z(k)为
{z(k)=Ψs⊥[us(k)ys(k)],us(k)=[uT(k−s)uT(k−s+1)⋯uT(k)]T,ys(k)=[yT(k−s)yT(k−sH)⋯yT(k)]T, | (2) |
式中:s为正整数,表示截断数据的长度;us(k)、ys(k)为长度为s+1的输入、输出数据矩阵;Ψs为基于数据驱动的系统信息矩阵[19];Ψs⊥为Ψs的正交补矩阵.
当系统出现异常时的状态检测量为
z(k)=Ψs⊥[us(k)ys(k)]=[0Hw,sWk,s+Vk,s+Hf,sFs], | (3) |
式中:Hw,s为过程噪声的Hankel矩阵;Wk,s为数据驱动的过程噪声矩阵;Vk,s为数据驱动的过程噪声矩阵;Hf,s为异常向量f(k)的Hankel矩阵;Fs为基于数据驱动的异常矩阵.
从式(3)可以看出,当系统出现异常时,z(k)不仅包含噪声信息,同时也包含了异常信息[20].
在基于模型的异常检测或故障检测领域,残差是通过比较实际系统的输出y(k)与估计输出ˆy(k)得到的,即定义r(k)=y(k)−ˆy(k). 残差数据中包含了系统扰动和系统异常的全部信息,残差产生器设计的主要目标是设法让异常检测系统对异常敏感,同时对其他信号具有鲁棒性.
结合Youla参数化和基于输入输出数据的状态检测量,残差生成器参数化形式为
r(k)=gTz(k)=gTΨs⊥[us(k)ys(k)], | (4) |
式中:g为参数向量,是残差生成器和异常检测的核心[21].
出于实际工程运营需要,通常在系统运行前就需要完成g与Ψs⊥的设计,并且在一段时间内都不能改变其参数.
由于残差数据中包含了系统扰动和系统异常的全部信息,因此,在系统运行时,可以通过对残差进行分析、处理,以及基于某种准则或阈值,来进一步掌握系统的运行状态信息,判断系统的运行状态.
当系统运行时,残差生成器实时获得残差数据r(k),并将残差数据送入残差评估单元,计算残差评估函数J(r)=Jk(r). 在设置评估函数阈值Jth后,根据式(5)判断系统是否发生异常
{Jk(r)>Jth,异常,Jk(r)⩽ | (5) |
传统的异常检测方案通常采用基于假设检验的方法进行残差评估,即假设残差的分布知识是准确已知的,或者可以从历史数据中估计出来[22]. 然而,对于磁浮列车悬浮系统而言,使用历史过程数据获得高度可信的概率分布估计比较繁琐,经验估计与真实概率分布的偏差可能会导致异常检测结果不可靠[23].
尽管难以获得悬浮系统的高度可信概率分布,但在无异常时的{{\boldsymbol{z}}}(k) 可以视为属于同一个集合. 考虑到悬浮系统运行数据的统计特性,可以依据数据的均值和协方差矩阵对其进行建模.
设无异常情况下,{{{\boldsymbol{z}}}}(k)满足
{S_{\rm{h}}} = \left\{ {{S_{\rm{z}}} \in {\zeta _{\rm{z}}}\left| {{{{\boldsymbol{E}}}} = {\bar {{\boldsymbol{z}}}_{\rm{h}}}{{,{\boldsymbol{V}}}} = {{\boldsymbol{\varSigma}} _{\rm{h}}}} \right.} \right\}, | (6) |
式中:{S_{\rm{h}}}为无异常数据的置信集;{\zeta _{{\rm{z}}}{\rm{}}}为所有{\boldsymbol{z}}(k)的分布集合;{\bar {\boldsymbol{z}}_{\rm{h}}}和{{\boldsymbol{\varSigma}} _{\rm{h}}}为无异常时状态检测变量的均值和方差,{\bar {\boldsymbol{z}}_{\rm{h}}}和{{\boldsymbol{\varSigma}} _{\rm{h}}}是通过对健康运行数据进行计算得来的,通常可以选择刚投入运行时无异常的运行数据,也可以是某一次无异常运行时的数据;E为{{{\boldsymbol{z}}}}(k)的均值;V为{{{\boldsymbol{z}}}}(k)的方差.
在悬浮系统运行初期,由于缺乏异常数据,难以分析出异常的数据分布特性. 因此,通常会依据正常运行时的数据特点预先设计一个异常检测系统,以实现对大部分潜在异常的检测. 假设系统运行时会出现M个典型异常,且发生第i个异常时的异常信号设为{{\boldsymbol{f}}_i}(k),i = 1,2, \cdots ,M. 为不失一般性,假设每个异常发生的概率相同. 采用均值和协方差矩阵对异常进行建模. 设异常情况时的 {{\boldsymbol{z}}}(k) 满足
{S _{{f_i}}} = \left\{ {{S_{\rm{z}}} \in {\zeta _{\rm{z}}}\left| {{{{\boldsymbol{E}}}}= {{\bar {{\boldsymbol{z}}}}_{{{\boldsymbol{f}}_i}}}{{,{\boldsymbol{V}}}} = {{\boldsymbol{\varSigma}} _{{\boldsymbol{f}}_i}}} \right.} \right\}, | (7) |
式中: {S _{{{\boldsymbol{f}}_i}}} 为异常数据的置信集;{\bar {\boldsymbol{z}}_{{{\boldsymbol{f}}_i}}}为{{{\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}}(k)的均值,{{{\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}}(k) 为第i个异常时的状态检测变量; {S _{{{\boldsymbol{f}}_i}}}为{{{\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}}(k)的方差.
当有可用的异常数据时,可以直接使用该数据建立置信集. 若无异常数据或异常数据少时,则可以自定义M个异常情况下的均值和协方差,进而构建置信集,当有新可用的异常数据时,记录数据并离线更新置信集.
针对第 i个异常,基于异常检测的残差信号为
{{{\boldsymbol{r}}}_i}(k) = {{\boldsymbol{g}}}_i^{\rm{T}}{{\boldsymbol{z}}}(k), | (8) |
式中:{{{\boldsymbol{g}}}_i}为第i个异常时的参数向量.
因此,第i个异常时的残差评估函数{\boldsymbol{J}}({\boldsymbol{r}}) 与阈值J_{{\rm{th}}}为
\left\{ \begin{gathered} {\boldsymbol{J}}({\boldsymbol{r}}) = {\left( {{{{\boldsymbol{r}}}_i}(k) - {{{\boldsymbol{r}}}_{i,{\rm{h}}}}} \right)^2}, \\ {{{\boldsymbol{r}}}_{i,{\rm{h}}}} = {{\boldsymbol{g}}}_i^{\rm{T}}{{\bar {\boldsymbol{z}}}_{\rm{h}}}, \\ {J_{{\rm{th}}}} = 1, \\ \end{gathered} \right. | (9) |
式中:{{{\boldsymbol{r}}}_{i,h}}为第i个无异常场景的残差.
定义异常误检率{P_{\rm{F}}} = \Pr \left\{ {J({\boldsymbol{r}}) > 1\left| {{{{\boldsymbol{f}}}_i}(k) = {\boldsymbol{0}}} \right.} \right\},定义异常漏检率为{P_{\rm{M}}} = \Pr \left\{ {J({\boldsymbol{r}}) \leqslant 1\left| {{{{\boldsymbol{f}}}_i}(k) \ne {\boldsymbol{0}}} \right.} \right\}. 理想情况下,异常检测系统应该对所有潜在的异常检测且不存在异常的误报,即{P_{\rm{F}}} = {P_{\rm{M}}} = 0,在实际中却难以实现. 当设计的{P_{\rm{M}}}越小时,往往会导致{P_{\rm{F}}}较高,反之亦然. 因此,异常检测系统通常会采取{P_{\rm{F}}}和{P_{\rm{M}}}的折中.
本文研究的悬浮系统的最低异常漏检率可被表述为{P_{\rm{F}}}固定时的最小化{P_{\rm{M}}},即考虑当{P_{\rm{F}}}固定且可接受时,如何使得{P_{\rm{M}}}最小化[24]. 根据式(8)可知,{{{\boldsymbol{g}}}_i}是{{\boldsymbol{r}}_i}(k)生成的关键, {{{\boldsymbol{g}}}_i} 的取值会影响异常检测的{P_{\rm{F}}}和{P_{\rm{M}}}. 因此,基于最低漏检率的残差生成器的设计最终转化为求得参数向量{{{\boldsymbol{g}}}_i}({{{\boldsymbol{g}}}_i} \ne {\boldsymbol{0}}),使得
\begin{gathered} \min {\alpha _i}, \\ {\rm{ s.t.}}\left\{ \begin{gathered} \sup \Pr \left\{ {{{\left( {{{\boldsymbol{g}}}_i^{\rm{T}}\left( {{{\boldsymbol{z}}}(k) - {{\bar {\boldsymbol{z}}}_{\rm{h}}}} \right)} \right)}^2} \leqslant 1} \right\} \leqslant {\alpha _i}, \\ \sup \Pr \left\{ {{{\left( {{{\boldsymbol{g}}}_i^{\rm{T}}\left( {{{\boldsymbol{z}}}(k) - {{\bar {\boldsymbol{z}}}_{\rm{h}}}} \right)} \right)}^2} > 1} \right\} \leqslant {\beta _i}, \\ \end{gathered} \right. \end{gathered} | (10) |
式中: {\alpha _i} \in (0,1) , {\beta _i} \in (0,1) ,分别为{P_{\rm{M}}}和{P_{\rm{F}}}的上界, {\beta _i} 在设计时应预先给出[25].
进一步,异常检测式(10)可以被重构为
\begin{gathered} \max a({{{\boldsymbol{g}}}_i}), \\ {\rm{ s.t.}}\left\{ \begin{gathered} {k_{\beta_i}}\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}\varSigma _{\rm{h}} {{{\boldsymbol{g}}}_i}} \leqslant 1, \\ a({{{\boldsymbol{g}}}_i}) = \frac{{\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}} {{\tilde {\boldsymbol{z}}}_{{\boldsymbol{f}}_i}} {\tilde {\boldsymbol{z}}}_{{\boldsymbol{f}}_i}^{\rm{T}} {{\boldsymbol{g}}_i}} - {k_{\beta_i}}\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}\varSigma_{\rm{h}} {{{\boldsymbol{g}}}_i}} }}{{\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}\varSigma_{{{\boldsymbol{f}}_i}} {{{\boldsymbol{g}}}_i}} }} > 0 , \\ \end{gathered} \right. \\ \end{gathered} | (11) |
式中:{k_{\beta_i}} = \sqrt {1/{\beta _i}};{\tilde{\boldsymbol{z}}}_{{\boldsymbol{f}}_i} = {\bar {\boldsymbol{z}}_{{\boldsymbol{f}}_i}} - {\bar {\boldsymbol{z}}_{\rm{h}}}.
当获得最优{{{\boldsymbol{g}}}_i}后,只需将{{{\boldsymbol{g}}}_i}乘以一个非零常数,即可满足{k_{\beta_i}}\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{\varSigma _{\rm{h}}}{{{\boldsymbol{g}}}_i}} \leqslant 1. 因此,先分析a({{{\boldsymbol{g}}}_i}) > 0时a({{{\boldsymbol{g}}}_i})的最大值.
令
\left\{ \begin{gathered} b({{{\boldsymbol{g}}}_i}) = \frac{{\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}} {{\boldsymbol{z}}_{{\boldsymbol{f}}_i}} {{\boldsymbol{z}}_{{\boldsymbol{f}}_i}^{\rm{T}}} {{{\boldsymbol{g}}}_i}} }}{{\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{{\boldsymbol{\varSigma }} _{{\boldsymbol{f}}_i}}{{{\boldsymbol{g}}}_i}} }}, \\ c({{{\boldsymbol{{\boldsymbol{g}}}}}_i}) = \frac{{{k_{\beta_i}}\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{{\boldsymbol{\varSigma }}_{\rm{h}}}{{{\boldsymbol{g}}}_i}} }}{{\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{{\boldsymbol{\varSigma }} _{{\boldsymbol{f}}_i}}{{{\boldsymbol{g}}}_i}} }}, \\ \end{gathered} \right. | (12) |
式中:{{\boldsymbol{\varSigma }} _{{\boldsymbol{f}}_i}} \ne 0;b({{{\boldsymbol{g}}}_i}) > c({{{\boldsymbol{g}}}_i}) > 0.
则d({{{\boldsymbol{g}}}_i}) = c({{{\boldsymbol{g}}}_i})/b({{{\boldsymbol{g}}}_i}) < 1,因此,可以得到式(13)的关系.
\left\{ \begin{gathered} a({{{\boldsymbol{g}}}_i}) = b({{{\boldsymbol{g}}}_i}) - c({{{\boldsymbol{g}}}_i}) = (1 - d({{{\boldsymbol{g}}}_i}))b({{{\boldsymbol{g}}}_i}), \\ e({{{\boldsymbol{g}}}_i}) = {b^2}({{{\boldsymbol{g}}}_i}) - {c^2}({{{\boldsymbol{g}}}_i}) = \frac{{1 + d({{{\boldsymbol{g}}}_i})}}{{1 - d({{{\boldsymbol{g}}}_i})}}{a^2}({{{\boldsymbol{g}}}_i}). \\ \end{gathered} \right. | (13) |
对e({{{\boldsymbol{g}}}_i})求关于a({{{\boldsymbol{g}}}_i})的导数,得
\frac{{\partial e({{{\boldsymbol{g}}}_i})}}{{\partial a({{{\boldsymbol{g}}}_i})}} = 2a({{{\boldsymbol{g}}}_i})\frac{{1 + d({{{\boldsymbol{g}}}_i})}}{{1 - d({{{\boldsymbol{g}}}_i})}} > 0. | (14) |
因此,e({{\boldsymbol{g}}_i})关于a({{{\boldsymbol{g}}}_i})是单调增加的,所以式(11)可以变为求{{{\boldsymbol{g}}}_i} = \arg \mathop {\max }\limits_{{{\boldsymbol{g}}_i} \ne 0} e({{{\boldsymbol{g}}}_i})的问题[26],即
\mathop {\max }\limits_{{{\boldsymbol{g}}_i} \ne 0} e({{{\boldsymbol{g}}}_i}) = \mathop {\max }\limits_{{{\boldsymbol{g}}_i} \ne 0} \frac{{\sqrt {{{\boldsymbol{g}}}_i^{\rm{T }} ({{\tilde {\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}} { \tilde{\boldsymbol{z}}^{\rm{T}}_{{{\boldsymbol{f}}_i}}} - k_{\beta_i}^2{{{\boldsymbol{\varSigma }}} _{\rm{h}}}) {{{\boldsymbol{g}}}_i}} }}{{\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{{{\boldsymbol{\varSigma }}} _{{{\boldsymbol{f}}_i}}}{{{\boldsymbol{g}}}_i}} }}. | (15) |
通过上述分析,基于最低漏检率的残差生成器的设计问题最终转换为求解式(15)的广义特征值和特征向量[27],具体求解步骤如算法1所示.
1) 算法 1:最优参数向量的求解
步骤 1 对\sqrt {{{{\boldsymbol{\varSigma }}} _{{{\boldsymbol{f}}_i}}}}做 SVD 分解,\sqrt {{{{\boldsymbol{\varSigma }}} _{{\boldsymbol{f}}_i}}} = {{{{\boldsymbol{U}}}}_i}\left[ { {{{{{\boldsymbol{S}}}}_i}}\;\; 0 } \right] {{{\boldsymbol{V}}}}_i^{\rm{T}},其中,\lambda _{{\rm{max}}} 为最大奇异值;{\boldsymbol{U}}_i、{\boldsymbol{S}}_i、{\boldsymbol{V}}_i 均为SVD分解后得出的矩阵;{\boldsymbol{v}}_i 为{\boldsymbol{V}}_i 中的元素向量.
步骤 2 通过求解方程组
\left\{ \begin{gathered} {{\boldsymbol{v}}}_i^{\rm{T}}({\lambda _{m,i}}{{{\boldsymbol{I}}}} - {{{\boldsymbol{S}}}}_i^{ - 1}{{{\boldsymbol{U}}}}_i^{\rm{T}}{{{\varXi}} _i}{{{{\boldsymbol{U}}}}_i}{{{\boldsymbol{S}}}}_i^{ - 1}){{{\boldsymbol{v}}}_i} = 0, \\ {\varXi _i} = {{\tilde {\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}} { {\boldsymbol{z}}^{\rm{T}}_{{{\boldsymbol{f}}_i}}} - k_{{\beta _i}}^2{{\boldsymbol{\varSigma}} _{\rm{h}}} , \\ {\lambda _{m,i}} = {\lambda _{\max }}\left\{ {{{{\boldsymbol{S}}}}_i^{ - 1}{{{\boldsymbol{U}}}}_i^{\rm{T}}{\varXi _i}{{{{\boldsymbol{U}}}}_i}{{{\boldsymbol{S}}}}_i^{ - 1}} \right\}, \\ {{\boldsymbol{v}}}_i^{\rm{T}}{{{\boldsymbol{v}}}_i} = 1, \\ \end{gathered} \right. |
得出{{{\boldsymbol{v}}}_i}.
步骤 3 依据 {{{\boldsymbol{g}}}_i} = {{{{\boldsymbol{U}}}}_i}{{{\boldsymbol{S}}}}_i^{ - 1}{{{\boldsymbol{v}}}_i}计算 {{{\boldsymbol{g}}}_i} .
步骤 4 依据 {{\hat {\boldsymbol{g}}}_i} = \dfrac{{{{\boldsymbol{g}}}_i}}{{{k_{{\beta _i}}}\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{{{{\boldsymbol{\varSigma}}}} _{\rm{h}}}{{{\boldsymbol{g}}}_i}} }}计算出满足式(10)的最优参数向量{{\hat {\boldsymbol{g}}}_i}.
在得出最优参数向量后,第i个异常的参数化残差生成器可以通过式(16)构建.
{{{\boldsymbol{r}}}_i}(k) = {{\hat {\boldsymbol{g}}}_i}^{\rm{T}}{{{{\boldsymbol{\varPsi}} }}_s}^ \bot \left[ {\begin{array}{*{20}{c}} {{{{\boldsymbol{u}}}_s}(k)} \\ {{{{\boldsymbol{y}}}_s}(k)} \end{array}} \right]. | (16) |
令{{\hat {\boldsymbol{r}}}_{i,{\rm{h}}}} = {{\hat {\boldsymbol{g}}}_i}^{\rm{T}}{\bar {\boldsymbol{z}}_{\rm{h}}},异常检测器的残差评估函数和异常检测阈值为
\left\{ \begin{gathered} {\boldsymbol{J}}({{{\boldsymbol{r}}}_i}) = {\left( {{{{\boldsymbol{r}}}_i}\left( k \right) - {{{\hat {\boldsymbol{r}}}}_{i,h}}} \right)^2}, \\ {J_{{\rm{th}}}} = 1. \\ \end{gathered} \right. | (17) |
完整的异常检测见算法 2
2) 算法 2:磁悬浮列车悬浮系统异常检测算法
① 离线辨识
步骤 1 采集无异常情况下的输入输出数据,基于子空间辨识法离线辨识出矩阵{{{{\boldsymbol{\varPsi}} }}_{\boldsymbol{s}}}^ \bot.
步骤 2 利用无异常的输入输出数据,计算得到{\bar {\boldsymbol{z}}_{\rm{h}}}和{{\boldsymbol{\varSigma}} _{\rm{h}}},构建健康状况下的置信集{S_{\rm{h}}}.
步骤 3 建立M个异常的置信集{S_{{{\boldsymbol{f}}_i}}} . 当有M个异常时的输入输出数据,计算每个异常的 {\bar {\boldsymbol{z}}_{{{\boldsymbol{f}}_i}}} 和{{{{\boldsymbol{\varSigma}} }}_{{\boldsymbol{f}}_i}},建立 {S _{{{\boldsymbol{f}}_i}}} . 如果无异常数据,预设每个异常的{\bar {\boldsymbol{z}}_{{{\boldsymbol{f}}_i}}}和{{{{\boldsymbol{\varSigma}}}} _{{{\boldsymbol{f}}}_i}},依此建立 {S _{{{\boldsymbol{f}}_i}}} .
步骤 4 确定每个异常下的{\beta _i},根据算法1计算出{{\hat {\boldsymbol{g}}}_i}和{{\hat {\boldsymbol{r}}}_{i,h}}.
② 在线异常检测
步骤 5 在线收集时间间隔为\left[ {k-s,k} \right]的输入输出数据,并构造数据矩阵\left[ {{{{\boldsymbol{u}}}_s}(k)}\quad{{{{\boldsymbol{y}}}_s}(k)} \right]^{\rm{T}}.
步骤 6 依据式(16)计算{{{\boldsymbol{r}}}_i}\left( k \right),依据式(17)执行异常检测并输出异常检测的结果.
步骤 7 令k = k + 1,重复步骤5和步骤6,直到中止时刻.
本节结合长沙磁浮快线某辆磁浮列车的悬浮系统的健康数据和异常数据,验证所提出的基于最低漏检率的残差生成器在悬浮系统异常检测中的有效性.
以长沙磁浮快线的运营磁浮列车为例,如图1所示,作为典型的常导电磁型磁浮列车[28],通常选用间隙作为其悬浮系统异常检测的检测变量,通过比较间隙值和异常阈值的大小来判断是否出现异常. 磁浮列车悬浮系统的工程异常阈值设置为标准间隙值 ±4 mm. 该线路的标准悬浮工作间隙为9 mm,工程异常阈值分别设为5 mm和13 mm.
磁浮列车悬浮系统经常会遭受到3类异常:间隙突变异常、砸轨异常和加速度传感器异常. 间隙突变异常的现象是指在运行中间隙有明显的波动,但没有超过工程异常阈值,此时加速度、电流也对应有一定的波动,出现两个尖峰,但电压几乎没有波动. 该异常通常是由于磁浮列车经过轨道接缝引起的,如图2所示.
砸轨异常的现象是间隙有剧烈波动,在一段时间内间隙值超过工程异常阈值并且需要经历5 s左右的恢复期,此时加速度、电流、电压均在一段时间内出现剧烈波动,如图3所示.
加速度计异常的现象为间隙没有超过工程异常阈值,也没有较为明显的波动,但加速度的波动十分明显且变为负值,电压也出现两个极小值尖峰. 由于实际加速度计存在着1g的偏置,当加速度计输出信号为负时,表明此刻的加速度已经低于 - 10\;{\rm{m/{s^2}}},如图4所示.
以磁浮列车悬浮系统的第13个悬浮点的运行数据作为检测对象,如图5所示.
列车正线行驶时存在着站内悬浮和站间行驶两个运行工况,考虑到上述3类异常均出现在站间行驶运行工况,因此,以速度为划分基准,分段提取站间行驶工况的运行数据. 本节选择一个无异常站间行驶数据段、3个包含间隙异常的数据段、2个包含砸轨异常的运行数据段和包含3个加速度传感器异常的数据段作为检测数据.
选择悬浮系统的电流数据作为输入数据,间隙数据和加速度数据作为输出数据,如图6所示. 依据正常运行时的输入输出数据建立{S_{\rm{h}}}并离线辨识出{{{{\boldsymbol{\varPsi}} }}_s}^ \bot.
依据3种异常的输入输出数据建立异常置信集 {S _{{{\boldsymbol{f}}_i}}} . 出于工程运营需求,设期望的故障误报率{P_{\rm{F}}} = 5{\text{%}} ,执行算法1可得
{{\hat {\boldsymbol{g}}}_i} =(0.257\;5, 0.257\;5,\cdots, 0.257\;5)_{1\times9}. | (18) |
基于最低漏检率的残差生成器为
{{{\boldsymbol{r}}}_i}(k) = 0.257\;5 {{{{\boldsymbol{\varPsi}} }}_s}^ \bot \left[ {\begin{array}{*{20}{c}} {{{{\boldsymbol{u}}}_s}(k)} \\ {{{{\boldsymbol{y}}}_s}(k)} \end{array}} \right]. | (19) |
依据式(17)算出9个检测数据段的J({{\boldsymbol{r}}_i})并送入检测单元. 图7为无异常数据段的检测结果. 从图中可以看出,当系统无异常时,J({{\boldsymbol{r}}_i})小于0.4,均在异常检测阈值{J_{{\rm{th}}}}之下,不存在异常误检.
含间隙突变异常的3个数据段的检测结果如图8所示. 在每个数据段,大部分数据都低于异常检测阈值{J_{{\rm{th}}}},但是在出现异常的时刻,存在着J({{{\boldsymbol{r}}}_i})大于{J_{{\rm{th}}}}的情况,个数分别为8、6和5,说明该方案能够实现对间隙异常的检测. 此外,间隙突变异常的J({{{\boldsymbol{r}}}_i})的峰值略高于检测阈值{J_{{\rm{th}}}},一般在2以下,且不存在第二个大于{J_{{\rm{th}}}}的次尖峰. 表1为间隙突变异常的检测位置与异常实际发生位置的对比,从表中可以看出,间隙突变异常1的检测时刻落后实际异常发生时刻 0.2 s,其他两个数据段可以及时检测出异常.
异常 | 检测位置/s | 实际位置/s |
间隙异常 1 | 5.9 | 5.7 |
间隙异常 2 | 568.7 | 568.5 |
间隙异常 3 | 21.1 | 21.1 |
含砸轨异常的两个数据段的检测结果如图9所示. 从图中可以看出,该方法能够实现对砸轨异常的检测. 图9(a)中有20个样本的J({{{\boldsymbol{r}}}_i})>{J_{{\rm{th}}}},图9(b)中有20个样本的J({{{\boldsymbol{r}}}_i})>{J_{{\rm{th}}}},此外,砸轨异常位置的J({{{\boldsymbol{r}}}_i})存在着两个高于{J_{{\rm{th}}}}的峰值. 表2为砸轨异常的检测位置与异常实际发生位置对比. 从表2可以看出,砸轨异常的检测时刻落后于实际异常的发生时刻,落后约0.2 s. 结合图3可以发现,造成这一现象的主要原因是砸轨异常发生时刻数据变化不明显,间隙和加速度都有一个缓慢减小的阶段,没有大范围变化.
异常 | 检测位置/s | 实际位置/s |
砸轨异常 1 | 7.1 | 7.0 |
砸轨异常 2 | 20.2 | 20.0 |
含加速度计异常的3个数据段的检测结果如图10所示. 在每个数据段,大部分数据都低于异常检测阈值{J_{{\rm{th}}}},但是在出现加速度计异常的时刻,存在着J({{{\boldsymbol{r}}}_i})>{J_{{\rm{th}}}}的情况,个数分别为5、7和4,且J({{{\boldsymbol{r}}}_i})的值较大,说明该方案能够实现对加速度计异常的检测. 表3为加速度计异常的检测位置与异常实际发生位置对比. 从表3中可以看出,一旦出现加速度计异常,本文所提出的方法可以立即检测出该异常,及时性较好.
异常 | 检测位置/s | 实际位置/s |
加速度计异常 1 | 52.2 | 52.2 |
加速度计异常 2 | 52.1 | 52.1 |
加速度计异常 3 | 40.5 | 40.5 |
经过上述分析,基于残差生成器的异常检测方法可以用于磁悬浮列车悬浮系统的异常检测. 基于最低漏检率的异常检测方法能够实现磁悬浮列车悬浮系统3类典型异常的完全检测,不存在漏检的情况且不存在误检.
本文针对电磁悬浮系统的异常检测问题开展了研究,以基于参数化残差的异常检测为核心,依据数据的均值和方差建立置信集,设计了基于最低异常漏检率的残差生成器,依据生成的残差数据在线进行异常检测. 最后,以长沙磁浮快线的运行数据为例,对悬浮系统常见的3类异常进行检测,实验结果表明,在给定异常误报率为5%时,本文所提出的方法能够实现对3类典型异常的全部检测,没有异常的漏检和误检. 检测时间对比结果进一步说明,本文提出的方法检测实时性较好.
[1] |
苏洋. 公铁两用双层桥梁风屏障气动机理及优化研究[D]. 成都: 西南交通大学, 2017.
|
[2] |
韩旭,彭栋,向活跃,等. 横风作用下高速铁路桥梁全封闭声屏障气动特性的风洞试验研究[J]. 铁道建筑,2019,59(7): 151-155. doi: 10.3969/j.issn.1003-1995.2019.07.35
HAN Xu, PENG Dong, XIANG Huoyue, et al. Research on wind tunnel tests for aerodynamic characteristics of closed noise barriers on high speed railway bridges under crosswinds[J]. Railway Engineering, 2019, 59(7): 151-155. doi: 10.3969/j.issn.1003-1995.2019.07.35
|
[3] |
孙延国,廖海黎,李明水. 基于节段模型试验的悬索桥涡振抑振措施[J]. 西南交通大学学报,2012,47(2): 218-223,264. doi: 10.3969/j.issn.0258-2724.2012.02.008
SUN Yanguo, LIAO Haili, LI Mingshui. Mitigation measures of vortex-induced vibration of suspension bridge based on section model test[J]. Journal of Southwest Jiaotong University, 2012, 47(2): 218-223,264. doi: 10.3969/j.issn.0258-2724.2012.02.008
|
[4] |
SIMIU E S R H. Wind effects on structures[M]. New York: Wiley, 1986.
|
[5] |
EHSAN F, SCANLAN R H. Vortex-induced vibrations of flexible bridges[J]. Journal of Engineering Mechanics, 1990, 116(6): 1392-411. doi: 10.1061/(ASCE)0733-9399(1990)116:6(1392)
|
[6] |
HARTLEN R T, CURRIE IAIN G. Lift-oscillator model of vortex-induced vibration[J]. Journal of the Engineering Mechanics Division, 1970, 96(5): 577-91. doi: 10.1061/JMCEA3.0001276
|
[7] |
周帅,陈克坚,陈政清,等. 大跨桥梁涡激共振幅值估算方法的理论基础与应用[J]. 高速铁路技术,2019,10(5): 25-31. doi: 10.12098/j.issn.1674-8247.2019.05.006
ZHOU Shuai, CHEN Kejian, CHEN Zhengqing, et al. Theoretical basis and practical applications of various vortex-induced vibration amplitudes estimation methods for large-span bridges[J]. High Speed Railway Technology, 2019, 10(5): 25-31. doi: 10.12098/j.issn.1674-8247.2019.05.006
|
[8] |
张志田,陈政清. 桥梁节段与实桥涡激共振幅值的换算关系[J]. 土木工程学报,2011,44(7): 77-82. doi: 10.15951/j.tmgcxb.2011.07.009
ZHANG Zhitian, CHEN Zhengqing. Similarity of amplitude of sectional model to that of full bridge in the case of vortex-induced resonance[J]. China Civil Engineering Journal, 2011, 44(7): 77-82. doi: 10.15951/j.tmgcxb.2011.07.009
|
[9] |
周奇,孟晓亮,朱乐东. 基于非线性涡激力广义模型的涡振幅值换算[J]. 土木工程学报,2020,53(10): 82-88. doi: 10.15951/j.tmgcxb.2020.10.008
ZHOU Qi, MENG Xiaoliang, ZHU Ledong. Amplitude conversion of vortex-induced vibration based on generalized model of nonlinear vortex-induced force[J]. China Civil Engineering Journal, 2020, 53(10): 82-88. doi: 10.15951/j.tmgcxb.2020.10.008
|
[10] |
SUN Y G, LI M S, LIAO H L. Nonlinear approach of vortex-induced vibration for line-like structures[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2014, 124: 1-6. doi: 10.1016/j.jweia.2013.10.011
|
[11] |
IRWIN P. Full aeroelastic model tests[M]. [S.l.]: Routledge, 2017: 125-35.
|
[12] |
HJORTH-HANSEN E. Section model tests[M]. [S.l.]: Routledge, 2017: 95-112.
|
[13] |
秦浩,廖海黎,李明水. 变截面连续钢箱梁桥典型施工阶段涡激振动[J]. 西南交通大学学报,2014,49(5): 760-765,786. doi: 10.3969/j.issn.0258-2724.2014.05.003
QIN Hao, LIAO Haili, LI Mingshui. Vortex-induced vibration of continuous steel box-girder bridge with variable cross-sections at typical erection stages[J]. Journal of Southwest Jiaotong University, 2014, 49(5): 760-765,786. doi: 10.3969/j.issn.0258-2724.2014.05.003
|
[14] |
DUAN J L, HUANG W P. CFD-based numerical analysis of a variable cross-section cylinder[J]. Journal of Ocean University of China, 2014, 13(4): 584-588. doi: 10.1007/s11802-014-2048-0
|
[15] |
陈政清. 工程结构的风致振动、稳定与控制[M]. 北京: 科学出版社, 2013.
|
[16] |
王新敏. ANSYS结构动力分析与应用 [M]. 北京: 人民交通出版社, 2014.
|
[1] | LI Qi, AI Yuxuan, SUN Cai, QIU Yibin, CHEN Weirong. Optimal Reconfiguration of Distribution Network Based on Backtracking Search Algorithm Under the Background of Non-cooperative Game Theory[J]. Journal of Southwest Jiaotong University, 2024, 59(2): 438-446. doi: 10.3969/j.issn.0258-2724.20210547 |
[2] | LIU Taofeng, LI Yiyuan, LI Wei, ZHOU Zhongrong. Influence of Surface Feature Height of Deterministic Texture on Tactile Perception of Fingertip[J]. Journal of Southwest Jiaotong University, 2020, 55(2): 372-378. doi: 10.3969/j.issn.0258-2724.20180038 |
[3] | LI Tian, QIN Deng, AN Chao, ZHANG Jiye. Effect of Computational Grid on Uncertainty in Train Aerodynamics[J]. Journal of Southwest Jiaotong University, 2019, 54(4): 816-822. doi: 10.3969/j.issn.0258-2724.20180503 |
[4] | HU Guangdi, WANG Guohui, LUO Huiyu, ZHOU Ke, LANG Xiaoyue. Robust Control of the Air Conditioning System of an Electric Vehicle with Actuator Fault[J]. Journal of Southwest Jiaotong University, 2018, 53(2): 351-358. doi: 10.3969/j.issn.0258-2724.2018.02.018 |
[5] | SUN Yue, TAN Ruoxi, TANG Chunsen, WANG Zhihui, DAI Lin. Optimized Design of New Coupling Mechanism for Electric Vehicles[J]. Journal of Southwest Jiaotong University, 2018, 53(5): 1078-1086. doi: 10.3969/j.issn.0258-2724.2018.05.027 |
[6] | SONG Kai, LIU Dan, LIU Jian. Multiple Model Analysis for Studying Groundwater Uncertainties[J]. Journal of Southwest Jiaotong University, 2018, 53(3): 574-581. doi: 10.3969/j.issn.0258-2724.2018.03.019 |
[7] | SU Yugang, ZHANG Shuai, XU Yong, TANG Chunsen. Design and Switching Control of Power Supply Coils Applied to ICPT-Based Electric Vehicles[J]. Journal of Southwest Jiaotong University, 2016, 29(1): 168-176. doi: 10.3969/j.issn.0258-2724.2016.01.024 |
[8] | ZHANG Junfeng, GE Tengteng, CHEN Qiang, WANG Fei. 4D Trajectory Prediction and Uncertainty Analysis for Departure Aircraft[J]. Journal of Southwest Jiaotong University, 2016, 29(4): 800-806. doi: 10.3969/j.issn.0258-2724.2016.04.027 |
[9] | LI Shanmei, XU Xiaohao, WANG Fei. Risk Prediction Model and Methodology of Airport Congestion Based on Probabilistic Demand[J]. Journal of Southwest Jiaotong University, 2013, 26(1): 154-159. doi: 10.3969/j.issn.0258-2724.2013.01.024 |
[10] | p style, ZHANG Ying, XIE Hua, HU Minghua. Optimization Model for Area Traffic Flow Management under Constraint of Deterministic Airspace Capacity[J]. Journal of Southwest Jiaotong University, 2012, 25(2): 348-354. doi: 10.3969/j.issn.0258-2724.2012.02.028 |
[11] | TIAN Wen, HU Minghua. Airspace Sector Probabilistic Traffic Demand Prediction Model[J]. Journal of Southwest Jiaotong University, 2011, 24(2): 340-346. doi: 10.3969/j.issn.0258-2724.2011.02.027 |
[12] | ZHANG Yougang, LUO Zhiyong. LMI-Based Robust Stabilization of Linear Uncertain Systems with Input Time-Delay[J]. Journal of Southwest Jiaotong University, 2009, 22(5): 682-687. |
[13] | WANG Cuihong, HUANG Tianmin. Robust Strictly Passive Control of Uncertain Descriptor Systems[J]. Journal of Southwest Jiaotong University, 2009, 22(6): 860-864,911. |
[14] | DUAN Zhongdong, ZHOU Qingsheng. Hierarchical Decentralized Robust Control Method of Uncertainty Structure Vibration[J]. Journal of Southwest Jiaotong University, 2009, 22(2): 155-159. |
[15] | ZHANGKe-yue, CAO Deng-qing. Stability of Uncertain Linear Systems with Multiple Time-Delays[J]. Journal of Southwest Jiaotong University, 2004, 17(2): 176-180. |
[16] | ZHANG Shen-feng, WUYu-hua, GUOJun-peng. Robust Stability Analysis for Interval DEA[J]. Journal of Southwest Jiaotong University, 2004, 17(2): 253-236. |
[17] | HANLing, ZHAOLian-wen. Static Game Model of Incomplete Information[J]. Journal of Southwest Jiaotong University, 2003, 16(3): 359-362. |
[18] | SHI Zhi-ping, JIN Fan. Deterministic Interleaver Design for High-Rate Turbo Codes[J]. Journal of Southwest Jiaotong University, 2002, 15(5): 544-547. |
[19] | MAXian-min, CHENQuan-shi. The Design and Simulation of Electric Vehicle Drive System with High Frequency AC Pulse Density Modulation[J]. Journal of Southwest Jiaotong University, 2002, 15(1): 81-84. |
[20] | Cao Dengqing. Lateral Stability Analysis for Rail Vehicle Dynamic Models with Uncertain Parameters[J]. Journal of Southwest Jiaotong University, 1999, 12(3): 253-258. |
异常 | 检测位置/s | 实际位置/s |
间隙异常 1 | 5.9 | 5.7 |
间隙异常 2 | 568.7 | 568.5 |
间隙异常 3 | 21.1 | 21.1 |
异常 | 检测位置/s | 实际位置/s |
砸轨异常 1 | 7.1 | 7.0 |
砸轨异常 2 | 20.2 | 20.0 |
异常 | 检测位置/s | 实际位置/s |
加速度计异常 1 | 52.2 | 52.2 |
加速度计异常 2 | 52.1 | 52.1 |
加速度计异常 3 | 40.5 | 40.5 |