1、1A 题 制动器试验台的控制方法分析汽车的行车制动器（以下简称制动器）联接在车轮上，它的作用是在行驶时使车辆减速或者停止。制动器的设计是车辆设计中最重要的环节之一，直接影响着人身和车辆的安全。为了检验设计的优劣，必须进行相应的测试。在道路上测试实际车辆制动器的过程称为路试，其方法为：车辆在指定路面上加速到指定的速度；断开发动机的输出，让车辆依惯性继续运动；以恒定的力踏下制动踏板，使车辆完全停止下来或车速降到某数值以下；在这一过程中，检测制动减速度等指标。假设路试时轮胎与地面的摩擦力为无穷大，因此轮胎与地面无滑动。为了检测制动器的综合性能，需要在各种不同情况下进行大量路试。但是，车辆设计阶段无法
2、路试，只能在专门的制动器试验台上对所设计的路试进行模拟试验。模拟试验的原则是试验台上制动器的制动过程与路试车辆上制动器的制动过程尽可能一致。通常试验台仅安装、试验单轮制动器，而不是同时试验全车所有车轮的制动器。制动器试验台一般由安装了飞轮组的主轴、驱动主轴旋转的电动机、底座、施加制动的辅助装置以及测量和控制系统等组成。被试验的制动器安装在主轴的一端，当制动器工作时会使主轴减速。试验台工作时，电动机拖动主轴和飞轮旋转，达到与设定的车速相当的转速(模拟实验中，可认为主轴的角速度与车轮的角速度始终一致)后电动机断电同时施加制动，当满足设定的结束条件时就称为完成一次制2动。路试车辆的指定车轮在制动时承
3、受载荷。将这个载荷在车辆平动时具有的能量（忽略车轮自身转动具有的能量）等效地转化为试验台上飞轮和主轴等机构转动时具有的能量，与此能量相应的转动惯量(以下转动惯量简称为惯量)在本题中称为等效的转动惯量。试验台上的主轴等不可拆卸机构的惯量称为基础惯量。飞轮组由若干个飞轮组成，使用时根据需要选择几个飞轮固定到主轴上，这些飞轮的惯量之和再加上基础惯量称为机械惯量。例如，假设有 4 个飞轮，其单个惯量分别是：10、20、40、80 kgm2，基础惯量为 10 kgm2，则可以组成10，20，30，160 kgm2 的 16 种数值的机械惯量。但对于等效的转动惯量为 45.7 kgm2 的情况，就不能精确
4、地用机械惯量模拟试验。这个问题的一种解决方法是：把机械惯量设定为40 kgm2，然后在制动过程中，让电动机在一定规律的电流控制下参与工作，补偿由于机械惯量不足而缺少的能量，从而满足模拟试验的原则。一般假设试验台采用的电动机的驱动电流与其产生的扭矩成正比（本题中比例系数取为 1.5 A/Nm） ；且试验台工作时主轴的瞬时转速与瞬时扭矩是可观测的离散量。由于制动器性能的复杂性，电动机驱动电流与时间之间的精确关系是很难得到的。工程实际中常用的计算机控制方法是：把整个制动时间离散化为许多小的时间段，比如 10 ms 为一段，3然后根据前面时间段观测到的瞬时转速与/或瞬时扭矩，设计出本时段驱动电流的值，
5、这个过程逐次进行，直至完成制动。评价控制方法优劣的一个重要数量指标是能量误差的大小，本题中的能量误差是指所设计的路试时的制动器与相对应的实验台上制动器在制动过程中消耗的能量之差。通常不考虑观测误差、随机误差和连续问题离散化所产生的误差。现在要求你们解答以下问题：1. 设车辆单个前轮的滚动半径为 0.286 m，制动时承受的载荷为 6230 N，求等效的转动惯量。2. 飞轮组由 3 个外直径 1 m、内直径 0.2 m 的环形钢制飞轮组成，厚度分别为 0.0392 m、0.0784 m、0.1568 m，钢材密度为7810 kg/m3，基础惯量为 10 kgm2，问可以组成哪些机械惯量？设电动机
6、能补偿的能量相应的惯量的范围为 -30, 30 kgm2，对于问题 1 中得到的等效的转动惯量，需要用电动机补偿多大的惯量？、3. 建立电动机驱动电流依赖于可观测量的数学模型。在问题 1 和问题 2 的条件下，假设制动减速度为常数，初始速度为 50 km/h，制动 5.0 秒后车速为零，计算驱动电流。4. 对于与所设计的路试等效的转动惯量为 48 kgm2，机械惯量为 35 kgm2，主轴初转速为 514 转/分钟，末转速为 257 转/ 分钟，时间步长为 10 ms 的情况，用某种控制方法试验得到的数据见附表。请对该方法执行的结果进行评价。5. 按照第 3 问导出的数学模型，给出根据前一个时
7、间段观测到的瞬时转速与/或瞬时扭矩，设计本时间段电流值的计算机控制方法，并对该方法进行评价。6. 第 5 问给出的控制方法是否有不足之处？如果有，请重新设计一个尽量完善的计算机控制方法，并作评价。4制动器试验台的控制方法模型摘要制动器试验台用于模拟理想的路试环境下汽车的制动过程，从而检验制动器的性能。试验台上用飞轮组的转动惯量模拟路试中汽车的惯性，使得飞轮组的转动动能与汽车的平动动能相当。由于飞轮组只能组合出离散的转动惯量，无法精确地组合出等效的转动惯量，所以我们采用对试验台系统提供驱动电流的办法补偿由于转动惯量偏差带来的动能的偏差，使得试验中制动器做功模拟路试中制动器做功，试验中飞轮组转速模
8、拟路试中车轮转速。为了简化问题，我们做出了一些必要的假设，将路试中汽车的制动过程归结为匀减速过程，由于试验中轮轴的转速与制动扭矩都可以读出，我们假设制动扭矩和轮轴实际转速为已知函数。试验中，通过驱动电流产生驱动扭矩，使得在任一时间微元内，制动器对系统做功相当于对等效转动惯量的系统做功，同时，任意时刻系统满足动力学方程，而初始时刻的转速可以直接读出，是已知的，于是可以建立驱动扭矩控制飞轮组系统的微分方程模型。把该模型用于题目提供的一组试验数据，可以看出，模型与实际情况吻合很好。利用微分方程模型，可以解出驱动扭矩关于制动扭矩的函数关系，以及驱动扭矩关于轮轴角速度变化率的函数关系。试验中制动扭矩与轮
9、轴转速都可以直接读出，于是，我们可以利用前一个时间段观测出的制动扭矩或者轮轴转速，计算出本时间段的控制扭矩以及控制电流的大小。但是，通过对题目提供的一组试验数据的分析可以看出，在试验中，不可避免地会出现不确定因素的干扰，导致轮轴转速偏离理论值，进而导致整个制动过程的能量误差增大。为了尽量减少不确定因素的影响，我们通过轮轴转速的读数与理论值的比较，时刻监测系统的运转状态，调整驱动电流的大小，补偿下一个时间微元内制动扭矩实际做功与理论值之间的偏差，修正轮轴转速，从而建立新的微分方程模型，它在一定程度上能针对不确定因素带来的干扰进行调整，使得系统具有较强的抗干扰能力。由于缺乏实地试验的条件，我们利5
10、用动力学方程，用计算机模拟一次制动试验，对新的模型进行检验。我们从制动器做功模拟、飞轮组转速模拟、以及控制系统的抗干扰能力三个方面来理解所谓的“完善的控制方法” 。 由于不同的试验中飞轮组的动能变化不同，能量误差的大小不能反映控制方法的优劣，我们引进了相对能量误差，使得能量误差在实际评价过程中有着更加客观的价值。通过对题目提供的试验数据中相对能量误差和模拟实验的相对能量误差的比较可以看出，新的模型的确可以减小相对能量误差，提高控制台试验模拟的精度。最后，我们对两微分方程模型各自的优缺点进行了评价。关键词：微分方程、能量误差、能量补偿、转速补偿、图像对比、不确定因素一、 问题的提出在汽车设计阶段
11、中，通常利用制动器试验台代替路试来检验制动器的性能。模拟试验的原则是试验台上制动器的制动过程与路试车辆上制动器的制动过程尽可能一致。现假设轮胎与地面无滑动。将路试时车辆的平动动能等效地转化为试验台上飞轮和主轴等机构的转动动能，与此相应的转动惯量（忽略车轮自身转动具有的能量）称为等效的转动惯量（以下简称惯量） 。试验台上不可拆卸机构的惯量称为基础惯量。飞轮组由若干个飞轮组成，飞轮组合出的惯量加上基础惯量称为机械惯量。由于飞轮组的惯量为固定的几个值，故为了精确模拟等效惯量，人们通常让电动机在一定规律的电流控制下参与工作，补偿由于机械惯量不足而缺少的能量。 假设试验台采用的电动机的驱动电流与其产生的
12、扭矩成正比；且试验台工作时主轴的瞬时转速与瞬时扭矩是可观测的离散量。由于制动器性能的复杂性，电动机驱动电流与时间之间的精确关系是很难得到的。我们把把整个制动时间离散化为许多小的时间段，然后根据前面时间段观测到的瞬时转速与瞬时扭矩，设计出本时段驱动电流的值，这个过程逐次进行，直至完成制动。评价控制方法优劣的一个重要数量指标是能量误差的大小。计算时我们不考虑观测误差、随机误差和连续问题离散化所产生的误差。 现在我们要解决的问题如下：1. 设车辆单个前轮的滚动半径为 0.286 m，制动时承受的载荷为 6230 N，求得等效的转动惯量。2. 飞轮组由 3 个外直径 1 m、内直径 0.2 m 的环形
13、钢制飞轮组成，厚度分别为 0.0392 m、0.0784 m、0.1568 m，钢材密度为 7810 kg/m3，基础惯量为 10 kgm2，问可以组成哪些机械惯量？设电动机能补偿的能量相应的惯量的范围为 -30, 30 kgm2，对于问题 1 中得到的等效的转动惯量，解释需要用电动机补偿惯量的大小。63. 建立电动机驱动电流依赖于可观测量的数学模型。在问题 1 和问题 2 的条件下，假设制动减速度为常数，初始速度为 50 km/h，制动 5.0 秒后车速为零，计算驱动电流。4. 对于与所设计的路试等效的转动惯量为 48 kgm2，机械惯量为 35 kgm2，主轴初转速为 514 转/分钟，末
14、转速为 257 转/ 分钟，时间步长为 10 ms的情况，根据附表给出的数据，对该方法执行的结果进行评价。5. 按照第 3 问导出的数学模型，给出根据前一个时间段观测到的瞬时转速与/或瞬时扭矩，设计本时间段电流值的计算机控制方法，并对该方法进行评价。6. 评价第 5 问给出的控制方法的不足之处。且重新设计一个尽量完善的计算机控制方法，并作评价。二、合理的假设根据制动器试验台的一些具体要求，并为了达到简化问题的目的，除问题中已经给出的假设外，我们进一步作如下假设：1.假设汽车路试时地面的滚动摩擦力很小，在我们讨论的范围内可以忽略不计。2.假设汽车在路试时受到恒定的制动力，因而其在理想路试的状态下
15、作匀减速运动。3.由于制动扭矩 (t)不可控制，但可从试验中直接读出，我们认为它是一0个已知函数，并且不受轮轴转速微小变动的影响。4.试验中轮轴的转速 (t)也可以直接读出，所以作为已知函数。04.忽略试验数据的获得过程的主观误差和仪器误差。5.假设试验过程中没有任何仪器损坏现象。6.假设计算机能够完全实现所需要进行的控制方法。三、符号系统 ： 主轴标准（路试）角速度(rad/s)： 主轴观测（试验）角速度(rad/s)0： 各时刻的角速度(rad/s)： 主轴转速（ r/m） ： 各时刻的转速（r/m）： 等效转动惯量（kg* ）0 2J： 机械转动惯量（ kg* ）2t： 时间（s）： 制
16、动扭矩（N*m）07： 各时刻的制动扭矩（N*m）L： 驱动扭矩（N*m ）I： 驱动电流（A）k： 驱动电流关于驱动扭矩的比例系数（A/Nm）： 控制台模拟过程中制动扭矩做功（J） E1： 理想路试情况下制动扭矩做功（J） E2： 能量误差（J ）： 相对能量误差（J）四、问题的分析在理想的路试中，以恒定的力踏下制动踏板，汽车以恒定的加速度减速。但是在模拟实验中，由于飞轮组无法精确模拟实际状况，必须由驱动电流补偿不足的转动惯量，驱动电流产生驱动扭矩，使得实验室中制动扭矩做功模拟路试时的制动扭矩做功，并且试验室中轮轴转速模拟路试时的实际转速。 】在模拟试验中，任意时间段内制动扭矩做功应该等于路
17、试的理想情况下制动扭矩做的功，据此可以建立微分方程模型。对于题目给出的实验数据，通过作图可以发现，转速数据以及制动扭矩数据在初期都有较大变化，而后转速数据趋于线性，而制动扭矩数据在小范围内波动。两组数据的变化与波动应该与随机因素有关，最终造成能量误差。虽然随机因素的出现是不可预测的，但是，通过计算机监测可以及时发现不确定因素的出现，然后控制驱动电流，在接下来的一段时间内调整转速，补偿由于不确定因素带来的做功模拟以及转速模拟的偏差。据此可以建立针对不确定因素干扰修正系统的微分方程模型。 转动惯量定义(惯性力矩)：定义 1：构件中各质点或质量单元的质量与其到给定轴线的距离平方乘积的总和。 所属学科
18、： 机械工程(一级学科) ；机构学(二级学科) ；机构动力学(三级学科) 定义 2：面积或刚体质量与一轴线位置相关联的量，是面积微元或组成刚体的质量微元到某一指定轴线距离的二次方的乘积之积分。 所属学科： 水利科技(一级学科) ；工程力学、工程结构、建筑材料(二级学科) ；工程力学（水利）(三级学科) 五、模型的建立与解答路试中汽车的惯性可以等效地转化为试验台上飞轮和主轴等机构转动时具有的转动惯量，使得试验台飞轮系统的转动动能与路试时汽车的平动动能相当，8等效转动惯量由飞轮组的转动惯量与驱动电流的补偿共同构成，驱动电流通过产生驱动扭矩，使制动扭矩做功等效于对一个转动惯量等于等效惯量的系统做功。
19、在试验台上，等效惯量 由两部分组成，一部分 是飞轮组和基础惯量之JMJ和的机械惯量,另一部分 由电模拟产生。E问题一：汽车平动所具有的动能为E = m ；12 2飞轮组系统转动时所具有的动能为E= ；120 2两者应该相等，所以m = ；12 2 120 2由=得等效惯量为：= m = = kg*0 2 2 62309.80.2862= 52 kg*2问题二：记圆环的外半径为 R，内半径为 r，厚度为 h，密度为 ，圆环的转动惯量为根据本题中所给物理量，可得飞轮质量的计算式为：（1） 21()4mVrh再根据飞轮即圆筒的转动惯量的公式为 1：（2） )(21rJ（3）J=12（ 44）9R=1
20、 m，r=0.2 m 已知，h 共有三种选择：当 h=0.0392m 时，J=30 kg*2当 h=0.0784m 时，J=60 kg*2当 h=0.1568m 时，J=120 kg*2考虑基础惯量为 10 kgm2，所以可以组成的机械惯量有：10 kg* ，40 kg* ，70 kg* ，100 kg* ，130 kg* ，160 kg* ，2 2 2 2 2 2190 kg* ，220 kg*2 2方案一：组合出 40 kgm2 的机械惯量，并用电动机补偿 12 kgm2 的惯量。方案二：组合出 70 kgm2 的机械惯量，并用电动机补偿 -18 kgm2 的惯量。主要误差是由于电流补偿的
21、那部分，所以由电动机补偿的惯量的大小越小越好！对于问题一中得到的转动惯量 52 kg* ，选择 40 kg* 的机械惯量最为接近，2 2也最为合理，用电动机补偿 12 kg* 的转动惯量。2问题三：假设 t 时刻驱动扭矩为 L(t)，主轴的角速度为 ，经过 dt 时间增加()d ， dt 时间内制动扭矩做功为L(t) dt在试验中制动扭矩做的功！12 (+)2 12 2 ()模拟等效系统的动能增量为：在理想路试上制动扭矩所做的功！120 (+)2 120 2为了模拟等效系统，两者应该相等：L(t) dt = 12 (+)2 12 2 () 120 (+)2 120 2得L(t)=(0)制动扭矩
22、为 ，再根据动量矩定理 1，制动力矩为：0()公式： dtLJ0()()=已知，于是微分方程模型为：(0)=010L(t)= (0)0()()=初始条件 (0)已知 ()在问题 1 和问题 2 的条件下，假设制动减速度为常数，初始速度为 50 km/h，制动 5.0 秒后车速为零，这样计算得：汽车的角速度为：rad/s，56.482.0635rv因为本题中假设制动减速度为常数，所以实验台上的飞轮角减速度为：rad/s2， 所以：的 dw/dt=9.71 rad/271.9tL(t)= (0)=116.55 驱动电流为： I = kL(t) =174.8A问题四：将表格绘制成曲线图：图一：制动扭
23、矩随时间变化曲线11图二：转速随时间变化曲线从图中可以看出，制动扭矩在 01 秒内呈现递增特性，而在 1 秒之后在一个恒定值附近小范围的波动。可以看出，这是控制系统的调节作用。与之对应的，转速曲线在刚开始一段时间内变动并不规律，在 1 秒附近有一次突然地减小，而 1 秒之后则基本符合线性规律。在理想的路试情况下，从初始转速达到末转速的变化过程中车轮角速度应该匀速递减。下图是实验测得的单个车轮的转速数据与理想路试情况下其转速数据的对比。图三：试验转速数据与理想路试情况转速数据的对比图（图中虚线表示理想路试情况下转速随时间变化情况，实线表示实验测得的车轮转速数据）12从图中可见，除去开始的一段车轮
24、转速的不规律变化情况外（01 秒内） ，其余时间实际数据和理想结果之间吻合程度还是比较好的。下面再计算这个控制方法的能量误差：将 4.67 秒的减速时间分成长度为 0.01 秒的时间小段。认为每段时间内制动扭矩为常量，等于该时段初测得的数值。则全过程中制动扭矩做功为：E1=0()()= 49304 J在理想路试情况下中制动器做功等于系统动能变化量，为：E2=120(0)2120()2=52209 所以能量误差为： = E2 E1=2905 为度量能量误差大小，我们引入相对的能量误差，它是能量误差与理想情况下制动器做功的比值：= E2=5.56%于是可以看出此处的能量误差较大。现在我们分析能量误
25、差产生的主要原因。根据问题三中的微分方程模型(*)，可解得：(#)()=(0)1000()考虑到实际情况是离散的，=0101=0本实验中 ， 可直接读出，利用上述方程可以计算0=300=53.86 /出 的理论值， 的理论值与观测值的对比见下图： 13图四： 的理论值与测量值的对比图（虚线为 的理论值，实线为 的测量值） 以 1s 时刻突变后测得的转速为初始转速，计算 1s 后转速的理论值，得到下图：图五从图中可以看出，01s 内理论值与测量值吻合很好，但在 1s 时，测量值突然减小，而后理论值与测量值以相同的速率线性递减，如果以 1s 时的测量值作为初始值计算 1s 后的理论值，则与测量情况
26、完全吻合。 微分方程模型没有考虑摩擦力及不确定因素造成的转速的随机变化，这说明而在 1 秒的时刻，试验中系统受到突然阻力，比如碰到某阻碍物，使得转速突然减小，而其余时间模型14与实际情况吻合很好，这说明摩擦力可以忽略不计，试验中的能量误差主要来源于 1s 时刻的突发阻力，使得转速在其余时间均小于理论值，进而导致制动扭矩做功小于理论数值。综合分析，我们得出的结论是这种控制方法对于转速的模拟情况较好，而应对不确定因素导致的突变时容易产生较大的能量误差，因此需要进行改进。问题五：第一种算法：根据微分方程(*)，得：()=000()所以，得到控制电流为：()=()=000()于是，可以根据某时刻测得的
27、制动扭矩 控制下一个时间段的驱动电流0()I(t)。第二种算法：在 n 到 n+1 所对应的时间内，通过补偿扭矩使得制动器的实际做功相当于对等效系统的做功，据此可以得到方程：t = 12 2 12 2+1+ 120 2 120 2+1解得：( )= 12t(0) 22+1但是问题在于此时 未知，所以我们 利用上一个时间段的数据来近似控+1制本时间段，即把上述公式修改为如下：( )+1= 12t(0) 22+1电流控制式为：( )I+1=+1= 12t(0) 22+1这样，我们同样可以根据某时间段初、末测得的轮轴转动角速度 控1n,制下一个时间段的驱动电流 In+1。对两种算法的评价：第一种算法
28、通过扭矩来对驱动电流进行预测，比较直观方便；而第二种算法由于要利用到仅间隔 0.01 秒的转速数据，会导致误差较大，从原表中的数据可以看出，有时测得的间隔 0.01 秒的两转速甚至相等。从这个角度看，第一种方法较优。？15两种方法的共同缺陷在于均利用第 n 个时间段的驱动扭矩理论数值去控制下一个时间段，这样就产生了滞后误差。问题六：从第四问中可以看出，实际试验过程中常常会有不确定因素出现而导致实际转速偏离理论值，而上题中控制方法没有考虑随机因素的影响，所以还是存在不足。下面讨论对于这种方法的改进。假设某次试验的主轴初始角速度为 ，将问题四解答中（#）式离散化，(0)得角速度的理论值为：()=(
29、0)101=00()但实际试验中由于不确定因素的影响，测得的 ，于是要通过调0()()整驱动扭矩 L(t)补偿做功偏差，使得下一个时间段内制动扭矩做功的偏差等于驱动扭矩补偿的能量：0()()0()0()=()000()0()解得()=0()()0()0此即为修正后的驱动扭矩控制表达式。在 t 时刻成立：()0()=于是，建立转速修正模型：()=0()()0()0()0()=初始条件 (0)已知 ()下面我们对问题四给出的数据验证这种控制方法确实能够减小能量误差。由修正模型解得修正后的角速度为：dt()=(0)100()()对其离散化，得到：()=(0)11=00()()而 )(3)(ftt已经
30、假设制动扭矩 (t)不受轮轴转速微小变动的影响，所以如果初始角速016度 取问题四中的数据，则制动扭矩 (i)也可以取问题四中的相应数据。计)0(0算数据见附表，由于数据量过多，附表中只列了间隔为 0.05 秒的各时刻的相关数据。将这些数据绘制成下面的曲线图，图中虚线为修正后的转速，实线为未修正的原始数据：图六：修正后的转速与原始数据的对比全过程中制动扭矩做功为：1=1=00()()=50184.10 在理想路试情况下中制动器做功为：E2=120(0)2120()2=49099.28 能量误差为： = E2 E1=1984.82 相对的能量误差为：= E2=4.04%误差相对问题四中的相对能量
31、误差 5.56%减小，同时，观察作图可以发现，此控制方法对于突变情况具有良好的反映和“自我修复”的能力，说明此方法相对比较完善。六、模型评价171、从图五可以看出， 式的常微分方程模型可以很好的模拟没有不确定()因素干扰的试验，但出现不确定因素干扰时，模拟偏差较大。2、问题六中建立的修正模型(*)，通过时刻监测主轴转速，调整驱动电流值修正由不确定因素导致的转速偏差，能较好的模拟路试时制动器的做功。3、从图六中可以看出，虽然通过调整驱动电流可以很好的补偿转速偏差和做功偏差，但是修正后的转速曲线斜率出现少量偏差，仍待改进。参考文献：1 姜启源，谢金星，叶俊 数学模型 （第三版） 北京 高等教育出版
32、社 20032 钟锡华，周岳明 力学（大学物理通用教程） 北京 北京大学出版社 20043 楼红卫，林伟 常微分方程 上海 复旦大学出版社 2007附表：模拟试验数据时间(s)L0(t)(N.m)L(t)(N.m)f(t)(rpm)w(t)(rad/s)原始转速(rpm)制动力矩做功(累计)(J)0 40 10.833 514.33 53.862 514.33 21.540.05 45 12.199 513.92 53.82 513.79 134.60.1 57.5 15.716 513.42 53.767 512.15 276.470.15 67.5 18.368 512.81 53.703
33、 512.15 447.73180.2 91.25 24.793 512.05 53.624 511.6 665.70.25 115 31.266 511.05 53.519 510.51 948.240.3 150 40.728 509.78 53.386 509.42 1308.950.35 181.25 49.421 508.17 53.217 507.23 1757.210.4 208.75 56.801 506.26 53.017 505.6 2283.530.45 226.25 61.828 504.12 52.793 502.87 2862.60.5 242.5 66.47 50
34、1.8 52.551 500.14 3485.150.55 241.25 65.745 499.38 52.297 498.5 4125.730.6 262.5 71.091 496.9 52.037 496.86 4787.920.65 266.25 72.454 494.27 51.761 493.58 5476.510.7 268.75 73.781 491.63 51.485 489.76 6164.020.75 277.5 76.192 488.92 51.201 487.03 6870.890.8 267.5 74.675 486.23 50.919 482.12 7561.730
35、.85 282.5 78.561 483.54 50.637 479.93 8262.710.9 277.5 76.816 480.76 50.346 477.75 8971.720.95 276.25 76.474 478.04 50.062 475.02 9661.631 283.75 82.882 475.29 49.773 465.19 10363.811.05 277.5 82.745 472.6 49.492 459.73 11053.311.1 283.75 84.955 469.91 49.211 456.46 11750.621.15 273.75 81.294 467.23
36、 48.929 454.82 12434.151.2 283.75 84.988 464.56 48.65 451 13117.421.25 282.5 84.602 461.85 48.366 448.27 13805.231.3 276.25 82.407 459.18 48.087 446.08 14473.981.35 285 85.378 456.49 47.805 442.81 15152.021.4 276.25 81.705 453.79 47.522 441.71 15821.31.45 280 82.846 451.13 47.244 438.98 16476.571.5
37、281.25 82.487 448.42 46.96 437.35 17140.311.55 273.75 81.002 445.72 46.678 433.52 17791.281.6 283.75 83.993 443.07 46.399 430.79 18438.361.65 281.25 83.604 440.35 46.114 427.52 19091.921.7 275 81.05 437.68 45.836 425.88 19726.561.75 286.25 84.754 434.99 45.554 422.6 20370.461.8 278.75 82.138 432.26
38、45.268 420.42 21012.621.85 278.75 82.92 429.61 44.99 416.6 21633.871.9 281.25 83.676 426.93 44.709 413.87 22261.931.95 275 82.202 424.24 44.427 410.59 22883.832 282.5 83.711 421.6 44.151 408.95 23497.952.05 283.75 84.838 418.86 43.865 405.13 24125.22.1 270 79.994 416.19 43.584 403.49 24729.86192.15
39、282.5 84.519 413.53 43.307 399.67 25336.082.2 282.5 84.897 410.82 43.022 396.4 25949.692.25 275 83.053 408.14 42.741 393.12 26547.512.3 285 85.693 405.49 42.464 390.94 27144.082.35 281.25 84.152 402.77 42.18 388.75 27743.062.4 273.75 81.929 400.11 41.9 386.02 28325.982.45 283.75 85.784 397.45 41.622
40、 382.2 28910.252.5 282.5 85.396 394.75 41.339 379.47 29498.872.55 275 83.576 392.09 41.061 376.19 30069.112.6 282.5 86.335 389.48 40.787 372.92 30636.032.65 285 86.214 386.76 40.502 371.28 31218.372.7 276.25 83.569 384.06 40.22 368.55 31785.062.75 282.5 85.503 381.43 39.944 365.82 32342.812.8 291.25
41、 87.635 378.66 39.654 363.64 32924.472.85 277.5 84.776 375.94 39.37 359.27 33483.662.9 280 85.16 373.32 39.095 357.08 34026.152.95 287.5 87.444 370.64 38.814 354.35 34579.883 278.75 84.32 367.93 38.53 352.17 35127.183.05 275 83.25 365.31 38.256 349.44 35652.793.1 287.5 87.999 362.66 37.979 345.62 36
42、189.853.15 285 85.797 359.93 37.693 344.53 36732.863.2 275 84.62 357.25 37.412 339.61 37257.733.25 282.5 86.508 354.63 37.137 337.43 37777.333.3 288.75 88.402 351.91 36.853 334.7 38308.733.35 281.25 86.575 349.21 36.57 331.42 38829.643.4 276.25 84.157 346.6 36.297 329.78 39330.673.45 286.25 88.222 3
43、43.94 36.019 325.96 39840.563.5 282.5 86.553 341.24 35.736 323.78 40352.313.55 277.5 85.539 338.59 35.458 320.5 40850.283.6 282.5 87.633 335.96 35.183 317.23 41343.493.65 288.75 89.565 333.27 34.902 314.5 41845.953.7 283.75 87.495 330.58 34.619 312.31 42340.883.75 273.75 84.95 327.95 34.344 309.04 4
44、2819.353.8 285 88.497 325.32 34.068 306.31 43299.993.85 286.25 89.41 322.62 33.785 303.03 43786.443.9 280 86.934 319.94 33.505 300.85 44262.543.95 278.75 87.16 317.34 33.232 297.57 44723.484 283.75 88.766 314.7 32.957 294.84 45190.144.05 291.25 91.094 312.01 32.675 292.11 45663.92204.1 281.25 87.408
45、 309.33 32.394 289.93 46125.54.15 275 86.645 306.73 32.122 286.1 46569.474.2 287.5 90.058 304.1 31.846 283.92 47022.064.25 287.5 90.018 301.39 31.562 281.19 47480.974.3 282.5 89.036 298.72 31.282 277.91 47927.954.35 272.5 85.971 296.12 31.011 275.18 48356.614.4 285 90.006 293.53 30.739 272.45 48788.874.45 291.25 92.587 290.85 30.459 269.18 49229.094.5 281.25 89.401 288.17 30.178 266.45 49661.864.55 273.75 87.08 285.56 29.905 263.72 50081.314.6 281.25 89.593 282.99 29.636 260.99 50492.894.65 288.75 92.013 280.35 29.359 258.26 50914.314.67 288.75 91.999 279.27 29.246 257.17 51084.1
