Processing math: 100%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus
  • Indexed by Core Journals of China, Chinese S&T Journal Citation Reports
  • Chinese S&T Journal Citation Reports
  • Chinese Science Citation Database
LYU Biao, LIU Yumeng. Resilience Assessment Based on Bayesian Network for on-Board Subsystem of CTCS-3 Train Control System[J]. Journal of Southwest Jiaotong University, 2022, 57(5): 949-959. doi: 10.3969/j.issn.0258-2724.20210102
Citation: LIU Pengfei, LIU Hongjun, GAO Hao, REN Zunsong, ZHANG Kailong. Elastic Vibration of Wheelset and Its Dynamic Effect on Railway Heavy-Haul Freight Wagon[J]. Journal of Southwest Jiaotong University, 2022, 57(1): 90-98, 234. doi: 10.3969/j.issn.0258-2724.20210024

Elastic Vibration of Wheelset and Its Dynamic Effect on Railway Heavy-Haul Freight Wagon

doi: 10.3969/j.issn.0258-2724.20210024
  • Received Date: 11 Jan 2021
  • Accepted Date: 17 Nov 2021
  • Rev Recd Date: 04 May 2021
  • Available Online: 22 Nov 2021
  • Publish Date: 12 May 2021
  • To study the characteristics of wheelset elastic vibration and their effect on the dynamic performance of heavy-haul freight wagon, a railway wagon with a 30 t axle load is considered as the research object in this work. The running stability and the curving performance of vehicle are investigated in a comparative manner. First, the mathematical modeling method of elastic body based on multi-body dynamics is presented. Then, a finite element model of wheelset’s flexible body is established, and the elastic vibration modes are analyzed. The flexible wheelset is integrated in to the multi-rigid-body system, which then forms the rigid-flexible coupling dynamic model of heavy-haul wagon. Finally, for the multi-rigid body and rigid-flexible coupling modeling methods used for the wagon, the differences in wheelset vibration responses, running stability, and curving performance characteristics determined by these methods for a heavy-haul wagon were compared when using the main-line track irregularity with short-wave irregularity as the excitation source. The research results indicate that the deformation of elastic wheelset can alleviate the wheel-rail impact and also weaken the wheel-rail rigid restraining effect when compared with rigid wheelset. Therefore, the vibration amplitude for the flexible wheelset is reduced and the vehicle’s nonlinear critical speed is decreases by 9%. In sharp curves, the wheel-rail lateral force of flexible wheelset is reduced by 13.7%. The influence of wheelset’s elastic vibration on the dynamic performance of the heavy-haul wagon thus can not be ignored.

     

  • CTCS-3 (China train control system-3)级列控系统是保障时速300 km以上动车组安全、高效运行的神经中枢. 车载子系统是CTCS-3级列控核心组成部分[1]. 已有研究普遍使用可用性或可靠性作为车载子系统运行稳定性的测度指标. 文献[2]组合可靠性框图与马尔可夫随机过程方法计算了列控系统的可靠性、可用性及维修性指标. 文献[3]基于动态故障树方法研究了CTCS-3级列控车载设备的可靠性. 文献[4-5]运用贝叶斯网络(Bayesian network, BN)评估了CTCS-3级列控车载子系统的可靠性. 文献[6]以可靠度收益和维修成本为优化目标,提出了列控系统维修周期的确定方法. 文献[7]使用动态BN评估了CTCS-300T列控车载子系统的运行可靠性与可用性. 上述研究[2-7]在建模工具、考虑因素和侧重点上虽各有不同,但均遵从传统可靠性工程研究的基本思路,即不考虑特定扰动事件影响、以统计平均性能确定部件寿命分布和维修时间分布;同时,基于列控系统部件由电子器件构成的现实,假定部件寿命与维修时间均服从指数分布,即具有恒定的失效率和维修率. 传统方法用来评价列车运行环境稳定、无重大扰动发生时的系统性能是合适的. 但“7·23”事故表明,当面临雷电等非常态扰动时,传统可用性、可靠性指标不再适用. 这不仅因为面临非常态扰动时列控系统失效概率会大幅增加,而且失效概率与扰动类型直接相关,不能再视为常数. 同时,运营企业不仅关心列控系统能否抵御非常态扰动,也关心系统一旦失效能否快速有效恢复到可用状态. 因此,设计合理指标、专门针对非常态扰动下列控系统性能研究,有利于提升系统应对雷电等重大扰动事件的能力.

    不同于可用性指标,韧性通常以非常态扰动为背景,可同时描述系统抵御扰动和从扰动中快速恢复的能力. 韧性最初应用于生态学领域. 1973年,生态学家Holling[8]首次使用韧性评价生态系统稳定性,随后韧性被引入工程领域. 随着研究不断深入,交通运输领域的韧性问题也逐步引起研究者关注[9]. Murray-Tuite[10]首次使用韧性评价交通系统的性能. Henry等[11]将韧性定义为恢复期内系统性能恢复程度与扰动事件所造成的系统性能损失程度之比. Twumasi-Boakye等[12]认为有韧性的系统应是吸收扰动能力强、失效恢复速度快的系统. 文献[13]通过航空服务从破坏性事件中恢复的速度来衡量韧性. 文献[14]通过乘客延误等指标量化铁路运输系统韧性. 文献[15]提出了一种基于OD (origin-destination)数据集评估城市交通韧性的方法. 文献[16]通过对10个地区的智能交通系统进行随机和有针对性的破坏来研究网络效率和韧性,证明了韧性对智能交通系统的适用性. 文献[17]利用弹性三角模型对列控车载设备在不同故障情况下的恢复能力进行了评估. 综合来看,交通运输领域韧性研究目前处于起步阶段,仅有的少量研究主要针对航空网络、道路网络、城市公交网络等展开.

    基于上述分析,引入韧性作为评价非常态扰动下CTCS-3级列控车载子系统运行稳定性的指标. 首先,界定了韧性的内涵,定义了适用于列控车载子系统的韧性量化方法和5种基于韧性的部件重要度指标,构建了基于BN的列控车载子系统韧性评估模型,并开展算例研究.

    CTCS-3级列控车载子系统(简称车载子系统)是保证列车安全、高效运行的核心装置,其结构如图1所示[4-5].

    图  1  CTCS-3级列控车载子系统结构
    Figure  1.  On-board subsystem structure of CTCS-3

    车载子系统采用分布式结构,由安全计算机(vital computer, VC)、轨道电路信息接收单元(track circuit receiver, TCR)、列车接口单元(train interface unit, TIU)、无线传输单元(radio transmission unit, RTU)、GSM-R (global system for mobile communi- cations-railway)、无线电台(radio station, RS)、测速测距单元(speed and distance processing unit, SDU)、应答器传输模块(balise transmission module, BTM)及应答器天线(BTM ant)、人机界面(driver machine interface, DMI)、司法记录仪(juridical recorder unit, JRU)以及PROFIBUS (process fIeld bus)接口等部件组成. 为保证高可靠、高可用性,关键部件均采用冗余配置,例如,BTM、BTM ant、DMI等采用冷备结构,ATP-CU (automatic train protection-core unit)、C2-CU (China train control system-2 core unit)、速度传感器、雷达、TCR、GSM-R单元及其天线等采用热备结构[4-5].

    由于学科视角与对象不同,不同研究者对韧性内涵界定存在较大差异[18]. 部分研究将韧性表述为可靠性、脆弱性、恢复性等传统指标的简单延伸或扩展,与韧性内涵严重不匹配[18]. 文献[9]指出,韧性内涵应包括两方面:一是系统抵御扰动的能力,二是系统偏离稳定状态后的快速恢复能力. 对车载子系统而言,实际应用中遭遇雷电等非常态扰动时,可能导致部件损坏,由于部分部件采用冗余配置,单套设备损坏后系统仍能正常工作,即系统具有吸收或抵御扰动的能力;而非冗余配置部件的损坏则可能导致系统失效,但通过采取一定的技术手段(如替换失效部件)可尽快使系统恢复到正常状态,即系统具有从扰动中快速恢复的能力. 车载子系统的韧性是上述两种能力的综合体现.

    基于上述分析,将车载子系统或部件的韧性定义为:系统或部件受到非常态扰动时具有的维持正常运行或失效后快速恢复正常运行的能力.

    图2描述了非常态扰动事件下韧性与可用性指标的差异. 图中:t为时间;t0tdtr分别为扰动事件发生时刻、系统性能退化到最低程度时刻、系统性能恢复时刻;F(t)为系统性能函数;红色实线表示可用性指标,黑色实线表示韧性指标. 为便于分析,将扰动事件下系统性能划分为4种状态,分别为:1) 正常状态:0<t<t0,出现扰动前,系统保持稳定工作状态. 2) 降级状态:t0<t<td,扰动事件发生导致系统性能退化;3) 恢复状态:td<t<tr,通过维修受损部件将系统性能恢复到预期水平. 4) 新稳定状态:t>tr,完成维修后,系统性能恢复到新的稳定状态. 可以看出:扰动事件发生后,可用性指标与韧性指标出现明显分化,由于可用性不考虑特定事件影响,因此指标稳定;而韧性存在事件依赖,出现明显的降级然后恢复的过程.

    图  2  非常态扰动事件下韧性与可用性指标的差异
    Figure  2.  Differences in resilience and availability indexes under abnormal disturbance events

    车载子系统是由多部件构成的复杂系统,由于扰动事件通常导致部件失效进而引发系统失效或性能退化,因此,车载子系统韧性评估应以部件韧性评估为基础. 为便于描述,令Xei(t)表示扰动事件e发生后时刻t部件i的状态变量,具体定义如下:

    Xei(t)={0,扰动事件e发生后时刻t部件i不可用 1, 扰动事件 e发生后时刻t部件i可用(部件 i未失效或虽失效但已恢复).
    (1)

    ψei(t)表示扰动事件e发生后时刻t部件i的韧性值. 根据前文所述,韧性内涵为抵御扰动的能力和无法抵御扰动时的快速恢复能力. 根据这一界定,可将ψei(t)具体表示如下:

    ψei(t) = P(Xei(t)=1),
    (2)

    式中:P(·)为概率计算.

    从式(2)可以看出:扰动事件e发生后部件i在时刻t的韧性值等于该时刻部件处于可用状态的概率. 存在两种可能性,一是该部件能抵御扰动事件影响并未失效,二是该部件虽失效但在时刻t已恢复正常.

    ρei(t0)表示时刻t0发生扰动事件e时部件i失效的概率,φei(tt0)表示扰动事件 e发生后部件i在时刻t0未失效且在时段tt0仍正常工作的概率,ϕei(tt0)表示时刻t0发生扰动事件e后部件i失效但在时段tt0恢复正常的概率.

    根据式(1)、(2)可得

    ψei(t) = (1ρei(t0))φei(tt0)+ρei(t0)ϕei(tt0),t>t0.
    (3)

    由于列控系统部件由电子器件构成,工程计算中通常以指数分布描述部件寿命与维修时间分布. 因此,假定φei(tt0)ϕei(tt0)分别服从部件i的失效率λi和部件i的维修率μi 的指数分布. 在此条件下,φei(tt0)ϕei(tt0)可分别展开为

    φei(tt0) = eλi(tt0),t>t0,
    (4)
    ϕei(tt0) = 1eμi(tt0),t>t0.
    (5)

    根据式(4)、(5),可将式(3)重新描述如下:

    ψei(t) = (1ρei(t0))eλi(tt0)+ρei(t0)(1eμi(tt0)),t>t0.
    (6)

    式(6)中:(1ρei(t0))eλi(tt0)反映部件i抵御扰动的能力,ρei(t0)(1eμi(tt0))反映部件i无法抵御扰动失效后的恢复能力. 因此,式(6)集中反映了韧性内涵的两个方面. 同时可以看出:部件i的韧性与扰动事件紧密相关(不同类型扰动事件造成部件失效的概率一般不同,即ρei(t0)随事件e变化而取值可能不同),而可用性指标与具体扰动事件无关,即ρei(t0)不影响可用性指标,这一点在图2中可直观反映.

    显然,车载子系统的韧性取决于所有部件的韧性. 因此,可将车载子系统的韧性视为所有部件韧性的函数. 但由于车载子系统结构复杂,很难建立解析的函数表达式. 基于此,本文利用BN强大的推理能力,建立基于BN的车载子系统韧性评估模型.

    BN是一种用于表达和推理不确定信息的图形化概率模型[19],可表示为B = <G,P> = <<V,E>,P>. 其中:G = <V,E>为有向无环图,V为图中节点集合,集合中每一个元素代表一个随机变量;E为图中有向边集合,集合中元素表示变量之间的关系;概率参数P包括先验概率和条件概率表两部分. BN中,两节点间如存在一条有向边,则箭头连接的节点称为子节点,箭尾连接的节点称为父节点. 网络中,没有父节点的节点称为根节点,没有子节点的节点称为叶节点,其余节点称为中间节点. 先验概率表示根节点的边缘分布情况,而条件概率表用于表示有向边相连节点的依赖关系和依赖强度.

    BN具有强大推理能力,包括正向推理和反向推理. 其中,正向推理又称因果推理,可根据根节点状态推理非根节点状态;反向推理又称诊断推理,可根据叶节点状态反向推理非叶节点状态.

    当前,国内外有近10种BN建模软件,其中GeNIe软件以其简单直观、功能完善而被广泛使用. 本文使用GeNIe 2.0作为建模工具,考虑到韧性指标的动态时变性,使用图3所示的动态BN构建车载子系统韧性评估模型,图中:VC-C2、VC-C3分别为CTCS-2、CTCS-3级列控系统的安全计算机;VC-C21、VC-C22分别为VC-C2的冗余配置部件1、部件2;VC-C31、VC-C32分别为VC-C3的冗余配置部件1、部件2.

    图  3  基于GeNIe 2.0的车载子系统韧性评估模型
    Figure  3.  Resilience evaluation model of on-board subsystem based on GeNIe 2.0

    建模过程中,考虑到车载子系统的核心功能是安全防护,因此,对不影响行车安全的部件如司法记录仪JRU在建模过程中不予考虑[4-5]. 同时,为简化模型,将冷备结构当作单系处理,热备结构按并联处理[4-5].

    韧性评估的根本目的在于提升系统韧性,而要提升系统韧性,必须从提升部件韧性入手. 现实条件下,由于条件限制,不大可能对所有部件都采取改进措施,这就需要按照某种规则确定优先次序,即定义部件重要度指标确定部件重要度排序,选择排序靠前的部件采取改进措施. 参照可靠性工程领域[20]部件重要度定义思路,提出如下5个基于韧性的部件重要度指标:

    1) BIRNBAUM重要度(BIRNBAUM importance, BI)

    扰动事件e发生后时刻t部件i的BI为

    IBI,ei(t)=ψs,e(ψe(t))ψei(t)=ψs,e(ψe(t))|ψei(t)=1ψs,e(ψe(t))|ψei(t)=0,
    (7)

    式中:ψe(t) = (ψe1(t),ψe2,(t),,ψei(t),)为所有部件韧性值向量;ψs,e(ψe(t))为部件状态为ψe(t)时的系统韧性;ψs,e(ψe(t))|ψei(t)=1ψs,e(ψe(t))|ψei(t)=0 分别为当部件i可用和失效时的系统韧性.

    从式(7)可以看出:部件的BI反映了该部件可用与否对系统韧性的影响程度. 该指标的不足在于部件自身处于可用或失效的概率对重要度没有影响.

    2) 关键重要度(criticality importance, CI)

    扰动事件e发生后时刻t部件i的CI为

    ICI,ei(t)=IBI,ei(t)(1ψei(t))1ψs,e(ψe(t)) = (ψs,e(ψe(t))|ψei(t)=1ψs,e(ψe(t))|ψei(t)=0)(1ψei(t))1ψs,e(ψe(t)).
    (8)

    从式(8)可以看出:CI是以BI为基础构建的. CI考虑了部件自身失效的概率对系统韧性的潜在影响,可以弥补BI的不足.

    3) 改进潜力(improvement potential, IP)

    扰动事件e发生后时刻t部件i的改进潜力为

    IIP,ei(t)=ψs,e(ψe(t))|ψei(t)=1ψs,e(ψe(t)).
    (9)

    从式(9)可以看出:改进潜力的本质是确保部件处于可用状态对系统韧性的改进程度.

    4) 风险增加值(risk achievement worth, RAW)

    扰动事件e发生后时刻t部件i的风险增加值为

    IRAW,ei(t)=1ψs,e(ψe(t))|ψei(t)=01ψs,e(ψe(t)).
    (10)

    式(10)中,分母表示扰动事件e发生后时刻t系统处于不可用状态的概率,分子表示扰动事件e发生后时刻t部件i不可用时系统处于不可用状态的概率,其实质为部件i不可用对系统处于不可用状态的影响程度.

    5) 风险减少值(risk reduction worth, RRW)

    扰动事件e发生后时刻t部件i的风险减少值为

    IRRW,ei(t)=1ψs,e(ψe(t))1ψs,e(ψe(t))|ψei(t)=1.
    (11)

    式(11)中,分子表示扰动事件e发生后时刻t系统处于不可用状态的概率,分母表示扰动事件e发生后时刻t时如果部件i可用系统处于不可用状态的概率,其实质为部件i可用对减少系统处于不可用状态的贡献程度.

    我国地域辽阔,列车运行可能遭遇不同气象条件影响. 车载子系统由大量电子器件构成,对雷电、磁暴、冰雪天气等比较敏感. 以雷电干扰为例,2011年7月,雷电引起列控系统故障、轨道电路发码异常,导致甬温线特大铁路交通事故发生,据铁路部门统计,每10次雷电事故中就有3次是在雷电破坏铁路信号设备的情况下发生的,因此,雷电对铁路信号设备性能具有重要影响. 再以磁暴干扰为例,1989年发生超强磁暴引起俄罗斯高尔基铁路信号集中闭塞系统功能异常;2003年磁暴期间,瑞典Vladimir和Arzamasskaya铁路部分信号灯多次显示错误信号. 冰雪可能导致安装在车体表面的电子设备冻结进而发生故障,例如,2015年辽沈地区出现因大量降雪导致安装在车体表面的信号感应线圈、测速雷达等设备冻结,继而出现车载设备功能异常. 基于上述分析,选取雷电、磁暴和冰雪作为非常态扰动事件,开展车载子系统韧性评估.

    由于雷电、磁暴、冰雪等非常态扰动事件发生频度低,缺乏足够的历史数据,因此,遭受上述事件影响导致车载子系统部件失效的概率很难准确估计. 鉴于此,根据行业专家意见并结合不同事件对不同部件的潜在影响程度,将部件失效概率设置为区间取值. 将扰动事件对部件的影响程度划分为3类,即重大、重要和一般,3类影响程度下部件的失效概率取值范围依次为10−2 ~ 10−1、10−3 ~ 10−2和10−4 ~ 10−3. 表1描述了雷电、磁暴、冰雪3类扰动事件下部件失效概率ρei(t0)、部件的失效率λi和维修率μi的取值. 其中,部件失效概率ρei(t0)综合行业专家意见后确定,参数λiμi依据文献[5]确定. 从表1可以看出:雷电主要影响TCR、RS、BTM ant,磁暴主要影响TCR、GSM-R、RS、BTM、BTM ant,冰雪主要影响TCR、RS、BTM ant.

    表  1  不同扰动情景下部件参数取值
    Table  1.  Component parameter values under different disturbance scenarios
    编号名称雷电磁暴冰雪λi/h−1μi/h−1
    影响
    程度
    ρei(t0)影响
    程度
    ρei(t0)影响
    程度
    ρei(t0)
    E1/E2C2-CU1/C2-CU2重要10−3 ~ 10−2重要10−3 ~ 10−2一般10−4 ~ 10−31.20 × 10−52.0000
    E3/E4TCR1/TCR2重大10−2 ~ 10−1重大10−2 ~ 10−1重要10−3 ~ 10−22.30 × 10−62.0000
    E5/E6TIU1/TIU2重要10−3 ~ 10−2重要10−3 ~ 10−2一般10−4 ~ 10−32.10 × 10−52.0000
    E7/E8RTU1/RTU2重要10−3 ~ 10−2重要10−3 ~ 10−2一般10−4 ~ 10−31.80 × 10−52.0000
    E9/E10GSM-R1/GSM-R2重要10−3 ~ 10−2重大10−2 ~ 10−1一般10−4 ~ 10−31.45 × 10−82.0000
    E11/E12RS1/RS2重大10−2 ~ 10−1重大10−3 ~ 10−2重大10−2 ~ 10−11.20 × 10−52.0000
    E13/E14ATP-CU1/
    ATP-CU2
    重要10−3 ~ 10−2重要10−3 ~ 10−2一般10−4 ~ 10−31.49 × 10−52.0000
    E15/E16SDU1/SDU2重要10−3 ~ 10−2重要10−3 ~ 10−2一般10−4 ~ 10−32.50 × 10−90.2500
    E17/E18PROFIBUS1/
    PROFIBUS2
    重要10−3 ~ 10−2重要10−3 ~ 10−2一般10−4 ~ 10−36.00 × 10−60.0625
    E19BTM重要10−3 ~ 10−2重大10−2 ~ 10−1一般10−4 ~ 10−32.00 × 10−62.0000
    E20BTM ant重大10−2 ~ 10−1重大10−2 ~ 10−1重要10−3 ~ 10−27.00 × 10−80.2500
    E21DMI重要10−3 ~ 10−2重要10−3 ~ 10−2一般10−4 ~ 10−35.00 × 10−62.0000
    下载: 导出CSV 
    | 显示表格

    假定扰动发生在t = 5000 h,以扰动发生后2 h为测试期,每15 min取样一次. 基于保守思想,选取失效概率区间最大值作为扰动事件下部件失效概率. 根据前述韧性计算方法,将表1中参数取值代入式(6),依次计算得到3种扰动情景下不同时刻部件的韧性指标,将其设置为根节点先验概率,再利用BN正向推理,计算相应扰动情景及时间点非根节点的韧性指标.

    图4比较了韧性与可用性指标的差异. 可看出:在扰动发生后的较短时间内可用性指标几乎不发生变化,而韧性指标随扰动事件变化而变化. 例如,在扰动发生时刻(t = 5000 h),车载子系统可用性为0.999996,扰动发生后2 h (t = 5002 h),其值仍保持在0.999996 (由于数值变化很小,受计算精度限制,看不出数值变化);而伴随着扰动事件的发生,韧性指标发生了明显的变化. 在t = 5000 h时,面临雷击、磁暴、冰雪扰动时,系统韧性指标分别为0.8819、0.8017和0.9880;当t = 5002 h时,系统韧性指标相应变化为0.9389、0.9373和0.9939. 这是因为可用性指标反映的是统计意义上系统的平均性能,不受特定扰动事件影响;而韧性指标与扰动事件类型紧密相关,不同扰动事件下系统韧性指标不同.

    图  4  韧性指标与可用性指标的比较
    Figure  4.  Comparison of resilience index and availability index

    图5描述了不同扰动情景下部件与车载子系统的韧性. 其中,t = 5000 h时的韧性值反映系统抵御扰动的能力,t > 5000 h时的韧性值变化反映系统的恢复能力. 可以看出:韧性可全面描述部件与车载子系统抵御扰动和从扰动中恢复的能力;不同扰动情景下车载子系统韧性明显不同. 其中,磁暴影响最为显著,其次是雷电,冰雪影响最小. 这是因为磁暴、雷电会对多个部件造成重大、重要影响,而冰雪仅对少量部件造成重大、重要影响,对其余部件仅造成一般性影响. 而车载子系统的韧性是关于部件韧性的单调函数,因此出现图5所示的结果.

    图  5  扰动情景对韧性的影响
    Figure  5.  Effects of disturbance scenarios on resilience

    利用BN反向推理,可以计算车载子系统在扰动发生后某时刻失去韧性(不可用)时各部件在该时刻失去韧性(不可用)的概率. 通过BN反向推理可以找出造成车载子系统失去韧性的主要原因. 以雷击扰动为例,从图6可以看出:若车载子系统在扰动发生时(t = 5000 h)失去韧性,此时部件E20 (BTM ant)、E3/E4 (TCR)、E11/E12 (RSS)失去韧性的概率分别为0.84549、0.10075和0.10068,说明此时车载子系统不可用的最主要原因是部件E20 (BTM ant)不可用,其次是E3/E4 (TCR)和E11/E12 (RSS)不可用. 此外,需要注意的是,若扰动发生后2 h车载子系统仍不可用,此时部件E20 (BTM ant)、E3/E4 (TCR)、E11/E12 (RSS)失去韧性的概率分别变化为0.99256、0.00183和0.00184,几乎可以认定是由于部件E20 (BTM ant)不可用造成的. 出现这一结果的原因是E20 (BTM ant)的维修率明显低于部件E3/E4 (TCR)和E11/E12 (RSS)的维修率.

    图  6  BN反向推理
    Figure  6.  BN backward reasoning

    为识别不同扰动事件下不同部件的重要程度,按照前述重要度指标定义,计算得到的3种扰动情景下5个重要度指标排序结果如表2 ~ 4所示. 可以看出:同一扰动事件下,使用不同指标得到的部件重要度排序结果不完全一致. 例如,雷电干扰下,扰动发生时(t = 5000 h)部件E5的BI、CI、RAW排序为7,而IP、RRW排序为5. 这是因为不同重要度指标从不同角度反映部件的重要程度. 在实际应用中,应根据实际需要合理选择重要度指标.

    表  2  雷电扰动下部件重要度排序
    Table  2.  Component importance rankings under lightning disturbance
    排序t = 5000 ht = 5001 ht = 5002 h
    BICIIPRAWRRWBICIIPRAWRRWBICIIPRAWRRW
    1E7E7E7E7E7E7E7E1E7E1E7E7E1E7E1
    2E9E9E9E9E9E9E9E5E9E5E9E9E3E9E3
    3E13E13E13E11E13E13E13E7E11E7E11E13E5E11E5
    4E11E1E1E13E1E11E1E9E13E9E13E1E7E13E7
    5E1E11E5E1E5E1E11E13E1E13E1E11E9E1E9
    6E3E3E11E3E11E3E3E3E3E3E3E3E11E3E11
    7E5E5E3E5E3E5E5E11E5E11E5E5E13E5E13
    8E15E15E15E15E15E15E19E15E15E15E15E19E15E15E15
    9E17E17E17E17E17E17E21E17E17E17E17E21E17E17E17
    10E19E19E19E19E19E21E15E21E19E21E19E15E19E19E19
    11E21E21E21E20E21E19E17E19E20E19E21E17E21E20E21
    12E20E20E20E21E20E20E20E20E21E20E20E20E20E21E20
    下载: 导出CSV 
    | 显示表格
    表  3  磁暴扰动下部件重要度排序
    Table  3.  Component importance rankings under magnetic storm disturbance
    排序t = 5000 ht = 5001 ht = 5002 h
    BICIIPRAWRRWBICIIPRAWRRWBICIIPRAWRRW
    1E7E7E7E7E7E7E7E1E7E1E7E7E1E7E1
    2E13E13E13E9E13E13E13E7E9E7E9E13E3E9E3
    3E9E1E1E11E1E9E1E13E11E13E11E1E5E11E5
    4E11E9E5E13E5E11E9E5E13E5E13E9E7E13E7
    5E1E11E9E1E9E1E11E9E1E9E1E11E9E1E9
    6E3E3E11E3E11E3E3E11E3E11E3E3E11E3E11
    7E5E5E3E5E3E5E5E3E5E3E5E5E13E5E13
    8E15E15E15E15E15E1512E15E15E15E15E21E15E15E15
    9E17E17E17E17E17E17E15E17E17E17E17E19E17E17E17
    10E21E21E21E19E21E21E17E21E19E21E21E15E21E19E21
    11E19E19E19E20E19E19E19E19E20E19E19E17E19E20E19
    12E20E20E20E21E20E20E20E20E21E20E20E20E20E21E20
    下载: 导出CSV 
    | 显示表格

    面临的扰动类型不同,部件重要度排序也存在差异. 例如,雷电、磁暴扰动发生时(t = 5000 h)部件E9 (GSM-R)的IP排序分别为2和5. 这一结论有很强的现实指导意义. 我国地域辽阔,不同地域列车面临的主要扰动事件类型不同,如南方容易遭受雷电干扰,而北方容易遭受冰雪灾害. 不同区域的铁路运营企业应根据面临的主要扰动事件类型,有针对性地确定车载子系统中的重要部件,加强对重要部件的检修和防护,以提高系统韧性.

    除此之外,部件的重要度排序可能随时间动态变化,并且此变化不具有单调性. 例如,磁暴干扰下,当t = 5000 h时,部件E13的RAW排序为2,而当t = 5001 h和t = 5002 h时,E13的RAW排序分别变为3和7.

    表  4  冰雪扰动下部件重要度排序
    Table  4.  Component importance rankings under snow and ice disturbances
    排序t = 5000 ht = 5001 ht = 5002 h
    BICIIPRAWRRWBICIIPRAWRRWBICIIPRAWRRW
    1E7E7E7E7E7E7E7E1E7E1E7E7E1E7E1
    2E9E9E9E9E9E9E9E3E9E3E9E9E3E9E3
    3E13E13E13E11E13E11E13E5E11E5E11E13E5E11E5
    4E11E11E1E13E1E13E11E7E13E7E13E11E7E13E7
    5E1E1E5E1E5E1E1E9E1E9E1E1E9E1E9
    6E3E3E3E3E3E3E3E11E3E11E3E3E11E3E11
    7E5E5E11E5E11E5E5E13E5E13E5E5E13E5E13
    8E15E15E15E15E15E15E19E15E15E15E15E19E15E15E15
    9E17E17E17E17E17E17E20E17E17E17E17E20E17E17E17
    10E19E19E19E19E19E19E15E19E19E19E19E15E19E19E19
    11E20E20E20E20E20E20E17E20E20E20E20E17E20E20E20
    12E20E20E20E20E20E20E20E20E20E20E20E20E20E20E20
    下载: 导出CSV 
    | 显示表格

    由于表1中设定的不同扰动情景下部件失效概率不是一个确定值,而是一定范围的数值区间. 因此,有必要通过灵敏度分析,验证部件失效概率的取值变化对系统韧性的潜在影响. 假定部件失效概率在取值范围内服从均匀分布,依据模特卡罗仿真思想,随机抽样2000次得到车载子系统韧性的频率直方图如图7所示. 可以看出:雷电、磁暴、冰雪扰动发生时(t = 5000 h),车载子系统的韧性取值范围分别为 [0.8800,0.9900]、[0.8400,0.9900] 和 [0.9900, 1.0000]. 这一结果反映了不同扰动事件对车载子系统韧性的潜在影响程度.

    图  7  不同扰动下车载子系统韧性的频率直方图
    Figure  7.  Frequency histogram of on-board subsystem resilience under different disturbances

    引入韧性作为非常态事件下车载子系统运行稳定性测度指标,构建了车载子系统韧性量化评估方法和基于贝叶斯网络的韧性评估模型,并定义了5种基于韧性的部件重要度指标. 算例结果表明:

    1) 韧性可全面描述车载子系统抵御扰动和从扰动中恢复的能力,并与可用性指标存在明显差异. 由于不考虑特定扰动事件影响,雷电等非常态事件发生时可用性指标稳定,而韧性与扰动事件紧密相关,且随扰动事件不同而取值不同.

    2) 不同扰动情境下车载子系统韧性明显不同. 其中,磁暴影响最为显著,其次是雷电,冰雪影响最小.

    3) 部件重要度与扰动情景相关,同一部件在不同扰动情景下重要度排序可能不同;同时,部件重要度随时间动态变化,扰动发生后的不同时间点,同一部件同一指标的重要度排序可能不同.

  • [1]
    钟硕乔. 车辆轨道耦合系统中的旋转柔性轮对建模研究[D]. 成都: 西南交通大学, 2017.
    [2]
    PENG B, IWNICKI S, SHACKLETON P, et al. The influence of wheelset flexibility on polygonal wear of locomotive wheels[J]. Wear, 2019, 432/433: 102917.1-102917.11.
    [3]
    BAEZA L, GINER-NAVARRO J, THOMPSON D J, et al. Eulerian models of the rotating flexible wheelset for high frequency railway dynamics[J]. Journal of Sound and Vibration, 2019, 449: 300-314. doi: 10.1016/j.jsv.2019.03.002
    [4]
    崔潇,姚建伟,胡晓依,等. 欧拉坐标系下柔性轮对旋转效应对轮轨力的影响[J]. 中国铁道科学,2019,40(4): 120-128. doi: 10.3969/j.issn.1001-4632.2019.04.15

    CUI Xiao, YAO Jianwei, HU Xiaoyi, et al. Rotation effect of flexible wheelset on wheel-rail force in Euler coordinate system[J]. China Railway Science, 2019, 40(4): 120-128. doi: 10.3969/j.issn.1001-4632.2019.04.15
    [5]
    FOURIE D, FRÖHLING R, HEYNS S. Railhead corrugation resulting from mode-coupling instability in the presence of veering modes[J]. Tribology International, 2020, 152: 106499. doi: 10.1016/j.triboint.2020.106499
    [6]
    CRUCEANU I C, SOROHAN Ș. Determination of the harmonic response of a railway wheelset using the finite element analysis method[J]. Procedia Manufacturing, 2020, 46: 173-179. doi: 10.1016/j.promfg.2020.03.026
    [7]
    李国芳,岳鹏,丁旺才,等. 轮对柔性对车辆动态曲线通过性能的影响研究[J]. 铁道标准设计,2019,63(9): 173-179.

    LI Guofang, YUE Peng, DING Wangcai, et al. Research on the influence of flexibility wheelset on vehicle dynamic curve passing performance[J]. Railway Standard Design, 2019, 63(9): 173-179.
    [8]
    张波,罗光兵,蒋忠城,等. 柔性结构对车辆运动稳定性的影响[J]. 技术与市场,2019,26(6): 5-8.

    ZHANG Bo, LUO Guangbing, JIANG Zhongcheng, et al. The influence of flexible structure on vehicle dynamic stability[J]. Technology and Market, 2019, 26(6): 5-8.
    [9]
    杨云帆,周青,巩磊,等. 轮对柔性对直线电机车辆动态响应的影响分析[J]. 西南交通大学学报,2020,55(6): 1313-1319. doi: 10.3969/j.issn.0258-2724.20180866

    YANG Yunfan, ZHOU Qing, GONG Lei, et al. Influence of wheelset flexibility on dynamic response of linear induction motor vehicles[J]. Journal of Southwest Jiaotong University, 2020, 55(6): 1313-1319. doi: 10.3969/j.issn.0258-2724.20180866
    [10]
    金学松,吴越,梁树林,等. 车轮非圆化磨耗问题研究进展[J]. 西南交通大学学报,2018,53(1): 1-14. doi: 10.3969/j.issn.0258-2724.2018.01.001

    JIN Xuesong, WU Yue, LIANG Shulin, et al. Mechanisms and countermeasures of out-of-roundness wear on railway vehicle wheels[J]. Journal of Southwest Jiaotong University, 2018, 53(1): 1-14. doi: 10.3969/j.issn.0258-2724.2018.01.001
    [11]
    杨润芝,曾京. 高阶车轮多边形对轮轨系统振动影响分析[J]. 振动与冲击,2020,39(21): 101-110.

    YANG Runzhi, ZENG Jing. Influences of higher order wheel polygon on vibration of wheel-rail system[J]. Journal of Vibration and Shock, 2020, 39(21): 101-110.
    [12]
    王相平,王红兵,贾文慧,等. 车轮扁疤对高速车辆动态曲线通过性能的影响[J]. 铁道科学与工程学报,2020,17(9): 2198-2207.

    WANG Xiangping, WANG Hongbing, JIA Wenhui, et al. The influence of wheel flat on the performance of high-speed vehicle dynamic on railway curve negotiation[J]. Journal of Railway Science and Engineering, 2020, 17(9): 2198-2207.
    [13]
    刘国云,曾京,邬平波,等. 车轮扁疤所引起的车辆系统振动特性分析[J]. 机械工程学报,2020,56(8): 182-189. doi: 10.3901/JME.2020.08.182

    LIU Guoyun, ZENG Jing, WU Pingbo, et al. Vibration characteristic analysis of vehicle systems due to wheel flat[J]. Journal of Mechanical Engineering, 2020, 56(8): 182-189. doi: 10.3901/JME.2020.08.182
    [14]
    杨广雪,赵方伟,李秋泽,等. 高速列车轮轨接触几何参数对轮轨磨耗的影响研究[J]. 铁道学报,2019,41(2): 50-56.

    YANG Guangxue, ZHAO Fangwei, LI Qiuze, et al. Study of influences of high-speed train wheel-rail contact geometric parameters on wheel-rail wear[J]. Journal of the China Railway Society, 2019, 41(2): 50-56.
    [15]
    魏世明. 30 t轴重货车HFZ915铸钢车轮的研发[J]. 电力机车与城轨车辆,2016,39(1): 58-60.

    WEI Shiming. Research and development of HFZ915 cast steel wheel for 30 t axle-load freight car[J]. Electric Locomotives & Mass Transit Vehicles, 2016, 39(1): 58-60.
    [16]
    吴毅,刘鑫贵,项彬,等. 重载铁路货车用LZ45CrV车轴钢综合性能的试验研究[J]. 中国铁道科学,2015,36(2): 68-72. doi: 10.3969/j.issn.1001-4632.2015.02.10

    WU Yi, LIU Xingui, XIANG Bin, et al. Experimental study on the comprehensive performance of LZ45CrV axle steel for heavy haul freight car[J]. China Railway Science, 2015, 36(2): 68-72. doi: 10.3969/j.issn.1001-4632.2015.02.10
    [17]
    国家铁路局. 铁道车辆轮对及轴承型式与基本尺寸: TB/T 1010—2016 [S]. 北京: 中国铁道出版社, 2016.
    [18]
    翟婉明. 车辆-轨道耦合动力学[M]. 4版. 北京: 科学出版社, 2015.
    [19]
    刘鹏飞,王开云,翟婉明. 长大列车通过弹性曲线轨道仿真求解方法研究[J]. 西南交通大学学报,2018,53(1): 31-37. doi: 10.3969/j.issn.0258-2724.2018.01.004

    LIU Pengfei, WANG Kaiyun, ZHAI Wanming. Simulation solution method for long and heavy-haul trains negotiating elastic curved tracks[J]. Journal of Southwest Jiaotong University, 2018, 53(1): 31-37. doi: 10.3969/j.issn.0258-2724.2018.01.004
    [20]
    李亨利,李芾,傅茂海,等. 斜楔摩擦减振器建模及动力学分析[J]. 铁道科学与工程学报,2015,12(5): 1191-1199. doi: 10.3969/j.issn.1672-7029.2015.05.031

    LI Hengli, LI Fu, FU Maohai, et al. Dynamic modeling and simulation of wedges friction damper[J]. Journal of Railway Science and Engineering, 2015, 12(5): 1191-1199. doi: 10.3969/j.issn.1672-7029.2015.05.031
    [21]
    王勇,曾京,张卫华. 铁道货车非线性稳定性[J]. 交通运输工程学报,2002,2(2): 36-40. doi: 10.3321/j.issn:1671-1637.2002.02.009

    WANG Yong, ZENG Jing, ZHANG Weihua. Nonlinear stability of railway freight cars[J]. Journal of Traffic and Transportation Engineering, 2002, 2(2): 36-40. doi: 10.3321/j.issn:1671-1637.2002.02.009
    [22]
    国家铁路局. 铁路线路设计规范: TB 10098—2017[S].北京: 中国铁道出版社, 2017.
  • Relative Articles

    [1]YANG Chunlei, WANG Kaiyun. Simulation Calculation Method and Application for Relative Displacement of Heavy Hall Freight Suspension on Curved Track[J]. Journal of Southwest Jiaotong University, 2022, 57(5): 1008-1016. doi: 10.3969/j.issn.0258-2724.20210158
    [2]YU Zhixiang, HU Guanghua, LI Tongmei, FENG Dubei. Analysis of Vibration Reduction and Vibration Measurement for Long-Span Railway Station Floor Slab[J]. Journal of Southwest Jiaotong University, 2019, 54(2): 296-303, 342. doi: 10.3969/j.issn.0258-2724.20170746
    [3]LI Jincheng, DING Junjun, WU Pengpeng, YANG Yang, LI Fu. Analysis of Dynamics and Wheel Wear of Low Floor Vehicle Based on Different Patterns[J]. Journal of Southwest Jiaotong University, 2019, 54(1): 14-22. doi: 10.3969/j.issn.0258-2724.20170757
    [4]YANG Chunlei, HUANG Yunhua, LI Fu. Comparative Study on the Curving Performance of Chinese Railway Heavy Haul Freight Bogies[J]. Journal of Southwest Jiaotong University, 2019, 54(3): 619-625. doi: 10.3969/j.issn.0258-2724.20180178
    [5]WANG Boming, CHEN Hongyu, LI Zhize, TAN Hongyuan. Curve-Passing Dynamic Performance of Portal Bogie with Forced Steering Mechanism[J]. Journal of Southwest Jiaotong University, 2016, 29(1): 91-97. doi: 10.3969/j.issn.0258-2724.2016.01.014
    [6]YANG Chun-Lei, LI Fei, HUANG Yun-Hua. Optimization ofP rimary Vertical Suspension of Heavy Haul Freight Car[J]. Journal of Southwest Jiaotong University, 2011, 24(5): 820-825.
    [7]ZHOU Ning, ZHANG Weihua. Dynamic Performances of Pantograph-Catenary System with Double Pantographs[J]. Journal of Southwest Jiaotong University, 2009, 22(4): 552-557.
    [8]MA Weihua, LUO Shihui, SONG Rongrong. Improvement of Dynamic Performance of Speed-Raising Locomotive with Motor Mounted on Frame[J]. Journal of Southwest Jiaotong University, 2007, 20(1): 84-88.
    [9]WANG Yong-guan, CHENKang. Effects ofCrosswinds on Curve Negotiation ofH igh-Speed Power Cars[J]. Journal of Southwest Jiaotong University, 2005, 18(2): 224-227.
    [10]WANG Yong, L Ke-wei, GUO Xiao-feng, CHEN Bao-shun. Dynamic Behaviors and Parameters of Freight Cars with Frame-Bogies[J]. Journal of Southwest Jiaotong University, 2003, 16(2): 127-131.
    [11]WANG Kai-yun, ZHAI Wan-ming, CAI Cheng-biao. Simulation of Lateral Dynamic Performances of Locomotives and Vehicles—Comparison of Vehicle/Track Coupling Model with Traditional Vehicle Model[J]. Journal of Southwest Jiaotong University, 2003, 16(1): 17-21.
    [12]LIUHong-you, ZENG Jing, LI Wen-xue. Dynamic Performances of Meter-Gauge Passenger Car with Forced-Steering Bogies[J]. Journal of Southwest Jiaotong University, 2002, 15(5): 579-583.
    [13]CHI Mao-ru, WANG Kai-wen, ZHANG Wei-hua. Dynamic Performances of the Bogie with Magnetic Fluid Coupled Wheelsets[J]. Journal of Southwest Jiaotong University, 2002, 15(6): 664-668.
    [14]WANG Ping, LIUXue-yi, LI Cheng-hui. Dynamic Analysis ofTrack Wave Propagation Supported on Continuous Elastic Foundations[J]. Journal of Southwest Jiaotong University, 2001, 14(1): 17-22.
    [15]ZENGJing, WANG Yong. Nonlinear Dynamic Analysis for Railway Freight Cars[J]. Journal of Southwest Jiaotong University, 2000, 13(4): 399-403.
  • Cited by

    Periodical cited type(7)

    1. 王鹏,孟鑫,赵军,刘宇,杨翠萍. 基于Faster R-CNN的铁路货车心盘螺栓智能检修. 自动化与仪器仪表. 2024(05): 237-240+245 .
    2. 贾彬,王志伟,张铭涛,刘鹏飞,刘国云,张卫华. 考虑闸瓦制动的重载列车动力学建模与车辆曲线制动安全分析. 重庆理工大学学报(自然科学). 2024(06): 155-163 .
    3. 贾小平,戴焕云,朱程,魏来. 考虑柔性轮对的地铁车辆系统振动特性研究. 铁道车辆. 2024(06): 17-24+31 .
    4. 孙小娟,孙洪亮,李幸人,王宇. 多边形车轮对六轴重载货运机车动力学性能影响分析. 机械强度. 2023(03): 668-676 .
    5. 刘嘉兴,吕大立,吕可维,张琪昌,李玉龙. 基于BP神经网络的转向架斜楔参数优化. 铁道学报. 2023(10): 52-59 .
    6. 朱红霞,刘红军,张成成,王鹏,王翠艳. 不同轨道激扰下柔性轮对重载货车垂向振动特性分析. 机械设计与研究. 2022(04): 87-93 .
    7. 和风,冯洋,佟来生,舒瑶,赵春发. 高速磁浮车辆通过平面曲线时悬浮架和电磁铁的弹性变形分析. 机车电传动. 2022(04): 1-8 .

    Other cited types(15)

  • Created with Highcharts 5.0.7Amount of accessChart context menuAbstract Views, HTML Views, PDF Downloads StatisticsAbstract ViewsHTML ViewsPDF Downloads2024-062024-072024-082024-092024-102024-112024-122025-012025-022025-032025-042025-050510152025
    Created with Highcharts 5.0.7Chart context menuAccess Class DistributionFULLTEXT: 43.4 %FULLTEXT: 43.4 %META: 53.2 %META: 53.2 %PDF: 3.4 %PDF: 3.4 %FULLTEXTMETAPDF
    Created with Highcharts 5.0.7Chart context menuAccess Area Distribution其他: 6.4 %其他: 6.4 %其他: 0.4 %其他: 0.4 %China: 0.6 %China: 0.6 %Germany: 0.3 %Germany: 0.3 %Seattle: 0.3 %Seattle: 0.3 %上海: 0.7 %上海: 0.7 %东莞: 0.3 %东莞: 0.3 %临汾: 0.3 %临汾: 0.3 %乌鲁木齐: 0.1 %乌鲁木齐: 0.1 %伊利诺伊州: 0.1 %伊利诺伊州: 0.1 %保定: 0.3 %保定: 0.3 %信阳: 0.1 %信阳: 0.1 %兰州: 0.1 %兰州: 0.1 %内江: 0.1 %内江: 0.1 %北京: 2.7 %北京: 2.7 %十堰: 0.4 %十堰: 0.4 %南京: 0.9 %南京: 0.9 %南昌: 0.1 %南昌: 0.1 %南通: 0.3 %南通: 0.3 %台州: 0.1 %台州: 0.1 %合肥: 0.1 %合肥: 0.1 %哈尔滨: 0.3 %哈尔滨: 0.3 %哥伦布: 0.6 %哥伦布: 0.6 %嘉兴: 0.3 %嘉兴: 0.3 %大连: 0.3 %大连: 0.3 %天津: 0.7 %天津: 0.7 %安康: 0.1 %安康: 0.1 %宣城: 0.4 %宣城: 0.4 %岳阳: 0.1 %岳阳: 0.1 %常州: 0.3 %常州: 0.3 %广州: 0.3 %广州: 0.3 %张家口: 1.2 %张家口: 1.2 %成都: 1.6 %成都: 1.6 %扬州: 0.1 %扬州: 0.1 %昆明: 0.6 %昆明: 0.6 %杭州: 0.9 %杭州: 0.9 %池州: 0.3 %池州: 0.3 %洛阳: 0.6 %洛阳: 0.6 %淄博: 0.1 %淄博: 0.1 %深圳: 0.1 %深圳: 0.1 %温州: 0.4 %温州: 0.4 %漯河: 1.3 %漯河: 1.3 %盐城: 0.1 %盐城: 0.1 %石家庄: 4.3 %石家庄: 4.3 %秦皇岛: 0.1 %秦皇岛: 0.1 %芒廷维尤: 25.1 %芒廷维尤: 25.1 %芝加哥: 0.1 %芝加哥: 0.1 %苏州: 0.1 %苏州: 0.1 %西宁: 38.6 %西宁: 38.6 %诺沃克: 0.1 %诺沃克: 0.1 %贵阳: 0.3 %贵阳: 0.3 %运城: 0.9 %运城: 0.9 %重庆: 0.1 %重庆: 0.1 %长春: 0.3 %长春: 0.3 %长沙: 2.4 %长沙: 2.4 %青岛: 0.7 %青岛: 0.7 %马鞍山: 0.3 %马鞍山: 0.3 %其他其他ChinaGermanySeattle上海东莞临汾乌鲁木齐伊利诺伊州保定信阳兰州内江北京十堰南京南昌南通台州合肥哈尔滨哥伦布嘉兴大连天津安康宣城岳阳常州广州张家口成都扬州昆明杭州池州洛阳淄博深圳温州漯河盐城石家庄秦皇岛芒廷维尤芝加哥苏州西宁诺沃克贵阳运城重庆长春长沙青岛马鞍山

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(18)

    Article views(425) PDF downloads(24) Cited by(22)
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return