FLAC3D学习笔记.doc

gen zone copy 0 1.8 0 ran x 0.9 2.7 ;将单元复制 到（0,1.8,0）以原点为基点gen zone reflect normal 0 1 0 ori 0 1.8 0 ran y 1.8 2.7 ;normal为对称面法线（或倾向倾角 dip,dd）ori 0 1.8 0面过这点ini z add -21.0 range group zhuang ;group zhuang单元下移 21mint 1 del ;删除接触面pl block state ;查看塑 性屈服单元ini state 0 ;认为不存在 塑性状态，计算时重新从弹性状态开始 ini xdis ydis zdis 0 ;将位移重 新设为 0，不将平衡阶段发生的位移记录在内pl block group ran pile ;调整 scale 和 grid point ID's查看节点号或历史记录查看hist unbal ;记 录不平衡力，默认 ID为 1hist id 2 gp zd 0 0 0 ;历史记录并 在图上展示计算过程pl ad hist 2pl showpl gpfix red sk ;查看 模型边界条件pl set plane ori norm ;切片（过 ori 点和 norm 法向）pl con sz plane ;显 示 Z 方向应力的切片pl fos con ssi outline on vel red ;查看安 全系数，剪应力增量云图hist write 1 vs 2 file 1.xls ;将历史记录 ID 为 1,2 的对应关系输入到 Excel 中expgrid 5-1.flac3dimpgrid name ;网格输出hist id 2 gp zdis id 3662set large ;计算过程 网格发生变形set logfile zdisp_1.dat ;命名日 志文件set log on print gp disp ran id 9844 ;输出 ID=9844点的所有位移del zon ran gr d1 any gr d2 any gr d3 anyset log offset movie avi step 1 file 6-5.avi ;求 解过程录制影片movie startprint sel node link range id 2 ;获得 id 为2 的节点的链接的 idsel delete link range id 29 ;删除 ID 为 29的链接sel link id=100 2 targe zone ;在 ID为 2 的节点上新建 ID为 100的链接print sel link attach ;查看链接 情况interface 1 wrap Material Bin range plane ori 0 0 0 normal 1 -1 0 aboveinterface 2 wrap Material Bin range plane ori 0 0 0 normal 1 -1 0 below;生成庄周接触面interface 1 face range cylinder end1 (0,1.8,0.1) end2 (0,1.8,-20.1) radius .31 &cylinder end1 (0,1.8,0.1) end2 (0,1.8,-20.1) radius .29 not;生成庄端接触面interface 2 face range cylinder end1 (0,1.8,-19.9) end2 (0,1.8,-20.1) radius .31;桩与承台刚接----之所以 radius<0.3 是防止与土刚接（桩径 0.3）0.29 not 表示大于0.29 的第一个边界gen merge 1e-3 ran cylinder end1 (0,0,-0.01) end2 (0,0,0.01) radius 0.29modle mohr property K G friction cohesion tension dilation(摩尔 -库伦模型有6 个参数，不赋值的默认 0）modle elasticproperty K G (弹性模型只有两个 参数）K=E/3（1-2v）G=E/2（1+v ）fish 保留字pnt = gp_near(0,0,10) ;靠近坐标（0， 0，10）的节点的地址赋予 pntpnt = gp_find(id) ;ID 为 id 节点的地 址赋予 pnt。