We demonstrate the potential of collocation methods for efficient higher-order analysis on standard nodal finite element meshes. We focus on a collocation method that is variationally consistent and geometrically flexible, converges optimally, embraces concepts of reduced quadrature, and leads to symmetric stiffness and diagonal consistent mass matrices. At the same time, it minimizes the evaluation cost per quadrature point, thus reducing formation and assembly effort significantly with respect to standard Galerkin finite element methods. We provide a detailed review of all components of the technology in the context of elastodynamics, that is, weighted residual formulation, nodal basis functions on Gauss-Lobatto quadrature points, and symmetrization by averaging with the ultra-weak formulation. We quantify potential gains by comparing the computational efficiency of collocated and standard finite elements in terms of basic operation counts and timings. Our results show that collocation is significantly less expensive for problems dominated by the formation and assembly effort, such as higher-order elastostatic analysis. Furthermore, we illustrate the potential of collocation for efficient higher-order explicit dynamics. Throughout this work, we advocate a straightforward implementation based on simple modifications of standard finite element codes. We also point out the close connection to spectral element methods, where many of the key ideas are already established.
AbstractWe demonstrate the potential of collocation methods for efficient higher-order analysis on standard nodal finite element meshes. We focus on a collocation method that is variationally consistent and geometrically flexible, converges optimally, embraces concepts of reduced quadrature, and leads to symmetric stiffness and diagonal consistent mass matrices. At the same time, it minimizes the evaluation cost per quadrature point, thus reducing formation and assembly effort significantly with respect to standard Galerkin finite element methods. We provide a detailed review of all components of the technology in the context of elastodynamics, i.e. weighted residual formulation, nodal basis functions on Gauss-Lobatto quadrature points, and symmetrization by averaging with the ultra-weak formulation. We quantify potential gains by comparing the computational efficiency of collocation and standard finite elements in terms of basic operation counts and timings. Our results show that collocation is significantly less expensive for problems dominated by the formation and assembly effort, such as higher-order elastostatic analysis. Furthermore we illustrate the potential of collocation for efficient higher-order explicit dynamics. Throughout this work, we advocate a straightforward implementation based on simple modifications of standard finite element codes. We also point out the close connection to spectral element methods, where many of the key ideas are already established.