自由下落猫姿态最优控制的混合优化策略

研究猫自由下落时的转体运动对探索宇航员在太空失重状态下的空间运动规律具有重要的参考价值.针对猫自由下落时四肢先着地现象的姿态最优控制问题,提出一种由Gauss伪谱法求解可行解与直接打靶法求解最优精确解相结合的混合优化策略.首先,根据猫在自由下落过程中的角动量守恒原理,推导出简化后的两对称刚体系统的非完整姿态运动方程;然后基于Gauss伪谱法将此无漂移系统的姿态非完整运动规划问题离散为非线性规划问题,并在不考虑实际性能指标函数的条件下利用序列二次规划算法求解此非线性规划问题在较少节点时对应的控制变量可行解,再通过三次样条插值获取较多节点时的控制变量值;最后基于直接打靶法将插值得到的控制值作为序列二次规划算法的初始猜测值,从而求解得到最优的控制输入,再代入系统姿态运动方程,通过数值积分得到落体猫的转体姿态运动曲线.通过数值仿真,求解得到的姿态运动曲线是光滑的,能以较高的精度到达预定的目标姿态;最优控制输入也能满足预先设计的零边界控制要求以及最大控制要求;结果表明了此混合优化策略具有较强的鲁棒性和有效性.     

Chebyshev spectral variational integrator and applications

The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.     

基于对偶变量变分原理的完整约束系统保辛算法

基于对偶变量变分原理,选择积分区间两端位移为独立变量,构造了求解完整约束哈密顿动力系统的高阶保辛算法。首先,利用拉格朗日多项式对作用量中的位移、动量及拉格朗日乘子进行近似;然后,对作用量中不包含约束的积分项采用Gauss积分近似,对作用量中包含约束的积分项采用Lobatto积分近似,从而得到近似作用量;最后,在此近似作用量的基础上,利用对偶变量变分原理,将求解完整约束哈密顿动力系统问题转化为一组非线性方程组的求解。算法具有保辛性和高阶收敛性,能够在位移的插值点处高精度地满足完整约束。算法的收敛阶数及数值性质通过数值算例验证。     

Analytical solutions for dynamics of dual-spin spacecraft and gyrostat-satellites under magnetic attitude control in omega-regimes

The attitude dynamics of a dual-spin spacecraft (a gyrostat with one rotor) with magnetic actuators attitude control is considered in the constant external magnetic field at the presence of the spacecraft’s own magnetic dipole moment, which is created proportionally to the angular velocity components (this motion regime can be called as “the omega-regime” or “the omega-maneuver”). The research of the dual-spin spacecraft angular motion under the action of the magnetic restoring torque is fulfilled in the generalized formulation close to the classical mechanics’ task of the heavy body/gyrostat motion in the Lagrange top. Analytical exact solutions of differential equations of the motion are obtained for all parameters in terms of elliptic integrals and the Jacobi functions. New obtained analytical solutions can be classified as results developing the classical fundamental problem of the rigid body and gyrostat motion around the fixed point. The technical application of the omega-regime to the angular reorientation of the spacecraft longitudinal axis along the angular momentum vector is considered.     

不规则区域平面弹性问题的正则区域重心插值配点法

将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10^-14以上。数值算例验证了所建立方法的有效性和计算精度。     

航天器太阳阵伸展过程最优控制的遗传算法

本文讨论航天器太阳阵伸展中航天器姿态的最优控制问题。利用动量矩守恒原理导出带太阳阵航天器的动力方程,指出了系统的非」完整约束性质。将航天器太阳阵展开过程中主体姿态控制问题转化为非线性系统最优控制问题。厚优控制中引入遗传算法,替代传统的牛顿迭代方法。提出基于遗传算法的非完整运动规划最优控制算法。通过数值仿真表明,该方法对太阳阵伸展过程航天器姿态控制是有效的。     

大型圆形贮液池力学分析中的重心插值拟谱法

针对传统大型圆形贮液池力学分析方法计算复杂、精度不高的问题,用拟谱法对圆形贮液池的线弹性静力问题进行了研究.采用重心Lagrange插值多项式建立未知函数的微分矩阵,将池壁的控制方程表示为代数方程组.通过求解代数方程组,求得池壁各个离散点挠度,进而采用微分矩阵直接求得池壁内力.算例表明,该方法原理简单,易于程序实现和数值计算精度高.     

拟谱-微分求积混合方法求解一类双曲电报方程

拟谱方法和微分求积法是两类重要的无网格法,二者都已在科学和工程计算中获得了广泛应用。采用拉格朗日插值多项式作为二者的试函数,且采用同一种网格点分布,指出了在空间域上,微分求积法是拟谱方法的一种特殊形式。在此基础上,结合二者各自的特点,提出了拟谱-微分求积混合方法用于求解一类双曲电报方程。理论分析和数值测试表明,新方法在空间域上具有谱精度收敛性,在时间域上是A-稳定的,比较适合于求解多维电报方程。     

最优轨迹规划方法及其平动点航天器编队均衡耗能重构应用

由于均衡耗能航天器编队能够提高整体航天器编队服役时间,针对平动点航天器编队重构的均衡耗能最优轨迹规划问题,提出一种以状态、协态和控制三类变量插值为核心的求解非线性最优控制问题的新方法。基于连续时间表达的非线性最优控制问题通过变分原理转化为非线性方程组的求解,并进一步推导非线性方程组显示格式的Jacobi矩阵提高非线性方程组的计算效率。本文方法既满足最优控制理论的一阶必要条件又具有较大的收敛域;同时,不需要对协态初值准确猜测,避免了大规模非线性规划问题的求解。通过对中心航天器位置固定和无中心航天器两种情况的数值模拟,结果表明,本文方法对航天器编队重构轨迹规划问题能够达到均衡耗能的目标,具有一定的应用价值。     

非线性MEMS微梁的重心有理插值迭代配点法分析

文章利用重心有理插值迭代配点法分析计算非线性MEMS微梁问题。通过处理MEMS微梁的几何通过假设初始函数,将微梁非线性控制方程转换为线性化微分方程,建立逼近非线性微分方程的线性化迭代格式。采用重心有理插值配点法求解线性化微分方程,提出了数值分析MEMS微梁非线性弯曲问题的重心插值迭代配点法。给出了非线性微分方程的直接线性化和Newton线性化计算公式,详细讨论了非线性积分项的计算方法和公式。利用重心有理插值微分矩阵,建立了矩阵-向量化的重心插值迭代配点法的计算公式。数值算例结果表明,重心插值迭代配点法求解微梁非线性弯曲问题,具有计算公式简单、程序实施方便和计算精度高的特点。     

多储液腔航天器刚液耦合动力学与复合控制

采用复合控制方法对充液航天器的姿态和轨道机动进行高精度控制.通过傅里叶-贝塞尔级数展开法,将低重力环境下液体的弯曲自由表面的动态边界条件转化为简单的微分方程,其中耦合液体晃动方程的状态向量由相对势函数的模态坐标和波高的模态坐标组成.通过广义准坐标下的拉格朗日方程得到航天器刚体部分运动和液体燃料晃动的耦合动力学方程,提出了自适应快速终端滑模策略和输入整形技术相结合的复合控制器,并分别用于控制携带有一个燃料腔和四个燃料腔航天器的轨道机动和姿态机动.通过数值模拟来验证控制器的效率和精度.结果表明,对于多储液腔航天器,如果在设计航天器的姿态和轨道控制器时没有充分考虑燃料晃动效应,那么在受控航天器系统中将会出现刚-液-控耦合问题并导致航天器姿态不稳定.而本研究中的复合自适应终端滑模控制器可以实现航天器机动的高精度控制并有效抑制液体燃料晃动.     

New way to construct high order Hamiltonian variational integrators

This paper develops a new approach to construct variational integrators. A simplified unconventional Hamilton’s variational principle corresponding to initial value problems is proposed, which is convenient for applications. The displacement and momentum are approximated with the same Lagrange interpolation. After the numerical integration and variational operation, the original problems are expressed as algebraic equations with the displacement and momentum at the interpolation points as unknown variables. Some particular variational integrators are derived. An optimal scheme of choosing initial values for the Newton-Raphson method is presented for the nonlinear dynamic system. In addition, specific examples show that the proposed integrators are symplectic when the interpolation point coincides with the numerical integration point, and both are Gaussian quadrature points. Meanwhile, compared with the same order symplectic Runge-Kutta methods, although the accuracy of the two methods is almost the same, the proposed integrators are much simpler and less computationally expensive.     

非自旋航天器混沌姿态运动及其参数开闭环控制

研究万有引力场中受大气阻力且存在结构内阻尼的非自旋航天器在椭圆轨道上平面天平动的混沌及其参数开闭环控制问题．在建立数学模型的基础上确定出现混沌的必要条件并数值验证混沌的存在性，提出非线性振动系统混沌运动的参数开闭环控制并应用于控制航天器的混沌姿态运动．     

燃料大幅晃动等几何分析仿真及航天器耦合动力学研究

现代航天器肩负许多周期长且复杂的航天任务,通常需要携带大量的液体燃料.贮箱中液体燃料大幅晃动会严重影响航天器的姿态稳定性和控制精度,是现代航天器耦合动力学建模和精确控制研究的重要问题.本文提出了一种新的液体大幅晃动数值仿真方法,采用等几何分析方法对贮箱内气体和液体整体进行建模和空间离散,采用压力修正的分步法对控制方程进行时间离散,结合水平集方法划分气体和液体区域并且实时追踪液体晃动自由面.提出了一种质量修正方法以消除水平集函数演化产生的液体质量误差.基于燃料大幅晃动等几何分析仿真方法,对携带太阳能帆板的充液航天器进行动力学建模和耦合运动数值仿真.对于太阳能帆板的振动问题则采用Kirchhoff-Love板理论建模和模态分析法数值求解.通过将数值仿真结果与解析解对比,证明了本文给出方法的正确性.本文还对燃料大幅晃动下的航天器刚-液-柔耦合运动进行了数值仿真,发现液体晃动对航天器的姿态变化和结构振动的幅值和频率具有不可忽视的影响.     

A mixture differential quadrature method for solving two-dimensional incompressible Navier-Stokes equations

Ｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ（ＤＱＭ）ｐｒｏｐｏｓｅｄｂｙＲ．Ｂｅｌｌｍａｎ［１,２］ｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｅｍｐｌｏｙｅｄｉｎｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｓｏｆｐｒｏｂｌｅｍｓｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｈｙｓｉｃａｌｓｃｉｅｎｃｅ．ＢｅｃａｕｓｅｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｎａｌｌｇｒｉｄｐｏｉｎｔｓｉｓｕｓｅｄｔｏｆｉｔｔｈｅｄｅｒｉｖａｔｉｖｅｓａｔｇｒｉｄｐｏｉｎｔｓｉｎｔｈｅＤＱＭ,ｉｔｉｓｅｎｏｕｇｈｔｏｏｂｔａ…     

带空间机械臂航天器姿态运动规划的数值算法研究

自由漂浮空间机械臂系统在无外力矩作用时,系统的动量矩守恒而成为非完整系统。利用这一特性本文研究了自由漂浮空间机械臂系统的三维姿态运动规划问题。导出带空间机械臂的航天器三维姿态运动数学模型,将系统的非完整运动规划问题转化为非线性系统最优控制问题,在最优控制中利用小波逼近控制输入规律,提出基于遗传算法的最优控制数值算法。通过数值仿真,表明该方法对带空间机械臂航天器系统的非完整姿态运动规划是有效的。     

HYBRID OPTIMIZATION STRATEGY FOR MOTION PLANNING OF SPACE ROBOT SYSTEM WITH PRISMATIC JOINT

研究了自由漂浮带滑移铰空间机器人非完整运动规划的最优控制问题,提出一种由高斯伪谱法求解可行解与直接打靶法求解最优解相结合的混合优化策略.首先,根据多体系统动力学理论建立空间机器人的动力学模型,给定系统的初始和目标位形,将空间机器人运动规划问题描述成博尔察（Bolza）型最优控制问题;然后,利用高斯伪谱法将最优控制问题离散为非线性规划问题,求解在较少勒让德-高斯（Legendre-Gauss,LG）点时状态变量和控制变量对应的可行解;最后,在LG点处离散控制变量,作为直接打靶法的初值,利用序列二次规划算法求解空间机器人系统的优化运动轨迹和最优控制输入.通过数值仿真,系统优化运动轨迹光滑平稳,最优控制输入也能很好地满足各种约束条件,仿真结果验证了该混合优化策略的鲁棒性和有效性.     

A numerical investigation of a spectral‐type nodal collocation discontinuous Galerkin approximation of the Euler and Navier–Stokes equations

In this paper, we focus on the applicability of spectral‐type collocation discontinuous Galerkin methods to the steady state numerical solution of the inviscid and viscous Navier–Stokes equations on meshes consisting of curved quadrilateral elements. The solution is approximated with piecewise Lagrange polynomials based on both Legendre–Gauss and Legendre–Gauss–Lobatto interpolation nodes. For the sake of computational efficiency, the interpolation nodes can be used also as quadrature points. In this case, however, the effect of the nonlinearities in the equations and/or curved elements leads to aliasing and/or commutation errors that may result in inaccurate or unstable computations. By a thorough numerical testing on a set of well known test cases available in the literature, it is here shown that the two sets of nodes behave very differently, with a clear advantage of the Legendre–Gauss nodes, which always displayed an accurate and robust behaviour in all the test cases considered.Copyright © 2012 John Wiley & Sons, Ltd.     

勒让德伪谱法求解三维刚体摆姿态运动最优控制问题

研究伪谱法求解三维刚体摆姿态运动最优控制问题. 针对三维刚体摆这类含有约束的力学模型,提出了基于勒让德伪谱法的三维刚体摆姿态最优控制方法. 利用插值逼近设计了三维刚体摆姿态运动最优控制算法,得到了三维刚体摆的姿态最优控制轨迹,并结合松弛参数来控制插值点的取值,寻找满足的可行解. 仿真结果表明,基于勒让德伪谱法的最优控制算法使得三维刚体摆能以较小的误差运动到期望的末端姿态,且计算速度快,能够获得精度较高的控制输入量.     

Chaotic attitude and reorientation maneuver for completely liquid-filled spacecraft with flexible appendage

The present paper investigates the chaotic attitude dynamics and reorientation maneuver for completely viscous liquid-filled spacecraft with flexible appendage. All of the equations of motion are derived by using Lagrangian mechanics and then transformed into a form consisting of an unperturbed part plus perturbed terms so that the system＇s nonlinear characteristics can be exploited in phase space. Emphases are laid on the chaotic attitude dynamics produced from certain sets of physical parameter values of the spacecraft when energy dissipation acts to derive the body from minor to major axis spin. Numerical solutions of these equations show that the attitude dynamics of liquid-filled flexible spacecraft possesses characteristics common to random, non- periodic solutions and chaos, and it is demonstrated that the desired reorientation maneuver is guaranteed by using a pair of thruster impulses. The control strategy for reorientation maneuver is designed and the numerical simulation results are presented for both the uncontrolled and controlled spins transition.