Flac3D基础知识
1. Generate 命令产生网格,简写为 GE。
2. zone 参数代表了三维网格体(单元体),简写为 zo。
3. (1)size 参数指示单元体各边小单元体的个数; (2)dimension 参数指示单元体内部巷道各边边长,如果没给出就默认为两点间距的20%;(3) edge 参数指示网格体各边的尺寸 Ge zone brick p0(a,b,c) edge x y z 边长为x y z 顶点为p0;(边长也可以用确定点坐标来利用点与点的间距来确定p0 p1 p2……)。
4.各种网格体名称及关键字:矩形网格(Brick),退化矩形网格(dbrick),楔形网格(wedge),金字塔网格(pyramid),四面体形网格(tetrahedron),圆柱体网格(cylinder),矩形体外环绕放射状网格(radbrick),平行六面体外环绕放射状网格(radtunnel),圆柱体外环绕放射状网格(radcylinder),柱形壳体网格(Cshell),交叉圆柱体网格(cylint),交叉平行六面体网格(tunint)。
5.绘制网格体语句:
①Generate zone brick size a,b,c (其中的brick可以被其他的网格体替换)。
②Generate zone brick p0=(,,) &
P1=(,,) &
P2=(,,) & (边长由点点间距确定)
P3=(,,) & → ( &符号为继续行的符号)
Size a b c → (各边小单元体个数)
③如果在p1 p2 p3 前面加 add 则各点坐标是相对于p0点的坐标。
④如果只给出p0 则由 edge a & 给出边长。(见第 3 条);利用ratio a,b,c 给出小单元体长度比例依次变化。
⑤内部有空单元体或者巷道的网格体,如果内部巷道需要填充为单元体,使用 fill group groupname 语句。(其中groupname为巷道内填充单元体的组名)。
皇帝诏曰
⑥显示图形:plot t rotation a b c 使图形可以围绕一点旋转
Plot magnification a 设置图形比例大小
Plot add axes red 设置坐标系及其颜色
Plot show 显示图形
Plot add surface red 设置图形表面颜色,设置坐标轴某一范围
的表面颜色时利用range参数(课本39页4-1),add surface red range z= a b 或者 课本42页的 plot block group yellow red
Plot blo gro 自动为各组分层,用不同颜色显示
马丁路德改革
⑦Generate zone brick p0=(,,) size a b 1 一层陪衬的底网格体。
↗第4条中的网格体名称(课本37页)
⑧Generate zone后面的关键字
↘copy复制,reflect镜像两个
⑨reflect镜像参数、copy复制参数
(1)Genertae zone reflect dip # dd # origin =(,,)
黑体部分表示(dip 与xy平面的夹角即倾角);(dd 镜像面的倾向,倾向线与y轴的夹角);38页手绘。Dip走向—正Y轴开始顺时针转动,dd倾向—xy平面顺时针转动为正。
描写雪景的古诗
(2)Generate zone reflect normal x y z origin=(,,)以normal和origin所示点连线作为镜像面的法线,找到镜像面。46页例题中以(-1,0,0康乃馨的寓意)即X轴为法线,找到镜像面作镜像的。
(3)Generate zone鲁滨逊全名 copy x y z range # 距离中心点(x,y,z)向量的距离复制规定范围内的实体。
6.Fix命令参数
Fix x range x= a b 表示a、b两个平面间的所有网格节点指定参数值不变。
7.求解,保存,等值线(应力、位移),range参数的使用。(见22)√
8.pau 暂停 continue继续 plot blo gro 显示分组
9.t magnification a 设置比例大小
10.fill group inner 填充内部网格体(见文件夹2)
12. 定义小组
group 小组1
group 小组2 range 小组1 not
13.shift +M 缩小 M 放大
14.plot 与语句 之间不需要用 & 继续符号
15.矩形体外环绕放射状网格 p0 p1 p2 p3 位置的确定 (见文件夹2)
16.Ctrl + R 恢复原来状态
17.一个单元只能有一种颜色 ,定义颜色时过半就能显示出颜色,不过半不显示。(见文件夹4)
18.对于一个界面,只能定义一个坐标轴。Size,dimension标注见图。
(这个在文件夹6,有体现出来,在定义坐标时)
19. 请问怎么把hist记录的数据全部输出到txt文件那? his write 1v2 file 文件名.txt
20.考虑孔隙水压力:ini pp XXX
21.记录采样(x,y,z)处某方向的节点位移: hist gp ydisp (x,y,z)
记录采样(x,y,z)节处某方向的单元体应力: hist zone syy (x,y,z)
22 hist 1 位移
Hist 2应力
输出2和1的hist关系图用:plot his -2 vs 1 如果应力值都为负数,则加-2是为了应力坐标轴为正!
23 摩尔库仑模型压缩测试
应变硬化/软化模型软化测试: prop XXXtab
Table 1 赋值含义 给谁赋值 ??没搞懂 你们研究下吧 给好评哈·
24 材料参数附加关键字
Gauss_dev A 输入值为平均数 A为偏差 随机进行高斯正态分布计算
Gradient gx gy gz 值=输入值+gx*X+gy*Y+gz*Z
Uniform_dev A 输入值为平均数 A为偏差数 随机进行均匀分布计算
25.创建桩结构单元:l pile id=1 begin x y z end x y z ng M 该桩整体包括M个桩构件
设置桩相关参数命令:l pile pro ……
26.Apply命令对任何边界或者部分边界施加力。
Apply szz =-1e5 sxz=-.5e5 range z=(-.1 , .1) 在范围-0.1~0.1内对模型边界施加应力分量。-----nstress(法向力) dstress(局部面倾向的切应力) sstress(局部面走向的切应力)
Apply nstress -1e6 range plane dip # dd# origin x y z 对过点x y z 走向逆反心理教育 倾向的平面施加-1e6的法向应力---肥胖纹严重图片见文件夹示例 plot surface red fap yellow 红色表面 、黄色应力线。
渐变应力施加—Apply szz -1e6 gradient 0,0,1e5 range z -1000,0 按照公式S=S0——gx*x+gy*y+gz*z 进行应力渐变
27. 记录网格节点速度历时曲线:
(1)跟踪记录网格中一些关键点的速度(命令:HIST gp vel <x y z>)并绘制其历时曲线。
(2)采用PLOT vel命令绘制完整的速度矢量场图。
在命令app nstress ……与solve间插入下面的命令:
hist unbal ;记录计算过程中的体系最大不平衡力
plot hist 1 ;绘制体系最大不平衡力历时曲线 四眼鱼
根据上述命令即可观察计算过程中体系最大不平衡力随计算时步的演化过程,,体系最大不平衡力随着计算的进行,逐渐趋近于0,表明体系最终达到了力平衡状态。不过,力平衡状态仅表示所有网格节点的合力为0,并非表明体系处于真实的物理平衡状态,因为在力平衡状态下,体系也有可能正在发生稳定的塑性流动。这时,就需要借助其它方式如观察网格节点速度、塑性区标识等来进一步评估模型所处的状态。
29.plot bound black 显示网格结构
30.对模型进行部分挖空---model null range x a,b y a,b z a,b
plot t plane dip 90 dd 180 origin 0,2.5,0 设置前视平面
30.显示体积模量(bulk)参数赋值情况。plot block prop bulk
31. print gp position range id=14647 输出节点坐标
32. 有些外部显示的命令 都可以在 点击图形——plotitems——modify里面设置。