铁路车辆轮轨接触几何与力学分析仿真系统

资 源 简 介

本项目旨在建立一套基于MATLAB的高精度轮轨接触分析平台，用于解决铁路车辆动力学中核心的轮对与钢轨相互作用问题。项目功能涵盖了从基础数据处理到复杂接触力学计算的全过程：1. 轮轨几何关系解算：采用迹线法或距离函数法，精确计算轮对在不同横移量、摇头角及侧滚角下的空间几何位置，自动搜索并确定轮轨接触点（包括单点及多点接触情况），计算接触角、滚动圆半径等几何参数。2. 等效锥度计算：依据UIC519标准，通过积分法或微分法计算轮对的等效锥度，分析不同轮轨匹配下的车辆运行稳定性。3. 接触力学分析：集成Hertz接触理论用于法向力计算，并结合Kalker线性理论、简化非线性理论（FASTSIM算法）或Polach模型进行切向蠕滑力与自旋力矩的求解，实现对轮轨界面应力分布的精确模拟。4. 廓形数据管理：支持标准廓形（如LM、LMA、S1002、CHN60等）及实测磨耗廓形的导入，内置三次样条插值与平滑算法以消除测量噪声对计算结果的影响。5. 接触斑可视化：提供二维及三维视图，直观展示接触斑形状、面积及内部应力云图分布。该系统可用于指导轮轨廓形优化设计、磨耗预测以及车辆脱轨安全性评估。

详 情 说 明

铁路车辆轮对轮轨接触关系综合分析仿真系统

项目简介

本项目建立了一套基于MATLAB的轮轨接触分析平台，用于解决铁路车辆动力学中核心的轮对与钢轨相互作用问题。该系统通过数值仿真手段，实现了从轮轨廓形数据生成、几何接触关系解算、等效锥度估算到接触斑内部力学特性分析的全流程计算。

系统的核心价值在于能够量化分析不同横移量下的轮轨几何参数变化，并通过集成的接触力学算法（Hertz + FASTSIM）求解轮轨界面的微观应力分布与蠕滑力，为车辆运行稳定性评估提供数据支撑。

主要功能特性

• 廓形数据生成与处理：内置LMA车轮踏面与CN60钢轨廓形的数学生成模型，支持廓形数据的三次样条插值和平滑滤波处理，以及基于轨距和轨底坡的坐标变换。
• 轮轨几何关系解算：基于最小距离法自动搜索左右轮轨接触点，计算不同横移量下的滚动圆半径差、接触角差及接触点空间位置。
• 等效锥度计算：依据UIC519标准思路，通过对小幅横移范围内的滚动圆半径差进行线性回归，快速估算名义滚动圆处的等效锥度。
• 接触力学分析：

* 法向力计算：基于Hertz接触理论，估算接触斑长短半轴及最大Hertz压力。 * 切向力求解：实现了简化的FASTSIM算法，考虑纵向、横向及自旋蠕滑率，求解接触斑内的剪切应力分布及总蠕滑力。

• 结果可视化：提供接触几何参数变化曲线及接触斑内部应力分布的直观展示。

系统要求

• MATLAB R2016a 或更高版本
• 无需额外工具箱（主要使用基础数学与插值函数）

使用方法

1. 确保MATLAB当前工作路径包含主脚本文件。
2. 直接运行主程序函数。
3. 程序将依次执行以下步骤：

* 生成并预处理轮轨廓形数据。 * 执行-10mm至+10mm范围内的横移工况循环，解算几何参数。 * 计算并输出名义等效锥度。 * 选取对中工况（横移量接近0），基于Hertz和FASTSIM算法计算接触斑应力。 * 绘制结果图表并在控制台输出关键力学指标（总纵向/横向蠕滑力、等效锥度值）。

代码实现逻辑与算法细节

1. 全局参数与工况设定

系统通过全局结构体定义了材料属性（弹性模量、泊松比、摩擦系数）和车辆几何参数（名义半径、轨距、轨底坡）。计算工况涵盖了横移量从-10mm到10mm的变化，并设定了单轮轮重。

2. 廓形数据管理

代码未直接读取外部文件，而是通过内部子函数模拟生成廓形数据：

• 车轮廓形：采用LMA磨耗型踏面模型，由R450圆弧、轮缘段和轮背段拼接而成。
• 钢轨廓形：采用CN60钢轨模型，由轨顶圆弧和侧面复合圆弧构成。
• 预处理：对生成的原始数据进行三次样条插值（Spline）以提高分辨率，并应用简单的移动平均滤波器（Window Size = 5）消除数据噪声。
• 坐标变换：钢轨数据利用旋转矩阵进行了坐标变换，考虑了左右轨距偏移（+/- L0/2）和1/40的轨底坡旋转。

3. 几何接触搜索算法

核心几何解算流程在主循环中实现：

• 横移模拟：将左/右车轮廓形在水平方向平移当前横移量。
• 接触点搜索：

* 确定轮轨横坐标的重叠区域。 * 在重叠区域内计算车轮与钢轨垂直方向的高度差（距离向量）。 * 搜索距离向量的极值点（最小垂直间隙）作为接触点。

• 几何参数计算：

* 接触角：在接触点位置对钢轨廓形进行微分（有限差分法），利用斜率的反正切值求得接触角。 * 滚动圆半径：基于接触点垂直坐标相对于名义半径的偏差进行估算。 * 几何偏差：计算左右轮的滚动圆半径差和接触角差。

4. 等效锥度计算

采用简化的差分/回归法计算：

• 选取横移量绝对值小于3mm的数据范围。
• 利用 polyfit 对该范围内的滚动圆半径差随横移量的数据进行一次多项式拟合。
• 取拟合直线斜率的一半作为等效锥度值（依据 $Delta r = 2 lambda y$ 的线性化假设）。

5. 接触力学模型

• Hertz接触：使用工程简化公式，根据主要接触区的等效曲率半径（设定的常数形式）和轮重，计算接触椭圆的长半轴 $a$、短半轴 $b$ 和最大中心压力 $P_0$。
• FASTSIM算法：

* 将接触椭圆离散化为 $31 times 31$ 的网格。 * 基于Kalker系数计算弹性柔度系数。 * 在网格节点上计算法向分布力（半椭球分布）。 * 计算刚性滑移（包含纵向、横向和自旋分量）。 * 应用线性弹性理论计算未饱和区的切向应力，并依据库伦摩擦定律（$mu P_z$）对饱和区应力进行截断修正，最终通过数值积分获得总蠕滑力 $F_x, F_y$ 和回转力矩 $M_z$。

6. 可视化

程序最后调用绘图函数（代码中引用但具体实现为了简洁未完全展开），传入廓形数据、几何解算结果结构体以及应力分布矩阵，实现数据的图形化展示。