The primary computational bottle-neck in implicit structural dynamics is the repeated inversion of the underlying stiffness matrix. In this paper, a fast inversion technique is proposed by merging four distinct but complementary concepts: (1) voxelization with adaptive local refinement, (2) assembly-free (a.k.a. matrix-free or element-by-element) finite element analysis (FEA), (3) assembly-free deflated conjugate gradient (AF-DCG), and (4) multicore parallelization. In particular, we apply these concepts to the well-known Newmark-beta method, and the resulting AF-DCG is well-suited for large-scale problems. It can be easily ported to many-core central processing unit (CPU) and multicore graphics-programmable unit (GPU) architectures, as demonstrated through numerical experiments.
Issue Section:
Research Papers
References
1.
Cook
, R. D.
, Malkus
, D. S.
, Plesha
, M. E.
, and Witt
, R.
, 2002
, Concepts and Applications of Finite Element Analysis
, 4th ed., John Wiley & Sons
, New York.2.
Cook
, R. D.
, 2007
, Concepts and Applications of Finite Element
, John Wiley & Sons
, New York
.3.
Newmark
, N. M.
, 1959
, “A Method of Computation for Structural Dynamics
,” J. Eng. Mech. Div.
, 85
(EM3
), pp. 67
–94
.4.
5.
Arbenz
, P.
, van Lenthe
, G. H.
, Mennel
, U.
, Müllerand
, R.
, and Sala
, M.
, 2008
, “A Scalable Multi-Level Preconditioner for Matrix-Free μ-Finite Element Analysis of Human Bone Structures
,” Int. J. Numer. Methods Eng.
, 73
(7
), pp. 927
–947
.10.1002/nme.21016.
Saad
, Y.
, 2003
, Iterative Methods for Sparse Linear Systems
, SIAM
, Philadelphia, PA.10.1137/1.97808987180037.
Adams
, M.
, 2002
, “Evaluation of Three Unstructured Multigrid Methods on 3D Finite Element Problems in Solid Mechanics
,” Int. J. Numer. Methods Eng.
, 55
(2
), pp. 519
–534
.10.1002/nme.5068.
Bell
, N.
, 2008
, “Efficient Sparse Matrix-Vector Multiplication on CUDA
,” Nvidia Corporation, Nvidia Technical Report No. NVR-2008-004.9.
Hughes
, T. J. R.
, Levit
, I.
, and Winget
, J.
, 1983
, “An Element-by-Element Solution Algorithm for Problems of Structural and Solid Mechanics
,” Comput. Methods Appl. Mech. Eng.
, 36
(2
), pp. 241
–254
.10.1016/0045-7825(83)90115-910.
Augarde
, C. E.
, Ramage
, A.
, and Staudacher
, J.
, 2006
, “An Element-Based Displacement Preconditioner for Linear Elasticity Problems
,” Comput. Struct.
, 84
(31–32
), pp. 2306
–2315
.10.1016/j.compstruc.2006.08.05711.
Bellavia
, S.
, Gondzio
, J.
, and Morini
, B.
, 2013
, “A Matrix-Free Preconditioner for Sparse Symmetric Positive Definite Systems and Least-Squares Problems
,” SIAM J. Sci. Comput.
, 35
(1
), pp. A192
–A211
.10.1137/11084081912.
Suresh
, K.
, and Yadav
, P.
, 2012
, “Large-Scale Modal Analysis on Multi-Core Architectures
,” ASME
Paper No. DETC2012-70281. 10.1115/DETC2012-7028113.
Mueller
, E.
, Guo
, X.
, Scheichl
, R.
, and Shi
, S.
, 2013
, “Matrix-Free GPU Implementation of a Preconditioned Conjugate Gradient Solver for Anisotropic Elliptic PDEs
,” e-print arXiv:1302.7193[cs, math].14.
Akbariyeh
, A.
, Carrigan
, T. J.
, Dennis
, B. H.
, Chan
, W. S.
, Wang
, B. P.
, and Lawrence
, K. L.
, 2012
, “Application of GPU-Based Computing to Large Scale Finite Element Analysis of Three-Dimensional Structures
,” Proceedings of the Eighth International Conference on Engineering Computational Technology
, Stirlingshire, UK, Paper No. 6.15.
Saad
, Y.
, Yeung
, M.
, Erhel
, J.
, and Guyomarc'h
, F.
, 2000
, “A Deflated Version of the Conjugate Gradient Algorithm
,” SIAM J. Sci. Comput.
, 21
(5
), pp. 1909
–1926
.10.1137/S106482959833976116.
Duster
, A.
, Parvizian
, J.
, Yang
, Z.
, and Rank
, E.
, 2008
, “The Finite Cell Method for 3D problems of Solid Mechanics
,” Comput. Methods Appl. Mech. Eng.
, 197
(45–48), pp. 3768
–3782
.10.1016/j.cma.2008.02.03617.
Karabassi
, E. A.
, Papaioannou
, G.
, and Theoharis
, T.
, 1999
, “A Fast Depth-Buffer-Based Voxelization Algorithm
,” J. Graphics Tools
, 4
(4
), pp. 5
–10
.10.1080/10867651.1999.1048751018.
Zienkiewicz
, O. C.
, and Taylor
, R. L.
, 2005
, The Finite Element Method for Solid and Structural Mechanics
, Elsevier
, Oxford, UK.19.
Wilson
, E. L.
, Taylor
, R. L.
, Doherty
, W. P.
, and Ghaboussi
, J.
, 1973
, “Incompatible Displacement Models
,” Numerical and Computer Methods in Structural Mechanics
, Vol. A
, Academic Press, New York, pp. 43
–57
.20.
Taiebat
, H. H.
, and Carter
, J. P.
, 2001
, “Three-Dimensional Non-Conforming Elements
,” Centre for Geotechnical Research, The University of Sydney, Sydney, Australia, Report No. R808.21.
Goddeke
, D.
, Strzodka
, R.
, and Turek
, S.
, 2007
, “Performance and Accuracy of Hardware-Oriented Native-Emulated- and Mixed-Precision Solvers in FEM Simulations
,” Int. J. Parallel Emergent Distrib. Syst.
, 22
(4
), pp. 221
–256
.10.1080/1744576060112207622.
Melanz
, D.
, Khude
, N.
, Jayakumar
, P.
, and Negrut
, D.
, 2014
, “A Matrix-Free Newton–Krylov Parallel Implicit Implementation of the Absolute Nodal Coordinate Formulation
,” ASME J. Comput. Nonlinear Dyn.
, 9
(1
), p. 011006
.10.1115/1.402528123.
Aubry
, R.
, Mut
, F.
, Dey
, S.
, and Lohner
, R.
, 2011
, “Deflated Preconditioned Conjugate Gradient Solvers for Linear Elasticity
,” Int. J. Numer. Methods Eng.
, 88
(11
), pp. 1112
–1127
.10.1002/nme.320924.
Briggs
, W. L.
, Henson
, V. E.
, and McCormick
, S. F.
, 2000
, A Multigrid Tutorial
, SIAM
, Philadelphia, PA.25.
NVIDIA Corporation
, 2008
, NVIDIA CUDA: Compute Unified Device Architecture, Programming Guide
, NVIDIA Corporation
, Santa Clara, CA
.26.
Copyright © 2015 by ASME
You do not currently have access to this content.
