A symplectic finite element method based on Galerkin discretization for solving linear systems

We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.

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.

Symplectic analysis for regulating wave propagation in a one-dimensional nonlinear graded metamaterial

An analytical method,called the symplectic mathematical method,is proposed to study the wave propagation in a spring-mass chain with gradient arranged local resonators and nonlinear ground springs.Combined with the linearized perturbation approach,the symplectic transform matrix for a unit cell of the weakly nonlinear graded metamaterial is derived,which only relies on the state vector.The results of the dispersion relation obtained with the symplectic mathematical method agree well with those achieved by the Bloch theory.It is shown that wider and lower frequency bandgaps are formed when the hardening nonlinearity and incident wave intensity increase.Subsequently,the displacement response and transmission performance of nonlinear graded metamaterials with finite length are studied.The dual tunable effects of nonlinearity and gradation on the wave propagation are explored under different excitation frequencies.For small excitation frequencies,the gradient parameter plays a dominant role compared with the nonlinearity.The reason is that the gradient tuning aims at the gradient arrangement of local resonators,which is limited by the critical value of the local resonator mass.In contrast,for larger excitation frequencies,the hardening nonlinearity is dominant and will contribute to the formation of a new bandgap.

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.

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.

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.

THEORETIC SOLUTION OF RECTANGULAR THIN PLATE ON FOUNDATION WITH FOUR EDGES FREE BY SYMPLECTIC GEOMETRY METHOD

The theoretic solution for rectangular thin plate on foundation with four edges free is derived by symplectic geometry method. In the analysis proceeding, the elastic foundation is presented by the Winkler model. Firstly, the basic equations for elastic thin plate are transferred into Hamilton canonical equations. The symplectic geometry method is used to separate the whole variables and eigenvalues are obtained simultaneously. Finally, according to the method of eigen function expansion, the explicit solution for rectangular thin plate on foundation with the boundary conditions of four edges frees are developed. Since the basic elasticity equations of thin plate are only used and it is not need to select the deformation function arbitrarily. Therefore, the solution is theoretical and reasonable. In order to show the correction of formulations derived, a numerical example is given to demonstrate the accuracy and convergence of the current solution.

On the symplectic superposition method for free vibration of rectangular thin plates with mixed boundary constraints on an edge

A novel symplectic superposition method has been proposed and developed for plate and shell problems in recent years.The method has yielded many new analytic solutions due to its rigorousness.In this study,the first endeavor is made to further developed the symplectic superposition method for the free vibration of rectangular thin plates with mixed boundary constraints on an edge.The Hamiltonian system-based governing equation is first introduced such that the mathematical techniques in the symplectic space are applied.The solution procedure incorporates separation of variables,symplectic eigen solution and superposition.The analytic solution of an original problem is finally obtained by a set of equations via the equivalence to the superposition of some elaborated subproblems.The natural frequency and mode shape results for representative plates with both clamped and simply supported boundary constraints imposed on the same edge are reported for benchmark use.The present method can be extended to more challenging problems that cannot be solved by conventional analytic methods.

Explicit multi-symplectic method for the Zakharov-Kuznetsov equation

We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler method for temporal discretization. It is verified that the proposed method has corresponding discrete multi-symplectic conservation laws. Numerical simulations indicate that the proposed scheme is characterized by excellent conservation.

Numerical experiment of anharmonic oscillators by using the symplectic scheme-shooting method

Symplectic scheme-shooting method (SSSM) is applied to solve the energy eigenvalues of anharmonic oscillators characterized by the potentials V(x)=λx 4 and V(x)=(1/2)x 2+λx 2α with α=2,3,4 and doubly anharmonic oscillators characterized by the potentials V(x)=(1/2)x 2+λ 1x 4 +λ 2x 6, and a high order symplectic scheme tailored to the "time"-dependent Hamiltonian function is presented. The numerical results illustrate that the energy eigenvalues of anharmonic oscillators with the symplectic scheme-shooting method are in good agreement with the numerical accurate ones obtained from the non-perturbative method by using an appropriately scaled basis for the expansion of each eigenfunction; and the energy eigenvalues of doubly anharmonic oscillators with the sympolectic scheme-shooting method are in good agreement with the exact ones and are better than the results obtained from the four-term asymptotic series. Therefore, the symplectic scheme-shooting method, which is very simple and is easy to grasp, is a good numerical algorithm.

Multi-symplectic methods for membrane free vibration equation

In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws-a multi-symplectic conservation law （CLS）, a local energy conservation law （ECL） as well as a local momentum conservation law （MCL） --is constructed to discrete the PDEs that are derived from the membrane free vibration equation. The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior,

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

This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a ＂major compo- nent＂ and a ＂high order small quantity＂ and using perturba- tion transformation technique, then the solution to the orig- inal 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 sym- plectic 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, ef- ficiency 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 or- der, the perturbation solutions tend to exact solutions rapidly.

Multi-Symplectic Splitting Method for Two-Dimensional Nonlinear Schrodinger Equation

Using the idea of splitting numerical methods and the multi-symplectic methods,we propose a multisymplectic splitting(MSS) method to solve the two-dimensional nonlinear Schrodinger equation(2D-NLSE) in this paper.It is further shown that the method constructed in this way preserve the global symplecticity exactly.Numerical experiments for the plane wave solution and singular solution of the 2D-NLSE show the accuracy and effectiveness of the proposed method.

Double Symplectic Eigenfunction Expansion Method of Free Vibration of Rectangular Thin Plates

The free vibration problem of rectangular thin plates is rewritten as a new upper triangular matrixdifferential system.For the associated operator matrix,we find that the two diagonal block operators are Hamiltonian.Moreover,the existence and completeness of normed symplectic orthogonal eigenfunction systems of these two blockoperators are demonstrated.Based on the completeness,the general solution of the free vibration of rectangular thinplates is given by double symplectic eigenfunction expansion method.

Multi-symplectic Runge-Kutta methods for Landau-Ginzburg-Higgs equation

Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations （PDEs） derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.

Multi-symplectic wavelet splitting method for the strongly coupled Schrodinger system

We propose a multi-symplectic wavelet splitting equations. Based on its mu]ti-symplectic formulation, method to solve the strongly coupled nonlinear SchrSdinger the strongly coupled nonlinear SchrSdinger equations can be split into one linear multi-symplectic subsystem and one nonlinear infinite-dimensional Hamiltonian subsystem. For the linear subsystem, the multi-symplectic wavelet collocation method and the symplectic Euler method are employed in spatial and temporal discretization, respectively. For the nonlinear subsystem, the mid-point symplectic scheme is used. Numerical simulations show the effectiveness of the proposed method during long-time numerical calculation.

A relativistic canonical symplectic particlein-cell method for energetic plasma analysis

A relativistic canonical symplectic particle-in-cell(RCSPIC)method for simulating energetic plasma processes is established.By use of the Hamiltonian for the relativistic Vlasov-Maxwell system,we obtain a discrete relativistic canonical Hamiltonian dynamical system,based on which the RCSPIC method is constructed by applying the symplectic temporal discrete method.Through a 106-step numerical test,the RCSPIC method is proven to possess long-term energy stability.The ability to calculate energetic plasma processes is shown by simulations of the reflection processes of a high-energy laser(1?×?1020 W cm-2)on the plasma edge.

A Symplectic Method of Numerical Simulation on Local Buckling for Cylindrical Long Shells under Axial Pulse Loads

In this paper,the local buckling of cylindrical long shells is discussed under axial pulse loads in a Hamiltonian system.Using this system,critical loads and modes of buckling of shells are reduced to symplectic eigenvalues and eigensolutions respectively.By the symplectic method,the solution of the local buckling of shells can be employed to the expansion series of symplectic eigensolutions in this system.As a result,relationships between critical buckling loads and other factors,such as length of pulse load,thickness of shells and circumferential orders,have been achieved.At the same time,symmetric and unsymmetric buckling modes have been discuss.Moreover,numerical results show that modes of post-buckling of shells can be Bamboo node-type,bending type,concave type and so on.Research in this paper provides analytical supports for ultimate load prediction and buckling failure assessment of cylindrical long shells under local axial pulse loads.

A symplectic eigensolution method in transversely isotropic piezoelectric cylindrical media

This paper reports establishment of a symplectic system and introduces a 3D sub-symplectic structure for transversely isotropic piezoelectric media. A complete space of eigensolutions is obtained directly. Thus all solutions of the problem are re- duced to finding eigenvalues and eigensolutions, which include zero-eigenvalue solutions and all their Jordan normal form of the corresponding Hamiltonian matrix and non-zero-eigenvalue solutions. The classical solutions are described by zero-eigen- solutions and non-zero-eigensolutions show localized solutions. Numerical results show some rules of non-zero-eigenvalue and their eigensolutions.

Second order conformal multi-symplectic method for the damped Korteweg–de Vries equation

A conformal multi-symplectic method has been proposed for the damped Korteweg–de Vries(DKdV) equation, which is based on the conformal multi-symplectic structure. By using the Strang-splitting method and the Preissmann box scheme,we obtain a conformal multi-symplectic scheme for multi-symplectic partial differential equations(PDEs) with added dissipation. Applying it to the DKdV equation, we construct a conformal multi-symplectic algorithm for it, which is of second order accuracy in time. Numerical experiments demonstrate that the proposed method not only preserves the dissipation rate of mass exactly with periodic boundary conditions, but also has excellent long-time numerical behavior.