本发明公开了一种两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，包括步骤：一、确定被测试制动器摩擦系数-温度间的函数关系μ＝f(T)：通过建立三维几何模型及三维热-机耦合有限元模型，对所建立三维热-机耦合有限元模型进行耦合温度场仿真分析得出；二、根据I型制动性能热衰退试验方法对被测试两轴汽车两轴汽车制动系统进行重复制动温升及各轴制动器摩擦性能虚拟测试；三、采用I型制动性能热衰退试验方法对被测试两轴汽车重复制动之后的热衰退制动效能进行虚拟测试。本发明设计合理、投资成本低、模拟效果好且操作简便，能解决现有I型制动性能热衰退虚拟系统及方法存在的成本高、模拟性较差、试验方法复杂、试验周期长等缺陷。
1.一种两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征在于该方法包括以下步骤：
步骤一、确定被测试两轴汽车制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)，被测试制动器为鼓式制动器或盘式制动器，其确定过程如下：
101、建立三维几何模型：根据被测试制动器的装配结构，利用CAD软件建立被测试制动器的三维几何模型；
102、建立三维热-机耦合有限元模型，其建模过程如下：
1021、将步骤101中所述三维几何模型导入多场耦合分析软件中；
1022、在多场耦合分析软件环境下，选定分析单元并以所选定的分析单元对所述三维几何模型进行网格划分，所选定的分析单元包括被测试制动器的三维实体结构单元、摩擦接触单元和对流辐射单元三种类型单元；
1023、输入根据工程材料手册所查询出的或者通过常规材料性能测试试验所测得的被测试制动器的材料物性参数，选择一种收敛快速的摩擦接触算法，同时选定摩擦系数-温度模型，最终建立三维热-机耦合有限元模型，所述摩擦系数-温度模型为被测试制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)，式中μ为被测试制动器制动摩擦副的摩擦系数，T为被测试制动器制动摩擦副的摩擦接触面温度；所述材料物性参数包括被测试制动器制动摩擦副所采用材料的弹性模量、密度、导热系数、比热容、热膨胀系数和泊松比；
103、对所建立的三维热-机耦合有限元模型进行耦合温度场仿真分析，在被测试制动器制动摩擦副上选定一个或多个测试点，将对所述三维热-机耦合有限元模型在各测试点上的温度场仿真结果与采用常规制动器温升试验方法对各测试点上所测试的温度测试结果进行比较，从而对步骤1023中所述函数关系μ＝f(T)进行验证，实现对函数关系μ＝f(T)中的系数进行准确辨识，最终得出被测试制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)；
步骤二、利用运算处理器且根据I型制动性能热衰退试验方法对被测试两轴汽车各轴制动器的温升及摩擦性能进行虚拟测试，其测试过程包括以下步骤：
201、初始参数设定、计算与存储：根据汽车I型制动性能热衰退试验方法规定的进行重复制动的试验工况，确定对被测试制动器进行虚拟试验的试验参数，包括被测试制动器所处车辆的初始速度V0且此时V0＝常量v、汽车制动减速度a0、被测试制动器的初始温度T0、循环制动次数和制动周期，并且根据所确定的循环制动次数和制动周期对总循环次数m和循环周期t进行设定；之后，调用参数设置单元输入所述试验参数和汽车及制动器的结构参数，同时根据测试精度及效率需求设定时间步长Δt，并计算出每一次循环制动过程的时间总步数k＝t/Δt，且将所输入的V0、a0、m、t、Δt、k和T0均同步存入存储单元内；同时调用参数计算单元并根据公式 计算得出被测试制动器中旋转部件的初始角速度ω0，并将初始角速度ω0分别转存为ωf和ωr后同步存入存储单元内；被测试制动器所处车辆为两轴汽车，则所述被测试制动器包括对两轴汽车的前轮进行制动的前轴制动器和对两轴汽车的后轮进行制动的后轴制动器；所述旋转部件为所述鼓式制动器的制动鼓或盘式制动器的制动盘；
202、建模：所述运算处理器根据汽车五自由度制动动力学方程组：
建立被测试制动器所处两轴汽车的动力学模型，所述两轴汽车的前轮和后轮的半径相同，式中∑FX、∑FZ、∑MY、∑Mwf和∑Mwr分别为汽车纵向合力、悬挂质量垂向合力、悬挂质量俯仰合力矩、前轮合力矩和后轮合力矩；Ffb和Frb为根据轮胎模型、滑移率和车轮载荷计算出来的汽车前、后轮纵向地面作用力；Fe、Ffk、Ffd、Frk、Frd、Mbf和Mbr分别为发动机制动力、前悬架弹性元件作用力、前悬架阻尼元件作用力、后悬架弹性元件作用力、后悬架阻尼元件作用力、前轴制动器输出力矩和后轴制动器输出力矩，且上述变量均为在预先设置初始值的基础上根据制动动力学方程组动态解算；M、Ms、Is、Iwf和Iwr分别为汽车质量、悬挂质量、悬挂质量绕质心转动惯量、前轮转动惯量和后轮转动惯量且在步骤201中进行设定；a、b和r分别为汽车质心到前轴距离、质心到后轴距离和车轮的滚动半径且在步骤201中进行设定；aX、asZ、εsY、εwf和εwr分别为汽车纵向加速度、悬挂质量垂向加速度、悬挂质量俯仰角加速度、前轮角加速度和后轮角加速度，且上述变量根据制动动力学方程组动态解算；
建立所述两轴汽车动力学模型的同时，将上述参数分别对应存入存储单元内，且此时
203、被测试两轴汽车制动热衰退循环过程中持续制动时制动器的管路压力值确定，其确定过程如下：
2031、调用管路压力值确定模块进行第一个时间步长Δt内制动器管路压力值的确定，其确定过程如下：
20311、先按照管路压力非线性递增规律确定前轴制动器和后轴制动器的管路压力递增量Δpf和Δpr，再根据公式pif＝p0f+Δpf和pir＝p0r+Δpr，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的管路压力pif和pir，并将pif和pir分别转存为pf和pr，对存储单元内的管路压力参数进行实时更新，其中p0f和p0r为被测试两轴汽车前轴制动器和后轴制动器的推出压耗；
20312、根据步骤20311中计算得出的pf和pr，计算得出前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并对存储单元内的输出力矩参数进行实时更新；
20313、结合步骤202中所述的发动机制动力Fe以及汽车前、后轮纵向地面作用力Ffb和Frb，并根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf和 计算得出被测试两轴汽车前后轮的角减速度εwf和εwr，并对存储单元内的角减速度参数进行实时更新；
20314、调用角速度计算模块且根据公式ωif＝ωf-εwf×Δt和ωir＝ωr-εwr×Δt，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir，并转存为ωf和ωr以对存储单元内的角速度参数进行实时更新；
20315、根据公式∑FX＝-2(Ffb+Frb)-Fe＝-MaX，计算得出此时间步长Δt结束时被测试两轴汽车的减速度aX；再根据公式Vi0＝V0-aX·Δt计算得出此时间步长Δt结束时被测试两轴汽车的车速Vi0并转存为V0，以对存储单元内的车速参数进行实时更新；
20316、结合ωf、ωr和V0，并根据公式 和 计
算得出前轮和后轮的滑移率sif和sir，并分别转存为sf和sr以对存储单元内的滑移率参数进行实时更新；
20317、根据轮胎模型且结合sf和sr，计算得出此时间步长Δt结束时汽车前、后轮纵向地面作用力Ffb和Frb并同步存入存储单元内，且相应对步骤202中所建立的动力学模型进行实时更新；
20318、根据车速V0并结合经常规试验测试得出的发动机制动力Fe与车速V0之间的函数关系，计算得出此时间步长Δt结束时发动机制动力Fe并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
20319、判断|aX-3|是否小于δ，其中δ＝0.01～0.001：当|aX-3|＜δ时，将步骤
20311中计算得出的pf和pr输出并另存为pmf和pmr，并将pmf和pmr作为制动热衰退循环过程中持续制动时前轴和后轴制动器的管路压力值，进入步骤204；否则，进入步骤2032；
2032、调用管路压力值确定模块进行下一时间步长Δt内制动器管路压力值的确定，其确定过程如下：
20321、先按照管路压力非线性递增规律确定前轴制动器和后轴制动器的管路压力递增量Δpf和Δpr，再根据公式pif＝pf+Δpf和pir＝pr+Δpr，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的管路压力pif和pir，并将pif和pir分别转存为pf和pr，以对存储单元内的管路压力参数进行实时更新；
20322、再根据步骤20321中计算得出的pf和pr，计算得出前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并对存储单元内的输出力矩参数进行实时更新；
20323、结合步骤202中所述的发动机制动力Fe以及汽车前、后轮纵向地面作用力Ffb和Frb，并根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf和 计算得出被测试两轴汽车前后轮的角减速度εwf和εwr，并对存储单元内的角减速度参数进行实时更新；
20324、调用角速度计算模块且根据公式ωif＝ωf-εwf×Δt和ωir＝ωr-εwr×Δt，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir，并转存为ωf和ωr以对存储单元内的角速度参数进行实时更新；
20325、根据公式∑FX＝-2(Ffb+Frb)-Fe＝-MaX，计算得出此时间步长Δt结束时被测试两轴汽车的减速度aX；再根据公式Vi0＝V0-aX·Δt计算得出此时间步长Δt结束时被测试两轴汽车的车速Vi0并转存为V0，以对存储单元内的车速参数进行实时更新；
20326、结合ωf、ωr和V0，并根据公式 和 计
算得出前轮和后轮的滑移率sif和sir，并分别转存为sf和sr，以对存储单元内的滑移率参数进行实时更新；
20327、根据轮胎模型且结合sf和sr计算得出此时间步长Δt结束时汽车前、后轮纵向地面作用力Ffb和Frb并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
20328、根据车速V0并结合经常规试验测试得出的发动机制动力Fe与车速V0之间的函数关系，计算得出此时间步长Δt结束时发动机制动力Fe并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
20329、判断|aX-3|是否小于δ，式中δ＝0.01～0.001：当|aX-3|＜δ时，将步骤
20321中计算得出的pf和pr输出并另存为pmf和pmr，并将pmf和pmr作为制动热衰退循环过程中持续制动时前轴和后轴制动器的管路压力值；否则，返回步骤20321；
204、循环制动：按照步骤201所设定的参数与步骤202中所建立的两轴汽车的动力学模型，在保持前轴制动器和后轴制动器的最大管路压力pf＝pmf且pr＝pmr的前提下，分m次分别对所述前轴制动器与后轴制动器的制动热衰退过程进行循环虚拟测试，其循环虚拟测试的测试过程如下：
2041、调用循环测试模块进行制动热衰退试验过程中第一个时间步长Δt内的试验，其试验过程如下：
20411、按照步骤20311至步骤20318，分别计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir、汽车前、后轮纵向地面作用力Ffb和Frb以及发动机制动力Fe，并将前轴制动器和后轴制动器的角速度ωif和ωir分别转存为ωf和ωr后，相应对步骤202中所建立的动力学模型进行实时更新；
20412、制动器温升计算，调用制动器温升计算模块对此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)分别进行计算，其计算过程如下：
20412a、根据公式pbf＝Mbfωf和pbr＝Mbrωr，计算出此时间步长Δt内的前轴制动器和和后轴制动器的吸收制动功率pbf和pbr，并将此时所求得的前轴制动器和后轴制动器的吸收制动功率pbf和pbr同步存入所述存储单元内，对所述存储单元内的吸收制动功率参数进行实时更新；
20412b、根据温 度变化微分 方程Pbf×dt-h×A[T(t)-TE]dt＝ m×c×dT(t)和
Pbr×dt-h×A[T(t)-TE]dt＝m×c×dT(t)，分别求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)，并将此时所求得的前轴制动器和后轴制动器的瞬时动态温度T(t)同步存入所述存储单元内，对所述存储单元内的前轴制动器和后轴制动器的瞬时动态温度参数进行实时更新；式中，h为前轴制动器或后轴制动器的对流换热系数，m为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热，TE为环境温度；
20413、热衰退摩擦性能计算：调用热衰退摩擦性能计算单元，根据步骤103中验证的函数关系μ＝f(T)且结合步骤20412b中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，相应分别计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr并相应存入存储单元内；
20414、调用所述参数计算单元，且根据制动器效能因数的常规计算方法且结合步骤
20413中计算得出的热态摩擦系数μf和μr，计算得出此时间步长Δt内前轴制动器和后轴制动器的效能因数BFf和BFr；再根据制动器输出力矩的常规计算方法，且结合此时所计算出来的效能因数BFf和BFr与前轴制动器管路压力pf和后轴制动器管路压力pr，相应推算得出此时间步长Δt结束时，前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并将所计算出来的效能因数BFf和BFr以及输出力矩Mbf和Mbr同步存入所述存储单元内，对所述存储单元内的制动器效能因数和输出力矩参数进行实时更新；同时，相应对步骤202中所建立的两轴汽车动力学模型进行实时更新；
20415、利用步骤20414中实时更新后的两轴汽车动力学模型，且根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf，调用参数计算单元计算得出此时间步长Δt结束时的前轮角加速度εwf和后轮角加速度εwr，并将此时所计算出来的前轮角加速度εwf和后轮角加速度εwr同步存入所述存储单元内，对所述存储单元内的前轮角加速度参数和后轮角加速度参数进行实时更新；
2042、调用循环测试模块进行制动测试过程中下一时间步长Δt内的测试，其测试过程如下：
20421、按照步骤20321至步骤20328，分别计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir、汽车前、后轮纵向地面作用力Ffb和Frb以及发动机制动力Fe，并将前轴制动器和后轴制动器的角速度ωif和ωir转存为ωf和ωr，并相应对步骤
20422、制动器温升计算，调用制动器温升计算模块对此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)分别进行计算，其计算过程如下：
20422a、根据公式pbf＝Mbfωf和pbr＝Mbrωr，计算出此时间步长Δt内的前轴制动器和和后轴制动器的吸收制动功率pbf和pbr，并将此时所求得的前轴制动器和后轴制动器的吸收制动功率pbf和pbr同步存入所述存储单元内，对所述存储单元内的吸收制动功率参数进行实时更新；
20422b、根据温 度变化微分 方程Pbf×dt-h×A[T(t)-TE]dt＝ m×c×dT(t)和
Pbr×dt-h×A[T(t)-TE]dt＝m×c×dT(t)，分别求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)，并将此时所求得的前轴制动器和后轴制动器的瞬时动态温度T(t)同步存入所述存储单元内，对所述存储单元内的前轴制动器和后轴制动器的瞬时动态温度参数进行实时更新；式中，h为前轴制动器或后轴制动器的对流换热系数，m为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热，TE为环境温度；
20423、热衰退摩擦性能计算：调用热衰退摩擦性能计算单元，根据步骤103中验证的函数关系μ＝f(T)且结合步骤20422b中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，相应分别计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr并相应存入存储单元内；
20424、调用所述参数计算单元，且根据制动器效能因数的常规计算方法且结合步骤
20423中计算得出的热态摩擦系数μf和μr，计算得出此时间步长Δt内前轴制动器和后轴制动器的效能因数BFf和BFr；再根据制动器输出力矩的常规计算方法，且结合此时所计算出来的效能因数BFf和BFr与前轴制动器管路压力pf和后轴制动器管路压力pr，相应推算得出此时间步长Δt结束时，前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并将所计算出来的效能因数BFf和BFr以及输出力矩Mbf和Mbr同步存入所述存储单元内，对所述存储单元内的制动器效能因数和输出力矩参数进行实时更新；同时，相应对步骤202中所建立的两轴汽车动力学模型进行实时更新；
20425、利用步骤20424中实时更新后的两轴汽车动力学模型，且根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf，调用参数计算单元计算得出此时间步长Δt结束时的前轮角加速度εwf和后轮角加速度εwr，并将此时所计算出来的前轮角加速度εwf和后轮角加速度εwr同步存入所述存储单元内，对所述存储单元内的前轮角加速度参数和后轮角加速度参数进行实时更新；
20426、本时间步长Δt结束时，对当前所述存储单元内更新后的被测试两轴汽车的车速V0进行判断：当 时，返回步骤2042，继续进行制动试验过程中的下一时间步长Δt内的试验测试，如此不断循环；当 时，说明已经完成本次热衰退试验虚拟测试过程中的制动试验过程，则进入步骤2043，此时本次热衰退试验中制动试验过程所用时间步数为n1；
2043、制动器提速散热阶段的制动器温度计算，其计算过程如下：
20431、调用制动器提速散热模块进行循环提速散热计算过程中第一个时间步长Δt内的计算，其计算过程如下：
20431a、根据公式-h×A[T(t)-TE]dt＝m×c×dT(t)，求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)并同步存入所述存储单元内，且对所述存储单元内的温度参数进行实时更新，式中TE为环境温度，h为前轴制动器或后轴制动器的对流换热系数，mbf为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热；
20431b、提速散热摩擦性能计算：根据步骤103中验证的函数关系μ＝f(T)且结合步骤20431a中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr，并相应存入存储单元内；
20432、调用制动器提速散热模块进行循环提速散热计算过程中下一个时间步长Δt内的计算，其计算过程如下：
20432a、根据公式-h×A[T(t)-TE]dt＝m×c×dT(t)，求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)并同步存入所述存储单元内，且对所述存储单元内的温度参数进行实时更新；
20432b、提速散热摩擦性能计算：根据步骤103中验证的函数关系μ＝f(T)且结合步骤20432a中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr，并相应存入存储单元内；
20433、返回步骤20432，继续进行循环提速散热计算过程中下一个时间步长Δt的计算，不断重复且经m1＝(t-n1×Δt)/Δt个时间步长计算后进入步骤2044，且本次热衰退试验中提速散热计算过程所用时间步数为m1；
2044、重复步骤2041至步骤2043，直至完成m次热衰退性能循环虚拟测试，则完成所述被测试两轴汽车前轴制动器和后轴制动器的温升及摩擦性能虚拟测试过程；
步骤二进行虚拟测试过程中，通过与所述运算处理器相接的显示器输出所述存储单元内所存储的所有计算数据进行同步显示，此时被测试两轴汽车各轴制动器的温升及摩擦性能虚拟测试过程结束。
2.按照权利要求1所述的两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征在于：步骤2044中完成所述被测试两轴汽车前轴制动器和后轴制动器的温升及摩擦性能虚拟测试过程后，还需利用运算处理器且根据I型制动性能热衰退试验方法对被测试两轴汽车重复制动之后的热衰退制动效能进行虚拟测试，其测试过程包括以下步骤：
301、初始参数设定与存储：根据I型制动性能热衰退试验方法对被测试两轴汽车重复制动之后的热衰退制动效能进行虚拟测试的试验工况，确定对被测试两轴汽车进行虚拟试验的试验参数，包括被测试制动器所处车辆的初始速度V2并同步存入存储单元内；
302、利用步骤202中所建立的两轴汽车动力学模型，计算得出在初始速度为V2且紧急停车状态下所述两轴汽车的汽车制动距离和制动减速度；
步骤302中进行虚拟测试过程中，通过与所述运算处理器相接的显示器输出所述存储单元内所存储的所有计算数据，汽车热衰退制动效能测试过程结束。
3.按照权利要求1或2所述的两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征在于：步骤1021中所述的多场耦合分析软件为ADINA软件。
4.按照权利要求1或2所述的两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征在于：步骤103中所述的常规制动器温升试验方法为采用底盘综合试验台对制动器后桥总成进行的定速温升试验方法或者采用制动器试验台对所述制动器后桥总成进行定速温升试验方法。
两轴汽车行车制动系I型制动性能热衰退虚拟试验方法
技术领域
[0001] 本发明涉及一种制动性能热衰退虚拟试验方法，尤其是涉及一种两轴汽车行车制动系I型制动性能热衰退虚拟试验方法。
背景技术
[0002] 汽车制动器的实际工作过程比较复杂，而制动器作为汽车的一个重要安全部件，它的实际工作性能及稳定性不仅受到“人-车-环境”系统中其它非车因素影响，同时还与汽车其它系统和总成相互耦合与作用。对于汽车行车制动系I型制动性能热衰退虚拟试验系统的建立，不但需要基于制动器台架试验及其有限元仿真相结合所验证得到的制动器的摩擦系数温度特性，同时还需要综合考虑汽车重量、几何和其它有关部件性能参数作用下对制动器动态输出力矩的影响，并设法进行误差修正。
[0003] 汽车整车制动试验主要为验证是否符合相应的法规、标准的要求而进行的，汽车整车制动方面主导性的国家强制性标准主要是GB12676-1999《汽车制动系统结构、性能和试验方法》和GB7258-1997《机动车运行安全技术条件》，其中1999年修订后的GB12676-1999等效采用了ECE R13法规和ISO 6597-1991，ISO7634-1995，ISO7635-1991等国际标准；GB7258是根据我国机动车运行中的安全要求，由公安部、交通部提出并于1987年开始实施，1997年作了修订，其第6节“制动系”的要求与GB12676-1999的内容大体相似。GB12676-1999《汽车制动系统结构、性能和试验方法》涉及的整车热衰退试验包括：I型试验(紧急制动热衰退和恢复试验)和II型试验(下长坡试验)。汽车实车道路制动试验特点是真实可靠，一般在汽车试验场进行，可以模拟一些车辆的实际使用工况，是汽车性能研究不可缺少的基础和重要环节，但它需要投入较多的人力、物力和各种测量仪器，费用及试验周期仍然相对较长；另外有时由于试验条件和试验安全方面考虑也有所受限。而虚拟试验为汽车试验提供了特殊高效解决手段。虚拟试验就是借助于计算机的高速解算，按实际试验要求对基于描述汽车系统及其动态过程的数学模型的仿真模型进行模拟试验：不仅可以作为真实试验的前期准备工作或在一定程度上替代传统的试验(如一些极限工况)；大幅减少真实试验次数，降低试验费用，缩短试验周期；较好的交互性，使得各种试验信息及时反馈；不受气象条件、场地、时间和次数的限制，试验过程可方便实现回放、再现和重复。
发明内容
[0004] 本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其设计合理、试验投资成本低、模拟效果好且使用操作简便，能有效解决现有制动系I型制动性能热衰退虚拟测量系统及方法所存在的投资成本高、模拟性较差、试验方法复杂、试验周期长等多种缺陷和不足。
[0005] 为解决上述技术问题，本发明采用的技术方案是：一种两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征在于该方法包括以下步骤：
[0006] 步骤一、确定被测试两轴汽车制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)，被测试制动器为鼓式制动器或盘式制动器，其确定过程如下：
[0007] 101、建立三维几何模型：根据被测试制动器的装配结构，利用CAD软件建立被测试制动器的三维几何模型；
[0008] 102、建立三维热-机耦合有限元模型，其建模过程如下：
[0009] 1021、将步骤一中所述三维几何模型导入多场耦合分析软件中；
[0010] 1022、在多场耦合分析软件环境下，选定分析单元并以所选定的分析单元对所述三维几何模型进行网格划分，所选定的分析单元包括被测试制动器的三维实体结构单元、摩擦接触单元和对流辐射单元三种类型单元；
[0011] 1023、输入根据工程材料手册所查询出的或者通过常规材料性能测试试验所测得的被测试制动器的材料物性参数，选择一种收敛快速的摩擦接触算法，同时选定摩擦系数-温度模型，最终建立三维热-机耦合有限元模型，所述摩擦系数-温度模型为被测试制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)，式中μ为被测试制动器制动摩擦副的摩擦系数，T为被测试制动器制动摩擦副的摩擦接触面温度；所述材料物性参数包括被测试制动器制动摩擦副所采用材料的弹性模量、密度、导热系数、比热容、热膨胀系数和泊松比；
[0012] 103、对所建立的三维热-机耦合有限元模型进行耦合温度场仿真分析，在被测试制动器制动摩擦副上选定一个或多个测试点，将对所述三维热-机耦合有限元模型在各测试点上的温度场仿真结果与采用常规制动器温升试验方法对各测试点上所测试的温度测试结果进行比较，从而对步骤3中所述函数关系μ＝f(T)进行验证，实现对函数关系μ＝f(T)中的系数进行准确辨识，最终得出被测试制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)；
[0013] 步骤二、利用运算处理器且根据I型制动性能热衰退试验方法对被测试两轴汽车各轴制动器的温升及摩擦性能进行虚拟测试，其测试过程包括以下步骤：
[0014] 201、初始参数设定、计算与存储：根据汽车I型制动性能热衰退试验方法规定的进行重复制动的试验工况，确定对被测试制动器进行虚拟试验的试验参数，包括被测试制动器所处车辆的初始速度V0且此时V0＝常量v、汽车制动减速度a0、被测试制动器的初始温度T0、循环制动次数和制动周期，并且根据所确定的循环制动次数和制动周期对总循环次数m和循环周期t进行设定；之后，调用参数设置单元输入所述试验参数和汽车及制动器的结构参数，同时根据测试精度及效率需求设定时间步长Δt，并计算出每一次循环制动过程的时间总步数k＝t/Δt，且将所输入的V0、a0、m、t、Δt、k、ω0和T0均同步存入存储单元内；同时调用参数计算单元并根据公式 计算得出被测试制动器中旋转部件的初始角速度ω0，并将初始角速度ω0分别转存为ωf和ωr后同步存入存储单元内；被测试制动器所处车辆为两轴汽车，则所述被测试制动器包括对两轴汽车的前轮进行制动的前轴制动器和对两轴汽车的后轮进行制动的后轴制动器；所述旋转部件为所述鼓式制动器的制动鼓或盘式制动器的制动盘；
[0015] 202、建模：所述运算处理器根据汽车五自由度制动动力学方程组：
[0016] 建立被测试制动器所处两轴汽车的动力学模型，所述两轴汽车的前轮和后轮的半径相同，式中∑FX、∑FZ、∑MY、∑Mwf和∑Mwr分别为汽车纵向合力、悬挂质量垂向合力、悬挂质量俯仰合力矩、前轮合力矩和后轮合力矩；Ffb和Frb为根据轮胎模型、滑移率和车轮载荷计算出来的汽车前、后轮纵向地面作用力；Fe、Ffk、Ffd、Frk、Frd、Mbf和Mbr分别为发动机制动力、前悬架弹性元件作用力、前悬架阻尼元件作用力、后悬架弹性元件作用力、后悬架阻尼元件作用力、前轴制动器输出力矩和后轴制动器输出力矩，且上述变量均为在预先设置初始值的基础上根据制动动力学方程组动态解算；M、Ms、Is、Iwf和Iwr分别为汽车质量、悬挂质量、悬挂质量绕质心转动惯量、前轮转动惯量和后轮转动惯量且在步骤201中进行设定；a、b和r分别为汽车质心到前轴距离、质心到后轴距离和车轮的滚动半径且在步骤201中进行设定；aX、asZ、εsY、εwf和εwr分别为汽车纵向加速度、悬挂质量垂向加速度、悬挂质量俯仰角加速度、前轮角加速度和后轮角加速度，且上述变量根据制动动力学方程组动态解算；
[0017] 建立所述两轴汽车动力学模型的同时，将上述参数分别对应存入存储单元内，且此时
[0018] 203、被测试两轴汽车制动热衰退循环过程中持续制动时制动器的管路压力值确定，其确定过程如下：
[0019] 2031、调用管路压力值确定模块进行第一个时间步长Δt内制动器管路压力值的确定，其确定过程如下：
[0020] 20311、先按照管路压力非线性递增规律确定前轴制动器和后轴制动器的管路压力递增量Δpf和Δpr，再根据公式pif＝p0f+Δpf和pir＝p0r+Δpr，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的管路压力pif和pir，并将pif和pir分别转存为pf和pr，对存储单元内的管路压力参数进行实时更新，其中p0f和p0r为被测试两轴汽车前轴制动器和后轴制动器的推出压耗；
[0021] 20312、根据步骤20311中计算得出的pf和pr，计算得出前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并对存储单元内的输出力矩参数进行实时更新；
[0022] 20313、结合步骤202中所述的发动机制动力Fe以及汽车前、后轮纵向地面作用力Ffb和Frb，并根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf和计算得出被测试两轴汽车前后轮的角减速度εwf和εwr，并对存储单元内的角减速度参数进行实时更新；
[0023] 20314、调用角速度计算模块且根据公式ωif＝ωf-εwf×Δt和ωir＝ωr-εwr×Δt，计算得出此时间步长Δt结束时前轴制动器或后轴制动器的角速度ωif和ωir并转存为ωf和ωr以对存储单元内的角速度参数进行实时更新；
[0024] 20315、根据公式∑FX＝-2(Ffb+Frb)-Fe＝-MaX，计算得出此时间步长Δt结束时被测试两轴汽车的减速度aiX并转存为aX；再根据公式Vi0＝V0-aX·Δt计算得出此时间步长Δt结束时被测试两轴汽车的车速Vi0并转存为V0，以对存储单元内的车速参数进行实时更新；
[0025] 20316、结 合 ωf、ωr 和 V0，并 根 据 公 式 和计算得出前轮和后轮的滑移率sif和sir并分别转存为sf和sr以对存
储单元内的滑移率参数进行实时更新；
[0026] 20317、根据轮胎模型且结合sf和sr，计算得出此时间步长Δt结束时汽车前、后轮纵向地面作用力Ffb和Frb并同步存入存储单元内，且相应对步骤202中所建立的动力学模型进行实时更新；
[0027] 20318、根据车速V0并结合经常规试验测试得出的发动机制动力Fe与车速V0之间的函数关系，计算得出此时间步长Δt结束时发动机制动力Fe并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
[0028] 20319、判断|aX-3|是否小于δ，其中δ＝0.01～0.001：当|aX-3|＜δ时，将步骤20311中计算得出的pf和pr输出并另存为pmf和pmr，并将pmf和pmr作为制动热衰退循环过程中持续制动时前轴和后轴制动器的管路压力值；否则，进入步骤2032；
[0029] 2032、调用管路压力值确定模块进行下一时间步长Δt内制动器管路压力值的确定，其确定过程如下：
[0030] 20321、先按照管路压力非线性递增规律确定前轴制动器和后轴制动器的管路压力递增量Δpf和Δpr，再根据公式pif＝pf+Δpf和pir＝pr+Δpr，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的管路压力pif和pir并将pif和pir分别转存为pf和pr，以对存储单元内的管路压力参数进行实时更新；
[0031] 20322、再根据步骤20321中计算得出的pf和pr，计算得出前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并对存储单元内的输出力矩参数进行实时更新；
[0032] 20323、结合步骤202中所述的发动机制动力Fe以及汽车前、后轮纵向地面作用力Ffb和Frb，并根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf和计算得出被测试两轴汽车前后轮的角减速度εwf和εwr，并对存储单元内的角减速度参数进行实时更新；
[0033] 20324、调用角速度计算模块且根据公式ωif＝ωf-εwf×Δt和ωir＝ωr-εwr×Δt，计算得出此时间步长Δt结束时前轴制动器或后轴制动器的角速度ωif和ωir，并转存为ωf和ωr以对存储单元内的角速度参数进行实时更新；
[0034] 20325、根据公式∑FX＝-2(Ffb+Frb)-Fe＝-MaX，计算得出此时间步长Δt结束时被测试两轴汽车的减速度aiX并转存为aX；再根据公式Vi0＝V0-aX·Δt计算得出此时间步长Δt结束时被测试两轴汽车的车速Vi0并转存为V0，以对存储单元内的车速参数进行实时更新；
[0035] 20326、结 合 ωf、ωr 和 V0，并 根 据 公 式 和计算得出前轮和后轮的滑移率sif和sir并分别转存为sf和sr，以对存
储单元内的滑移率参数进行实时更新；
[0036] 20327、根据轮胎模型且结合sf和sr计算得出此时间步长Δt结束时汽车前、后轮纵向地面作用力Ffb和Frb并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
[0037] 20328、根据车速V0并结合经常规试验测试得出的发动机制动力Fe与车速V0之间的函数关系，计算得出此时间步长Δt结束时发动机制动力Fe并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
[0038] 20329、判断|aX-3|是否小于δ，式中δ＝0.01～0.001：当|aX-3|＜δ时，将步骤20321中计算得出的pf和pr输出并另存为pmf和pmr，并将pmf和pmr作为制动热衰退循环过程中持续制动时前轴和后轴制动器的管路压力值；否则，返回步骤20321；
[0039] 204、循环制动：按照步骤201所设定的参数与步骤202中所建立的两轴汽车的动力学模型，在保持前轴制动器和后轴制动器的最大管路压力pf＝pmf且pr＝pmr的前提下，分m次分别对所述前轴制动器与后轴制动器的制动热衰退过程进行循环虚拟测试，其循环虚拟测试的测试过程如下：
[0040] 2041、调用循环测试模块进行制动热衰退试验过程中第一个时间步长Δt内的试验，其试验过程如下：
[0041] 20411、按照步骤20311至步骤20318，分别计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir、汽车前、后轮纵向地面作用力Ffb和Frb以及发动机制动力Fe，并将前轴制动器和后轴制动器的角速度ωif和ωir分别转存为ωf和ωr后，相应对步骤202中所建立的动力学模型进行实时更新；
[0042] 20412、制动器温升计算，调用制动器温升计算模块对此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)分别进行计算，其计算过程如下：
[0043] 20412a、根据公式pbf＝Mbfωf，计算出此时间步长Δt内的前轴制动器和和后轴制动器的吸收制动功率pbf和pbr，并将此时所求得的前轴制动器和后轴制动器的吸收制动功率pbf和pbr同步存入所述存储单元内，对所述存储单元内的吸收制动功率参数进行实时更新；
[0044] 20412b、根据温度变化微分方程Pbf×dt-h×A[T(t)-TE]dt＝m×c×dT(t)和Pbr×dt-h×A[T(t)-TE]dt＝m×c×dT(t)，分别求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)，并将此时所求得的前轴制动器和后轴制动器的瞬时动态温度T(t)同步存入所述存储单元内，对所述存储单元内的前轴制动器和后轴制动器的瞬时动态温度参数进行实时更新；式中，h为前轴制动器或后轴制动器的对流换热系数，m为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热，TE为环境温度；
[0045] 20413、热衰退摩擦性能计算：调用热衰退摩擦性能计算单元，根据步骤103中验证的函数关系μ＝f(T)且结合步骤20412b中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，相应分别计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr并相应存入存储单元内；
[0046] 20414、调用所述参数计算单元，且根据制动器效能因数的常规计算方法且结合步骤20413中计算得出的热态摩擦系数μf和μr，计算得出此时间步长Δt内前轴制动器和后轴制动器的效能因数BFf和BFr；再根据制动器输出力矩的常规计算方法，且结合此时所计算出来的效能因数BFf和BFr与前轴制动器管路压力pf和后轴制动器管路压力pr，相应推算得出此时间步长Δt结束时，前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并将所计算出来的效能因数BFf和BFr以及输出力矩Mbf和Mbr同步存入所述存储单元内，对所述存储单元内的制动器效能因数和输出力矩参数进行实时更新；同时，相应对步骤203中所建立的两轴汽车动力学模型进行实时更新；
[0047] 20415、利用步骤20414中实时更新后的两轴汽车动力学模型，且根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf，调用参数计算模块计算得出此时间步长Δt结束时的前轮角加速度εwf和后轮角加速度εwr，并将此时所计算出来的前轮角加速度εwf和后轮角加速度εwr同步存入所述存储单元内，对所述存储单元内的前轮角加速度参数进行实时更新；
[0048] 2042、调用循环测试模块进行制动测试过程中下一时间步长Δt内的测试，其测试过程如下：
[0049] 20421、按照步骤20321至步骤20328，分别计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir、汽车前、后轮纵向地面作用力Ffb和Frb以及发动机制动力Fe，并将前轴制动器和后轴制动器的角速度ωif和ωir转存为ωf和ωr，并相应对步骤202中所建立的动力学模型进行实时更新；
[0050] 20422、制动器温升计算，调用制动器温升计算模块对此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)分别进行计算，其计算过程如下：
[0051] 20422a、根据公式pbf＝Mbfωf，计算出此时间步长Δt内的前轴制动器和和后轴制动器的吸收制动功率pbf和pbr，并将此时所求得的前轴制动器和后轴制动器的吸收制动功率pbf和pbr同步存入所述存储单元内，对所述存储单元内的吸收制动功率参数进行实时更新；
[0052] 20422b、根据温度变化微分方程Pbf×dt-h×A[T(t)-TE]dt＝m×c×dT(t)和Pbr×dt-h×A[T(t)-TE]dt＝m×c×dT(t)，分别求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)，并将此时所求得的前轴制动器和后轴制动器的瞬时动态温度T(t)同步存入所述存储单元内，对所述存储单元内的前轴制动器和后轴制动器的瞬时动态温度参数进行实时更新；式中，h为前轴制动器或后轴制动器的对流换热系数，m为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热，TE为环境温度；
[0053] 20423、热衰退摩擦性能计算：调用热衰退摩擦性能计算单元，根据步骤103中验证的函数关系μ＝f(T)且结合步骤20412b中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，相应分别计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr并相应存入存储单元内；
[0054] 20424、调用所述参数计算单元，且根据制动器效能因数的常规计算方法且结合步骤20423中计算得出的热态摩擦系数μf和μr，计算得出此时间步长Δt内前轴制动器和后轴制动器的效能因数BFf和BFr；再根据制动器输出力矩的常规计算方法，且结合此时所计算出来的效能因数BFf和BFr与前轴制动器管路压力pf和后轴制动器管路压力pr，相应推算得出此时间步长Δt结束时，前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并将所计算出来的效能因数BFf和BFr以及输出力矩Mbf和Mbr同步存入所述存储单元内，对所述存储单元内的制动器效能因数和输出力矩参数进行实时更新；同时，相应对步骤203中所建立的两轴汽车动力学模型进行实时更新；
[0055] 20425、利用步骤20424中实时更新后的两轴汽车动力学模型，且根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf，调用参数计算模块计算得出此时间步长Δt结束时的前轮角加速度εwf和后轮角加速度εwr，并将此时所计算出来的前轮角加速度εwf和后轮角加速度εwr同步存入所述存储单元内，对所述存储单元内的前轮角加速度参数和后轮角加速度参数进行实时更新；
[0056] 20426、本时间步长Δt结束时，对当前所述存储单元内更新后的被测试两轴汽车的车速V0进行判断：当 时，返回步骤2042，继续进行制动试验过程中的下一时间步长Δt内的试验测试，如此不断循环；当 时，说明已经完成本次热衰退试验虚拟测试过程中的制动试验过程，则进入步骤2043，此时本次热衰退试验中制动试验过程所用时间步数为n1；
[0057] 2043、制动器提速散热阶段的制动器温度计算，其计算过程如下：
[0058] 20431、调用制动器提速散热模块进行循环提速散热计算过程中第一个时间步长Δt内的计算，其计算过程如下：
[0059] 20431a、根据公式-h×A[T(t)-TE]dt＝m×c×dT(t)，求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)并同步存入所述存储单元内，且对所述存储单元内的温度参数进行实时更新，式中TE为环境温度，h为前轴制动器或后轴制动器的对流换热系数，mbf为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热；
[0060] 20431b、提速散热摩擦性能计算：根据步骤103中验证的函数关系μ＝f(T)且结合步骤20431a中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr，并相应存入存储单元内；
[0061] 20432、调用制动器提速散热模块进行循环提速散热计算过程中下一个时间步长Δt内的计算，其计算过程如下：
[0062] 20432a、根据公式-h×A[T(t)-TE]dt＝m×c×dT(t)，求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)并同步存入所述存储单元内，且对所述存储单元内的温度参数进行实时更新；
[0063] 20432b、提速散热摩擦性能计算：根据步骤103中验证的函数关系μ＝f(T)且结合步骤20432a中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr，并相应存入存储单元内；
[0064] 20433、返回步骤20432，继续进行循环提速散热计算过程中下一个时间步长Δt的计算，不断重复且经m1＝(t-n1×Δt)/Δt个时间步长计算后进入步骤20434，且本次热衰退试验中提速散热计算过程所用时间步数为m1；
[0065] 2044、重复步骤2041至步骤2043，直至完成m次热衰退性能循环虚拟测试，则完成所述被测试两轴汽车前轴制动器和后轴制动器的温升及摩擦性能虚拟测试过程；
[0066] 步骤二进行虚拟测试过程中，通过与所述运算处理器相接的显示器输出所述存储单元内所存储的所有计算数据进行同步显示，此时被测试两轴汽车各轴制动器的温升及摩擦性能虚拟测试过程结束。
[0067] 上述两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征是：步骤2044中完成所述被测试两轴汽车前轴制动器的温升及摩擦性能虚拟测试过程后，还需利用运算处理器且根据I型制动性能热衰退试验方法对被测试两轴汽车重复制动之后的热衰退制动效能进行虚拟测试，其测试过程包括以下步骤：
[0068] 301、初始参数设定与存储：根据I型制动性能热衰退试验方法对被测试两轴汽车重复制动之后的热衰退制动效能进行虚拟测试的试验工况，确定对被测试两轴汽车进行虚拟试验的试验参数，包括被测试制动器所处车辆的初始速度V2并同步存入存储单元内；
[0069] 302、利用步骤202中所建立的两轴汽车动力学模型，计算得出在初始速度为V2且紧急停车状态下所述两轴汽车的汽车制动距离和制动减速度；
[0070] 步骤302中进行虚拟测试过程中，通过与所述运算处理器相接的显示器输出所述存储单元内所存储的所有计算数据，汽车热衰退制动效能测试过程结束。
[0071] 上述两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征是：步骤1021中所述的多场耦合分析软件为ADINA软件。
[0072] 上述两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，其特征是：步骤103中所述的常规制动器温升试验方法为采用底盘综合试验台对制动器后桥总成进行的定速温升试验方法或者采用制动器试验台上对所述制动器总成进行定速温升试验方法。
[0073] 本发明与现有技术相比具有以下优点：
[0074] 1、设计合理、操作简便且实现方便，明显简化了两轴汽车行车制动系I型制动性能热衰退试验的试验过程，大大了缩短试验周期，大幅减小试验成本。
[0075] 2、计算精度较高。
[0076] 3、实用价值高，能简便投入工业生产应用中，可以在制动器生产过程中很方便地用于进行制动性能的评估与监控，可经济、高效地对制动器的设计、生产、鉴定、销售等环节进行实时跟踪测试，不仅保证产品质量，还可及时调整产品各项性能指标，增强产品生产的柔性。
[0077] 4、推广应用前景广泛，符合当前我国两轴汽车制动系I型制动性能热衰退虚拟测试行业对测试方法所要求的最有效、最精准、最节约的目标要求。
[0078] 5、适用面广，可方便地用于对汽车制动热衰退和制动器摩擦性能热衰退的单因素或多因素分析、极限工况分析等方面的虚拟试验检测研究。
[0079] 综上所述，本发明设计合理、试验投资成本低、模拟效果好且使用操作简便，能有效解决现有制动系I型制动性能热衰退虚拟测量系统及方法所存在的投资成本高、模拟性较差、试验方法复杂、试验周期长等多种缺陷和不足。
[0080] 下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0081] 图1为本发明的试验方法流程图。
具体实施方式
[0082] 如图1所示的一种两轴汽车行车制动系I型制动性能热衰退虚拟试验方法，包括以下步骤：
[0083] 步骤一、确定被测试两轴汽车制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)，被测试制动器为鼓式制动器或盘式制动器，其确定过程如下：
[0084] 101、建立三维几何模型：根据被测试制动器的装配结构，且采用CAD软件建立被测试制动器的三维几何模型；
[0085] 102、建立三维热-机耦合有限元模型，其建模过程如下：
[0086] 1021、将步骤一中所述三维几何模型导入多场耦合分析软件中；本实施例中，所述多场耦合分析软件为ADINA软件；
[0087] 1022、在多场耦合分析软件环境下，选定分析单元并以所选定的分析单元对所述三维几何模型进行网格划分，所选定的分析单元包括被测试制动器的三维实体结构单元、摩擦接触单元和对流辐射单元三种类型单元；
[0088] 1023、输入根据工程材料手册(具体为由马之庚和任陵柏主编的《现代工程材料手册》一书，也可参考常规使用的相关材料手册)所查询出的或者通过常规材料性能测试试验所测得的被测试制动器的材料物性参数，选择一种收敛快速的摩擦接触算法，同时选定摩擦系数-温度模型，最终建立三维热-机耦合有限元模型，所述摩擦系数-温度模型为被测试制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)，式中μ为被测试制动器制动摩擦副的摩擦系数，T为被测试制动器制动摩擦副的摩擦接触面温度；所述材料物性参数包括被测试制动器制动摩擦副所采用材料的弹性模量、密度、导热系数、比热容、热膨胀系数和泊松比；
[0089] 103、对所建立的三维热-机耦合有限元模型进行耦合温度场仿真分析，在被测试制动器制动摩擦副上选定一个或多个测试点，将对所述三维热-机耦合有限元模型在各测试点上的温度场仿真结果与采用常规制动器温升试验方法对各测试点上所测试的温度测试结果进行比较，从而对步骤3中所述函数关系μ＝f(T)进行验证，实现对函数关系μ＝f(T)中的系数进行准确辨识，最终得出被测试制动器制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)；
[0090] 本实施例中，所述常规制动器温升试验方法为采用底盘综合试验台对制动器后桥总成进行的定速温升试验方法或者采用制动器试验台上对所述制动器总成进行定速温升试验方法。
[0091] 本实施例中，被测试制动器所处车辆为两轴汽车，且被测试制动器包括前轴制动器和后轴制动器，则步骤一中需分别对前轴制动器和后轴制动器的制动摩擦副的摩擦系数-温度之间的函数关系μ＝f(T)分别进行确定。
[0092] 步骤二、利用运算处理器且根据I型制动性能热衰退试验方法对被测试两轴汽车各轴制动器的温升及摩擦性能进行虚拟测试，其测试过程包括以下步骤：
[0093] 201、初始参数设定、计算与存储：根据汽车I型制动性能热衰退试验方法规定的进行重复制动的试验工况，确定对被测试制动器进行虚拟试验的试验参数，包括被测试制动器所处车辆的初始速度V0且此时V0＝常量v、汽车制动减速度a0、被测试制动器的初始温度T0、循环制动次数和制动周期，并且根据所确定的循环制动次数(即重复制动次数)和制动周期对总循环次数m和循环周期t进行设定；之后，调用参数设置单元输入所述试验参数和汽车及制动器的结构参数，同时根据测试精度及效率需求设定时间步长Δt，并计算出每一次循环制动过程的时间总步数k＝t/Δt，且将所输入的V0、a0、m、t、Δt、k、ω0和T0均同步存入存储单元内；同时调用参数计算单元并根据公式 计算得出被测试制动器中旋转部件的初始角速度ω0，并将初始角速度ω0分别转存为ωf和ωr后同步存入存储单元内；被测试制动器所处车辆为两轴汽车，则所述被测试制动器包括对两轴汽车的前轮进行制动的前轴制动器和对两轴汽车的后轮进行制动的后轴制动器；所述旋转部件为所述鼓式制动器的制动鼓或盘式制动器的制动盘。
[0094] 实际使用过程中，可以根据具体需要，通过调用参数设置单元输入所述试验参数和汽车及制动器的结构参数。
[0095] 202、建模：所述运算处理器根据汽车五自由度制动动力学方程组：
[0096] 建立被测试制动器所处两轴汽车的动力学模型，所述两轴汽车的前轮和后轮的半径相同，式中∑FX、∑FZ、∑MY、∑Mwf和∑Mwr分别为汽车纵向合力、悬挂质量垂向合力、悬挂质量俯仰合力矩、前轮合力矩和后轮合力矩；Ffb和Frb为根据轮胎模型计算出来的的汽车前、后轮纵向地面作用力；Fe、Ffk、Ffd、Frk、Frd、Mbf和Mbr分别为发动机制动力(其大小与被检测汽车的车速、档位等参数有关)、前悬架弹性元件作用力、前悬架阻尼元件作用力、后悬架弹性元件作用力、后悬架阻尼元件作用力、前轴制动器输出力矩和后轴制动器输出力矩，且上述变量均为在预先设置初始值的基础上根据制动动力学方程组动态解算；M、Ms、Is、Iwf和Iwr分别为汽车质量、悬挂质量、悬挂质量绕质心转动惯量、前轮转动惯量和后轮转动惯量且在步骤201中进行设定；a、b和r分别为汽车质心到前轴距离、质心到后轴距离和车轮的滚动半径且在步骤201中进行设定；
aX、asZ、εsY、εwf和εwr分别为汽车纵向加速度、悬挂质量垂向加速度、悬挂质量俯仰角加速度、前轮角加速度和后轮角加速度，且上述变量根据制动动力学方程组动态解算。
[0097] 建立所述两轴汽车动力学模型的同时，将上述参数分别对应存入存储单元内，且此时
[0098] 实际对转动惯量I计算时，对于两轴车(含双后轴的三轴车)，其前轴制动器和驱动轴制动器所承受的转动惯量I分别为：
[0099]若为双后轴的三轴车，其后制动器应承受的转
动惯量I为
[0100]
[0101] 式中If、Ir、Idr-分别为前、后轴制动器和双后轴的三轴车后制动器应承受的转动2
惯量，kg·m·s ；
[0102] β-前后轴制动力比，对于双后轴的三轴车，两后轴制动力之和作为后轴制动力；
[0103] G0、Ga、r、g-分别为汽车空车总质量(kg)、汽车满载总质量(kg)、车轮滚动半径2
(m)、重力加速度(m/s)。
[0104] 203、被测试两轴汽车制动热衰退循环过程中持续制动时制动器的管路压力值确定，其确定过程如下：
[0105] 2031、调用管路压力值确定模块进行第一个时间步长Δt内制动器管路压力值的确定，其确定过程如下：
[0106] 20311、先按照管路压力非线性递增规律确定前轴制动器和后轴制动器的管路压力递增量Δpf和Δpr，再根据公式pif＝p0f+Δpf和pir＝p0r+Δpr，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的管路压力pif和pir，并将pif和pir分别转存为pf和pr，对存储单元内的管路压力参数进行实时更新，其中p0f和p0r为被测试两轴汽车前轴制动器和后轴制动器的推出压耗；
[0107] 20312、根据步骤20311中计算得出的pf和pr，计算得出前轴制动器和后轴制动器的输出力矩Mbf和Mbr(具体按照制动器输出力矩与管路压力之间的函数关系进行计算)，并对存储单元内的输出力矩参数进行实时更新；
[0108] 20313、结合步骤202中所述的发动机制动力Fe以及汽车前、后轮纵向地面作用力Ffb和Frb，并根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf和计算得出被测试两轴汽车前后轮的角减速度εwf和εwr，并对存储单元内的角减速度参数进行实时更新；
[0109] 20314、调用角速度计算模块且根据公式ωif＝ωf-εwf×Δt和ωir＝ωr-εwr×Δt，计算得出此时间步长Δt结束时前轴制动器或后轴制动器的角速度ωif和ωir，并转存为ωf和ωr以对存储单元内的角速度参数进行实时更新；
[0110] 20315、根据公式∑FX＝-2(Ffb+Frb)-Fe＝-MaX，计算得出此时间步长Δt结束时被测试两轴汽车的减速度aiX并转存为aX；再根据公式Vi0＝V0-aX·Δt计算得出此时间步长Δt结束时被测试两轴汽车的车速Vi0并转存为V0，以对存储单元内的车速参数进行实时更新；
[0111] 20316、结 合 ωf、ωr 和 V0，并 根 据 公 式 和计算得出前轮和后轮的滑移率sif和sir并分别转存为sf和sr以对存
储单元内的滑移率参数进行实时更新；
[0112] 20317、根据轮胎模型且结合sf和sr，计算得出此时间步长Δt结束时汽车前、后轮纵向地面作用力Ffb和Frb并同步存入存储单元内，且相应对步骤202中所建立的动力学模型进行实时更新；
[0113] 20318、根据车速V0并结合经常规试验测试得出的发动机制动力Fe与车速V0之间的函数关系，计算得出此时间步长Δt结束时发动机制动力Fe并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
[0114] 20319、判断|aX-3|是否小于δ，其中δ＝0.01～0.001：当|aX-3|＜δ时，将步骤20311中计算得出的pf和pr输出并另存为pmf和pmr，并将pmf和pmr作为制动热衰退循环过程中持续制动时前轴和后轴制动器的管路压力值；否则，进入步骤2032；
[0115] 2032、调用管路压力值确定模块进行下一时间步长Δt内制动器管路压力值的确定，其确定过程如下：
[0116] 20321、先按照管路压力非线性递增规律确定前轴制动器和后轴制动器的管路压力递增量Δpf和Δpr，再根据公式pif＝pf+Δpf和pir＝pr+Δpr，计算得出此时间步长Δt结束时前轴制动器和后轴制动器的管路压力pif和pir，并将pif和pir分别转存为pf和pr，以对存储单元内的管路压力参数进行实时更新；
[0117] 20322、再根据步骤20321中计算得出的pf和pr，计算得出前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并对存储单元内的输出力矩参数进行实时更新；
[0118] 20323、结合步骤202中所述的发动机制动力Fe以及汽车前、后轮纵向地面作用力Ffb和Frb，并根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf和计算得出被测试两轴汽车前后轮的角减速度εwf和εwr，并对存储单元内的角减速度参数进行实时更新；
[0119] 20324、调用角速度计算模块且根据公式ωif＝ωf-εwf×Δt和ωir＝ωr-εwr×Δt，计算得出此时间步长Δt结束时前轴制动器或后轴制动器的角速度ωif和ωir，并转存为ωf和ωr以对存储单元内的角速度参数进行实时更新；
[0120] 20325、根据公式∑FX＝-2(Ffb+Frb)-Fe＝-MaX，计算得出此时间步长Δt结束时被测试两轴汽车的减速度aiX并转存为aX；再根据公式Vi0＝V0-aX·Δt计算得出此时间步长Δt结束时被测试两轴汽车的车速Vi0并转存为V0，以对存储单元内的车速参数进行实时更新；
[0121] 20326、结 合 ωf、ωr 和 V0，并 根 据 公 式 和计算得出前轮和后轮的滑移率sif和sif，并分别转存为sf和sr，以对
存储单元内的滑移率参数进行实时更新；
[0122] 20327、根据轮胎模型且结合sf和sr计算得出此时间步长Δt结束时汽车前、后轮纵向地面作用力Ffb和Frb并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
[0123] 20328、根据车速V0并结合经常规试验测试得出的发动机制动力Fe与车速V0之间的函数关系，计算得出此时间步长Δt结束时发动机制动力Fe并同步存入存储单元内，并相应对步骤202中所建立的动力学模型进行实时更新；
[0124] 20329、判断|aX-3|是否小于δ，式中δ＝0.01～0.001：当|aX-3|＜δ时，将步骤20321中计算得出的pf和pr输出并另存为pmf和pmr，并将pmf和pmr作为制动热衰退循环过程中持续制动时前轴和后轴制动器的管路压力值；否则，返回步骤20321。
[0125] 具体而言：根据制动器结构参数和摩擦系数μ，可求得制动器的效能因数BF：不同型式的制动器，计算公式也不同。对于常见典型的钳盘式制动器，制动器效能因数BF为：
[0126] BF＝2μ，式中μ——盘与制动衬块间的摩擦系数。
[0127] 支点固定的领从蹄鼓式制动器，领蹄、从蹄的效能因数BF分别为
[0128]
[0129]
[0130] 若两蹄张开力F1、F2相同，则其制动器效能因数BF为
[0131] BF＝BF1+BF2
[0132] α——最大压力线与包角平分线夹角；
[0133] l0——压力中心圆直径，
[0134] γ——摩擦角，γ＝arctanμ，μ为摩擦副摩擦系数；
[0135] λ＝γ+β-α；λ′＝γ-β+α；
[0136]
[0137] 制动管路压力p的计算过程如下：对于气压制动系统，可利用制动器输出力矩表达式Mb＝(p-p0)·Ac·ηm·BF·ρ·r′(即管路压力p，由Mb、BF和p间关系表达式)来计算得到，式中：
[0138] p、p0——管路压力、推出压耗；
[0139] ρ、Ac、ηm、r′、BF——制动气室至制动器杠杆比、制动气室面积、制动气室至制动器驱动机构机械效率、制动鼓或制动盘的作用半径、制动器效能因数；
[0140] 对于液压制动系统，根据管路压力p，由Mb、BF和p间关系表达式，即其制动器输出力矩为：
[0141] Mb＝(p-p0)·Awc·ηa·BF·r′(此公式为制动器输出力矩与管路压力之间的函数关系)
[0142] 式中
[0143] Awc、ηa——制动轮缸工作面积、制动驱动机构中制动主缸之后机械效率，r′为被测试制动器的制动鼓或制动盘的作用半径。
[0144] 204、循环制动：按照步骤201所设定的参数与步骤202中所建立的两轴汽车的动力学模型，在保持前轴制动器和后轴制动器的最大管路压力pf＝pmf且pr＝pmr的前提下，分m次分别对所述前轴制动器与后轴制动器的制动热衰退过程进行循环虚拟测试，其循环虚拟测试的测试过程如下：
[0145] 2041、调用循环测试模块进行制动热衰退试验过程中第一个时间步长Δt内的试验，其试验过程如下：
[0146] 20411、按照步骤20311至步骤20318，分别计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir、汽车前、后轮纵向地面作用力Ffb和Frb以及发动机制动力Fe，并将前轴制动器和后轴制动器的角速度ωif和ωir分别转存为ωf和ωr后，相应对步骤202中所建立的动力学模型进行实时更新；
[0147] 20412、制动器温升计算，调用制动器温升计算模块对此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)分别进行计算，其计算过程如下：
[0148] 20412a、根据公式pbf＝Mbfωf，计算出此时间步长Δt内的前轴制动器和和后轴制动器的吸收制动功率pbf和pbr，并将此时所求得的前轴制动器和后轴制动器的吸收制动功率pbf和pbr同步存入所述存储单元内，对所述存储单元内的吸收制动功率参数进行实时更新；
[0149] 20412b、根据温度变化微分方程Pbf×dt-h×A[T(t)-TE]dt＝m×c×dT(t)和Pbr×dt-h×A[T(t)-TE ]dt＝m×c×dT(t)，分别求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)，并将此时所求得的前轴制动器和后轴制动器的瞬时动态温度T(t)同步存入所述存储单元内，对所述存储单元内的前轴制动器和后轴制动器的瞬时动态温度参数进行实时更新；式中，h为前轴制动器或后轴制动器的对流换热系数，m为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热，TE为环境温度；
[0150] 20413、热衰退摩擦性能计算：调用热衰退摩擦性能计算单元，根据步骤103中验证的函数关系μ＝f(T)且结合步骤20412b中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，相应分别计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr并相应存入存储单元内；
[0151] 20414、调用所述参数计算单元，且根据制动器效能因数的常规计算方法且结合步骤20413中计算得出的热态摩擦系数μf和μr，计算得出此时间步长Δt内前轴制动器和后轴制动器的效能因数BFf和BFr；再根据制动器输出力矩的常规计算方法，且结合此时所计算出来的效能因数BFf和BFr与前轴制动器管路压力pf和后轴制动器管路压力pr，相应推算得出此时间步长Δt结束时，前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并将所计算出来的效能因数BFf和BFr以及输出力矩Mbf和Mbr同步存入所述存储单元内，对所述存储单元内的制动器效能因数和输出力矩参数进行实时更新；同时，相应对步骤203中所建立的两轴汽车动力学模型进行实时更新；
[0152] 20415、利用步骤20414中实时更新后的两轴汽车动力学模型，且根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf，调用参数计算模块计算得出此时间步长Δt结束时的前轮角加速度εwf和后轮角加速度εwr，并将此时所计算出来的前轮角加速度εwf和后轮角加速度εwr同步存入所述存储单元内，对所述存储单元内的前轮角加速度参数进行实时更新；
[0153] 2042、调用循环测试模块进行制动测试过程中下一时间步长Δt内的测试，其测试过程如下：
[0154] 20421、按照步骤20321至步骤20328，分别计算得出此时间步长Δt结束时前轴制动器和后轴制动器的角速度ωif和ωir、汽车前、后轮纵向地面作用力Ffb和Frb以及发动机制动力Fe，并将前轴制动器和后轴制动器的角速度ωif和ωir转存为ωf和ωr，并相应对步骤202中所建立的动力学模型进行实时更新；
[0155] 20422、制动器温升计算，调用制动器温升计算模块对此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)分别进行计算，其计算过程如下：
[0156] 20422a、根据公式pbf＝Mbfωf，计算出此时间步长Δt内的前轴制动器和和后轴制动器的吸收制动功率pbf和pbr，并将此时所求得的前轴制动器和后轴制动器的吸收制动功率pbf和pbr同步存入所述存储单元内，对所述存储单元内的吸收制动功率参数进行实时更新；
[0157] 20422b、根据温度变化微分方程Pbf×dt-h×A[T(t)-TE]dt＝m×c×dT(t)和Pbr×dt-h×A[T(t)-TE]dt＝m×c×dT(t)，分别求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)，并将此时所求得的前轴制动器和后轴制动器的瞬时动态温度T(t)同步存入所述存储单元内，对所述存储单元内的前轴制动器和后轴制动器的瞬时动态温度参数进行实时更新；式中，h为前轴制动器或后轴制动器的对流换热系数，m为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热，TE为环境温度；
[0158] 20423、热衰退摩擦性能计算：调用热衰退摩擦性能计算单元，根据步骤103中验证的函数关系μ＝f(T)且结合步骤20412b中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，相应分别计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr并相应存入存储单元内；
[0159] 20424、调用所述参数计算单元，且根据制动器效能因数的常规计算方法且结合步骤20423中计算得出的热态摩擦系数μf和μr，计算得出此时间步长Δt内前轴制动器和后轴制动器的效能因数BFf和BFr；再根据制动器输出力矩的常规计算方法，且结合此时所计算出来的效能因数BFf和BFr与前轴制动器管路压力pf和后轴制动器管路压力pr，相应推算得出此时间步长Δt结束时，前轴制动器和后轴制动器的输出力矩Mbf和Mbr，并将所计算出来的效能因数BFf和BFr以及输出力矩Mbf和Mbr同步存入所述存储单元内，对所述存储单元内的制动器效能因数和输出力矩参数进行实时更新；同时，相应对步骤203中所建立的两轴汽车动力学模型进行实时更新；
[0160] 20425、利用步骤20424中实时更新后的两轴汽车动力学模型，且根据公式∑Mwf＝Mbf-Ffb·r＝Iwfεwf，调用参数计算模块计算得出此时间步长Δt结束时的前轮角加速度εwf和后轮角加速度εwr，并将此时所计算出来的前轮角加速度εwf和后轮角加速度εwr同步存入所述存储单元内，对所述存储单元内的前轮角加速度参数和后轮角加速度参数进行实时更新；
[0161] 20426、本时间步长Δt结束时，对当前所述存储单元内更新后的被测试两轴汽车的车速V0进行判断：当 时，返回步骤2042，继续进行制动试验过程中的下一时间步长Δt内的试验测试，如此不断循环；当 时，说明已经完成本次热衰退试验虚拟测试过程中的制动试验过程，则进入步骤2043，此时本次热衰退试验中制动试验过程所用时间步数为n1；
[0162] 2043、制动器提速散热阶段的制动器温度计算，其计算过程如下：
[0163] 20431、调用制动器提速散热模块进行循环提速散热计算过程中第一个时间步长Δt内的计算，其计算过程如下：
[0164] 20431a、根据公式-h×A[T(t)-TE]dt＝m×c×dT(t)，求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)并同步存入所述存储单元内，且对所述存储单元内的温度参数进行实时更新，式中TE为环境温度，h为前轴制动器或后轴制动器的对流换热系数，mbf为前轴制动器或后轴制动器的旋转部件的质量，c为前轴制动器或后轴制动器的旋转部件所采用材料的比热；
[0165] 20431b、提速散热摩擦性能计算：根据步骤103中验证的函数关系μ＝f(T)且结合步骤20431a中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr，并相应存入存储单元内；
[0166] 20432、调用制动器提速散热模块进行循环提速散热计算过程中下一个时间步长Δt内的计算，其计算过程如下：
[0167] 20432a、根据公式-h×A[T(t)-TE]dt＝m×c×dT(t)，求解出此时间步长Δt内前轴制动器和后轴制动器的瞬时动态温度T(t)并同步存入所述存储单元内，且对所述存储单元内的温度参数进行实时更新；
[0168] 20432b、提速散热摩擦性能计算：根据步骤103中验证的函数关系μ＝f(T)且结合步骤20432a中计算出来的前轴制动器和后轴制动器的瞬时动态温度T(t)，计算出前轴制动器和后轴制动器在此时间步长Δt内的热态摩擦系数μf和μr，并相应存入存储单元内；
[0169] 20433、返回步骤20432，继续进行循环提速散热计算过程中下一个时间步长Δt的计算，不断重复且经m1＝(t-n1×Δt)/Δt个时间步长计算后进入步骤20434，且本次热衰退试验中提速散热计算过程所用时间步数为m1；
[0170] 2044、重复步骤2041至步骤2043，直至完成m次热衰退性能循环虚拟测试，则完成所述被测试两轴汽车前轴制动器和后轴制动器的温升及摩擦性能虚拟测试过程；
[0171] 步骤二进行虚拟测试过程中，通过与所述运算处理器相接的显示器输出所述存储单元内所存储的所有计算数据进行同步显示，此时被测试两轴汽车各轴制动器的温升及摩擦性能虚拟测试过程结束。
[0172] 本实施例中，步骤2044中完成所述被测试两轴汽车前轴制动器的温升及摩擦性能虚拟测试过程后，还需利用运算处理器且根据I型制动性能热衰退试验方法对被测试两轴汽车重复制动之后的热衰退制动效能进行虚拟测试，其测试过程包括以下步骤：
[0173] 301、初始参数设定与存储：根据I型制动性能热衰退试验方法对被测试两轴汽车重复制动之后的热衰退制动效能进行虚拟测试的试验工况，确定对被测试两轴汽车进行虚拟试验的试验参数，包括被测试制动器所处车辆的初始速度V2并同步存入存储单元内；
[0174] 302、利用步骤202中所建立的两轴汽车动力学模型，计算得出在初始速度为V2且紧急停车状态下所述两轴汽车的汽车制动距离和制动减速度，也就是说，本步骤中需对所述两轴汽车动力学模型的动力学参数进行修改；
[0175] 步骤302中进行虚拟测试过程中，通过与所述运算处理器相接的显示器输出所述存储单元内所存储的所有计算数据，汽车热衰退制动效能测试过程结束。
[0176] 以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。
