WKBJ近似保辛吗?

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

非线性方程的保辛近似算法

保守体系的微分方程可用Hamilton体系的方法描述,其特点是保辛。两个辛矩阵之和不能保辛,两个辛矩阵的乘积仍是辛矩阵。最常用的小参数摄动法用的是加法,因此对辛矩阵不能保辛。从保辛的角度,要用正则变换。本文针对非线性微分方程,运用自变量坐标变换,对原系统进行变换。由此推导出变换后系统的变分原理。引入Hamilton对偶变量,通过数学变换,得到变系数非线性方程。针对该方程,本文提出了保辛摄动算法。通过数值算例,对不同步长下,保辛摄动法、多尺度摄动法、龙格库塔法和精确解的结果做了比较。数值例题表明,对于非线性方程,本文提出的保辛摄动算法有良好的精度。在步长增大的情况下,保辛摄动保持了良好的稳定性。     

对一般Hamilton体系近似解保辛条件的讨论

文献[1]给出的一般Hamilton体系近似解保辛的条件,尚需讨论。其中(3.5~6)的证明适用于v·=H(z)v的系统,要求H(z)是Hamilton矩阵,齐次方程。即使Hamilton矩阵是与位移有关的H(z,q),仍可适用。但一般Hamilton体系未必能表示为v·=H·v,例如存在有势外力的情况,此时是非齐次线性方程。故在一般情况下,近似解是否保辛的原则尚需明确。一般的Hamilton正则方程体系非线性,通常用数值积分近似求解,故时间坐标离散而成为长η的序列,离散坐标体系。分析结构力学[2]考虑了离散坐标的情况,其中证明了区段两端状态向量之间关系ζ=ζ(v)的微商S(v)是…     

计算力学学报第25卷(2008年)总目次

第1期短波近似的保辛算法……………………………………………………………………钟万勰,孙雁(1-7)复杂机构的控制与结构同步优化设计……………………………………保宏,段宝岩,杜敬利,等(8-13)基于宏观调控策略的结构损伤检测PSO改进算法………………………余岭,万祖勇,朱宏平     

基于哈密顿动力系统新变分原理的保辛算法之一：变分原理和算法构造

提出了哈密顿动力系统的一个新变分原理,并基于此变分原理构造了四类保辛算法。通过新的变分原理定义修正作用量,然后将位移和动量采用拉格朗日多项式近似,并采用高斯积分对时间近似积分得到近似的修正作用量。在修正作用量的基础上,通过选择时间步两端不同的位移或动量作为独立变量,可构造四种不同类型的保辛算法。     

非线性问题的变分迭代方法及其应用

本文应用变分的概念，提出了求解非线性问题的加速迭代方法．这一方法的基本思想是：先给出问题的近似解，再引进一乘子校正其近似解；而乘子可用变分的概念最优确定，几个例子说明这种方法是有效的．     

动力学平衡方程的Euler中点辛差分求解格式

给出了动力学方程${\pmb M}\ddot {\pmb x} + {\pmb C}\dot {\pmb x} + {\pmb K \pmb x} = {\pmb R}$的二阶Euler中点隐式差分求解格式，分保守系统、无 阻尼受迫振动系统和阻尼系统3种情况, 讨论了算法中Jacobi矩阵${\pmb A}$的性质，譬 如${\pmb A}$是否为辛矩阵以及谱半径等. 对于无阻尼系统，证明了无论是否存在外 载荷，Jacobi 矩阵都是辛矩阵. 证明了辛矩阵的所有本征值的模为1，其谱半径永远 为1, 以及$\delta = 0.5$和$\alpha = 0.25$的Newmark算法就是Euler中点隐式差 分格式，对保守系统它们都是辛算法. 严格证 明了Euler中点辛格式是严格保持系统能量的. 通过算例详细讨论了保辛算法用于求解非保 守系统动态特性的优越性，如广义保结构特性等；分析了保辛算法的相位误差以及由其引起 的系统的附加能量特性；分析了保辛算法和$\delta \ne 0.5$的Newmark算法的精度随着激励频率与系统固有频率比的变化情况等     

基于哈密顿动力系统新变分原理的保辛算法之二：算法保辛性质证明

文献[1]给出了哈密顿系统的一个新的变分原理,并基于此变分原理,通过选择一个时间步长两端不同广义位移或广义动量为独立变量,给出了四种不同类型的求解哈密顿动力系统的数值方法。本文将分别证明这四类数值方法都是保辛的数值方法。     

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

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

基于哈密顿动力系统新变分原理的保辛算法之三：数值算例

文献[1,2]给出了四种不同类型的求解哈密顿动力系统的数值方法,并证明了它们的保辛特性。本文将讨论这四类算法的具体数值性能,包括算法的线性稳定性,精度和效率等。     

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 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.     

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

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

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

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

粘弹性厚壁筒问题的辛本征解方法

将哈密顿体系引进到粘弹性力学厚壁筒问题中,在辛体系下重新描述了基本问题,即建立了正则方程组。借助于积分变换,得到了拉伸、扭转和弯曲等问题的解以及有边界局部效应的解。将原问题归结为辛几何空间中的零本征值本征解和非零本征值本征解问题,从而建立了一种有效的分析问题方法和数值方法。为解决同类问题提供了一条可行的路径。     

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.