pfcFISH.ppt

目 录 FISH使用技巧若干示例 朱永生 依泰斯卡(武汉)咨询有限公司 ITASCA数值技术研讨会 中国.西安 2009.5.24 FLAC3D to Tecplot UDEC to Tecplot Dxf to UDEC/PFC2D Stl to PFC3D CAD to 3DEC 开始 搜索节点、单元总数，并定义、初始化后期所需要的数据交换数组 读入FLAC3D结果文件，完成相关辅助设置，如开挖单元删除和开挖指标设定等 FISH：_Null_Switch函数 FISH：E_and_GP_Num函数 FISH：_Arr_Para函数 FISH：Array_Define函数 数据交换数组装配，主要完成功能： 节点id 、位移、速度和孔隙水压力计算 单元主应力、正应力和剪应力计算 单元正应变、剪应变计算 基于Mohr-Coulomb模型的单元安全系数计算 弹性应变能密度计算 单元主应力方向计算 如果开挖指示为真，则进行节点、单元编号调整 FISH：_Gp_Var函数 FISH：_El_Stress函数 FISH：_El_Strain函数 FISH：_El_Fos函数 FISH：_El_Energy函数 FISH：_Sig_Dir函数 FISH：_Gp_Swap函数 将数据交换数组输出至外部后处理程序可执行文件 结束 Fig. FLAC3D TO Tecplot程序流程 数据输出交换数据初始化：标题、变量名称、节点和单元结构、各节点和单元变量数据值(位移、速度、应力、应变等) FISH：_Tec_Arr_Assembly函数 FISH：_Tec_Write函数 FLAC3D to Tecplot Fig. 隧洞开挖示例 F2 F1 Fig. 边坡示例 Fig. UDEC TO Tecplot程序流程 开始 搜索节点、单元总数，数据交换数组维数定义 读入UDEC结果文件，完成相关辅助设置 FISH：_Model_Setup函数 FISH：_Info_Get函数 数据交换数组变量装配，包括： 节点类 节点指针 计算速度、位移 单元类 单元指针 单元应力：主应力、正应力、剪应力 单元应变：正应变和剪应变 单元节点信息 FISH：_Gp_Pointer函数 FISH：_GP_Var函数 FISH：_El_Pointer函数 FISH：_El_Stress函数 FISH：_El_Strain函数 FISH：_El_Gp函数 将数据交换数组输出至外部后处理程序可执行文件 结束 数据输出交换数据初始化：标题、变量名称、节点和单元结构、各节点和单元变量数据值(位移、速度、应力、应变等) FISH：_Tec_UDEC_Data函数 FISH：_Tec_Write函数 UDEC to Tecplot Fig. 边坡示例 开始 读入Model.dxf文件,装配数组 参数初始化，包括层面、断层、节理、层间带等结构面总数目和UDEC/PFC2D标识号(id)，另外需说明输出类型，即UDEC或PFC2D FISH：Setup函数 FISH：_Para函数 FISH：_Dxf_Read函数(_Database) 遍历数据(_Database)，搜索各类型结构面总数，提取polyline类Entities，并初始化结构面表，最终将坐标数据装配至各结构面对象数组表中 FISH：_Dxf_Operation_1函数 FISH：_Dxf_Operation_2函数 输出类型？ UDEC PFC2D 按照UDEC程序命令格式整理数据交换数组_Arr_Geom 按照UDEC程序命令格式整理数据交换数组_Arr_Geom 将数据交换数组输出至外部UDEC/PFC2D可执行文件 结束 FISH：_Geom_INI函数 FISH：_Geom函数 Fig. DXF TO UDEC/PFC2D程序流程 DXF to UDEC/PFC2D Fig. DXF-UDEC Fig. DXF-UDEC c1 c2 c3 c4 c5 c6 注释： C1-C6表示层面编号 Fig. DXF-PFC2D C1 C2 C3 C4 C5 C6 注释： C1-C6表示层面编号 Fig. DXF-PFC2D 开始 初始条件设定，包括阻尼条件和输出文件名称 完成模型相关设定，如颗粒(Ball)密度、颗粒数、墙体(Wall)刚度、模型尺寸规模等 FISH：ModelSetup函数 FISH：Simini函数 读入Model.stl文件，主要处理对象： 模型三角网处理 三角网轮廓线处理 结果输出结束 Fig. STL TO PFC3D程序流程 生成颗粒集、墙体集： 生成几何模型 初始化颗粒集，完成颗粒、墙体属性标定 初始化颗粒、墙体速度 FISH：ReadStlFile函数 FISH：MakeAssemb