2017-05-30 142 views
0

我正在尝试生成一个多维数据集的四面体网格。我遇到的问题是,用于生成三角面的几个节点不会用于生成任何四面体元素。有没有办法避免这种情况? 这是代码。GMSH中2D和3D元素之间的节点不一致性

cl__1 = 0.01; 
    Point(1) = {0, 0, 0, cl__1}; 
    Point(2) = {1, 0, 0, cl__1}; 
    Point(3) = {1, 1, 0, cl__1}; 
    Point(4) = {0, 1, 0, cl__1}; 
    Line(1) = {1, 2}; 
    Line(2) = {2, 3}; 
    Line(3) = {3, 4}; 
    Line(4) = {4, 1}; 
    Line Loop(5) = {4, 1, 2, 3}; 
    Plane Surface(6) = {5}; 
    Extrude {0, 0, 1} { 
     Surface{6}; 
    } 
    Coherence; 
    Transfinite Line {1, 2, 3, 4, 9, 10, 11, 8, 18, 22, 13, 14} = 2 Using Progression 1; 
    Transfinite Surface {19}; 
    Transfinite Surface {23}; 
    Transfinite Surface {27}; 
    Transfinite Surface {15}; 
    Transfinite Surface {28}; 
    Transfinite Surface {6}; 
    Transfinite Volume{1} = {1, 2, 3, 4, 6, 10, 14, 5}; 
    Physical Surface("top") = {28}; 
    Physical Surface("bottom") = {6}; 
    Physical Surface("x_min") = {15}; 
    Physical Surface("x_max") = {23}; 
    Physical Surface("y_min") = {19}; 
    Physical Surface("y_max") = {27}; 
    Physical Volume("bottom_volume") = {1}; 

由于 蹇

+0

欢迎来到Stack Overflow。如果可能的话,尝试添加你的问题的一些代码和/或你已经完成的工作。好运 –

+0

谢谢。我更新了我的代码。 – Jian

+0

谢谢,如果可能的话,给你添加标签问题,指出使用的语言和其他提示,以便用户可以找到你的问题。 GL –

回答

0

该输入的结果给出了这些节点:

$Nodes 
8 
1 0 0 0 
2 1 0 0 
3 1 1 0 
4 0 1 0 
5 0 1 1 
6 0 0 1 
7 1 0 1 
8 1 1 1 
$EndNodes 

与这些节点可以定义所有的三角形面和定义该网格中的四面体单元(有在我看来没有不一致)。

相关问题