有限元网格自适应加密

有限元一阶四面体单元自适应加密程序

程序说明

运行fem_mesh.exe程序执行加密程序

该程序实现了对部分四面体单元进行加密，需要导入 .xyz 和 .lv 格式的文件。.xyz 文件中储存了点的坐标信息，第一行表示所有点的个数。.lv 文件储存了四面体单元上的点的全局编号信息（按局部编号排列），第一行表示所有四面体单元的个数。

程序对四面体的边进行加密，因此需要导入 .le 文件和 .re 文件。.le 文件储存了每个四面体单元的六条边的全局编号。.re 文件储存了每条边上的误差大小，并按降序排列。

程序运行结束后会输出.xyz 和.lv两个新文件。

关于单元点和边的全局编号说明

• 单元上点的局部编号是按右手螺旋法则进行编号。

• 加密结束后，单元上新的边的全局编号需要重新排序。

• 四面体上的边的局部编号按照以下法则进行编号，左边表示边的局部编号，右边则是连接的两个点的局部编号：

  1. 边1：（1，2）
  2. 边2：（1，3）
  3. 边3：（1，4）
  4. 边4：（2，3）
  5. 边5：（2，4）
  6. 边6：（3，4）

加密算法的说明

1. 读取网格几何加密信息
   2.依次遍历所有单元
2. 遍历该单元的所有边，判断是否需要加密
   4.若某条边需要加密，则新生成一个中点（如果已经生成过中点了，则不需要再生成）
3. 新生成一个单元，并修改旧单元的编号
4. 继续遍历所有单元，直到遍历结束