gmsh生成非结构网格详解(之二).docVIP

Gmsh 生成非结构网格详解（之二） 利用Gmsh软件可以快速生成三维非结构网格，我们下面可以看一个例子，例子在gmsh的tutorial目录中： 下面是 t5.geo的完整说明： // 先定义几个特征长度： lcar1 = .1; lcar2 = .0005; lcar3 = .075; // 为了能够不修改文件，而全局的改变这些长度，我们能够在 // 命令行指定一个因子或者指定 Mesh.CharacteristicLengthFactor // 选项来达成。例如 // // $ gmsh t5 -clscale 1 // // 将会产生一个大约两千个节点和一万个四面体的网格（在一台 // 666M 主频的 Alpha Workstation 上需要约 3 秒钟）。而 // // $ gmsh t5 -clscale 0.2 // // 会将所有的特征长度都缩小 5 倍，得到的网格将有大约 17 万 // 个节点，和大约 100 万个四面体，在同样的机器上，这需要花 // 费大约 16 分钟的时间，Gmsh 还需要很多改进来达到更高的效 // 率。 Point(1) = {0.5,0.5,0.5,lcar2}; Point(2) = {0.5,0.5,0,lcar1}; Point(3) = {0,0.5,0.5,lcar1}; Point(4) = {0,0,0.5,lcar1}; Point(5) = {0.5,0,0.5,lcar1}; Point(6) = {0.5,0,0,lcar1}; Point(7) = {0,0.5,0,lcar1}; Point(8) = {0,1,0,lcar1}; Point(9) = {1,1,0,lcar1}; Point(10) = {0,0,1,lcar1}; Point(11) = {0,1,1,lcar1}; Point(12) = {1,1,1,lcar1}; Point(13) = {1,0,1,lcar1}; Point(14) = {1,0,0,lcar1}; Line(1) = {8,9}; Line(2) = {9,12}; Line(3) = {12,11}; Line(4) = {11,8}; Line(5) = {9,14}; Line(6) = {14,13}; Line(7) = {13,12}; Line(8) = {11,10}; Line(9) = {10,13}; Line(10) = {10,4}; Line(11) = {4,5}; Line(12) = {5,6}; Line(13) = {6,2}; Line(14) = {2,1}; Line(15) = {1,3}; Line(16) = {3,7}; Line(17) = {7,2}; Line(18) = {3,4}; Line(19) = {5,1}; Line(20) = {7,8}; Line(21) = {6,14}; Line Loop(22) = {11,19,15,18}; Plane Surface(23) = {22}; Line Loop(24) = {16,17,14,15}; Plane Surface(25) = {24}; Line Loop(26) = {-17,20,1,5,-21,13}; Plane Surface(27) = {26}; Line Loop(28) = {4,1,2,3}; Plane Surface(29) = {28}; Line Loop(30) = {7,-2,5,6}; Plane Surface(31) = {30}; Line Loop(32) = {6,-9,10,11,12,21}; Plane Surface(33) = {32}; Line Loop(34) = {7,3,8,9}; Plane Surface(35) = {34}; Line Loop(36) = {10,-18,16,20,-4,8}; Plane Surface(37) = {36}; Line Loop(38) = {-14,-13,-12,19}; Plane Surface(39) = {38}; // 除了使用文件包含机制以外，我们还能定义函数。在下面的函数中， // 我们使用了保留词 'newp'，这样能够自动的选择一个新的点数。这 // 个数是当前最大数加一得到的。相似的，'newreg' 保留词是所有除 // 去点以外的实体数中最大的加一。