一种高速列车制动系统热机耦合摩擦学行为预测方法

技术特征：
1.一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述预测方法包括：s1、获取磨合阶段完成后的接触压力分布，在参考接触压力分布和非均匀热流密度的前提下，进行正式服役阶段的制动系统热机耦合响应状态分析，计算制动系统的热机耦合响应；s2、读取某段中间时刻的响应状态作为有限元模型的初始状态，对正式服役阶段的热机耦合摩擦磨损行为进行预测，并实时更新有限元模型，当达到截止条件时记录并保存所有结果；s3、根据步骤s2的结果提取服役过程中某一磨损状态下的有限元结果并进一步更新有限元模型，对几何形貌变化的有限元模型进行振动加速度分析，从而对制动系统的摩擦振动噪声行为进行分析预测。2.根据权利要求1所述的一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述步骤s1具体包括以下内容：s11、将制动盘接触区域划分为n
h
个接触环带；s12、计算第j个接触环带的平均摩擦半径r
j
和面积s
j
；s13、读取正式服役阶段开始时的接触应力场p
s
，计算第j个接触环带的接触压力f
j
；s14、基于摩擦功率法，计算第j个接触环带的热流密度q
j
；s15、基于盘块热流分配系数，将摩擦热流按比例施加在制动盘与摩擦块表面；s16、根据温度位移分析程序完成正式服役阶段的热机耦合响应状态求解，记录并保存每个增量步的热机耦合响应状态结果。3.根据权利要求1所述的一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述步骤s2具体包括以下内容：s21、读取热机耦合响应状态结果，将正式服役阶段等分为l段，取第k段中间时刻的响应状态作为有限元模型的初始状态，计算器k的取值为1到l；s22、建立制动盘和摩擦块接触对，并进行短时间的热机耦合摩擦磨损行为计算；s23、更新计算器k＝k 1，如果k<l，则读取步骤s1求解的热机耦合响应状态结果，并再次提取下一段中间时刻的响应状态作为有限元模型的初始状态；s24、重复步骤s22和s23，直至正式服役阶段的热机耦合摩擦磨损行为预测完成，记录并保存所有结果。4.根据权利要求3所述的一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述进行短时间的热机耦合摩擦磨损行为计算包括：a1、将整段时间离散为t个增量步，初始化增量步计数器y＝1，根据温度位移分析程序进行热机耦合摩擦磨损分析；a2、读取有限元模型的节点坐标、节点温度、接触应力及相对滑移距离；a3、根据节点温度修正磨损系数并更新求解公式中的磨损系数，计算节点磨损深度和磨损方向；a4、反馈节点的磨损深度和磨损方向，根据任意拉格朗日-欧拉技术移动节点并更新坐标，同时在保证网格拓扑形式不变的前提下重绘网格，生成新的有限元模型；a5、判断增量步计数器y的大小，如果y<t，则更新计数器y＝y 1，生成新的有限元模型；
a6、重复步骤a2-a5，直至模拟结束，记录并输出摩擦块磨损状态w
k
，生成新的有限元模型。5.根据权利要求1所述的一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述步骤s3具体包括以下内容：s31、读取正式服役阶段的计算结果，提取服役过程中某一磨损状态下的有限元结果m
x
；s32、将有限元结果m
x
导入主程序生成新的有限元模型，并根据实际工况重新施加载荷，设置相互作用关系和边界条件；s33、启动模态分析程序，基于子空间投影法对几何形貌变化的有限元模型进行复特征值分析，或者利用隐式/显式动态分析程序对几何形貌变化的有限元模型进行振动加速度分析，从而对制动系统的摩擦振动噪声行为进行分析预测。6.根据权利要求1-5任意一项所述的一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述预测方法还包括有限元模型建立步骤和磨合阶段模拟计算步骤；所述预测方法在第一次运行时所述有限元模型建立步骤和磨合阶段模拟计算步骤执行于步骤s1之前，所述有限元模型建立步骤执行与所述磨合阶段模拟计算步骤之前；在所述预测方法第二次运行时以及之后运行不在执行所述有限元模型建立步骤和磨合阶段模拟计算步骤。7.根据权利要求6所述的一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述有限元模型建立步骤包括：根据制动系统设计需求，通过有限元软件对原始或者改进后的制动系统建立有限元模型；根据实际情况为有限元模型赋予材料属性，施加载荷并设置相互作用关系和边界条件。8.根据权利要求6所述的一种高速列车制动系统热机耦合摩擦学行为预测方法，其特征在于：所述磨合阶段模拟计算步骤包括：b1、将磨合过程离散为n个增量步，初始化增量步计数器i＝1；b2、开始第i个增量步的制动系统响应状态计算，记录并保存接触应力场和接触节点相对滑移距离；b3、进行第i个增量步的磨损计算，读取节点坐标和第i个增量步的节点接触应力及相对滑移距离，利用磨损求解公式计算节点磨损深度，设置磨损方向；b4、向有限元主程序反馈接触节点的磨损深度和磨损方向，根据任意拉格朗日-欧拉技术移动节点并更新坐标，同时在保证网格拓扑形式不变的前提下重绘网格，生成新的有限元模型；b5、更新计算器i＝i 1，开始下一个增量步的制动系统响应状态计算，记录并保存接触应力和接触节点相对滑移距离；b6、判断接触应力场是否达到相对稳定状态，如果没有达到，则重复步骤b3-b5，如果达到，则磨合阶段磨损计算完成，记录并输出磨合结束后，即正式服役阶段开始时的接触应力场p
s
和摩擦块磨损状态w
s
。

技术总结
本发明涉及一种高速列车制动系统热机耦合摩擦学行为预测方法，包括：S1、进行正式服役阶段的制动系统热机耦合响应状态分析，计算制动系统的热机耦合响应；S2、对正式服役阶段的热机耦合摩擦磨损行为进行预测，并实时更新有限元模型，当达到截止条件时记录并保存所有结果；S3、对几何形貌变化的有限元模型进行振动噪声行为预测。本发明充分考虑温度对材料摩擦磨损性能的影响，弥补了现有预测方法在进行磨损预测时忽略温度影响的缺陷，同时考虑了摩擦表面形貌变化对摩擦振动噪声的影响，提高了制动系统热机耦合摩擦学行为预测结果的准确性和有效性，为制动系统关键零部件的服役状态及使用寿命预测提供了有效手段，降低了高速列车的运维成本。的运维成本。的运维成本。


技术研发人员：卢纯 尹家宝 张庆贺 莫继良
受保护的技术使用者：西南交通大学
技术研发日：2021.11.25
技术公布日：2022/4/15