发明内容
本发明所要解决的技术问题是提供一种高速列车车头优化设计方法,能够减少高速列车的阻力,并且可以缩短研制周期,减轻设计人员的劳动强度,实现高速列车车头的快速高效的设计。
为解决上述技术问题,本发明高速列车车头优化设计方法的技术方案是,包括如下步骤:
(1)建立高速列车车头的参数化模型,其中包括描绘列车的头部形状的三条主控制型线和两条辅助控制型线,所述三条主控制型线包括纵向对称面最大控制型线、俯视最大控制型线和车体截面外廓型线,所述两条辅助控制型线使车头背部和侧面光滑连接,所述纵向对称面最大控制型线主要控制列车的阻力,所述俯视最大控制型线主要控制列车的会车压力波,所述车体截面外廓型线保持与现有的高速列车相同;
(2)纵向对称面最大控制型线和俯视最大控制型线用二次曲线来表示,每条二次曲线用3个参数来控制,所述3个参数分别是列车头部长度、角度和形状参数,纵向对称面最大控制型线和俯视最大控制型线共由6个参数控制,根据现有高速列车车头的外形,对这6个参数选取合理的变化范围,并取具体数值,分别对纵向对称面最大控制型线和俯视最大控制型线的参数排列组合,列出数据设计表;
(3)根据所述设计表中的参数,利用计算流体动力学(CFD)软件分别对纵向对称面最大控制型线和俯视最大控制型线的试验方案进行二维流动计算,对纵向对称面最大控制型线的试验方案计算阻力,选取使得阻力系数最小的参数值作为纵向对称面最大控制型线的参数,对俯视最大控制型线的试验方案计算会车压力波的幅值,选取使得会车压力波的幅值最小的参数值作为俯视最大控制型线的参数;
(4)根据所选取的纵向对称面最大控制型线和俯视最大控制型线的参数,结合车体截面外廓型线通过无限插值得到Coons曲面,根据需要增加辅助控制型线数量对Coons曲面进行必要的修改,确定三维高速列车车头模型;
(5)对三维高速列车车头模型开展三维流动计算,计算得到的阻力与现有的高速列车相比,验证阻力是否减小,如果是,则设计完成;如果不是,再回到第(2)步,更改6个参数的变化范围,重复第(3)步到第(5)步。
本发明通过采用上述方法,减少了高速列车的阻力,并且缩短了研制周期,减轻了设计人员的劳动强度,实现高速列车车头的快速高效的设计。
具体实施方式
本发明通过参数建模,建立磁悬浮列车车头的参数化模型,然后采用计算流体力学(CFD)软件进行流动计算,利用正交设计表或者均匀设计表分别对纵向型面和水平型面进行二维的优化。随后根据纵向型面和水平型面通过无限插值得到Coons曲面,按照需要增加辅助控制型线对Coons曲面进行必要的修改,确定三维列车车头模型,最后对三维列车车头模型开展三维流动计算,与现有的磁悬浮列车进行比较,以确认是否实现减阻的目标。具体包括如下步骤:
(1)建立高速列车车头的参数化模型,通过对现有的磁悬浮列车车头的分析,列车的头部形状由三条主控制型线和两条辅助控制型线来控制。三条主控制型线为纵向对称面最大控制型线、俯视最大控制型线和车体截面外廓型线。两条辅助控制型线主要使车头背部和侧面光滑连接。如图1所示,A线为纵向对称面最大控制型线,D线为俯视最大控制型线,E线为车体截面外廓型线,B线和C线为辅助控制型线。所述纵向对称面最大控制型线主要控制列车的阻力,所述俯视最大控制型线主要控制列车的会车压力波,所述车体截面外廓型线保持与现有的高速列车相同。
(2)纵向对称面最大控制型线和俯视最大控制型线用二次曲线来表示,每条二次曲线用3个参数来控制,所述3个参数分别是列车头部长度L、角废θ和形状参数ρ,纵向对称面最大控制型线和俯视最大控制型线共由6个参数控制,根据现有高速列车车头的外形,对这6个参数选取合理的变化范围,并取具体数值,分别对纵向对称面最大控制型线和俯视最大控制型线的参数排列组合,列出数据设计表。如图2所示,其中AB是一条二次曲线,AC和BC分别是A点和B点的切线,M是AB的中点,S是MC和AB的交点,L为列车头部的长度,∠ACM为角度θ,形状参数ρ的定义为SM线段与CS线段之比,即ρ=MS/MC。二次曲线由预定的起点A和终点B以及这些点处所希望的切线方向(C点)决定;并通过二次曲线的形状参数ρ控制曲线形状:ρ大于0.5、等于0.5和小于0.5分别对应着双曲线、抛物线和椭圆。特别是当ρ→1时,曲线趋于折线ACB;当ρ→0时,曲线趋于直线AB;当ρ=0.4142且CM为AB的垂直平分线时为圆弧。
对于A线,所取的参数范围包括3个水平,如表1所示,
表1
参数 |
水平1 |
水平2 |
水平3 |
长度L(米) |
5.187 |
6 |
7 |
角度θ(度) |
57.6 |
50 |
65 |
形状参数ρ |
0.46 |
0.36 |
0.56 |
利用正交设计表对表1的数据进行组合,可以得到表2的数据组合:
表2
序号 |
长度L(米) |
角度θ(度) |
形状参数ρ |
1 |
5.187 |
57.6 |
0.46 |
2 |
5.187 |
50 |
0.36 |
3 |
5.187 |
65 |
0.56 |
4 |
6 |
57.6 |
0.36 |
5 |
6 |
50 |
0.56 |
6 |
6 |
65 |
0.46 |
7 |
7 |
57.6 |
0.56 |
8 |
7 |
50 |
0.46 |
9 |
7 |
65 |
0.36 |
对于D线,由于其角度参数为90°,可以为定值,所以变化范围只限于列车头部长度L和形状参数ρ,如表3所示,
表3
参数 |
水平1 |
水平2 |
水平3 |
头部长度L(米) |
5.187 |
6 |
7 |
形状参数ρ |
0.5 |
0.6 |
0.7 |
利用正交设计表对表3的数据进行组合,可以得到表4的数据组合:
表4
试验方案 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
形状参数ρ |
0.5 |
0.5 |
0.5 |
0.6 |
0.6 |
0.6 |
0.7 |
0.7 |
0.7 |
头部长度L(米) |
5.187 |
6.0 |
7.0 |
5.187 |
6.0 |
7.0 |
5.187 |
6.0 |
7.0 |
(3)根据所述设计表中的参数,利用计算流体动力学(CFD)软件分别对纵向对称面最大控制型线和俯视最大控制型线的试验方案进行二维流动计算,对纵向对称面最大控制型线的试验方案计算阻力,选取使得阻力系数最小的参数值作为纵向对称面最大控制型线的参数,对俯视最大控制型线的试验方案计算会车压力波的幅值,选取使得会车压力波的幅值最小的参数值作为俯视最大控制型线的参数。
所述步骤(2)中,列出的数据设计表为正交设计表或者均匀设计表,在步骤(3)中,如果选取的参数组合未在正交设计表或者均匀设计表中出现,则对所述选取的参数组合计算阻力,验证得到的阻力与其它参数组合计算得到的阻力相比是否是最小的。采用这种方法,可以大大减少运算量,提高设计效率。
对于A线,对表2的参数进行计算,得到如表5所示的数据:
表5
其中,表5中的序号与表2中的序号相对应,对表5中序号为1~9的数据进行分析,用K
1表示第一种水平试验的阻力系数之和,
为K
1的算术平均值,即
用K
2表示第一种水平试验的阻力系数之和,
为K
2的算术平均值,即
用K
3表示第一种水平试验的阻力系数之和,
为K
3的算术平均值,即
极差为
和
中最大值和最小值的差值,反映了因子变化时,目标函数变化的幅度,如表6所示,
表6
对表6进行分析可以得到以下结论:
从上表中可以看出阻力系数随三个参数的变化规律:
(1)随着长度的增加,阻力系数减少;
(2)随着角度θ的增加,阻力系数增加;
(3)随着形状参数ρ的增加,阻力系数增加。
根据极差的大小可以判断,这三个参数影响阻力系数的主次顺序为:角度θ、长度L和形状参数ρ。因为极差越大,说明与其对应的因子对阻力系数的影响越显著,阻力系数对因子的水平变化越敏感。
步骤(3)中对于纵向对称面最大控制型线除计算阻力外,还计算升力,确定根据选取的参数所计算的升力是否在安全范围之内;如果比其它大多数参数组合所得到的升力都小,则回到第(2)步重新确定参数范围。如果升力过大,会导致高速列车在运行时失去控制。但是,在保证安全的范围内,较大的升力可以降低对高速列车支承力的要求,所以还是有好处的。
因此,最优的试验方案是:L3-θ2-ρ2。把这种方案命名为第10种,计算结果如表5中序号10,确实这种方案的阻力系数最小,同时这种方案的升力系数也较大,仅次于第7种和第8种。表5中,TR08为现有高速列车的参数计算值,与TR08的纵向型面相比,第10种方案的阻力系数下降了8%,升力增加了3.3%。
对于D线,首先需要在列车上的头车、车厢和尾车上取3个测点,如图3所示,对于这3个测点采用表4的数据分别计算会车压力波,其中Pmax为会车压力波的最大值,Pmin为会车压力波的最小值,ΔP为会车压力波的幅值,M
1、M
2和M
3分别表示3个水平所对应的会车压力波的幅值之和,
和
分别表示3个水平所对应的会车压力波的算术平均值。
测点1的会车压力波的数据如表7所示,
表7
序号 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
TR08 |
Pmax |
9244 |
9160 |
8973 |
9370 |
9317 |
9111 |
9548 |
9251 |
9262 |
9288 |
Pmin |
-8380 |
-7340 |
-7370 |
-7906 |
-7609 |
-7268 |
-7243 |
-6893 |
-6378 |
8185 |
ΔP |
17624 |
16500 |
16343 |
17276 |
16926 |
16379 |
16791 |
16144 |
15640 |
17473 |
测点1的会车压力波幅值的直观分析数据如表8所示,
表8
测点2的会车压力波的数据如表9所示,
表9
序号 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
TR08 |
Pmax |
7756 |
7773 |
7806 |
7761 |
7689 |
7610 |
7853 |
7681 |
7670 |
7730 |
Pmin |
-12929 |
-11919 |
-11597 |
-13520 |
-13571 |
-13543 |
-11517 |
-9371 |
-9471 |
-12192 |
ΔP |
20685 |
19692 |
19403 |
21281 |
21260 |
21153 |
19370 |
17052 |
17141 |
19922 |
测点2的会车压力波幅值的直观分析数据如表10所示,
表10
测点3的会车压力波的数据如表11所示,
表11
序号 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
TR08 |
Pmax |
5153 |
5450 |
5598 |
5157 |
5464 |
5541 |
4246 |
5502 |
5717 |
5035 |
Pmin |
-20415 |
-20251 |
-20907 |
-21940 |
-24192 |
-23580 |
-16189 |
-17294 |
-17763 |
-19002 |
ΔP |
25568 |
25701 |
26505 |
27097 |
29656 |
29121 |
20435 |
22796 |
23480 |
24037 |
测点3的会车压力波幅值的直观分析数据如表12所示,
表12
对表7~表12的数据进行分析可以得到,形状参数对于会车压力波幅值的影响更大,而且形状参数越大,会车压力波的幅值越小;而列车头部长度越长,会车压力波也越小,因此L3-ρ3的数据最为合理。
(4)根据所选取的纵向对称面最大控制型线和俯视最大控制型线的参数,结合车体截面外廓型线通过无限插值得到Coons曲面,根据需要增加辅助控制型线数量对Coons曲面进行必要的修改,确定三维高速列车车头模型,如图4所示;
(5)对三维高速列车车头模型开展三维流动计算,计算得到的阻力与现有的高速列车相比,验证阻力是否减小,如果是,则设计完成;如果不是,再回到第(2)步,更改6个参数的变化范围,重复第(3)步到第(5)步。
根据本发明优化的车头参数进行计算,得到的结果于现有的车头数据进行比较,如表13所示,
表13
列车 |
压阻(N) |
摩阻(N) |
总阻力(N) |
阻力系数 |
TR08 |
21270.78 |
29542.75 |
50813.53 |
0.426 |
优化列车 |
7850.14 |
38234.42 |
46084.54 |
0.386 |
变化比例 |
63.1% |
-29.4% |
9.4% |
9.4% |
表13中的变化比例指:(TR08的值-优化列车的值)/TR08的值。
由表13可以看出,优化的车头的阻力比现有的车头大大降低。采用本发明得到的车头与现有的车头外型的比较可参见图5所示。
与现有的技术相比,本发明具有如下优点:
1.车头表面建模简单方便,本发明选取了三条主控制型线和两条辅助控制型线,并且将三条主控制型线用二次曲线来表示,避免了采用手工设计型线的繁琐;由于二次曲线的控制参数比较少,避免了采用NURBS曲线控制参数多的缺点;采用两条辅助控制型线和Coons曲面生成,可以适用于一些复杂外形,避免了用某个统一的数学函数来描述车头适用范围较窄的不足。
2.由于控制型线的参数少,可以方便地生成所需要的外形,使各种优化方法(如正交设计表等)可以有效的开展,满足工程实用的要求。
3.控制型线的设计与CFD软件的计算密切相关,互相耦合,使生成的外形可以满足减少阻力的目标,以提高速度,减少能耗。同时本发明纵向对称面最大控制型线主要控制列车的阻力和升力,俯视最大控制型线主要控制列车的会车压力波,分别对这两种型线采用二维流动计算,可以在较短的时间内得到优化的方案。避免了完全采用三维流动计算进行优化的费时费力。