plot fip outside on 显示力的箭头
建模完事后就要给模型赋值，赋参数了，关键词，model model mor（赋摩尔库伦模型，想要什么模型改相应的词）
pro den 2500（密度）bulk2e8体积模量，shear 2e8 （剪切模量）co 2e6（内聚力）fric45（摩擦角）ten1e6 抗拉强度。range...........(我们要赋值的范围) 给模型赋值就是hi这么简单，难在参数的选取。
现在看边界，如果是不动面，用fix x range。。。。。如果边界带速率，或位移用ini。。。 set grav 0 0，-9.81 ini dens 1000
设置重力就一定要设置密度
不考虑重力就不要设计密度。两者共存的，设计程序的时候一定要记的 set mech force 50 solve
用SET命令设定了最大不平衡力的上限值为50N,即：计算的最大不平衡力一旦小于此值，我们便人为模型达到了平衡状态。求解过程终止。设置这个目的是在一定精度要求下，加快计算（值越大，算的量越少，花的时间也越少）
说到这里，add与set的区别 add原来有的加新的用add set原来没有的，设置用set plot create GravV
plot set plane dip=90 dd=0 origin=3,4,0 plot set rot 15 0 20 plot set center 2.5 4.2 4.0 plot add bound behind plot add bcont szz plane plot add axes
plot show 。如这个命令，显示图像的时候命令有很多， 那么我们想显示别的图像，又不想吧这个弄好的图像丢掉 则用plot create name建立一个图像，在建立新的图像时 plot create name1（新的名） 显示以前的图像
那么输入如下命令：plot current name。另一种方法，输入：plot show name 这里有一点最好开始的时候将网格剖分的相对稀疏一些（1000～1500个单元 左右），先求解一下，看看结果，然后再根据大致的结果有针对性的增加模型的 网格数，以得到更精确的解答。
更快的检验一下自己 的模型是否建错了 不同的区域定义不同的本构模型 model mohr range z=0,5 prop ClayeyGravel range z=0,5 在定义参数后面加个range范围
为了避免在输入命令时总要通过range＋全局坐标域的方法来限定该命令执 行的区域，我们可以先通过range命令来给不同的区域命名，然后就可以用range ＋区域名的方法来限定范围了。 这样就避免了每次输入命令要限定在该区域执行 时，都要输入一遍坐标范围的麻烦。具体如何给区域命名，以及日后如何用区域 名来限定该区域，可以看下面这个例子： range name trench x=-1,1 y=0,4 z=-2,2 model null range trench
今天做方案，没时间，就讲几点得了，1：右手法则，建立模型时，的点的排放时有规律的， 方法如下，拿出右手 伸直了
4指合并指向就是p0到p1的方向 掌心指向就是p0到p2的方向
拇指与4指垂直的方向就是p0到p3的方向
很多人说网格后的云图有带齿什么的，改进的方法就是重复用下面的命令 Gen ADJust
用此命令来微调网格, 使之离散平滑。该命令可以连续使用以增进平滑效果的。
plo blo gro可以看到不同的group的颜色分布（例不同的土层定义组后，就可以用此命令显示不同的色了） 如何查看剖面上任一点的位移值？ print gp disp range x() y( ) z( )或 print gp disp range id
plot gp disp range x() y( ) z( )或 plot gp disp range id
很多人都是这么做的，我也是学习， initial 与 apply 有何区别？
initial初始化命令，如初始化计算体的应力状态等； apply边界条件限制命令，如施加边界的力、位移等约束等。
initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化
如果想看模型在计算过程中的变化情况
Set movie avi step 1 file 1.avi；定义动画，动画记录的频率为1时步一频率（计算量大的时候就没必要1时步的记录，可以10,100时步的记录，减少电脑的使用） Movie start；动画开始
hist n= 5 ：N(Nstep)=5指定了每迭代计算5次 记录一次相关的值（如果不指定，则系统默认值为10） hist unbal 命令的第二行指定记录
最大不平衡力（the maximum unbalanced force） hist gp zdisp 4,4,8 记录坐标值为
（4，4，8）的节点Z方向的位移（这个用的多，查看收敛量时就是这个命令，上面两个个人人认为用处不大，也很少用过）
如何在程序中进行暂停，并可恢复计算？
答：在命令中加入pause命令，用continue进行继续 个人认为这个命令非常好
用在我们分步求解中想得到某一个过程中的结果，不用等到全求完，还可以在分布求解错误的时候就进行改正。而不是等到结果出来 关于dd、dip的理解：
不同人有不同的理解，我的理解就是dip面与水平面的夹角 dd是面与水平面的交线的垂线与y轴的夹角
1) 在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负);5 L0 D# _1 ~1 C- L 2) \表示对其值作\如果大家使用AUTOCAD的话就明白\的含义了,比如上面的 hist 1 gp ydisp的值是0~100,那么vs -1就变成-100~0.以此类推.
3) 为什么要这么做呢?主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内.
最近看有很多对后处理中的收敛有很多不同的想法， 自己也总结一下
plot sk dis scale 0.07显示位移矢量箭头的后面的scale默认的是0.05 3越大箭头越大
magf也是放大变形的意思，但是不能放在sk dis的后面，计算机会出错
回到主题，收敛数据，我们想得到一条收敛量随时间变化，或者不同应力直间关系的图线 方法
1用set log on然后把数据导入excel或者其他的处理软件进行处理
2.用hist id=1（不写id的话，电脑会自己按123456的顺序进行排序）gp（点，如果是体的话 用zone）zdis xyz（xyz点处z向位移）
这里应力应变之类的都是用zone，位移之类的都是用gp 最后相要出相应的结果用命令 plot hist 1；
如果相要的不同的id之间的关系（比如说竖向应力szz与剪应力sxz的关系）命令 plot hist 2 v 3 （2,3是对应应力开始设的id号）
flac3d程序里面自动设置的步长是10，也就是说每隔10步记录一次，这样如果模型小得出的结果就数据少图线不光滑，如果模型大，得到的数据多计算文件又太大，浪费空间。设置步长的命令 set hist_rep n(n就是步长) 上面说的plot多是flac里显示的 如果想把结果存到别处去 用命令
hist write 2 v 3 file namehist.txt(其实2相当于x轴，3相当于y轴) 设置影片格式上面见过，温故而知新
set movie avi step 1 （存储步长，可以任意设，就是每隔几步存储的意思）file name.avi movie start
分步开挖就是分布计算
每多一步（多添一块土，或多开挖，多掘进等）就计算一下，加一个solve即可
如果想保存每一步的结果和图，那就在每一步的solve后面添加保存命令设立相应的文件，以后调用即可 1.设定初始地应力场是十分必要的 没有地应力场进行计算会使模型在重力作用下发生过大的变形 从而不能与实际情况相符
2。apply是施加的外力 ini是初始的内力。
3. 有人问如果改参数这样做的话地应力不是被人为改变了吗？为什么要人为加大它的参数呢，很不合理？ 回答：不能改变bulk和shear模量，因为这两个参数是用来计算应力的。而con和fric是用来判断破坏状态的。如果将bulk和shear模量改变了，那你在重力条件下计算的应力就不真实了。咱们改大参数的目的是让模型在不破坏的条件下得到正确的应力，所以只有将con和fric变大才能达到目的。
当然如果在弹性模型下计算初始应力就没有这种问题了。因为弹性模型是没有破坏准侧的当然也没有描述达到破坏状态的参数了。
4. 大家正好讨论地应力的问题，我也有个疑问，平衡地应力后，水平方向的地应力如何确定的？是由本构模型来确定的吗？
这个问题现在我也没答案也是我的问题
5. gen zone brick size 10,10,10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohr
prop bulk 5e9 shear 3e9 fric 35 ini density 2500 set gravity 0,0,-10 fix x y z range z -0.1 0.1 ini szz = -5.0e6 grad 0,0,2.5e4 ini sxx = -2.5e6 grad 0,0,1.25e4 ini syy = -2.5e6 grad 0,0,1.25e apply szz = -4.5e6 range z 19.9 20.1
