AVL（Athena Vortex Lattice）如何编写与应用（一）

AVL，全称为Athena Vortex Lattice，是MIT的Mark Drela教授开发的一个气动分析程序。程序最初由Harold Youngren 1988年为MIT Athena TODOR航空软件集编写。经过Mark Drela和Harold Youngren大量修改后，现在已经发展到3.36版了。AVL的启动分析部分用FORTRAN编写而成，图形显示则是基于XWindows，用C语言编写，正是因为基于这两种以速度著称的语言，AVL运算起来特别快。而且它的跨平台性也比较好，在UNIX和Windows系统下都能运行。

原文链接：https://blog.csdn.net/weixin_41767350/article/details/79564733

感觉没有一个很详细的文章介绍怎么使用，研究了一下说明书和学长，浅浅写一下相关参数的作用。

一、xxx.avl文件

这个文件是用来建立飞机模型的文件，是一定要有的文件，其他的.run和.mass并不是必须要有的，甚至mass文件并没有看到有人在使用。

基本的文件格式如下：

Plane 1
0.0                   !   Mach
0     0     0.0       !   iYsym  iZsym  Zsym
0.190 0.239 0.800      !   Sref   Cref   Bref   reference area, chord, span
0.05934 0.0 -0.0004924 !   Xref   Yref   Zref   moment reference location (arb.)

#
#==============================================================
#
SURFACE
Main Wing  
# Main wing
5  1.0  17  1.0   ! Nchord   Cspace   Nspan  Sspace
#
# reflect image wing about y=0 plane
YDUPLICATE
     0.00000 
#
# twist angle bias for whole surface
ANGLE
     0.00000    

#
# x,y,z bias for whole surface
TRANSLATE
    0.0  0.00000     0.00000

#Here the Section Start 
#----------------------- Outter section------------------------
#
#    Xle         Yle         Zle         chord       angle   
SECTION
    0       0.400        0       0.20       0.000  

AFIL
naca0012E.dat


CONTROL
aileron  1.0   0.658    0.0 1.0 0.0   -1.0

#------------------------Inner section-------------------------
#--------------------------------------------------------------
#    Xle         Yle         Zle         chord       angle   
SECTION
     0        0.0        0.0         0.275        0.000   

AFIL
naca0012E.dat

CONTROL
aileron  1.0   0.60    0.0 1.0 0.0   -1.0



#Wing Section defined

#==============================================================
#Horizontal Stabilizer
SURFACE
H-Stab
4  1.0  9  -1.5  !  Nchord   Cspace
#
# reflect image wing about y=0 plane
YDUPLICATE
     0.00000 
#
# twist angle bias for whole surface
ANGLE
     0.00000    
#
# x,y,z bias for whole surface
TRANSLATE
   0.56    0.00000     0.00000
#--------------------Section 2--------------------------------
SECTION
   0.015      0.22        0.0         0.13        0.000   

AFIL
naca0012E.dat

CONTROL
elevator  1.0    0.45    0.0 1.0 0.0     1.0

#Horizontal Stabilizer Defined

#--------------------Inner 1---------------------------------
#    Xle         Yle         Zle         chord       angle  
SECTION
   0          0.0        0.0         0.145        0.000   

AFIL
naca0012E.dat

CONTROL
elevator  1.0    0.41    0.0 1.0 0.0     1.0

#==============================================================
#Vertial Stabilizer
SURFACE
V-Stab
4  1.0  15  0.75  ! Nchord   Cspace
#
# x,y,z bias for whole surface
TRANSLATE
    0.65    0.00000     0.00000
#--------------------Section 1---------------------------------
#    Xle         Yle         Zle         chord       angle   
SECTION
   0           0.0       0        0.15        0.000   

AFIL
naca0012.dat  

#--------------------Section 2---------------------------------
SECTION
   0.040       0.0        0.18         0.09        0.000   

AFIL
naca0012.dat
#==============================================================
#Vertial Stabilizer Defined

首先第一行是飞机的名称，并没有什么讲究就是要写上而已，和文件名并没有关联。

第二行是mach 默认自由流马赫数，在这里默认为零，对气动计算影响不大。

第三行的三个参数：实际上也可以都设为0，不用改变

        iYsym：在-1,01之间选择，等于1是飞机关于y=0对称等等

Sref、Cref、Bref是飞机的参考面积参考弦长和参考的展长，我的理解就是所设计飞机的S，c，b这三个参数。

XYZref：力矩和转速被定义的相关默认位置，就是重心位置，且配平计算时必须是CG位置。

需要注意的是这里定义的Mach，XYZref值可以在运行中改变。

———————————————以上是飞机的基本信息———————————————

在此之后设置飞机的机翼：

SURFACE
Main Wing  
# Main wing
5  1.0  17  1.0   ! Nchord   Cspace   Nspan  Sspace
#
# reflect image wing about y=0 plane
YDUPLICATE
     0.00000 
#
# twist angle bias for whole surface
ANGLE
     0.00000    

#
# x,y,z bias for whole surface
TRANSLATE
    0.0  0.00000     0.00000

#Here the Section Start 
#----------------------- Outter section------------------------
#
#    Xle         Yle         Zle         chord       angle   
SECTION
    0       0.400        0       0.20       0.000  

AFIL
naca0012E.dat


CONTROL
aileron  1.0   0.658    0.0 1.0 0.0   -1.0

#------------------------Inner section-------------------------
#--------------------------------------------------------------
#    Xle         Yle         Zle         chord       angle   
SECTION
     0        0.0        0.0         0.275        0.000   

AFIL
naca0012E.dat

CONTROL
aileron  1.0   0.60    0.0 1.0 0.0   -1.0

每一个SURFACE对应一个翼面直到下一个surface或body出现

SURFACE关键字下面一行是这个翼面的名字

在此之后设置的是马蹄涡的个数，这个数量没有一个固定的要求，不需要太多。

前两个是弦向的涡个数和以cos函数分布的间距参数（设为1即可）后两个可以在这里声明也可以写到后面的每一个section中，是展向的涡个数和以sin函数分布的间距参数。

当声明了YDUPLICATE这个关键字，相当于飞机是轴对称的，它的对称面就是y=关键词下的数字。

需要注意的是，只有当定义了iYsym = 0的时候才可以使用，否则会出错。这也是为什么在最开始不用额外设计那三个量的原因。

ANGLE是入射角对应关键字，更改入射角而无需更改每一个section的Ainc值，参数表示添加到所有定义的Ainc上的偏移量。

在此之后最重要的就是SECTION关键词

它定义了这个型号的机翼的位置弦长入射角等。

SECTION                             |  (keyword)
0.0 5.0 0.2   0.50  1.50   5 -2.0   | Xle Yle Zle   Chord Ainc   [ Nspan Sspace ]

需要注意的是定义翼面的时候一定要从左到右来定义

垂尾从上到下排序是最方便的。

另外所有的翼型文件一定要和avl文件保存在一个目录下才能使用。

之后的control和操作步骤问题之后再说，先说这么多吧~