技术领域
[0001] 本
发明属于无砟轨道检测领域，更具体地，涉及一种无砟轨道弹性扣件损伤识别方法。
背景技术
[0002] 国内外的高速
铁路轨道结构主要采用两种形式：有砟轨道与无砟轨道。其中高速铁路无砟轨道结构，由于其良好的整体性能，营运时极高的平顺性和
稳定性，已逐步在高速
铁路轨道结构中得到广泛应用。
[0003] 在高速铁路无砟轨道结构中，弹性扣件连接
钢轨和轨道板，起到固定钢轨的作用，是轨道结构中的重要组成部分。随着列车高速重载化的发展趋势，轨道结构的轨下
基础将受到很大的冲击和振动，导致弹性扣件产生不同程度的损伤甚至松脱。扣件松脱会加剧车辆和钢轨的振动，改变车辆-轨道耦合系统的动
力学性能，极大地影响列车运营的舒适性和安全性，特别是列车行驶时速高达300km/h。因此，对扣件损伤进行识别并及时的进行维护显得尤为重要。
[0004] 目前，在许多城市的铁路实际运营中，铁路结构损坏的检测通常依赖于目视检查。对弹性扣件的检测主要依靠人工沿着线路进行排查，这样传统的人工巡检方法效率较低，并且工作时间受
天窗时间的限制，而且现有扣件损伤检测技术存在一定的
缺陷，该方法只能检测扣件的外观状态，即是否缺失和断裂，无法同时精确地识别出损伤
位置和程度。
发明内容
[0005] 针对
现有技术的以上缺陷或改进需求，本发明提供了一种无砟轨道弹性扣件损伤识别方法，其中结合实际激励力数据和无损有限元模型，能够消除不确定因素的影响，准确识别弹性扣件的损伤位置和损伤程度，因而尤其适用于无砟轨道的检测。
[0006] 为实现上述目的，本发明提出了一种无砟轨道弹性扣件损伤识别方法，该方法包括如下步骤：
[0007] S1针对无损状态下的无砟轨道建立无损有限元模型并优化，然后确定弹性扣件的无损
刚度系数；
[0008] S2对待测无砟轨道施加激励，测得实际激励力数据和各
数据采集点的实际
加速度时程数据；
[0009] S3将所述实际激励力数据代入所述无损有限元模型中，获得模拟加速度时程数据；
[0010] S4根据所述步骤S2获得的实际加速度时程数据和所述步骤S3获得的模拟加速度时程数据，利用MCMC方法计算各个所述弹性扣件的实际刚度系数；
[0011] S5判断各个所述弹性扣件的实际刚度系数与无损刚度系数的比值是否为1，若是，则该弹性扣件没有损坏，若否，则该弹性扣件发生损坏。
[0012] 作为进一步优选地，所述步骤S1包括如下子步骤：
[0013] S11在无损状态下对无砟轨道施加激励，获得无损激励力数据和各数据采集点的无损加速度时程数据；
[0014] S12针对无损状态下的无砟轨道建立无损有限元模型；
[0015] S13将所述步骤S11获得的无损激励力数据代入所述无损有限元模型中获得模拟加速度时程数据；
[0016] S14将所述模拟加速度时程数据与所述子步骤S11获得的无损加速度时程数据进行比较，判断误差是否在预设范围内，若是，则直接转入所述子步骤S15，若否，则修正所述无损有限元模型，直至所述模拟加速度时程数据与所述无损加速度时程数据的误差在预设范围内，然后转入所述步骤S15；
[0017] S15根据所述步骤S11获得的无损加速度时程数据确定所述弹性扣件的无损刚度系数。
[0018] 作为进一步优选地，所述子步骤S11的具体过程为：在无砟轨道上确定数据采集点的位置，并在相应位置安装
传感器，在无损状态下选取任意一个数据采集点施加激励，通过所述传感器获得无损激励力数据和各数据采集点的无损加速度时程数据。
[0019] 作为进一步优选地，所述子步骤S11中传感器安装在钢轨和轨道板上。
[0020] 作为进一步优选地，所述子步骤S12中建立无损有限元模型的方法为：针对无损状态下无砟轨道的钢轨、弹性扣件、底座板、
水泥
沥青砂浆、轨道板以及基础进行仿真模拟，对所述轨道板与
水泥沥青砂浆之间以及所述水泥沥青砂浆与轨道板之间的层间关系进行仿真模拟，用
弹簧模拟所述弹性扣件，建立所述无损有限元模型。
[0021] 作为进一步优选地，所述步骤S4包括如下子步骤：
[0022] S41利用下式计算第一次取样即g＝1时的核
密度函数κg：
[0023]
[0024] 式中，g为取样次数，Nss为样本总数，W(j)为第j个样本的权重，N(θj(i),σ(j))为高斯函数；
[0025] S42根据下式对核密度函数κg进行更新，获得g＝gr时的核密度函数κgr：
[0026]
[0027] A为计算核密度函数时的系数，gr为
迭代次数；
[0028] S43利用下式计算目标拟合函数J(θ)：
[0029]
[0030] 式中，J(θ)为目标拟合函数，N0为总的数据点，k为数据采集点的编号，Nd为数据采集点的总数，j为数据中的序列号，NS为序列号的最大值，n为敲击的编号，Nf为敲击的总次数， 为第k个数据采集点的第j时间步长的实际加速度时程数据， 为第k个数据采集点的第j时间步长的模拟加速度时程数据，θ为弹性扣件的刚度系数，|| ||表示范数；
[0031] S44根据所述步骤S42获得的κgr和所述步骤S43获得的J(θ)，利用下式计算g＝gr时的边缘概率密度p(θ|D)：
[0032]
[0033] 式中，θ为各弹性扣件的刚度系数，D为实际加速度时程数据，c为归一化常数；
[0034] S45根据所述步骤S44获得的p(θ|D)反算所述弹性扣件的实际刚度系数。
[0035] 作为进一步优选地，所述迭代次数gr的取值优选为10次～20次。
[0036] 总体而言，通过本发明所构思的以上技术方案与现有技术相比，主要具备以下的技术优点：
[0037] 1.本发明通过建立无损有限元模型，能够综合考虑无砟轨道的各种不确定因素，使得在对待测无砟轨道进行测试时能够消除不确定因素的影响，仅考虑弹性扣件的损伤情况，从而有效提高无砟轨道弹性扣件损伤识别的精确度，通过MCMC法确定弹性扣件的损伤位置和损伤程度，对于无砟轨道后期的运营以及弹性扣件损伤识别与维护加固等工作具有前瞻性作用；
[0038] 2.同时，本发明在构建无损有限元模型时，将模拟加速度时程数据与无损加速度时程数据进行比较，保证误差在预设范围内，能够有效提高无损有限元模型的准确性；
[0039] 3.尤其是，本发明采用基于蒙特卡洛的贝叶斯修正方法对弹性扣件的实际损伤刚度系数进行优化，通过对核密度函数进行迭代优化，能够获得较为准确的联合概率密度，从而保证损伤识别的准确性，同时不需要进行前期的模型类选择，便可有效快速地进行检测。
附图说明
[0040] 图1是本发明提供的无砟轨道弹性扣件损伤识别方法的
流程图；
[0041] 图2是本发明步骤S4的流程图；
[0042] 图3是本发明优选
实施例中基于时程加速度数据的模型验证示意图；
[0043] 图4是本发明优选实施例中MCMC方法计算的实际刚度系数的分布示意图，其中(a)为弹性扣件1，(b)为弹性扣件2，(c)为弹性扣件3，(d)为弹性扣件4；
[0044] 图5是本发明优选实施中实际激励力数据与模型输入力的对比示意图；
[0045] 图6是本发明优选实施例中弹性扣件的边缘概率密度分布图，其中(a)为弹性扣件1，(b)为弹性扣件2，(c)为弹性扣件3，(d)为弹性扣件4；
[0046] 图7是本发明优选实施例中弹性扣件及轨道板的平面安装示意图；
[0047] 图8是本发明优选实施例中弹性扣件及轨道板的轴视图；
[0048] 图9是本发明优选实施例中识别结果示意图；
[0049] 图10是本发明优选实施例中MATLAB、Python语言以及ABAQUS交互式使用的代码流程图。
具体实施方式
[0050] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0051] 如图1所示，本发明提出了一种无砟轨道弹性扣件损伤识别方法，该方法包括如下步骤：
[0052] S1针对无损状态下的无砟轨道建立无损有限元模型并优化，然后确定弹性扣件的无损刚度系数；
[0053] S2对待测无砟轨道施加激励，测得实际激励力数据和各数据采集点的实际加速度时程数据；
[0054] S3将实际激励力数据代入无损有限元模型中，获得模拟加速度时程数据；
[0055] S4根据步骤S2获得的实际加速度时程数据和步骤S3获得的模拟加速度时程数据，利用MCMC(蒙特卡洛)方法计算各个所述弹性扣件的实际刚度系数；
[0056] S5判断各个弹性扣件的实际刚度系数与无损刚度系数的比值是否为1，若是，则该弹性扣件没有损坏，若否，则该弹性扣件发生损坏。
[0057] 进一步，步骤S1包括如下子步骤：
[0058] S11综合考虑无砟轨道结构的尺寸、识别精确度以及经济因素，确定钢轨上安装传感器的位置和个数，为保证结果的准确性，可将传感器安装在钢轨和轨道板上，然后在无损状态下选取钢轨上某数据采集点A进行动力特性测试，对其进行激励，通过设置的传感器获得无损激励力数据和各数据采集点的无损加速度时程数据；
[0059] S12针对无损状态下无砟轨道的钢轨、弹性扣件、底座板、水泥沥青砂浆、轨道板以及基础进行仿真模拟，对轨道板与水泥沥青砂浆之间以及水泥沥青砂浆与轨道板之间的层间关系进行仿真模拟，用弹簧模拟弹性扣件，建立无损有限元模型；
[0060] S13将所述步骤S11获得的无损激励力数据代入所述无损有限元模型中获得模拟加速度时程数据；
[0061] S14将所述模拟加速度时程数据与所述子步骤S11获得的无损加速度时程数据进行比较，判断误差是否在预设范围内，若是，则直接转入所述子步骤S15，若否，则修正所述无损有限元模型，直至所述模拟加速度时程数据与所述无损加速度时程数据的误差在预设范围内，然后转入所述步骤S15；
[0062] S15根据所述步骤S11获得的无损加速度时程数据确定所述弹性扣件的无损刚度系数。
[0063] 进一步，如图2所示，步骤S4包括如下子步骤：
[0064] S41利用式(1)计算第一次取样即g＝1时的核密度函数κg：
[0065]
[0066] 式中，g为取样次数，Nss为样本总数，W(j)为第j个样本的权重，N(θj(i),σ(j))为高斯函数，θj(i)为第j个样本，σ(j)为第j个样本高斯函数的方差；
[0067] S42根据式(2)对核密度函数κg进行更新，获得g＝gr时的核密度函数κgr：
[0068]
[0069] gr为迭代次数，根据试验经验确定gr优选为10次～20次；
[0070] S43利用式(3)计算目标拟合函数J(θ)：
[0071]
[0072] 式中，J(θ)为目标拟合函数，N0为总的数据点，k为数据采集点的编号，Nd为数据采集点的总数，j为数据中的序列号，NS为序列号的最大值，n为敲击的编号，Nf为敲击的总次数， 为第k个数据采集点的第j时间步长的实际加速度时程数据， 为第k个数据采集点的第j时间步长的模拟加速度时程数据，θ为弹性扣件的刚度系数，|| ||表示范数；
[0073] S44根据步骤S42获得的κgr和步骤S43获得的J(θ)，利用式(4)计算g＝gr时的边缘概率密度p(θ|D)：
[0074]
[0075] 式中，p(θ|D)为在现场测量响应数据D的情况下刚度系数为θ的概率。θ为各弹性扣件的刚度系数，D为实际加速度时程数据，c为归一化常数，使得概率密度在参数空间上的积分等于1；
[0076] S45根据步骤S44获得的p(θ|D)反算弹性扣件的实际刚度系数。
[0077] 下面针对具体实施例对本发明作进一步说明。
[0078] S1在无砟轨道的钢轨和铁道板上安装传感器，如图7、8所示，在无损状态下对无砟轨道施加激励，获得无损激励力数据和各数据采集点的无损加速度时程数据，针对无损状态下的无砟轨道建立无损有限元模型，图3是无损有限元模型中模拟加速度时程数据与无损加速度时程数据的对比，误差在预设范围内，说明该模型具有较高的精确度，然后根据无损加速度时程数据确定弹性扣件的无损刚度系数；
[0079] S2对待测无砟轨道施加激励，测得实际激励力数据和各数据采集点的实际加速度时程数据；
[0080] S3将实际加速度时程数据作为目标数据导入数学商用
软件MATLAB中，将实际激励力数据代入无损有限元模型中，获得模拟加速度时程数据；
[0081] S4根据所述步骤S2获得的实际加速度时程数据和所述步骤S3获得的模拟加速度时程数据，利用MCMC方法计算各个所述弹性扣件的实际刚度系数；
[0082] 更具体地， 误差满足均值为0，方差为κg的高斯分布，图10是MATLAB、Python语言以及ABAQUS交互式使用的代码流程图，详述如下：
[0083] 在已有的无砟轨道模型上，构建含弹性扣件损伤的CRTS板式无砟轨道的模型，并生成建模
数据库INP文件，在INP文件中寻找所述损伤弹性扣件刚度系数的行数；
[0084] 针对弹性扣件
弹性模量范围，设置一定数量初始值，根据迭代初始值，
修改INP文件中已寻找到的对应模拟弹性扣件的弹簧的参数，并生成新的INP文件；
[0085] 针对新生成的INP文件，提交ABAQUS计算，生成结果数据库ODB文件；
[0086] 编写Python语言，提取结果数据库ODB文件中预先设置的数据采集点的加速度时程数据，并存入TXT文件中；
[0087] 将数据采集点的加速度数据与预先设置的目标数据进行如式(3)计算，J(θ)为目标拟合函数，再根据MCMC贝叶斯概率理论方法，可求得弹性扣件识别结果的概率密度分布概率；
[0088] 图4是MCMC方法计算的实际刚度系数的分布示意图，基于上述MCMC贝叶斯模型修正方法所识别出来的无砟轨道弹性扣件损伤位置和程度，使用有限元软件ABAQUS构建具有对应损伤情况的有限元模型，施加对应现场的激励力，如图5现场试验力与模型输入力对比，并利用Python语言提取对应数据采集点的加速度时程；
[0089] S5判断各个弹性扣件的实际刚度系数与无损刚度系数的比值是否为1，若是，则该弹性扣件没有损坏，若否，则该弹性扣件发生损坏，结果如图6、9所示，弹性扣件4的比值为0.65，其损伤程度为35％。
[0090] 本发明中，将基于MCMC的贝叶斯概率理论方法与工程进行深度结合，并采用商业数学软件MATLAB与Python语言以及通用有限元软件ABAQUS的交互式使用，可以有效快速地测出无砟轨道结构弹性扣件的损伤位置和程度，并求出对应的概率分布，对于轨道结构后期的运营以及无砟轨道结构损伤维护加固，具有前瞻性作用，具有较强的工程应用价值。
[0091] 本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
