一、前沿

随着分析的深入，对我们的理解越加具有挑战性。本节介绍的概念比较抽象，构建的模型内部带有体积孔。但是我们在创建基本模型的时候，其思想并没有改变（点-线-线环-面-面环-体）。接下来我们将逐步解析本节将要介绍的代码功能。

【注】使用的软件Gmsh

二、代码（来自于Gmsh官方文档“t5.geo”）

三、解释

本节代码分四部分阐述：

（1）第一部分

这一部分比较简单，构建缺角立方体形状。此处不再赘述。

（2）第二部分

这一部分代码，可以理解为定义了一个函数，函数的作用就是创建小球。其中的细节，现解释如下：

a） 人为定义标签，不可能不重复的定义；所以，我们使用保留变量名，系统自动定义标签

newp：它自动选择一个新的点标签

newc:  它自动选择一个新的曲线标签

newcl: 它自动选择一个新的曲线循环标签

news: 它自动选择一个新的曲面标签

newsl: 它自动选择一个新的曲面循环标签

newv:  它自动选择一个新的体积标签

（3）第三部分

这一部分，可以理解为for循环结构；循环画出5个球体；其中需要传入参数（函数中定义的参数）；

【信息打印】

在终端上打印一些变量（注意，因为所有“.geo”文件中的变量在内部被视为浮点数，格式字符串应该只包含有效的浮点格式，诸如‘%g’、‘%f’、‘%e’等说明符）：

Printf("Hole %g (center = {%g,%g,%g}, radius = %g) has number %g!", t, x, y, z, l, thehole);

（4）第四部分

将整体加上外表面，构成几何体。

四、网格划分

五、总结

本节，在构建模型方面，使用了宏，循环等定义，简化了建模过程，在构建复杂模型的道路上更加前进了一步。

参考文献

[1]. 有限元分析（1）

[2]. 有限元分析（2）

[3]. 有限元分析（3）

[4]. 有限元分析（4）

[5]. Gmsh官方教程"t5.geo"