基于CFD_CSD耦合算法的机翼颤振分析_曾宪昂_图文

1、基于CFD /CSD 耦合算法的机翼颤振分析¹曾宪昂,徐敏,安效民,陈士橹(西北工业大学航天学院,陕西西安710072摘要:用计算流体力学/计算结构力学(CFD/CSD耦合算法对标准气动弹性模型AGARD 445.6机翼作了颤振分析,主要研究机翼的跨音速颤振求解问题。采用常体积转换法(CVT 进行流体与结构之间的数据交换并运用松耦合方法对气动弹性方程进行时域推进仿真。计算机翼在M a =0.4991.072的颤振边界,并将计算结果同偶极子格网法(DLM 的计算结果与试验结果比较,结果显示CFD/CSD 耦合计算结果较DLM 计算结果更接近于试验值,尤其是在非线性强的跨音速区域。可见,

2、CFD /CSD 耦合计算比DLM 具有很大的优越性。关键词:计算流体力学/计算结构力学(CFD/CSD耦合算法,AGARD 445.6机翼,颤振分析,常体积转换(CVT ,偶极子格网法(DLM 中图分类号:V 211.47文献标识码:A 文章编号:1000-2758(200801-0079-04颤振是飞行器飞行时常见的一种气动弹性现象,它对飞行器的破坏是灾难性的。在过去,广泛运用于飞行器的颤振计算方法是偶极子格网法(DLM,它是基于线化位流理论的一种颤振计算方法。在很多的商业软件中都应用DLM 进行气弹分析,如广泛应用于工业设计的M SC .NASTRAN 的气动弹性模块1等。但是DLM 是

3、一种基于线化理论的方法,无法解决非线性强的流场,而且它大多采用的是平面模型,无法计入机翼的厚度、迎角等。近年来,CFD 计算技术发展十分迅速,计算机的性能也有很大提高,因此CFD/CSD 耦合计算方法2,3也迅速发展起来。由于CFD 计算求解的是非线性方程组,因此这种方法可以计算飞行器在非线性强的流场中的运动,同时这种方法可更加直观、实时地显示飞行器结构的变形及流场的变化,并且它还可以计算较为复杂的外形。本文首先利用有限元方法对机翼结构进行模态分析;然后采用Euler 方程计算非定常气动力;结构与流体的数据交换采用常体积转换法(CVT 进行4。在具体计算中,给定一个来流速度,在时域中推进气动弹

4、性方程,观察各阶模态位移的时间历程。若幅值增长则来流速度大于颤振速度,若幅值衰减则来流速度小于颤振速度,当幅值保持不变时,此时的来流速度即为颤振速度。1机翼结构模型AGARD 445.6机翼是国际上用于检验颤振计算方法的一个标准模型,它有着较为完备的风洞试验数据5 所示。图1AG A RD 445.6机翼前四阶模态变形图2008年2月第26卷第1期西北工业大学学报Jour nal o f No rt hw ester n P olyt echnical U niv ersity F eb.2008V ol.26N o.1¹收稿日期:2007-01-17基金项目:国家自然科学基金(90

5、405002资助作者简介:曾宪昂(1983-,西北工业大学硕士生,主要从事气动伺服弹性的研究。表1AGARD445.6机翼的模态频率与试验值的比较一阶模态9.60239.6150.62694.9572计算流体动力学及动网格方法在惯性坐标系下,非定常积分形式的Euler 方程可写为55t 8U d 8+SF d S =0(1式中,U 为守恒型的流动变量;F 为无黏的流通矢量;8为运动控制体积;S 是运动控制体的表面积。采用非结构网格,全场网格数共253563个(机翼表面网格见图2,控制方程采用有限体积法离散求解。计算中动网格采用弹性系数算法对网格进行重构。弹性系数算法是将任意2个网格节点之间的边

6、等效为1根弹簧。结点的位移将产生与位移成正比的弹性力。在平衡状态下,每个结点所受的所有与它相连边上的弹性力之和为零。通过平衡方程可得结点位移的计算方程$S m +1i=njkij smj/njk ij(2式中,s i 是结点i 的位移;s j 是与结点i 相连的结点j 的位移;n 是与结点i 相连的结点数目;k ij 是结点i 与结点j 之间的弹性系数,其计算式为k ij =ûs i -s j û-1/2( 3那么已知边界结点的位移就可通过(2式进行迭代,收敛后就得到了新的内部结点位置。图2机翼表面网格3CFD /CSD 数据交换方法在CFD/CSD 耦合计算中,流体与结构

7、之间的数据交换方法十分重要。它直接关系到计算结果的好坏。在数据交换过程为了保证力守恒、力矩守恒以及虚功守恒。本文采用CV T 方法进行数据交换:对每个气动网格点a (t ,首先要在结构网格上找出离它最近的3个结构网格点,其顶点用s i (t 、s j (t 、s k (t 表示,随时间的推移,它们有如下关系式a (t =A s i (t +B s j (t +C s k (t +v (t (s j (t -s i (t ×(s k (t -s i (t (4式中,A 、B 、C 是常数,且满足A +B +C =1,再利用变形时保持有向体积不变以及气动网格点在结构三角形中的投影位置保持

8、不变关系可求得变形后的气动网格点位置。与CVT 方法配套的气动力转换关系式为Q i =A Q a ,Q j =B Q a ,Q k =C Q a 。其中Q a 为气动网格点的气动力,Q i 、Q j 、Q k 为Q a 转换到3个结构网格点的气动力。图3显示了气动网格点对机翼前四阶振型插值得到的网格变形。图3CV T 插值得出的机翼前四阶模态变形图4耦合算法基于模态的气动弹性方程有如下形式q ¨i +2Ni X i q õi +X 2i q i =F i (5式中,q i 为广义模态位移;Ni 为模态阻尼;X i 为模态频率;F i 为广义模态气动力;下标i 表示模态的阶数

9、。本文采用松耦合方法6在时域内对气动弹性方程进行推进。具体步骤为:¹计算t 时刻气动网格点80西北工业大学学报第26卷上的气动力。º将气动网格点上的气动力通过数据交换转换为广义气动力。»用龙格-库塔方法求解气动弹性方程求得t +$t 时刻的广义位移。¼将广义位移转换成实际位移,并通过插值得到气动网格的变形。重复步骤¹¼步就能得到机翼的响应历程。5计算结果计算状态同参考文献5,给定第一阶模态一个小的位移作为初始扰动,为保证计算精度,取时间步长为0.0002s 。计算马赫数范围为0.4991.072。颤振速度和颤振频率都以无量纲的形式给出

10、V f =U /(b s X A L (6X f =X /X A(7式中,U 为来流速度;b s 为机翼半根弦长;X 为振动频率;X A 为机翼固有一阶扭转频率;L 为质量比,L =m /(Q V ,m 为机翼质量,Q 为来流密度,V 为机翼体积。计算结果如图4图8所示。CFD /CSD 耦合计算、DLM 方法计算得到的颤振边界以及试验结果的比较如图9图10 所示。图4M a =0.499,A =0°,V f =0. 4467图5M a =0.678,A =0°,V f =0.4156图6M a =0.901,A =0°,V f =0.3608图7M a =0.9

11、6,A =0°,V f =0.307 8图8M a =1.072,A =0°,V f =0.3534图9A =0°时的颤振速度边界图10A =0°时的颤振频率边界6结论从图9、图10可以很直观地看出CFD /CSD 耦合算法得到的颤振边界比DLM 方法计算结果更接近于试验值,可见CFD/CSD 耦合算法的精度比DLM 高。由于机翼厚度很小,在亚音速区域DLM计算结果是比较精确的,完全适用于工程计算。但在跨音速区(Ma =0.96计算结果与试验值有较大出入,这是由于DLM 采用的线化理论不适合于非线性强的跨音速流场计算,因此它所计算出的颤振速度是不可靠的。

12、CFD/CSD 耦合计算方法可以求解非线性的流场,因此在跨音速区域的计算结果比DLM 方法要精确。在超音速区(Ma =1.072CFD/81第1期曾宪昂等:基于CFD /CSD 耦合算法的机翼颤振分析CSD 耦合计算方法与试验结果还存在相当差距,其主要原因是本文未考虑流体粘性,可以肯定粘性对计算结果有一定影响,但还有其它一些原因。CFD/CSD 耦合计算的精度高,但是计算时间长,尤其是在搜索颤振速度时要耗费大量时间。在后续工作中作者将进一步对基于N -S 方程的CFD /CSD 耦合计算方法以及如何提高CFD/CSD 耦合计算效率方面做深入研究,使其具有更高的实际工程应用价值。参考文献:1R

13、odden W P,Jo hnson E H.M SC/NA ST RA N A er oelastic A na ly sis U ser s Guide Ver sio n 68.T he M acN eal SchwendlerCo rp,19942L ee -R ausch E M ,Batina J T .Calcula tio n of A GA RD W ing 445.6Flut ter U sing N avier -Stokes A er odynamics .AI AA 93-34763R amji Kam ako ti,Wei Shyy ,Siddhar th T ha

14、kur ,Bhavani Sankar.T ime Dependent R A NS Com putation for an A er oelas-t ic Wing.A IA A -2004-08864徐敏,陈士橹.CF D /CSD 耦合计算研究.应用力学学报,2004,21(2:3336X u M ,Chen S L .CF D /CSD Coupling Calculation R esear ch .Chinese Jo urnal o f A pplied M echanics ,2004,21(2:3336(in Chinese X u M ,A n X M ,Chen S L.

15、CF D/CSD Coupling N umerical Co mputat ional M et ho do lo gy.Acta A ero nautica et A str onautica Sinica ,2006,27(1:3337(in ChineseWing Flutter Analysis Using CFD /CSD AlgorithmZeng Xian ang ,Xu M in ,A n Xiaomin ,Chen Shilu(College o f Astr onautics,No rt hw ester n P olyt echnical U niv ersity ,X i an 710072,ChinaAbstract :The flutter analysis of the AGARD wing 445.6,w hich is a standar d a