Variational discretization of constrained Birkhoffian systems

In this paper, we derive a variational characterization of constrained Birkhoffian dynamics in both continuous and discrete settings. When additional algebraic constraints appear, derivation of the necessary conditions under which the Pfaff action is extremized gives constrained Birkhoffian equations. Inspired by this continuous framework, we directly discretize the constraints as well as the Pfaff action and consequently formulate the discrete constrained Birkhoffian dynamics. Via this discrete variational approach which is parallel with the continuous case, the resulting discrete constrained Birkhoffian equations automatically preserve the intrinsic symplectic structure when identified as numerical algorithms. Considering that the obtained algorithms require not only the specification of an initial configuration but also a second configuration to operate, we present a natural, reasonable, and efficient method of initialization of simulations. While retaining the structure-preserving property, the obtained discrete schemes exhibit excellent numerical behaviors, demonstrated by numerical examples dealing with the mathematical pendulum and the 3D pendulum.     

Birkhoff力学的研究进展

Birkhoff力学是Hamilton力学的一个自然发展,是分析力学发展的一个新阶段,它广泛应用于力学、物理学和工程.本文总结Birkhoff力学的形成和发展,特别是近二十年所取得的成就.首先,从Birkhoff的《动力系统》中的有关段落开始,叙述Birkhoff力学的起源.其次,叙述这个力学的基本原理——Pfaff-Birkhoff原理以及这个力学的基本方程——Birkhoff方程的形成和发展.第三,简述Birkhoff力学的一些专门问题,包括约束Birkhoff系统,Birkhoff方程的积分方法,Birkhoff动力学逆问题,Birkhoff方程的运动稳定性,Birkhoff系统的几何方法,Birkhoff系统的全局分析等.最后,对Birkhoff力学的未来研究提出一些建议.      

Birkhoff 系统稳定性的动力学控制

本文研究 Birkhoff 系统和广义 Birkhoff 系统平衡稳定性的动力学控制. 首先建立系统的运动方程和平衡方程. 其次,研究 Birkhoff 系统中控制参数出现在 Birkhoff 函数中平衡稳 定性的动力学控制. 方法是通过选取控制参数使得 Birkhoff 函数 $B$ 成为定号函数,而其时间导数 $\dot {B}$ 为与 $B$ 反号的常号函数. 再次,研究广义 Birkhoff 系统平衡稳定性的动力学控制,通过选取 Birkhoff 函数或附加项中包含控制参数的方法,使得 Birkhoff 函数是定号函数,而其时间导数为反号的常号函数,从而控制系统的平衡稳定性. 最后举例说明结果的应用.      

Optimal flow control for Navier–Stokes equations: drag minimization

Optimal control and shape optimization techniques have an increasing role in Fluid Dynamics problems governed by partial differential equations (PDEs). In this paper, we consider the problem of drag minimization for a body in relative motion in a fluid by controlling the velocity through the body boundary. With this aim, we handle with an optimal control approach applied to the steady incompressible Navier–Stokes equations. We use the Lagrangian functional approach and we consider the Lagrangian multiplier method for the treatment of the Dirichlet boundary conditions, which include the control function itself. Moreover, we express the drag coefficient, which is the functional to be minimized, through the variational form of the Navier–Stokes equations. In this way, we can derive, in a straightforward manner, the adjoint and sensitivity equations associated with the optimal control problem, even in the presence of Dirichlet control functions. The problem is solved numerically by an iterative optimization procedure applied to state and adjoint PDEs which we approximate by the finite element method. Copyright © 2007 John Wiley & Sons, Ltd.     

DYNAMICAL CONTROL OF STABILITY FOR BIRKHOFFIAN SYSTEM$^{\bf 1)}$

本文研究 Birkhoff 系统和广义 Birkhoff 系统平衡稳定性的动力学控制. 首先建立系统的运动方程和平衡方程. 其次,研究 Birkhoff 系统中控制参数出现在 Birkhoff 函数中平衡稳 定性的动力学控制. 方法是通过选取控制参数使得 Birkhoff 函数 $B$ 成为定号函数,而其时间导数 $\dot {B}$ 为与 $B$ 反号的常号函数. 再次,研究广义 Birkhoff 系统平衡稳定性的动力学控制,通过选取 Birkhoff 函数或附加项中包含控制参数的方法,使得 Birkhoff 函数是定号函数,而其时间导数为反号的常号函数,从而控制系统的平衡稳定性. 最后举例说明结果的应用.     

Birkhoff系统的一类积分不变量的构造

分别建立了自由Birkhoff系统和约束Birkhoff系统的非等时变分方程,并且利用系统的Birkhoff方程及其非等时变分方程证明,可由第一积分直接构造该系统基于非等时变分的一类积分不变量。文中,举例说明结果的应用。     

动态规划求解空间双臂机器人非完整运动最优控制问题

通过自适应动态规划研究自由漂浮空间双臂机器人运动的最优控制问题．针对空间双臂机器人的非完整性运动,采用自适应动态规划(Adaptive Dynamic Programming, ADP)方法求解其最优控制问题．根据多体动力学理论,推导了载体位置、姿态均无控制条件下,双臂空间机器人满足的系统动量守恒关系的非完整约束方程,并将其转化为控制系统的状态方程,从而将双臂空间机器人的非完整运动规划问题转化为对非线性系统的控制问题．文中根据自适应动态规划网络结构,利用神经网络来近似性能指标函数,进而用龙格库塔法求解状态变量．并给出了适合该类问题的一种效用函数具体表达式,保证了空间双臂机器人到达期望位置后不再继续运动．实现了对空间双臂机器人非完整运动规划的最优控制．数值仿真实验验证了ADP对求解空间双臂机器人非完整运动规划最优控制问题的有效性．     

一种Birkhoff形式下结构动响应问题的保辛中点格式

结构动响应预测是结构设计的基础,是结构振动控制、载荷识别的前提。本文在辛体系下针对结构动响应问题,提出了一种Birkhoff形式下的保辛中点格式。首先引入状态变量,并基于摄动方法将结构动响应方程转化为线性自治Birkhoff方程的形式,进一步利用中心差分推导出线性自治Birkhoff方程的中点格式,其证明是保辛的。该格式不要求Birkhoff方程系数矩阵非奇异,因此适用于奇数维系统。两个不同数值算例的结果充分验证了本文方法的卓越性,也凸显了相对于传统算法在计算精确度和稳定性方面的明显优势。     

Noether symmetry method for Birkhoffian systems in terms of generalized fractional operators

Compared with the Hamiltonian mechanics and the Lagrangian mechanics, the Birkhoffian mechanics is more general. The Birkhoffian mechanics is discussed on the basis of the generalized fractional operators, which are proposed recently. Therefore, differential equations of motion within generalized fractional operators are established. Then, in order to find the solutions to the differential equations, Noether symmetry, conserved quantity, perturbation to Noether symmetry and adiabatic invariant are investigated. In the end, two applications are given to illustrate the methods and results.     

Generalized canonical transformation for second-order Birkhoffian systems on time scales

The theory of time scales, which unifies continuous and discrete analysis, provides a powerful mathematical tool for the study of complex dynamic systems. It enables us to understand more clearly the essential problems of continuous systems and discrete systems as well as other complex systems. In this paper, the theory of generalized canonical transformation for second-order Birkhoffian systems on time scales is proposed and studied, which extends the canonical transformation theory of Hamilton canonical equations. First, the condition of generalized canonical transformation for the second-order Birkhoffian system on time scales is established.Second, based on this condition, six basic forms of generalized canonical transformation for the second-order Birkhoffian system on time scales are given. Also, the relationships between new variables and old variables for each of these cases are derived. In the end, an example is given to show the application of the results.     

The method of variation on parameters for integration of a generalized Birkhoffian system

This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedure for solving the problem can be divided into two steps:the first step,a system of auxiliary equations is constructed and its general solution is given;the second step,the parameters are varied,and the solution of the problem is obtained by using the properties of generalized canonical transformation.The method of variation on parameters for the generalized Birkhoffian system is of universal significance,and we take a nonholonomic system and a nonconservative system as examples to illustrate the application of the results of this paper.     

A new Lie symmetrical method of finding conserved quantity for Birkhoffian systems

For a Birkhoffian system, a new Lie symmetrical method to find a conserved quantity is given. Based on the invariance of the equations of motion for the system under a general infinitesimal transformation of Lie groups, the Lie symmetrical determining equations are obtained. Then, several important relationships which reveal the interior properties of the Birkhoffian system are given. By using these relationships, a new Lie symmetrical conservation law for the Birkhoffian system is presented. The new conserved quantity is constructed in terms of infinitesimal generators of the Lie symmetry and the system itself without solving the structural equation which may be very difficult to solve. Furthermore, several deductions are given in the special infinitesimal transformations and the results are reduced to a Hamiltonian system. Finally, one example is given to illustrate the method and results of the application.     

NUMERICAL SOLUTIONS OF OPTIMAL CONTROL FOR THERMALLY CONVECTIVE FLOWS

We study the numerical solution of optimal control problems associated with two-dimensional viscous incompressible thermally convective flows. Although the techniques apply to more general settings, the presentation is confined to the objectives of minimizing the vorticity in the steady state case and tracking the velocity field in the non-stationary case with boundary temperature controls. In the steady state case we develop a systematic way to use the Lagrange multiplier rules to derive an optimality system of equations from which an optimal solution can be computed; finite element methods are used to find approximate solutions for the optimality system of equations. In the time-dependent case a piecewise-in-time optimal control approach is proposed and the fully discrete approximation algorithm for solving the piecewise optimal control problem is defined. Numerical results are presented for both the steady state and time-dependent optimal control problems. © 1997 John Wiley & Sons, Ltd.     

判定非自治Birkhoff系统稳定性的广义组合梯度方法

研究判定非自治Birkhoff系统稳定性的广义组合梯度方法．首先,给出非自治Birkhoff系统和非自治广义Birkhoff系统的运动微分方程;其次,给出一类将广义梯度系统和广义斜梯度系统组合而成的广义组合梯度系统,并讨论广义组合梯度系统的一些性质;最后,将非自治Birkhoff系统和非自治广义Birkhoff系统在一定条件下表示成广义组合梯度系统,并用广义组合梯度系统的性质研究了这两类Birkhoff系统的稳定性．举例说明结果的应用．     

Use of the energy integral in optimal control of the spacecraft spatial attitude

In this paper, we study the problem of spacecraft optimal rotation into the position with given attitude. The rotation takes a fixed time. The control program is optimized by minimizing the energy integral. We present the analytic solution of this problem and obtain formalized equations and computational expressions for constructing the optimal rotation program. We also present an example and the results of mathematical simulation of the spacecraft motion dynamics under optimal control, which demonstrate the practical realizability of our method for the spacecraft spatial attitude control.     

A shape optimisation method of a body located in adiabatic flows

The purpose of this study is to derive an optimal shape of a body located in adiabatic flow. In this study, we use the equation of motion, the equation of continuity and the pressure–density relation derived from the Poisson’s law as the governing equation. The formulation is based on an optimal control theory in which a performance function of fluid force is taken into consideration. The performance function should be minimised satisfying the governing equations. This problem can be solved without constraints by using the adjoint equation with adjoint variables corresponding to the state equation. The performance function is defined by the drag and lift forces acting on the body. The weighted gradient method is applied as a minimisation technique, the Galerkin finite element method is used as a spatial discretisation and the implicit scheme is used as a temporal discretisation to solve the state equations. The mixed interpolation, the bubble function for velocity and the linear function for density, is employed as the interpolation. The optimal shape is obtained for a body in adiabatic flows.     

一类多体系统非完整运动最优规划的拟牛顿算法

研究带有非完整约束的一类多体系统运动规划问题。多体系统中的非完整约束通常是由不可积的速度约束或不可积的守恒律引起。在系统动量和动量矩守恒情况下,动力学方程降阶为非完整形式约束方程,系统的控制问题可转化为无漂移系统的非完整运动规划问题。文中首先导出具有多体开链系统的非完整运动模型。利用最优控制理论和最优化技术,采用输入参数化的方法将连续的最优控制问题转化为离散的最优控制问题,提出一种非完整多体系统运动规划的拟牛顿算法。最后将该方法用于自由漂浮的空间三连杆机构,仿真结果验证了该方法的有效性。     

Based on wavelet analysis to optimal control of motion planning of space manipulator

ＩｎｔｒｏｄｕｃｔｉｏｎＲｅｃｅｎｔｌｙ ,ｍｏｒｅａｎｄｍｏｒｅｐｅｏｐｌｅｐａｙａｔｔｅｎｔｉｏｎｔｏｔｈｅｎｏｎｈｏｌｏｎｏｍｉｃｍｏｔｉｏｎｐｌａｎｎｉｎｇｏｆｓｐａｃｅｍａｎｉｐｕｌａｔｏｒ[1～4].Ｓｐａｃｅｍａｎｉｐｕｌａｔｏｒｉｓａｔｙｐｉｃａｌｍｕｌｔｉｂｏｄｙｓｙｓｔｅｍｗｉｔｈｏｕｔｒｏｏｔ.Ｗｈｅｎｔｈｅｍａｎｉｐｕｌａｔｏｒｉｓｏｐｅｒａｔｅｄ ,ａｎｄａｔｔｉｔｕｄｅｃｏｎｔｒｏｌｓｙｓｔｅｍｏｆｔｈｅｓｐａｃｅｖｅｈｉｃｌｅｉｓｓｗｉｔｃｈｅｄｏｆｆ,ｒｅｌａｔｉｖｅｍｏｔ…     

Birkhoff系统的广义正则变换及其基本形式

动力学方程的积分问题是分析力学研究的一个重要方面．由于求解一般的动力学方程往往会遇到很大困难，因此可利用变量变换，使方程变得容易求解．文章研究Birkhoff系统的广义正则变换．首先，建立Birkhoff系统的广义正则变换的充分必要条件；其次，基于该条件，给出Birkhoff系统的广义正则变换的六种基本形式，导出每一种情况下新旧变量之间的变换关系．作为特例，文中给出Hamilton方程的正则变换．文末，给出算例以说明结果的应用．     

基于Gauss伪谱法的3D刚体摆姿态最优控制

研究Gauss伪谱法求解3D刚体摆姿态最优控制问题．针对其最优姿态控制问题,既要满足由任意位置运动到平衡位置姿态运动规划问题,又要满足系统含有动力学约束的力学模型问题,提出基于四元数来描述3D刚体摆的数学模型,建立3D刚体摆姿态的动力学和运动学方程,为了解决3D刚体摆在平衡位置处的姿态最优控制问题,设计基于Gauss伪谱算法的最优姿态开环控制器,得到了3D刚体摆的姿态最优控制轨迹,得到满足的可行解,通过仿真实验验证了其开环解在平衡位置的控制姿态最优性．