CN116663455B - 一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法 - Google Patents
一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法 Download PDFInfo
- Publication number
- CN116663455B CN116663455B CN202310683954.6A CN202310683954A CN116663455B CN 116663455 B CN116663455 B CN 116663455B CN 202310683954 A CN202310683954 A CN 202310683954A CN 116663455 B CN116663455 B CN 116663455B
- Authority
- CN
- China
- Prior art keywords
- wing
- lift
- lift wing
- speed
- aerodynamic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 32
- 238000010187 selection method Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 60
- 238000013461 design Methods 0.000 claims abstract description 48
- 238000011160 research Methods 0.000 claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 24
- 238000009434 installation Methods 0.000 claims description 22
- 101100457838 Caenorhabditis elegans mod-1 gene Proteins 0.000 claims description 13
- 101150110972 ME1 gene Proteins 0.000 claims description 13
- 238000004088 simulation Methods 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000005477 standard model Effects 0.000 claims description 6
- 238000010835 comparative analysis Methods 0.000 claims description 5
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000010205 computational analysis Methods 0.000 claims 1
- 230000002093 peripheral effect Effects 0.000 claims 1
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000011161 development Methods 0.000 description 5
- 230000001965 increasing effect Effects 0.000 description 5
- 238000005265 energy consumption Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 239000000725 suspension Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法，以布设在高速列车车顶的升力翼为研究对象，以高速铁路限界为约束条件，建立多方案升力翼三维几何模型，采用计算流体动力学计算方法，从流线型外观设计、气动力特性、流场结构等多目标参数综合评价和选择，确定适应下一代时速400+km高速列车升力翼翼型的基本三维结构方案。该基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法能够给现阶段高速列车装配升力翼的科学选型提供一个典型参考方案，填补了行业该方面的技术空白，能够有效解决针对现有高速列车升力翼选型适用性低、实车上线试验困难及试验成本高等问题。
Description
技术领域
本发明涉及轨道交通装备制造及列车空气动力学领域，具体涉及一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法。
背景技术
在分析世界高速列车及我国高速列车发展现状的基础上，下一代高速动车组应该在更高速、更安全、更环保、更经济、更舒适、更友好6个方面具有明显的代际特征，其中，更高速体现在深入优化改善轮轨关系及弓网关系，进一步降低列车运行阻力及全方面的轻量化设计，实现时速达400+km的商业运营，形成谱系化产品；更安全体现在持续提升全生命周期的智能化列车运行安全保障及全方位的预警监测；更环保体现在进一步提高环保性指标，其中运行阻力减少5-10％，能耗减少10-15％，车厢内环境噪声降低3-5dB，同时加大高效能、轻量化等领域新材料及结构设计应用，实现轴重减少5-10％；更经济体现在深化智能化技术，提升故障预测与健康管理，实现预期全生命周期成本比现有高速动车组降低10-15％，可用性提升5-10％。
随着预期设计运行速度的进一步提升，高速列车运行牵引能耗将持续增大、轮轨磨耗将进一步加剧、车体振动及轮轨动力学问题将更加突出、列车空气动力学效应愈加明显，此过程不可避免的面临牵引供电资源紧张、列车关键走行部件维修周期压缩与使用寿命降低、气动噪声与气动阻力骤增，以及复杂风环境下等工况下的列车行车安全问题。为了在更高速度运行条件下降低列车全寿命周期成本，研究提出了带有升力翼的高速列车概念，突破传统高速列车气动外形设计理念，结合高速列车和飞行器各自优势，以期通过增加列车气动升力，实现高速列车整体节能降耗。
20世纪末，日本东北大学最早提出了“气动悬浮列车”设计概念，通过在地面附近布置地效翼，利用地面效应增升，给列车提供升力，同时对气动悬浮列车所用翼型进行了初步设计研究，认为其运载经济效率要高于磁浮列车和高速民航客机，并制作出了气动悬浮列车实验车型，提出了一种添加升力翼的概念设计方案，在车顶和车底侧面布置“仿机翼”，并给出了一些可供选择的翼型。
对于翼型的选择及运用方面，现阶段也做过大量研究，从对风洞环境下的NACA0015翼型数值仿真研究及试验结果对比中，发现两者的差异主要是由翼型的气流分离位置差异所导致的；在针对NACA0015翼型不同迎角下的气动特性数值仿真研究中，发现该翼型在10°以上迎角下即出现失速现象；通过风洞试验对二维NACA4412翼型在不同迎角、不同飞行高度（升力翼与车顶之间的距离）时的升力和阻力变化规律进行了研究，发现翼型在各迎角下靠近地面都会导致上表面的吸力损失；在针对低雷诺数条件下多种翼型的气动特性的试验研究中，得到了气动特性较优的翼型及不同迎角下的翼型升力变化规律。
上述相关研究表明，翼型的气动设计直接关系着高速列车气动特性的优劣，设计出具备良好气动特性的升力翼是升力翼列车技术的关键。截止目前，围绕这一目标，研究提出了多种设计方案，比如授权公告号为CN113602299B，发明名称为高速列车气动力调控的伸缩翼装置、高速列车及控制方法，授权公告号为CN210133111U，发明名称为高速轨道交通列车侧翼升力控制机构的中国实用新型专利，授权公告号为CN202175052U和CN202175053U公开的一种高速列车车翼装置等。但总体来看并未充分与高速列车的发展实际相结合，尤其表现在升力翼选型方法、结构设计方案、安装布置形式、控制方式及实车应用等方面基本处于空白状态。
基于此，在现阶段我国围绕高速智能绿色铁路装备大发展的背景下，研发选择一种适应现阶段铁路限界约束的的增升效果明显、阻力系数小、气动噪声小、安装空间小的高速列车升力翼是目前高速列车提速运行及贯彻节能降耗发展绿色铁路装备亟待解决的问题之一。
发明内容
针对中国现阶段高速铁路及标准动车组标准限界，为能够科学地研发选择一种适应现阶段铁路限界约束的的增升效果明显、阻力系数小、气动噪声小、安装空间小的高速列车升力翼，本发明提出一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法。
为实现上述技术目的，本发明采用如下技术方案予以实现：
一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法，所述方法以布设在高速列车车顶的升力翼为研究对象，以高速铁路限界为约束条件，建立多方案升力翼三维几何模型，采用计算流体动力学计算方法，从流线型外观设计、气动力特性、流场结构等多目标参数综合评价和选择，确定适应下一代时速400+km高速列车升力翼翼型的基本三维结构方案。具体基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法包括以下步骤：
1）明确限界内可供升力翼布设的有效安装空间：
11）根据铁路建筑限界、桥隧限界条件及具体适用车型车辆限界条件，确定高速列车升力翼装置在限界内车顶上部最大有效安装空间（S1），以及车体两侧有效安装空间（S2）；
12）考虑下一代时速400+km高速列车车侧流线型表面设计、车门及车窗位置、风挡及车顶设备等对升力翼两侧布设位置的无序性影响，布设选择限界内车顶上部最大有效安装空间（S1），舍弃车体两侧有效安装空间（S2）；
2）建立多方案升力翼数字化三维几何模型：
21）确定高速列车升力翼基本翼剖面：参考飞机等飞行器的机翼设计理念与方法，提出适合对于时速400~450km范围运行的轮轨列车的翼剖面备选方案：平凸型、凹凸型和双凸型3种常见翼型，在该速度范围内确定相同气动设计尺度，采用计算流体动力学方法，模拟计算并优先考虑小阻力大升力的翼剖面设计，同时结合可升降升力翼与高速列车车顶面接触空气动力学性能，选用适配车顶局部表面的平凸型翼剖面进行高速列车升力翼的设计研究；
22）确定平凸型升力翼数字化三维几何模型：以平凸型翼剖面为雏形，根据高速列车升力翼工作气动环境、尺度空间和多角度姿态特征等，初步选定升力翼弦长、最大厚度、前缘半径、后缘角为固定参数，主要针对升力翼横向流线型曲面下偏角和纵向前后掠角对升力翼气动性能的主要影响，给出零偏角平凸型（Mod1）、带偏角平凸型（Mod2）和带掠角平凸型（Mod3）3种升力翼设计方案；其中，零偏角平凸型（Mod1）升力翼为横向任意位置翼剖面几何参数均一致设计；带偏角平凸型（Mod2）升力翼为横向从中间向两侧翼剖面逐渐过渡减薄；带掠角平凸型（Mod3）为在所述带偏角平凸型（Mod2）的基础上两侧向后同时在水平面内前缘设计存在前掠角，后缘设计存在后掠角；
3）建立升力翼流体动力学计算模型：
31）设定数学模型：设定计算采用有粘可压缩的Navier-Stokes方程，湍流模拟采用标准的涡黏性模型，其中所述标准的涡黏性模型将湍流黏性系数和湍流动能与湍流动能耗散率相关联，求解模型采用有限体积法，对流项采用二阶迎风和有界中心格式，离散项采用二阶隐式格式；
32）确定带升力翼高速列车计算模型：以中国标准动车组外观构造为雏形，选择包含1节头车、1节中间车及1节尾车组成的1比1的三编组标准模型，其中升力翼按单组设置，布设在头车车顶纵向中心处；
33）确定计算域：以带单排升力翼的所述三编组标准模型的列车总长（L）为参考对象，创建长为4倍列车总长（L），宽为2倍列车总长（L），高为1倍列车总长（L）的长方体外流场计算域；所述带升力翼高速列车计算模型位于所述长方体外流场计算域对称边界中部，其中高速列车头车鼻尖距离外流场前侧面距离为1倍列车总长（L），尾车鼻尖距离外流场后侧面距离为2倍列车总长（L），所述带升力翼高速列车计算模型车体底面距离外流场下底面为0.4m；
34）计算设置：采用有限体积法进行计算域内非结构化网格划分及求解，其中对于高速列车头车及尾车、流固接触面和升力翼外围等流场变化突出的区域，采用全局和局部网格相配合叠加的处理模式，计算湍流强度设为0.5%，残差取值范围为10-6~10-4；
4）多方案带升力翼高速列车气动力特性比较分析：
41）气动力比较分析：分别以零偏角平凸型（Mod1）、带偏角平凸型（Mod2）和带掠角平凸型（Mod3）3种不同升力翼设计方案为模拟计算对象，设定计算风速为250、300、350、400、450km/h，升力翼工作迎角为16°，升力翼起升高度为400mm，进行气动特性分析，分析计算不同速度等级条件下升力翼所受气动升力F L 与气动阻力值F D ，根据Davis公式采用过原点开口向上的抛物线拟合升力翼所受气动升力F L 与气动阻力F D 随速度变化函数；根据拟合曲线判断确定小气动阻力、大气动升力变化趋势的升力翼翼型方案I；
42）气动系数比较分析：在此基础上采用统一量纲升力系数C L 及阻力系数C D 分别计算升力翼所受气动升和阻力大小，其中升力系数C L 计算公式为C L =F L /(0.5ρA L V 2 )，阻力系数C D 计算公式为C D =F D /(0.5ρA D V 2 )，式中：A L 为升力翼垂向投影面积，所述垂向投影面积计算公式为A L =ckcosγ；其中，c为升力翼弦长，k为升力翼在限界内的横向展长，k在限界内的约束条件为K＞k≥1600mm，K为铁路建筑限界内可安装升力翼空间的最大横向尺寸；A D 为升力翼纵向投影面积，其中所述升力翼纵向投影面积A D 采用三维几何最大轮廓投影计算的方式在数字化三维辅助设计平台上计算；ρ为15℃标准大气压下的空气密度；V为计算风速；经流体动力学计算分析，采用列表方式比对各型升力翼气动力系数，样式为：
| 翼型 | 零偏角平凸型（Mod1） | 带偏角平凸型（Mod2） | 带掠角平凸型（Mod3） |
| 升力系数C L | 计算值n1 | 计算值n2 | 计算值n3 |
| 阻力系数C D | 计算值n4 | 计算值n5 | 计算值n6 |
根据气动力系数计算表，判断确定小阻力系数、大升力系数的升力翼翼型方案II；
43）升力翼表面气动压力系数比较分析：根据步骤3）所建立升力翼流体动力学计算模型，在同一速度等级条件下，模拟计算的升力翼中线翼剖面轮廓线上每个计算单元节点的静压值P，计算每个所述计算单元节点的压力系数C P ，其中所述压力系数计算公式为C P =(P-P ∞)/(0.5ρV 2 )，式中：P为计算单元节点的静压值；P ∞为模拟远环境静压值，也即为15℃环境温度对应的标准大气压值；ρ为15℃标准大气压下的空气密度；V为计算风速；
以升力翼前缘到后缘纵向长度为参考值，将其标度归一化处理，在同一坐标系中依次导入多方案升力翼纵向归一化标度-气动压力系数散点，分别以前缘和后缘点为分界拟合出升力翼上弧线和下弦线对应的气动压力曲线；分别对比分析升力翼翼剖面的上弧线和下弦线受压波动情况，重点以上弧线压力波动变化幅度、频率和峰值差为主要对比指标，分析确定波动变化幅度小、频率小和峰值差小的升力翼翼型方案III；
44）升力翼纵向对称面上流场结构比较分析：根据步骤3）所建立升力翼流体动力学计算模型，在同一速度等级条件下，模拟计算各型升力翼表面和纵向对称面上稳态压力分布及流动迹线，以升力翼上弧面尾部负压强度小、下弧面底部正压强度小、流场干扰效应小、涡量小及影响范围小为对比指标，确定升力翼翼型方案IV；
5）确定适配升力翼翼型结构：
分别依据步骤41）速度-气动拟合曲线比较分析优选的升力翼翼型方案I，步骤42）气动系数比较分析优选的升力翼翼型方案II，步骤43）升力翼表面气动压力系数比较分析优选的升力翼翼型方案III，步骤44）升力翼纵向对称面上流场结构比较分析优选的升力翼翼型方案IV，根据高速列车400~600km/h速度等级范围内升力翼的主要应用指标次序综合分析确定；其中升力翼的主要应用指标重要程度为：气动力＞气动系数＞升力翼表面气动压力系数＞升力翼纵向对称面上流场结构。
作为优选地，所述铁路建筑限界、桥隧限界条件为中国高速铁路标准限界范围，所述适用车型车辆限界为依据中国标准动车组横向截面轮廓的空间限界。
作为优选地，高速列车升力翼装置在限界横断内车顶上部所述最大有效安装空间（S1）为长方形空间，所述长方形空间关于纵向中心对称，其最大横向尺寸为3400mm。
本发明的有益效果为：该基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法能够给现阶段高速列车装配升力翼科学选型提供一个典型参考方案，有效填补该方面的技术空白，能够有效解决针对现有高速列车升力翼选型适用性低、实车试验困难及成本高等问题。
附图说明
图1为本发明的一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法流程图；
图2为本发明的升力翼在适应高速列车车顶上有效安装空间的安装布设示意图；
图3为本发明的升力翼在铁路限界内车顶上部有效布置空间选择示意图；
图4为本发明的升力翼在铁路限界内车体两侧有效布置空间选择示意图；
图5为本发明的多方案升力翼模型及参数设计示意图；
图6为本发明的多方案升力翼模型表面气动压力系数比较分析图。
具体实施方式
下面结合附图对本发明作进一步说明：
如图1所示，一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法，所述方法以布设在高速列车车顶的升力翼为研究对象，以高速铁路限界为约束条件，建立多方案升力翼三维几何模型，采用计算流体动力学计算方法，从流线型外观设计、气动力特性、流场结构等多目标参数综合评价和选择，确定适应下一代时速400+km高速列车升力翼翼型的基本三维结构方案，具体基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法包括以下步骤：
1）明确限界内可供升力翼布设的有效安装空间：（参考图2、3、4）
11）根据铁路建筑限界、桥隧限界条件及具体适用车型车辆限界条件，确定高速列车升力翼装置在限界内车顶上部最大有效安装空间（S1），以及车体两侧有效安装空间（S2）；
12）考虑下一代时速400+km高速列车车侧流线型表面设计、车门及车窗位置、风挡及车顶设备等对升力翼两侧布设位置的无序性影响，所述方法布设选择限界内车顶上部最大有效安装空间（S1），舍弃车体两侧有效安装空间（S2）；
2）建立多方案升力翼数字化三维几何模型：（参考图5）
21）确定高速列车升力翼基本翼剖面：参考飞机等飞行器的机翼设计理念与方法，提出适合对于时速400~450km范围运行的轮轨列车的翼剖面备选方案：平凸型、凹凸型和双凸型3种常见翼型，在该速度范围内确定相同气动设计尺度，采用计算流体动力学方法，模拟计算并优先考虑小阻力大升力的翼剖面设计，同时结合可升降升力翼与高速列车车顶面接触空气动力学性能，选用适配车顶局部表面的平凸型翼剖面进行高速列车升力翼的设计研究；
22）确定平凸型升力翼数字化三维几何模型：以平凸型翼剖面为雏形，根据高速列车升力翼工作气动环境、尺度空间和多角度姿态特征等，初步选定升力翼弦长、最大厚度、前缘半径、后缘角为固定参数，主要针对升力翼横向流线型曲面下偏角和纵向前后掠角对升力翼气动性能的主要影响，给出零偏角平凸型（Mod1）、带偏角平凸型（Mod2）和带掠角平凸型（Mod3）3种升力翼设计方案；其中，零偏角平凸型（Mod1）升力翼为横向任意位置翼剖面几何参数均一致设计；带偏角平凸型（Mod2）升力翼为横向从中间向两侧翼剖面逐渐过渡减薄；带掠角平凸型（Mod3）为在所述带偏角平凸型（Mod2）的基础上两侧向后同时在水平面内前缘设计存在前掠角，后缘设计存在后掠角；
3）建立升力翼流体动力学计算模型：
31）设定数学模型：设定计算采用有粘可压缩的Navier-Stokes方程，湍流模拟采用标准的涡黏性模型，其中所述标准的涡黏性模型将湍流黏性系数和湍流动能与湍流动能耗散率相关联，求解模型采用有限体积法，对流项采用二阶迎风和有界中心格式，离散项采用二阶隐式格式；
32）确定带升力翼高速列车计算模型：以中国标准动车组外观构造为雏形，选择包含1节头车、1节中间车及1节尾车组成的1比1的三编组标准模型，其中升力翼按单组设置，布设在头车车顶纵向中心处；
33）确定计算域：以带单排升力翼的所述三编组标准模型的列车总长（L）为参考对象，创建长为4倍列车总长（L），宽为2倍列车总长（L），高为1倍列车总长（L）的长方体外流场计算域；所述带升力翼高速列车计算模型位于所述长方体外流场计算域对称边界中部，其中高速列车头车鼻尖距离外流场前侧面距离为1倍列车总长（L），尾车鼻尖距离外流场后侧面距离为2倍列车总长（L），所述带升力翼高速列车计算模型车体底面距离外流场下底面为0.4m；
34）计算设置：采用有限体积法进行计算域内非结构化网格划分及求解，其中对于高速列车头车及尾车、流固接触面和升力翼外围等流场变化突出的区域，采用全局和局部网格相配合叠加的处理模式，计算湍流强度设为0.5%，残差取值范围为10-6~10-4；
4）多方案带升力翼高速列车气动力特性比较分析：
41）气动力比较分析：分别以零偏角平凸型（Mod1）、带偏角平凸型（Mod2）和带掠角平凸型（Mod3）3种不同升力翼设计方案为模拟计算对象，设定计算风速为250、300、350、400、450km/h，升力翼工作迎角为16°，升力翼起升高度为400mm，进行气动特性分析，分析计算不同速度等级条件下升力翼所受气动升力F L 与气动阻力值F D ，根据Davis公式采用过原点开口向上的抛物线拟合升力翼所受气动升力F L 与气动阻力F D 随速度变化函数；根据拟合曲线判断确定小气动阻力、大气动升力变化趋势的升力翼翼型方案I；
42）气动系数比较分析：在此基础上采用统一量纲升力系数C L 及阻力系数C D 分别计算升力翼所受气动升和阻力大小，其中升力系数C L 计算公式为C L =F L /(0.5ρA L V 2 )，阻力系数C D 计算公式为C D =F D /(0.5ρA D V 2 )，式中：A L 为升力翼垂向投影面积，所述垂向投影面积计算公式为A L =ckcosγ；其中，c为升力翼弦长，k为升力翼在限界内的横向展长，k在限界内的约束条件为K＞k≥1600mm，K为铁路建筑限界内可安装升力翼空间的最大横向尺寸；A D 为升力翼纵向投影面积，其中所述升力翼纵向投影面积A D 采用三维几何最大轮廓投影计算的方式在数字化三维辅助设计平台上计算；ρ为15℃标准大气压下的空气密度；V为计算风速；经流体动力学计算分析，采用列表方式比对各型升力翼气动力系数，样式为：
| 翼型 | 零偏角平凸型（Mod1） | 带偏角平凸型（Mod2） | 带掠角平凸型（Mod3） |
| 升力系数C L | 计算值n1 | 计算值n2 | 计算值n3 |
| 阻力系数C D | 计算值n4 | 计算值n5 | 计算值n6 |
根据气动力系数计算表，判断确定小阻力系数、大升力系数的升力翼翼型方案II；
43）升力翼表面气动压力系数比较分析（参考图6）：根据步骤3）所建立升力翼流体动力学计算模型，在同一速度等级条件下，模拟计算的升力翼中线翼剖面轮廓线上每个计算单元节点的静压值P，计算每个所述计算单元节点的压力系数C P ，其中所述压力系数计算公式为C P =(P-P ∞)/(0.5ρV 2 )，式中：P为计算单元节点的静压值；P ∞为模拟远环境静压值，也即为15℃环境温度对应的标准大气压值；ρ为15℃标准大气压下的空气密度；V为计算风速；
以升力翼前缘到后缘纵向长度为参考值，将其标度归一化处理，在同一坐标系中依次导入多方案升力翼纵向归一化标度-气动压力系数散点，分别以前缘和后缘点为分界拟合出升力翼上弧线和下弦线对应的气动压力曲线；分别对比分析升力翼翼剖面的上弧线和下弦线受压波动情况，重点以上弧线压力波动变化幅度、频率和峰值差为主要对比指标，分析确定波动变化幅度小、频率小和峰值差小的升力翼翼型方案III；
44）升力翼纵向对称面上流场结构比较分析：根据步骤3）所建立升力翼流体动力学计算模型，在同一速度等级条件下，模拟计算各型升力翼表面和纵向对称面上稳态压力分布及流动迹线，以升力翼上弧面尾部负压强度小、下弧面底部正压强度小、流场干扰效应小、涡量小及影响范围小为对比指标，确定升力翼翼型方案IV；
5）确定适配升力翼翼型结构：
分别依据步骤41）速度-气动拟合曲线比较分析优选的升力翼翼型方案I，步骤42）气动系数比较分析优选的升力翼翼型方案II，步骤43）升力翼表面气动压力系数比较分析优选的升力翼翼型方案III，步骤44）升力翼纵向对称面上流场结构比较分析优选的升力翼翼型方案IV，根据高速列车400~600km/h速度等级范围内升力翼的主要应用指标次序综合分析确定；其中升力翼的主要应用指标重要程度为：气动力＞气动系数＞升力翼表面气动压力系数＞升力翼纵向对称面上流场结构。
如图2、3、4所示，所述铁路建筑限界、桥隧限界条件为中国高速铁路标准限界范围，所述适用车型车辆限界为依据中国标准动车组横向截面轮廓的空间限界，高速列车升力翼装置在限界横断内车顶上部所述最大有效安装空间（S1）为长方形空间，所述长方形空间关于纵向中心对称，其最大横向尺寸为3400mm。
针对利用本方法选定升力翼方案的计算特性说明：
与零偏角平凸型（Mod1）升力翼和带掠角平凸型（Mod3）升力翼相比，带偏角平凸型（Mod2）升力翼气动设计在整体增升效果、减阻特性、流场效应及与车顶流线型契合度等方面具备较好的综合优势，运行升力系数为1.15，其中在时速450km工作时，贡献气动升力达19.08kN，而气动阻力仅占升力的32.7％。
其中需要说明的是，本文献中提到的“左”、“右”、“前”、“后”、“内”、“外”、“上”、“下”等指示方位或位置关系为基于附图所示的位置关系，仅是为了便于描述本技术方案和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造或操作。因此不能理解为对技术方案的限制，所述连接关系可以指直接连接关系，也可以指间接连接关系。
显然，本领域的技术人员可以对本发明进行各种改动和变形而不脱离本发明的精神和范围，这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (3)
1.一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法，其特征在于：所述方法以布设在高速列车车顶的升力翼为研究对象，以高速铁路限界为约束条件，建立多方案升力翼三维几何模型，采用计算流体动力学计算方法，从流线型外观设计、气动力特性和流场结构多目标参数综合评价和选择，确定适应下一代时速400+km高速列车升力翼翼型的基本三维结构方案，具体基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法包括以下步骤：
1）明确限界内可供升力翼布设的有效安装空间：
11）根据铁路建筑限界、桥隧限界条件及具体适用车型车辆限界条件，确定高速列车升力翼装置在限界内车顶上部最大有效安装空间（S1），以及车体两侧有效安装空间（S2）；
12）考虑下一代时速400+km高速列车车侧流线型表面设计、车门及车窗位置、风挡和车顶设备对升力翼两侧布设位置的无序性影响，所述方法的升力翼布设选择限界内车顶上部最大有效安装空间（S1），舍弃车体两侧有效安装空间（S2）；
2）建立多方案升力翼数字化三维几何模型：
21）确定高速列车升力翼基本翼剖面：参考飞机和飞行器的机翼设计理念与方法，提出适合对于时速400~450km范围运行的轮轨列车的翼剖面备选方案：平凸型、凹凸型和双凸型3种常见翼型，在该速度范围内确定相同气动设计尺度，采用计算流体动力学方法，模拟计算并优先考虑小阻力大升力的翼剖面设计，同时结合可升降升力翼与高速列车车顶面接触空气动力学性能，选用适配车顶局部表面的平凸型翼剖面进行高速列车升力翼的设计研究；
22）确定平凸型升力翼数字化三维几何模型：以平凸型翼剖面为雏形，根据高速列车升力翼工作气动环境、尺度空间和多角度姿态特征，初步选定升力翼弦长、最大厚度、前缘半径、后缘角为固定参数，主要针对升力翼横向流线型曲面下偏角和纵向前后掠角对升力翼气动性能的主要影响，给出零偏角平凸型（Mod1）、带偏角平凸型（Mod2）和带掠角平凸型（Mod3）3种升力翼设计方案；其中，零偏角平凸型（Mod1）升力翼为横向任意位置翼剖面几何参数均一致设计；带偏角平凸型（Mod2）升力翼为横向从中间向两侧翼剖面逐渐过渡减薄；带掠角平凸型（Mod3）为在所述带偏角平凸型（Mod2）的基础上两侧向后同时在水平面内前缘设计存在前掠角，后缘设计存在后掠角；
3）建立升力翼流体动力学计算模型：
31）设定数学模型：设定计算采用有粘可压缩的Navier-Stokes方程，湍流模拟采用标准的涡黏性模型，其中所述标准的涡黏性模型将湍流黏性系数和湍流动能与湍流动能耗散率相关联，求解模型采用有限体积法，对流项采用二阶迎风和有界中心格式，离散项采用二阶隐式格式；
32）确定带升力翼高速列车计算模型：以中国标准动车组外观构造为雏形，选择包含1节头车、1节中间车及1节尾车组成的1比1的三编组标准模型，其中升力翼按单组设置，布设在头车车顶纵向中心处；
33）确定计算域：以带单排升力翼的所述三编组标准模型的列车总长（L）为参考对象，创建长为4倍列车总长（L），宽为2倍列车总长（L），高为1倍列车总长（L）的长方体外流场计算域；所述带升力翼高速列车计算模型位于所述长方体外流场计算域对称边界中部，其中高速列车头车鼻尖距离外流场前侧面距离为1倍列车总长（L），尾车鼻尖距离外流场后侧面距离为2倍列车总长（L），所述带升力翼高速列车计算模型车体底面距离外流场下底面为0.4m；
34）计算设置：采用有限体积法进行计算域内非结构化网格划分及求解，其中对于高速列车头车及尾车、流固接触面和升力翼外围流场变化突出的区域，采用全局和局部网格相配合叠加的处理模式，计算湍流强度设为0.5%，残差取值范围为10-6~10-4；
4）多方案带升力翼高速列车气动力特性比较分析：
41）气动力比较分析：分别以零偏角平凸型（Mod1）、带偏角平凸型（Mod2）和带掠角平凸型（Mod3）3种不同升力翼设计方案为模拟计算对象，设定计算风速为250、300、350、400、450km/h，升力翼工作迎角为16°，升力翼起升高度为400mm，进行气动特性分析，分析计算不同速度等级条件下升力翼所受气动升力F L 与气动阻力值F D ，根据Davis公式采用过原点开口向上的抛物线拟合升力翼所受气动升力F L 与气动阻力F D 随速度变化函数；根据拟合曲线判断确定小气动阻力、大气动升力变化趋势的升力翼翼型方案I；
42）气动系数比较分析：在此基础上采用统一量纲升力系数C L 及阻力系数C D 分别计算升力翼所受气动升和阻力大小，其中升力系数C L 计算公式为C L =F L /(0.5ρA L V 2 )，阻力系数C D 计算公式为C D =F D /(0.5ρA D V 2 )，式中：A L 为升力翼垂向投影面积，所述垂向投影面积计算公式为A L =ckcosγ；其中，c为升力翼弦长，k为升力翼在限界内的横向展长，k在限界内的约束条件为K＞k≥1600mm，K为铁路建筑限界内可安装升力翼空间的最大横向尺寸；A D 为升力翼纵向投影面积，其中所述升力翼纵向投影面积A D 采用三维几何最大轮廓投影计算的方式在数字化三维辅助设计平台上计算；ρ为15℃标准大气压下的空气密度；V为计算风速；经流体动力学计算分析，采用列表方式比对各型升力翼气动力系数，样式为：
根据气动力系数计算表，判断确定小阻力系数、大升力系数的升力翼翼型方案II；
43）升力翼表面气动压力系数比较分析：根据步骤3）所建立升力翼流体动力学计算模型，在同一速度等级条件下，模拟计算的升力翼中线翼剖面轮廓线上每个计算单元节点的静压值P，计算每个所述计算单元节点的压力系数C P ，其中所述压力系数计算公式为C P =(P-P ∞)/(0.5ρV 2 )，式中：P为计算单元节点的静压值；P ∞为模拟远环境静压值，也即为15℃环境温度对应的标准大气压值；ρ为15℃标准大气压下的空气密度；V为计算风速；
以升力翼前缘到后缘纵向长度为参考值，将其标度归一化处理，在同一坐标系中依次导入多方案升力翼纵向归一化标度-气动压力系数散点，分别以前缘和后缘点为分界拟合出升力翼上弧线和下弦线对应的气动压力曲线；分别对比分析升力翼翼剖面的上弧线和下弦线受压波动情况，重点以上弧线压力波动变化幅度、频率和峰值差为主要对比指标，分析确定波动变化幅度小、频率小和峰值差小的升力翼翼型方案III；
44）升力翼纵向对称面上流场结构比较分析：根据步骤3）所建立升力翼流体动力学计算模型，在同一速度等级条件下，模拟计算各型升力翼表面和纵向对称面上稳态压力分布及流动迹线，以升力翼上弧面尾部负压强度小、下弧面底部正压强度小、流场干扰效应小、涡量小及影响范围小为对比指标，确定升力翼翼型方案IV；
5）确定适配升力翼翼型结构：
分别依据步骤41）速度-气动拟合曲线比较分析优选的升力翼翼型方案I，步骤42）气动系数比较分析优选的升力翼翼型方案II，步骤43）升力翼表面气动压力系数比较分析优选的升力翼翼型方案III，步骤44）升力翼纵向对称面上流场结构比较分析优选的升力翼翼型方案IV，根据高速列车400~600km/h速度等级范围内升力翼的主要应用指标次序综合分析确定；其中升力翼的主要应用指标重要程度为：气动力＞气动系数＞升力翼表面气动压力系数＞升力翼纵向对称面上流场结构。
2.根据权利要求1所述的一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法，其特征在于：所述铁路建筑限界、桥隧限界条件为中国高速铁路标准限界范围，所述适用车型车辆限界为依据中国标准动车组横向截面轮廓的空间限界。
3.根据权利要求1所述的一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法，其特征在于：高速列车升力翼装置在限界横断内车顶上部所述最大有效安装空间（S1）为长方形空间，所述长方形空间关于纵向中心对称，其最大横向尺寸为3400mm。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310683954.6A CN116663455B (zh) | 2023-06-10 | 2023-06-10 | 一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310683954.6A CN116663455B (zh) | 2023-06-10 | 2023-06-10 | 一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116663455A CN116663455A (zh) | 2023-08-29 |
| CN116663455B true CN116663455B (zh) | 2023-10-13 |
Family
ID=87720448
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310683954.6A Active CN116663455B (zh) | 2023-06-10 | 2023-06-10 | 一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116663455B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115468732A (zh) * | 2022-07-26 | 2022-12-13 | 华设设计集团股份有限公司 | 一种高速列车升力翼安装布置及协同控制方法 |
| CN115476887A (zh) * | 2022-11-03 | 2022-12-16 | 兰州交通大学 | 一种高速列车升力翼 |
| CN115570983A (zh) * | 2022-10-26 | 2023-01-06 | 中南大学 | 一种列车用自适应升力受电弓 |
| CN115795642A (zh) * | 2022-09-08 | 2023-03-14 | 西北工业大学 | 一种高速列车专用翼型优化设计方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11295045B2 (en) * | 2019-04-03 | 2022-04-05 | GM Global Technology Operations LLC | Tools and methods for aerodynamically optimizing the geometry of vehicle bodies |
-
2023
- 2023-06-10 CN CN202310683954.6A patent/CN116663455B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115468732A (zh) * | 2022-07-26 | 2022-12-13 | 华设设计集团股份有限公司 | 一种高速列车升力翼安装布置及协同控制方法 |
| CN115795642A (zh) * | 2022-09-08 | 2023-03-14 | 西北工业大学 | 一种高速列车专用翼型优化设计方法及系统 |
| CN115570983A (zh) * | 2022-10-26 | 2023-01-06 | 中南大学 | 一种列车用自适应升力受电弓 |
| CN115476887A (zh) * | 2022-11-03 | 2022-12-16 | 兰州交通大学 | 一种高速列车升力翼 |
Non-Patent Citations (1)
| Title |
|---|
| 王瑞东 等.高速列车串列升力翼翼型优化设计.空气动力学学报.2021,130-136. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116663455A (zh) | 2023-08-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | A study of the influence of bogie cut outs' angles on the aerodynamic performance of a high-speed train | |
| Zhai et al. | Dynamics of high-speed train in crosswinds based on an air-train-track interaction model | |
| Guo et al. | Influences of affiliated components and train length on the train wind | |
| CN108583847B (zh) | 一种适用于长航时无人机的低雷诺数高功率因子翼型 | |
| CN104442859A (zh) | 带翼准悬浮与悬浮高速列车系统 | |
| Tian et al. | Flow structure around high-speed train in open air | |
| He et al. | Detached eddy simulation of a closely running lorry platoon | |
| CN108733914A (zh) | 基于人工神经网络的跨音速翼型自然层流延迟转捩设计方法 | |
| CN102951166B (zh) | 基于高速铁路的车体底盘带有仿机翼的空气动力悬浮列车 | |
| Liu et al. | Effects of bottom deflectors on aerodynamic drag reduction of a high-speed train | |
| Wu et al. | Aerodynamic drag performance analysis of different types of high-speed train pantograph fairing | |
| CN116663455B (zh) | 一种基于限界约束的高速列车升力翼翼型的计算流体动力学选择方法 | |
| Adamu et al. | An investigation of influence of windshield configuration and train length on high-speed train aerodynamic performance | |
| Xue et al. | Influence of variable cross-section on pressure transients and unsteady slipstream in a long tunnel when high-speed train passes through | |
| CN107264774B (zh) | 一种采用前缘支撑翼的m形翼高亚声速飞行器气动布局 | |
| CN115421426B (zh) | 一种高速列车侧翼升力调控装置安装布置及协同控制方法 | |
| CN116882311B (zh) | 高速列车升力翼常态化工作迎角的计算流体动力学确定方法 | |
| CN111859545A (zh) | 一种考虑升力匹配的宽速域高升阻比机翼优化设计方法 | |
| Xu et al. | Effect of two bogie cavity configurations on the underbody flow and near wake structures of a high-speed train | |
| CN111874029A (zh) | 一种降低气动噪声的高速列车排障器结构 | |
| CN105550465A (zh) | 一种基于橫风非定常理论的车厢横断面优化方法 | |
| CN115510710A (zh) | 尾部吹/吸气减阻式列车的气动阻力计算方法及系统 | |
| CN118133432B (zh) | 一种带翼稍小翼的高速列车升力翼优化设计方法及升力翼 | |
| Render et al. | Aerodynamics of battle-damaged finite-aspect-ratio wings | |
| CN112810639A (zh) | 一种高速列车减阻降噪的控制装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |