旁路式双喉道喷管气动矢量特性数值研究

旁路式双喉道喷管(BDTN)通过设置旁路通道即可实现稳定的推力矢量,无需增加额外的高压气源和次流系统,研究其推力矢量特性对于提高飞行器的机动性和敏捷性具有重要意义。采用数值方法对BDTN在不同落压比下的内外流场及其在高温高压环境下的内流特性进行数值仿真。结果表明:BDTN的流场结构与普通双喉道喷管(DTN)相仿,但矢量效果更好,可在推力系数为0.966的同时获得24.6°的推力矢量角;落压比(NPR)大于4时,随着入口总压的增加,喷管推力矢量角和推力系数均逐渐减小;在入口总温为3 000K、落压比为100时,喷管的推力矢量角减小至12.75°,推力系数减小至0.882,即高温高压的内流环境对喷管矢量性能具有很大的消极影响。     

喷管气动参数对推力矢量影响的数值模拟

以二元收扩喷管为对象,开展了基于二次流喷射的流体推力矢量技术研究。基于CFD技术,分析了激波矢量控制技术实现推力矢量的机理。重点分析了喷管落压比NPR,二次流总压比SPR,自由流马赫数Ma∞对流体推力矢量性能的影响。数值模拟结果表明:推力矢量的大小与斜激波的位置、角度以及流动分离区的大小有关。在所计算的参数范围内,推力矢量随喷管主次流的变化规律是:推力矢量角随NPR的增大逐渐减小;随SPR的增加,推力矢量角单调增加;在大落压比时,自由流马赫数对推力矢量的影响是有限的,而在低落压比时,自由流马赫数增加,推力矢量角减小。     

双喉道推力矢量喷管的内流特性研究

为了研究双喉道推力矢量喷管(DTN)在非推力矢量和推力矢量情况下的内流特性,基于数值模拟的方法,计算分析了不同几何参数和气动参数对DTN的影响。结果表明,DTN在非推力矢量时,仅在落压比(ZNPR)为3～4之间才具有较高的内流性能(推力系数达0.97,流量系数为0.94),当落压比增加时,推力系数迅速下降。在推力矢量时,DTN可以获得很大的推力矢量效率(当落压比为4,引射量为3%时达到4),且推力系数也较高(0.94以上),其综合性能优于单喉道偏移和激波操纵式矢量喷管。二次流量、落压比、凹腔扩张角和收敛角、引射角度都对推力矢量状态下的DTN内流性能有着不同的影响。为了实现DTN在推力矢量和非推力矢量下都有较好的内流综合性能,所建议的设计参数为:落压比为3～4,引射量为3%,凹腔扩张角为10°左右,收敛角在20°～30°,引射角度为30°逆流引射角(β=30°)。     

二次流喷射位置对流体推力矢量喷管气动性能影响的数值模拟

数值模拟了二次流喷射位置对基于激波控制的二维收-扩(2DCD)喷管的流体推力矢量气动性能,结果表明喷射位置对矢量角有较大影响,在喷管落压比NPR=4.5时,二次流喷射位置相对位于喷管发散段中部时,矢量角最大;在NPR≥7.2时,二次流喷射位置越靠近出口,矢量角越大.喷射位置对推力系数及流量系数的影响不大.     

基于激波控制的S 弯二元矢量喷管数值模拟

为了掌握S弯二元矢量喷管的气动性能，采用CFD数值模拟方法研究了有无二次流喷射状态下S弯二元矢量喷管激波 诱导的工作机理，以及二次流喷射位置、主流落压比和二次流与主流总压比对S弯二元矢量喷管推力系数、矢量角、壁面静压的影 响。结果表明：在二次流流通面积不变、次流与主流流量比Ws/Wp≤6%的情况下，喷管上、下壁面分别喷射二次流产生的最大矢量 角分别为22.9°和15.9°；喷射位置对矢量角有较大影响，对推力系数影响不大，随着二次流喷射位置逐渐靠近出口，矢量角先增大 后减小；射流位置固定，随着主流落压比的增大，推力系数增大，当主流落压比从2增大到6时，推力系数最多提高17.9 %，矢量角 先增大后减小；随着二次流与主流总压比的增大，推力系数整体呈单调减小趋势，矢量角先增大后减小。     

扩张段射流对旁路式双喉道喷管矢量特性的影响研究

为探究腔体扩张段射流对旁路式双喉道喷管气动矢量特性的影响,采用数值模拟方法对喷管在不同次流入射位置和次流压比下的内流情况进行仿真研究.结果表明:在扩张段引入次流能够改善喷管内流性能,随着次流入射位置后移,推力矢量角先增大后减小,推力系数逐渐增大且增幅渐缓;随着次流压比增加,喷管推力矢量角逐渐增加后基本保持不变,推力系数先增加后快速下降,而矢量效率先急剧上升后趋于平稳;改进后的喷管在最佳算例中得到推力矢量角为27.59°,推力系数为0.956,矢量效率3°/1％次流流量.     

次流通道对双喉道气动矢量喷管的性能影响研究

为了进一步了解次流通道对双喉道气动矢量喷管(Dual Throat Nozzle,DTN)气动矢量性能的影响,基于二维DTN喷管构型进行了详细的数值模拟研究。结果表明:收扩的次流通道相对于平直的次流通道不会改善DTN喷管的推力矢量性能;随着次流流量增加,DTN的推力矢量角存在一个最大值,超过该值各推力矢量性能均会降低;次流通道宽度对DTN喷管推力矢量性能影响显著:(1)次流宽度不影响所能达到的最大推力矢量角;(2)小的次流宽度在进口总压较高的条件下具有更高的推力矢量角;(3)增大次流宽度能降低次流总压,但达到的矢量角会有所减小。     

固定几何气动矢量喷管气动性能数值仿真

为掌握固定几何气动矢量喷管气动性能,通过CFD数值模拟的方法,研究了主流落压比、扩张段二次流落压比、扩张段二次流角度和引射对固定几何气动矢量喷管轴向推力系数的影响;主流落压比、扩张段二次流落压比和扩张段二次流角度对矢量角的影响;主流落压比、喉道二次流落压比和喉道二次流角度对喉道控制率的影响。结果表明:随主流落压比增大轴向推力系数增大,矢量角减小,喉道控制率减小;随扩张段二次流落压比增大推力系数减小,矢量角增大;随喉道二次流落压比增大,喉道控制率增大;随扩张段二次流角增大轴向推力系数减小,矢量角略有减小;随喉道二次流角增大喉道控制率增大;随引射方式增加喷管推力系数增大。     

轴对称射流矢量喷管非矢量状态下推力特性研究

通过数值模拟,研究了非矢量状态下（射流流流量为零）轴对称射流矢量喷管在低落压比时（NPR≤3.5）,射流缝内形成的空腔流动对喷管推力特性的影响;利用这种影响来控制低落压比时喷管的气流分离位置,提高了喷管处于低落压比时的喷管推力系数;推力系数可以提高到0.90以上。因此合理布置双连通射流缝的位置,不仅有利于提高射流的推力矢量效率,而且有利于提高低落压比无次流时的推力系数,从而弥补射流矢量喷管几何固定带来的不足。     

几何构型对逆流推力矢量喷管特性影响

采用数值模拟方法研究了几何构型对逆流推力矢量喷管内部流动结构及性能的影响.研究结果表明:主流未附体时,增大缝宽可提高推力矢量角,但同时却降低了合成推力系数,抽吸二次流流向由同向转变为逆向且流量逐渐增大;增大斜切角可致使推力矢量角减小,合成推力系数先增大后减小.随着斜切角的增大,抽吸二次流流向由逆向逐渐转变为同向二次流,流量比在2%以内,且呈先减小后增大趋势.     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

μ理论在电液负载模拟器中的应用

 研究μ理论在电液负载模拟器中的应用，提出电液负载模拟器的鲁棒控制策略。电液负载模拟器是一个复杂的机-电-液复合系统，且是一个强耦合、时变受控对象。综合考虑参数变化、模型变动和干扰等不确定性的影响，利用μ综合控制理论设计电液负载模拟器的鲁棒力控制器，并通过μ分析使用鲁棒力控制器时系统的鲁棒稳定性和鲁棒性能。给出使用鲁棒力控制器和经典力控制器时的实验结果，实验结果表明所设计鲁棒力控制器的有效性和优越性。     

基于圆型限制性三体问题模型的行星卫星逃逸能量研究

 基于由飞行器、行星及其卫星组成圆型限制性三体问题模型，通过庞加莱映射的方法，研究了飞行器从行星卫星附近逃逸的问题。在Jacobi常数确定的前提下，通过逆向积分，飞行器从L1或L2点附近返回近月点，得到近月点速度出发速度。研究结果表明绕飞L1点和L2点逃逸行星卫星需要的最低能量是不同的，从月球表面逃逸所需的速度脉冲分别比开普勒算法节省46.5m/s和42.3m/s，且均小于Villac等人在Hill模型下得到38.9m/s，从而改进了Villac等人的相关工作，同时也给出了从太阳系主要行星卫星表面逃逸所需的最小能量。     

Model Development and Adaptive Imbalance Vibration Control of Magnetic Suspended System

对星载磁悬浮单框架控制力矩陀螺,建立其转子的平动与转动的动力学方程。所建模型反映了转子的动、静不平衡特性,以及转子运动与框架、星体运动的耦合关系。模型是主动磁悬浮控制与动框架效应分析的基础。对动、静不平衡量未知的转子,设计了自适应对中控制器。在分散PD、比例交叉反馈控制器与高低通滤波器的基础上,通过动、静不平衡量的辨识与补偿控制,使转子绕其惯量主轴旋转,衰减角动量交换执行机构高速转子因不平衡而对基座产生的干扰力,消除星上主要颤振源,达到降低卫星姿态抖动的目的。     

Love Waves in Layered Graded Composite Structures with Imperfectly Bonded Interface

Theoretical analysis and numerical calculations of Love wave propagation in layered graded composites with imperfectly bonded interface are described in this paper. On the basis of WKB method, the approximate analytic solutions for Love waves are obtained. By the interface shear spring model, the dispersion relations for Love waves in layered graded composite structures with rigid, slip, and imperfectly bonded interfaces are given, and the effects of the interface conditions on the phase velocities of Love waves in SiC/Al lay- ered graded composites are discussed. Numerical analysis shows that the phase velocity decreases when the defined flexibility parameter is greater. For the general imperfectly bonded interface, the phase velocity changes in the range of the velocities for the rigid and slip interface conditions.     

Compression Creep Behavior of High Volume Fraction of SiC Particles Reinforced Al Composite Fabricated by Pressureless Infiltration

The compression creep deformation of the high volume fraction of SiC particles reinforced Al-Mg-Si composite fabricated by pressure-less infiltration was investigated. The experimental results show that the creep stress exponents are very high at temperatures of 673 K, 723 K and 773 K, and if taking the threshold stress into account, the true stress exponent of minimum creep strain rate is still approximately 5, although the volume fraction of reinforcements is very high. The creep strain rate in the high volume fraction rein- forced aluminum alloy matrix composites is controlled by matrix lattice diffusion. It is found that the creep-strengthening effect of high volume fraction of silicon carbide particles is significant, although the particles do not form effective obstacles to dislocation motion.     

Experimental Investigation on Performance of Pulse Detonation Rocket Engine Model

The PDRE test model used in these experiments utilized kerosene as the fuel, oxygen as oxidizer, and nitrogen as purge gas. The solenoid valves were employed to control intermittent supplies of kerosene, oxygen and purge gas. PDRE test model was 50 mm in inner diameter by 1.2 m long. The DDT （deflagration to detonation transition） enhancement device Shchelkin spiral was used in the test model. The effects of detonation frequency on its time-averaged thrust and specific impulse were experimentally investigated. The obtained results showes that the time-averaged thrust of PDRE test model was approximately proportional to the detonation frequency. For the detonation frequency 20 Hz, the time-averaged thrust was around 107 N, and the specific impulse was around 125 s. The nozzle experiments were conducted using PDRE test model with three traditional nozzles. The experimental results obtained demonstrated that all of those nozzles could augment the thrust and specific impulse. Among those three nozzles, the convergent nozzle had the largest increased augmentation, which was approximately 18%, under the specific condition of the experiment.     

高超声速飞行器高温流场数值模拟面临的问题

随着高超声速飞行器目标光辐射和电磁散射特性研究的发展和深入,高温流场特性日益引起人们的关注。由于高温流场特性研究中涉及到非常多的复杂气动现象,如气动加热、烧蚀、辐射、燃烧、化学反应以及湍流等,因此其数值模拟面临着诸多挑战。这里基于连续流计算流体力学(CFD)技术和稀薄气体蒙特卡罗直接仿真(DSMC)方法,从化学物理模型建模、方法稳定性与数值求解效率出发,分析了高超声速飞行器外部绕流、尾迹和发动机喷焰三方面的流场特性数值模拟在不同弹道、热防护手段和飞行流域环境下所面临的问题。在此基础上提出了数值求解技术和化学物理模型建模今后需要发展的方向,为有效提高高超声速高温流场特性数值模拟效率、增加流场特性预测精度提供了指导,从而为研究流场对高超声速飞行器目标光辐射和电磁散射特性影响提供有效的基础数据。     

Adaptive Control of Helicopter Ground Resonance with Magnetorheological Damper

The methodology for adaptive control of helicopter ground resonance with magnetorheological （MR） damper is presented. The adaptive inverse control method is used to control the output damping force of MR damper and the range of the damping force is given. Through the adaptive inverse control, the damping force of MR damper is fit to a desired damping force. With the background of applying MR damper to control of helicopter ground resonance, a model of loss force and an adaptive arithmetic for stabilization of the coupled rotor/fuselage system are presented. The simulation shows that the controller presented in this paper can stabilize the rotor/fuselage coupling system quickly and control the helicopter ground resonance effectively.     

The Full Flowpath Analysis of a Hypersonic Vehicle

对一种类X43-A吸气式高超飞行器全流道开展了M7一级的三维数值仿真研究,深入探索了进气道处于起动状态和不起动状态时全流道的冷流流场结构和和气动力特性,且部分结果与实验数据进行了对比。研究结果表明:(1)前体横截面上存在显著的展向压强梯度,使得经过预压缩的气流偏离了进气道进口,但同时也减少了进入内通道的边界层气流,提高了进口流场的品质;(2)后体喷流股的膨胀过程受到了周围外流的显著干扰,因而沿流动方向其截面形状不断发生变化,如在喷口附近为近似矩形,而在后体末端附近则演化为近似三角形;(3)当进气道处于不起动状态时,其全流道流动结构发生了显著变化,进气道的外部压缩波系往复振荡,尾喷管出口的喷流股也在不停的膨胀和收缩,具有极强的非定常特征;(4)当进气道处于不起动状态时,全机的气动力特性呈周期性变化,升阻比的变化幅度较大,在最大、最小值分别可达起动状态的2倍和1/4倍。另外,升力、阻力系数的变化曲线之间存在一定的相位差;(5)与实验结果的对照表明,所采用的数值仿真方法具有较高的精度。