气动力辅助变轨的数值模拟

采用了航天器在行星上层大气中进行高超声速飞行时的轨道动力学方程,针对航天器从地球静止轨道转移到一个共面圆形低地轨道的变轨过程,进行了气动力辅助变轨过程的模拟.在变轨过程中,航天器从地球静止轨道开始,经过8次大气路径,耗时43.7小时,到达圆形低地轨道,与霍曼转移进行对比,其所消耗的推进剂质量仅为霍曼转移的41%.研究结果表明:气动力辅助变轨技术能够在降低推进剂消耗的情况下实现航天器的轨道转移.     

超声速流中有限长柱壳蒙皮颤振的精确线性理论

本文讨论了有限长柱壳超声速颤振的线性问题的数学基础,包括问题的提法以及气动力和颤振振型的精确求法;所采用的方法是作者在前文[11]中提出的方法的推广. 考虑到壳面气动力的非定域性以及壳体的精确力短理论,根据严格的线性理论,将柱壳蒙皮颤振问题化成了讨论一个新型的八阶非自共轭微分积分方程以及相应的八个边界条件所构成的复固有值问题:     

膜扑翼飞行器的变形研究

最近昆虫翼的变形成了研究热点,而扑翼飞行器的变形力学研究却几乎无人问津.然而,无论昆虫、鸟类还是扑翼飞行器在飞行时,翼的变形都是存在的,要精确计算翼扑动产生的气动力,必须考虑其变形.本文比较了导致变形产生的膜扑翼飞行器的惯性力和气动力在一个周期中的变化情况,发现它们的峰值比值为2左右,然后提出了在随体坐标系中的固支边界条件,采用有限元法计算了惯性力和气动力分别对变形的影响,发现扑翼飞行器的气动力对变形的影响是不可忽略的重要因素,而惯性力与气动力的合力引起的最大正变形发生在下扑初始阶段,最大负变形发生在上扑初始阶段.本文为扑翼飞行器的设计提供了力学分析基础.     

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

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

《力学进展》总目录2005年第35卷第1(总第141)期—第4(总第144)期

(括弧外数字为期数，括弧内数字为页数) 综述评论 扩展有限元法(XFEM)及其应用.……，…、，…，……卜...……，.…，.............……李录贤王铁军1(5) ALE迎风有限元法研究进展.，.......……，.........……，..……、........……岳宝增彭武王照林1(21) 气动力辅助变轨与最优控制研究进展..…，........……，……，.，……，...……‘.吴德隆彭伟斌王小军1(30) 齿轮系统非线性振动研究进展.，.……，...……，.............................……王建军李其汉李润方l(37) 材料力学性能退化的超声无损检测与评价...................…     

文章导读

正飞机结构气动弹性分析与控制研究(1–27,doi:10.6052/0459-1879-15-423)胡海岩,赵永辉,黄锐评述飞机结构气动弹性研究进展并重点介绍作者的如下研究成果:在亚音速气动弹性研究中,揭示多个控制面间隙、反馈时滞对颤振的影响机理,提出时滞最优控制方法,经风洞实验证明可在较大流速范围内抑制颤振;在跨音速气动弹性研究中,提出非线性气动力的降阶模型,进而有效提高了跨音速颤振分析和颤振主动抑制的效率.     

高超声速风洞短时气动力智能辨识算法研究

风洞测力试验是高超声速飞行器研发的重要环节,随着这项技术的发展,试验模型的大尺度化成为高超声速风洞试验的趋势.在几百毫秒的有效测试时间内,大尺度测力系统刚度减弱等问题会严重导致气动力辨识精度变差,试验模型大尺度化对短时脉冲燃烧风洞精确气动力辨识带来了挑战.对此本文提出了一种新的基于传统信号处理结合深度学习的智能气动力辨...     

平流层飞艇空气动力估算

本文采用有限基本解方法与工程估算方法相结合的气动力工程计算方法,用以计算平流层飞艇的气动力。将飞艇所受的气动力分成飞艇艇身和尾翼所受气动力两部分,每一部分的气动力按照无粘性流产生的线性气动力和粘性引起的非线性气动力分别进行计算。根据势流理论对飞艇艇身线性气动力进行分析计算,由于飞艇艇体是旋成体,故根据Allen的横流阻力理论对其所受的非线性气动力进行计算;尾翼的线性气动力采用有限基本解方法进行计算,非线性气动力用Polhamus-Lamar吸力比拟方法估算。该方法中考虑了由于尾翼安装在体上后,处于艇体产生的上洗流场中,尾翼气动力的变化和尾翼对艇身气动力的干扰作用。通过算例的计算与实验结果比较得出该方法可以快速、准确的计算飞艇所受的气动力。     

侧向多喷口干扰复杂流动数值模拟研究

采用具有高分辨率的NND格式,通过数值求解N-S方程对典型外形多喷口侧向喷流复杂干扰流动进行了数值模拟. 为了提高计算效率,采用了LU-SGS隐式算法. 采用分块对接网格技术,生成高质量的贴体计算网格,精确模拟喷口截面. 对比分析了不同计算格式、限制器形式、网格拓扑及流动形态(层流与湍流)对喷流干扰流场结构和压力分布特性的影响,研究和分析了喷口附近流场的涡系结构、波系结构和喷流干扰引起的气动力特性. 在上述研究的基础上,针对典型飞行器外形的侧向喷流干扰特性进行了详细的数值模拟,得到了喷口参数(喷口位置、数目等)及来流条件对喷流干扰流场结构、气动力特性的影响规律,并对其流动机理进行了相应的分析. 研究表明,发展的针对多喷口侧喷干扰的数值计算方法是成功的,可以应用于飞行器侧向喷流干扰的流场结构分析及气动力特性数值预测.      

基于ROM技术的阵风响应分析方法

阵风响应分析是大型飞机设计过程中必不可少的环节. 现有的阵风响应分析主要采用基于线化升力面理论的气动力模型,不能考虑到各种非线性效应,不适合于跨音速气动弹性的分析. 基于CFD技术,采用系统辨识方法,在状态空间内建立了降阶的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程、非定常气动力模型(结构运动)、外激阵风的气动力模型,建立了基于CFD技术的阵风响应分析模型.算例研究了某一典型机翼在方波激励下的阵风响应问题,对比了各阶模态位移的响应以及翼根弯矩的响应. 基于ROM技术的计算结果与CFD/CSD直接耦合仿真结果吻合,证明了该方法的正确性和精度.      

三体轨道动力学研究进展

三体系统轨道动力学问题是航天动力学领域中的经典问题, 具有丰富的理论与工程意义, 并将在人类由近地延伸到深空的航天活动过程中起到至关重要的作用. 本文回顾并总结了三体系统轨道动力学相关研究进展, 并结合未来的深空探测的发展趋势, 展望了三体系统轨道动力学研究中的热点与挑战. 首先阐述了三体问题的研究背景及意义, 简要回顾了三体系统动力学模型的发展历程. 其次, 系统概述了三体系统平衡点附近的局部运动特性, 介绍了平衡点附近周期轨道解析与数值求解方法, 给出了拟周期运动的最新进展. 同时总结了共振轨道、循环轨道、自由返回轨道等三类三体系统全局周期运动的动力学特性与研究进展. 再次, 从不变流形理论和弱稳定边界理论两个方面综述了三体系统中低能量转移与捕获轨道设计的研究进展. 最后, 综述了三体系统轨道动力学在编队飞行、导航星座设计两方面的应用, 并展望了全月面覆盖轨道设计、三体系统下的小推力轨道优化和三体系统的三角平衡点开发利用中值得关注的轨道动力学与控制问题.      

State-of-the-art and prospects for orbital dynamics and control near small celestial bodies

小天体探测是未来深空探测的重点领域之一, 而小天体附近轨道动力学与控制问题是小天体探测任务迫切需要解决的关键问题. 该问题涉及形状不规则小天体附近的动力学环境建模与小天体附近轨道动力学机理. 本文从不规则形状小天体引力场的建模、小天体附近的自然轨道动力学、小天体附近的受控轨道动力学3 个方面综述了小天体附近轨道动力学与控制的研究现状与发展趋势, 并分析了小天体附近轨道动力学所面临的挑战与难题, 最后对我国未来小天体探测任务可能涉及的轨道动力学与控制问题的发展方向进行了展望.     

Stability and control of tethered satellite with chemical propulsion in orbital plane

The tethered satellite with chemical propulsion has broad application prospects in the disposal of abandoned satellites, the orbital rescue of spacecrafts, and the transportation of space supplies, which is completely different from the traditional applications of tethered satellites. Therefore, new research on its dynamics, stability, and control becomes useful and interesting. In this article, based on a dumbbell model of tethered satellite, the dynamics equations of tethered system in orbital maneuvering are established. Furthermore, according to the definitions of transversal and radial propulsive coefficients, analytical solutions of the equilibrium position for librational angle are derived during maneuvering in orbital plane; meanwhile, the effects of propulsive coefficients on librational stability are analyzed, which provides a basis for a selection of expected attitude trajectory. Then, a method of hierarchical sliding-mode tension control is presented to track the expected in-plane angle. This method can address the underactuated problem of tethered systems without either complex coordinate transformation for the system state model or constraint equation restrictions. During orbital flight, in-plane and out-of-plane angles are decoupled, so the tether tension control cannot be conducted to inhibit the out-of-plane angle. To solve this problem, the binormal component of thrust acceleration normal to the orbital plane is adopted as a control variable, and a feedback linearization-based thrust controller is designed to damp out the out-of-plane angle. Afterwards, orbital transfer cases between two circular orbits are studied to demonstrate the effectiveness of the tethered satellite with chemical propulsion. Numerical simulation results indicate that the stability of librational angles has a close relation to propulsive coefficients, and distributions of stable centers and unstable saddle points are totally different on both sides of bifurcation point. In addition, tracking control requirements for tethered satellite are guaranteed by designed controllers, which ensure flight safety in orbital maneuvering.     

A new method to retrieve aerosol optical thickness from satellite images on a parallel system

A wide variety of algorithms have been developed to monitor aerosol burden from satellite images. Still, few solutions currently allow for real-time and efficient retrieval of aerosol optical thickness （AOT）, mainly due to the extremely large volume of computation necessary for the numeric solution of atmospheric radiative transfer equations. Taking into account the efforts to exploit the SYNergy of Terra and Aqua Modis （SYNTAM, an AOT retrieval algorithm）, we present in this paper a novel method to retrieve AOT from Moderate Resolution Imaging Spectroradiometer （MODIS） satellite images, in which the strategy of block partition and collective communication was taken, thereby maximizing load balance and reducing the overhead time during inter-processor communication. Experiments were carried out to retrieve AOT at 0.44, 0.55, and 0.67μm of MODIS/Terra and MODIS/Aqua data, using the parallel SYNTAM algorithm in the IBM System Cluster 1600 deployed at China Meteorological Administration （CMA）. Results showed that parallel implementation can greatly reduce computation time, and thus ensure high parallel efficiency. AOT derived by parallel algorithm was validated against measurements from ground-based sun-photometers; in all cases, the relative error range was within 20%, which demonstrated that the parallel algorithm was suitable for applications such as air quality monitoring and climate modeling.     

The present status and prospects in the research of orbital dynamics and control near small celestial bodies

Small celestial body exploration is of great significance to deep space activities. The dynamics and control of orbits around small celestial bodies is of top priority in the exploration research. It includes the modeling of dynamics environment and the orbital dynamics mechanism. This paper introduced state-of-the-art researches, major challenges, and future trends in this field. Three topics are mainly discussed: the gravitational field modeling of irregular-shaped small celestial bodies, natural orbital dynamics and control, and controlled orbital dynamics. Finally, constructive suggestions are made for China's future space exploration missions.     

自由漂浮空间机器人路径优化的Legendre伪谱法

基于Legendre 伪谱法研究自由漂浮空间机器人非完整路径规划的最优控制问题. 自由漂浮是空间机器人执行任务常用的工作模式,其路径优化是空间机器人完成复杂空间任务的基础. 由于空间机器人不具有固定基座,机械臂和载体之间存在非完整约束,使得自由漂浮空间机器人路径规划完全不同于地面机器人而变得具有挑战性. 本文提出自由漂浮空间机器人路径规划的最优控制伪谱方法. 首先,利用多体动力学理论建立自由漂浮空间机器人动力学模型,给定系统的初始和目标位形,选取机械臂关节耗散能最小为性能指标,并考虑实际控制输入受限,建立其路径规划的Bolza 问题. 然后,应用Legendre 伪谱法,将状态和控制变量在Legendre-Gauss-Lobatto (LGL) 点上离散,并构造Lagrange 插值多项式逼近系统状态和控制变量,将连续路径优化问题离散化为非线性规划问题求解. 最后通过数值仿真表明,应用Legendre 伪谱法求解自由漂浮空间机器人非完整路径规划问题,得到的机械臂和载体最优运动轨迹,较好地满足各种约束条件,且计算精度高、速度快,并具有良好的实时性.      

Natural convection of an alumina-water nanofluid inside an inclined wavy-walled cavity with a non-uniform heating using Tiwari and Das’ nanofluid model

The present study is devoted to numerical analysis of natural convective heat transfer and fluid flow of alumina-water nanofluid in an inclined wavy-walled cavity under the effect of non-uniform heating. A single-phase nanofluid model with experimental correlations for the nanofluid viscosity and thermal conductivity has been included in the mathematical model. The considered governing equations formulated in dimensionless stream function, vorticity, and temperature have been solved by the finite difference method. The cavity inclination angle and irregular walls (wavy and undulation numbers) are very good control parameters for the heat transfer and fluid flow. Nowadays, optimal parameters are necessary for the heat transfer enhancement in different practical applications. The effects of the involved parameters on the streamlines and isotherms as well as on the average Nusselt number and nanofluid flow rate have been analyzed. It has been found that the heat transfer rate and fluid flow rate are non-monotonic functions of the cavity inclination angle and undulation number.     

Method of Optimizing Maneuvers of an Orbital Transfer Vehicle in a Strong Central Field of Gravity

International Applied Mechanics - A method of constructing the optimal trajectories of interorbital flights of an orbital transfer vehicle with a low thrust propulsion system is described. Testing...     

Low-thrust Lambert transfer based on two-stage constant-vector thrust control method

The low-thrust Lambert transfer refers to that the spacecraft achieves the orbital transfer whose boundary conditions are represented by two sets of orbital elements at initial and final time by the low-thrust propulsion system. The modulus and direction of the low-thrust solutions in previous methods change with time, which leads to high control requirements for the engine. In this paper, to reduce the requirements of the engine, a practical two-stage constant-vector thrust control method is proposed, in which the magnitude and direction of the thrust are deemed as segmental constant value in TNH frame, where three components of the thrust are f_t, f_n, and f_h. First, the mathematical model of the two-stage constant-vector thrust is formulated, and a rapid algorithm is presented to obtain the solution based on the linearized sensitivity matrix, which describes the relationship between the constant-vector thrust and the change of the orbital elements approximately. Furthermore, two low-thrust Lambert strategies based on the two-stage constant-vector thrust are presented for cases of short-time transfer and long-time transfer. A sequence of numerical simulations demonstrated the efficiency of the proposed approaches. The proposed control strategies are solved rapidly, and they are also suitable for different types of orbits with J₂ perturbation, which are practical options for engineering applications.

     

Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.