参考文献-中国力学学会.docVIP

基于动量源方法的直升机旋翼/机身流场数值模拟 王博 招启军南京航空航天大学[1]。这类方法不太适合于分析粘性影响较大、甚至存在气流分离现象的旋翼/机身干扰流场；另一类方法是新兴的计算流体力学（CFD）分析方法[2][3]。由于Euler方程和Navier-Stokes（N-S）方程从反映流体本质的基本方程出发对旋翼/机身流场进行求解，因而能充分地捕捉该干扰流场的细节特征和反映粘性对流场特性的影响，但由于方法本身的复杂性及受计算机资源的限制，其研究还远未成熟。发展有效的、低成本和短周期的旋翼/机身流场的分析方法是当前直升机CFD的一个重要发展方向。本文拟从研究复杂的直升机旋翼流场入手，建立一个基于动量源方法的有效和工程实用的旋翼/机身组合流场的数值模拟方法，开展直升机旋翼/机身组合流场数值模拟及其设计参数的影响分析。 国内关于旋翼/机身流场模拟的研究主要是基于涡流理论，例如文献[4][5]，CFD方法的研究还主要集中于孤立旋翼流场的模拟[6][7]。 2动量源方法 为了降低计算代价和提高模拟旋翼下洗流场的效率Rajagopalan[8] 和Chaffin[9]采用以桨盘代替旋翼的办法对旋翼流场进行了计算。采用动量源方法最大的特点是把桨叶对气流的作用以动量源的形式来代表，并忽略桨叶附近的流动细节特征，把周期性的流动通过时间平均的办法转化为“准定常”流动。采用上述方法在保证旋翼下洗流场的本质属性的前提下，舍去求解旋翼桨叶周围流场的流动细节给计算精度提出的要求。同时，由于用围绕整个桨盘的网格来取代围绕桨叶的贴体网格，减小了网格生成的难度和网格数目，从而有效地节省了计算时间。因此，在许多不需要关注桨盘表面流动细节和对计算效率要求较高的情况下，采用动量源方法可以获得很好的效果。 2.1控制方程 由于旋翼下洗流场的速度不高，将整个流场看作不可压有粘流动，控制方程采用定常不可压层流N-S方程，其中X向动量方程如下 X向动量方程： （2） 其中、是动量源项，分别表示某一点上的旋翼诱导力与单位体积之比。Y向、Z向动量方程与此类似。 2.2桨盘气动模型，桨盘直角坐标系为，二者关系如图1所示，转换关系如下 （3） 其中，A、B分别表示桨盘左倾角和前倾角，是旋翼旋转中心在计算域内的坐标，当悬停时，A和B都为0。 图1 计算域直角坐标系和桨盘直角坐标系 图2 桨盘直角坐标系与桨盘圆柱坐标系 为了便于计算，进一步定义桨盘圆柱坐标系，与基于桨盘的直角坐标系的关系如图2所示，二者矢量转换通过(4)式完成。 （4） 考虑到桨叶将会发生变形，还要建立一个桨叶偏转坐标系，它与桨盘圆柱坐标系的关系如图3所示。 图3 桨叶偏转坐标系 图4桨叶剖面气动力和来流示意图 桨叶偏转坐标系与桨盘圆柱坐标系之间的转换式为 （5） 在图3中，由下式确定 （6） 桨盘上动量源项可以通过下式计算 （7） 其中，b为桨叶片数，为桨叶经过动量源项对应的网格所需时间与旋转周期的比例，是坐标系下桨叶微段所受的空气作用力，由通过坐标转换得到。为坐标系下桨叶微段所受的空气作用力，如图4所示，可以由下式计算 （8） 其中 （9） （10） 结合上述关于坐标系变换的内容，可以得出动量源项的计算过程，具体计算步骤为： 第一步：读入单元几何参数及计算域坐标下的来流速度。由坐标判断方位角及到桨毂的距离，结合已知的转速和周期变距参数就可以得到桨叶在这个位置上的桨距、转动牵连速度等。通过坐标变换获得桨叶剖面坐标系下的气流速度。结合转动牵连速度和剖面来流速度可以获得诱导角。 第二步：由第一步获得的参数进一步得出桨叶剖面攻角和来流速度。空气密度、弦长及翼型升、阻力线斜率等。采用公式(8)、(9)及(10)可以得出桨叶在某一单元上所受的升力和阻力。 第三步：通过坐标转换