FLAC如何划分网格.ppt

GEN命令的一种最简单的形式可以在建立网格时给网格定义坐标系统，比如，现假定开挖例题的模型在坐标系中的空间范围是这样的：x方向，-10～10；y方向，-10～10；z方向,-20～0。 gen zone brick size 6,8,8 p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) p0 -10, -10, -20 p1 10, -10, -20 p2 -10, 10, -20 p3 -10, -10, 0 gen zone brick size 6,8,8 plot surf GEN(GENERATE)+关键词 fix x range x 5.9 6.1 fix y range y -0.1 0.1 fix y range y 7.9 8.1 fix z range z -0.1 0.1 如果我们想在迭代计算过程中记录一些有用的变量或参数的变化，从而来判断分析是否已经达到平衡状态或者模型已经发生垮塌，就需要使用到HISTORY（HIST）命令 hist n= 5 hist unbal hist gp zdisp 4,4,8 命令的第一行：指定了参数记录的步频：N(Nstep)=5指定了每迭代计算5次记录一次相关的值（如果不指定，则系统默认值为10） 。命令的第二行指定记录最大不平衡力（the maximum unbalanced force），命令的第三行指定记录坐标值为（4，4，8）的节点Z方向的位移。在计算时记录最大不平衡力是一个很好的习惯，如果最大不平衡力接近一个很小的值，并且位移记录值不再发生变化，那就 表明计算已经到达了平衡状态。 set mech force=50 solve 这样，软件便开始了问题的求解，在求解过程中，会在程序窗口中显示计算 的时间步数以及对应的最大不平衡力的值， 当最大不平衡力小于所设定的上限值 （本例题中我们用SET命令设定了最大不平衡力的上限值为50N,即：计算的最大 不平衡力一旦小于此值，我们便人为模型达到了平衡状态。 ）求解过程终止* 。因为我们之前没有关闭绘图窗口，你会发现显示在窗口中图像不断的被更新。如果我们先关闭了绘图窗口，那么求解所需的时间会更短。 计算停止在351步**，可以通过查看最大不平衡力以及所记录节点Z方向的位移采样图来分析，模型是否达到了平衡状态 plot hist 1 Hist 2 最大不平衡力采样图中显示最大不平衡力接近于零，位移图中可以看出位移已经趋于一个固定值。这两个采样结果均表明系统已经达到了平衡状态。 输入的采样记录，系统按输入顺序从“1”开始一直往下编号，这就是为什么输入hist 1命令显示的是最大不平衡力的采样记录图，输入hist 2命令可显示Z方向的位图采样图的原因了。当采样参数较多，自己又不记得输入的顺序时，可以在“FLAC3D”提示符状态下输入：print hist命令来显示所有采样参数对应的编号 接着我们通过下面几行PLOT命令来建立一个比较复杂的新视图，以检查重力场是否已经施加到模型中 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命令，进入“plot”相关的提示符状态，那么命令中所有的plot均可删除。执行以上命令后，便会出现一个名叫GravV的视窗，并且系统指定它为当前显示窗口。我们在命令中已经定义了当前视图的剖面： 倾角为90度 （dip=90） ， 倾向为0度 （dd=0） ，通过点（3，4，0）。在命令中我们还定了视角（plot set rot 15 0 20 ；plot set center （2.5 4.2 4.0）并且使用bound（boundary）命令来增加边界线框，线框只在平面的后面显示（由关键字behind定义，注可简写为BE）。plot add bcont szz plane命令的作用是在当前平面中绘出垂直应力分布图 倾角、倾向：假定了x轴方向为“东”，y轴方向为“北”，z轴方向指“上”。倾角指的是，在x-y平面下半空间区域部分与x-y平面的夹角。倾向指投影在x-y 平面上，以y轴为基准，绕着顺时针方向转，与y轴之间的夹角。指向坐标轴方向的矢量，比如力，位移，速度，只要指向坐标正方向即为正。 现在，最好保存以上得出的初始状态，便于在将来任何时候返回到该状态重新设定相关参数来做相应