
Citation: | ZHOU Linli, TAN Ping, TENG Xiaofei. Lateral Performance of Semicircular Corrugated Steel Plate Shear Wall with Edge Stiffeners[J]. Journal of Southwest Jiaotong University, 2021, 56(3): 550-557. doi: 10.3969/j.issn.0258-2724.20190319 |
随着中速磁悬浮列车的广泛商业化应用[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)⩽Jth,健康. | (5) |
传统的异常检测方案通常采用基于假设检验的方法进行残差评估,即假设残差的分布知识是准确已知的,或者可以从历史数据中估计出来[22]. 然而,对于磁浮列车悬浮系统而言,使用历史过程数据获得高度可信的概率分布估计比较繁琐,经验估计与真实概率分布的偏差可能会导致异常检测结果不可靠[23].
尽管难以获得悬浮系统的高度可信概率分布,但在无异常时的z(k)可以视为属于同一个集合. 考虑到悬浮系统运行数据的统计特性,可以依据数据的均值和协方差矩阵对其进行建模.
设无异常情况下,z(k)满足
Sh={Sz∈ζz|E=ˉzh,V=Σh}, | (6) |
式中:Sh为无异常数据的置信集;ζz为所有z(k)的分布集合;ˉzh和Σh为无异常时状态检测变量的均值和方差,ˉzh和Σh是通过对健康运行数据进行计算得来的,通常可以选择刚投入运行时无异常的运行数据,也可以是某一次无异常运行时的数据;E为z(k)的均值;V为z(k)的方差.
在悬浮系统运行初期,由于缺乏异常数据,难以分析出异常的数据分布特性. 因此,通常会依据正常运行时的数据特点预先设计一个异常检测系统,以实现对大部分潜在异常的检测. 假设系统运行时会出现M个典型异常,且发生第i个异常时的异常信号设为fi(k),i=1,2,⋯,M. 为不失一般性,假设每个异常发生的概率相同. 采用均值和协方差矩阵对异常进行建模. 设异常情况时的z(k)满足
Sfi={Sz∈ζz|E=ˉzfi,V=Σfi}, | (7) |
式中:Sfi为异常数据的置信集;ˉzfi为zfi(k)的均值,zfi(k)为第i个异常时的状态检测变量;Sfi为zfi(k)的方差.
当有可用的异常数据时,可以直接使用该数据建立置信集. 若无异常数据或异常数据少时,则可以自定义M个异常情况下的均值和协方差,进而构建置信集,当有新可用的异常数据时,记录数据并离线更新置信集.
针对第 i个异常,基于异常检测的残差信号为
ri(k)=gTiz(k), | (8) |
式中:gi为第i个异常时的参数向量.
因此,第i个异常时的残差评估函数J(r)与阈值Jth为
{J(r)=(ri(k)−ri,h)2,ri,h=gTiˉzh,Jth=1, | (9) |
式中:ri,h为第i个无异常场景的残差.
定义异常误检率PF=Pr,定义异常漏检率为{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类典型异常的全部检测,没有异常的漏检和误检. 检测时间对比结果进一步说明,本文提出的方法检测实时性较好.
张彬彬,曹万林,潘毅,等. 双向单排配筋T形剪力墙抗震性能试验研究[J]. 土木建筑与环境工程,2011,33(增刊1): 203-208.
ZHANG Binbin, CAO Wanlin, PAN Yi, et al. Experimental study on seismic performance of T-shear wall structure with bidirectional single row of steel bars[J]. Journal of Civil Architectural and Environmental Engineering, 2011, 33(S1): 203-208.
|
张彬彬,曹万林,张建伟,等. 双向单排配筋Z形剪力墙抗震性能试验研究[J]. 工程抗震与加固改造,2011,33(5): 45-52. doi: 10.3969/j.issn.1002-8412.2011.05.008
ZHANG Binbin, CAO Wanlin, ZHANG Jianwei, et al. Experimental study on seismic performance of Z-shear wall structure with bidirectional single row of steel bars[J]. Earthquake Resistant Engineering and Retrofitting, 2011, 33(5): 45-52. doi: 10.3969/j.issn.1002-8412.2011.05.008
|
American Institute of Steel Construction. Specification for structural steel buildings: ANSI/AISC360-10[S]. Chicago: [s.n.], 2010.
|
BERMAN J, BRUNEAU M. Plastic analysis and design of steel plate shear walls[J]. Journal of structural Engineering, 2003, 129(11): 1448-1456.
|
张哲,李国强,孙飞飞. 波纹腹板钢梁受弯稳定性研究[J]. 工程力学,2011,28(8): 77-82.
ZHANG Zhe, LI Guoqiang, SUN Feifei. Research on flexural stability of the H-beams with corrugated webs[J]. Engineering Mechanics, 2011, 28(8): 77-82.
|
李国强,范昕,孙飞飞,等. 波纹腹板H型钢梁柱铰接节点静力性能试验研究及有限元分析[J]. 建筑结构学报,2013,34(4): 107-114.
LI Guoqiang, FAN xin, SUN Feifei, et al. Investigation of static behavior of pinned connections of H-beams with corrugated webs[J]. Journal of Building Structures, 2013, 34(4): 107-114.
|
清华大学, 巴特勒(上海)有限公司. 波浪腹板钢结构应用技术规程: CECS290—2011[S]. 北京: 中国计划出版社, 2011.
|
中国建筑标准设计研究院有限公司. 门式刚架轻型房屋钢结构技术规范: GB 51022—2015[S]. 北京: 中国建筑工业出版社, 2015.
|
BOTROS R B G. Nonlinear finite element analysis of corrugated steel plate shear walls[R]. Calgary: University of Calgary, 2006.
|
兰银娟. 折板钢板剪力墙抗侧力结构理论研究[D]. 西安: 西安建筑科技大学, 2006.
|
李峰,薛港,问晓朋,等. 斜置单层折板钢板剪力墙的屈曲分析[J]. 水利与建筑工程学报,2011,9(4): 89-92. doi: 10.3969/j.issn.1672-1144.2011.04.021
LI Feng, XUE Gang, WEN Xiaopeng, et al. Elastic buckling analysis for inclined and single fold plate steel shear wall[J]. Journal of Water Resources and Architectural Engineering, 2011, 9(4): 89-92. doi: 10.3969/j.issn.1672-1144.2011.04.021
|
EMAMI F, MOFID M, VAFAI A. Experimental study on cyclic behavior of trapezoidally corrugated steel shear walls[J]. Engineering Structures, 2013, 48: 750-762.
|
魏瑶,谭平,李洋,等. 侧边加劲波纹钢板墙弹性剪切屈曲分析[J]. 地震工程与工程振动,2015,35(3): 199-209.
WEI Yao, TAN Ping, LI Yang, et al. Elastic shear buckling analysis of corrugated steel plate wall with edge stiffeners[J]. Earthquake Engineering and Engineering Dynamics, 2015, 35(3): 199-209.
|
谭平,魏瑶,李洋,等. 波纹钢板剪力墙抗震性能试验研究[J]. 土木工程学报,2018,51(5): 8-15.
TAN Ping, WEI Yao, LI Yang, et al. Experimental investigation on performance of corrugated steel plate shear wall[J]. China Civil Engineering Journal, 2018, 51(5): 8-15.
|
孙军浩. 波纹钢板剪力墙的抗侧及抗震性能研究[D]. 天津: 天津大学, 2016.
|
赵秋红,李楠,孙军浩. 波纹钢板剪力墙结构的抗侧性能分析[J]. 天津大学学报(自然科学与工程技术版),2016,49(增刊1): 152-160.
ZHAO Qiuhong, LI Nan, SUN Junhao. Analysis on lateral performance of sinusoidally corrugated steel plate shear walls[J]. Journal of Tianjin University (Science and Technology), 2016, 49(S1): 152-160.
|
王威,张龙旭,苏三庆,等. 波形钢板剪力墙抗震性能试验研究[J]. 建筑结构学报,2018,39(5): 36-44.
WANG Wei, ZHANG Longxu, SU Sanqing, et al. Experimental research on seismic behavior of corrugated steel plate shear wall[J]. Journal of Building Structures, 2018, 39(5): 36-44.
|
郭彦林, 童精中, 姜子钦. 波形腹板钢结构设计原理与应用[M]. 北京: 科学出版社, 2015: 15-34.
|
龙驭球, 包世华. 结构力学[M]. 北京: 高等教育出版社, 2001: 10-26.
|
孙训方, 方孝淑, 关来泰. 材料力学[M]. 5版. 北京: 高等教育出版社, 2009: 5-21.
|
哈尔滨工业大学. 钢板剪力墙技术规程: JGJ/T 380—2015[S]. 北京: 中国建筑工业出版社, 2016.
|
[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 |