航天器运输包装箱仿真验证平台技术研究

针对数量有限的物理跑车试验无法满足减振与保温性能测试需求的问题，提出一套航天器运输包装箱动力学与热学仿真验证方法，包括：建立适用于包装箱系统的刚柔耦合多体动力学系统，通过结合线路条件测试生成的动力学系统外部激励，实现减振性能虚拟跑车测试；建立基于计算流体力学的包装箱热学模型，通过模拟自然对流和空调控制，实现包装箱保温性能虚拟跑车测试；基于C/S架构和导航式流程设计思想，建立航天器运输包装箱仿真验证平台，通过实际案例证明该平台仿真结果与实际跑车测试数据具有较高的一致性。     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

A progressive approach to predict shot peening process parameters for forming integral panel of Al7050-T7451

In this paper, a progressive approach to predict the multiple shot peening process parameters for complex integral panel is proposed. Firstly, the invariable parameters in the forming process including shot size, mass flow, peening distance and peening angle are determined according to the empirical and machine type. Then, the optimal value of air pressure for the whole shot peening is selected by the experimental data. Finally, the feeding speed for every shot peening path is predicted by regression equation. The integral panel part with thickness from 2 mm to 5 mm and curvature radius from 3200 mm to 16000 mm is taken as a research object, and four experiments are conducted. In order to design specimens for acquiring the forming data, one experiment is conducted to compare the curvature radius of the plate and stringer-structural specimens, which were peened along the middle of the two stringers. The most striking finding of this experiment is that the outer shape error range is below 3.9%, so the plate specimens can be used in predicting feeding speed of the integral panel. The second experiment is performed and results show that when the coverage reaches the limit of 80%, the minimum feeding speed is 50 mm/s. By this feeding speed, the forming curvature radius of the specimens with different thickness from the third experiment is measured and compared with the research object, and the optimal air pressure is 0.15 MPa. Then, the plate specimens with thickness from 2 mm to 5 mm are peened in the fourth experiment, and the measured curvature radius data are used to calculate the feeding speed of different shot peening path by regressive analysis method. The algorithm is validated by forming a test part and the average deviation is 0.496 mm. It is shown that the approach can realize the forming of the integral panel precisely.     

Development and performance evaluation of tropospheric scintillation model on Ku-band satellite link over Akure,Nigeria

Existing amplitude scintillation prediction models often perform less satisfactorily when deployed outside the regions where they were formulated. This necessitates the need to evaluate the performance of scintillation models developed in one region using data data from other regions while documenting their relative errors. Due to its variation with elevation angle, frequency, other link parameters and meteorological factors, we employed three years (January 2016 to December 2018) of concurrently measured satellite radio beacons and tropospheric weather parameters to develop a location-based amplitude scintillation prediction model over the Earth-space path of Akure (7.17^oN, 5.18^oE), South-western Nigeria. The satellite beacon measurement used Tektronix Y400 NetTek Analyzer at 1 s integration time while meteorological parameters, namely; temperature, pressure and relative humidity were measured using Davis Vantage Vue weather station at 1 min integration time. Comparative study of the model’s performance with nine (9) existing scintillation prediction models indicates that the best and worst performing models, in terms of root mean square error (RMSE), are the Statistical Temperature and Refractivity (STN) and direct physical and statistical prediction (DPSP) models with values 11.48 and 51.03 respectively. Also, worst month analysis indicates that April, with respective enhancement and fade values of 0.88 and 0.90 dB for 0.01% exceedance, is the overall worst calendar month for amplitude scintillation.     

Study of highly perturbed spacecraft formation dynamics via approximation

This paper explores methods for approximating and analyzing the dynamics of highly perturbed spacecraft formations with an emphasis on computationally efficient approaches. This facilitates on-board computation or rapid preliminary mission design analysis. Perturbed formation dynamics are often approximated as linear time-varying (LTV) systems, for which Floquet theory can be used to analyze the degree of system instability. Furthermore, the angular momentum of the relative orbital state can be computed with the approximate dynamics to provide additional insight. A general methodology is developed first and then applied to the problem of unstable formation dynamics in asteroid orbits. Here the dominant perturbative effects due to low-order gravitational harmonics and solar radiation pressure are modeled. Numerical simulations validate the approach and illustrate the approximation accuracy achieved.     

Fault mode probability factor based fault-tolerant control for dissimilar redundant actuation system

This paper presents a Fault Mode Probability Factor (FMPF) based Fault-Tolerant Control (FTC) strategy for multiple faults of Dissimilar Redundant Actuation System (DRAS) composed of Hydraulic Actuator (HA) and Electro-Hydrostatic Actuator (EHA). The long-term service and severe working conditions can result in multiple gradual faults which can ultimately degrade the system performance, resulting in the system model drift into the fault state characterized with parameter uncertainty. The paper proposes to address this problem by using the historical statistics of the multiple gradual faults and the proposed FMPF to amend the system model with parameter uncertainty. To balance the system model precision and computation time, a Moving Window (MW) method is used to determine the applied historical statistics. The FMPF based FTC strategy is developed for the amended system model where the system estimation and Linear Quadratic Regulator (LQR) are updated at the end of system sampling period. The simulations of DRAS system subjected to multiple faults have been performed and the results indicate the effectiveness of the proposed approach.     

Topology optimization in lightweight design of a 3D-printed flapping-wing micro aerial vehicle

Topology optimization is an effective method to obtain a lightweight structure that meets the requirements of structural strength. Whether the optimization results meet the actual needs mainly depends on the accuracy of the material properties and the boundary conditions, especially for a tiny Flapping-wing Micro Aerial Vehicle (FMAV) transmission system manufactured by 3D printing. In this paper, experimental and numerical computation efforts were undertaken to gain a reliable topology optimization method for the bottom of the transmission system. First, the constitutive behavior of the ultraviolet (UV) curable resin used in fabrication was evaluated. Second, a numerical computation model describing further verified via experiments. Topology optimization modeling considering nonlinear factors, e.g. contact, friction and collision, was presented, and the optimization results were verified by both dynamic simulation and experiments. Finally, detailed discussions on different load cases and constraints were presented to clarify their effect on the optimization. Our methods and results presented in this paper may shed light on the lightweight design of a FMAV.     

Accurate approximation of in-ecliptic trajectories for E-sail with constant pitch angle

Propellantless continuous-thrust propulsion systems, such as electric solar wind sails, may be successfully used for new space missions, especially those requiring high-energy orbit transfers. When the mass-to-thrust ratio is sufficiently large, the spacecraft trajectory is characterized by long flight times with a number of revolutions around the Sun. The corresponding mission analysis, especially when addressed within an optimal context, requires a significant amount of simulation effort. Analytical trajectories are therefore useful aids in a preliminary phase of mission design, even though exact solution are very difficult to obtain. The aim of this paper is to present an accurate, analytical, approximation of the spacecraft trajectory generated by an electric solar wind sail with a constant pitch angle, using the latest mathematical model of the thrust vector. Assuming a heliocentric circular parking orbit and a two-dimensional scenario, the simulation results show that the proposed equations are able to accurately describe the actual spacecraft trajectory for a long time interval when the propulsive acceleration magnitude is sufficiently small.