中文

基于 B 样条加权高斯求积规则的高效质量与刚度矩阵组装

数值分析 2017-10-04 v1

摘要

Calabro 等人 (2017) 将质量与刚度计算的范式从传统的逐单元组装转变为逐行概念,表明后者提供的积分速度可快数个数量级。将 B 样条基函数视为非负测度,每个质量矩阵行通过其自身关于该测度的求积规则进行积分。每条规则易于计算,因为它归结为一个线性方程组,然而这些求积规则属于 Newton-Cotes 型,即所需求积节点数等于样条空间的维数。在本工作中,我们提出高斯型加权求积规则,它以最少的求积节点数保证关于权函数的积分精确性。加权高斯规则作为非线性方程组的解出现。我们为均匀 C1C^1 二次和 C2C^2 三次等几何离散推导了质量矩阵和刚度矩阵的规则。与 Calabro 等人 (2017) 相比,我们的规则将求积节点数进一步减少了 (p+12p+1)d(\frac{p+1}{2p+1})^d 倍(pp 为多项式次数,dd 为问题维数),从而将质量与刚度矩阵组装的计算成本降低了类似的倍数。

关键词

引用

@article{arxiv.1710.01048,
  title  = {Efficient mass and stiffness matrix assembly via weighted Gaussian quadrature rules for B-splines},
  author = {Michael Bartoň and Vladimir Puzyrev and Quanling Deng and Victor Calo},
  journal= {arXiv preprint arXiv:1710.01048},
  year   = {2017}
}