一种轮轨滚动接触斑分析方法技术

本发明专利技术公开了一种轮轨滚动接触斑分析方法，具体为：根据实际轮轨滚动接触关系，建立轮轨滚动接触关系多体动力学仿真模型，输出对应轮轨滚动接触斑的轮轨滚动接触变量；根据轮轨滚动接触斑的相关轮轨滚动接触变量，基于Kalker线性简化理论计算轮轨滚动接触斑内总滑动矢量的分布情况及界定接触斑内的黏滑区情况；基于Archard磨耗模型计算轮轨滚动接触斑内垂直磨耗量分布情况。本发明专利技术磨耗计算过程贴近实际现场磨耗情况，计算结果精确，且对于轮轨滚动接触斑的滑动效应、磨耗效应分析过程高效、简单，可作为SIMPACK、GENSYS等多体动力学仿真软件的轮轨滚动接触斑磨耗算法接口，实现对列车运行过程中轮对磨耗情况的仿真与计算。

An analysis method of wheel rail rolling contact spot

【技术实现步骤摘要】
一种轮轨滚动接触斑分析方法
本专利技术涉及轮轨滚动接触分析
，特别是一种轮轨滚动接触斑分析方法。
技术介绍
随着地铁系统的迅猛发展，地铁系统也面临着诸多问题。以地铁列车系统为例，轮对作为列车车体重要组成部件之一，在列车运行过程中起着引导列车车体沿着钢轨正确前进与转弯的作用，此外还要承受来自轨道和列车车体的全部动、静载荷，严峻的工作环境使得轮对磨耗问题已然成为地铁系统最常碰到的问题之一。实时掌握轮对磨耗情况，保障轮对正常的服役性能，对保证地铁列车正常安全运营具有重要意义。而研究轮对磨耗问题的核心在于分析轮轨滚动接触关系，轮轨滚动接触关系反映轮轨之间的匹配状态，是分析轮对磨耗问题的基础，其表现形式为轮轨滚动接触斑。而现有针对轮轨滚动接触斑的研究将接触斑视为一个整体，对其进行分析滑动效应及磨耗研究，导致最终分析结果准确性不高。
技术实现思路
本专利技术的目的在于提供一种能实现轮轨之间相对滑动效应分析、磨耗效应分析的高效、准确的轮轨滚动接触斑分析方法。实现本专利技术目的的技术解决方案为：一种轮轨滚动接触斑分析方法，包括以下步骤：步骤1、根据现场实际的轮轨滚动接触关系，建立轮轨滚动接触关系多体动力学仿真模型，输出对应轮轨滚动接触斑的轮轨滚动接触变量；步骤2、根据轮轨滚动接触斑的相关轮轨滚动接触变量，基于Kalker线性简化理论计算轮轨滚动接触斑内总滑动矢量w的分布情况及界定接触斑内的黏滑区情况；步骤3、基于Archard磨耗模型计算轮轨滚动接触斑内垂直磨耗量分布情况。进一步地，步骤1所述的根据现场实际的轮轨滚动接触关系，建立轮轨滚动接触关系多体动力学仿真模型，输出对应轮轨滚动接触斑的轮轨滚动接触变量，具体如下：步骤1.1，根据现场轮对基本接触参数，包括车轮直径、轮对横移量、轮对摇头角、轮缘内侧距，建立轮对多体动力学仿真模型；步骤1.2，根据现场轨道基本接触参数及轨道几何线性情况，建立轨道多体动力学仿真模型；步骤1.3，建立轮轨滚动接触关系多体动力学仿真模型，实时模拟列车运行过程中的实际轮轨滚动接触关系，并输出该过程中轮轨滚动接触斑的轮轨滚动接触变量，包括：接触斑数量np，横向η、纵向ε、自旋蠕滑率，轮上接触斑起始yws、终止ywe位置，接触斑半轴长a、b，接触斑垂向Qp接触应力，Kalker系数C11、C22、C23。进一步地，步骤2所述的根据轮轨滚动接触斑的相关轮轨滚动接触变量，基于Kalker线性简化理论计算轮轨滚动接触斑内总滑动矢量w的分布情况及界定接触斑内的黏滑区情况，具体如下：步骤2.1，基于轮轨滚动接触斑的轮轨滚动接触变量，得到轮轨稳态滚动接触状态下的滑动方程表达式；步骤2.2，对滑动方程表达式进行无量纲化操作，实现原轮轨滚动接触斑到单位圆接触斑的转换并得到总滑动矢量w的方程表达式；步骤2.3，将单位圆接触斑离散化为10×10个离散点，引入库伦摩擦定律分析，计算各离散点的黏滑情况及总滑动矢量w。进一步地，步骤2.1所述的基于轮轨滚动接触斑的轮轨滚动接触变量，得到轮轨稳态滚动接触状态下的滑动方程表达式，具体如下：步骤2.1.1、分别计算接触斑纵向、横向切向力px、py：其中，C11、C22、C23为Kalker系数，G为轮轨材料合成剪切模量，η、ε、分别为横向、纵向、自旋蠕滑率，a、b为椭圆接触斑半轴长；步骤2.1.2、分别计算接触斑纵向、横向上的材料柔度系数Lx、Ly：步骤2.1.3、基于Kalker线性简化理论得到轮轨稳态滚动接触状态下的滑动方程表达式：其中，vx、vy分别为纵向、横向滑动速度，vv为机车运行速度，Lx、Ly分别为纵向、横向材料柔度系数，px、py分别为纵向、横向切向力，x、y为基于轮轨坐标系接触斑内的坐标点。进一步地，步骤2.2所述的对滑动方程表达式进行无量纲化操作，实现原轮轨滚动接触斑到单位圆接触斑的转换并得到总滑动矢量w的方程表达式，具体如下：步骤2.2.1、对滑动方程表达式进行无量纲化操作，实现原轮轨滚动接触斑到单位圆接触斑的转换，令：其中，N为接触斑内法向合力，f为轮轨摩擦系数，x'为无量纲化后的横坐标，y'为无量纲化后的纵坐标，z0为接触斑内法向力最大值，p′x、p′y为无量纲化后接触斑内切向力，nx、ny为无量纲化后接触斑内的位移量，L'y为无量纲化后横向方向上的材料柔度系数，为无量纲化后的自旋蠕滑率分量，wi为滑动量；步骤2.2.2、得到新的滑动方程表达式为：步骤2.2.3、将新的滑动方程表达式转换为矩阵矢量形式：其中，w＝(wxwy)为总滑动矢量，为刚性滑动矢量，为弹性滑动矢量。进一步地，步骤2.3所述的将单位圆接触斑离散化为10×10个离散点，引入库伦摩擦定律分析，计算各离散点的黏滑情况及总滑动矢量w，具体如下：步骤2.3.1、取单位圆上任意一条平行于x’轴，宽为dy’的带，从带上的一点x1’＝x0’-h到x0’作积分，当h取值趋向于0时，将矩阵矢量表达式转换为：w1/2＝s1/2+p′1-p′0其中，为点x1’与点x0’之间的总滑动矢量，为点x1’与点x0’之间的刚性滑动矢量，p'1＝p'(x'1)为点x1’处的弹性滑动矢量，p'0＝p'(x'0)为点x0’处的弹性滑动矢量，h为点x1’与点x0’之间的长度，x'1/2为点x1’与点x0’之间的中点；步骤2.3.2、h取带长的1/10，从带右边界开始，在边界上的力理想化为0，故p’(x0’)＝(00)，向后按照依次对整条带的刚性滑动矢量进行求解，实现将单位圆接触斑离散化为10×10个离散点；步骤2.3.3、引入库伦摩擦定律分析，计算各离散点的黏滑情况及总滑动矢量w，即单位圆接触斑内任意一点都满足以下情况：定义z’＝p/z0，pH＝p0’-s│xo’-h/2，1)如果|pH|≤fz’，则该点为黏着状态，p1’＝pH且w│xo’-h/2＝(00)；2)如果|pH|>fz’，则该点为滑动状态，p1’＝(fz’/|pH|)×pH且w│xo’-h/2＝-λp1’，其中λ＝|pH|/(fz’)-1且λ>0。进一步地，步骤3所述的基于Archard磨耗模型计算轮轨滚动接触斑内垂直磨耗量分布情况，具体如下：步骤3.1，针对轮轨滚动接触斑内的滑动区，基于Archard磨耗模型计算各滑动离散点垂直磨耗量Δz(x,y)为：其中，Qp为法向接触力，H为两接触材料中较软材料的硬度，kw为无量纲的磨耗系数，Δx为离散点在横坐标上的长度；步骤3.2，针对轮轨滚动接触斑内的黏着区，默认Δz(x,y)为0。本专利技术与现有技术相比，其显著优点为：(1)将轮轨滚动接触斑的磨耗离散化为10×10个离散点的磨耗分布，磨耗计算过程更贴近实际现场磨耗情况，计算结果更为精确；(2)轮轨滚动接触斑的滑动效应、磨耗效本文档来自技高网...

【技术保护点】
1.一种轮轨滚动接触斑分析方法，其特征在于，包括以下步骤：/n步骤1、根据现场实际的轮轨滚动接触关系，建立轮轨滚动接触关系多体动力学仿真模型，输出对应轮轨滚动接触斑的轮轨滚动接触变量；/n步骤2、根据轮轨滚动接触斑的相关轮轨滚动接触变量，基于Kalker线性简化理论计算轮轨滚动接触斑内总滑动矢量w的分布情况及界定接触斑内的黏滑区情况；/n步骤3、基于Archard磨耗模型计算轮轨滚动接触斑内垂直磨耗量分布情况。/n

【技术特征摘要】
1.一种轮轨滚动接触斑分析方法，其特征在于，包括以下步骤：
步骤1、根据现场实际的轮轨滚动接触关系，建立轮轨滚动接触关系多体动力学仿真模型，输出对应轮轨滚动接触斑的轮轨滚动接触变量；
步骤2、根据轮轨滚动接触斑的相关轮轨滚动接触变量，基于Kalker线性简化理论计算轮轨滚动接触斑内总滑动矢量w的分布情况及界定接触斑内的黏滑区情况；
步骤3、基于Archard磨耗模型计算轮轨滚动接触斑内垂直磨耗量分布情况。


2.根据权利要求1所述的轮轨滚动接触斑分析方法，其特征在于，步骤1所述的根据现场实际的轮轨滚动接触关系，建立轮轨滚动接触关系多体动力学仿真模型，输出对应轮轨滚动接触斑的轮轨滚动接触变量，具体如下：
步骤1.1，根据现场轮对基本接触参数，包括车轮直径、轮对横移量、轮对摇头角、轮缘内侧距，建立轮对多体动力学仿真模型；
步骤1.2，根据现场轨道基本接触参数及轨道几何线性情况，建立轨道多体动力学仿真模型；
步骤1.3，建立轮轨滚动接触关系多体动力学仿真模型，实时模拟列车运行过程中的实际轮轨滚动接触关系，并输出该过程中轮轨滚动接触斑的轮轨滚动接触变量，包括：接触斑数量np，横向η、纵向ε、自旋蠕滑率，轮上接触斑起始yws、终止ywe位置，接触斑半轴长a、b，接触斑垂向Qp接触应力，Kalker系数C11、C22、C23。


3.根据权利要求1所述的轮轨滚动接触斑分析方法，其特征在于，步骤2所述的根据轮轨滚动接触斑的相关轮轨滚动接触变量，基于Kalker线性简化理论计算轮轨滚动接触斑内总滑动矢量w的分布情况及界定接触斑内的黏滑区情况，具体如下：
步骤2.1，基于轮轨滚动接触斑的轮轨滚动接触变量，得到轮轨稳态滚动接触状态下的滑动方程表达式；
步骤2.2，对滑动方程表达式进行无量纲化操作，实现原轮轨滚动接触斑到单位圆接触斑的转换并得到总滑动矢量w的方程表达式；
步骤2.3，将单位圆接触斑离散化为10×10个离散点，引入库伦摩擦定律分析，计算各离散点的黏滑情况及总滑动矢量w。


4.根据权利要求3所述的轮轨滚动接触斑分析方法，其特征在于，步骤2.1所述的基于轮轨滚动接触斑的轮轨滚动接触变量，得到轮轨稳态滚动接触状态下的滑动方程表达式，具体如下：
步骤2.1.1、分别计算接触斑纵向、横向切向力px、py：



其中，C11、C22、C23为Kalker系数，G为轮轨材料合成剪切模量，η、ε、分别为横向、纵向、自旋蠕滑率，a、b为椭圆接触斑半轴长；
步骤2.1.2、分别计算接触斑纵向、横向上的材料柔度系数Lx、Ly：



步骤2.1.3、基于Kalker线性简化理论得到轮轨稳态滚动接触状态下的滑动方程表达式：



其中，vx、vy分别为纵向、横向滑动速度，vv为机车运行速度，Lx、Ly分别为纵向、横向材料柔度系数，px、py分别为纵向、横向切向力，x、y为基于轮轨坐标系接触斑内的坐标点。


5.根据权利要求3所述的轮轨滚动接触斑分...

【专利技术属性】
技术研发人员：叶正君，徐筱筱，刘英舜，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏;32