涡格法.docVIP

螺旋桨涡格法的初步实现 螺旋桨升力面理论有涡分布法，偶极子分布法和加速度势法等。经过多年的应用实践，现在普遍采用的是涡分布法。虽然面元法的发展比升力面方法在理论上更为完善，尤其是压力分布的预报计算显著优于后者，但就螺旋桨合力的计算而言，两者基本具有同样的计算精度，因此升力面理论方法目前仍是螺旋桨理论设计所采用的主要工具。随着计算机速度和内存的不断提高和扩大，涡格法（Vortex Lattice Method，即VLM）得到越来越广泛的应用。至今，各学者所发表的各种涡升力面理论的论文，就原理而言都是一致的，而主要差别仅在于涡模型包括尾涡处理和数值方法两方面。 作为螺旋桨理论的初学者，为了深入理解升力面理论以及应用方式，我们将通过自编程序在计算机上进行简单的涡格法实现。由于螺旋桨的几何形状较为复杂，而我们的目的在于对升力面涡格法有初步的认识，因此我们对螺旋桨作一定程度上的简化，既能反映涡格法的计算特征，又能使计算相对简便。在对涡格法理论模型进行简单的介绍之后，本文将给出涡格法的实现程序，并对不同攻角以及不同拱弧面下的计算结果进行分析比较。 一、涡格法理论模型 1、螺旋桨几何特征 本文使用一空间曲面（拱弧面，即升力面）代替螺旋桨桨叶，桨叶厚度可由离散的源汇分布来模拟，但在本文中不予考虑。该拱弧面在xoy平面的投影为展长2.0，弦长1.0的矩形平面，在xoz平面上的投影为一剖面拱弧。为了比较不同的拱弧线对桨叶性能的影响，本文使用三种不同的剖面拱弧线，它们分别为NACA剖面拱弧、圆弧以及抛物线拱弧。为了比较不同的最大拱高值对桨叶性能的影响，本文取两种不同的最大拱高，分别为fmax=0.02*c和fmax=0.04*c，其中c为弦长1.0。 2、控制方程 涡格法的思想是将升力面设置在拱弧面上，然后在升力面上进行涡格划分，并在每个涡格上布置马蹄涡和控制点，而每个马蹄涡的强度是未知的，可以通过在控制点处满足拱弧面上法向速度为零的边界条件进行求解。本文中涡格为均匀网格划分，展向划分m个涡格，弦向划分n个涡格，在每个涡格中将马蹄涡的展向涡段布置在每格的1/4弦向格长位置处，而控制点置于3/4弦向格长及1/2展向格长位置，Weissinger证明了这样做已隐含着库塔条件的满足。对于弦向涡线则沿涡格线分成直线涡段，亦即说整个弦向涡的几何形状用折线代替，尾涡区涡线另行处理。这样，对于每一个控制点，满足物面不可穿透条件可得： 其中和分别表示第个涡格控制点的诱导速度和法向矢量，表示流场的入流速度。为所有的升力面马蹄涡在第个涡格控制点的诱导速度的叠加，即： 其中，，分别表示单位涡强的第个马蹄涡的展向涡段、弦向涡段以及尾涡段对第个涡格控制点的诱导速度，可根据Biot-Savart定理求得，而表示第个马蹄涡的涡强，为待求的未知量。 对所有的控制点满足方程联立求解，即可得每个涡格内马蹄涡的涡强。 3、叶面区的涡系模型 在叶面区内，展向均匀划分15个格，弦向均匀划分20格，而每一个马蹄涡按涡格划分为展向涡段和弦向涡段。由于NACA的剖面拱弧型值为一系列散点坐标，所以在求解涡段端点及控制点坐标时需要进行插值处理。考虑到最大拱高值相对于弦长较小，使用线性插值和三次样条插值的结果差别不大，因此选用较为简单的线性插值即可。对于控制点的法向矢量，通过求取与控制点相邻两点所组成的直线斜率，然后再作简单的三角处理即可求得。 对于圆弧以及抛物线拱弧，可直接求得各点的坐标值和法向矢量。但为了保持程序的统一性，本文依然先求得拱弧线的一系列散点坐标，然后按与NACA拱弧相类似的处理方法求取各点的坐标值和法向矢量。 4、尾流区的涡系模型 尾流区内尾涡的处理较为复杂，且对计算结果有较大的影响。尾流区是自由空间，故涡的走向须沿着流场内该点的速度方向。而流场速度是需问题解出后才知道，故尾流区是事先未知的，这就构成了解升力面问题的非线性。 为了使求解问题简化，本文事先假定了尾流区的范围和形状。本文将尾流区分为过渡尾流区及远尾流区，以便更好的符合实际情况。在过渡尾流区，尾涡在随边处沿剖面拱弧的切线方向泄出，由于受来流速度的影响，尾涡在过渡尾流区结束的地方与来流方向一致。过渡尾流区的长度为一倍弦长。在过渡尾流区内，尾涡的走向成二次抛物线型。与叶面区类似，过渡尾流区内的尾涡按弦向共划分n个涡段。同时，考虑到展向两端的尾涡卷曲，尾涡在随边处的下泄方向不仅与剖面拱弧的切线方向有关，而且与来流速度相关，因此本文中尾涡的下泄方向为： 其中，，分别表示尾涡下泄方向，来流速度方向和随边处剖面拱弧的切线方向，与展向坐标有关，其表达式为： 其中为尾涡下泄点的展向坐标，为展长。 在远尾流区，尾涡方向与来流速度方向一致，为半无限长涡线，处理较为简单。 二、数值结果的比较分析 本文计算了NACA剖面拱弧、圆弧