基于FW-H方程的旋翼气动声学计算研究

由流体力学N S方程导出的非齐次波动方程———FfowcsWilliams Hawkings方程(简称FW H方程),可以精确地描述在静止流体中运动的物体与流体相互作用的发声问题。以FW H方程为理论模型,将旋翼桨叶运动发声问题等效为包含桨叶的任意运动控制面(声源面)的声辐射问题,并在旋翼绕流Euler方程数值模拟的基础上,在时域内计算了悬停旋翼和前飞旋翼的声场。应用于UH 1H和AH 1/OLS两种旋翼模型的气动声学计算表明:计算结果与噪声实验值符合良好;所研制的程序不仅能够较准确地计算单极子噪声和偶极子噪声,而且具有较强的跨音速四极子噪声预测能力。     

关于用影响系数法处理准静弹性飞机动态特性问题的一些意见

 <正> 一、引言近几年来,国外刊物发表了一些用影响系数法分析和计算弹性飞机动态特性问题的文章（见文献[1～7]）。这些文章在准静弹性假设下处理飞机运动问题时,引入了一组惯性导数和方位角导数,从而使运动方程在形式上与习惯用的刚性飞机运动方程有了差别。这样,除增加了导数数目外,还使过去对刚性飞机建立起来的、关于某些导数的概念以及导数与运     

昆虫悬停飞行的动稳定性:理论分析与数值模拟

首先,从昆虫(身体和拍动翅)的动力学方程和N-S方程出发,在一定假设下,将运动方程简化为6自由度刚体的方程,并用线化理论给出了反映昆虫悬停飞行纵向动稳定性的理论解.然后,用完全的动力学方程和N-S方程偶合的数值解对上述理论的简化假设进行了检验.     

惯性平台离心试验的力学分析

惯性平台离心试验是平台例行试验项目之一。本文叙述:(1)惯性平台组合的运动分解;(2)惯性力的分布;(3)惯性力对外框架轴之矩Mx的表达式;(4)Mx与离心机臂长成反比;(5)达朗贝尔原理的应用。     

跨音速机翼结构动力响应的计算研究

通过交替求解流体运动方程和结构动力方程来计算时域内二维机翼的结构动力响应,并研究了结构参数对系统稳定性的影响,使用控制中的相平面图对计算结果进行分析,来判断系统的结构稳定性。计算中流体运动方程采用二维非定常Ｅｕｌｅｒ方程,并用Ｊａｍｅｓｏｎ提出的有限体积方案,四步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进求解;结构模型采用二自由度ＰＡＰ模型。经分析,计算结果合理。     

旋流加力室火焰前锋的计算

基本方程 1.微分方程 计算旋流加力室火焰前锋的位置,基于气泡理论。由于旋转力场,冷气流离心向外运动,迫使径向点燃的火焰微团向心运动。火焰微团的位置由轴向和径向速度来确定。假设这两个运动是独立的。      

弱可压缩流体与可压缩流体模型的管路水击研究

为提升液体推进系统模型对可压缩流体的仿真能力,针对供应管路的瞬变流动,采用一维有限体积法,使用Modelica语言基于MWorks平台分别开发了适用于补燃循环发动机仿真模型库的弱可压缩流体与可压缩流体的管路模型,不同于控制方程中显含声速的弱可压缩流体模型,可压缩流体模型通过将守恒方程与真实流体的SRK型状态方程耦合求解来计算压力波的传播。在管路供应系统中,分别以液氢、液氧及液甲烷为工质进行关阀水击的仿真研究,结果表明:使用SRK状态方程避免了可压缩流体模型对内能拟合公式的依赖,减轻了建模难度;可压缩流体模型对液氧及液甲烷的仿真结果较为准确,最大误差不超过2%,对液氢的最大计算误差约9%;相对于可压缩流体模型,弱可压缩流体模型计算得到的水击压力在振幅及频率上均偏大;可压缩流体模型计算得到的流体密度波动不超过3%,因此工程实践中处理水击问题时将液体假设为弱可压具有一定的合理性。     

非惯性非正交曲线坐标系中的边界层方程

应用张量分析方法将 N-S方程写于作非惯性运动的非正交曲线坐标系中 ,然后用匹配渐近展开法导出此坐标系中的边界层方程。消去压强项后 ,与惯性坐标系中的边界层方程相比 ,非惯性系中的边界层方程在 x1和 x2 方向的动量方程中分别只多了一哥氏力项 ,它们是由物体转动角速度的物面法向分量产生的 ,其它惯性力项均不进入边界层方程。方程中所含曲线坐标系的拉梅系数及其导数只需在物面上取值 ,物面曲率影响不进入边界层方程     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

翼型跨音速小扰动势流-边界层计算

翼型绕流的外部无粘流,采用跨音速小扰动速势方程和改进小扰动速势方程对三种差分格式进行对比计算:(1)Murman-Cole非守恒式;(2)守恒式;(3)Engquist-Osher式。在相同的初值条件下,它们的计算稳定性、收敛性和收敛解基本相同,但激波位置(1)在前,(2)在后,(3)在中间。对人工时间阻尼项εφ_(xt)的作用进行了线化分析,所得结论与数值实验结果相符。二级精度格式计算表明:采用这种方案能在几乎不增加计算机时的前提下,大大提高差分计算精度。本文最后采用边界层积分方程法,计算翼型边界层,并与势流进行迭代计算。在弱激波条件下,迭代收敛解与实验值接近。     

发动机故障诊断经验影响系数的定量化问题

本文提出了根据发动机故障样本建立故障方程的方法。发动机经验故障方程是不同于发动机小偏差方程的另一类故障方程,它为发动机故障诊断提供了一条简便易行的途径,并且使发动机故障诊断范围扩展到气路分析方法难以适应的场合。文中讨论了经验故障方程的建立方法及其在发动机故障诊断中的应用,给出了利用经验故障进行发动机故障诊断的实例。     

发动机经验故障方程的建立与应用

本文提出了根据发动机故障样本建立故障方程的方法。发动机经验故障方程是不同于发动机小偏差方程的另一类故障方程,它为发动机故障诊断提供了一条简便易行的途径,并且使发动机故障诊断范围扩展到气路分析方法难以适用的场合。文中讨论了经验故障方程的建立方法及其在发动机故障诊断中的应用,给出了利用经验故障方程进行发动机故障诊断的实例。      

N-S方程及其简化形式在圆柱绕流分离点邻域的数学性质

证明了在圆柱绕流分离点的邻域上 ,N-S方程和高智的简化 N-S方程是正则的 ,即存在级数解 ;而 Davis和 ГКП的简化 N-S与经典的边界层方程一样存在 Goldstein奇性     

旋转坐标系中可压缩湍流的基本方程和湍流模型

本文在可压缩湍流瞬时量基本方程的基础上,用质量加权平均的概念,推导了旋转坐标系中可压湍流流动的连续方程、动量方程和能量方程。进而导出了旋转坐标系中可压湍流的应力输运方程、脉动平均动能方程和湍流耗散率输运方程。这些方程中包括了由旋转引起的 Coriolis力和离心力的作用项。这些方程是在旋转坐标系中研究平均 N-S方程求解湍流流场和研究湍流模型的基础。      

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

四边简支压电热弹性体层合开口壳的精确解

根据压电热弹性体的本构关系和热传导关系的相似性,通过增维方式将热传导关系写进压电热弹性体的本构关系中,直接得到了一个新的增维的压电热弹性体本构关系。依据此增维的本构关系,结合压电热弹性体平衡方程和热平衡方程,成功地导出了压电热弹性体层合开口壳机、电、热耦合问题的齐次状态方程。这类方法大大简化了人们在分析压电热弹性体耦合问题时通常要求解的非齐次微分方程和关于温度的平衡方程和导热方程的二阶微分方程的繁琐方法,同时,齐次状态方程的提出大大减少了数值计算的工作量。     

Motion equations of hemispherical resonator and analysis of frequency split caused by slight mass non-uniformity

The mass non-uniformity of hemispherical resonator is one of reasons for frequency split, and frequency split can cause gyroscope to drift. Therefore, it is of great significance to analyze the relationship between mass non-uniformity and frequency split, which can provide a theoretical basis for mass balance of imperfect resonator. The starting point of error mechanism analysis for gyroscope is the motion equations of resonator. Firstly, based on the Kirchhoff-Love hypothesis in the elastic thin shell theory, the geometric deformation equations of resonator are deduced. Secondly, the deformation energy equation of resonator is derived according to the vibration mode and relationship between the stress and strain of hemispherical thin shell. Thirdly, the kinetic energy equation of resonator is deduced by the Coriolis theorem. Finally, the motion equations of resonator are established by the Lagrange mechanics principle. The theoretical values of precession factor and natural frequency are calculated by the motion equations, which are substantially consistent with the ones by the finite element method and practical measurement, the errors are within a reasonable range. Simultaneously, the varying trend of natural frequency with respect to the geometrical and physical parameters of resonator by the motion equations is consistent with that by the finite element analysis. The above conclusions prove the correctness and rationality of motion equations. Similarly, the motion equations of resonator with mass non-uniformity are established by the same modeling method in case of ignoring the input angular rate and damping, and the state equations with respect to the velocity and displacement of vibration system are derived, then two natural frequencies are solved by the characteristic equation. It is concluded that one of reasons for frequency split is the 4th harmonic of mass non-uniformity, and thus much attention should be paid to minimizing the 4th harmonic of mass non-uniformity in the course of mass balancing for imperfect resonator.     

锥形薄壳的弹性稳定性

 本文用变分原理推导了锥形薄壳稳定性的基本微分方程式。在推导过程中,考虑了壳体微元在中面内的变形和绕各个轴的转动。所得的方程最后化为用三个位移表达的微分方程式,它适用于各种载荷情况,并避免了“唐纳尔型”方程的限制和出于物理想像来建立锥形壳微分方程式的错觉。作为本文理论应用的一个例子,最后用伽辽金方法近似解答了全锥在外压作用下的稳定性问题,获得了简洁的临界压力计算公式。所得结果与实验一致。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。