0　引言

普遍存在的自然风以及汽车会车、超车及通过隧道等特殊行驶工况的叠加和耦合常常会恶化汽车外流场，导致汽车所受气动力的变化更加复杂.汽车在侧风环境下会车时，随着两车距离的减小，周围的气流产生严重的气动干扰，引起空气绕流的变化.从而产生瞬时变化的气动力，直接影响汽车的阻力特性、行驶稳定性和安全性.因此，有必要对汽车在侧风环境下会车的气动特性开展相关研究.目前国内外对于车辆交会的空气动力特性研究主要是集中于轨道车辆，而对于汽车这一类公路运输车辆的研究较少.高磊等［1］对会车过程进行三维瞬态数值模拟，得到了会车过程的瞬态空气动力特性，并且研究了不同相对车速对会车过程汽车瞬态空气动力特性的影响.张英朝等［2］使用商业的CFD软件，对简单外形汽车隧道中会车过程的瞬态空气动力学特性进行了模拟.徐国英等［3］运用计算流体力学数值模拟方法对在一级公路和四级公路上行驶的汽车迎面会车过程进行了仿真研究.上述研究主要是针对非侧风情况下汽车会车的气动特性，而对侧风下会车工况的气动特性研究很少.

1　数值模型及计算工况

为简化计算，侧风下两车等速会车的气动特性研究只考虑二维情形.对于两车交会的相对运动采用滑移网格技术处理，因此，在图1所示计算域中将整个流动区域划分为4个区域（图中L为车长）.共划分为两个移动区域和两个静止区域，移动区域与相邻区域通过交界面来交换数据.移动区域01（面EFJK）中包含第一辆汽车car01，其移动速度为V1；移动区域02（面CDLM）中包含第二辆汽车car02，其移动速度为V2.流动计算的car01和car02相隔1.2 m，且外形一致（图2），其中L、W、R分别为车长、车宽、圆角半径（mm）.考虑为等速（V1=V2=25 m/s）会车，侧风风速为10 m/s，两车相对速度计算的流动雷诺数和马赫数分别为5.6E6、0.147.因此，可将汽车交会的外流场简化为二维、粘性、不可压缩、不考虑热量交换的非定常湍流流动.在计算模拟中选用可实现的k-ε湍流模型，通过求解连续性方程、动量方程、湍动能方程和耗散率方程得到整个流场的数值解［4-6］.整个计算域采用结构化四边形网格进行离散（图3），车体表面加有细密的边界层网格.

  

图1　流动计算域Fig.1 Computing domain of flow field

  

图2　汽车二维CAD模型（cm）Fig.2 Two-dimension CAD model of vehicle（cm）

  

图3　车体周围网格Fig.3 Mesh around the body

为了得到具体流动问题的特解，表1中列出了相应的流动边界条件.在非定常计算中，其时间步长为0.001 s，单步中迭代次数为20，计算时间与真实的流动时间一致.

 

表1　流动边界条件Tab.1 Boundary conditions of the flowing

  

几何边边界条件HI速度进口CM、DL、EK、FJ交界面AB、AN、BC、HI、GH、IJ、ML、EF压力出口

2　仿真结果及分析

2.1　交会气动特性

在数值计算处理中，参考面积为1.8 m2，参考速度为51 m/s，特征长度取车宽1.8 m，垂向长度按1 m计算.通过非定常计算可以获得侧风下交会两车的气动特性时程变化趋势.图4、图5分别交会时两车的阻力系数及侧向力系数随时间变化的历程，由于气动干扰，交会时两车的气动特性表现出较明显的非定常性.由图4的阻力系数变化可以看出，处于交会侧风上游的汽车car01的阻力系数变化较平缓，而处于下游的汽车car02的阻力系数变化较剧烈.上游的car01的阻力系数随着两车逼近逐渐减小，当会车后两车头部间隔△d=4.0 L时，阻力系数基本保持不变.下游的car02的阻力系数经历先增加，后减小，再增加，再减小的变化过程.图5中也反映出汽车car02的侧向力系数较car01的侧向力系数变化更复杂.当会车的两车头部平齐时，car02的侧向力达到第一个峰值点，而后逐渐下降，再继续上升到第二个峰值点.随着会车过程的继续进行，侧向力系数逐渐减小，当△d=2.0 L时，侧向力出现了负值，即由原来所受的排斥力变化为吸引力.处于侧风上游的汽车car01的侧向力系数首先逐渐增加，再减小，而后缓慢增加直至趋于基本不变.

施工前，应先对施工现场地形进行测量，划分出基坑支护高度分段；在分段高度范围内，最深点对应剖面设计参数，深度变浅时，土钉需按照先减上、后减下的方式逐步减少。

图11、图12分别为t=2.00 s时刻会车车辆car01和car02沿x方向、y方向的静压系数比较.由图中曲线的积分面积可以看出：t=2.00 s时刻car01和car02所受到的侧向力相差较大，而阻力的大小很接近.图13、图14分别为t=2.25 s时刻会车车辆car01和car02沿x方向、y方向的静压系数比较.在t=2.25 s时刻，car02沿x方向的积分面积明显大于car01，但沿y方向的则相反.因此，该时刻下，car01受到的阻力略大于car02，而侧向力则明显小于car02.

  

图4　阻力系数变化Fig.4 Variation of dragging force coefficient

  

图5　侧向力系数变化Fig.5 Variation of side force coefficient

2.2　特定时刻流场分析

通过对侧风下会车的两车气动特性的数值研究，可以获得以下几点主要结论：

  

图6　各时刻两车相对位置Fig.6 Relative place of two passing vehicles at the different moment

1）在侧风环境下，等速会车的两车均表现出较强的非定常性，阻力系数和侧向力系数均呈现出一定的非周期性变化.

  

图7　t=1.90 s时刻Fig.7 At the moment of t=1.90 s

  

图8　t=2.05 s时刻Fig.8 At the moment of t=2.05 s

  

图9　t=2.20 s时刻Fig.9 At the moment of t=2.20 s

  

图10　t=2.30 s时刻Fig.10 At the moment of t=2.30 s

世界第四大铁矿生产商，澳大利亚的福蒂丘集团也在厄瓜多尔设立了办公室，该公司承诺今年投资1亿美元进行勘探。

  

图11　沿x方向的静压系数比较Fig.11 Comparison of pressure coefficient in the x direction

  

图12　沿y方向的静压系数比较Fig.12 Comparison of pressure coefficient in the y direction

  

图13　沿x方向的静压系数比较Fig.13 Comparison of pressure coefficient in the x direction

  

图14　沿y方向的静压系数比较Fig.14 Comparison of pressure coefficient in the y direction

相会的两车在不同时刻（t=2.00 s、t=2.15 s），其气动特性也发生明显变化.图15、图16分别为汽车car01沿x方向和y方向的静压系数比较.图17、图18分别为汽车car02沿x方向和y方向的静压系数比较.在t=2.00 s时刻，两车车头正好重叠，随着会车过程的继续，直到t=2.15 s时刻，此时两车头部间距1.5倍车长，一半的车身产生重叠.相比t=2.00 s时刻，t=2.15 s时刻的气动干扰更加剧烈，大大削弱了两车所受的气动力.因此，无论是car01还是car02，后一时刻所受到的阻力和侧向力均小于前一时刻.

  

图15　沿x方向的静压系数比较Fig.15 Comparison of pressure coefficient in the x direction

  

图16　沿y方向的静压系数比较Fig.16 Comparison of pressure coefficient in the y direction

  

图17　沿x方向的静压系数比较Fig.17 Comparison of pressure coefficient in the x direction

  

图18　沿y方向的静压系数比较Fig.18 Comparison of pressure coefficient in the y direction

3　结论

为了理解以上的汽车在侧风下交会时气动干扰特性，可以针对特定时刻的流场进行单独分析.图6为相会两车的几种典型位置关系，其中△d表示两车头部的间隔距离，L为车长，t为给定时刻.

图7—图10为4个特定时刻t=1.90 s、t=2.05 s、t=2.20 s、t=2.30 s的流场压力云图和流线图.从这些压力分布和流线图中可以看出，在侧风和汽车行驶风的联合作用下，在车体背风侧和尾部不同程度地出现了负压区.特别是在两车相隔间距较近的时刻，复杂的压力分布直接导致作用在车体上的气动力发生显著的变化.在t=2.05 s时刻，两车头部相隔0.5倍车长，车身的一半重叠在一起，在car02的迎风侧正压区直接影响到car01背风侧的负压区，因此，该时刻car01所受到的侧向力明显小于car02.在t=2.20 s时刻，两车头部间距为2倍间距，车身正好错开.在car02的背风侧出现了一个很明显的负压区，由于car02尾流的影响，在car01的背风侧尚未出现明显的负压区，因此，该时刻car02所受到的侧向力明显大于car01.

2）侧风会车时，两车的气动干扰对两车的气动特性产生明显的影响，但其中对于侧向力的影响更加剧烈.

作为老牌的教学医院，齐鲁医院有自己的胸怀和担当。陈玉国表示，在齐鲁医院的文化氛围中，学生会在住培过程中被熏陶出良好的行为和素养，将来无论在哪里，都会携带齐鲁医院优秀的文化和基因。

从表2中的公式所得纵向抗剪承载力可以看出，波形PBL连接件的抗剪承载力随着贯穿钢筋直径的增加而提高.但根据试验结果来看，波形PBL连接件的抗剪承载力虽在提高但之后的贯穿钢筋直径的变化却对承载力提高的不大.分析原因可能是公式中并未考虑到开孔中粗骨料的情况，当开孔孔径较小，而贯穿钢筋直径较大时，粒径较大的粗骨料不能进入孔径中，这种情况大大降低了混凝土榫的承载能力.因此在采用公式(1)计算PBL连接件抗剪承载力应该先考虑贯穿钢筋截面面积在开孔面积的占比再考虑PBL连接件抗剪承载力的计算公式.

3）简化的二维数值模拟与真实的会车会有一定的偏差，但仍可为后续的真实车体会车绕流流动的气动分析提供借鉴和参考.

参考文献

［1］高磊，谢金法.车速对会车车辆瞬态气动特性的影响［J］.农业装备技术，2009，35（6）：15-17.

［2］张英朝，傅立敏.简单外形汽车隧道中会车过程的瞬态空气动力学数值模拟［J］.吉林大学学报（工学版），2006，36（3）：302-306.

［3］徐国英，卢军，薛劲橹.基于CFD的汽车迎面会车过程气动干扰特性仿真研究［J］.测试技术学报，2010，24（6）：502-506.

［4］罗建斌，胡爱军.高架桥声屏障高度对列车气动特性影响的数值模拟［J］.计算物理，2012，29（1）：65-72.

［5］王福军.计算流体动力学分析［M］.北京：清华大学出版社，2006.

［6］罗建斌，胡爱军，粟腾超.声屏障高度影响列车气动特性的数值研究［J］.广西科技大学学报，2014，25（3）：14-19.