基于最小值原理的冗余度机器人动力学控制

为了解冗余度机器人全局法优化中数值求解的困难，本文讨论了动力学方程的建立，无约束和有约束最优控制问题之间的内在联系，重点分析了求解最优控制问题的数值方法，文中提出了双向异步积分迭代的求解正则方程组的直接迭代法，较好解决了状态方程和协态方程稳定相逆给求解两点这值问题带来的困难。     

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

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

基于图剖分的多块结构网格负载平衡方法

负载平衡是影响并行计算性能的重要因素。针对多块结构网格,给出了一种改进的多层次图剖分负载平衡方法。该方法设计了新的网格剖分算法,采用改进的子块分裂方法与图剖分算法的循环调用实现结构对接网格剖分,并通过建立不同物体重叠网格间的连接关系,实现了结构重叠网格的负载平衡。采用2个典型算例对方法进行了对比验证,数值结果表明,子块分裂方法对剖分结果具有重要影响,采用循环调用算法及改进的子块分裂方法能有效地实现计算负载均衡及通信量优化,同时显著减少了网格块数及因虚网格导致的内存需求,有利于提高并行效率。该负载平衡方法与网格拓扑无关,适用于多块结构对接网格及重叠网格,且整体型剖分方式对于多块结构重叠网格具有更好的剖分效果。     

Propagation and interaction of interplanetary transient disturbances. Numerical simulations

We study the heliocentric evolution of ICME-like disturbances and their associated transient forward shocks (TFSs) propagating in the interplanetary (IP) medium comparing the solutions of a hydrodynamic (HD) and magnetohydrodynamic (MHD) models using the ZEUS-3D code [Stone, J.M., Norman, M.L., 1992. Zeus-2d: a radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. i – the hydrodynamic algorithms and tests. Astrophysical Journal Supplement Series 80, 753–790]. The simulations show that when a fast ICME and its associated IP shock propagate in the inner heliosphere they have an initial phase of about quasi-constant propagation speed (small deceleration) followed, after a critical distance (deflection point), by an exponential deceleration. By combining white light coronograph and interplanetary scintillation (IPS) measurements of ICMEs propagating within 1 AU [Manoharan, P.K., 2005. Evolution of coronal mass ejections in the inner heliosphere: a study using white-light and scintillation images. Solar Physics 235 (1–2), 345–368], such a critical distance and deceleration has already been inferred observationally. In addition, we also address the interaction between two ICME-like disturbances: a fast ICME 2 overtaking a previously launched slower ICME 1. After interaction, the leading ICME 1 accelerates and the tracking ICME 2 decelerates and both ICMEs tend to arrive at 1 AU having similar speeds. The 2-D HD and MHD models show similar qualitative results for the evolution and interaction of these disturbances in the IP medium.     

二维三角晶格的Peierls相变前后的声子色散曲面的研究

利用晶格动力学理论计算了二维三角晶格的Peierls相变前后的声子能谱。并且在第一布里渊区绘制了声子色散曲面。通过对Peierls相变前后的声子能谱比较可知，随其对称性降低，声子能谱的简并度被破缺。二维三角晶格的Peierls相变前只有两支声频支声子能谱，Peierls相变后有两支声频支声子能谱与两支光频支声子能谱，并且声频支声子能谱与光频支声子能谱不相交，两支光频支声子沿ky轴平缓地变化。     

Numerical prediction and wind tunnel experiment for a pitching unmanned combat air vehicle

The low-speed flowfield for a generic unmanned combat air vehicle (UCAV) is investigated both experimentally and numerically. A wind tunnel experiment was conducted with the Boeing 1301 UCAV at a variety of angles of attack up to 70 degrees, both statically and with various frequencies of pitch oscillation (0.5, 1.0, and 2.0 Hz). In addition, pitching was performed about three longitudinal locations on the configuration (the nose, 35% MAC, and the tail). Solutions to the unsteady, laminar, compressible Navier–Stokes equations were obtained on an unstructured mesh to match results from the static and dynamic experiments. The computational results are compared with experimental results for both static and pitching cases. Details about the flowfield, including vortex formation and interaction, are shown and discussed, including the non-linear aerodynamic characteristics of the vehicle.     

Comparison of mesopause region meteor radar winds,medium frequency radar winds and low frequency drifts over Germany

During 2004 and 2005 measurements of mesospheric/lower thermospheric (80–100 km) winds have been carried out in Germany using three different ground-based systems, namely a meteor radar (36.2 MHz) at the Collm Observatory (51.3°N, 13°E), a MF radar (3.18 MHz) at Juliusruh (54.6°N, 13.4°E) and the LF D1 measurements using a transmitter (177 kHz) at Zehlendorf near Berlin and receivers at Collm with the reflection point at 52.1°N, 13.2°E. This provides the possibility of comparing the results of different radar systems in nearly the same measuring volume. Meteor radar winds are generally stronger than the winds observed by MF and especially by LF radars. This difference is small near 80 km but increases with height. The difference between meteor radar and medium frequency radar winds is larger during winter than during summer, which might indicate an indirect influence of gravity waves on spaced antenna measurements.     

Determination of an instantaneous Laplace plane for Mercury’s rotation

Mercury is the target of two space missions: MESSENGER, which carried out its first and second flybys of Mercury on January 14, 2008 and October 6, 2008, and the ESA/JAXA space mission BepiColombo, scheduled to arrive at Mercury in 2020. The preparation of these missions requires a good knowledge of the rotation of Mercury.     

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

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

轮控小行星探测器建模及跳跃移动仿真

小行星探测器在弱引力环境下无法采用传统的轮式机构进行移动。为解决该问题,采用反作用飞轮对探测器跳跃移动进行控制,并分析了该方案的可行性。根据Hertz碰撞定律及简化的Karnopp切向摩擦力模型,建立了探测器与地面的接触力模型。分析了轮控小行星探测器的起跳过程,给出了探测器静止起跳所需要的最小飞轮力矩关系。考虑到反作用飞轮存在惯性、粘滞、摩擦等情况,建立了轮控探测器的姿态动力学模型,并对探测器在均匀重力场下的连续跳跃过程进行了控制策略设计及仿真。结果表明:基于飞轮控制的小行星探测器跳跃移动在微重力环境下是可行的,且可以通过施加合适的控制力矩维持探测器跳跃的方向及跳跃过程的稳定性。