Efficient mass and stiffness matrix assembly via weighted Gaussian quadrature rules for B-splines

DOI10.1016/J.CAM.2019.112626zbMATH Open1493.65054arXiv1710.01048OpenAlexW2990351226MaRDI QIDQ2297117FDOQ2297117

Authors: Michael Bartoň, Vladimir Evgenievich Puzyrev, Quanling Deng, Victor M. Calo 

Publication date: 18 February 2020

Published in: Journal of Computational and Applied Mathematics (Search for Journal in Brave)

Abstract: Calabro et al. (2017) changed the paradigm of the mass and stiffness computation from the traditional element-wise assembly to a row-wise concept, showing that the latter one offers integration that may be orders of magnitude faster. Considering a B-spline basis function as a non-negative measure, each mass matrix row is integrated by its own quadrature rule with respect to that measure. Each rule is easy to compute as it leads to a linear system of equations, however, the quadrature rules are of the Newton-Cotes type, that is, they require a number of quadrature points that is equal to the dimension of the spline space. In this work, we propose weighted quadrature rules of Gaussian type which require the minimum number of quadrature points while guaranteeing exactness of integration with respect to the weight function. The weighted Gaussian rules arise as solutions of non-linear systems of equations. We derive rules for the mass and stiffness matrices for uniform $C^1$ quadratic and $C^2$ cubic isogeometric discretizations. Our rules further reduce the number of quadrature points by a factor of $(fracp+12p+1{)}^d$ when compared to Calabro et al. (2017), $p$ being the polynomial degree and $d$ the dimension of the problem, and consequently reduce the computational cost of the mass and stiffness matrix assembly by a similar factor.

Full work available at URL: https://arxiv.org/abs/1710.01048

Recommendations

• Fast formation and assembly of finite element matrices with application to isogeometric linear elasticity
• Efficient assembly based on B-spline tailored quadrature rules for the IgA-SGBEM
• Fast formation of isogeometric Galerkin matrices by weighted quadrature
• Weighted quadrature for hierarchical B-splines
• Quadrature rules in the isogeometric Galerkin method: state of the art and an introduction to weighted quadrature

zbMATH Keywords

B-splinesisogeometric analysismass and stiffness matrix assemblyweighted Gaussian quadrature

Mathematics Subject Classification ID

Numerical computation using splines (65D07) Numerical quadrature and cubature formulas (65D32) Approximate quadratures (41A55)

Cites Work

• A nonlinear optimization procedure for generalized Gaussian quadratures
• Isogeometric analysis. Toward integration of CAD and FEA
• Efficient quadrature for NURBS-based isogeometric analysis
• Reduced Bézier element quadrature rules for quadratic and cubic splines in isogeometric analysis
• Selective and reduced numerical integrations for NURBS-based isogeometric analysis
• A simple algorithm for obtaining nearly optimal quadrature rules for NURBS-based isogeometric analysis
• Weighted quadrature rules for finite element methods
• Gauss quadrature for refineable weight functions
• Moment Theory for Weak Chebyshev Systems with Applications to Monosplines, Quadrature Formulae and Best One-Sided $L^1 $-Approximation by Spline Functions with Fixed Knots
• Isogeometric analysis with geometrically continuous functions on two-patch geometries
• Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization
• The variational collocation method
• Optimal quadrature for univariate and tensor product splines
• Optimal quadrature rules for odd-degree spline spaces and their application to tensor-product-based isogeometric analysis
• Gaussian quadrature for splines via homotopy continuation: rules for \(C^2\) cubic splines
• Dispersion-minimized mass for isogeometric analysis
• Dispersion optimized quadratures for isogeometric analysis
• Spline functions, convex curves and mechanical quadrature
• Integration by interpolation and look-up for Galerkin-based isogeometric analysis
• Fast formation of isogeometric Galerkin matrices by weighted quadrature
• Optimal and reduced quadrature rules for tensor product and hierarchically refined splines in isogeometric analysis
• Efficient assembly based on B-spline tailored quadrature rules for the IgA-SGBEM
• Gaussian quadrature rules for \(C^1\) quintic splines with uniform knot vectors
• Exploring matrix generation strategies in isogeometric analysis
• Dispersion-optimized quadrature rules for isogeometric analysis: modified inner products, their dispersion properties, and optimally blended schemes
• Reduced integration at superconvergent points in isogeometric analysis
• Dispersion-minimizing quadrature rules for \(C^1\) quadratic isogeometric analysis

Cited In (13)

• Efficient assembly based on B-spline tailored quadrature rules for the IgA-SGBEM
• Refined isogeometric analysis for generalized Hermitian eigenproblems
• Machine learning discovery of optimal quadrature rules for isogeometric analysis
• Efficient Quadrature Rules for Computing the Stiffness Matrices of Mass-Lumped Tetrahedral Elements for Linear Wave Problems
• Quadrature rules for \(C^1\) quadratic spline finite elements on the Powell-Sabin 12-split
• Tuning three-dimensional tumor progression simulations on a cluster of GPGPUs
• Quadrature rules in the isogeometric Galerkin method: state of the art and an introduction to weighted quadrature
• Concurrent algorithm for integrating three-dimensional B-spline functions into machines with shared memory such as GPU
• High‐fidelity tensor‐decomposition based matrix formation for isogeometric buckling analysis of laminated shells with solid‐shell formulation
• Efficient matrix computation for isogeometric discretizations with hierarchical B-splines in any dimension
• Fast formation and assembly of finite element matrices with application to isogeometric linear elasticity
• Weighted quadrature for hierarchical B-splines
• Meshless physics-informed deep learning method for three-dimensional solid mechanics

This page was built for publication: Efficient mass and stiffness matrix assembly via weighted Gaussian quadrature rules for B-splines

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2297117)