一种新的空间飞行器再入制导律研究

基于阻力与能量关系的优化剖面,以该再入剖面作为基准轨迹,给出了跟踪这一基准轨迹的空间飞行器再入制导律与控制律,将基准轨迹转化为相应的迎角剖面和滚转角剖面,设计了相应的结构配置。仿真结果表明,基于阻力与能量关系剖面设计的制导律满足了各项约束条件,并能成功地控制飞行器到达所要求的位置。     

THE RESEARCH OF INCREASING LIFT ANDREDUCING DRAG FOR WING-TIP SAILS

ＴＨＥＲＥＳＥＡＲＣＨＯＦＩＮＣＲＥＡＳＩＮＧＬＩＦＴＡＮＤＲＥＤＵＣＩＮＧＤＲＡＧＦＯＲＷＩＮＧ－ＴＩＰＳＡＩＬＳＣｈｅｎＭｉｎｇｙａｎ，ＱｉＭｅｎｇｂｕ（ＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕｔｉｃｓ，...     

数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

某型飞机阻力伞抛放系统故障模式分析

针对某型飞机阻力伞系统放伞、抛伞故障,通过分析系统的组成、工作原理、常见故障模式以及故障原因,提出了改进建议,为该型飞机阻力伞抛放系统故障的排除提供思路。     

基于改进准则法的复合壳阻尼结构拓扑减振动力学优化

针对复合壳阻尼结构的拓扑减振优化问题，以约束阻尼层的有限单元为设计变量，采用体积比、模态频率和振型为优化约束条件，构建以多模态权重系数的结构模态损耗因子数值关系为优化目标函数的拓扑减振优化模型。为了拓展优化目标灵敏度具有不局限于某一变密度法插值模型的形式，推导了数值表达式的一般函数式。动力学优化中优化目标灵敏度正、负数集共存，使得非凸性的目标函数设计变量出现负值或优化函数寻优于局部极值点。为此，推导出复合壳阻尼结构的全域灵敏度改进优化准则法迭代格，以确保每次迭代域均为全域设计变量集。结合有限单元法编程实现了复合壳阻尼结构改进准则法，并对复合壳结构进行拓扑减振优化分析。结果表明：在敷设体积减为全覆盖的50%时，复合壳结构的模态损耗因子增减偏差为10%，具有提升减振的轻量化设计目的；各阶目标函数和拓扑构型所需的迭代次数少，中间密度区域较小，多阶优于单阶模态优化函数，易于获得全域寻优的有效减振。     

Laminar flow control research at TsAGI: Past and present

This paper presents a brief review of activities in laminar flow control being performed at the Central Aerohydrodynamic Institute named after Prof. N.E. Zhukovsky (TsAGI). These efforts are focused on the improvement of the existing laminar flow control methods and on the development of new ones. The investigations have demonstrated the effectiveness of aircraft surface laminarization applications with the aim of friction drag reduction. The opportunity of considerable delaying of laminar-turbulent transition due to special wing profile geometry and using boundary layer suction and surface cooling has been verified at sub- and supersonic speeds through various wind tunnel testing at TsAGI and during flying laboratory experiments at the Flight Research Institute (LII). The investigations on using hybrid laminar flow control systems for friction drag reduction were also carried out. New techniques of laminar flow control were proposed, in particular, the method of local heating of the wing leading edge, boundary layer laminarization by means of receptivity control, and electrohydrodynamic methods of boundary layer stability control.     

细长杆降低超声速客机气动噪声的数值分析

在超声速客机机头加装适当的细长杆可以降低飞行噪声。本文借助基于AUSM+格式的准三维数值模拟手段对7组13种细长杆方案做了分析比较,总结出了近场气体参数与远场气体参数之间的关系的规律。文章还从气体动力学理论和激波理论出发,对数值模拟结果作了解释,并揭示了细长杆降噪效果与细长杆外形的内在联系,以及近场、远场两者关系的经验公式。文章从各种方案中选取了降噪效果较好的细长杆方案做三维流场分析,并通过流固耦合计算对其结构强度做了校核,验证了方案的可行性。     

水气动量交换的实验研究

在风浪槽中进行水气动量交换实验，测量了不同风速和风区下的波浪及空气湍流流场。测量数据表明水面上风速分布符合对数律，水气界面附近在小范围的动量交换常通量层。阻力系数ＣＤ随风速增加线性增长，随波龄增加而减小。这一趋向与理论分析及外海帝测结果一致。对风浪发展过程中波浪和空气湍流谱结构的分析证明了水气运动的耦合关系。     

Global aerodynamic design optimization based on data dimensionality reduction

In aerodynamic optimization, global optimization methods such as genetic algorithms are preferred in many cases because of their advantage on reaching global optimum. However, for complex problems in which large number of design variables are needed, the computational cost becomes prohibitive, and thus original global optimization strategies are required. To address this need, data dimensionality reduction method is combined with global optimization methods, thus forming a new global optimization system, aiming to improve the efficiency of conventional global optimization. The new optimization system involves applying Proper Orthogonal Decomposition (POD) in dimensionality reduction of design space while maintaining the generality of original design space. Besides, an acceleration approach for samples calculation in surrogate modeling is applied to reduce the computational time while providing sufficient accuracy. The optimizations of a transonic airfoil RAE2822 and the transonic wing ONERA M6 are performed to demonstrate the effectiveness of the proposed new optimization system. In both cases, we manage to reduce the number of design variables from 20 to 10 and from 42 to 20 respectively. The new design optimization system converges faster and it takes 1/3 of the total time of traditional optimization to converge to a better design, thus significantly reducing the overall optimization time and improving the efficiency of conventional global design optimization method.     

Refining DORIS atmospheric drag estimation in preparation of ITRF2008

In preparation of ITRF2008, all geodetic technique services (VLBI, SLR, GPS and DORIS) are generating new solutions based on combination of individual analysis centers solutions. These data reprocessing are based on a selection of models, parameterization and estimation strategy unique to each analysis center and to each technique. While a good agreement can be found for models between groups, thanks to the existence of the IERS conventions, a great diversity still exist for parameter estimation, allowing possible future improvements in this direction. The goal of this study is to focus on the atmospheric drag estimation used to generate the new DORIS/IGN ignwd08 time series prepared for ITRF2008. We develop here a method to inter-compare different processing strategies. In a first step, by analyzing single-satellite solutions for a few weeks of data but for a large number of possible analysis strategies, we demonstrate that estimating drag coefficient more frequently (typically every 1–2 h instead of previously every 4–8 h) for the lowest DORIS satellites (SPOTs and Envisat) provides better geodetic results for station coordinates and polar motion. This new processing strategy also solved earlier problem found when processing DORIS data during intense geomagnetic events, such as geomagnetic storms. Differences between drag estimation strategies can mostly be found during these few specific periods of extreme geomagnetic activity (few days per year). In such a case, when drag coefficient is only estimated every 6 h or less often for single-satellite solution, a significant degradation in station coordinate accuracy can be observed (120 mm vs. 20 mm) and significant biases arose in polar motion estimation (5 mas vs. 0.3 mas). In a second step, we reprocessed a full year of DORIS data (2003) in a standard multi-satellite mode. We were able to provide statistics on a more reliable data set and to strengthen these conclusions. Our proposed DORIS analysis is easy to implement in all software packages and is now already used by several analysis centers of the International DORIS Service (IDS) when submitting reprocessed solutions for ITRF2008.