gen zone brick size 6 8 8 ;建立有限差分网格，X方向6个，Y方向8个，Z方向8个
plot ;进入绘图模式
creat Trench ;创建一个叫“Treach”的视图窗口
add surface yellow ;是表面显示黄色
add axes black ;使坐标轴显示黑色
show ;在屏幕上显示
;/t ;执行回车命令(我想在这一步让程序自动执行回车命令)
model mohr ;定义本构模型为摩尔-库伦本构模型
;print mem sav mem sys ;查看此问题需要占用的内存
;定义材料的各项参数
prop bulk=1e8 ;体积模量=1e8
prop shear=0.3e8 ;剪切模量=0.3e8
prop fric=35 ;内摩擦角=35度 注：以上三个还可以写成：prop bulk=1e8 shear=0.3e8 fric=35
prop coh=1e10 tens=1e10 ;定义粘聚力=1e10，抗拉强度=1e10
set grav 0,0,-9.81 ;加载重力，单位m/s2
ini dens=1000 ;定义初始网格质量密度(后面会显示：384 zones modified)
;定义模型边界条件
fix x rang x -0.1 0.1 ;固定边界面X=0边界面上所有X方向节点的位移
fix x rang x 5.9 6.1 ;固定边界面X=6边界面上所有X方向节点的位移
fix y rang y -0.1 0.1 ;固定边界面Y=0边界面上所有Y方向节点的位移
fix y rang y 7.9 8.1 ;固定边界面Y=8边界面上所有Y方向节点的位移
fix z rang z -0.1 0.1 ;固定边界面Z=0方向上所有Z方向节点的位移
;记录某些参数的变化
hist n=5 ;指定参数的记录频数
hist unbal ;记录最大不平衡力
hist gp zdisp 4,4,8 ;记录坐标值为（4，4，8）节点在Z方向的位移(输入后会显示：Displacement Mag.Gp(537)(4.0000e+000,4.0000e+000,8.0000e+000))
set mech force=50 ;设定最大不平衡力
solve ;按一定的精度自动求解平衡状态
;显示最大不平衡力以及节点（4,4,8）在Z方向上的位移图
;plot ;进入自动绘图模式
;hist 1 ;显示最大不平衡力变化曲线
;hist 2 ;显示节点（4,4,8）在Z方向的位移变化曲线
;注：当采样的参数较多时，并且自己又记不住输入的次序，可以在“FLAC3D>”提示符下输入“print hist”
;检查应力场是否已经加到模型中
plot
creat GravV ;创建一个名叫GravV的视图窗口
set plane dip=90 dd=0 origin=3,4,0 ;设置一平面，倾角为90度，倾向为0度，经过节点（3,4,0）
;设定视角
plot ;此处出现bug，因为不能利用上面的plot直接调用"FLAC3D>"
set rot 15 0 20
set center 2.5 4.2 4.0
add boundary behind ;增加边界线框，并且线框只在平面的后面显示
add bcont szz plane ;在当前面绘制出垂直应力分布图
add axes ;添加坐标系
show ;显示图形
;到此可以利用save Trench.sav命令来保存一个名叫Trench的.sav类型的文件
;-----------------------------------------------------------------------开挖阶段-----------------------------------------------------------------------
prop coh=1e3 tens=1e3 ;重新定义粘聚力和抗拉强度，之所以取如此小的值是为了让土体发生
垮塌
model null range x=2,4 y=2,6 z=5,10 ;开挖部分（通过设置模型的材料为null）通过限定xyz三个方向的区域来实现
set large ;垮塌属于大变形，将计算设置成大变形
ini xdis=0 ydis=0 zdis=0 ;将所有网格点位移清零（因为需要模拟的是开挖过程引起的位移的变化，而不是从加载重力开始到开挖整个过程的位移变化）
step 2000 ;执行2000步时间
;查看位移等值曲线图
plot ;进入绘图设置
creat DispCont ;创建名为DispCont的视图窗口
copy GravV DispCont settings
add cont disp plane behind shad on
add axes
show