如何获得网格的面积？

Question

如果细分为三角形，如何获得网格的面积？

我找到了每个三角形的叉积并计算了面积。根据搅拌机，面积出错了。我查看了堆栈流其他帖子，但它们没有任何帮助。你能帮我弄清楚为什么我会得到一个我的网格的低面积为 16.3。

 for i in f:
   # print("i",i)
    for k in i:
    # print("k",k)
     for j in z:
     # print("z",z)

      if (k==z.index(j)):


        f_v.append(j)
   # print(f_v)

    v0=np.array(f_v[0])
    v1=np.array(f_v[1])
    v2=np.array(f_v[2])

    ax=np.subtract(f_v[1],f_v[0])
    ax=np.subtract(f_v[2],f_v[1])
    ay=np.subtract(f_v[3],f_v[1])
   ..

    cxx=np.power(cx,2)
 #   cyy=np.power(ay,2)
    #czz=np.power(cz,2)

Answer 1

我创建了meshplex来帮助您完成此类任务，并且可以快速完成。 只需加载网格并总结单元体积：

import meshplex

points = numpy.array([[0.0, 0.0, 0.0], [1.0, 0.0, 0.0], [0.0, 1.0, 0.0]])
cells = numpy.array([[0, 1, 2]])
mesh = meshplex.MeshTri(points, cells)
# or read it from a file
# mesh = meshplex.read("circle.vtk")

print(numpy.sum(mesh.cell_volumes))