一种实测轮轨蠕滑曲线模拟方法技术

本发明专利技术涉及铁道工程领域，尤其涉及一种实测轮轨蠕滑曲线模拟方法。本发明专利技术引入四个修正参数μ

【技术实现步骤摘要】
一种实测轮轨蠕滑曲线模拟方法


[0001]本专利技术涉及铁道工程领域，尤其涉及一种实测轮轨蠕滑曲线模拟方法。

技术介绍

[0002]现有车辆
‑
轨道耦合动力学中黏着
‑
蠕滑曲线的模拟大多采用Polach的模型或其扩展形式，但无法获取轮轨局部接触解；科研人员为准确分析低黏着状态下轮轨接触斑内的接触应力以及滑移，应考虑黏着
‑
蠕滑特性、接触力学行为的耦合作用关系。Maksym Spiryagin等在Kalker简化理论的基础上提出了改进的FASTSIM算法，E.A.H.Vollebregt在CONTACT中扩展开发了考虑与滑移相关的摩擦系数和第三介质层的赫兹算法。
[0003]Polach模型、改进FASTSIM算法和CONTACT的扩展程序均可以模拟实测的蠕滑曲线，但存在着无法获取切向应力分布、切向应力分布误差大、参数选取缺乏科学性等问题。Polach在Kalker简化理论和线性理论的基础上，通过引入经验参数发展的计算模型可取得与实测蠕滑曲线相吻合的数据，但是Polach模型假定切向力在黏着区呈线性分布直接获取蠕滑力解析解，无法获取切向应力分布；虽然Polach模型应用较为广泛，但曲线拟合参数的获取并没有一套标准、高效的流程，多为研究人员主观选取，缺乏科学性；Maksym Spiryagin等在Kalker简化理论的基础上提出了改进的FASTSIM算法，由于算法中采用了不准确的抛物线切向接触边界，虽然可以得到轮轨接触的切向接触应力分布，但误差较大，约为5r/>‑
25％；E.A.H.Vollebregt在CONTACT中扩展开发了考虑与滑移相关的摩擦系数和第三介质层的赫兹算法，计算准确但计算成本高，无法应用于车辆
‑
轨道耦合动力学中。Sichani改进的FaStrip算法较大地提高了切向解的精度，在存在横向蠕滑和自旋的情况下，切向解的误差较于FASTSIM从18％降低到2％，同时计算成本较低，但是无法获得与实测蠕滑曲线相吻合的结果，需要进行扩展。

技术实现思路

[0004]针对
技术介绍
中存在的问题，提出一种实测轮轨蠕滑曲线模拟方法。
[0005]本专利技术引入四个修正参数μ
s
、A、B和f，对FaStrip算法进行修正，使其可考虑实测蠕滑曲线；提出拟合参数初值的计算方法，利用遗传算法获得与实测数据相吻合的拟合参数，使拟合蠕滑曲线与实测数据之间的误差降到最小，可输出较为准确的切向接触解，提高了拟合稳定性，减小了参数的误差，且通用性强。另外，此方法与原FaStrip算法对比，局部接触解更准确，更能反应出实测线路的黏着
‑
蠕滑特性。
[0006]本专利技术提出一种实测轮轨蠕滑曲线模拟方法，该方法基于修正的FaStrip算法，步骤如下：
[0007]S1、处理实测数据，获取待拟合曲线的折减系数控制区间和摩擦参数控制区间；
[0008]S2、引入四个修正参数μ
s
、A、B和f，确定参数的取值范围：
[0009]S21、斜率因子f的取值范围：根据实测试验确定Kalker系数c
11
，接触椭圆长半轴长度a，轮轨材料剪切模量G，第三介质层的厚度h
(3)
，第三介质层剪切模量G
(3)
的数值；将上述
参数值代入公式得到f的数值；将f的初步拟合数值的
±
10％作为遗传算法寻优的取值范围；
[0010]S22、变摩擦参数μ
s
、A、B的取值范围：提取摩擦系数控制区间的数据点；依据变摩擦公式μ＝μ
s
((1
‑
A)e
‑
Bω
+A)对摩擦系数控制区间的数据点进行非线性拟合；对拟合公式进行参数辨识，获得μ
s
、A、B的初步拟合数值；将μ
s
、A、B的初步拟合数值的
±
10％作为遗传算法寻优的取值范围；
[0011]S3、基于遗传算法进行参数拟合：
[0012]S31、确定目标函数；
[0013]S32、参数拟合：按照一定规则将参数即种群的个体转换为计算机可识别的数据串；随后计算目标函数，得到种群每个个体的适应度值，由适应度值选择个体，被选的个体会发生类似于自然界中染色体上的交叉和变异的现象；
[0014]S33、经过上述步骤完成一次进化，得到一个新的更符合要求的种群，然后再不断重复该过程，直到遗传代数达到设定值，精度达到要求时停止拟合过程。
[0015]优选的，变摩擦公式中μ
s
为最大摩擦系数，A和B为摩擦调节参数，ω是滑移速度矢量。
[0016]优选的，忽略自旋，ω可以表示为其中，V为机车速度，m/s；ξ、η分别为纵向和横向蠕滑率。
[0017]优选的，S1中需对实测数据进行分段，取实测蠕滑曲线的峰值作为分界点，对应蠕滑率称为滑行临界蠕滑率，小于滑行临界蠕滑率的部分数据在拟合时主要由折减系数控制，称为折减系数控制区间；大于滑行临界蠕滑率的部分数据在拟合时主要有摩擦参数控制，称为摩擦参数控制区间。
[0018]优选的，在S31中目标函数为其中μ
real
,μ
fit
分别代表拟合蠕滑曲线上的黏着系数和实测黏着系数，实验数据点的个数为n个。
[0019]优选的，在拟合蠕滑曲线上取与实测数据的点相同蠕滑率处的对应黏着系数，两者的差值R
i
即反应在该点上拟合曲线与实测数据的误差大小；此目标函数把每个点的误差平方相加得到整条曲线总的误差值R。
[0020]优选的，在S32中，交叉是参数之间的部分数据交换，变异是改变个体即参数编码数据串上某一位或某些位的值。
[0021]与现有技术相比，本专利技术具有如下有益的技术效果：
[0022]本专利技术通过引入四个修正参数μ
s
、A、B和f改进FaStrip算法，并基于遗传算法实现实测数据的拟合，得到的拟合曲线能达到很好的拟合效果，使拟合蠕滑曲线与实测数据之间的误差降到最小，可输出较为准确的切向接触解，提高了拟合精度，减小了参数的误差。另外，此方法与原FaStrip算法对比，局部接触解更准确，更能反应出实测线路的黏着
‑
蠕滑特性，具有更高的拟合精度和稳定性，且通用性强。
附图说明
[0023]图1为本专利技术一种实施例中实测数据示意图；
[0024]图2为本专利技术一种实施例中方法流程框图；
[0025]图3为本专利技术一种实施例中拟合效果示意图；
[0026]图4为本专利技术一种实施例中切向接触结构示意图(纵向蠕滑率为0.02)。
具体实施方式
[0027]本专利技术提出的一种实测轮轨蠕滑曲线模拟方法，该方法基于修正的FaStrip算法，步骤如下：
[0028]一、处理实测数据，获取待拟合曲线的折减系数控本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种实测轮轨蠕滑曲线模拟方法，其特征在于，该方法通过引入修正参数改进FaStrip算法，基于遗传算法实现对实测蠕滑曲线的模拟，步骤如下：S1、处理实测数据，获取待拟合曲线的折减系数控制区间和摩擦参数控制区间；S2、引入四个修正参数μ
s
、A、B和f，确定参数的取值范围：S21、斜率因子f的取值范围：根据实测试验确定Kalker系数c
11
，接触椭圆长半轴长度a，轮轨材料剪切模量G，第三介质层的厚度h
(3)
，第三介质层剪切模量G
(3)
的数值；将上述参数值代入公式得到f的数值；将f的初步拟合数值的
±
10％作为遗传算法寻优的取值范围；S22、变摩擦参数μ
s
、A、B的取值范围：提取摩擦系数控制区间的数据点；依据变摩擦公式μ＝μ
s
((1
‑
A)e
‑
Bω
+A)对摩擦系数控制区间的数据点进行非线性拟合；对拟合公式进行参数辨识，获得μ
s
、A、B的初步拟合数值；将μ
s
、A、B的初步拟合数值的
±
10％作为遗传算法寻优的取值范围；S3、基于遗传算法进行参数拟合：S31、确定目标函数；S32、参数拟合：按照一定规则将参数即种群的个体转换为计算机可识别的数据串；随后计算目标函数，得到种群每个个体的适应度值，由适应度值选择个体，被选的个体会发生类似于自然界中染色体上的交叉和变异的现象；S33、经过上述...

【专利技术属性】
技术研发人员：安博洋，田春香，王平，张洪吉，李粮余，莫宏愿，孙耀亮，宋娟，杨春凯，陈嵘，徐井芒，
申请(专利权)人：中铁二院工程集团有限责任公司，
类型：发明
国别省市：