Automated design of gravity-assist trajectories to Mars and the outer planets

In this paper, a new approach to planetary mission design is described which automates the search for gravity-assist trajectories. This method finds all conic solutions given a range of launch dates, a range of launch energies and a set of target planets. The new design tool is applied to the problems of finding multiple encounter trajectories to the outer planets and Venus gravity-assist trajectories to Mars. The last four-planet grand tour opportunity (until the year 2153) is identified. It requires an Earth launch in 1996 and encounters Jupiter, Uranus, Neptune, and Pluto. Venus gravity-assist trajectories to Mars for the 30 year period 1995–2024 are examined. It is shown that in many cases these trajectories require less launch energy to reach Mars than direct ballistic trajectories.Assistant Professor, School of Aeronautics and AstronauticsGraduate Student, School of Aeronautics and Astronautics     

Orbit Determination of Binary Asteroids

In addition to the detection of an asteroid moon or a binary asteroid, the knowledge of the satellite’s true orbit is of high importance to derive fundamental physical parameters of the binary system such as its mass and to shed light on its possible formation history and dynamical evolution (prograde/retrograde orbit, large/small eccentricity or inclination, etc.). A new methodology for preliminary orbit determination of binary asteroids – and visual binaries in general – is proposed. It is based on Thiele–Innes method combined with a ‘trial and error’ Monte-Carlo technique. This method provides the full set of solutions (bundle of orbits, with the 7 orbital elements) even for a reduced number of observations. The mass is a direct by-product of this orbit determination, from which one can next infer the bulk-density and porosity. In addition to the bundle of orbits, the method provides the marginal probability densities of the foreseen parameters. Such error analysis – since it avoids linear approximation – can be of importance for the prediction of the satellite’s position in the plane-of-sky during future stellar occultations or subsequent observations, but also for the analysis of the orbit’s secular evolution. After briefly describing the method, we present the algorithm and its application to some practical cases, with particular emphasis on asteroids binaries and applications on orbital evolution.     

On the possibility of investigation of the magnetic field structure in subphotospheric layers of the Sun according to observations of sunspot torsional oscillations

A method of investigation of the magnetic field structure in subphotospheric layers of the Sun has been developed. The method is based on observations of the torisonal oscillations of single sunspots. Characteristics of the torsional oscillations have been obtained from observations of the longitudinal magnetic field and radial velocities of seven single sunspots in the photospheric line Fe I λ5253 Å. The parameters of the torsional oscillations and magnetic tubes in the deep layers have been determined. The radius of the cross section of a magnetic flux tube forming a sunspot is greatest near the Sun’s surface and is approximately equal to the radius of a sunspot umbra. Down to the deeper layers, it decreases quite quickly. The longitudinal electric current appearing in the magnetic tube changes direction. The typical time of the current changes is determined by the period of the torsional oscillations. The intensity of the longitudinal magnetic field in the tube increases with depth. The Alfven wave velocity averaged over the length of a magnetic tube is tens or hundreds of times less than this velocity in a sunspot umbra. It decreases with an increase in the period of oscillations. A decrease in the Alfven wave velocity leads to an increase in the twisting angle of magnetic field lines.     

月球重力场对绕月低轨卫星的长期影响

从解析形式出发,利用月球重力场模型JGL165P1,分析了月球重力场(带谐项)对绕月低轨卫星的长期影响。为了减少计算误差,保证计算精度,在分析解中使用循环公式来计算倾角函数。结果指出对于一个高度为100km的极月轨道卫星,冻结轨道存在的可能性不大,但是当轨道倾角在i=90°附近或者高度再高一些,则有可能存在冻结轨道;对于100km高的初始圆轨道,卫星在无控的情况下半年内将会坠落到月球表面,如果高度增加到200km,则不进行轨道控制也不会坠落到月面上。利用仿真软件GEODYN解算出来的结果证实了上述结论。     

The characteristics of motion of thin magnetic tubes

Starting from the equations of motion of a thin magnetic tube, the characteristic curves and velocities and compatibility relations are derived as basis for investigating its motion and for correctly formulating the problem of stationary solution. It is shown that the characteristic velocity of transverse waves is related to the Alfvén Mach number of the flow in the tube. When the flow velocity exceeds the critical value for the Kelvin-Helmholtz instability, transverse waves cease to exist.     

Homoclinic and heteroclinic orbits in the photogravitational restricted three-body problem

We study numerically the asymptotic homoclinic and heteroclinic orbits around the hyperbolic Lyapunov periodic orbits which emanate from Euler's critical points L ₁ and L ₂, in the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these Lyapunov orbits, are also presented. Poincaré surface of sections of these manifolds on appropriate planes and several homoclinic and heteroclinic orbits for the gravitational case as well as for varying radiation factor q ₁, are displayed. Homoclinic-homoclinic and homoclinic-heteroclinic-homoclinic chains which link the interior with the exterior Hill's regions, are illustrated. We adopt the Sun-Jupiter system and assume that only the larger primary radiates. It is found that for small deviations of its value from the gravitational case (q ₁ = 1), the radiation pressure exerts a significant impact on the Hill's regions and on these asymptotic orbits.     

GLONASS卫星位置的动力学改进算法

为了减弱由于广播星历拟合引起的误差,本文提出了一种基于动力学轨道改进的GLONASS卫星位置算法。该算法仅采用GLONASS广播星历信息就能够提高GLONASS卫星位置坐标的计算精度。算例表明,该方法能够有效减弱GLONASS卫星广播星历的拟合误差,卫星位置的精度损失从0.5m降到mm级;外推2-3h,精度损失也仅在cm级水平。     

内约束方钢管混凝土柱抗震性能研究

由于方钢管对混凝土约束作用较弱,地震作用下方钢管混凝土柱底部钢管易出现屈曲,因此本文提出一种新型内约束方钢管混凝土柱。基于ABAQUS有限元软件,本文采用合理的材料本构模型建立内约束方钢管混凝土柱三维实体精细有限元模型,该模型能准确反应钢管、混凝土以及拉筋之间的相互作用,又能反应拟静力作用下混凝土的塑性损伤和钢材的循环硬化规律。有限元结果与试验结果吻合良好。首先,在此基础上笔选出最佳内约束形式,对拉箍筋方钢管混凝土柱的抗震性能明显优于圆环箍筋;其次,提出在不同轴压比下内约束方钢管混凝土柱的焊接拉筋最佳布置长度和合理体积配箍率;再次,探讨不同参数对内约束方钢管混凝土柱滞回性能的影响,结果表明:提高截面含钢率和长细比能有效改善组合柱的极限承载力,而轴压比在一定范围内有利于能提高柱的承载力;最后,讨论了约束措施对内约束方钢管混凝土柱耗能性能的影响。     

高钢级绳索取心钻杆管材残余应力的初步研究

有观点认为管材的残余应力是绳索取心钻杆非正常损坏的诱因之一。据此采用环切法对不同热处理和矫直工艺的冷拔绳索取心管材制造过程中的残余应力进行了检测比对,认为冷拔管材制造环节产生的残余应力是存在的,也是可控的。只要制造工艺措施得当,可以将高钢级合金管材的残余应力控制在可接受范围内;合格管材出厂时存在较高残余应力造成的钻杆折断等事故较为罕见。应从管材制造、钻杆加工、现场使用等方面入手,努力降低钻杆断裂事故。     

高轨卫星轨道预报中神经网络模型优化设计

高轨卫星是我国卫星导航系统的重要组成部分。提升该类卫星的轨道预报精度有利于用户定位精度的提高。提出了一种改进高轨卫星轨道预报精度的新方法。该方法避开了精化动力学模型的困难,尝试从轨道预报误差的规律中寻找突破。利用神经网络作为建立预报模型的工具,将某历史时刻的轨道预报误差作为训练样本,利用训练好的神经网络模型补偿当前时刻的预报轨道以提高轨道预报精度。对影响神经网络模型补偿效果的各因素进行了详细分析,制定了适应于高轨卫星短期、中期和长期预报的神经网络最优模型。利用实测数据进行了试验分析,结果表明:预报8,15及30 d应选择的训练步长分别为10,20及25 min;轨道预报8~30 d时,训练噪声均选取0.01。神经网络模型有效地改进了高轨卫星的轨道预报精度,预报4~30 d,轨道精度提高幅度为34.67%~82.37%不等。