调整网格简析

1、FLAC3D 5.0培训教程（武汉）,工程师 李振 2014.3.27-3.28 Itasca(武汉)咨询有限公司,FLAC3D 5.0培训日程安排,FLAC3D V5.0界面操作 FLAC3D基本操作方法vs应用流程；FISH调整网格 FLAC3D内置Fish语言的应用； FLAC3D结构单元vs接触单元； FLAC3D渗流模块 其他,2014.3.27 2014.3.28,利用fish调整规整网格建立交叉隧道,一、建模思路 二、建模过程fish语句的分析 三、小结,一、建模思路,1 首先用内置brick生成砖型规整网格 2 在现有网格的基础上调整gridpoint的坐标， 坐标的调整依据建

2、模人员期望的网格形 状来进行 3 规划好网格的调整次序（如先调整端面 网格，再调整中间部分的网格）,镜像两次生成最终网格,调用fish调整网格,二、建模过程fish语句的分析,2.1 为了分析整个网格调整的过程，将手册中的命令进行分解简化，从整个命令来看，最先调整的两个砖型网格依次为brick_3, brick_1（此名称代表两个类组）；为了清晰明了地理解调整过程，因此，在建立brick时，只建立此两个砖型网格，即只保留如下的命令行（前面的参数设置命令不变）。,只分析brick_1和brick_3,2.2 将shape函数进行分解简化，将与brick_1和brick_3无关的fish语句全部删

3、除或处理为注释行。为了弄清楚每个调整步的作用，将每个调整步骤重新定义为一个新的函数，依次运行并在运行后使用plot显示调整后的效果，同时做出深入的分析，弄清每个语句的作用。,Fish语句中，函数retrieve_pointer， 函数set_pointer，函数get_parray， 函数create_and_fill_pointers， 函数find_extremes，保持不变。 将函数create_all_bricks更名为函数 create_13_bricks，只保留与两个brick相关的 语句，如下所示。,将shape_tunnel_b函数分解成如下的单个小的函数， 分别运行这些小函数

4、，分析这些小函数各自的作用。,对块体顶点进行编号,Brick_3即为tunnel-b调整前 的母体网格,分析create_13_bricks这个函数，此函数中的关键部分是create_and_fill_pointers函数,def create_13_bricks ii = out(Indexing Brick 1) pxsize = zone_x_a ; Brick 1 pysize = zone_x_b pzsize = zone_y_tun xmin = x0brick xmax = x1brick ymin = y0brick ymax = y1brick zmin = z1brick

5、 zmax = z2brick xrat = 1.0 yrat = 1.0 zrat = 1.0 create_and_fill_pointers parray01 = parray,ii = out(Indexing Brick 3) pxsize = zone_z_b ; Brick 3 pysize = zone_x_b pzsize = zone_y_tun xmin = x1brick xmax = x2brick ymin = y0brick ymax = y1brick zmin = z1brick zmax = z2brick xrat = zone_rat_b yrat =

6、1.0 zrat = 1.0 create_and_fill_pointers parray03 = parray end create_13_bricks,图中所示建模所用的部分参数， 除brick_3的x方向网格比率为1.2以外， 其余网格变化率皆为1.0， 即定居等分。,x0brick,zone_rat_b=1.2,y1brick,x1brick,x2brick,y0brick,z0brick,z1brick,重点分析create_and_fill_pointers函数,函数get_parray的作用便是计算当前的网格体的gridpoint的数 目，然后利用内置函数get_mem(isi

7、ze)向计算机内存来申请 isize个fish变量对象。 用这些个fish变量对象来存储网格体中 每个gridpoint的地址。以brick_1和brick_3为例，每个块体有 (pxsize+1) * (pysize+1) * (pzsize+1) =6*6*6个gridpoint。 继续执行create_and_fill_pointers函数，下面的语句便是寻找 每个gridpoint的地址，并将其值（pointer类型）传递给相应 的fish变量对象，而这个过程通过函数set_pointer来实现。此 函数中ioffset用来控制这isize个fish变量对象首地址的偏移量。,手册中对内

8、存函数用法的解释,lose_mem(n,ad) This function returns n FISH -variable objects to the host. The parameter ad is the address of the start of the array of objects; there is no checking done to verify that ad is a valid address. The returned value is undefined. get_mem(n)This function gets n FISH -variable obj

9、ects from the hosts memory space and returns the address of the start of the contiguous array of objects mem(ad) This function may be used either as a source or destination (i.e., on the right- or left-hand side of an expression): var = mem(ad) The value of the FISH -variable at address ad is transm

10、itted to the regular FISH variable var (together with the type). mem(ad)=var The value and type of the regular FISH variable var is placed in the FISH -variable at address ad. It is the users responsibility to check that the addresses given contain valid FISH variables. During execution, an error message is given if ad is not of type pointer.,X0 x1 x2 x3 x4 x5,y5 y4 y3 y2 y1 y0,z0 z1 z2 z3 z4 z5,Gridpoint遍历的次序，循环从内到外,另外要注意的就是寻找块体网格节点附近坐标的方法， if xrat # 1.0 then fxpos = (xratixpos)-1.0) / (xratpxsize)-1.0) else fxpos = float(ixpos) / float(pxs