基于混合Lie算子辛算法的不变流形计算

平动点附近周期轨道的不变流形因其在低能轨道转移中起着重要作用而受到广泛关注.在设计低能轨道过程中不变流形要实时进行能量匹配,但利用传统数值积分方法进行积分时能量会耗散.显式辛算法具有比隐式辛算法计算效率高的优势,但其要求Hamilton系统必须分成两个可积的部分,而旋转坐标系下的圆型限制性三体问题是不可分的,因而显式辛算法难以用于求解旋转坐标系下的圆型限制性三体问题.本文通过引入混合Lie算子,成功实现了带三阶导数项的力梯度辛算法对圆型限制性三体问题的求解,并将基于混合Lie算子的带三阶导数项的辛算法与Runge-Kutta78算法和Runge-Kutta45算法进行仿真对比,仿真结果表明基于混合Lie算子的含有三阶导数项的辛算法位置精度高、能量误差小且计算效率高.利用基于混合Lie算子的带三阶导数项的辛算法计算不变流形,可以实现低能轨道转移过程中轨道拼接点的能量精准匹配.     

日地系统L_2点Halo轨道绳系卫星编队动力学

研究平动点附近周期轨道上旋转多体绳系卫星编队系统的非线性耦合动力学问题。编队系统为各卫星质量接近的轮辐状结构,位于日地系统第二平动点附近,整个系统的旋转保持系绳处于张紧状态,建立Hill限制性三体问题的绳系编队系统动力学模型。针对处于留位阶段的典型对称三星编队,在位于较大Halo轨道上无控制力作用的情况下,进行母卫星轨道运动与系绳摆动耦合运动的动力学模拟,分析轨道方向、卫星质量比、系绳长度以及初始旋转速度对编队系统整体稳定性的影响。     

平动点Halo轨道航天器保持的多目标优化

针对日地系统平动点附近Halo轨道航天器保持任务,考虑航天器的能量消耗与轨道保持精度需求,采用多目标优化方法设计了改进时变控制器用于航天器Halo轨道保持任务。首先,基于圆形限制性三体模型推导了航天器的相对动力学方程并基于此设计了线性时变控制器。然后,采用多目标优化方法对时变控制器参数进行优化,得到满足航天器能量消耗与轨道保持精度之间平衡的Pareto最优解。最后,通过对考虑模型与环境干扰情况的数值模拟,结果表明多目标优化方法对平动点Halo轨道航天器保持任务达到低能能耗与高精度目标,具有一定的应用价值。     

卫星编队飞行中相对轨道的J2摄动分析

详细分析了$J_{2}$摄动对编队卫星相对轨道构形的影响. $J_{2}$ 摄动对相对轨道的影响分为相对轨道构形的漂移、相对轨道平面的章动和进动. 首先, 分析了相对轨道构形漂移速度、章动角速度和进动角速度的一阶近似表达式的数量级及其影 响因素. 其次,给出一个准则,来判断同一相对轨道的漂移和转动之间的关系. 最后, 利用该准则,分析了主、从星的轨道根数差对相对轨道的漂移和转动的影响.     

THE EVOLUTION CHARACTERIZATION OF LIBRATIN POINT ORBITS AND APPLICATION RESEARCH

三体问题中, 轨道的受力和运动规律非常复杂. 对于特定的任务, 如何选择轨道的初始解是一大难题.针对平面三体问题, 利 用近拱点庞加莱映射, 对平动点顺行轨道和逆行轨道的长期和短期演化性质进行分析.根据轨道的初始状态将其分为逃逸轨道和捕获轨道.对于逃逸轨道, 给出了同宿轨道和异宿轨道的设计方法, 并利用两级微分修正法消除了拼接点处的位置不连续问题.对于捕获轨道, 得到了几类典型的周期和准周期轨道.对逆行轨道的演化性质进行分析时发现, 逆行轨道通常为准周期轨道, 比顺行轨道更加稳定.利用近拱点庞加莱映射可以快速确定不同类型轨道对应的初始状态, 为特定任务需求下的轨道设计提供了一种快速而有效的选择方案.     

COLLINEAR LIBRATION POINTS WITH CONTINUOUS LOW THRUST

研究了圆型限制性三体问题的平动点在连续小推力作用下的具体位置和动力学特征的变化. 研究表明,随着小推力方向在空间中的变化,平动点的具体位置也会发生相应的变化,文章详细阐述了这些变化的特征.针对航天任务中应用较多的共线平动点L1和L2,研究了其附近运动的稳定性状态,给出了线性条件稳定解,并在此基础上,构造了条件稳定解的高阶形式,将其结果与数值积分轨道进行了比对,两者符合得很好. 最后,进一步研究了共线平动点附近周期轨道族的演化状态,由于连续小推力引入的非对称性,周期轨道族会发生分岔现象.     

一类二维非自治离散系统中的复杂簇发振荡结构1）

簇发振荡是多时间尺度系统复杂动力学行为的典型代表,簇发振荡的动力学机制与分类问题是簇发研究的重要问题之一,但当前学者们所揭示的簇发振荡的结构大多较为简单.研究以非自治离散Duffing系统为例,探讨具有复杂分岔结构的新型簇发振荡模式,并将其分为两大类,一类经由Fold分岔所诱发的对称式簇发,另一类经由延迟倍周期分岔所诱发的非对称式簇发.快子系统的分岔表现为典型的含有两个Fold分岔点的S形不动点曲线,其上、下稳定支可经由倍周期（即Flip）分岔通向混沌.当非自治项（即慢变量）穿越Fold分岔点时,系统的轨线可以向上、下稳定支的各种吸引子（例如,周期轨道和混沌）进行转迁,因此得到了经由Fold分岔所诱发的各种对称式簇发;而当非自治项无法穿越Fold分岔点,但可以穿越Flip分岔点时,系统产生了延迟Flip分岔现象.基于此,得到了经由延迟Flip分岔所诱发的各种非对称簇发.特别地,文中所报道的簇发振荡模式展现出复杂的反向Flip分岔结构.研究结果表明,这与非自治项缓慢地反向穿越快子系统的Flip分岔点有关.研究结果丰富了离散系统簇发的动力学机理和分类.     

非均匀水流中非线性波传播的数值模拟

以一种考虑波流相互作用的新型{Boussinesq}型方程为控制方程组, 采用五阶{Runge}-{Kutta}-{England}格式离散时间积分,采用七点 差分格式离散空间导数,并通过采用恰当的出流边界条件,从而建立了非均匀水流中非线性 波传播的数值模拟模型. 通过对均匀水流与水深水域内和潜堤地形上存在弱流或强流时波浪 传播的数值模拟,说明模型能有效地反映水流对波浪传播的影响.     

双小行星探测轨道动力学研究进展

双小行星系统由在万有引力作用下彼此环绕的两颗小行星组成, 对研究太阳系起源、行星系统演化和行星防御都具有重要的价值, 近年来成为行星科学和航天动力学研究的热门对象, 对双小行星系统的原位探测也即将迎来热潮. 双小行星系统的独特构型和附近的复杂动力学环境为探测器轨道动力学和任务设计带来了全新的挑战, 为应对这些挑战所进行的研究也推动了轨道动力学基础理论的发展. 本文对双小行星探测轨道动力学的研究进展进行综述, 首先介绍了双小行星研究和探测的背景及意义, 简要阐述了双小行星系统形成理论及其附近轨道动力学的研究概况. 其次, 介绍了双小行星系统不规则引力场和相互引力势的建模方法, 进而展示了双星的姿态轨道耦合动力学, 即完全二体问题, 包括双星相对运动的平衡构型和稳定性. 接着, 介绍了描述双星附近探测器轨道运动的限制性完全三体问题的动力学模型, 以及该模型下的平动点、平动点周期轨道、大范围周期轨道、转移轨道和轨道维持等方面的研究进展. 第四部分综述了环绕双小行星系统单颗星的受摄二体问题, 以轨道摄动理论和行星系统中受摄二体问题的研究现状为背景, 介绍了环绕双小行星系统主星的半解析轨道动力学建模与轨道稳定性分析. 之后, 介绍了目前面向探测任务需求和考虑实际约束的轨道动力学研究和轨道设计. 最后, 基于目前研究进展, 分析了面临的若干问题, 对未来双小行星探测轨道动力学及相关技术的发展进行了讨论和展望.      

一类二维非自治离散系统中的复杂簇发振荡结构

簇发振荡是多时间尺度系统复杂动力学行为的典型代表,簇发振荡的动力学机制与分类问题是簇发研究的重要问题之一,但当前学者们所揭示的簇发振荡的结构大多较为简单.研究以非自治离散Duffing系统为例,探讨具有复杂分岔结构的新型簇发振荡模式,并将其分为两大类,一类经由Fold分岔所诱发的对称式簇发,另一类经由延迟倍周期分岔所诱发的非对称式簇发.快子系统的分岔表现为典型的含有两个Fold分岔点的S形不动点曲线,其上、下稳定支可经由倍周期(即Flip)分岔通向混沌.当非自治项(即慢变量)穿越Fold分岔点时,系统的轨线可以向上、下稳定支的各种吸引子(例如,周期轨道和混沌)进行转迁,因此得到了经由Fold分岔所诱发的各种对称式簇发;而当非自治项无法穿越Fold分岔点,但可以穿越Flip分岔点时,系统产生了延迟Flip分岔现象.基于此,得到了经由延迟Flip分岔所诱发的各种非对称簇发.特别地,文中所报道的簇发振荡模式展现出复杂的反向Flip分岔结构.研究结果表明,这与非自治项缓慢地反向穿越快子系统的Flip分岔点有关.研究结果丰富了离散系统簇发的动力学机理和分类.     

Some basic solutions for nonlinear creep

The aim of the paper is to derive the exact analytical expressions for torsion and bending creep of rods that obey the Norton–Bailey, Prandtl–Garofalo and Naumenko–Altenbach–Gorash constitutive models. The common secondary creep constitutive model is the Norton–Bailey law which gives a power law relationship between creep rate and stress. The closed form solutions for fractional Norton–Bailey creep law are derived. The analytical formulas express the torque and bending moment as functions of the time for the period of relaxation. Other formulas express the twist rate and curvature as functions of the time for the duration of engineering creep experiment. The derived formulas are suitable for the practically important problems of machinery. Namely, the formulas are relevant for calculation of hereditary effects for helical, leaf and disk springs and twisted shafts.     

Convective heat transfer from two rotating circular cylinders in tandem arrangement by using lattice Boltzmann method

The numerical investigation of the two-dimensional laminar flow past two ro- tating circular cylinders in the tandem arrangement is conducted by the lattice Boltzmann method. The numerical strategy is used for dealing with curved and moving boundaries of the second-order accuracy for velocity and temperature fields. The effects of various rotational speed ratios and gap spacing are studied with the Reynolds number of 100 and the Prandtl number of 0.71. A varied range of rotational speed ratios are investigated for four different gap spacing, i.e., 3.0, 1.5, 0.7, and 0.2. The results show that, for the first cylinder, the lift and drag coefficients for large gap spacing are similar to those for a single cylinder; for the second cylinder, the lift coefficient descends with the increase in the angular velocity for all gap spacing, while the drag coefficient ascends except for the gap spacing of 3.0. The results of the averaged periodic Nusselt number on the surface of the cylinders show that, for small distances between the cylinders and low angular velocities, conduction is a dominant mechanism of heat transfer, but for large distances and high angular velocities, convection is the main mechanism of heat transfer.     

Analytical solutions for long waves on a circular island with combined topographies

We derived analytical solutions for long waves on a circular island with combined topographies. The inner region is a circular island on a flat bottom and the outer region is a cylindrical island on a shoal. The solution for the circular island on a flat bottom was developed previously by Jung et al. (2010) [13]. The solution for the cylindrical island on a shoal was developed in the present study using the methods of separation of variables, Taylor series expansion, and Frobenius series. The present solution for the cylindrical island on a shoal was verified by comparing with previous solutions which was developed for limited cases. Also, some cases were investigated with different values of wave period, shoal vertex level, and power of the radius. The solutions for two different types of real islands were obtained using the solution for the combined topographies. The present solution for the combined topographies can be applied more accurately to the real island than the previous solution for the single bottom topography.     

An efficient coupled layerwise theory for static analysis of piezoelectric sandwich beams

Summary An efficient one-dimensional model is developed for the statics of piezoelectric sandwich beams. Third-order zigzag approximation is used for axial displacement, and the potential is approximated as piecewise linear. The displacement field is expressed in terms of three primary displacement variables and the electric potential variables by satisfying the conditions of zero transverse shear stress at the top and bottom and its continuity at layer interfaces. The deflection field accounts for the piezoelectric transverse normal strain. The governing equations are derived using a variational principle. The present results agree very well with the exact solution for thin and thick highly inhomogeneous simply supported hybrid sandwich beams. The developed theory can accurately model open and closed circuit boundary conditions. The first author is grateful to DST, Government of India, for financial support for this work.     

Feedback minimization of first-passage failure of quasi non-integrable Hamiltonian systems

The non-linear stochastic optimal control of quasi non-integrable Hamiltonian systems for minimizing their first-passage failure is investigated. A controlled quasi non-integrable Hamiltonian system is reduced to an one-dimensional controlled diffusion process of averaged Hamiltonian by using the stochastic averaging method for quasi non-integrable Hamiltonian systems. The dynamical programming equations and their associated boundary and final time conditions for the problems of maximization of reliability and of maximization of mean first-passage time are formulated. The optimal control law is derived from the dynamical programming equations and the control constraints. The dynamical programming equations for maximum reliability problem and for maximum mean first-passage time problem are finalized and their relationships to the backward Kolmogorov equation for the reliability function and the Pontryagin equation for mean first-passage time, respectively, are pointed out. The boundary condition at zero Hamiltonian is discussed. Two examples are worked out to illustrate the application and effectiveness of the proposed procedure.     

网格自适应技术在复杂外形流场模拟中的应用

建立了一套适用于非结构混合网格自适应方法,针对激波和涡的不同特征采用不同加密探测器,各向异性加密棱柱单元并沿物面法向方向剖分所有棱柱层,各向异性剖分四面体单元,并保证四面体与棱柱交界面上网格协调。构造Hermit插值近似投影物面新加网格点和基于Laplacian光滑方法对空间网格进行优化。通过网格自适应加密,使用Roe格式计算高超声速球头绕流的红玉现象得到明显减轻。F16飞机含激波和脱体涡的流场自适应计算表明,网格加密集中在激波面和涡核附近区域,激波和涡计算更准确。     

A coupled zigzag theory for the dynamics of piezoelectric hybrid cross-ply plates

A recently developed coupled third-order zigzag theory for the statics of piezoelectric hybrid cross-ply plates is extended to dynamics. The theory combines a third-order zigzag approximation for the in-plane displacements and a sub-layerwise linear approximation for the electric potential, considering all components of the electric field. The nonuniform variation of the transverse displacement due to the piezoelectric field is accounted for. The conditions for the absence of shear traction at the top and bottom surfaces and continuity of transverse shear stresses in the presence of electromechanical loading are satisfied exactly, thereby reducing the number of displacement variables to five, which is the same as in a first- or third-order equivalent single-layer theory. The governing equations of motion are derived from the extended Hamilton's principle. The theory is assessed by comparing the Navier solutions for the free and forced harmonic vibration response of simply supported plates with the exact three-dimensional piezoelasticity solutions. Comparisons for hybrid test, composite and sandwich plates establish that the present theory is quite accurate for the dynamic response of moderately thick plates.     

任意挠曲输液管道瞬时变分原理

给出了弹性系统瞬时势能泛函的一般表达式，得到了其关于任意挠曲输液管道的具体表达式，根据瞬时最小势能原理建立了该结构的有限元动力方程，并首次推导出其离心力载荷公式，讨论了管道结构静平衡位置的发散流速的确定，最后，用简单的实验结果对理论分析作了验证。     

Long-wavelength dispersion of acoustic waves in transversely inhomogeneous anisotropic plates

Long-wavelength onset of the fundamental branches is described for a free anisotropic plate with arbitrary through-plate variation of material properties. Main attention is given to the flexural branch. Closed-form expressions for the leading-order dispersion coefficient of the velocity and displacement are derived for a generic case and exemplified for the various types of either continuous, or discrete, or periodic inhomogeneity and for the monoclinic symmetry. The relevance of the static averaging is examined in detail. The bounds for the slope of the flexural velocity branch are established. The upper fundamental branches are considered for the case when these are uncoupled inplane and shear horizontal ones.     

Analysis of shock wave reflection from fixed and moving boundaries using a stabilized particle method

In the present paper, the efficiency of an enhanced formulation of the stabilized corrective smoothed particle method （CSPM） for simulation of shock wave propagation and reflection from fixed and moving solid boundaries in compressible fluids is investigated. The Lagrangian nature and its accuracy for imposing the boundary conditions are the two main reasons for adoption of CSPM. The governing equations are further modified for imposition of moving solid boundary conditions. In addition to the traditional artificial viscosity, which can remove numerically induced abnormal jumps in the field values, a velocity field smoothing technique is introduced as an efficient method for stabilizing the solution. The method has been implemented for one- and two-dimensional shock wave propagation and reflection from fixed and moving boundaries and the results have been compared with other available solutions. The method has also been adopted for simulation of shock wave propagation and reflection from infinite and finite solid boundaries.