FLAC3D学习笔记(自己总结版)

  • 格式:docx
  • 大小:69.80 KB
  • 文档页数:6

下载文档原格式

  / 6
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

一、建立模型笔记

1键盘上的键的功能:X,Y,Z分别控制所建立的模型围绕X,Y,Z三个轴做逆时针旋转.如果打开大小写锁定键,分边按XYZ时建立的模型会做顺时针旋转.+号的功能是加大模型旋转的角度,-号的功能是减小模型旋转的角度。

2键盘上的上下左右四个键可以实现图形的向上向下,向左和向右的移动.

3当在一个区域建立多个模型的时候,有时候为了镜像单个模块儿,需要预先编组,编组的方法如下:

gen zon tunint p0 0 0 0 p1 0 10 0 p2 5 5 0 p3 0 0 -5 p4 5 10 0 &

p5 5 5 -5 p6 0 10 -5 p7 5 10 -5 dim sq2 sq2 1 sq2 1 1 1 group(编组) z1(组名为:z1) &

fill group tunnel

4镜像整个模型的方法:首先判断即将镜像的模型在原模型的哪个方向(即三个轴的方向),然后在轴所在的方向上去两个点.其中一个是原点(origin(xyz)).

5镜像部分模型的方法: gen zone reflect normal (x,y,z) origin (x,y,z) range group z1就可以了.

6特别要注意的一点:交叉巷道的巷道充填和巷道的网格是两个组,所以在镜像的时候不要忘了给充填部分镜像.

7 group radcyl range group cylint not(编组的时候不包括cylint)

7 快捷键总结

Ctrl+Z 局部缩放的功能。Ctrl+R的功能是恢复到初始状态。

8模型建立的方法:

方法1:利用默认值生成网格,

各个默认值:(1)p0的默认值是(0,0,0),网格的每边的单元格数size默认值是10,网格每边的长度默认值为size的值.

方法2:利用4个点的坐标来生成矩形网格.p0~p3 size

方法3:利用edge来确定边长生成矩形网格.

方法4:利用参数ratio来确定单元体几何边画边生成矩形网格.

方法5:利用参数add(相对坐标)来生成矩形网格.

方法6:利用八个点的左边来生成矩形网格.

9、用户自定义模型的方法

我来贡献一点自己的成果:

FLAC3D的二次开发环境提供了开放的用户接口,在软件安装文件中包含了软件自带所有本构模型的源代码,且给出了Mohr-Coulomb模型和应变软化模型的编译示例,因此可以方便地进行本构模型的修改与开发。为了方便起见,下面的说明以建立UserModel模型为例。

(1) 在模型头文件(usermodel.h)中进行新的本构模型派生类的声明,修改模型的ID(为避免与已有模型冲突,一般要求大于100)、名称和版本,修改派生类的私有成员,主要包括模型的基本参数及程序执行过程中主要的中间变量。( ? p6 u' J5 Q3 y( a

(2) 在程序C++文件(usermodel.cpp)中修改模型结构(UserModel::UserModel(bool bRegister): Constit -utive Model)的定义,这是一个空函数,主要功能是给(1)中定义的所有私有成员赋初值,一般均赋值为0.0。(3) 修改const char **UserModel: roperties()函数,该函数包含了给定模型的参数名称字符串,在FLAC3D的计算命令中需要用到这些字符串进行模型参数赋值。

(4) const char **UserModel::States()函数是单元在计算过程中的状态指示器,可以按照需要进行修改指示器的内容。

(5) 按照派生类中定义的模型参数变量修改double UserModel::GetProperty()和void UserModel:: SetProperty()函数,这两个函数共同完成模型参数的赋值功能。' U; e G' W" Q# R4 q/ @9 G" h

(6) const char * UserModel::Initialize()函数在执行CYCLE命令或大应变模式下对于每个模型单元(zone)调用一次,主要执行参数和状态指示器的初始化,并对派生类声明中定义的私有变量进行赋值。值得注意的是,Initialize()函数调用时没有定义应变分量,但可以调用应力分量,但不能对应力进行修改。9 n# e8 |' c- B/ q, B B5 m

(7) const char * UserModel::Run()是整个模型编制过程中最主要的函数,它对每一个字单元(sub-zone)在每次循环时均进行调用,由应变增量计算得到应力增量,从而获得新的应力。在计算过程中,要根据单元应力情况对单元状态指示器进行赋值。当进行塑性模型编制时,需对达到塑性的应力状态进行修正。2 [/ ~. U( A"

G. b3 @6 ]. `

(8) 修改const char * UserModel::SaveRestore()中的变量,修改方法同(2)和(5),该函数的主要功能是对计算结果进行保存。

(9) 程序的调试有两种方法。①在VC++的工程设置中将FLAC3D软件中的EXE文件路径加入到程序的调试范围中,并将FLAC3D自带的DLL文件加入到附加动态链接库(Additional DLLs)中,然后在Initialize()或Run()函数中设置断点,进行调试;②在程序文件中加入return()语句,这样可以将希望得到的变量值以错误提示的形式在FLAC3D窗口中得到。

10、网格剖面的画法

方法:首先在命令行里输入edit

11、想知道有多少单元体和节点?

方法:用命令print info

12、绘制塑形区域命令是什么?

方法:Plot block state shear ;获得剪切屈服的单元,包括shear-n和shear-p

Plot block state tension-p ;获得过去拉伸屈服的单元

Plot block state now ;获得当前处于塑性状态的单元,包括shear-n 和tension-n

Plot block state past ;获得过去处于塑性状态的单元,包括shear-p 和tension-p

13、step后面的数字是怎么确定的?

答:

14、开挖模型的命令

Model null range group 组名

Model null range group x x1,x2

画圈的正负号是啥意思?

答:正号表示10的正的几次方.而负号表示负的几次方.

11、关于ratio

ratio代表单元体尺寸变化率,关于只需要记住一点:不同模型的两个尺寸相等的面,采用相同的集合变化率,网格才会重合,ratio对于减少单元的数目有很好的效果。

12、地震波的模拟命令流

set dyn on

set dyn multi on

set large

free x y z range z -35.1 -34.9 ;去掉模型底部原有的静力条件

apply nquiet squiet dquiet ran z -35.1 -34.9 ;静态边界条件

table 1 read crvel20.txt

apply szz -1.524293e7 hist table 1 ran z -35.1 -34.9 ;加动力载荷