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 perturbation series methodology for non-conservative linear Hamiltonian system with damping

In this paper,the symplectic perturbation series methodology of the non-conservative linear Hamiltonian system is presented for the structural dynamic response with damping.Firstly,the linear Hamiltonian system is briefly introduced and its conservation law is proved based on the properties of the exterior products.Then the symplectic perturbation series methodology is proposed to deal with the non-conservative linear Hamiltonian system and its conservation law is further proved.The structural dynamic response problem with eternal load and damping is transformed as the non-conservative linear Hamiltonian system and the symplectic difference schemes for the non-conservative linear Hamiltonian system are established.The applicability and validity of the proposed method are demonstrated by three engineering examples.The results demonstrate that the presented methodology is better than the traditional Runge–Kutta algorithm in the prediction of long-time structural dynamic response under the same time step.

VARIATIONS ON A THEME BY EULER

The oldest and simplest difference scheme is the explicit Euler method. Usually, it is not symplectic for general Hamiltonian systems. It is interesting to ask: Under what conditions of Hamiltonians, the explicit Euler method becomes symplectic? In this paper, we give the class of Hamiltonians for which systems the explicit Euler method is symplectic. In fact, in these cases, the explicit Euler method is really the phase how of the systems, therefore symplectic. Most of important Hamiltonian systems can be decomposed as the summation of these simple systems. Then composition of the Euler method acting on these systems yields a symplectic method, also explicit. These systems are called symplectically separable. Classical separable Hamiltonian systems are symplectically separable. Especially, we prove that any polynomial Hamiltonian is symplectically separable.