| Line 1: | Line 1: | ||
== Abstract == | == Abstract == | ||
| − | We examine a residual and matrix-free Jacobian formulation of compressible and nearly incompressible ( 0.5) displacement-only linear isotropic elasticity with high-order hexahedral finite elements. A matrix-free p-multigrid method is combined with algebraic multigrid on the assembled sparse coarse grid matrix to provide an effective preconditioner. The software is verified with the method of manufactured solutions. We explore convergence to a predetermined | + | We examine a residual and matrix-free Jacobian formulation of compressible and nearly incompressible (v → 0.5) displacement-only linear isotropic elasticity with high-order hexahedral finite elements. A matrix-free p-multigrid method is combined with algebraic multigrid on the assembled sparse coarse grid matrix to provide an effective preconditioner. The software is verified with the method of manufactured solutions. We explore convergence to a predetermined L<sub>2</sub> error of 10<sup>-4</sup>, 10<sup>-5</sup> and 10<sup>-6</sup> for the compressible case and 10<sup>-4</sup>, 10<sup>-5</sup> for the nearly-incompressible cases, as the Poisson's ratio approaches 0.5, based upon grid resolution and polynomial order. We compare our results against results obtained from C3D20H mixed/hybrid element available in the commercial finite element software ABAQUS that is quadratic in displacement and linear in pressure. We determine, for the same problem size, that our matrix-free approach for displacement-only implementation is faster and more efficient for quadratic elements compared to the C3D20H element from ABAQUS that is specially designed to handle nearly-incompressible and incompressible elasticity problems. However, as we approach the near incompressibility limit, the number of Conjugate Gradient iterations required to achieve the desired solution increases significantly. |
== Full document == | == Full document == | ||
<pdf>Media:Draft_Content_929931938p3551.pdf</pdf> | <pdf>Media:Draft_Content_929931938p3551.pdf</pdf> | ||
Latest revision as of 09:27, 12 March 2021
Abstract
We examine a residual and matrix-free Jacobian formulation of compressible and nearly incompressible (v → 0.5) displacement-only linear isotropic elasticity with high-order hexahedral finite elements. A matrix-free p-multigrid method is combined with algebraic multigrid on the assembled sparse coarse grid matrix to provide an effective preconditioner. The software is verified with the method of manufactured solutions. We explore convergence to a predetermined L2 error of 10-4, 10-5 and 10-6 for the compressible case and 10-4, 10-5 for the nearly-incompressible cases, as the Poisson's ratio approaches 0.5, based upon grid resolution and polynomial order. We compare our results against results obtained from C3D20H mixed/hybrid element available in the commercial finite element software ABAQUS that is quadratic in displacement and linear in pressure. We determine, for the same problem size, that our matrix-free approach for displacement-only implementation is faster and more efficient for quadratic elements compared to the C3D20H element from ABAQUS that is specially designed to handle nearly-incompressible and incompressible elasticity problems. However, as we approach the near incompressibility limit, the number of Conjugate Gradient iterations required to achieve the desired solution increases significantly.
Full document
Document information
Published on 10/03/21
Submitted on 10/03/21
Volume 1400 - Software, High Performance Computing, 2021
DOI: 10.23967/wccm-eccomas.2020.302
Licence: CC BY-NC-SA license
Share this document
Keywords
claim authorship
Are you one of the authors of this document?