平头增阻再入体俯仰动态特性计算与流动机理分析

采用自由振荡法数值模拟了一种平头增阻再入体的高超声速俯仰动态特性,沿用稳定性导数模型作为动态特性的指标,通过把稳定性导数视为攻角的函数来考虑非线性影响,并应用奇异分解线性最小二乘法辨识稳定性导数。对于本文研究对象,计算结果表明,平头外形在较高马赫数和小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的。文中对平头外形出现这种动态特性现象的空气动力学原因进行了分析,有必要作进一步的确证与研究。     

一种新型三维仿生扑翼机构设计与分析

以昆虫翅膀运动机理为基础,提出了一种复合3种运动的仿生扑翼新机构。由于能使翼尖按8字形或香蕉形运动及机翼扭转,同时可保证两翼的运动对称,使扑翼运动更符合昆虫运动形态。为了实现生物运动轨迹,建立了机构参数的优化模型,并据此设计了扑翼微型飞行器样机,其扑动规律能更逼真地模拟昆虫翼的扑动形态,为新型微型飞行器的设计提供样机。     

Integrated design of topology and material for composite morphing trailing edge based compliant mechanism

The morphing trailing edge based compliant mechanism is a developing technology which can increase lift-drag ratio for variable flight modes by bending down the trailing edge. Composite material design is integrated into topology optimization for the morphing trailing edge based compliant mechanism in the paper. A two-step optimization strategy is established to solve the integrated design problem. Initially, lamination parameters are introduced and viewed as a bridge between structure stiffness and fiber angles for composite material. Design variables include the lamination parameters and element density. The least-squares between actual and desired displacements at output points along trailing edge is adopted to evaluate the deformed capability of the trailing edge. An integrated optimization model for the composite morphing trailing edge is established with the volume constraints. The optimal topologic shape and lamination parameters are initially obtained. Subsequently, a least-squares optimization between fiber angles and the optimal lamination parameters is implemented to obtain optimal fiber angles. Finally, morphing capability of composites trailing edge based compliant mechanism is investigated by simulation and experiments. The results indicate the composites trailing edge based compliant mechanism can approximately bend down 8 degrees and satisfies the design requirement.     

Strong interactions of incident shock wave with boundary layer along compression corner

The coherent structure and instability of the interaction of incident shock wave with boundary layer developing on a compression corner are experimentally studied. The experiments are carried out in a supersonic wind tunnel of Mach number 2. Particular attention is paid to shock patterns and unsteady shock motions induced by the separation bubble. The high-speed schlieren is used to visualize the flowfield evolution and to characterize the instability. The snapshot proper orthogonal decomposition of schlieren sequences is applied to investigate the primary coherent structure in the flowfield. Fast Fourier transform and continuous wavelet transformation are applied to characterize the instability. The results show that there are large-scale low-frequency oscillations of the shock waves and small-scale high-frequency pulsations in the separation region. The peak frequency of shock oscillation is mainly concentrated in the range of 100–1000 Hz. The pulsation of the small flow structure in the separation bubble is mainly concentrated above 12.5 kHz. Based on the results of experimental analysis, the preliminary mechanism of the large-scale instability of such interaction is obtained.     

Positioning error compensation for parallel mechanism with two kinematic calibration methods

Compared with serial mechanisms, the parallel mechanism (PM) theoretically exhibited higher positioning accuracy, dynamic performance, strength-to-weight ratio, and lower manufacturing cost, but they had not been widely used in the practical application. One key issue, positioning accuracy, which directly affected their performance and was greatly influenced by the errors of kinematic structure parameters was analyzed. To effectively enhance the positioning precision of PMs, a novel modeless kinematic calibration method, namely the split calibration, was presented and its compensation effect of the positioning error was comprehensively compared with that of an integrated method on two different types of PMs. A strange phenomenon-correct and incorrect identified results were derived from two different PMs by the same integrated method, respectively-which had not been reported yet was discovered, and the origin of it was revealed utilizing numerical simulations. Finally, respective merits and drawbacks of these two methods obtained in this paper provided underlying insights to guide the practical application of the kinematic calibration for PMs.     

基于力位协同控制的大飞机机身壁板装配调姿方法

为校正中机身壁板由于重力和调姿内力产生的变形，提高中机身壁板装配调姿精度，提出了一种基于力位协同控制的装配调姿方法。通过将调姿机构等效为并联机构，推导了调姿机构的解析正反解模型；根据螺旋理论，建立了力传感器测量值与重力、调姿内力之间的映射关系，实现重力补偿值的动态计算，基于局部刚体-弹性连接假设，通过多元线性回归方法构建了调姿内力转化为位置补偿量的模型；根据Clamped-Free变形协调原理，简化了定位器调姿内力之间的协调关系，在此基础上提出了重力前馈补偿和调姿内力转化为位置补偿的力位协同控制策略，并对其进行了理论分析与设计。最后，对所提出的控制策略进行了仿真分析，结果表明采用力位协同控制方法，调姿定位精度提高35.3%，调姿内力降低77.8%，通过应用实验，说明了该方法的可行性和有效性。     

非圆截面变体头锥机构设计与运动仿真分析

设计具有变形能力的飞行器结构,通过主动改变其气动外形是提高空天飞行器飞行能力和环境自适应性的有效途径。本文提出了一种可实现单侧弯曲与多级伸展变形的空天飞行器非圆截面变体头锥机构设计方法,可独立进行头锥的多级伸展及弯曲运动,实现变体形式多样化;通过在伸展机构级间设计锁定装置,保证了头锥变形的展开稳定性。结合飞行工况,计算了头锥伸展过程中所受载荷及所需驱动力。头锥伸展及弯曲运动仿真结果论证了方案的可行性。在给定驱动模式下头锥机构的轴向最大伸展位移与体长比为0.43,伸展弯曲比可达3.5。通过分析变形过程的位移、速度及加速度变化规律发现,所设计的变体头锥机构可实现较好的运动特性和运动平稳性。     

基于目标加速度方向观测的微分对策制导律

针对延迟信息条件下的机动目标拦截问题,研究了一种考虑目标加速度方向观测的微分对策制导律。首先,针对引入目标方向观测信息末制导博弈问题,建立了微分对策数学模型。随后,利用状态依赖黎卡提方程求解微分对策问题,得到包含无延迟目标机动项的微分对策制导律,可以更为有效地拦截机动目标。最后,在延迟信息条件下,利用目标加速度方向观测的方法补偿延迟的目标加速度,再将补偿后的目标加速度信息取代无延迟的目标加速度,得到延迟信息条件下考虑目标加速度方向观测的微分对策制导律。仿真结果表明,新的制导律可以在延迟信息条件下有效地拦截机动目标。     

Combining MHD and kinetic modelling of solar flares

Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving thermal plasma (normally considered using magnetohydrodynamic (MHD) approach) and non-thermal plasma (requiring a kinetic approach). In this mini-review we consider different approaches aimed at bridging the gap between fluid and kinetic modelling of solar flares. Two types of approaches are discussed: combined MHD/test-particle (MHDTP) models, which can be used for modelling the flaring corona with relatively small numbers of energetic particles, and hybrid fluid-kinetic methods, which can be used for modelling stronger events with higher numbers of energetic particles. Two specific examples are discussed in more detail: MHDTP models of magnetic reconnection and particle acceleration in kink-unstable twisted coronal loops, and a novel reduced-kinetic model of particle transport in converging magnetic fields.     

Multi precursors analysis associated with the powerful Ecuador (MW = 7.8) earthquake of 16 April 2016 using Swarm satellites data in conjunction with other multi-platform satellite and ground data

After DEMETER satellite mission (2004–2010), the launch of the Swarm satellites (Alpha (A), Bravo (B) and Charlie (C)) has created a new opportunity in the study of earthquake ionospheric precursors. Nowadays, there is no doubt that multi precursors analysis is a necessary phase to better understand the LAIC (Lithosphere Atmosphere Ionosphere Coupling) mechanism before large earthquakes. In this study, using absolute scalar magnetometer, vector field magnetometer and electric field instrument on board Swarm satellites, GPS (Global Positioning System) measurements, MODIS-Aqua satellite and ECMWF (European Centre for Medium-Range Weather Forecasts) data, the variations of the electron density and temperature, magnetic field, TEC (Total Electron Content), LST (Land Surface Temperature), AOD (Aerosol Optical Depth) and SKT (SKin Temperature) have been surveyed to find the potential seismic anomalies around the strong Ecuador (M_w = 7.8) earthquake of 16 April 2016. The four solar and geomagnetic indices: F10.7, D_st, K_p and a_p were investigated to distinguish whether the preliminary detected anomalies might be associated with the solar-geomagnetic activities instead of the seismo-ionospheric anomalies. The Swarm satellites (A, B and C) data analysis indicate the anomalies in time series of electron density variations on 7, 11 and 12 days before the event; the unusual variations in time series of electron temperature on 8 days preceding the earthquake; the analysis of the magnetic field scalar and vectors data show the considerable anomalies 52, 48, 23, 16, 11, 9 and 7 days before the main shock. A striking anomaly is detected in TEC variations on 1 day before earthquake at 9:00 UTC. The analysis of MODIS-Aqua night-time images shows that LST increase unusually on 11 days prior to main shock. In addition, the AOD variations obtained from MODIS measurements reach the maximum value on 10 days before the earthquake. The SKT around epicentral region presents anomalous higher value about 40 days before the earthquake. It should be noted that the different lead times of the observed anomalies could be acknowledged based on a reasonable LAIC earthquake mechanism. Our results emphasize that the Swarm satellites measurements play an undeniable role in progress the studies of the ionospheric precursors.