【FLAC3D】一段找出某个面分组下所有单元的代码

问题来源于二群小伙伴，见下图，本文演示找出面分组为‘Top'的单元，然后将其分组为’test‘。读者可根据需要自行编写其他相关操作的代码。

图1-1 问题来源

2.代码

如果采用python实现，那么需要用到face_in_group()函数，函数用法见下图

import itasca as itit.command("python-reset-state off")it.command("""    model new    zone import 'mesh.f3grid'    zone face skin""")for z in it.zone.list():    for i in range(6):        flag = z.face_in_group(i,'Top')        if flag is True:           z.set_group('test')

如果采用fish实现，则需要用到zone.face.isgroup()函数，函数用法如下：

model newzone import 'mesh.f3grid'zone face skinfish def test    loop foreach z zone.list        loop i(1,6)            flag = zone.face.isgroup(z,i,'Top')if flag = True                zone.group(z) = 'test'            endif        endloop   endloopend@test

图3-1 运行结果