FLAC3D建模方法探讨
来源:用户上传
作者:
摘要:针对FLAC3D前处理方面的功能不足,本文列举了目前FLAC3D建模的常用方法以及其建模思路,供需要者参考和选择,具有一定的实用价值。
关键词:FLAC3D;建模方法
1.引言
随着计算机技术的发展,数值模拟方法已广泛应用于岩土及地下工程的研究和设计中,FLAC3D是三维快速拉格朗日法(Fast Lagrangian Analysis of Continua-3D)的缩写,它是一种基于三维显式有限差分法的数值分析方法,可以模拟岩土或其他材料的三维力学行为[1]。这种分析方法非常适合于模拟大变形问题,尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点,但FLAC3D非完全可视化的建模方式和网格划分方法过于繁琐,限制了其通用性[2]。
2.建模方法探讨
随着FLAC3D的应用逐步普遍,出现了FLAC3D的一些建模的方法,在此列举出来,以供需要者参考和选择。
2.1 FLAC3D中INTERFACE建模方法
INTERFACE单元在数字模拟中应用相当广泛,很多数值模拟软件都包含有这种单元。它可以模拟结构的接缝,不同材料、不同块体的接触面,以及岩石的节理面等。
在FLAC3D中使用INTERFACE建模,目前来说基本上有三种方法,分离法、移动组合法和依次建模法。
分离法的思路是先把块体单元分组建好,然后将界面上两部分的接触网格分离为两部分,最后设置INTERFACE单元。
移动组合法的实现步骤是先在不同的位置分别建立两组块体,然后在一个块体的表面建立INTERFACE,最后通过移动使两组块体在INTERFACE处组合起来。
依次建模法是建立一组块体,然后在其相应表面建立INTERFACE,再建立与之接触的另外一组块体。对于复杂的模型,交替建立块体和INTERFACE[3]。
2.2 基于Midas/Gts的FLAC3D的建模方法
在Midas/Gts中建立复杂地质模型,可以应用Midas/Gts内嵌工具地形数据生成器(TGM),或者从AutoCAD导人线段建模等方法建立模型,图形建模法可以实现“所见即所得”,强大的布尔运算功能及网格划分功能也能使使用者很方便快捷地建立所需要的功能.模型建立完毕,使用菜单导出单元节点信息及节点坐标信息.通过FISH语言,将节点信息和单元信息转换为FLAC3D的Command格式。最后应用FLAC3D内嵌FISH语言,将由Midas/Gts建立的模型导入FLAC3D平台,进行赋值及后处理[4]。
2.3 基于ANSYS软件的FLAC3D建模方法
利用已有的建模与网格划分功能强大的ANSYS软件对复杂岩土工程建立相应数值模型(包括网格划分),根据ANSYS初始单元体和FLAC3D初始单元体的关系以及ANSYS的节点文件(NLIST.LIS)和单元文件(ELIST.LIS)的特点,再通过优化算法将ANSYS数据转化,直接生成FLAC3D可以识别的单元节点格式文件(*.FLAC3D)[5]。
2.4 基于Surfer平台的FLAC3D建模方法
AutoCAD生成的地表、基岩等高线(已赋高程)及地形地质图里则包含了地层界面、地质构造的三维坐标信息,用件将这些三维坐标信息提取后输入Surfer软件即可生成经过插值后的储存了地层界面信息的*.grid文件,此后将这些经过插值的、网格对应的*.grid文件另存为*.dat文件即生成网格点对应的各地层界面控制点三维坐标,然后用FLAC3D里面内嵌的Fish语言编辑的简单指令即可生成FLAC3D可以读取的*.dat三维地质模型文件,然后用Call命令读取这些*.dat文件即可生成复杂三维FLAC3D数值模型[6]。
2.5 应用MATLAB接口程序建立FLAC3D模型方法
用一定大小(可根据计算的精确度调整)的方格网覆盖于地形图之上,读出网格点的坐标 (x、y、z),即拟建山体模型的山顶控制点坐标。将其以一定格式记入数据文件中,如在*.txt文件中输入3列数据,形如“x、y、z”表示控制点坐标,以便MATLAB读取数据文件将其中包含的数据赋值给指定变量。在MATLAB环境中,建立一个能包含所有单元节点坐标的数组S(x,y,z),通过循环和MATLAB的打印命令生成FISH格式的脚本文件。
在建立较为复杂模型时,可以将该模型分成几块,编写不同的脚本文件,在FLAC3D中用CALL命令依次调用,完成建模过程[7]。
3.结论
针对FLAC3D前处理方面的功能不足,列举了目前比较常用的FLAC3D前期建模方法,供需要者参考和选择,具有一定的实用价值。
转载注明来源:https://www.xzbu.com/2/view-585138.htm