WKBJ近似保辛吗?

WKBJ短波近似是最常用的有效求解方法之一。保守体系的微分方程可用Hamilton体系的方法描述，其特点是保辛。保辛给出保守体系结构最重要的特性。但WKBJ短波近似却未曾考虑保辛的问题。本文给出验证近似解保辛的条件，并指出WKBJ近似难于保辛。然后给出正则变换的摄动保辛方法。数值例题展示了提出的保辛算法的有效性。     

短波近似的保辛算法

WKBJ短波近似是最常用的有效求解方法之一。保守体系的微分方程可用Hamilton体系的方法描述,其特点是保辛。保辛给出保守体系结构最重要的特性。但WKBJ短波近似却未曾考虑保辛的问题。WKBJ近似可用自变量坐标变换,然后再给出其保辛摄动。数值例题展示了本文变换保辛算法的有效性。     

Hamilton体系下介电弹性体圆形薄膜的动力学建模与辛求解

采用辛算法研究了Hamilton体系下介电弹性体圆形薄膜的动力学响应。首先,将该问题引入Hamilton对偶变量体系,借助Legendre变换,给出系统的广义动量和Hamilton函数,通过对Hamilton函数作用量的变分,得到Hamilton体系下的正则方程。其次,对于得到的正则方程给出了辛Runge-Kutta的计算格式。最后,采用二级四阶辛Runge-Kutta算法对动力学系统进行了数值求解,和四级四阶经典Runge-Kutta算法进行对比,结果表明,二级四阶辛Runge-Kutta算法具有保能量以及长时间数值稳定的优势,同时说明四级四阶经典Runge-Kutta算法对于步长依赖的局限性。     

Hamilton体系下介电弹性体圆形薄膜的动力学建模与辛求解

采用辛算法研究了Hamilton体系下介电弹性体圆形薄膜的动力学响应。首先，将该问题引入Hamilton对偶变量体系，借助Legendre变换，给出系统的广义动量和Hamilton函数，通过对Hamilton函数作用量的变分，得到Hamilton体系下的正则方程。其次，对于得到的正则方程给出了辛Runge-Kutta的计算格式。最后，采用二级四阶辛Runge-Kutta算法对动力学系统进行了数值求解，和四级四阶经典Runge-Kutta算法进行对比，结果表明，二级四阶辛Runge-Kutta算法具有保能量以及长时间数值稳定的优势，同时说明四级四阶经典Runge-Kutta算法对于步长依赖的局限性。     

计算最优控制辛数值方法

针对最优控制问题（OCP）的辛数值方法研究及应用进行综述。主要涉及内容包括,动力学系统为常微分方程描述的一般无约束、含不等式约束和状态时滞的最优控制问题,微分代数方程描述的一般无约束、含不等式约束和含切换系统的最优控制问题,以及闭环最优控制问题。从间接法和直接法两个求解框架出发,重点介绍本课题组在保辛算法方面的研究工作。在间接法框架下,首先基于生成函数和变分原理,将OCP保辛离散为非线性方程组,再数值求解方程组。在直接法框架下,将OCP保辛离散为有限维的非线性规划问题（NLP）,再数值求解。针对闭环最优控制问题,提出了保辛模型预测控制、滚动时域估计和瞬时最优控制算法。研究表明,保辛算法具有高精度和高效率的特点,在航空航天和机器人等领域有着广泛应用前景和价值。     

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

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

辛体系下含对称破缺因素动力学系统的近似守恒律

自冯康先生创立Hamilton系统辛几何算法以来,诸如辛结构和能量守恒等守恒律逐渐成为动力学系统数值分析方法有效性的检验标准之一。然而,诸如阻尼耗散、外部激励与控制和变参数等对称破缺因素是实际力学系统本质特征,影响着系统的对称性与守恒量。因此,本文在辛体系下讨论含有对称破缺因素的动力学系统的近似守恒律。针对有限维随机激励Hamilton系统,讨论其辛结构;针对无限维非保守动力学系统、无限维变参数动力学系统、Hamilton函数时空依赖的无限维动力学系统和无限维随机激励动力学系统,重点讨论了对称破缺因素对系统局部动量耗散的影响。上述结果为含有对称破缺因素的动力学系统的辛分析方法奠定数学基础。     

FPU方程的多尺度保辛摄动积分

FPU问题是一个经典非线性问题,其计算涉及多尺度分析。本文针对FPU问题,提出多尺度保辛摄动算法,该方法具有多尺度效应,可以按不同尺度显示计算结果,长时间计算保真,可以克服刚性问题,采用较大的积分步长,可以克服数值共振现象。数值算例显示了本文算法的有效性。     

非线性Schrdinger方程的保辛数值求解

首先将非线性Schrdinger方程化为Hamilton正则方程形式,而后建立Hamilton体系下的变分原理。再用有限元法离散空间坐标,同时对时间坐标进行精细积分,最后运用混合能变分原理,提出非线性Schrdinger方程保辛数值解法。这种解法在保辛的同时,可以让能量和质量在积分格点上亦全部达到守恒。数值算例验证了该方法的有效性。     

离散精细时程积分的自适应求积算法研究

基于Hamilton体系下的精细时程积分方法,通过对载荷项进行离散,应用中值法使载荷项在时间步长内为常值,从而将非齐次动力方程转化为齐次动力方程,避免了矩阵的求逆运算;基于积分区间逐次半分的思想实现了任意时间步长的自适应求积。数值算例结果表明:在同等时间步长的非齐次系统中,精细时程积分的最大误差为中心差分法的2.8%,为Newmark法的2.2%,最大求解误差仅为0.029%。这充分说明了本文的离散精细时程积分的自适应求积算法具有很好的收敛性。     

多体动力学的几何积分方法研究进展

动力系统的几何积分研究是近20年来工程计算领域非常活跃的方向. 多体动力学方程(微分方程, 微分代数方程)是一类典型的动力系统, 将其从Lagrange体系向Hamilton系统过渡, 目的在于从欧氏几何过渡到辛几何形态, 将对偶变量引入到力学研究中, 然后利用辛几何的数学框架对多体系统动力学方程进行数值计算, 可以预知多体动力学系统的一些定性信息, 并在数值离散时能保持这些定性性质特征, 尤其在表示关键的物理意义时需要强调保持这些几何性质. 简要介绍多体系统(无约束多刚体系统、完整约束多刚体系统和柔性多体系统)的Hamilton正则方程的建立和几何积分方法的构造, 着重介绍了在多体动力学计算中非常有应用前景的高阶辛算法(合成辛算法、分裂合成辛算法和辛精细积分法)、多辛算法, 以及广义Hamilton 系统与Lie 群积分方法等计算几何力学方法, 并对Lie 群积分的投影方法、流形局部坐标法等方法进行了阐述.     

Numerical solutions of linear quadratic control for time-varying systems via symplectic conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented. It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation. The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

Numerical solutions of linear quadratic control for time-varying systems via symplectic par conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented.It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation.The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

High order symplectic conservative perturbation method for time-varying Hamiltonian system

This paper presents a high order symplectic conservative perturbation method for linear time-varying Hamiltonian system.Firstly,the dynamic equation of Hamiltonian system is gradually changed into a high order perturbation equation,which is solved approximately by resolving the Hamiltonian coefficient matrix into a "major component" and a "high order small quantity" and using perturbation transformation technique,then the solution to the original equation of Hamiltonian system is determined through a series of inverse transform.Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes,the transfer matrix is a symplectic matrix;furthermore,the exponential matrices can be calculated accurately by the precise time integration method,so the method presented in this paper has fine accuracy,efficiency and stability.The examples show that the proposed method can also give good results even though a large time step is selected,and with the increase of the perturbation order,the perturbation solutions tend to exact solutions rapidly.     

NUMERICAL METHOD BASED ON HAMILTON SYSTEM AND SYMPLECTIC ALGORITHM TO DIFFERENTIAL GAMES

The resolution of differential games often concerns the difficult problem of two points border value (TPBV), then ascribe linear quadratic differential game to Hamilton system. To Hamilton system, the algorithm of symplectic geometry has the merits of being able to copy the dynamic structure of Hamilton system and keep the measure of phase plane. From the viewpoint of Hamilton system, the symplectic characters of linear quadratic differential game were probed; as a try, Symplectic-Runge-Kutta algorithm was presented for the resolution of infinite horizon linear quadratic differential game. An example of numerical calculation was given, and the result can illuminate the feasibility of this method. At the same time, it embodies the fine conservation characteristics of symplectic algorithm to system energy.     

二维非局部线弹性平面问题的辛分析

将二维非局部线弹性理论引入到Hamilton体系下,基于变分原理推导得出了二维线弹性理论的对偶方程和相应的边界条件.在分析验证对偶方程的准确性的基础上,该套方法被应用于二维弹性平面波问题的求解.将精细积分与扩展的W-W算法相结合在Hamilton体系下建立了求解平面Rayleigh波的数值算法.从推导到计算的保辛性确保了辛体系非局部理论与算法的准确性.通过对不同算例的数值计算,分析和对比了非局部理论方法与传统局部理论方法的差别,并进一步指出了该套算法的适用性和优势所在.     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.