Flac3D命令--完整经典版.doc

实例分析命令： X，Y，Z 旋转 Shift+ X，Y，Z反向旋转 Gen zone……；model……；prop……（材料参数）；set grav 0,0,-9.81（重力加速度） plot add block group red yellow把在group 中的部分染成红色和黄色 plot add axes black坐标轴线为黑色；print zone stress% K单元应力结果输出 ini dens 2000 ran z a b（设置初始密度，有时不同层密度不同）；ini……（设置初始条件）；fix……（固定界面） set plot jpg ；set plot quality 100 ；plot hard file 1.jpg图像输出（格式、像素、名称） plot set magf 1.0视图的放大倍数为1.0；plo con szz z方向应力云图 ini z add -1 range group one群one的所有单元，在z方向上向下移动1m；然后合并命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体，1e-5是容差 (基坑开挖步骤)：Step 1: create initial model state（建立初始模型）Step 2: excavate trench（开挖隧道） group Top range group Base not 定义（群组Base以外的为）群组Top plot blo gro使得各个群组不同颜色显示 （两个部分间设置界面；切割法）：gen separate Top使两部分的接触网格分离为两部分；interface 1 wrap Base Top在（Base和Top）这两部分之间添加接触单元；plot create view_int显示，并创建标题view_int；plot add surface显示表面；plot add interface red界面颜色红色 （简单的定义函数及运行函数）new；def setup定义函数setup；numy = 8定义常量numy 为8；depth = 10.0 定义depth为10；end结束对函数的定义；setup运行函数setup （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。 模拟模型的材料问题时为什么要去定义某个方向上的初始速度？— 渐变应力施加：apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60 dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 0 dip dd确定平面位置使用：（纠结） print gp position range id=14647 输出节点坐标 apply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0在这个求解方程中，z为变量，所以为： ；原点（0,0,0） free x range x -.1 .1 z 6.9 10.1放松x=0 平面上，z=7，10 这一部分在x 方向的约束（可以在此处产生破坏） 体积模量K和剪切模量G与杨氏模量及泊松比v之间的转换关系如下： 一般而言，大多数问题可以采用FLAC3D默认的收敛标准（或称相对收敛标准），即当体系最大不平衡力与典型内力的比率R小于定值10-5；（也可由用户自定义该值，命令：SET mech ratio value） 所谓体系最大不平衡力，是指每一个计算循环(或称计算时步)中，外力通过网格节点传递分配到体系各节点时，所有节点的外力与内力之差中的最大值； 所谓典型内力，则是指计算模型所有网格点力的平均值。 要遵循由简单到复杂、由少网格到多网格的思路。 混合离散化法的基本原理是通过适当调整四面体应变率张量中的第一不变量，来给予单元更多体积变形方面的灵活性。 gen zone ref以z=0平面为对称面镜像生成网格 ; gen zone ref gen zone ref dip 90以y=0面为对称面镜像生成网格 ; gen zone ref dip 90 dd 90以x=0平面为对称面镜像生成网格 球体模型建立：make_sphere函数循环使用。计算球心至节点距离、获取节点坐标、调整坐标值得到球面至球心距离、完成八分之一球模型。变换节点内存地址！ ①变形云图：（位移、应力、速度等等）；plot add cont szz out on shade on打开szz的等值线图（outline、effective网格、有效的…显示）；PLOT con szz out