随着汽车速度的不断提高，随之出现的与汽车气动性能有关的问题也越来越多，如空气阻力等。优化汽车外形，减少空气阻力，有限元分析是一个行之有效的方法。

从时速不超过50km/h的马车型汽车，发展到现在的流线型汽车时速可以达到200～300km/h，百年的汽车发展给人类社会带来了巨大而深刻的变化。随着汽车速度的不断提高，随之出现的与汽车气动性能有关的问题也越来越多，如空气阻力（当车速达到250～300km/h时，空气阻力占总阻力的50%以上）、气动噪声以及对周围环境的影响等，均与汽车外形的流线程度有关。为此，还必须进一步优化车体外形，减少空气阻力等对汽车的影响。有限元分析计算是研究汽车气动性的一种手段，本文运用ANSYS公司的有限元分析软件ANSYS对三种汽车外形进行了分析研究，针对三种汽车外形气动性能做了比较，研究车外形对气动性能的影响。

数学以及有限元模型
根据流体力学理论，计算的数学模型采用三维粘性不可压缩雷诺平均应力方程和工程上应用较广的k-∈双方程湍流模型。
连续性方程：

动量（矢量）方程：

式中me??为有效粘性系数，且me?? =m+mt，m为动力粘性系数，mt为湍流粘性系数，其中mt =Cmρκ2/∈，Cm为常数，一般取为Cm=0.09；P'为修正压力，且P'=P+—ρκ；k和∈值1可由k-∈双方程得出。

湍流动能k方程：

湍流动能耗散率方程： 

其中C∈1， C∈2， σκ， σ∈为常数，其取值一般为：C∈1=1.44，C∈2=1.92，σκ=1.0，σ∈=1.3。剪切生成项G的取值一般按下式得出

由于本文所使用的是ANSYS教学版软件，在计算规模上有很大限制，因此对计算模型进行了必要的简化，考虑汽车结构的对称性，在忽略侧风等非对称性因素的影响，考虑汽车的对称面，同时对车底、车头等部位均进行了必要的简化，建立了平面有限元模型。使用ANSYS软件的2D FLOUNT 141单元进行离散，用有限元的二维分析模型进行求解。

计算与分析

在此，针对不同时代的车型选用了三种有代表性的车型进行分析。

1.车型1，（如图1～4）

图1  速度为50km/h

 
图2  流体总体压力分布图

 
图3  速度为250km/h

图4  流体总体压力分布图

2.车型2，（如图5～8）

图5  速度为50km/h

图6  流体总体压力分布图

图7  速度为250km/h

图8  流体总体压力分布图

3.车型3，（图9～12）

图9  速度为50km/h

 
图10 流体总体压力分布图

 
图11 速度为250km/h

图12 流体总体压力分布图

4.速度分布

图1、3、5、7、9、11给出的是3种汽车分别在高低速时纵向对称面流场速度分布情况。可以看出，由于三种头部形状流线型都较好，当气流流过汽车头部时基本上未出现大的流动分离。但车型3流线型最佳且外形曲面变化最小，所以效果最好。

气流在经过汽车顶部时达到最高速度，然后气流回落速度也逐渐降低。在汽车尾部，三种车型都产生了涡流情况。随着离汽车尾部距离的增大，旋涡的涡强在减弱，涡核向外移动，并向地面靠近。这说明涡流在不断地扩散，能量也随之在耗散。车型3有较好的表现，尾部有较好的流线型，且经过了车型的改进使得涡流较小。而在高速时车型3的表现更突出，其他两种车在高速时，尾部涡流有变大的趋势，且流线型变差速度变化较大。

空气压力分布

图2、4、6、8、10为气流总体压力分布图。由于汽车底部靠近地面，把气流进行了压缩增速，所以，3种车都出现了不同层度的高压区域。相比较而言，车型3稍微好于其他两种。由于头部气流开始与汽车接触变化，气压也开始变大。车型1较差，在汽车前窗部分出现较大的高压区域，而其他两种较好。在尾部，由于涡流的产生等，出现了负压区域。在后上方，由于车顶部高速气流的作用，产生了一条较长的较高的压力带，且逐渐变小。

结论

由ANSYS分析计算得知：在汽车头部气流开始分流，而汽车外形的变化使流速开始加快。在汽车的顶部和底部产生高速区，且在车顶部分产生最高流速，然后开始逐渐回落。汽车尾流中有纵向的涡流产生，随尾流的发展，涡流由强到弱，涡心降低（逐渐接近轨面），涡核向外移动。通过比较，第三种车型优于其他两种车型，且在低速和高速时都有较好的空气动力性能。故在汽车外形发展设计研究时，应尽量将汽车外形设计成流线型，外形曲面曲率不应出现突变或变化过大，避免产生流动分离。