Symplectic partitioned Runge-Kutta method based onthe eighth-order nearly analytic discrete operator and its wavefield simulations

We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.

New exact solutions for free vibrations of rectangular thin plates by symplectic dual method

The separation of variables is employed to solve Hamiltonian dual form of eigenvalue problem for transverse free vibrations of thin plates, and formulation of the natural mode in closed form is performed. The closed-form natural mode satisfies the governing equation of the eigenvalue problem of thin plate exactly and is applicable for any types of boundary conditions. With all combinations of simplysupported （S） and clamped （C） boundary conditions applied to the natural mode, the mode shapes are obtained uniquely and two eigenvalue equations are derived with respect to two spatial coordinates, with the aid of which the normal modes and frequencies are solved exactly. It was believed that the exact eigensolutions for cases SSCC, SCCC and CCCC were unable to be obtained, however, they are successfully found in this paper. Comparisons between the present results and the FEM results validate the present exact solutions, which can thus be taken as the benchmark for verifying different approximate approaches.

Highly accurate symplectic element based on two variational principles

For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process is simple and straightforward. In this paper, based on the seminal idea of the generalized mixed methods, a simple, stable, and highly accurate 8-node noncompatible symplectic element(NCSE8) was developed by the combination of the modified Hellinger-Reissner mixed variational principle and the minimum energy principle. To ensure the accuracy of in-plane stress results, a simultaneous equation approach was also suggested. Numerical experimentation shows that the accuracy of stress results of NCSE8 are nearly the same as that of displacement methods, and they are in good agreement with the exact solutions when the mesh is relatively fine. NCSE8 has advantages of the clearing concept, easy calculation by a finite element computer program, higher accuracy and wide applicability for various linear elasticity compressible and nearly incompressible material problems. It is possible that NCSE8 becomes even more advantageous for the fracture problems due to its better accuracy of stresses.

Symplectic system based analytical solution for bending of rectangular orthotropic plates on Winkler elastic foundation

This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish dual variables and dual equations in the symplectic space.The operator matrix of the equation set is proven to be a Hamilton operator matrix.Separation of variables and eigenfunction expansion creates a basis for analyzing the bending of rectangular orthotropic plates on Winkler elastic foundation and obtaining solutions for plates having any boundary condition.There is discussion of symplectic eigenvalue problems of orthotropic plates under two typical boundary conditions,with opposite sides simply supported and opposite sides clamped.Transcendental equations of eigenvalues and symplectic eigenvectors in analytical form given.Analytical solutions using two examples are presented to show the use of the new methods described in this paper.To verify the accuracy and convergence,a fully simply supported plate that is fully and simply supported under uniformly distributed load is used to compare the classical Navier method,the Levy method and the new method.Results show that the new technique has good accuracy and better convergence speed than other methods,especially in relation to internal forces.A fully clamped rectangular plate on Winkler foundation is solved to validate application of the new methods,with solutions compared to those produced by the Galerkin method.

Particle swarm optimization-based algorithm of a symplectic method for robotic dynamics and control

Multibody system dynamics provides a strong tool for the estimation of dynamic performances and the optimization of multisystem robot design. It can be described with differential algebraic equations(DAEs). In this paper, a particle swarm optimization(PSO) method is introduced to solve and control a symplectic multibody system for the first time. It is first combined with the symplectic method to solve problems in uncontrolled and controlled robotic arm systems. It is shown that the results conserve the energy and keep the constraints of the chaotic motion, which demonstrates the efficiency, accuracy, and time-saving ability of the method. To make the system move along the pre-planned path, which is a functional extremum problem, a double-PSO-based instantaneous optimal control is introduced. Examples are performed to test the effectiveness of the double-PSO-based instantaneous optimal control. The results show that the method has high accuracy, a fast convergence speed, and a wide range of applications.All the above verify the immense potential applications of the PSO method in multibody system dynamics.

Symplectic eigenfunction expansion theorem for elasticity of rectangular planes with two simply-supported opposite sides

The eigenvalue problem of the Hamiltonian operator associated with plane elasticity problems is investigated.The eigenfunctions of the operator are directly solved with mixed boundary conditions for the displacement and stress in a rectangular region.The completeness of the eigenfunctions is then proved,providing the feasibility of using separation of variables to solve the problems.A general solution is obtained with the symplectic eigenfunction expansion theorem.

Symplectic Schemes for Birkhoffian System

A universal symplectic structure for a Newtonian system including nonconservative cases can be constructed in the framework of Birkhoffian generalization of Hamiltonian mechanics. In this paper the symplectic geometry structure of Birkhoffian system is discussed, then the symplecticity of Birkhoffian phase flow is presented. Based on these properties we give a way to construct symplectic schemes for Birkhoffian systems by using the generating function method.

Symplectic analysis for elastic wave propagation in two-dimensional cellular structures

On the basis of the finite element analysis, the elastic wave propagation in cellular structures is investigated using the symplectic algorithm. The variation principle is first applied to obtain the dual variables and the wave propagation problem is then transformed into two-dimensional （2D） symplectic eigenvalue problems, where the extended Wittrick-Williams algorithm is used to ensure that no phase propagation eigenvalues are missed during computation. Three typical cellular structures, square, triangle and hexagon, are introduced to illustrate the unique feature of the symplectic algorithm in higher-frequency calculation, which is due to the conserved properties of the structure-preserving symplectic algorithm. On the basis of the dispersion relations and phase constant surface analysis, the band structure is shown to be insensitive to the material type at lower frequencies, however, much more related at higher frequencies. This paper also demonstrates how the boundary conditions adopted in the finite element modeling process and the structures＇ configurations affect the band structures. The hexagonal cells are demonstrated to be more efficient for sound insulation at higher frequencies, while the triangular cells are preferred at lower frequencies. No complete band gaps are observed for the square cells with fixed-end boundary conditions. The analysis of phase constant surfaces guides the design of 2D cellular structures where waves at certain frequencies do not propagate in specified directions. The findings from the present study will provide invaluable guidelines for the future application of cellular structures in sound insulation.

THE HAMILTONIAN SYSTEM AND COMPLETENESS OF SYMPLECTIC ORTHOGONAL SYSTEM

In this paper, a new Banach space ZH is defined, and it is proved that there is completeness of eigenfunction system (symplectic orthogonal system) of a class of Hamiltonian system in ZH space. We have also proved the following results: ZH space can be continuously imbedded to L-2[0,1] X L-2[0,1], but ZH not equal L-2[0,1] X L-1[0,1].

The symplectic eigenfunction expansion theorem and its application to the plate bending equation

This paper deals with the bending problem of rectangular plates with two opposite edges simply supported. It is proved that there exists no normed symplectic orthogonal eigenfunction system for the associated infinite-dimensional Hamiltonian operator H and that the two block operators belonging to Hamiltonian operator H possess two normed symplectic orthogonal eigenfunction systems in some space. It is demonstrated by using the properties of the block operators that the above bending problem can be solved by the symplectic eigenfunction expansion theorem, thereby obtaining analytical solutions of rectangular plates with two opposite edges simply supported and the other two edges supported in any manner.

SYMPLECTIC ALGORITHM IN SOLVING OPTIMAL CONTROL PROBLEMS

A symplectic algorithm is used to solve optimal control problems. Linear and nonlinear examples aregiven. Numerical analyses show that the symplectic algorithm gives satisfactory performance in that it works inlarge step and is of high speed and accuracy. This indicates that the symplectic algorithm is more effective andreasonable in solving optimal control problems.

SYMPLECTIC SOLUTION SYSTEM FOR REISSNER PLATE BENDING

Based on the Hellinger_Reissner variatonal principle for Reissner plate bending and introducing dual variables,Hamiltonian dual equations for Reissner plate bending were presented.Therefore Hamiltonian solution system can also be applied to Reissner plate bending problem,and the transformation from Euclidian space to symplectic space and from Lagrangian system to Hamiltonian system was realized.So in the symplectic space which consists of the original variables and their dual variables,the problem can be solved via effective mathematical physics methods such as the method of separation of variables and eigenfunction_vector expansion.All the eigensolutions and Jordan canonical form eigensolutions for zero eigenvalue of the Hamiltonian operator matrix are solved in detail,and their physical meanings are showed clearly.The adjoint symplectic orthonormal relation of the eigenfunction vectors for zero eigenvalue are formed.It is showed that the all eigensolutions for zero eigenvalue are basic solutions of the Saint_Venant problem and they form a perfect symplectic subspace for zero eigenvalue.And the eigensolutions for nonzero eigenvalue are covered by the Saint_Venant theorem.The symplectic solution method is not the same as the classical semi_inverse method and breaks through the limit of the traditional semi_inverse solution.The symplectic solution method will have vast application.

TWO NEW CONSTRUCTIONS OF CARTESIAN AUTHENTICATION CODES FROM SYMPLECTIC GEOMETRY

In this paper, two new constructions of Cartesian authentication codes from symplectic geometry are presented and their size parameters are computed.

Symplectic solution for three dimensional couple stress problem and its variational principle

A new state vector is presented for symplectic solution to three dimensional couple stress problem. Without relying on the analogy relationship, the dual PDEs of couple stress problem are derived by a new state vector. The duality solution methodology in a new form is thus extended to three dimensional couple stress. A new symplectic orthonormality relationship is proved. The symplectic solution to couple stress theory based a new state vector is more accordant with the custom of classical elasticity and is more convenient to process boundary conditions. A Hamilton mixed energy variational principle is derived by the integral method.

Pseudospectral method with symplectic algorithm for the solution of time-dependent SchrSdinger equations

A pseudospectral method with symplectic algorithm for the solution of time-dependent Schrodinger equations （TDSE） is introduced. The spatial part of the wavefunction is discretized into sparse grid by pseudospectral method and the time evolution is given in symplectic scheme. This method allows us to obtain a highly accurate and stable solution of TDSE. The effectiveness and efficiency of this method is demonstrated by the high-order harmonic spectra of one-dimensional atom in strong laser field as compared with previously published work. The influence of the additional static electric field is also investigated.

Difference Discrete Variational Principles, Euler?Lagrange Cohomology and Symplectic, Multisymplectic Structures III: Application to Symplectic and Multisymplectic Algorithms

In the previous papers I and II, we have studied the difference discrete variational principle and the Euler?Lagrange cohomology in the framework of multi-parameter differential approach. We have gotten the difference discrete Euler?Lagrange equations and canonical ones for the difference discrete versions of classical mechanics and field theory as well as the difference discrete versions for the Euler?Lagrange cohomology and applied them to get the necessary and sufficient condition for the symplectic or multisymplectic geometry preserving properties in both the Lagrangian and Hamiltonian formalisms. In this paper, we apply the difference discrete variational principle and Euler?Lagrange cohomological approach directly to the symplectic and multisymplectic algorithms. We will show that either Hamiltonian schemes or Lagrangian ones in both the symplectic and multisymplectic algorithms are variational integrators and their difference discrete symplectic structure-preserving properties can always be established not only in the solution space but also in the function space if and only if the related closed Euler?Lagrange cohomological conditions are satisfied.

Difference Discrete Variational Principle,Euler—Lagrange Cohomology and Symplectic,Multisymplectic Structures II：Euler—Lagrange Cohomology

In this second paper of a series of papers, we explore the difference discrete versions for the Euler?Lagrange cohomology and apply them to the symplectic or multisymplectic geometry and their preserving properties in both the Lagrangian and Hamiltonian formalisms for discrete mechanics and field theory in the framework of multi-parameter differential approach. In terms of the difference discrete Euler?Lagrange cohomological concepts, we show that the symplectic or multisymplectic geometry and their difference discrete structure-preserving properties can always be established not only in the solution spaces of the discrete Euler?Lagrange or canonical equations derived by the difference discrete variational principle but also in the function space in each case if and only if the relevant closed Euler?Lagrange cohomological conditions are satisfied.

Two New Fourth-Order Three-Stage Symplectic Integrators

Two new fourth-order three-stage symplectic integrators are specifically designed for a family of Hamiltonian systems,such as the harmonic oscillator,mathematical pendulum and latticeφ4 model.When the nonintegrable latticeφ4 system is taken as a test model,numerical comparisons show that the new methods have a great advantage over the second-order Verlet symplectic integrators in the accuracy of energy,become explicitly better than the usual non-gradient fourth-order seven-stage symplectic integrator of Forest and Ruth,and are almost equivalent to a fourth-order seven-stage force gradient symplectic integrator of Chin.As the most important advantage,the new integrators are convenient for solving the variational equations of many Hamiltonian systems so as to save a great deal of the computational cost when scanning a lot of orbits for chaos.

Total Energy Conservation and the Symplectic Algorithm

Based on the principle of total energy conservation, we give two important algorithms, the total energy conservation algorithm and the symplectic algorithm, which are established for the spherical shallow water equations. Also, the relation between the two algorithms is analyzed and numerical tests show the efficiency of the algorithms.

A New Multi-Symplectic Integration Method for the Nonlinear Schrödinger Equation

We propose a new multi-symplectic integration method for the nonlinear Schrödinger equation.The new scheme is derived by concatenating spatial discretization of the multi-symplectic Fourier pseudospectral method with temporal discretization of a symplectic Euler scheme and it is semi-explicit in the sense that it does not need to solve the nonlinear algebraic equations at every time step.We verify that the multi-symplectic semi-discretization of the Schrödinger equation with periodic boundary conditions has N semi-discrete multi-symplectic conservation laws.The discretization in time of the semi-discretization leads to N full-discrete multi-symplectic conservation laws.Numerical results are presented to demonstrate the robustness and the stability.