Application of a modified generalized sail model on a solar sail with wrinkles

In this paper, we present an analysis of effect of wrinkles on the solar sail performance. We describe different analytical, semi-analytical and numerical approaches to the calculation of general large-scale curvature of a solar sail as well as parameters of so-called wrinkled domains, and introduce the impact of such wrinkles on the thrust and torque of the solar sail. Finally, we present a model of an optically-orthotropic surface for such non-ideal sail, providing a connection with the Generalized Sail Model, and other solar sail thrust models.     

空天飞行器返回滑翔段在线制导方法

针对空天往返飞行器的返回滑翔段在线制导问题，设计了一种新的滑翔段飞行剖面，实现了滑翔段终端交班高度、位置和倾角约束的自动满足，减少了在线制导算法中需处理的约束数量。推导了滑翔段运动状态、过程约束和性能指标的解析表达式，获得了剩余航程和终端速度间的函数关系。在此基础上，提出了一种双层在线制导方法：内层解析重构飞行剖面，同时通过解析确定路径点来改变剩余航程的变化率，进而对终端交班速度进行控制；外层借助解析表达式，使用粒子群优化算法(PSO)和改进共轭梯度法(CGM)优化飞行剖面，从而满足过程约束和指标要求。最后通过数学仿真验证了方法的正确性。     

多臂空间机器人的视觉伺服与协调控制

针对多臂空间机器人自主目标抓捕任务，首先建立多臂空间机器人的运动模型和其与目标的相对运动模型，采用Kane方法建立多臂空间机器人的动力学模型；其次，研究基于视觉伺服的机械臂在线轨迹规划算法，并引入零反作用机动，消除机械臂运动对平台姿态的扰动；再次，在不使用零反作用机动功能时，分别使用基于角动量前馈补偿的协调控制算法和逆动力学方法设计了协调控制器，在机械臂运动时保持平台姿态和相对目标的位置。最后，开发了基于Matlab的仿真软件MASS（多臂空间机器人仿真），仿真结果校验了上述方法的有效性。     

空间载荷机柜热控设计的仿真研究

空间站主动热控系统的作用是提供和保持航天员、仪器设备和结构部件所要求的温度、湿度和风速等热环境。本文以空间载荷机柜为热控系统的基本单元，建立了以单相流体循环回路为基础、采用液冷结合风冷的冷却方式的载荷机柜数学模型；介绍了模糊增量控制算法的工作原理；提出了采用模糊增量控制算法，以液冷支路进出口温差和风冷支路冷却量占总冷却量的比值为控制目标的反馈控制策略。同时文章采用数值模拟方法研究了控制系统的动态响应效果。仿真结果显示，模糊增量控制算法具有响应迅速、超调量小、无静态误差的特点，控制策略满足有效性、稳定性和高精确度的要求。     

Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

考虑测试性设计缺陷修正延时的测试性增长模型建模与评价方法

针对现有测试性增长模型忽略了测试性设计缺陷的修正延时过程，进而导致测试性增长模型（TGM）跟踪与预计精度低的问题，提出了一种考虑测试性设计缺陷修正延时的测试性增长模型建模理论与方法。首先分析测试性设计缺陷识别与修正之间的延时机理，得到剩余测试性设计缺陷（TDL）具有先增后减的铃形变化趋势；在此基础上，分别以Gamma、Raleigh和Delay-S 3种曲线拟合剩余测试性设计缺陷（RTDL）变化趋势，研究建立基于以上3种曲线考虑修正延时的测试性增长模型；最后，以某机载稳定跟踪平台测试性增长试验数据验证测试性增长模型拟合、跟踪及预计效能。研究结果表明：基于该测试性增长试验数据，Gamma曲线可以精确地拟合剩余测试性设计缺陷变化规律，测试性增长模型跟踪和预计精度可达到10^-2数量级。     

Geometric error analysis of an over-constrained parallel tracking mechanism using the screw theory

This paper deals with geometric error modeling and sensitivity analysis of an overconstrained parallel tracking mechanism. The main contribution is the consideration of overconstrained features that are usually ignored in previous research. The reciprocal property between a motion and a force is applied to tackle this problem in the framework of the screw theory. First of all, a nominal kinematic model of the parallel tracking mechanism is formulated. On this basis, the actual twist of the moving platform is computed through the superposition of the joint twist and geometric errors. The actuation and constrained wrenches of each limb are applied to exclude the joint displacement. After eliminating repeated errors brought by the multiplication of wrenches, a geometric error model of the parallel tracking mechanism is built. Furthermore,two sensitivity indices are defined to select essential geometric errors for future kinematic calibration. Finally, the geometric error model with minimum geometric errors is verified by simulation with SolidWorks software. Two typical poses of the parallel tracking mechanism are selected, and the differences between simulation and calculation results are very small. The results confirm the correctness and accuracy of the geometric error modeling method for over-constrained parallel mechanisms.     

基于Kriging模型梯度解析解的改进一次二阶矩方法

改进一次二阶矩（AFOSM）法是一种基于功能函数梯度的结构可靠性分析方法，鉴于其对隐式函数的梯度较难求解，提出了一种基于Kriging模型梯度解析解的AFOSM方法，利用Kriging代理模型的解析表达式推导得到功能函数对输入变量的梯度解析解，为AFOSM中设计点的确定提供高精度的梯度信息。通过Kriging与AFOSM的结合，很好地解决了基于有限元模型的隐式情况下梯度计算量相当大、可靠性分析难的问题。数值与工程算例验证了所提Kriging梯度解析解的较高精确性，同时验证了所提基于Kriging解析解的AFOSM结构可靠性分析方法的正确性与较高精度。     

基于自适应Kriging代理模型的交叉熵重要抽样法

对于复杂失效域和小失效概率耦合的可靠性分析问题，本文提出了一种交叉熵重要抽样（CE-IS）方法结合自适应Kriging （AK）代理模型的求解方法（CE-IS-AK）。所提方法基于交叉熵原理，用混合高斯模型逐步逼近最优重要抽样密度函数，并采用AK模型协助逼近过程中混合高斯模型的参数的更新，从而提高了CE-IS方法的计算效率。另外，本文还改进了CE-IS方法的收敛准则，避免了方法的冗余迭代，扩大了方法的适用范围。由于在CE-IS方法中引入了AK模型，因此，本文方法所构建的重要抽样函数在保证精度的基础上提高了效率。相较于AK-MCS方法，本文方法中引入了重要抽样的思想，因此在Kriging训练点数目基本相同的情况下，大幅缩减小失效概率计算时样本池规模，并且由于利用了混合高斯模型，因而对多失效域具有较好的适用性。算例分析也证明了本文所提方法的优越性。     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.