带落角约束的多导弹分布式协同制导律

为实现多枚导弹协同攻击机动目标,基于具有推力可控能力的导弹,提出了一种带落角约束的多导弹分布式协同制导律。将制导律的设计分离为视线方向和视线法向上两个部分:视线方向上基于多智能体协同控制理论和超螺旋控制算法,设计制导律控制导弹剩余时间在有限时间内趋于一致;视线法向上运用零化视线角速率思想和有限时间滑模控制理论,设计制导律控制导弹击中目标的同时满足落角约束。并针对两部分制导律中存在的目标机动信息,分别设计非齐次干扰观测器进行估计。仿真结果表明,提出的制导律能够有效完成协同攻击任务,脱靶量和落角误差分别控制在0.13 m和0.02°以内,并且有效抑制了抖振现象,有利于提高导弹自动驾驶仪的跟踪精度。     

大长细比导弹的气动弹性降阶模型

对于大长细比导弹,需要在设计阶段准确计算气动弹性/气动伺服弹性,但其复杂的气动力给计算带来困难,因此气动力降阶模型是突破大长细比导弹跨音速气动弹性分析与控制瓶颈的关键技术.虽然气动力模型降阶方法已在预测二维机翼结构的气动弹性方面取得重要进展,但几乎未见关于全机模型的气动力降阶模型研究报道.本文基于递归Wiener模型的气动力降阶方法,利用CFD计算的气动力作为模型辨识数据,用鲁棒子空间和Levenberg-Marquardt算法辨识降阶模型参数,建立了大长细比导弹气动力降阶模型.在此基础上与大长细比导弹有限元模型相结合,构造出气动弹性降阶模型,并在数值仿真中测试气动弹性降阶模型在不同马赫数下的适用性.数值仿真结果表明,该气动弹性降阶模型能够精确预测导弹模型在不同飞行条件下的非定常气动力和导弹模型的气动弹性频率响应特性.     

带有攻击角约束的多导弹协同制导律

针对带有攻击角约束的多导弹同时攻击机动目标问题,提出了一种带有攻击角约束的协同制导律。首先基于平面内的导弹-目标相对运动方程,建立了带有攻击角约束的协同制导模型;其次,把协同制导律的设计过程分离为两个部分:一是基于图论的有关内容,运用有限时间一致性理论设计沿着视线方向上的加速度指令来保证所有导弹与目标的相对距离在有限时间内到达一致,进而保证所有的导弹同时击中机动目标;二是利用非齐次干扰观测器对机动目标的加速度进行估计,并运用滑模控制设计视线法向上的加速度指令来保证每枚导弹与目标间的视线角速率收敛到零和视线角收敛到期望的终端视线角,即每枚导弹以期望的终端视线角成功击中目标;最后,对三枚导弹同时打击同一机动目标的情况进行仿真,仿真结果表明本文设计的带有攻击角约束的协同制导律的有效性和正确性。     

双导弹拦截角度协同的微分对策制导律

为提高命中高价值目标的概率,基于线性二次型微分对策理论,对两枚导弹协同拦截单个目标的制导律进行了研究。单枚导弹在最小化自身脱靶量的同时,与另一枚导弹实现拦截角度上的协同,从而构成特定的拦截态势,以提高拦截机动目标的性能和末制导尾端对目标的可观测性。所推导的微分对策制导律考虑到了对策三方的控制系统动态,且具有解析解,形式上为零控脱靶量和零控协同拦截角误差的线性组合。基于推导结果完成了微分对策制导律的制导增益和对策空间分析,给出了鞍点解的存在条件,并进行了分析。非线性系统仿真结果表明由于导弹间存在显式的协同关系,拦截目标所需的加速度较低,且在设定的协同拦截角度收敛后,加速度会进一步减小。     

导弹纵向俯仰气动特性数值分析

应用有限体积方法求解三维可压缩雷诺平均N-S方程,计算了某型巡航导弹的静、动态绕流流场和空气动力特性.湍流模型应用修正后的B-L代数湍流模型.以细长旋成体绕流流场为例进行数值计算,并和实验数据相比较,验证了本文方法的可靠性,在此基础上,开展了某型巡航导弹绕不同转轴和以不同频率进行俯仰振动的非定常气动力迟滞特性研究.计算结果表明,当导弹作快速俯仰振动时,在上仰和下俯过程中的同一迎角瞬间,绕导弹流场流动明显不同,表现出明显的非定常迟滞特性.并且,导弹的非定常气动力迟滞特性随振动频率的增大变化明显,且气动力迟滞曲线随着振动轴位置的变化而变化.     

火箭橇试验分离制导工具误差的有效性分析

火箭橇试验能在地面提供逼真的飞行环境,是制导系统可靠性检验与精度鉴定的有效手段,代表了未来制导系统测试的一个发展方向。利用火箭橇试验分离制导工具误差是制导系统精度鉴定的核心内容,对此进行了研究:给出了解析式平台制导系统的火箭橇试验方法;推导了制导工具误差系数分离模型;重点对火箭橇试验分离制导工具误差的有效性与精度进行了分析。结果表明,常规解析平台跟踪发射惯性坐标系的情况下,仅能有效辨识少数工具误差系数,系数辨识的精度很大程度上受制于外测精度;进一步分析表明其成因皆由火箭橇试验特性所致。最后给出了三种改善工具误差分离精度的途径。     

带有末端角度和轨迹路径点约束的MPSC制导律设计

针对多约束制导问题,给出了一种同时满足末端角度约束和飞行轨迹路径点约束的模型预测扩展控制制导方法,该制导方法通过满足飞行轨迹路径点约束实现灵活调节飞行轨迹,可以大大缩短目标防御反应时间。模型预测扩展控制制导方法是基于非线性最优控制理论,给出了控制量表达式以二次形式近似时制导律的设计过程。模型预测扩展控制制导方法只能对末端时刻输出量进行约束,通过对该制导方法进行扩展,使其还可以满足飞行轨迹路径点约束。仿真结果表明,考虑飞行轨迹路径点约束时,导弹经过设定的路径点并以给定的弹道倾角命中目标。     

有界控制导弹攻击防御性飞机的制导律

为提高攻击导弹同时面对目标飞机及其防御导弹情况下的命中概率,基于微分对策理论,对攻击导弹的制导律进行了设计。应对独立控制的多对象博弈问题,微分对策理论具有天然的优势,且相比于最优制导律,微分对策制导律对于目标机动估计误差和机动策略具有更强的鲁棒性。所推导的微分对策制导律进一步考虑了攻击导弹的控制有界性,且适用于攻击导弹、目标飞机和防御导弹具有高阶线性控制系统动态的情形。为验证制导律性能,进行了非线性系统仿真,结果表明该制导律在成功归避防御导弹的同时可实现趋于零脱靶量的目标拦截。攻击导弹为实现规避和攻击的双重任务,仅需要保持相比于防御导弹两倍左右的机动优势。     

飞行器纵向俯仰气动特性计算研究

应用有限体积方法求解三维可压缩雷诺平均N-S方程,计算了巡航导弹外形飞行器作小振幅俯仰运动时的动态绕流流场和空气动力特性,开展了导弹绕不同转轴、以不同频率和在不同迎角范围内进行俯仰运动的非定常气动力迟滞特性研究。计算结果表明,当导弹作快速俯仰运动时,在上仰和下俯过程中的同一迎角瞬间,绕导弹流场流动明显不同,表现出明显的非定常迟滞特性。导弹的非定常气动力迟滞特性随俯仰运动频率的增大明显增强,且气动力迟滞曲线随着俯仰轴位置的变化而变化。在同一减缩频率下,导弹在不同迎角范围内作周期俯仰运动时,相同的运动相位角所对应的升力系数对迎角的导数是一致的,而不同减缩频率下升力系数对迎角的导数随运动相位角变化曲线明显不同。     

风洞实验中翼型大迎角气动力分散性的机理分析

基于$k$-$\omega$的SST两方程湍流模型, 求解雷诺平均N-S方程获得定常和非定常气动力, 耦合 翼型弹性振动方程, 在时间域内模拟了不同的翼型非定常流动, 重点研究了大迎角下的分离 流问题, 研究结果表明: 在百万雷诺数条件下, 由于振动引起分离涡的不规律脱落, 可能导 致气动力平均值的变化; 而厚度大于20\%的翼型在一定大迎角范围内, 会出现分离涡流场平 衡态的转化, 从流体力学稳定性的角度, 解释了风洞实验中大迎角气动力数据的分散性, 为 大迎角气动力风洞实验的重复性和数据分散性给出了一种新解释.     

LES Study of Flow Between Shrouded Co-rotating Disks

Large eddy simulations of turbulent flow between shrouded co-rotating disks, representing a simplified model of a hard disk drive, are performed. The computation domain surrounds a complete disk and is bounded at top and bottom by half a disk. Therefore, it is possible to compute the fluctuating pressure field surrounding the middle disk. Also, the influence of the shroud geometry is taken into account by comparing a flat shroud wall and a wall with rib chambers. In the flat shroud case, the fluctuating pressure on the upper- and lower-surface of a disk indicates a strong correlation with fluid motion travelling across the disk-tip clearance region. However, in the ribbed shroud case the organized flow structure that is observed in the flat shroud case disappears and the fluctuating pressure acting on the surface of the disk is remarkably diminished.     

低雷诺数下附面层组合抽吸方案对压气机特性影响的研究

为了分析附面层抽吸流动控制对低雷诺数下压气机特性的影响,本文采用数值方法模拟了低雷诺数下附面层组合抽吸方案对NASA Rotor 37跨音速压气机性能和稳定性的影响特点及作用机理。通过在该压气机转子叶片吸力面和机匣上分别设计附面层抽吸槽,探讨了组合抽吸方案对低雷诺数下(H=20km)压气机性能和稳定性的影响。结果表明:采用组合抽吸方案后,压气机峰值效率提高约1.3%;压气机最大增压比提高约2.5%;压气机转子的近失速点流量减小约14.6%。进一步分析作用机理发现,组合抽吸槽有效抑制了附面层径向涡向叶顶的运动和聚集,使叶顶附面层分离区减少约70%从而有效改善了压气机的流场特性。     

侧压式道气道流场特性研究

应用表面油流和压力测量技术,结合数值模拟方法研究了超燃冲压发动机侧压式进气道无支板和有支板两种典型工况的基本流动特性.结果表明:由侧壁诱发的后掠激波与底板边界层相互作用,形成一对逆向旋转的漩涡,该漩涡在下游不断发展并偏离底板,在出口截面形成一个低马赫数低总压区;支板的引入增加了压缩率,但同时导致底板与侧壁边界层严重分离,使进气道性能偏离设计工况.     

Flow and performance characteristics of twin-entry radial turbine under full and extreme partial admission conditions

This paper presents numerical and experimental investigation of the performance and internal flow field characteristics of twin-entry radial inflow turbines at full and extreme partial admission conditions. The turbine is tested on a turbocharger test facility, which was developed for small and medium size turbochargers. Experimental results show that the lowest efficiency corresponds to extreme conditions. Therefore, flow field analyzing is employed to consider these conditions. The flow pattern in the volute and impeller of a twin-entry turbine is analyzed using an in-house fully three-dimensional viscous flow solver. The computational performance results are compared with the experimental results and good agreement is found. The flow field at the outlet of the turbine is investigated using a five-hole pressure probe; the numerical results are also compared with experimental measurements at the outlet of the rotor. For the volute, results show that lowest entropy gain factor corresponds to the extreme conditions, particularly when shroud side entry is fully closed. At the inlet of the rotor for equal admission conditions, the incidence angle is mostly in the optimum values. However, large variation in the incidence angle is seen in the extreme conditions, which lead to larger incidence losses and consequently a lower efficiency. In addition, entropy distribution contours corresponding to the exit plane are considered. For full admission, the location of low entropy gain factor at this plane occupies a region near the shroud side of suction surface as well as near the hub side of the pressure surface that corresponds to a region of high absolute flow angle. However, for the extreme cases, the low entropy gain factor occupies a relatively larger region near the shroud side than full admission. So, higher loss generation is noted at the extreme cases. Moreover, this entropy gain factor region is increased when shroud side entry is fully closed.     

Numerical simulation of rivulet evolution on a horizontal cable subject to an external aerodynamic field

On wet and windy days, the inclined cables of cable-stayed bridges may experience a large amplitude oscillation known as rain-wind-induced vibration (RWIV). It has previously been shown by in situ and wind-tunnel studies that the formation of rain-water accumulations or ‘rivulets’ at approximately the separation points of the external aerodynamic flow field and the resulting effect that these rivulets have on this field may be one of the primary mechanisms for RWIV. A numerical method has been developed to undertake simulations of certain aspects of RWIV, in particular, rivulet formation and evolution. Specifically a two-dimensional model for the evolution of a thin film of water on the outer surface of a horizontal circular cylinder subject to the pressure and shear forces that result from the external flow field is presented. Numerical simulations of the resulting evolution equation using a bespoke pseudo-spectral solver capture the formation of two-dimensional rivulets, the geometry, location and growth rate of which are all in good agreement with previous studies. Examinations of how the distribution and magnitude of aerodynamic loading and the Reynolds number influence the rivulet temporal evolution are undertaken, the results of which indicate that while all three affect the temporal evolution, the distribution of the loading has the greatest effect.     

Experimental study and numerical simulation of the characteristics of a percussive gas–solid separator

The presence of solid particles in the flow of hypersonic wind tunnels damages the appearance of the experiment models in the wind tunnel and influences the accuracy of experimental results. The design of a highly efficient gas–solid separator was therefore undertaken. Particle trajectory imaging methods were used to measure trajectories under different conditions. The flow field and particle movement characteristics for different head angles (HAs) and separation tooth angles (STAs), inlet velocities, and the exhaust gas outlet pressures in the separator, were calculated using simulations based on the discrete phase model. The particle separation efficiency, pressure loss, and flow loss resulting from different structural parameters were also studied. In line with experimental observations, the characteristic angle of particle movements in the separator and the separation efficiency of the separator were found to increase with decreasing HA and with increasing STA. Separation efficiency improves with increasing inlet velocity and with increasing negative pressure of the exhaust gas outlet; however, the corresponding pressure loss and the flow rate of the waste gas also increased.     

Investigation of flows past a plate in the presence of exoergic processes driven by the interaction between a nonequilibrium electric discharge and a propane-air mixture

Supersonic flow (M = 2) past a plate along which propane is injected is investigated within the framework of the solution of problems of combustion initiation and stabilization at low static temperatures and pressures in the presence of a nonequilibrium discharge and metal and dielectric interceptors mounted on the plate surface. The experiments show that two zones with exoergic reactions develop when a metal interceptor is mounted on the plate. One zone is located ahead of the leading separation zone and the other above and behind the interceptor edge, its boundary partially penetrating into the supersonic flow region. Using modern spectroscopic methods, the radiation intensity distributions of a series of plasmochemical reaction products are obtained in the neighborhood of the plate ahead of the interceptor, behind it, and above its edge. It is found that the fuel is intensively converted under the action of the discharge with the occurrence of a series of free radicals, atomic hydrogen and oxygen which are themselves chemically active.     

Effect of pressure-dependent slip on flow curve multiplicity

Various microstructural pictures for slip at polymer/solid interfaces lead to relations which have a region where multiple values of slip velocity are predicted for the same shear stress. This leads to the expectation of multivalued flow curves, which has been verified in specific cases by numerous researchers. We study the effect of pressure dependence on flow curve multiplicity using a simple multivalued slip relation to model the phenomena of hysteresis and spurt flow in polymer extrusion. A continuation technique is used to trace out the boundaries of the region of flow curve multiplicity as pressure drop and die length to diameter (L/D) ratio are changed. Results for Newtonian, shear thinning and viscoelastic constitutive equations show that, despite the multivalued nature of the slip model, multiplicity (and thus hysteresis) is absent at high L/D.  For the sake of completeness, we also carry out time-dependent simulations at constant piston speed taking fluid compressibility into account. These simulations show that oscillations in the pressure drop and exit volumetric flow rate result only if the system is operated in the multiplicity region of the steady state flow curve, in agreement with the results of similar simulations by researchers using various multivalued slip models without pressure dependence. The results demonstrate that a multivalued slip model does not guarantee multiplicity in the flow curve for the constant pressure drop operation, nor oscillations for constant piston speed operation. Received: 18 August 1997 Accepted: 30 March 1998     

Simultaneous measurements of disk vibration and pressure fluctuation in turbulent flow developing in a model hard disk drive

The complex flow features inside hard disk drive models are investigated in an axisymmetric and a semi-open shroud configurations. For the axisymmetric case, we have employed both experimental and computational approaches. The experiment focuses on both flow dynamics and the disk vibration, where measurements of the fluctuating pressure and velocity are undertaken at some representative points. The correlation between the disk vibration and the fluctuating pressure in the turbulent flow between disks is evident from the spectral analysis. The experimentally observed fluctuating pressure and velocity are partly due to the disk vibration and its contribution could be estimated by comparing the experiment with the results of a large eddy simulation. For the semi-open shroud case, although the characteristic peaks attributable to the large-scale vortical structure are still observed in the power spectra, the pressure fluctuation and the disk vibration are suppressed when the arm is inserted.     

Aerodynamics of laminar separation flutter at a transitional Reynolds number

Experimental observations of self-sustained pitch oscillations of a NACA 0012 airfoil at transitional Reynolds numbers were recently reported. The aeroelastic limit cycle oscillations, herein labelled as laminar separation flutter, occur in the range 5.0×10⁴≤Re_c≤1.3×10⁵. They are well behaved, have a small amplitude and oscillate about θ=0°. It has been speculated that laminar separation leading to the formation of a laminar separation bubble, occurring at these Reynolds numbers, plays an essential role in these oscillations. This paper focuses on the Re_c=7.7×10⁴ case, with the elastic axis located at 18.6% chord. Considering that the experimental rig acts as a dynamic balance, the aerodynamic moment is derived and is empirically modelled as a generalized Duffing–van-der-Pol nonlinearity. As expected, it behaves nonlinearly with pitch displacement and rate. It also indicates a dynamically unstable equilibrium point, i.e. negative aerodynamic damping. In addition, large eddy simulations of the flow around the airfoil undergoing prescribed simple harmonic motion, using the same amplitude and frequency as the aeroelastic oscillations, are performed. The comparison between the experiment and simulations is conclusive. Both approaches show that the work done by the airflow on the airfoil is positive and both have the same magnitude. The large eddy simulation (LES) computations indicate that at θ=0°, the pitching motion induces a lag in the separation point on both surfaces of the airfoil resulting in negative pitching moment when pitching down, and positive moment when pitching up, thus feeding the LCO.