改进的电热除冰系统的仿真方法

建立了电热除冰系统矩形微段的二维简化物理模型，采用焓法模型建立相变导热问题的微分控制方程，用控制容积法对该模型的控制方程及其边界、初始条件进行离散，使用块修正技术对原有的离散方程组的求解方法进行了改进，得到各个时刻的温度分布及冰的融解情况。计算结果与国外文献结果吻合，验证了改进的算法和程序的正确性。并计算分析了时间步长的选择、网格的划分对计算结果的影响。     

脉冲式燃烧风洞起动特性数值研究

利用计算流体力学软件Fluent,采用隐式双时间步方法模拟非定常过程,控制方程的离散采用二阶迎风格式,湍流模型选择剪切应力输运(SST)k-ω模型,并求解组分扩散方程,对脉冲式燃烧风洞的起动过程进行了数值模拟,分析了喷管起动时流场结构的变化情况和自由射流区超声速流场的建立过程.计算结果为改进风洞起动特性提供了一定的参考依据.      

CFD技术在动态压力校准中的应用研究

采用计算流体力学(CFD)流场可视化的手段,用数值模拟的方法研究了动态压力校准装置.研究工作主要集中于激波管动态压力校准装置和正弦动态压力校准装置的数值模拟,控制方程采用雷诺平均的N-S方程描述,方程的封闭采用精度相对较高的高雷诺数湍流模型,用有限体积法对控制方程进行空间离散和时间离散.并通过实验验证了控制方程和计算方法的正确性.     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

卡尔曼滤波在自旋卫星章动控制中的应用

细长体结构卫星的自旋运动是不稳定的,章动有发散趋势,须进行实时的章动监视,且章动角的测量值包含有误差.通过引入一个离散控制过程的系统,介绍了卡尔曼滤波器的基本原理,建立了自旋卫星自旋运动简化的动力学模型,给出了运动方程和观测方程,详细地推导了适用于自旋卫星章动控制的卡尔曼滤波处理过程,利用该模型对某自旋卫星章动进行仿真计算,并对姿态章动联合控制和应急章动控制情况下的章动控制应用情况进行分析,得出结论：基于卡尔曼滤波的卫星章动角计算的方法正确、结果准确.     

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

低速长航时无人机垂直突风过载分析

低速中空长航时无人机的升力非线性较强,利用标准离散突风模型和连续紊流模型计算得到的突风过载差别较大.针对该问题,本文采用改进的离散突风模型求解突风载荷,即直接求解突风条件下原始的运动学方程,并在运动历程中考虑非线性升力的影响.结果表明:改进的离散突风模型计算得到的突风载荷明显低于标准离散突风模型所得的载荷,更贴近无人机的实际情况,且改进的离散突风模型与连续紊流模型的计算结果也相近.     

耦合高效配平策略的旋翼气动特性分析方法

建立了一个耦合高效配平策略和计算流体力学（CFD）技术的旋翼气动特性分析方法.采用CFD求解器进行雅克比矩阵和旋翼流场的计算,建立了一个包含直接配平法的旋翼气动特性分析平台;引入差量法思想对配平策略进行简化,发展了一种差量配平法,解决了直接配平法中雅克比矩阵计算效率低下的问题,显著缩短了计算时间.旋翼流场求解中控制方程采用Navier-Stokes方程,空间离散格式选取2阶迎风Roe格式,时间推进采用隐式lower-upper symmetric Gauss-Seidel (LU-SGS)格式,湍流模型为Baldwin-Lomax (B-L)模型.分别采用建立的基于直接配平方法和差量配平方法的气动特性分析平台对AH-1G和Helishape 7A直升机的旋翼进行了配平计算和研究,对比分析了两种配平方法的计算效率和精度及背景网格的网格量对配平结果的影响.结果表明:与直接配平法相比,耦合差量配平法的气动特性分析方法计算圈数更少,可节省至少60%的计算时间;使用疏化的背景网格虽然不利于捕捉流场细节,但是可以在兼顾配平结果准确性的同时,降低配平计算代价,提高旋翼配平和流场计算的效率.      

直升机附件舱温度场的数值模拟

为考查某型直升机附件舱的热控制能力,采用有限元软件分析了其在强制冷却和自然冷却情况下的温度场分布.应用非结构化网格和有限体积法进行计算区域和控制方程的离散,采用标准k-ε模型封闭湍流控制方程.针对3种冷却系统设计方案,分析了不同冷却气流进口尺寸、相对位置以及不同附件散热量对附件舱温度分布的影响,计算了自然冷却情况下附件舱密闭时的温度场.计算和试验结果均表明,所采用的冷却方案可以满足附件舱冷却的需要.      

Structural Nonlinear Flutter Characteristics Analysis for an Actuator-fin System with Dynamic Stiffness

The flutter characteristics of an actuator-fin system are investigated with structural nonlinearity and dynamic stiffness of the electric motor. The component mode substitution method is used to establish the nonlinear governing equations in time domain and frequency domain based on the fundamental dynamic equations of the electric motor and decelerator. The existing describing function method and a proposed iterative method are used to obtain the flutter characteristics containing preload freeplay nonlinea...     

粘弹壁板颤振的非线性动力特性

研究粘弹材料壁板在超声速气流作用下颤振时的分岔及混沌等复杂动力学特性。采用von Karman大变形理论及Kelvin粘弹阻尼模型建立壁板的动力学方程,通过线性活塞理论建立气动力模型。利用迦辽金法将壁板颤振模型转化为常微分方程组,并使用Gear的BDF方法进行数值求解。通过数值模拟研究了该系统在粘弹阻尼作用下的动力学行为以及粘弹阻尼的影响。计算结果表明,粘弹壁板颤振系统表现出丰富的动力学行为,其二次分岔特性很复杂。随着粘弹性阻尼的增大,系统的稳定解区域在减小,而静态屈曲解几乎不受影响,同时发现混沌运动区域也随着粘弹阻尼的增大而减小。     

等温不可压流体系统稳态数值计算方法研究

航空工业领域广泛存在着流体系统,任何流体系统均可化为只有节点和分支组成的通用流体网络,本文研究了这种通用流体网络的数值计算方法.将连续方程施加到节点,将动量方程施加到分支,建立起等温不可压流体系统的非线性控制方程组.根据Newton-Raphson迭代思想,研究了基于高斯-约当消去法的流体系统控制方程组数值求解方法.在此基础上,采用压缩矩阵存储技术开发了流体系统仿真平台FOCUSS-FS.通过算例进行验证,结果表明:本文计算方法能够准确获得任意流体网络的压强和流量分布,收敛速度快且效率较高,内存消耗小.     

气流激励下叶片动力响应分析方法

利用在物理过程上的耦合代替数学方程上的耦合,完善了双向顺序耦合求解的理论,进一步给出了双向顺序耦合的数值求解方法,并描述了其在物理上对应的过渡过程.结合具体工程算例,分析得到某小型发动机二级静叶在气流激励下的流场分布形式和叶片动力响应的特点,指出了双向顺序耦合数值求解方法的优越性.      

多点敷设支撑层的粘弹性约束阻尼梁振动分析

针对夹层约束阻尼梁结构的减振特性进行深入研究,在其基础上扩展多处敷设增设支撑层的约束阻尼,建立了多点铺设增设支撑层约束阻尼梁的振动特性模型,利用粘弹性材料的Maxwell模型、假设模态法、有限元理论和力学原理,结合功能特性方程和Lagrange方程,导出了增设支撑层约束阻尼梁的振动位移响应方程.通过对敷设双支撑约束阻尼简支梁在支撑层厚度变化的参数研究,可以看到增设适当厚度的支撑层能够达到更好的减振效果.     

导弹垂直发射系统三维超声速射流的数值模拟

对导弹垂直发射系统(VLS,Vertical launching system)内的三维湍流流场进行了数值模拟.采用了单元中心有限体积法和三阶迎风TVD差分格式离散Navier-Stokes方程,采用双时间步和预处理过程提高计算收敛速度,引入对角化的双时间ADI算法求解控制方程.仿真结果准确地模拟了发射系统内的激波形成过程,真实地反映出VLS中的复杂流动.研究表明,压力室底部射流中心的波后压强以大约1 200Hz的频率振荡.      

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...