Plane elasticity in sectorial domain and the Hamiltonian system

ＰＬＡＮＥＥＬＡＳＴＩＣＩＴＹＩＮＳＥＣＴＯＲＩＡＬＤＯＭＡＩＮＡＮＤＴＨＥＨＡＭＩＬＴＯＮＩＡＮＳＹＳＴＥＭＺｈｏｎｇＷａｎ－ｘｉｅｔａｐ（钟万勰）（ＤａｌｉａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ．Ｄａｌｉａｎ）（ＲｅｃｅｉｖｅｄＤｅｃ．...     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper,a new analytical method of symplectic system.Hamiltonian system,is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain.In the system,the fundamental problem is reduced to all eigenvalue and eigensolution problem.The solution and boundary conditions call be expanded by eigensolutions using ad.ioint relationships of the symplectic ortho-normalization between the eigensolutions.A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space.The results show that fundamental flows can be described by zero eigenvalue eigensolutions,and local effects by nonzero eigenvalue eigensolutions.Numerical examples give various flows in a rectangular domain and show effectivenees of the method for solving a variety of problems.Meanwhile.the method can be used in solving other problems.     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper, a new analytical method of symplectic system, Hamiltonian system, is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain. In the system, the fundamental problem is reduced to an eigenvalue and eigensolution problem. The solution and boundary conditions can be expanded by eigensolutions using adjoint relationships of the symplectic ortho-normalization between the eigensolutions. A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space. The results show that fundamental flows can be described by zero eigenvalue eigensolutions, and local effects by nonzero eigenvalue eigensolutions. Numerical examples give various flows in a rectangular domain and show effectiveness of the method for solving a variety of problems. Meanwhile, the method can be used in solving other problems.     

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.     

Stokes流问题中的辛本征解方法

通过引入哈密顿体系,将二维Stokes流问题归结为哈密顿体系下的本 征值和本征解问题. 利用辛本征解空间的完备性,建立一套封闭的求解问题方法. 研究结果 表明零本征值本征解描述了基本的流动,而非零本征值本征解则显示着端部效应影响特点. 数值算例给出了辛本征值和本征解的一些规律和具体例子. 这些数值例子说明了端部非规则 流动的衰减规律. 为研究其它问题提供了一条路径.     

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

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

基于Hamilton体系的辛半解析法在各向异性电磁波导中的应用

力学中的Hamilton体系使用对偶变量来描述问题,而电磁场正好有电场和磁场这一对对偶变量。本文将力学中的Hamilton体系应用到电磁波导问题。根据电磁波导的Hamilton体系理论,辛几何可用于任意各向异性材料。将横向的电场和磁场构成对偶向量,基于Hamilton变分原理做半解析横向离散,并保持结构辛体系。本文以各向异性材料电磁波导为例,求解了问题的辛本征值,得到了镜像线的色散曲线。     

电磁波导辛体系下的周期皱波导分析

基于Hamilton变分原理的电磁波导辛体系自建立以来解决了传统电磁有限元所不能解决的一些问题。本文在介绍这一体系之后,经半解析横向离散及辛正则化,给出了类凝聚和协调质量阵。针对常见的周期皱波导问题,引入等效折射率概念,将皱波导转化为折射率周期变化的多层薄膜,并将其对应为力学分析中的条形域问题。最后给出的数值例子中所计算的通带辛本征值与解析解很接近,表明该理论方法有很高的计算效率。     

Hamiltonian formulation of anisotropic elasticity

The general expression of Hamiltonian partial differential equation for 3-D anisotropic elasticity is derived in this paper. It is a two-variable mixed formulation based on which new analytic and semi-analytic methods can be introduced.     

计算最优控制辛数值方法

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

弹性矩形板问题的Hamilton正则方程

为了采用辛算法求出弹性矩形板问题的解析解，中直接从弹性矩形板的控制方程出发推导了弹性矩形板，其中包括弹性矩形薄板和厚板问题以及弹性地基上矩形薄板和厚板问题的Hamilton正则方程，为利用辛几何方法求出任意边界条件下这类问题的理论解奠定了基础．     

Completeness of system of root vectors of upper triangular infinite-dimensional Hamiltonian operators appearing in elasticity theory

This paper deals with a class of upper triangular infinite-dimensional Hamiltonian operators appearing in the elasticity theory.The geometric multiplicity and algebraic index of the eigenvalue are investigated.Furthermore,the algebraic multiplicity of the eigenvalue is obtained.Based on these properties,the concrete completeness formulation of the system of eigenvectors or root vectors of the Hamiltonian operator is proposed.It is shown that the completeness is determined by the system of eigenvectors of the operator entries.Finally,the applications of the results to some problems in the elasticity theory are presented.     

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.     

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.     

辛方法和弹性圆柱壳在内外压和轴向冲击下的动态屈曲

针对有内压或外压的弹性圆柱壳在轴向冲击载荷耦合作用下的动态屈曲问题,构造哈密顿体系,在辛空间中将临界载荷和动态屈曲模态归结为辛本征值和本征解问题,从而形成一种辛方法。该方法直接得到非轴对称的屈曲模态。数值结果给出了圆柱壳问题的临界载荷和屈曲模态以及一些规律。     

EIGENVALUE PROBLEM OF A LARGE SCALE INDEFINITE GYROSCOPIC DYNAMIC SYSTEM

Gyroscopic dynamic system can be introduced to Hamiltonian system.Based on an adjoint symplectic subspace iteration method of Hamiltonian gyroscopic system, an adjoint symplectic subspace iteration method of indefinite Hamiltonian function gy- roscopic system was proposed to solve the eigenvalue problem of indefinite Hamiltonian function gyroscopic system.The character that the eigenvalues of Hamiltonian gyroscopic system are only pure imaginary or zero was used.The eigenvalues that Hamiltonian function is negative can be separated so that the eigenvalue problem of positive definite Hamiltonian function system was presented,and an adjoint symplectic subspace iteration method of positive definite Hamiltonian function system was used to solve the separated eigenvalue problem.Therefore,the eigenvalue problem of indefinite Hamiltonian function gyroscopic system was solved,and two numerical examples were given to demonstrate that the eigensolutions converge exactly.     

含裂纹纳米板振动问题的哈密顿体系方法

基于裂纹处范德华力效应,采用非局部弹性理论构造纳米板模型,并通过导入哈密顿体系建立含裂纹纳米板振动问题的对偶正则控制方程组。在全状态向量表示的哈密顿体系下,将含裂纹纳米板的固有频率和振型问题归结为广义辛本征值和本征解问题。利用哈密顿体系具有的辛共轭正交关系,得到问题解的级数解析表达式。结合边界条件,得到固有频率与辛本征值的代数方程关系式,进而直接给出固有频率的表达式。数值结果表明,非局部尺寸参数和裂纹长度对纳米板振动的各阶固有频率有直接的影响。对比表明,辛方法是准确且可靠的,可为工程应用提供依据。     

The Hamiltonian equations in some mathematics and physics problems

IntroductionHamiltoniansystemsareveryimportantinmathematics,physicsandmechanics.SpeciallyafterFengKang[1]firstproposedthesymplecticnumericalmethodforfinitedimensionalHamiltoniansystemin 1 980s,thestudyofHamiltoniansystemhasbeenfurtherimprovedinboththeoryandapplication .ItisthereforenaturaltoconsidertheproblemaboutwhichkindofpartialdifferentialequationscanberewrittenorrepresentedbyHamiltoniansystem .SuchkindofproblemhasnowbeenwarmlystudiedbymanymathematiciansandphysicianssuchasinRefs.[2 -5] …     

旋转运动中弹性梁耦合振动的辛方法

研究旋转梁结构的弹性耦合振动问题。通过引入对偶体系,建立了解决该类问题的辛方法。在辛体系中描述旋转梁纵向和横向耦合振动控制方程,即哈密顿正则方程。进一步求解得到结构的固有振动频率及相应的振动模态,发现固有振动频率随转动角速度先升后降以及模态之间的某种转化规律。