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 ...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.展开更多
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 te...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.展开更多
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 li...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.展开更多
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-f...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.展开更多
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 pa...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.展开更多
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...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.展开更多
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. Firs...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.展开更多
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 e...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.展开更多
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 ...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.展开更多
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 os...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.展开更多
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-...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,展开更多
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...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.展开更多
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 i...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.展开更多
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.M...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.展开更多
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 ...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.展开更多
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...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 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 rel...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.展开更多
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 eigen...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.展开更多
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 so...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.展开更多
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 Preissma...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.展开更多
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘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.
基金This research was supported by the National Natural Science Foundation of China (Nos. 41230210 and 41204074), the Science Foundation of the Education Department of Yunnan Province (No. 2013Z152), and Statoil Company (Contract No. 4502502663).
文摘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.
基金Project supported by the National Natural Science Foundation of China(Nos.12072266,12172297,11972287,and 12072262)the Open Foundation of the State Key Laboratory of Structural Analysis for Industrial Equipment of China(No.GZ22106)。
文摘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.
基金supported by the National Natural Science Foundation of China (10772014)
文摘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.
基金Project supported by the National Natural Science Foundation of China(Nos.91648101 and11672233)the Northwestern Polytechnical University(NPU)Foundation for Fundamental Research(No.3102017AX008)the National Training Program of Innovation and Entrepreneurship for Undergraduates(No.S201710699033)
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10374119 and 10674154), and The 0ne- Hundred-Talents Project of Chinese Academy of Science.Acknowledgments We gratefully acknowledge Professor Ding P Z and Professor Liu X S for their hospitality and help in symplectic algorithm.
文摘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.
文摘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.
基金the support from the National Natural Science Foundation of China(Grants 12022209,11972103,and 11825202)the Liaoning Revitalization Talents Program of China(Grant XLYC1807126)the Fundamental Research Funds for the Central Universities(Grant DUT21LAB124).
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 10971226, 91130013, and 11001270)the National Basic Research Program of China (Grant No. 2009CB723802)
文摘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.
文摘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.
基金the National Natural Science Foundation of China(Nos.10632030 and 10572119)Program for New Century Excellent Talents of Ministry of Education of China(No.NCET-04-0958)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment
文摘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,
基金supported by the National Natural Science Foun-dation of China (11172334)
文摘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.
基金Supported by the Natural Science Foundation of China under Grant No.0971226the 973 Project of China under Grant No.2009CB723802+1 种基金the Research Innovation Fund of Hunan Province under Grant No.CX2011B011the Innovation Fund of NUDT under Grant No.B110205
文摘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.
基金Supported by the National Natural Science Foundation of China under Grant No.10962004the Natural Science Foundation of Inner Mongolia under Grant No.2009BS0101+1 种基金the Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant No.20070126002the Cultivation of Innovative Talent of "211 Project"of Inner Mongolia University
文摘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.
基金supported by the National Natural Science Foundation of China (Nos. 10772147 and10632030)the Ph. D. Program Foundation of Ministry of Education of China (No. 20070699028)+2 种基金the Natural Science Foundation of Shaanxi Province of China (No. 2006A07)the Open Foundationof State Key Laboratory of Structural Analysis of Industrial Equipment (No. GZ0802)the Foundation for Fundamental Research of Northwestern Polytechnical University
文摘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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10971226,91130013,and 11001270)the National Basic Research Program of China(Grant No.2009CB723802)+1 种基金the Research Innovation Fund of Hunan Province,China (Grant No.CX2011B011)the Innovation Fund of National University of Defense Technology,China(Grant No.B120205)
文摘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.
基金supported by National Natural Science Foundation of China(Nos.11805203,11775222,11575185)the National Magnetic Confinement Fusion Energy Research Project of China(2015GB111003)the Key Research Program of Frontier Sciences CAS(QYZDB-SSW-SYS004)。
文摘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.
基金This research is funded by the grants from Dalian Project of Innovation Foundation of Science and Technology(No.2018J11CY005)Research Program of State Key Laboratory of Structural Analysis for Industrial Equipment(No.S18313).
文摘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.
基金Project (Nos. 19902014 and 10272024) supported by the NationalNatural Science Foundation of China
文摘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.
基金Project supported by the Program for Innovative Research Team in Science and Technology in Fujian Province University,China,the Quanzhou High Level Talents Support Plan,China(Grant No.2017ZT012)the Promotion Program for Young and Middle-Aged Teacher in Science and Technology Research of Huaqiao University,China(Grant No.ZQN-YX502)
文摘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.