基金项目:西北民族大学中央高校基本科研业务费资金资助项目(31920160065)
第一作者:卢永飞(1982-),男,博士生,讲师,主要从事桥梁结构动力学方面的教学与研究工作.E-mail: ban43954069@163.com
(1.西北民族大学 土木工程学院,甘肃 兰州,730030; 2.兰州交通大学 土木工程学院,甘肃 兰州,730070)
(1.School of Civil Engineering, Northwest Minzu University, Lanzhou 730030, China; 2.School of Civil Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China)
safety engineering; train; running behavior evaluation; set pair analysis; entropy; MATLAB
DOI: 10.15986/j.1006-7930.2020.03.014
针对过桥列车走行性评价中各评价指标数据测定的不确定性和随机性,在列车走行性综合评价中应用集对分析理论,首先分析测点实测数据与列车走行性评价指标之间的联系度,再根据联系度对各实测数据进行初步分类,然后进一步对初步分类的测点实测数据做综合联系度计算,判定列车走行性等级,并将信息论中的熵值理论引入评价模型确定各指标的权重,使权重的分配有了一定的理论依据.最后通过对一在不同车速下列车过桥的6组样本数据进行分析,用MATLAB编程计算得到,与传统等权重指标评价方法相比,基于熵权的集对分析模型概念清晰、准确度高、符合实际.
The theory of set pair analysis is applied to the comprehensive evaluation of train running performance in view of the uncertainty and randomness of data measurement of evaluation indexes in the evaluation of train running performance of bridge crossing.Firstly, the connection degree between the measured data of the measuring point and the evaluation index of train running was analyzed and then a preliminary classification, and a comprehensive connection calculation on the measured data of the preliminary classification of the measuring point,were also conducted in determining the train running level. The entropy value theory in information theory is introduced into the evaluation model to determine the weight of each index, so that the weight distribution has a certain theoretical basis.Finally, 6 sets of sample data of a vehicle passing through the bridge at different speeds were analyzed. Compared with the traditional equal weight index evaluation method, the evaluation result of the set pair analysis model based on entropy weight is more consistent with the reality, clearer and more accurate through the programming calculation of MATLAB.
随着列车行驶速度的不断提高和列车载荷的不断加大,列车通过桥梁的安全性和平稳性问题越来越得到人们的重视.然而由于列车过桥时列车-轨道-桥梁的动力相互作用是一个非常复杂的过程,许多相关因素都具有一定的随机性和模糊性,以致使实测数据不可避免的产生随机误差,对列车运行安全性和平稳性最终评价结果产生一定的影响[1-4].针对该情况采用模糊诊断方法评估该动力相互作用系统的安全性和平稳性是行之有效的.2000年,翟婉明等[5]首次从车辆-轨道相互作用整体系统的角度,运用车辆-轨道耦合动力学理论,对列车走行性受轨道结构刚度的影响规律进行了研究.2002年,曾庆元院士[6]第一次提出了列车脱轨能量随机分析理论,用此理论预报了一货物列车脱轨实例,预报结果与实际脱轨事故分析结果基本一致,并计算了该试验列车最大振动响应,计算结果与测试结果良好接近.在文献[6]的基础上,2002年,向俊等[7]进一步论述了列车脱轨能量随机分析的理论与实践,再次作了实例脱轨分析,提出了预防列车脱轨的措施,得到车轨系统抗脱轨安全系数的计算方法.2011年,周智辉等[8]人运用列车脱轨能量随机分析理论分析了金温线武义江桥列车走行性,并基于该理论,对铁路桥梁行车安全指标分析方法进行了研究.2015年,梁世庆等[9]把物元可拓数学方法、熵权理论和关联度函数运用到轨道交通建设中,建立了轨道交通建设时序决策的熵权物元可拓模型.2018年,刘清漪等[10]将交通、经济、社会和环境作为轨道交通的影响子系统,设计量化模型计算了外部效应对轨道交通线路的影响程度,并对深圳地铁3号线建立灰色聚类模型进行了实例分析.KUMAR K和GORTMAKER S L等人已经将集对分析理论已经广泛应用于多个领域[11-15].借助于指标数据测定的确定性和不确定性分析,本文运用集对分析理论,先对所测列车运行指标进行定性分析,通过计算联系度对列车过桥运行安全性和平稳性做出定量分析,在确定所测列车运行指标权重时,运用数学中的信息熵反映实测数据的效用值,对列车过桥运行安全性和平稳性做出全面、客观的评价,最后通过一实例应用对铁路桥上列车的走行安全性和平稳性做出评估.
集对概念由我国学者赵克勤提出,用来处理不确定性问题的一种数学分析方法.集对分析理论的核心思想:在具体问题背景下,对集合A、B组成的集对进行分析,两组集对得到N个特性,其中有S个为集对A、B所共同具有的元素数,有P个 为A、B中相互对立的元素数,在其余F个特性关系为既不共有也不对立[16],则两个集合的联系度
μ=S/N+F/Ni+P/Nj=a+bi+cj,(1)
其中,(a+b+c=1).式中:μ为集合A、B的联系度; S/N、F/N、P/N分别为集合A与集合B的同一度、差异度、对立度,分别简记为a、b、c; i为差异标记,在[-1,1]区间视不同情况取值; j为对立系数,一般取为-1,也可仅起标记作用.联系度μ与不确定系数i是本理论的核心,该理论对常见的一系列不确定现象均能够进行描述.
基于集对分析的列车过桥运行安全性和平稳性评价是将列车过桥走行性实测数据与规范或标准中的列车过桥走行性评价指标构成一个集对,分析构成的两个集对,就可以得到安全性和平稳性评价的量化指标.列车过桥运行安全性和平稳性评价常规评价指标有7个,分别为脱轨系数、轮重减载率、轮轨横向力、车体垂向加速度、车体横向加速度、Sperling舒适性垂向指标、Sperling舒适性横向指标,设一组实测样本中优于标准的有S个评价指标,劣于标准的有P个评价指标,未测或缺乏比较的有其余的F个评价指标构成,则各数据样本的联系度可通过式(1)计算出,公式中a、b、c三个数据的大小关系可初步反映该组样本数据所得评价指标的优劣.
列车过桥安全性和平稳性评价指标特性属于越小越优型,其联系度为
μst={1+0i+0j,x∈[0,S1]
(S2-x)/(S2-S1)+(x-S1)/(S2-S1)i+0j, x∈[S1, S2]
0+(S3-x)/(S3-S2)i+(x-S2)/(S3-S2)j, x∈[S2, S3]
0+0i+1j,x∈[S3,+SymboleB@](2)
式中S1、S2、S3分别为标准的界限值; t为第t项指标; s为第s个待测样本; x为测点s的第t项评价指标的实测值.
在目前已有的列车过桥运行安全性和平稳性评价中,各项评价指标均取等权重,没有考虑到不同指标的相对重要性,不能符合实际列车运行情况,使得走行性评价结果不够准确.指标权重的合理性与安全评价的准确度息息相关,在基于熵权-集对分析的列车过桥运行安全性和平稳性评价模型中,这里采用熵权法来确定各指标的权重,计算步骤如下[17].
(1)先构建m个评价对象,每个评价对象设有n个评价指标,构成判断矩阵R
R=(rst)m×n,(s=1,2,…,m; t=1,2,…,n)(3)
(2)R归一化处理,得到矩阵B.
越大越优型指标为
bst=(rst-rmin)/(rmax-rmin)
越小越优型指标为
bst=(rmax-rst)/(rmax-rmin)(4)
(3)各评价指标的熵为
Ht=-(∑ms=1fstlnfst)/lnm
(s=1,2,…,m; t=1,2,…,n)(5)
式中:fst=bst/∑ms=1bst
当fst=0时,lnfst无意义,可将其修正为
fst=(1+bst)/∑ms=1(1+bst)(6)
(4)计算各项指标的熵权
W=(ωt)1×n,ωt=(1-Ht)/(n-∑nt=1Ht),
且满足∑nt=1ωt=1(7)
设第s个待测样本的综合联系度为μ-s,首先根据式(1)初步计算出列车过桥走行性评价样本s的联系度μs; 再用式(2)对评价样本做进一步的集对分析,得到第k项指标联系度μsk; 再次,运用(3)-(7)式计算各项指标的权重W=(ω1,ω2,…,ωn)T,其应该满足归一化条件∑nt=1ωt=1,即熵权-集对分析模型如下:
μ-s'=μs×∑nt=1(ωt×μst)(8)
式中,n为评价指标的个数.
最后,对μ^-'s中各分量进行归一化处理,可得到平均联系度μ^-s.分析比较综合联系度μ^-s中a、b、c的相对大小,并结合实际工程,即可判断出测点样本所属等级,其评判标准则为:(1)若max[a,b,c]=b,则判断为良好;(2)若max[a,b,c]=a,且a+b≥0.7,则判为优秀,否则判为良好;(3)若max[a,b,c]=c,且b+c≥0.7,则判为合格,否则判为良好.
以文献[1]中武广高速铁路典型桥梁动力特性和列车走行型评估报告中梁家湾特大桥的列车动力响应数据为例,在不同车速下共有6组实测数据,每一组数据样本有7项评价指标分别为脱轨系数、轮重减载率、轮轨横向力、车体垂向加速度、车体横向加速度、Sperling舒适性垂向指标、Sperling舒适性横向指标见表1.依据相关规范将列车走行性状况分为3级[18-19],Ⅰ级为优秀,Ⅱ级为良好,Ⅲ三级为合格,其中轮轨横向力指标门限值根据测试车轮静荷载计算得到,集体分级标准见表2.
表2 列车过桥走行性评价指标分级标准
Tab.2 The grading standards of the train crossing bridge behavior evaluation index
首先根据式(1)以列车走行性分级Ⅰ、Ⅱ、Ⅲ的级别作为分析对象联系度表达式中的同一度、差异度、对立度的取值界限,计算6组测量样本数据的集对分析联系度分别为: μ1=7/7+0i+0j、 μ2=6/7+0i+1/7j、μ3=3/7+3/7i+1/7j、μ4=3/7+2/7i+2/7j、μ5=3/7+3/7i+1/7j、μ6=4/7+2/7i+1/7j.
比较6个列车走行性评价样本的联系度可知:样本3、5属于同一级别,列车走行性状况从优到劣依次为1、2、6、3和5、4.为进一步分析实测指标的具体数值与列车走行性分级标准之间的数量关系,需对各测点样本数值做同一、差异、对立的进一步分析.
以车速序号3为例,由式(2)可计算出7项实测指标相对于列车走行性分级标准的联系度分别为:μ31=(0.8-0.72)/(0.8-0.6)+(0.72-0.6)/(0.8-0.6)i+0j=0.4+0.6i+0j
μ32=0+(0.8-0.66)/(0.8-0.65)i+(0.66-0.65)/(0.8-0.65)
j=0+0.933i+0.067j
μ33=1+0i+0j
μ34=(2.95-2.556)/(2.95-2.45)+(2.556-2.45)/(2.95-2.45)
i+0j=0.788+0.212i+0j
μ35=1+0i+0j,
μ36=(2.75-2.73)/(2.75-2.5)+(2.73-2.5)/(2.75-2.5)
i+0j=0.08+0.92i+0j
μ37=1+0i+0j
用MATLAB对熵值法确定权重系数过程(3)-(7)式进行编程如图1,分别输入6种不同的车速对应的6组实测数值,计算得到列车过桥走行性评价指标的各自权重见表3.
按式(8)计算出各样本的平均联系度为μ^-'1、μ^-'2、μ^-'3、μ^-'4、μ^-'5、μ^-'6,再对计算出的平均联系度进行归一化可得各样本的综合联系度μ^-1、μ^-2、μ^-3、μ^-4、μ^-5、μ^-6.最后根据评价准则判断出各车速下实测数据的列车走行性等级见表4.
各种车速下列车安全性和平稳性状况从优到劣依次为1、2、6、5、3、4.车速序号1、2、6列车过桥车速最小,所得列车安全性和平稳性状况最优,样本4列车过桥车速最大,其安全性和平稳性最差,整个评价结果与实际相符.
(1)本文应用数学上的集对分析理论构造列车走行性评价模型,在确定列车过桥运行安全性和平稳性评价的7个常规评价指标权重时引入了信息论中熵值理论,通过信息熵理论计算列车过桥运行实测数据本身的效用值来确定各个指标的权重,使得指标权重的分配有理可依.
(2)与常规列车过桥运行安全性和平稳性的等权重的指标权重评价方法相比,本文的评价方法使轮轴横向力这项指标的权重增大为0.165,该结果更加符合实际情况.
(3)列车过桥运行安全性和平稳性评价是一项比较复杂的工程,需要综合多项测量指标的属性进行判定,本文仅仅考虑了列车过桥运行安全性和平稳性评价的7个常规评价指标,而且列车过桥时要受到风荷载气动影响和温度变化的动态影响,故对其进行综合评价仅仅从权重角度考虑还有待欠缺.