AVL教学

  • 格式:doc
  • 大小:895.00 KB
  • 文档页数:12

A VL使用说明1.什么是A VLA VL (Athena V ortex Lattice)是MIT的Mark Drela教授开发的一个气动分析程序。

程序最初由Harold Youngren 1988年为MIT Athena TODOR 航空软件集编写。

经过Mark Drela 和Harold Youngren大量修改后,现在已发展到3.15版了。

A VL 的气动分析部分用FORTRAN编写而成,图形显示则是基于XWindows,用C语言编写,正是因为基于这两种以速度著称的语言,A VL运算起来特别快。

而且它的的跨平台性也比较好,在UNIX和Windows系统下都能运行。

2.A VL能做什么A VL利用涡格法对给定的外形作气动分析,采用关键字描述几何外形,计算出数据并绘制图形,而且数据和图形都能够输出。

这在计算后利用其结果十分方便。

A VL具有为飞机外型快速分析的大量特性,主要的特征如下:气动分力升力面细长体外型描述关键字驱动的几何外型输入文件使用线形插值描述截面截面特性翼型可以是NACA xxxx, 或者来自于翼型输入文件操纵面偏转抛物线形剖面极线,尺度放缩对整个表面或机身做缩放、平移、旋转复制整个表面或机身特性马蹄涡(表面)源线+双合线(机身)有限主体选择离散化归一化正弦余弦组合操纵面偏转通过倾斜法向量前缘襟翼后缘襟翼铰轴线独立于离散过程普通自由流描述alpha, beta 气流角p,q,r 飞机旋转分力亚音速Prandtl-Glauert 可压缩流处理空气动力分力输出直接力和力矩特瑞夫兹平面(Trefftz-plane)力及力矩的导数,w.r.t自由流,旋转,控制力在体坐标或固定坐标下配平计算控制变量alpha, betap, q, r操纵面偏转约束对变量的直接约束通过指定升力系数、力矩来间接约束一次定义多个配平算例保存配平算例设置以备以后调用可选质量定义文件(仅在配平设置,模态分析中)用户选择单位逐条列出部件位置,质量和惯性约束下的配平设置平直或带斜坡的水平飞行固定速率的俯仰飞行模态分析准稳定模型下的刚体分析显示带参数的特征根级数实时显示模态运动输出动态系统矩阵3.涡格法模型使用的原则就如同其它数值算法一样,A VL的用处同样有所限制。

在使用时必须牢记。

1)外型涡格法主要适合于小攻角侧滑下的薄升力面组成的气动外形。

外型表面及其拖曳出的尾迹构成单层的涡流面,可以离散成马蹄涡线,这些流线平行于x轴。

A VL也能由源线和双合线模拟如机身或者发动机舱这类的细长体。

合力和合力矩符合细长体理论。

但对这种模型的认识相对有限,所以用它建模时尤其要注意。

如果要求机身对气动载荷影响很小,最好还是不要使用A VL模型。

2)非定常流A VL采取准稳定流计算,也就是非稳定的脱体涡会被忽略。

更准确地说,它假设减小的频率受到限制,摆动(例如倾斜)必须足够慢,使得摆动的周期比流过机翼弦线的时间要长得多。

这在实际上也跟任何期望的机动飞行相符。

同样,用于计算的滚转、俯仰、偏航的速率必须很小而使气流角很小。

它们可以用无量纲的滚转速率参数判断,这些量限制在如下经验范围:-0.10 < pb/2V < 0.10-0.03 < qc/2V < 0.03-0.25 < rb/2V < 0.25这些数据表明非常猛烈的机动飞行除了可能出现在低速飞行特技表演之外,一般都不会包括在任何典型的飞行任务中。

任何情况下,如果其中的任何参数超过上述范围,结果都需要特别注意小心处理。

3)可压缩性可压缩性使用Prandtl-Glauert (PG) 变换处理,其相对重要性由PG因子1/B = 1 / sqrt(1 - M^2) 。

其中“M”是自由流的马赫数。

下表给出期望的有效范围:M 1/B___ _____0.0 1.000 |0.1 1.005 |0.2 1.021 |0.3 1.048 |- PG有效值0.4 1.091 |0.5 1.115 |0.6 1.250 |0.7 1.400 PG怀疑范围(气流可能接近音速)0.8 1.667 PG不可靠范围(气流肯定是音速)0.9 2.294 PG无效范围对后掠翼的外型,PG模型的有效性适合用垂直翼马赫数衡量Mperp = M cos(sweep)只要Mperp < M,建模时后掠翼型可以比非后掠性适用于更高的马赫数。

例如,45度后掠翼在自由流M = 0.8时Mperp = 0.8 * cos(45) = 0.566仍然位于上面给出的PG有效范围内。

这种情况下一般能给出合理的结果。

如果速率参数就位于最低马赫数,即M = 0.2以下,最好就设M = 0。

这样可以加速计算,因为改变马赫数需要对VL影响线矩阵重新计算和重新分解因式,这些运算会抵消大部分的计算效果。

如果马赫数保持不变,计算仅一次就能完成。

4.输入文件A VL采用三种输入文件,它们都是纯文本格式。

每个算例三个输入文件都需要同样的文件名“xxx”,而每个的扩展名及说明如下:xxx.avl 必须的主要输入文件,定义几何外形。

xxx.mass 可选的输入文件。

给出质量和惯性及尺寸单位。

xxx.run 可选的输入文件。

定义算例的一些参数。

用户提供文件xxx.avl 和xxx.mass,这些文件一般由任何文本编辑器就可以创建。

给出的算例文件可以作为使用的模板。

xxx.run 文件是A VL在用户给的指令下创建的,它可以手动修改,但一般不需要这样做,更方便的办法是在A VL 中改变它的内容然后重新写这个文件。

输入文件数据的具体格式参见A VL的使用说明。

5.飞机的气动计算过程1)先建立飞机的外形文件,由于是低速飞机,所以M数取0。

下面给出我们飞机建立模型的数据。

文件:bs.avlPlane LanJian#Mach0.0#iYsym iZsym Zsym0 0 0#Sref Cref Bref19.3 1.7545 11.0#Xref Yref Zref0.0 0.0 0.0#==================================================================== #==================================================================== SURFACEMain Wing# Nchordwise Cspace10 1.0YDUPLICATE0.0ANGLE2.0# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace0.0 0.0 0.0 1.80 0.0 20 1.0AFILEepplear625.afoil# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace1.7 1.0 0.02.40 0.0 10 1.0AFILEepplear625.afoil# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace1.7 1.50 0.02.40 0.0 30 1.0AFILEepplear625.afoilCONTROL# Cname Cgain Xhinge HingeV ec SgnDupaileron -1.0 0.8 0.0 0.0 0.0 -1.0# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace4.055.5 0.0 1.2 0.0 0 0AFILEepplear625.afoilCONTROL# Cname Cgain Xhinge HingeV ec SgnDupaileron -1.0 0.8 0.0 0.0 0.0 -1.0#==================================================================== #==================================================================== SURFACEV-stab# Nchordwise Cspace Nspanwise Sspace10 1.0 20 1.0YDUPLICATE0.0TRANSLATE4.2 1.45 0.2# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace0.0 0.0 0.0 1.3 0. 0 0NACA0012CONTROL# Cname Cgain Xhinge HingeV ec SgnDuprudder 1.0 0.72 0.0 0.0 1.0 1.0# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace0.4 0.0 1.2 0.9 0. 0 0NACA0012CONTROL# Cname Cgain Xhinge HingeV ec SgnDuprudder 1.0 0.72 0.0 0.0 1.0 1.0#==================================================================== #==================================================================== SURFACEH-stab# Nchordwise Cspace Nspanwise Sspace10 1.0 20 1.0YDUPLICATE0.0TRANSLATE4.4 0.0 1.46# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace0.0 0.0 0.0 1.1 0. 0 0NACA0008CONTROL# Cname Cgain Xhinge HingeV ec SgnDupelevator 1.0 0.72 0.0 1.0 0.0 1.0# -------------------------------------------------------------SECTION# Xle Yle Zle Chord Ainc Nspanwise Sspace0.0 2.0 0.0 1.1 0. 0 0NACA0008CONTROL# Cname Cgain Xhinge HingeV ec SgnDupelevator 1.0 0.72 0.0 1.0 0.0 1.0# -------------------------------------------------------------# -------------------------------------------------------------关于长度的数据都是按分米取值。