A Hamiltonian particle method for non-linear elastodynamics

Yukihito Suzuki, Seiichi Koshizuka

Research output: Contribution to journalArticle

25 Citations (Scopus)

Abstract

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.

Original languageEnglish
Pages (from-to)1344-1373
Number of pages30
JournalInternational Journal for Numerical Methods in Engineering
Volume74
Issue number8
DOIs
Publication statusPublished - 2008 May 21
Externally publishedYes

Fingerprint

Hamiltonians
Particle Method
Elastodynamics
Hamiltonian Systems
Governing equation
Discrete Dynamics
Incompressibility
Symplectic Structure
Angular momentum
Meshless
Legendre
Large Deformation
Time Integration
Fragmentation
Angular Momentum
Conservation
Continuum
Strictly
Discretization
Flexibility

Keywords

  • Hamiltonian
  • Incompressible
  • Non-linear elastodynamics
  • Particle method
  • Symplectic

ASJC Scopus subject areas

  • Engineering (miscellaneous)
  • Applied Mathematics
  • Computational Mechanics

Cite this

A Hamiltonian particle method for non-linear elastodynamics. / Suzuki, Yukihito; Koshizuka, Seiichi.

In: International Journal for Numerical Methods in Engineering, Vol. 74, No. 8, 21.05.2008, p. 1344-1373.

Research output: Contribution to journalArticle

@article{7d1ac965ffb54ff7a21b0b8efaaaee3f,
title = "A Hamiltonian particle method for non-linear elastodynamics",
abstract = "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{\"o}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.",
keywords = "Hamiltonian, Incompressible, Non-linear elastodynamics, Particle method, Symplectic",
author = "Yukihito Suzuki and Seiichi Koshizuka",
year = "2008",
month = "5",
day = "21",
doi = "10.1002/nme.2222",
language = "English",
volume = "74",
pages = "1344--1373",
journal = "International Journal for Numerical Methods in Engineering",
issn = "0029-5981",
publisher = "John Wiley and Sons Ltd",
number = "8",

}

TY - JOUR

T1 - A Hamiltonian particle method for non-linear elastodynamics

AU - Suzuki, Yukihito

AU - Koshizuka, Seiichi

PY - 2008/5/21

Y1 - 2008/5/21

N2 - 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.

AB - 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.

KW - Hamiltonian

KW - Incompressible

KW - Non-linear elastodynamics

KW - Particle method

KW - Symplectic

UR - http://www.scopus.com/inward/record.url?scp=43449139410&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=43449139410&partnerID=8YFLogxK

U2 - 10.1002/nme.2222

DO - 10.1002/nme.2222

M3 - Article

VL - 74

SP - 1344

EP - 1373

JO - International Journal for Numerical Methods in Engineering

JF - International Journal for Numerical Methods in Engineering

SN - 0029-5981

IS - 8

ER -