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

gen zone ref以z=0平面为对称面镜像生成网格 ; gen zone ref gen zone ref dip 90以y=0面为对称面镜像生成网格 ; gen zone ref dip 90 dd 90以x=0平面为对称面镜像生成网格 球体模型建立 ①变形云图：〔位移、应力、速度等等〕；plot add cont szz out on shade on翻开szz的等值线图〔outline、effective网格、有效的…显示〕；PLOT con szz out on ；PLOT con zd ou on magf 20〔块云图bcon不能跟magf放大倍数〕；plo con zdisp outline on shade on 翻开zdisp的等值线云图 ②变形矢量图：〔变形方向、大小PLOT sk dis scale 0.07其中scale用来改变箭头大小，默认是0.05〕； ③塑性区分布print zone state都可以用来显示模型的塑性区 应力正位于屈服面，或说正处于破坏状态时以shear-n 或tension-n 标识； 曾进入过屈服状态，但现已经退出以shear-p 或tension-p 标识 根据塑性区标识判断破坏机制是否在起作用，可按下述步骤进行： ④变量监测〔历史跟踪〕：Hist gp zdisp 1；Hist gp szz 2输出2和1的hist关系图用：plot his -2 vs 1 如果应力值都为负数，那么加-2是为了应力坐标轴为正！默认id=1 2 3 4 …依次〔hist id=? Gp zdis x y z〕；hist write 7 v 8 file 6-3hist.txt 6-3hist.txt ⑤切片设置：剖面〔plot set plane ori 0 1.5 0 norm 0 1 0 输出信息：print zone stress 6个方向应力数值 print gp dis 全部节点三个方向的变形大小 SET log on flac3d.loglog状态关闭〔set log off〕为止 rest 6-3.sav set log on set logfile 6-2.log print zone stress print gp dis range id a any id b any set log off 之后重新翻开6-2.log 最大不平衡力不代表没有了塑性流动，只是不平衡力趋于零，并不完全归零。 记录不平衡力随步数的关系。 hist unbal监测不平衡力，并保存历史记录；save t1.sav hist gp zdisp 4,4,8监测网格坐标点〔4，4，8〕在z方向的位移，并保存历史纪录 plot set rot 20 0 30视图的旋转角度为〔20，0，30〕 PLOT block group查看计算模型 Ctrl+G变成灰色图像 采用2种方式观察网格节点速度：：HIST gpvelx y z记录某个节点的速度；PLOT vel力平衡时，由于网格节点力不为零就存在网格节点速度，仍有可能产生较大位移。 初始地应力场：弹性求解法、改变参数的弹塑性求解法以及分阶段弹塑性求解法。 由于为弹性求解，在体系到达平衡时，岩、土体中并未有产生屈服的区域model elas 更改强度参数的弹塑性求解法生成是指求解过程中始终采用塑性模型model mohr 〔先把基点处总的力算出，然后渐变，一般是递减〕，设置初始应力：ini szz 60e3 grad 0 0 20e3 ran z 0 3; ini sxx 30e3 grad 0 0 10e3 ran z 0 3; ini syy 30e3 grad 0 0 10e3 ran z 0 3;水平为竖向的一半。如果都在水下，那么孔隙水压力为ini pp 30e3 grad 0 0 -10e3 ran z 0 3渐变的过程 ：〔存在问题〕？？ 标题栏查看或者执行【Plotitems】/【1 Block group】/【Modify】’s 建立接触面单元〔桩 9.51。9.7 162页。 移来移去法：两个模型分开，一个建立接触面单元，然后移动使两个模型合并。 ini z add -1 range group 。 接触面单元interface () face range cylinder end1 (0,0,-4.9) end2 (0,0,-5.1) radius .31 ，—柱体两端面圆心及半径。 ② 导来导去法：del ran group 2 ，expgrid interface 1 face range x 1 y 1 2 z 1 2 interface 1 face ran