主要内容

网格数据作为[p,e,t]三元组

偏微分方程工具箱™使用带有三角形元素的网格用于二维几何图形,使用带有四面体元素的网格用于三维几何图形。偏微分方程工具箱的早期版本使用网格的形式[p, e t]三。的矩阵pe,t分别表示网格的点(节点)、元素和三角形或四面体。工具箱的后续版本支持[p, e t]网格的兼容性原因。

请注意

新特性可能与遗留工作流不兼容。有关推荐工作流中的网格数据的描述,请参见网格数据

二维网格的网格数据有以下组件:

  • p(点,网格节点)是一个2——- - - - - -Np节点矩阵,其中Np是网格中的节点数。每一列p (:, k)x点的-坐标kp (k)y点的-坐标kp (2 k)

  • e(边)是a7——- - - - - -边的矩阵,其中是网格中的边数。网格边缘向内e几何图形的边缘是一一对应的。的e矩阵表示几何的离散边,其方式与t矩阵表示离散面。的每一栏e矩阵表示一条边。

    • e (k)网格边第一个点的索引是多少k

    • e (2 k)第二点的索引是否在网格边缘k

    • e (k)参数值是否在边的第一个点k.该参数值与沿几何边缘的弧长有关。

    • e (4 k)参数值是否在边的第二个点k

    • e (5 k)是包含网格边的几何边的ID。可以使用该命令查看边缘idpdegplot(几何学,“EdgeLabels”,“上”)

    • e (k)是边缘左侧的子域名号。沿边缘的方向是通过增加参数值来给出的。的子域名0是几何图形的外部。

    • e (7 k)是边缘右侧的子域名号。

  • t(三角形)是a4——- - - - - -Nt三角形的矩阵7——- - - - - -Nt三角形的矩阵,取决于你是否调用generateMeshGeometricOrder名称-值对设置为“二次”“线性”,分别。initmesh创建只“线性”元素,具有大小4——- - - - - -NtNt是网格中三角形的数量。每一栏t中的点的索引p形成了三角形。例外的是列中的最后一个条目,它是子域名号。三角形点的顺序如下所示。

线性三角形元素的节点从最左边的节点开始逆时针编号1,2,3。二次三角形元素的节点也一样,每条边中间的附加节点编号为4,5,6。

三维网格的网格数据有以下组成部分:

  • p(点,网格节点)是一个3.——- - - - - -Np节点矩阵,其中Np是网格中的节点数。每一列p (:, k)x点的-坐标kp (k),y点的-坐标kp (2 k),以及z点的-坐标kp (3 k)

  • e是一个将网格面与几何边界关联起来的对象。在将几何边界上设置的边界条件转换为网格边界面时,偏微分方程工具箱函数使用这种关联。

  • t(四面体)要么是11乘的Nt四面体的矩阵或5 ×矩阵Nt四面体的矩阵,取决于你是否调用generateMeshGeometricOrder名称-值对设置为“二次”“线性”,分别。Nt是网格中四面体的个数。每一栏t中的点的索引p形成了四面体。例外的是列中的最后一个元素,它是子域名号。四面体点的排列顺序如图所示。

线性四面体单元的节点编号如下。底部的节点从最左边的节点开始逆时针依次为1、2、3。节点4在四面体的顶部。二次四面体单元的节点是相同的,只是在每条边的中间增加了节点。这些节点在四面体的底部编号为5,6,7,在四面体的两侧编号为8,9,10。

您可以创建[p, e t]通过使用以下方法之一进行网格划分:

  • 使用initmesh函数来创建二维图像[p, e t]网。

  • 使用generateMesh函数来创建二维或三维网格FEMesh对象。然后使用meshToPet函数将网格转换为a[p, e t]网。

Baidu
map