新手向——FLAC3D6.0如何输出节点位移

4月前浏览542

在利用FLAC3D计算完成后，有时需要提取单个或多个节点的位移，笔者在此提供常用的2种方法：①使用list功能提取；②编写fish函数提取。

======案例演示======

案例模型为6*6*6*的立方体，分组情况见图1。进行初始平衡后挖去分组‘2’，提取并输出模型各节点的位移。

图1 模型分组

命令流如下：

;建立模型并完成开挖计算model newzone create brick size 6 6 6 group '1'zone cmodel assign m-czone prop bulk 5e8 shear 3e8 coh 1e6 fric 20zone ini dens 2800zone face skinzone face apply v-n 0 range group 'West' or 'East'zone face apply v-n 0 range group 'North' or 'South'zone face apply v-n 0 range group 'Bottom'model gravity 9.81model solve elasticzone group '2' range p-x 1 5 p-z 1 5zone cmodel assign null range group '2'model solve

图2 Z向位移云图

①使用list命令提取位移，相应位移值会输出到软件console界面，见图3.

命令如下：

;提取全部节点位移zone gridpoint list displacement ;提取某个分组（包括面分组）的节点位移，则命令变为：zone gridpoint list displacement range group 'xxx';提取某一个节点(x,y,z)的位移，则命令变为：zone gridpoint list displacement range position (x,y,z);提取某个范围内的节点位移，则命令变为：zone gridpoint list displacement range position-x x1 x2 position-y y1 y2 position-x z1 z2

图3 list位移输出示意

Tips1：输出的位移看Displacement列，不看Small Displacement列，此项可以通过fish函数提取结果对比。

②使用fish函数并提取，本文以提取x向位移并输出到txt文件为例，命令如下：

;提取单个节点(x,y,z)结果并输出到consolefish def disp_out  gpnt = gp.near(x,y,z)  xdisp = gp.disp.x(gpnt)  ydisp = gp.disp.y(gpnt)  zdisp = gp.disp.z(gpnt)  io.out('The x-displacement of gridpoint is: \n' + string(xdisp))  io.out('The y-displacement of gridpoint is: \n' + string(ydisp))  io.out('The z-displacement of gridpoint is: \n' + string(zdisp)) end@disp_out ;提取模型所有节点的编号与x/y/z位移并输出到txt文件，以输出x位移为例：fish def disp_out    n = 0    loop foreach gpnt gp.list        n = n+1        table.x(1,n) = gp.id(gpnt)        table.y(1,n) = gp.disp.x(gpnt)    endloop    command        table '1' export 'xdisp.txt'     endcommandend@disp_out