基于绳系并联机器人支撑系统的SDM动导数试验可行性研究

详细给出了在低速风洞中,采用绳系并联机器人(WDPR)支撑模型,用强迫振荡法进行标准动态模型(SDM)动导数试验可行性的研究。试验中将杆式六分量应变天平内置入模型中以测量模型的气动力和气动力矩,建立了适用于绳系并联机器人支撑系统的模型运动控制子系统和数据采集子系统。采用绳拉力作为参考信号,对气动力矩信号与位姿信号进行数据的同步处理,解决了绳系并联机器人支撑系统应用于动导数试验时所测力矩信号与位姿信号之间的相位差确定问题,给出了WDPR支撑下模型动导数的计算方法。整个试验样机置于某开口式低速直流风洞中进行了俯仰、带偏航角的俯仰以及升沉的动导数试验,通过测量和计算得到各动导数。试验结果与参考文献相比较具有合理的一致性。研究结果表明,采用绳系并联机器人支撑模型进行动导数试验是可行的,至少对于SDM是这样的结果;使用一套绳系并联机器人支撑系统,可以完成多套硬式支撑系统才能完成的动导数试验,从而提高试验效率,降低试验成本。     

Numerical Simulation of Rotor Flow Field Based on Overset Grids and Several Spatial and Temporal Discretization Schemes

A numerical method based on solutions of Euler/Navier-Stokes (N-S) equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method (AUSM) scheme, upwind AUSM/van Leer scheme, AUSM+ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points (IGBPs). Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.     

运输机空投货物建模与任务性能影响因素研究

针对当前运输机空投货物模型可拓展性不强的问题,提出了分离体的建模方法。在无需分析货物与飞机整体质心变化的情况下,分别以货物和飞机为对象进行受力分析,推导了货物与飞机之间的相互作用力,建立了更加贴近实际的运输机单件货物空投和连续空投的动力学模型。最后,通过数值仿真,分析了飞行速度和牵引伞牵引比对空投任务性能的影响。结果表明,飞行速度和牵引比越大,飞机在货物干扰作用下的状态变化量越小。     

组合升力体动力学建模与仿真

设计了组合升力体总体布局，提出了新的组合升力体建模思路与方法，建立了完整的组合升力体飞行动力学数值仿真模型。建模时综合考虑了主要气动力部件的建模，详细讨论了旋翼／固定翼过渡飞行模式的建模方法。根据提出的建模方法进行了组合升力体的操纵响应分析，仿真结果证明了建模理论和方法的合理性和有效性。     

RELATIONSHIPS FOR MOTION DESCRIBED IN ANY TWO ROTATIONAL REFERENCE FRAMES

In this paper the definitions of the relative derivatives of a vector with respect to time in any reference frame by using the notation of vectrix are proposed. It is further shown that with the help of vectrix one can treat the relative derivatives more clearly and more rigorously when multiple reference frames are concerned. Very useful relationships between the first order relative derivatives described in any two reference frames are derived. Relationships for the second order relative derivatives are also presented in detail. A new concept called "rotation induced acceleration" which is caused by the rotation of two reference frames when motion is necessarily described in them is presented. The result obtained in this paper is the generalization of the acceleration theorem given in conventional theoretical mechanics.     

Structure and kinetic properties of slow-mode shocks associated with magnetic reconnection in the near-Earth magnetotail

We study the structure and kinetic properties of slow-mode shocks near the plasma sheet boundary layer (PSBL) associated with magnetic reconnection by Cluster observation. The presence of slow-mode shocks is confirmed by traditional Rankine–Hugoniot (RH) analysis and Monte-Carlo shock fitting method. The Walén analysis, applied to the tailward flow associated with slow-mode shocks, also supports that plasma was accelerated across a Petschek-type slow-mode shock connected to the diffusion region. Back-streaming ions were observed on the shock layer, and cold ions were accelerated and heated by slow-mode shocks. In addition, whistler and electrostatic solitary waves were observed around the slow-mode shocks. These waves might be excited by the observed field-aligned electron beams near the shocks.     

计算大型实对称特征问题的 Lanczos-QR 算法

为了计算大型实对称特征值问题Ｋｘ＝λＭｘ的少数低阶特征值对，本文给出Ｌａｎｃｚｏｓ－ＱＲ迭代方法。首先，给定初始迭代向量ｖ１，作ｍ步Ｌａｎｃｚｏｓ分解：ＫＶｍ＝ＭＶｍＴｍ＋ｈｍｅｍＴ。取Ｔｍ的ｄ个最大特征值为移位量，对Ｔｍ进行ｄ步带原点位移的ＱＲ分解。然后，修改初始迭代向量ｖ１。迭代地重新开始这一过程，迫使初始迭代向量ｖ１进入需求的特征子空间，从而使残量‖Ｋｘ－θＭｘ‖→０。数值例子表明，该方法收敛性强，且稳定、有效。     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

拖靶系统挂机安装角度计算

拖靶安全释放是拖靶系统正常执行任务的前提,拖靶挂机安装角度的优劣直接影响着拖靶的安全释放.依据拖靶释放特性建立了拖靶系统动力学模型,提出了拖靶安全释放的条件,给出了拖靶挂机安装角度的计算方法;通过实例对方法的正确性进行了验证,其结果直接应用于某新型拖机改装试验中,取得了满意的效果.     

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机，针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响，对全机带动力风洞试验模型进行了大规模并行非定常数值计算，再现了风洞试验现象，通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂，这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速，为保证飞行安全，继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示，在飞行雷诺数和螺旋桨转速下，相同侧风范围内风洞试验显示的流动不稳定因素基本消失，偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析，建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系，进行了偏航稳定性的虚拟试飞评估，研究成果也为AG600飞机的首飞和飞行试验所验证。