文件名称:用于 2D/3D FEM 矩阵的快速 MATLAB 组装函数:用于 2D/3D 泊松和线性弹性问题的 MATLAB 组装函数-matlab开发
文件大小:17KB
文件格式:ZIP
更新时间:2024-06-18 15:56:05
matlab
我们提出了 P 1 有限元方法的 MATLAB 实现,用于二维 (2D) 和三维 (3D) 中泊松问题和线性弹性问题的数值解。 该代码由矩阵(质量和刚度)和右侧的矢量化(和短)组装函数组成。 由于对于 P 1 有限元,单元质量矩阵和右侧很简单,因此实现仅在单元体积上使用 MATLAB 函数稀疏。 对于刚度矩阵,为了获得接近标准形式的 MATLAB 实现,使用元胞数组来存储元素基函数的梯度。 然后,组装过程可以在小尺寸单元阵列上使用矩阵/矢量乘积。 数值实验表明,我们的实现在时间上是快速且可扩展的。 引用Koko J.,使用元胞阵列方法在 2D 和 3D 中进行 FEM 矩阵的快速 MATLAB 组装,国际建模、模拟和科学计算杂志卷。 7、第3期(2016)
【文件预览】:
KPDE.zip