gmsh生成非结构网格详解.docVIP

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