Parallel computations in nonlinear solid mechanics using adaptive finite element and meshless methods

Zahur Ullah, C. E. Augarde, W. M. Coombs

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

Purpose:– A variety of meshless methods have been developed in the last 20 years with an intention to solve practical engineering problems, but are limited to small academic problems due to associated high computational cost as compared to the standard finite element methods (FEM). The purpose of this paper is to develop an efficient and accurate algorithms based on meshless methods for the solution of problems involving both material and geometrical nonlinearities.Design/methodology/approach:– A parallel two-dimensional linear elastic computer code is presented for a maximum entropy basis functions based meshless method. The two-dimensional algorithm is subsequently extended to three-dimensional adaptive nonlinear and three-dimensional parallel nonlinear adaptively coupled finite element, meshless method cases. The Prandtl-Reuss constitutive model is used to model elasto-plasticity and total Lagrangian formulations are used to model finite deformation. Furthermore, Zienkiewicz and Zhu and Chung and Belytschko error estimation procedure are used in the FE and meshless regions of the problem domain, respectively. The message passing interface library and open-source software packages, METIS and MUltifrontal Massively Parallel Solver are used for the high performance computation.Findings:– Numerical examples are given to demonstrate the correct implementation and performance of the parallel algorithms. The agreement between the numerical and analytical results in the case of linear elastic example is excellent. For the nonlinear problems load-displacement curve are compared with the reference FEM and found in a very good agreement. As compared to the FEM, no volumetric locking was observed in the case of meshless method. Furthermore, it is shown that increasing the number of processors up to a given number improve the performance of parallel algorithms in term of simulation time, speedup and efficiency.Originality/value:– Problems involving both material and geometrical nonlinearities are of practical importance in many engineering applications, e.g. geomechanics, metal forming and biomechanics. A family of parallel algorithms has been developed in this paper for these problems using adaptively coupled finite element, meshless method (based on maximum entropy basis functions) for distributed memory computer architectures.
LanguageEnglish
Pages1161-1191
JournalEngineering Computations
Volume33
Issue number4
DOIs
Publication statusAccepted/In press - 13 Jun 2016

Fingerprint

Mechanics
Finite element method
Parallel algorithms
Entropy
Geomechanics
Computer architecture
Biomechanics
Metal forming
Message passing
Constitutive models
Software packages
Error analysis
Interfaces (computer)
Plasticity
Data storage equipment
Costs

Keywords

  • Parallel computations
  • Adaptively coupled finite element-meshless method
  • Adaptivity
  • Error estimation
  • Maximum entropy shape functions
  • Meshless method

Cite this

Ullah, Zahur ; Augarde, C. E. ; Coombs, W. M. / Parallel computations in nonlinear solid mechanics using adaptive finite element and meshless methods. 2016 ; Vol. 33, No. 4. pp. 1161-1191.
@article{fa76487ce52b4cac951207c8451c0056,
title = "Parallel computations in nonlinear solid mechanics using adaptive finite element and meshless methods",
abstract = "Purpose:– A variety of meshless methods have been developed in the last 20 years with an intention to solve practical engineering problems, but are limited to small academic problems due to associated high computational cost as compared to the standard finite element methods (FEM). The purpose of this paper is to develop an efficient and accurate algorithms based on meshless methods for the solution of problems involving both material and geometrical nonlinearities.Design/methodology/approach:– A parallel two-dimensional linear elastic computer code is presented for a maximum entropy basis functions based meshless method. The two-dimensional algorithm is subsequently extended to three-dimensional adaptive nonlinear and three-dimensional parallel nonlinear adaptively coupled finite element, meshless method cases. The Prandtl-Reuss constitutive model is used to model elasto-plasticity and total Lagrangian formulations are used to model finite deformation. Furthermore, Zienkiewicz and Zhu and Chung and Belytschko error estimation procedure are used in the FE and meshless regions of the problem domain, respectively. The message passing interface library and open-source software packages, METIS and MUltifrontal Massively Parallel Solver are used for the high performance computation.Findings:– Numerical examples are given to demonstrate the correct implementation and performance of the parallel algorithms. The agreement between the numerical and analytical results in the case of linear elastic example is excellent. For the nonlinear problems load-displacement curve are compared with the reference FEM and found in a very good agreement. As compared to the FEM, no volumetric locking was observed in the case of meshless method. Furthermore, it is shown that increasing the number of processors up to a given number improve the performance of parallel algorithms in term of simulation time, speedup and efficiency.Originality/value:– Problems involving both material and geometrical nonlinearities are of practical importance in many engineering applications, e.g. geomechanics, metal forming and biomechanics. A family of parallel algorithms has been developed in this paper for these problems using adaptively coupled finite element, meshless method (based on maximum entropy basis functions) for distributed memory computer architectures.",
keywords = "Parallel computations, Adaptively coupled finite element-meshless method, Adaptivity, Error estimation, Maximum entropy shape functions, Meshless method",
author = "Zahur Ullah and Augarde, {C. E.} and Coombs, {W. M.}",
year = "2016",
month = "6",
day = "13",
doi = "10.1108/EC-06-2015-0166",
language = "English",
volume = "33",
pages = "1161--1191",
number = "4",

}

Parallel computations in nonlinear solid mechanics using adaptive finite element and meshless methods. / Ullah, Zahur; Augarde, C. E.; Coombs, W. M.

Vol. 33, No. 4, 13.06.2016, p. 1161-1191.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Parallel computations in nonlinear solid mechanics using adaptive finite element and meshless methods

AU - Ullah, Zahur

AU - Augarde, C. E.

AU - Coombs, W. M.

PY - 2016/6/13

Y1 - 2016/6/13

N2 - Purpose:– A variety of meshless methods have been developed in the last 20 years with an intention to solve practical engineering problems, but are limited to small academic problems due to associated high computational cost as compared to the standard finite element methods (FEM). The purpose of this paper is to develop an efficient and accurate algorithms based on meshless methods for the solution of problems involving both material and geometrical nonlinearities.Design/methodology/approach:– A parallel two-dimensional linear elastic computer code is presented for a maximum entropy basis functions based meshless method. The two-dimensional algorithm is subsequently extended to three-dimensional adaptive nonlinear and three-dimensional parallel nonlinear adaptively coupled finite element, meshless method cases. The Prandtl-Reuss constitutive model is used to model elasto-plasticity and total Lagrangian formulations are used to model finite deformation. Furthermore, Zienkiewicz and Zhu and Chung and Belytschko error estimation procedure are used in the FE and meshless regions of the problem domain, respectively. The message passing interface library and open-source software packages, METIS and MUltifrontal Massively Parallel Solver are used for the high performance computation.Findings:– Numerical examples are given to demonstrate the correct implementation and performance of the parallel algorithms. The agreement between the numerical and analytical results in the case of linear elastic example is excellent. For the nonlinear problems load-displacement curve are compared with the reference FEM and found in a very good agreement. As compared to the FEM, no volumetric locking was observed in the case of meshless method. Furthermore, it is shown that increasing the number of processors up to a given number improve the performance of parallel algorithms in term of simulation time, speedup and efficiency.Originality/value:– Problems involving both material and geometrical nonlinearities are of practical importance in many engineering applications, e.g. geomechanics, metal forming and biomechanics. A family of parallel algorithms has been developed in this paper for these problems using adaptively coupled finite element, meshless method (based on maximum entropy basis functions) for distributed memory computer architectures.

AB - Purpose:– A variety of meshless methods have been developed in the last 20 years with an intention to solve practical engineering problems, but are limited to small academic problems due to associated high computational cost as compared to the standard finite element methods (FEM). The purpose of this paper is to develop an efficient and accurate algorithms based on meshless methods for the solution of problems involving both material and geometrical nonlinearities.Design/methodology/approach:– A parallel two-dimensional linear elastic computer code is presented for a maximum entropy basis functions based meshless method. The two-dimensional algorithm is subsequently extended to three-dimensional adaptive nonlinear and three-dimensional parallel nonlinear adaptively coupled finite element, meshless method cases. The Prandtl-Reuss constitutive model is used to model elasto-plasticity and total Lagrangian formulations are used to model finite deformation. Furthermore, Zienkiewicz and Zhu and Chung and Belytschko error estimation procedure are used in the FE and meshless regions of the problem domain, respectively. The message passing interface library and open-source software packages, METIS and MUltifrontal Massively Parallel Solver are used for the high performance computation.Findings:– Numerical examples are given to demonstrate the correct implementation and performance of the parallel algorithms. The agreement between the numerical and analytical results in the case of linear elastic example is excellent. For the nonlinear problems load-displacement curve are compared with the reference FEM and found in a very good agreement. As compared to the FEM, no volumetric locking was observed in the case of meshless method. Furthermore, it is shown that increasing the number of processors up to a given number improve the performance of parallel algorithms in term of simulation time, speedup and efficiency.Originality/value:– Problems involving both material and geometrical nonlinearities are of practical importance in many engineering applications, e.g. geomechanics, metal forming and biomechanics. A family of parallel algorithms has been developed in this paper for these problems using adaptively coupled finite element, meshless method (based on maximum entropy basis functions) for distributed memory computer architectures.

KW - Parallel computations

KW - Adaptively coupled finite element-meshless method

KW - Adaptivity

KW - Error estimation

KW - Maximum entropy shape functions

KW - Meshless method

U2 - 10.1108/EC-06-2015-0166

DO - 10.1108/EC-06-2015-0166

M3 - Article

VL - 33

SP - 1161

EP - 1191

IS - 4

ER -