Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

基于剪压强度准则的UHPC键齿接缝抗剪承载力

徐略勤 徐意宏 陈亮 江林松 周建庭

陈维荣, 燕雨, 李奇. 基于状态机的燃料电池混合动力系统控制策略[J]. 西南交通大学学报, 2019, 54(4): 663-670. doi: 10.3969/j.issn.0258-2724.20170279
引用本文: 徐略勤, 徐意宏, 陈亮, 江林松, 周建庭. 基于剪压强度准则的UHPC键齿接缝抗剪承载力[J]. 西南交通大学学报. doi: 10.3969/j.issn.0258-2724.20230042
CHENG Weirong, YAN Yu, LI Qi. Control Strategy Based on State Machine for Fuel Cell Hybrid Power System[J]. Journal of Southwest Jiaotong University, 2019, 54(4): 663-670. doi: 10.3969/j.issn.0258-2724.20170279
Citation: XU Lueqin, XU Yihong, CHEN Liang, JIANG Linsong, ZHOU Jianting. Shear Bearing Capacity of Ultra-High Performance Concrete Keyed Joints Based on Shear-Compression Strength Criterion[J]. Journal of Southwest Jiaotong University. doi: 10.3969/j.issn.0258-2724.20230042

基于剪压强度准则的UHPC键齿接缝抗剪承载力

doi: 10.3969/j.issn.0258-2724.20230042
基金项目: 国家自然科学基金项目(51978113);重庆市英才计划(cstc2022ycjh-bgzxm0133)
详细信息
    作者简介:

    徐略勤(1983—),男,教授,博士,研究方向为桥梁工程,E-mail:xulueqin@163.com

  • 中图分类号: U443.3

Shear Bearing Capacity of Ultra-High Performance Concrete Keyed Joints Based on Shear-Compression Strength Criterion

  • 摘要:

    为探明超高性能混凝土(UHPC)键齿接缝的抗剪承载力及其计算方法,以接缝类型、侧压应力为参数开展5个UHPC键齿接缝和1个UHPC平接缝的足尺模型试验,研究UHPC接缝试件的破坏模式和抗剪承载力变化规律;基于混凝土剪压强度准则,采用八面体应力公式推导UHPC键齿接缝的抗剪承载力计算公式,并与试验结果进行对比验证;通过收集62个UHPC键齿接缝的已有试验数据,对本文所提的计算公式进行精度分析. 研究结果表明:UHPC键齿接缝因键齿根部出现竖向贯通主裂缝而破坏,具有明显的脆性特征;在相同侧压应力下,胶接缝试件的极限荷载比干接缝大21.3%,当侧压应力由3 MPa增大至12 MPa时,胶接缝试件的极限荷载增大了26.9%;本文计算公式可较准确地预测UHPC键齿接缝的抗剪承载力,且离散性较低,计算结果总体偏保守,其中本文5个键齿接缝试件的平均绝对误差为11%,62个已有试件的平均绝对误差为21%.

     

  • 随着中速磁悬浮列车的广泛商业化应用[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的系统输出;ABCD为一定维数的系统参数矩阵;BfDf为异常参数矩阵;f(k)为异常向量;w(k)为系统的过程噪声;v(k)为系统的输出噪声.

    则状态检测量z(k)

    {z(k)=Ψs[us(k)ys(k)],us(k)=[uT(ks)uT(ks+1)uT(k)]T,ys(k)=[yT(ks)yT(ksH)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是通过对健康运行数据进行计算得来的,通常可以选择刚投入运行时无异常的运行数据,也可以是某一次无异常运行时的数据;Ez(k)的均值;Vz(k)的方差.

    在悬浮系统运行初期,由于缺乏异常数据,难以分析出异常的数据分布特性. 因此,通常会依据正常运行时的数据特点预先设计一个异常检测系统,以实现对大部分潜在异常的检测. 假设系统运行时会出现M个典型异常,且发生第i个异常时的异常信号设为fi(k)i=1,2,,M. 为不失一般性,假设每个异常发生的概率相同. 采用均值和协方差矩阵对异常进行建模. 设异常情况时的z(k)满足

    Sfi={Szζz|E=ˉzfi,V=Σfi}, (7)

    式中:Sfi为异常数据的置信集;ˉzfizfi(k)的均值,zfi(k)为第i个异常时的状态检测变量;Sfizfi(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}) > 0a({{{\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 0b({{{\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.

    图  1  长沙磁浮快线
    Figure  1.  Changsha Maglev Express

    磁浮列车悬浮系统经常会遭受到3类异常:间隙突变异常、砸轨异常和加速度传感器异常. 间隙突变异常的现象是指在运行中间隙有明显的波动,但没有超过工程异常阈值,此时加速度、电流也对应有一定的波动,出现两个尖峰,但电压几乎没有波动. 该异常通常是由于磁浮列车经过轨道接缝引起的,如图2所示.

    图  2  间隙突变异常时间隙、电流、加速度和电压
    Figure  2.  Gap, current, acceleration, and voltage at gap mutation anomaly

    砸轨异常的现象是间隙有剧烈波动,在一段时间内间隙值超过工程异常阈值并且需要经历5 s左右的恢复期,此时加速度、电流、电压均在一段时间内出现剧烈波动,如图3所示.

    图  3  砸轨异常时的间隙、电流、加速度和电压
    Figure  3.  Gap, current, acceleration, and voltage at rail smashing anomaly

    加速度计异常的现象为间隙没有超过工程异常阈值,也没有较为明显的波动,但加速度的波动十分明显且变为负值,电压也出现两个极小值尖峰. 由于实际加速度计存在着1g的偏置,当加速度计输出信号为负时,表明此刻的加速度已经低于 - 10\;{\rm{m/{s^2}}},如图4所示.

    图  4  加速度计异常时的间隙、电流、加速度和电压
    Figure  4.  Gap, current, acceleration, and voltage at accelerometer anomaly

    以磁浮列车悬浮系统的第13个悬浮点的运行数据作为检测对象,如图5所示.

    图  5  第13个悬浮点的运行数据
    Figure  5.  Operational data of the 13th suspension point

    列车正线行驶时存在着站内悬浮和站间行驶两个运行工况,考虑到上述3类异常均出现在站间行驶运行工况,因此,以速度为划分基准,分段提取站间行驶工况的运行数据. 本节选择一个无异常站间行驶数据段、3个包含间隙异常的数据段、2个包含砸轨异常的运行数据段和包含3个加速度传感器异常的数据段作为检测数据.

    选择悬浮系统的电流数据作为输入数据,间隙数据和加速度数据作为输出数据,如图6所示. 依据正常运行时的输入输出数据建立{S_{\rm{h}}}并离线辨识出{{{{\boldsymbol{\varPsi}} }}_s}^ \bot.

    图  6  辨识所用的输入输出数据
    Figure  6.  Input and output data used for identification

    依据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}}}}之下,不存在异常误检.

    图  7  正常运行数据段的检测结果
    Figure  7.  Detection results of normal operational data segment

    含间隙突变异常的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,其他两个数据段可以及时检测出异常.

    图  8  间隙突变异常的检测结果
    Figure  8.  Detection results of gap mutation anomalies
    表  1  间隙突变异常的检测位置与异常实际发生位置
    Table  1.  Detection position and actual occurrence position of gap mutation anomalies
    异常检测位置/s实际位置/s
    间隙异常 15.95.7
    间隙异常 2568.7568.5
    间隙异常 321.121.1
    下载: 导出CSV 
    | 显示表格

    含砸轨异常的两个数据段的检测结果如图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可以发现,造成这一现象的主要原因是砸轨异常发生时刻数据变化不明显,间隙和加速度都有一个缓慢减小的阶段,没有大范围变化.

    图  9  砸轨异常的检测结果
    Figure  9.  Detection results of rail smashing anomalies
    表  2  砸轨异常的检测位置与异常实际发生位置
    Table  2.  Detection position and actual occurrence position of rail smashing anomalies
    异常检测位置/s实际位置/s
    砸轨异常 17.17.0
    砸轨异常 220.220.0
    下载: 导出CSV 
    | 显示表格

    含加速度计异常的3个数据段的检测结果如图10所示. 在每个数据段,大部分数据都低于异常检测阈值{J_{{\rm{th}}}},但是在出现加速度计异常的时刻,存在着J({{{\boldsymbol{r}}}_i})>{J_{{\rm{th}}}}的情况,个数分别为5、7和4,且J({{{\boldsymbol{r}}}_i})的值较大,说明该方案能够实现对加速度计异常的检测. 表3为加速度计异常的检测位置与异常实际发生位置对比. 从表3中可以看出,一旦出现加速度计异常,本文所提出的方法可以立即检测出该异常,及时性较好.

    图  10  加速度计异常的检测结果
    Figure  10.  Detection results of accelerometer anomalies
    表  3  加速度计异常的检测位置与异常实际发生位置
    Table  3.  Detection position and actual occurrence position of accelerometer anomalies
    异常检测位置/s实际位置/s
    加速度计异常 152.252.2
    加速度计异常 252.152.1
    加速度计异常 340.540.5
    下载: 导出CSV 
    | 显示表格

    经过上述分析,基于残差生成器的异常检测方法可以用于磁悬浮列车悬浮系统的异常检测. 基于最低漏检率的异常检测方法能够实现磁悬浮列车悬浮系统3类典型异常的完全检测,不存在漏检的情况且不存在误检.

    本文针对电磁悬浮系统的异常检测问题开展了研究,以基于参数化残差的异常检测为核心,依据数据的均值和方差建立置信集,设计了基于最低异常漏检率的残差生成器,依据生成的残差数据在线进行异常检测. 最后,以长沙磁浮快线的运行数据为例,对悬浮系统常见的3类异常进行检测,实验结果表明,在给定异常误报率为5%时,本文所提出的方法能够实现对3类典型异常的全部检测,没有异常的漏检和误检. 检测时间对比结果进一步说明,本文提出的方法检测实时性较好.

  • 图 1  试件构造及尺寸

    Figure 1.  Structure and dimensions of specimens

    图 2  试验加载设计

    Figure 2.  Loading setup of tests

    图 3  键齿接缝试件破坏形态

    Figure 3.  Failure patterns of specimens with keyed joints

    图 4  平接缝试件破坏形态

    Figure 4.  Failure pattern of specimen with flat joint

    图 5  荷载-滑移曲线

    Figure 5.  Load-slip curves

    图 6  荷载-侧压应力关系

    Figure 6.  Relationships between loads and lateral compressive pressure

    图 7  键齿接缝受力示意

    Figure 7.  Force of keyed joint

    图 8  计算值与实测值对比

    Figure 8.  Comparison of calculated and measured values

    图 9  计算精度分析

    Figure 9.  Analysis of calculation accuracy

    表  1  试件参数

    Table  1.   Specimen parameters

    试件
    编号
    接缝
    类型
    设计侧压
    应力/MPa
    实际侧压
    应力/MPa
    预估承
    载力/kN
    F3-P平接缝33.1386.40
    F3-G干接缝33.05263.17
    F3-J胶接缝32.99526.40
    F6-J胶接缝66.05589.53
    F9-J胶接缝99.01646.66
    F12-J胶接缝1212.04698.57
    下载: 导出CSV

    表  2  试件抗剪承载力对比

    Table  2.   Comparison of shear bearing capacity of specimens

    试件编号 Vt/kN Vc/kN λ \overline \lambda Sλ λcov
    F3-G 623.7 657.1 1.05 1.09 0.11 0.10
    F3-J 756.3 714.7 0.94
    F6-J 836.8 880.3 1.05
    F9-J 927.6 1045.9 1.13
    F12-J 960.1 1211.5 1.26
    下载: 导出CSV

    表  3  参数分析结果

    Table  3.   Analysis results of influencing parameters

    影响参数 试件数量/个 \overline \lambda Sλ λcov
    键齿数量 单键齿 35 0.92 0.20 0.25
    多键齿 27 1.08 0.21 0.25
    接缝形式 干接缝 38 0.96 0.18 0.22
    胶接缝 24 1.03 0.25 0.31
    约束应力比 ≤1% 9 0.93 0.15 0.16
    (1%,2%] 13 0.89 0.24 0.27
    (2%,5%] 16 0.87 0.17 0.20
    (5%,10%] 15 1.09 0.19 0.21
    >10% 10 1.22 0.29 0.25
    下载: 导出CSV
  • [1] 《中国公路学报》编辑部. 中国桥梁工程学术研究综述•2021[J]. 中国公路学报,2021,34(2): 1-97. doi: 10.3969/j.issn.1001-7372.2021.02.002

    Editorial Department of China Journal of Highway and Transport. Review on China’s bridge engineering research: 2021[J]. China Journal of Highway and Transport, 2021, 34(2): 1-97. doi: 10.3969/j.issn.1001-7372.2021.02.002
    [2] 欧智菁,陈伟隆,曹磊. UHPC 预制管混凝土组合柱抗震性能[J]. 西南交通大学学报,2025,60(1):63-71.

    OU ZHIjing, CHEN Weilong, CAO Lei. Seismic performance investigation of UHPC precast concrete composite columns [J]. Journal of Southwest Jiaotong University,2025,60(1):63-71.
    [3] 邓宗才,王天宇. 高强钢丝布增强UHPC 约束混凝土柱轴压性能[J]. 西南交通大学学报,2025,60(1):72-83.

    DENG Zongcai, WANG Tianyu. Axial compression behavior of concrete columns confined by ultra-high performance concrete reinforced with high-strength steel wire cloth [J]. Journal of Southwest Jiaotong University,2025,60(1):72-83.
    [4] BUYUKOZTURK O, BAKHOUM M M, BEATTIE S M. Shear behavior of joints in precast concrete segmental bridges[J]. Journal of Structural Engineering, 1990, 116(12): 3380-3401. doi: 10.1061/(ASCE)0733-9445(1990)116:12(3380)
    [5] TURMO J, RAMOS G, APARICIO A C. Shear strength of dry joints of concrete panels with and without steel fibres. Application to precast segmental bridges[J]. Engineering Structures, 2006, 28(1): 23-33. doi: 10.1016/j.engstruct.2005.07.001
    [6] ROBERTS C L, BREEN J, KREGER M. Measurement based revisions for segmental bridge design and construction criteria. final report[R]. Austin: University of Texas, 1993.
    [7] GOPAL B A, HEJAZI F, HAFEZOLGHORANI M, et al. Shear strength of dry and epoxy joints for ultra-high-performance fiber-reinforced concrete[J]. ACI Structural Journal, 2020, 117(1): 279-288.
    [8] 孙雪帅. 预制拼装桥梁节段间接缝抗剪性能试验研究[D]. 南京:东南大学,2015.
    [9] 刘桐旭. 节段预制拼装UHPC梁接缝抗剪性能试验与理论研究[D]. 南京:东南大学,2017.
    [10] JIANG H B, HUANG C W, FENG J H, et al. Direct shear behavior of castellated dry RPC joints in precast concrete segmental bridges[J]. Structures, 2021, 33: 4579-4595. doi: 10.1016/j.istruc.2021.07.013
    [11] 王景全,胡玉庆,刘桐旭,等. 超高性能混凝土大键齿干接缝受剪性能与承载力计算方法[J]. 建筑结构学报,2021,42(10): 177-185.

    WANG Jingquan, HU Yuqing, LIU Tongxu, et al. Shear performance and bearing capacity calculation method for UHPC dry joints with large-keys[J]. Journal of Building Structures, 2021, 42(10): 177-185.
    [12] 潘仁胜,何伟伟,程灵霄,等. UHPC大键齿干接缝直剪性能及尺寸参数分析[J]. 湖南大学学报(自然科学版),2021,48(7): 129-137.

    PAN Rensheng, HE Weiwei, CHENG Lingxiao, et al. Direct shear behavior and dimensional parameter analysis of UHPC dry joint with big shear key[J]. Journal of Hunan University (Natural Sciences), 2021, 48(7): 129-137.
    [13] 闫泽宇. 节段预制拼装UHPC接缝抗剪性能研究[D]. 长沙:湖南大学,2019.
    [14] 郑凡,江建,李有志. 节段预制UHPC桥梁接缝抗剪性能试验研究[J]. 公路工程,2022,47(6): 1-9,41. doi: 10.12361/2661-3522-04-06-93072

    ZHENG Fan, JIANG Jian, LI Youzhi. Experimental study on the shear behavior of joints in precast segmental UHPC bridges[J]. Highway Engineering, 2022, 47(6): 1-9,41. doi: 10.12361/2661-3522-04-06-93072
    [15] 中华人民共和国住房和城乡建设部,国家市场监督管理总局. 混凝土物理力学性能试验方法标准:GB/T 50081—2019[S]. 北京:中国建筑工业出版社,2019.
    [16] 蒋大骅. 国外混凝土强度理论述评[J]. 同济大学学报,1978,6(1): 145-156.

    TaHua Chiang. A brief review on the failure criterion for concrete[J]. Journal of Tongji University, 1978, 6(1): 145-156.
    [17] 蒋大骅. 钢筋混凝土梁抗剪强度的一个新解[J]. 土木工程学报,1984,17(3): 23-34.

    JIANG Dahua. A new solution for shear strength of reinforced concrete beams[J]. China Civil Engineering Journal, 1984, 17(3): 23-34.
    [18] KIM Y J, CHIN W J, JEON S J. Interface shear strength at joints of ultra-high performance concrete structures[J]. International Journal of Concrete Structures and Materials, 2018, 12(6): 767-780. doi: 10.1186/s40069-018-0298-8
    [19] AASHTO. Guide specifications for design and construction of segmental concrete bridges: 2003 interim revisions [S]. Washington, D. C.: AASHTO, 2003.
    [20] ACI Committee. Building code requirements for reinforced concrete:ACI 318-08 [S]. Farmington Hills: [s.n.], 2008.
  • 期刊类型引用(5)

    1. 赵超,付斌,林立. 基于改进樽海鞘算法的含电动汽车微电网经济优化调度. 控制理论与应用. 2025(01): 167-180 . 百度学术
    2. 马龙飞,张宝群,王立永,曾佳妮,焦然,宫成. 基于数字孪生混合储能的电动汽车参与微电网负荷功率波动平抑研究. 汽车工程. 2024(06): 1045-1053+1084 . 百度学术
    3. 罗世刚,刘正英,李威武,许青,妥建军. 新能源接入交直流混合微电网的设备调度方法研究. 自动化仪表. 2024(06): 111-115 . 百度学术
    4. 张懿璞,金雨哲,柯吉,茹锋,张成朔. 考虑电动汽车负荷的高速公路服务区离网微电网容量配置. 长安大学学报(自然科学版). 2024(05): 126-140 . 百度学术
    5. 何玉朝,郭枫,钱辉敏,姚胜红,张能. 微电网源网荷储分布鲁棒优化调度方法研究. 电子设计工程. 2023(17): 123-127 . 百度学术

    其他类型引用(4)

  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  199
  • HTML全文浏览量:  93
  • PDF下载量:  29
  • 被引次数: 9
出版历程
  • 收稿日期:  2023-02-07
  • 修回日期:  2023-09-02
  • 网络出版日期:  2024-07-06

目录

/

返回文章
返回