在FLAC3D6.0中使用fish语言批量处理结果文件

8月前浏览792

有时在计算完成之后，需要对多个结果文件进行重复的操作。若是直接在fish函数中调用结果文件，则会出现以下报错：

图1 fish调用结果文件报错

这是因为flac3d在执行model new或是model restore操作时，会清空当前的所有fish函数，从而出现错误，因此在fish函数中调用restore就可以认为是我杀死了我。

图2 是我杀了我

同时，我们知道FLAC3D是一款命令流驱动的软件，即软件所执行的所有操作都依靠发送给底层的命令流实现。因此，对于不同结果文件间的重复操作，我们虽不能直接使用fish函数去执行，但我们可以通过fish函数生成重复的调用命令，从而实现批量处理结果文件。看下面例子：

; 建立模型并进行初始平衡model newzone create brick size 6 6 6 group 'a'zone group 'b' range p-x 2 4 position-z 2 4zone splitzone cmodel assign elasticzone prop bulk 3e8 shear 1e8 dens 1800zone face skinzone face apply v-n 0 range group 'West' or 'East'zone face apply v-n 0 range group  'South' or 'North'zone face apply v-n 0 range group 'Bottom'model gravity 9.81model solve; 开挖6步并存储相应的结果文件fish def _kw    loop i(0,5)        j = i+1        filename = 'kw_' + string(j) + '.f3sav'        command            zone cmodel assign null range  group 'b' p-y @i @j            model solve            model save @filename        endcommand    endloopend@_kw; 使用fish语言实现结果文件的自动调用fish def _postProcess    array arr(1)    msg = ''    loop i(1,6)        filename = 'kw_' + string(i) + '.f3sav'        msg += "model restore " + "'" + filename + "'\n"        msg += "zone gridpoint list displacement range  position  3 3 6\n\n"    endloop    arr(1) = msg    oo = file.open('list.txt',1,1)    oo = file.write(arr,1)    oo = file.close()end@_postProcess; 调用生成的命令流文件call 'list.txt'

图3 模型图

图4 生成的命令流文件

======命令解释======

; 使用fish语言实现结果文件的自动调用fish def _postProcess    array arr(1)    msg = ''    loop i(1,6)        filename = 'kw_' + string(i) + '.f3sav'        msg += "model restore " + "'" + filename + "'\n"        msg += "zone gridpoint list displacement range position  3 3 6\n\n"    endloop    arr(1) = msg    oo = file.open('list.txt',1,1)    oo = file.write(arr,1)    oo = file.close()end@_postProcess

以上是生成重复命令流的代码，从中可以看出，使用fish批量后处理的流程为：首先通过将需要执行的命令流写作字符串msg，然后将字符串msg写入数组arr中，最后将数组arr写入文本文件list中。

在这个过程中，需要注意以下几点：

1、需要进行后处理的结果文件，在命名时应具有一定规律，参考本例的命名方式。

2、想要在字符串中加入引号时，推荐利用双引号""将单引号'括起来，见第7行的"'"。

3、+=为复合赋值运算符，其应用举例: a = a + 1 等效于 a += 1。

4、'\n'表示换一行，'\n\n'即为换两行。

5、以上案例只是简单示范，读者应根据实际需要自行构造msg的内容。

======总结======

使用fish语言进行计算结果间的后处理存在很大局限性，推荐使用python进行后处理操作，尤其是FLAC3D7.0中的python可以安装pandas等第三方库，极大地便利了批量后处理，相关内容择日发布。

来源：FLAC3D小技巧

著作权归作者所有，欢迎分享，未经许可，不得转载

首次发布时间：2023-12-26