A Hamiltonian particle method for non-linear elastodynamics

Yukihito Suzuki*, Seiichi Koshizuka


研究成果: Article査読

34 被引用数 (Scopus)


Particle methods are meshless simulation techniques in which motion of continua is approximated by discrete dynamics of a finite number of particles. They have a great degree of flexibility, for instance, in dealing with complex large deformations or the fragmentation of solids. In this paper, a particle method for non-linear elastodynamics of compressible and incompressible materials is developed based on a discretization of the Lagrangian from which the governing equations of elastodynamics are derived using the principle of least action. The discretized Lagrangian leads to a finite-dimensional Hamiltonian system via the Legendre transformation. If the material is incompressible, the Hamiltonian system is accompanied by holonomic constraints. Depending on whether the material is compressible or incompressible, the symplectic scheme adopted for numerical time integration is either the Störmer/ Verlet scheme or the RATTLE method, respectively. The resulting particle method inherits the symplectic structure possessed by the governing equations of elastodynamics. In the case of incompressible materials, incompressibility is strictly enforced at each time step. Some numerical tests indicate the excellence of the method for conservation of mechanical energy besides that of linear and angular momenta.

ジャーナルInternational Journal for Numerical Methods in Engineering
出版ステータスPublished - 2008 5月 21

ASJC Scopus subject areas

  • 工学(その他)
  • 応用数学
  • 計算力学


「A Hamiltonian particle method for non-linear elastodynamics」の研究トピックを掘り下げます。これらがまとまってユニークなフィンガープリントを構成します。