gmsh生成非结构网格详解(之一).doc

Gmsh 生成非结构网格详解（之一） 在流体计算中，我们经常要生成二维非结构网格。利用Gmsh软件可以快速生成，我们下面可以看一个例子，例子在gmsh的tutorial目录中： 下面是 t4.geo的完整说明： cm = 1e-02 ; e1 = 4.5*cm ; e2 = 6*cm / 2 ; e3 = 5*cm / 2 ; h1 = 5*cm ; h2 = 10*cm ; h3 = 5*cm ; h4 = 2*cm ; h5 = 4.5*cm ; R1 = 1*cm ; R2 = 1.5*cm ; r = 1*cm ; ccos = ( -h5*R1 + e2 * Hypot(h5,Hypot(e2,R1)) ) / (h5^2 + e2^2) ; ssin = Sqrt(1-ccos^2) ; Lc1 = 0.01 ; Lc2 = 0.003 ; // 下面是所有的操作符的列表，除了 ^ 以外，这些操作符和 C、C++ // 语言中都是一样的： // // - (in both unary and binary versions, i.e. as in -1 and 1-2) // ! (the negation) // + // * // / // % (the rest of the integer division) // // // = // = // == // != // (and) // || (or) // || (or) // ^ (power) // ? : (the ternary operator) // // 使用括号可以将表达式结合 // // 除了这些操作符以外，所有的 C 数学函数也能够使用（注意首字母 // 要大写）： // // Exp(x) // Log(x) // Log10(x) // Sqrt(x) // Sin(x) // Asin(x) // Cos(x) // Acos(x) // Tan(x) // Atan(x) // Atan2(x,y) // Sinh(x) // Cosh(x) // Tanh(x) // Fabs(x) // Floor(x) // Ceil(x) // Fmod(x,y) // // 我们还有增加的函数： // // Hypot(x,y) computes Sqrt(x^2+y^2) // Rand(x) generates a random number in [0,x] // // Gmsh 中唯一预定义的常数就是 Pi. Point(1) = { -e1-e2, 0.0 , 0.0 , Lc1}; Point(2) = { -e1-e2, h1 , 0.0 , Lc1}; Point(3) = { -e3-r , h1 , 0.0 , Lc2}; Point(4) = { -e3-r , h1+r , 0.0 , Lc2}; Point(5) = { -e3 , h1+r , 0.0 , Lc2}; Point(6) = { -e3 , h1+h2, 0.0 , Lc1}; Point(7) = { e3 , h1+h2, 0.0 , Lc1}; Point(8) = { e3 , h1+r , 0.0 , Lc2}; Point(9) = { e3+r , h1+r , 0.0 , Lc2}; Point(10)= { e3+r , h1 , 0.0 , Lc2}; Point(11)= { e1+e2, h1 , 0.0 , Lc1}; Point(12)= { e1+e2, 0.0 , 0.0 , Lc1}; Point(13)= { e2 , 0.0 , 0.0 , Lc1}; Point(14)= { R1 / ssin , h5+R1*ccos, 0.0 , Lc2}; Point(15)= { 0.0 , h5 , 0.0 , Lc2}; Point(16)= { -R1 / ssin , h5+R1*ccos, 0.0 , Lc2}; Point(17)= { -e2 , 0.0 , 0.0 , Lc1}; Point(18)= { -R2 , h1+h3 , 0.0 , Lc2}; Point(19)= { -R2 , h1+h3+h4, 0.0 , Lc2}; Point(20)= { 0.0 , h1+h3+h4