改进的万有引力公式

根据广义相对论给出的行星运动方程，导出如下改进的万有引力公式F＝－GMm/r^2-3G^2M^2mp/c^2r^4应用该公式可给出光线近日偏折问题的准确解。     

点火发动机形成的受限喷流流场数值模拟

从曲线坐标下的气体Ｎ－Ｓ方程出发，采用ＥＮＯ差分格式，对点火发动机工作初期，在固体火箭发动机内腔形成的受限冷喷流流场进行数值模拟，在此基础上，分析了喷流的流场结构及特性。     

一类特殊的一阶常微分方程的初等积分法

微分方程是与微积分同时诞生的，从诞生的时候起，人们就试图用初等积分法求微分方程的解。经过一百多年的努力，人们掌握了一些可以用初等积分法求解的微分方程类型，也发现有许多微分方程无法用初等积分法求解。但还有很多微分方程可用初等积分法求解，只是人们还没有发现它们。所以，用初等积分法求解微分方程仍是微分方程领域的一个基本研究内容。在这篇文章里讨论了一类类似黎卡提方程的一阶微分方程的初等积分法。     

基于SDRE方法的空间飞行器轨道交会状态估计

针对空间飞行器自主交会对接的最终逼近阶段,采用雷达作为相对测量设备,建立了目标飞行器轨道坐标系下的相对运动方程以及基于星载雷达的测量模型;运用基于状态相关黎卡提方程(SDRE)的方法设计滤波器,对目标飞行器运动状态进行估计,给出了状态相关黎卡提方程的数值解法,通过具体仿真算例分析了所设计滤波器的性能,并与扩展卡尔曼滤波器进行了比较,结果表明所提出的滤波器模型具有更高的相对状态确定精度.     

曲边单元间断有限元方法求解二维Euler方程

考虑到间断有限元方法对边界的敏感性,采用基于八节点曲边四边形单元的间断有限元方法求解了Euler方程的圆柱绕流问题。详细推导了八节点四边形单元的变换关系,给出了Jacobi矩阵行列式的具体表达式。对比直边单元和曲边单元的计算结果,采用曲边单元后,计算结果符合Euler方程的无粘假设,得出了八节点四边形单元间断有限元方法求解Euler方程是合适的结论。     

Stability-based SDRE controller for spacecraft momentum management

Momentum management of spacecraft aims to avoid the angular momentum accumulation of control momentum gyros through real-time attitude adjustment. An attitude control/momentum management controller based on state-dependent Riccati equation is developed for attitude-stabilized spacecraft. The governing equations of the system are formulated as three-axis coupled with full moment of inertia, which fully capture the nonlinearity of the system and are valid for systems with significant products of inertia or strong pitch to roll/yaw coupling. The state-dependent Riccati equation algorithm brings the nonlinear system to a linear structure having state dependent coefficients matrices and minimizing a quadratic-like performance index. The system equations are nondimensionalized, which avoid numerical problems at the same time make the weighting matrix more predictable. To guarantee closed-loop system stability, the state-dependent Riccati equation algorithm is also modified based on pole placement technique. The state-dependent Riccati equation is online calculated through the computational-efficient θ-D technique which reaches a tradeoff between control optimality and computation load. The dynamic characteristics of the system at torque equilibrium attitude are analyzed. Constraints on moment of inertia for successful momentum management are provided. Simulations demonstrate the excellent performance of the controller.     

利用二、三维嵌套技术数值模拟复杂边界下的流场

通过数值模拟复杂几何区域和边界条件下二、三维嵌套紊流模型,获得了计算区域内任意一点的流场信息。这一方法克服了常规的单纯模拟二维或三维紊流控制方程组而无法获得准确的自由表面和三维紊流场的缺点,通过两者有机的结合,并采用动边界处理等先进的计算技术,得出了精度更高、更为合理的紊流场。该方法曾被用于三峡工程大江截流完成之后导流明渠通航的流场模拟中,计算结果合理可信,部分成果的精度也已得到了实测资料的验证,成果为导流明渠通航路线的选择提供了依据     

无轴涵道旋翼气动特性数值研究

采用基于非结构网格的滑移网格技术对无轴涵道旋翼与传统涵道旋翼的气动特性进行了非定常Euler方程数值模拟。分别考察了无轴涵道旋翼与传统涵道旋翼气动特性差异,中心孔径、涵道扩散角、涵道翼型对无轴涵道旋翼气动特性的影响。涵道扩散角变化范围为-6°~10°,涵道翼型选择NACA66、NACA0018及NACA4415三种翼型进行研究。研究发现:涵道翼型类型对无轴涵道旋翼的拉力分配影响较大,对称翼型能减弱旋翼上方低压涡从而涵道能产生更大的拉力;无轴涵道旋翼比传统涵道旋翼具有更优的拉力性能,转速为18 000 r/min时其总拉力是后者的1.185倍;减小无轴涵道旋翼的中心孔径能提高总拉力值,但整体耗能也将随之提高;涵道拉力占比越高的翼型,其最佳涵道扩散角越大,该状态下涵道拉力占比最高。      

一种统一惯性导航方法

捷联式和平台式惯性导航系统的导航方程不同:同一惯性测量系统选取不同导航坐标系时,导航方程也不同;各类力学编排形式多样、结构复杂且输出结果具有局限性.针对这些问题,提出一种统一惯性导航方法.惯性导航试验结果表明,基于统一惯性导航方程与基于传统惯性导航方程解算的结果一致,验证了统一惯性导航方程的有效性.     

三维仿生悬停扑翼的时间非对称扑动气动特性

基于对现有黄蜂飞行观察实验的分析,建立了不同下扑和上挥时间非对称扑动悬停扑翼运动轨迹模型,并采用数值求解三维层流Navier-Stokes(N-S)方程的方法,研究了采用不同下扑和上挥时间比的仿生悬停扑翼气动特性.结果表明采用适当的时间非对称扑动可以增强悬停扑翼的气动性能.进一步通过对不同时间非对称扑动扑翼流场分析得出,采用时间非对称扑动不但可以增强使扑翼产生高升力的旋转环流机制,而且还可以增强扑翼表面涡流和展向流的强度,从而使扑翼具有更好的气动性能.