利用欧拉与NS方程模拟气动弹性过程中精度影响的研究.pdfVIP

第八膈(2003年)全国空气弹性学术交流会论文 利用欧拉和NS方程模拟气动弹性过程中 精度影响的研究+ 叶正寅、田八林 两北下业人学飞机系(西囊，邮编710072) 摘要用中心差分的有限体积法离散Navier-Stokes方榉的通量项，川取时间方法离散 时间导数项，为提高网格效率．采Ⅲ无限插值理论生成每一时间步的O型网格。NS方 祥与篼型运动方程进行耦合求解，在翼型运动方程中，不考虑结构的非线性项，以便重 点研究流动1E线性引起的非线性力学现象。数值计算表明，在气动弹性动态过程的模拟 中，使用烈时间方法时在每一步时间步长中，精度标准需要达到10一7左右。 关键词非线性；气动弹眭；数值模拟； 气动弹性特性预测和分析是飞行器设计中的一项必不可少的上=作，在传统的数值计算 方法(包括进口的商业软件包)中，使用的气动力模型都是以势流、线化方程为理论基础的 菲定常升力面方法，此技术途径的计算效率较高，但其主要缺点有：其空气动力理论是以零 迎角为基础的数值方法，只能适用丁-传统的附着流型，不能反映出迎角的影响．尤其是不能 应IL【』于流动出现分离的情况，而现代E行器已经广泛采用分离流型的气动设计思想：其次， 在很重要的跨音速阶段传统方法却不适用，而且它无法模拟激波随结构变形所产生的运动变 化。另外，传统方法不能考虑实际飞行中有迎角时的静载对弹性结构引起的静变形，而静变 形实际上又政变了机冀各部分的迎角分布，甚至静载引起的静变形会超出弹性结构的线性范 围。 近年来国内外从求解欧拉、Ns方程出发发展了计算气动弹性的计算方法”11，但这种 途径还处在初步发展阶段，很多影响因素需要通过人量的计算研究得出规律性的认识，只有 这样大量的经验积累，新的数值方法才能成功地应用于工程实际型号。本文主要研究双时问 方法的每一时间步长中精度标准对气动弹性计算结果的影响。 1非定常NS方程求解的双时间方法 +本文得到航空科学基金(编号：00A53001)资助 7 1 1网格生成 为了提高计岸效率，本文采_L}』改进的无限插值理论生成翼型的O型网格，其改进的融 合函数莆l网格生成方法具体见文献【4】。 1．2N—s方程的求解 本文所用的流动基本方科为可压缩、非定常N—S方程，其积分形式为 未妒们扩釉=。 其中： 4 Q=【npu，pv，pw，el z方向的速度分量和单位体积的总内能，再为面积分的法 P．“，v，w，e分别为空气密度、x，yt 向单If奇：向量，V为体积分域，ay为包含体积分域的边界，JFl为通量项，它包括无粘项|F：和 枯性项，．两部分： F=％十‘ (2) 戈于舣时间方法求解非定常N—s方程，详见文献【4】。 2结构运动方程的求解 本文认为翼型的变形量仍在线性范围之内，所以采用传统振动模态的概念来描述翼型的 变形： Z(x，Y，f)=∑中．(上，y)·q冉) (3) 其中．N为所选取的振型阶数，中i(互，Y)为机翼的第1阶振型．qj【f)为与第1阶振型对应 的广义坐标。 应用拉格朗日方程．翼型的运动方程可以写为下矩阵形式： [J=}，Jf叮}+[Gl{q}+【K】fq】=fF} (4) 其中，[MIN广义质量矩阵，【G】为结构阻尼项，【KIY,J广义刚度矩阵，{F)为广义气动力。 为了便于用龙格一库塔方法求解方程(4)，引入状态变量{El，即 (E)=[{q}7，㈨7]7=‰，q2，q3，…，qⅣ，ql*q2，…，q～】7 1jlI|，方程(4)变成： cs， t亡，=[：≥；：j品≯：：二；!；：：G，]tE，+[【：r]t，， 其中．【OJ为零方阵，田为单位阵。 8