Isogeometric analysis for explicit elastodynamics using a dual-basis diagonal mass formulation

Cosmin Anitescu, Chuong Nguyen, Timon Rabczuk, Xiaoying Zhuang

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

We propose a method to obtain diagonal mass matrices for NURBS-based approximation spaces by a “dual lumping” method. The use of lumped mass matrices is of great importance in elastodynamics problems, as they can be employed in explicit time integration schemes which do not require the solution of a linear system. In finite elements, several well-established methods, such as row-sum, diagonal scaling, or nodal quadrature methods have been used to obtain lumped mass matrices for different applications. However, for higher-order and higher continuity approximation spaces such as those derived from NURBS, these approaches have only limited (second-order) accuracy. In this work, we derive a dual basis which has optimal approximation and dispersion properties, while maintaining local support. The dual space has discontinuities at the element boundaries (knots) and it is used to provide the test functions in the context of a Petrov–Galerkin method. This results in a general framework for the study of lumped mass matrices which can be employed in explicit time integration schemes with high-order accuracy. Numerical experiments are presented to demonstrate the applicability of the method to problems with smooth solutions as well as to wave propagation problems with reduced regularity.

Original languageEnglish
Pages (from-to)574-591
Number of pages18
JournalComputer Methods in Applied Mechanics and Engineering
Volume346
DOIs
Publication statusPublished - 2019 Apr 1

Fingerprint

elastodynamics
formulations
matrices
lumping
approximation
linear systems
regularity
quadratures
continuity
Wave propagation
Linear systems
wave propagation
discontinuity
scaling
Experiments

Keywords

  • Bézier extraction
  • Dual-basis
  • Elastodynamics
  • Isogeometric analysis
  • Lumped mass

ASJC Scopus subject areas

  • Computational Mechanics
  • Mechanics of Materials
  • Mechanical Engineering
  • Physics and Astronomy(all)
  • Computer Science Applications

Cite this

Isogeometric analysis for explicit elastodynamics using a dual-basis diagonal mass formulation. / Anitescu, Cosmin; Nguyen, Chuong; Rabczuk, Timon; Zhuang, Xiaoying.

In: Computer Methods in Applied Mechanics and Engineering, Vol. 346, 01.04.2019, p. 574-591.

Research output: Contribution to journalArticle

Anitescu, Cosmin ; Nguyen, Chuong ; Rabczuk, Timon ; Zhuang, Xiaoying. / Isogeometric analysis for explicit elastodynamics using a dual-basis diagonal mass formulation. In: Computer Methods in Applied Mechanics and Engineering. 2019 ; Vol. 346. pp. 574-591.
@article{ccdcdcb5d8324c43a26be948affa5ac3,
title = "Isogeometric analysis for explicit elastodynamics using a dual-basis diagonal mass formulation",
abstract = "We propose a method to obtain diagonal mass matrices for NURBS-based approximation spaces by a “dual lumping” method. The use of lumped mass matrices is of great importance in elastodynamics problems, as they can be employed in explicit time integration schemes which do not require the solution of a linear system. In finite elements, several well-established methods, such as row-sum, diagonal scaling, or nodal quadrature methods have been used to obtain lumped mass matrices for different applications. However, for higher-order and higher continuity approximation spaces such as those derived from NURBS, these approaches have only limited (second-order) accuracy. In this work, we derive a dual basis which has optimal approximation and dispersion properties, while maintaining local support. The dual space has discontinuities at the element boundaries (knots) and it is used to provide the test functions in the context of a Petrov–Galerkin method. This results in a general framework for the study of lumped mass matrices which can be employed in explicit time integration schemes with high-order accuracy. Numerical experiments are presented to demonstrate the applicability of the method to problems with smooth solutions as well as to wave propagation problems with reduced regularity.",
keywords = "B{\'e}zier extraction, Dual-basis, Elastodynamics, Isogeometric analysis, Lumped mass",
author = "Cosmin Anitescu and Chuong Nguyen and Timon Rabczuk and Xiaoying Zhuang",
year = "2019",
month = "4",
day = "1",
doi = "10.1016/j.cma.2018.12.002",
language = "English",
volume = "346",
pages = "574--591",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0045-7825",
publisher = "Elsevier",

}

TY - JOUR

T1 - Isogeometric analysis for explicit elastodynamics using a dual-basis diagonal mass formulation

AU - Anitescu, Cosmin

AU - Nguyen, Chuong

AU - Rabczuk, Timon

AU - Zhuang, Xiaoying

PY - 2019/4/1

Y1 - 2019/4/1

N2 - We propose a method to obtain diagonal mass matrices for NURBS-based approximation spaces by a “dual lumping” method. The use of lumped mass matrices is of great importance in elastodynamics problems, as they can be employed in explicit time integration schemes which do not require the solution of a linear system. In finite elements, several well-established methods, such as row-sum, diagonal scaling, or nodal quadrature methods have been used to obtain lumped mass matrices for different applications. However, for higher-order and higher continuity approximation spaces such as those derived from NURBS, these approaches have only limited (second-order) accuracy. In this work, we derive a dual basis which has optimal approximation and dispersion properties, while maintaining local support. The dual space has discontinuities at the element boundaries (knots) and it is used to provide the test functions in the context of a Petrov–Galerkin method. This results in a general framework for the study of lumped mass matrices which can be employed in explicit time integration schemes with high-order accuracy. Numerical experiments are presented to demonstrate the applicability of the method to problems with smooth solutions as well as to wave propagation problems with reduced regularity.

AB - We propose a method to obtain diagonal mass matrices for NURBS-based approximation spaces by a “dual lumping” method. The use of lumped mass matrices is of great importance in elastodynamics problems, as they can be employed in explicit time integration schemes which do not require the solution of a linear system. In finite elements, several well-established methods, such as row-sum, diagonal scaling, or nodal quadrature methods have been used to obtain lumped mass matrices for different applications. However, for higher-order and higher continuity approximation spaces such as those derived from NURBS, these approaches have only limited (second-order) accuracy. In this work, we derive a dual basis which has optimal approximation and dispersion properties, while maintaining local support. The dual space has discontinuities at the element boundaries (knots) and it is used to provide the test functions in the context of a Petrov–Galerkin method. This results in a general framework for the study of lumped mass matrices which can be employed in explicit time integration schemes with high-order accuracy. Numerical experiments are presented to demonstrate the applicability of the method to problems with smooth solutions as well as to wave propagation problems with reduced regularity.

KW - Bézier extraction

KW - Dual-basis

KW - Elastodynamics

KW - Isogeometric analysis

KW - Lumped mass

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

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

U2 - 10.1016/j.cma.2018.12.002

DO - 10.1016/j.cma.2018.12.002

M3 - Article

AN - SCOPUS:85059037477

VL - 346

SP - 574

EP - 591

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0045-7825

ER -