基于ANSYS平台复杂地质体FLAC3D模型的自动生成 (1)

第 24卷 第 6 期 岩石力学与工程学报 Vol.24 No.6 2005年 3月 Chinese Journal of Rock Mechanics and Engineering March 2005 收稿日期 2003–11–08 修回日期 2003–12–01 基金项目 国家重点基础研究发展规划(973)项目(2002CB412702) 中国科学院知识创新工程资助项目(KJCX2–SW–L1–1) 作者简介 廖秋林(1977–) 男 现为博士研究生 主要从事工程地质与地质灾害方面的研究工作 E-mail qiulinliao@yahoo.com.cn 基于 ANSYS平台复杂地质体 FLAC3D模型的自动生成 廖秋林 1 曾钱帮 1 刘 彤 1 路世豹 2 侯哲生 1 (1. 中国科学院 地质与地球物理研究所工程地质力学重点实验室 北京 100029 2. 中国矿业大学 建筑工程学院 北京 100083) 摘要 由于 FLAC3D软件建模难度大 提出了一种快速建模 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 即以 ANSYS 有限元程序完成的复杂地质体建 模 网格划分为基础 采用 Visual Basic 语言编写了 FLAC3D-ANSYS接口程序 实现了 FLAC3D软件建模的直观 快速和自动化 通过新疆下坂地水库坝址区和北京地铁四号线车站三维建模实例检验了该方法的有效性和可行 性 该方法是运用不同程序优点解决复杂工程地质问题的典型范例 关键词 数值 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 接口程序 复杂地质体 自动生成 中图分类号 O 241 文献标识码 A 文章编号 1000–6915(2005)06–1010–04 AUTOMATIC MODEL GENERATION OF COMPLEX GEOLOGIC BODY WITH FLAC3D BASED ON ANSYS PLATFORM LIAO Qiu-lin1 ZENG Qian-bang1 LIU Tong1 LU Shi-bao2 HOU Zhe-sheng1 (1. Key Laboratory of Engineering Geomechanics Institute of Geology and Geophysics The Chinese Academy of Sciences Beijing 100029 China 2. School of Architecture China University of Mining and Technology Beijing 100083 China) Abstract For some complex geologic bodies numerical models with FLAC3D are difficult to construct. While finite element programs such as ANSYS can easily fulfill since it has perfect pre-processors such as entity modeling Boolean operation of geometric body and free mesh. In this paper a procedure is proposed to build numerical model of FLAC3D by using ANSYS and the interface program of FLAC3D-ANSYS is compiled with visual basic language since there are some differences between the element data conducted by the two programs. So the automatic model generation of complex geologic body can be conducted as follows the construction of numerical model by ANSYS the transferring of element data and calling FLAC3D. Then three-dimensional numerical models of a dam and a subway are developed and the case study shows that the intuition fastness and automation of modeling can be realized. The procedure gives an example to build complicated engineering geological model with the advantages of different programs. Key words numerical analysis interface program complex geologic body automatic generation 1 引 言 FLAC3D 软件是率先将连续体的快速拉格朗日 分析应用于岩土工程问题的计算软件 该软件在解 决岩土工程问题上具有许多优越性[1] 已逐渐成为 工程技术人员理想的三维数值模拟工具 然而 FLAC3D 软件在模型建立以及单元网格划分等前处 理问题上却存在以下不足 造成了其建模的不便性 (1) 模型的建立只能靠数据文件来实现 不是 万方数据 第 24卷 第 6期 廖秋林等. 基于ANSYS平台复杂地质体FLAC3D模型的自动生成 • 1011 • 很直观 不能像ANSYS或 ALGOR等有限元软件 可以直接进行图形的处理 (2) 对于比较复杂的工程模型 在建模时需要 各控制点详细的数据 容易出错 检查起来也不是 很容易 (3) 建模工作量大 花费时间长 直接造成了 三维模拟计算的周期长 难度大 为解决 FLAC3D 软件建模的不足 文[2]采用 FORTRAN语言专门编写了 FLAC3D的前处理程序 对于地表形态复杂 岩层和地质结构较单一的地质 体实现了快速 便捷的建模 文[3 4]研究了 FLAC3D 二维平面应变模型的快速建立 但是对于采动影响 下的矿山工程和地质结构错综复杂的边坡等地质 体 其 FLAC3D模型的建立及网格划分仍然非常不 便 然而 许多有限元程序对于复杂工程地质体数 值模型的建立有明显优势 由于在所研究对象的前 处理上 有限元法和有限差分法基本一致 即将研 究对象划分为许多具有一定形状且有一定编制 规则 编码规则下载淘宝规则下载天猫规则下载麻将竞赛规则pdf麻将竞赛规则pdf 的节点组成的单元体(面) 因此 作者试图利用已 有的建模与网格划分功能强大的 ANSYS 软件对 复杂工程地质体建立相应数值模型(包括网格划 分) 再通过数据转化实现 FLAC3D 模型的自动生 成 2 FLAC3D模型的自动生成 FLAC3D和ANSYS所采用的单元体形状大都相 同 但其每一单元节点编制的规则和节点坐标 即 单元数据 却有一定的差别 因此 根据这 2种软 件单元形状及其单元数据的关系编写 FLAC3D- ANSYS 接口程序是本文建议建模方法的关键 FLAC3D 模型的自动生成主要由以下步骤组成 ANSYS 模型的建立 ANSYS 和 FLAC3D的数据转 换以及 FLAC3D调用生成的模型数据文件 2.1 ANSYS模型的建立 ANSYS作为广泛应用的有限元软件 充分综合 了 CAD CAE CAM 等图像处理工具[5] 是建立 复杂计算模型有效而又方便快捷的平台 ANSYS 可以自上而下直接建立实体模型 还可通过自下而 上依次生成点 线 面和体 从而建立实体模型 强大的布尔运算工具可以实现实体之间加 减 分 类 搭接 粘接和分割等复杂运算 大大提高了建 立复杂地质体三维模型的效率 对于实体模型的网 格划分 ANSYS提供了功能强大的控制工具 如单 元大小和形状的控制 网格的划分类型(自由和映射) 以及网格的清除和细化 此外 还可对实体模型图 直接划分网格 最后 输出各单元节点坐标及单元 信息 NODE. DAT 和 ELE. DAT文件 2.2 FLAC3D与 ANSYS单元数据关系 要将ANSYS所生成的单元数据文件为FLAC3D 所利用 有必要掌握 FLAC3D与 ANSYS单元数据之 间的关系 在模拟对象的单元处理上 FLAC3D 与 ANSYS都提供了丰富的单元形状 根据地质体的 特征 计算精度要求以及单元形状的空间展布特 点 仅考虑以下 5种单元体 六面体 五面楔形体 五面锥形体 四面体和圆柱体 这 5种单元体基本 能满足各种地质体数值模型的建立 这 2种软件所 采用单元节点编制对应关系如表 1所示 此外 由 ANSYS 存在单元退化和二次单元等问题 而 FLAC3D则只能通过对 ANSYS单元退化节点的判断 用低节点的单元替换退化的高节点单元 例如 六 面体可退化为五面楔形体 五面锥形体和四面体 如表 1所示 表 1 ANSYS 与 FLAC3D单元数据关系对照 Table 1 Computation of element data between ANSYS and FLAC3D 单元数据 单元类型 ANSYS FLAC3D 备注 六面体 Brick单元 五面楔形体 其 ANSYS 单元 为Brick退化单元 五面锥形体 其 ANSYS 单元 为Brick退化单元 四面体(I) 其 ANSYS 单元 为Brick退化单元 四面体(II) 根据以上对 FLAC3D与 ANSYS 单元数据关系 的分析 作者利用 Visual Basic语言编写了 FLAC3D- ANSYS接口程序包 首先 该程序将单元节点坐标 万方数据 1012 岩石力学与工程学报 2005 转化成 FLAC3D单元节点坐标 其次 根据 ANSYS 提供的单元信息 ELE. DAT文件的格式特点 该程序 自动判断其每一单元的形状(也考虑了退化单元的 转换) 并生成相应的 FLAC3D单元 该程序除实现 了 2 种软件的单元数据的转换之外 还将 ANSYS 定义的不同实体遗传到 FLAC3D中 并形成相应的 Group 方便了计算参数的赋值 最后 该程序直接 产生 FLAC3D计算所需的数据文件 其主要内容包 括 (1) 产生节点的命令 主要格式为 GEN POINT ID 1 0.000 0.000 414.075 (2) 产生单元的命令 主要格式为 GEN ZONE BRICK(WEDGE)SIZE 1 1 1 P0 and POINT 5532 P1 POINT 5393 P2 POINT 719 P3 and POINT 5274 P4 POINT 580 P5 POINT 461 P6 and POINT 5511 P7 POINT 698 (3) 产生实体的命令 主要格式为 GROUP 5 RANGE ID 3591 2.3 模型的自动生成 最后 通过 FLAC3D命令 call 调入由接口程 序输出的数据文件 并加入边界条件 初始条件以 及岩土体的力学参数 即可生成数值模型 3 应用实例 3.1 新疆下坂地水库坝址区三维数值模型 新疆下坂地水库坝址区 峡谷长约 1 200 m 峡 谷区岸坡陡峭 比高近 600 m 峡谷及其两侧斜坡 累计宽度达 2 500 m 其岩层之间的接触关系比较复 杂 对于这样一个地表形态和地质结构复杂的坝 址 采用本文提出的建模方法可以很理想地建立其 三维模型 具体操作如下 (1) 准备几何图形数据 首先 将复杂地表划 分为一系列三角形面 以三角形的顶点作为主要地 表数据点 然后 根据岩层的产状 定义岩层接触 面控制点和模型边界点 最后 将这些数据做成 ANSYS的脚本文件 (2) 生成实体模型 先根据点数据生成图 1(a) 的表面形态 再沿三角形面向下生成一系列柱体 利用布尔运算相加成 4个体(每一岩层为一个实体) 之后 以模型底面为基础向上生成一个长方体 最 后 将前面生成的 4 个体与长方体相减 即得到 图 1(a) 图 1 新疆下坂地水库坝址区三维数值模型 Fig.1 Numerical model in the Xiabandi reservoir area Xinjiang Autonomous Region (3) 划分网格 该模型选择了四面体单元划分 网格并根据岩性定义了 3种材料 在网格划分时 根据计算需要进行网格密度和重点研究区网格细化 的控制 并通过 ANSYS 网格检验达到网格最优 化 最终 该模型共划分 2 357个节点 9 987个单 元 (4) 生成 FLAC3D三维模型 利用 ANSYS生成 的NODE. DAT和ELE. DAT文件和作者开发的接口 程序生成 FLAC3D 可调用的脚本文件 最后生成 FLAC3D三维数值模型 如图 1(b)所示 整个建模从数据采集到模型开始计算 共用时 不到 2 d 且实现了单元的最优化 3.2 北京地铁四号线车站数值模型 北京地铁四号线车站为一大型地下硐室群结 构 其两边隧道与车站的衔接处是三维数值模型建 立的难点 该模型若用 FLAC3D所提供的基本单元 (如 radcylinder或 cylinder)来处理这一复杂硐室群 而且要保证圆弧之间的衔接和实体之间网格的数目 匹配(FLAC3D 要求数目一致或成倍数关系 以免计 算得出结果在此处不连续) 难度很大 即使利用 FLAC3D中的 FISH语言来处理这一复杂硐室群 其 结果也并不理想 且耗时达一周之久 而根据本文所提出的建模方法 在获得该模型 的任一横截面的几何数据坐标后 直接利用 ANSYS 绘制出模型的横截面图 再对该截面划分网格形成 (a) 河谷 基岩 斜坡 (b) 万方数据 第 24卷 第 6期 廖秋林等. 基于ANSYS平台复杂地质体 FLAC3D模型的自动生成 • 1013 • 二维模型 之后利用 ANSYS 中的 Vdrag命令将该 截面沿指定方向拉成体 即完成了 ANSYS的建模 并生成了相应的数据文件 最后 利用接口程序生 成的 FLAC3D数据文件 建立三维数值模型 如图 2 所示 该模型建立的每一环节都简单 易操作 整 个建模时间不到半天 而且 该模型既美观 又符 合 FLAC3D的要求 其计算结果也令人满意 图 2 北京地铁四号线车站三维数值模型 Fig.2 Numerical model of a station of Beijing subway No.4 4 结 论 本文基于分析功能强大的岩土工程计算软件 FLAC3D 提出了基于 ANSYS 平台复杂地质体 FLAC3D数值模型的自动生成方法 该方法不仅克 服了建立复杂计算模型的困难 而且实现了建模的 自动化 其流程的每一个操作均简单 易行 所要 求的基本数据仅为几何图形数据点 大大减少了建 模所需的时间 精力 提高了数值模拟的频度 有 利于充分发挥 FLAC3D解决工程地质问题的强大功 能 也促进了 FLAC3D在岩土工程领域的普及和应 用 该方法的主要思路着眼于发挥和结合 FLAC3D 与 ANSYS软件的优点 并开发了 FLAC3D-ANSYS 接口程序 解决了复杂工程地质体的建模难题 在 科技手段发达且种类繁多的当今社会 对于复杂工 程地质问题的数值模拟 本文提出的方法值得借 鉴 参考文献(References) [1] Itasca Consulting Group Inc.. Fast Language Analysis of continua in three dimensions[R]. [s.l.] Itasca Consulting Group Inc. 1997. [2] 胡 斌 张倬元 黄润秋等. FLAC3D前处理程序的开发及仿真效 果检验[J]. 岩石力学与工程学报 2002 21(9) 1 387–1 391.(Hu Bin Zhang Zhuoyuan Huang Runqiu et al. Development of pre-processing package for FLAC3D and verification of its simulating effects[J]. Chinese Journal of Rock Mechanics and Engineering 2002 21(9) 1 387–1 391.(in Chinese)) [3] 油新华. 土石混合体随机结构模型及其应用研究[博士学位论 文][D]. 北方交通大学 2001.(You Xinhua. Stochastic structural model of the earth-rock aggregate and its application[Ph. D. Thesis][D]. Beijing Northern Jiaotong University 2001.(in Chinese)) [4] 吴剑波. 三峡库区白衣庵滑坡稳定性研究[硕士学位 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 ][D]. 北 京 中国矿业大学 2002.(Wu Jianbo. Study of the stability on Baiyi¢an ancient landslide in the Three Gorges District [M. S. Thesis][D]. Beijing China University of Mining and Technology 2002.(in Chinese)) [5] 边 萌. ANSYS 5.7 有限元实例分析教程[M]. 北京 机械工业出 版社 2002.(Bian Meng. The case-analysing tutorial by the finite element program ANSYS 5.7[M]. Beijing China Mechanical Industry Press 2002.(in Chinese)) 万方数据 基于ANSYS平台复杂地质体FLAC3D模型的自动生成 作者： 廖秋林， 曾钱帮， 刘彤， 路世豹， 侯哲生， LIAO Qiu-lin， ZENG Qian-bang， LIU Tong， LU Shi- bao， HOU Zhe-sheng 作者单位： 廖秋林,曾钱帮,刘彤,侯哲生,LIAO Qiu-lin,ZENG Qian-bang,LIU Tong,HOU Zhe-sheng(中国科学院,地质与 地球物理研究所工程地质力学重点实验室,北京,100029)， 路世豹,LU Shi-bao(中国矿业大学,建筑工程学院 ,北京,100083) 刊名： 岩石力学与工程学报 英文刊名： CHINESE JOURNAL OF ROCK MECHANICS AND ENGINEERING 年，卷(期)： 2005,24(6) 被引用次数： 59次 参考文献(5条) 1.边萌 ANSYS 5.7 有限元实例分析教程 2002 2.吴剑波 三峡库区白衣庵滑坡稳定性研究[学位论文] 2002 3.油新华 土石混合体的随机结构模型及其应用研究[学位论文] 2001 4.胡斌;张倬元;黄润秋 FLAC3D前处理程序的开发及仿真效果检验[期刊论文]-岩石力学与工程学报 2002(09) 5.Itasca Consulting Group Inc Fast Language Analysis of continua in three dimensions 1997 引证文献(61条) 1.杨宇宏.许胜 基于Surfer和ANSYS复杂地质体的FLAC3D实体建模[期刊论文]-山西建筑 2011(5) 2.陈星.郭永成.赵二平 综合考虑安全性和经济性的料场边坡开挖 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 优选[期刊论文]-三峡大学学报（自然科学版） 2010(1) 3.韦玮.严明 某水电站右岸缆机平台边坡稳定性数值分析[期刊论文]-山西建筑 2010(1) 4.王树仁.张海清 MIDAS/GTS-FLAC~(3D)耦合建模新方法及其应用[期刊论文]-土木建筑与环境工程 2010(1) 5.梁元贵.刘关锋 易门铜厂铜矿露天开采最终边坡角分析[期刊论文]-采矿技术 2010(1) 6.荆永滨.王李管.贾明涛.陈建宏 基于TetGen的复杂FLAC3D模型可视化建模方法[期刊论文]-岩土力学 2010(8) 7.成建.康永红.吴东 采矿方法对地表影响的数值模拟预测研究[期刊论文]-云南冶金 2010(3) 8.余超.陈星.邓华锋.李俊美 大山口二级水电站引水隧洞围岩稳定性研究[期刊论文]-水电能源科学 2010(9) 9.徐雪飞.许劲松.赵鹏.迟颖 龙王塘水库大坝预应力锚固工程数值模拟研究[期刊论文]-东北水利水电 2010(11) 10.刘心庭.唐辉明 FLAC3D复杂网格模型的构建及其工程应用[期刊论文]-金属矿山 2010(11) 11.吕韬.杨建国.吴孟.何飞 建设用地地灾危险性评估的数值模拟研究[期刊论文]-地下空间与工程学报 2010(6) 12.钟正强.彭振斌.彭文祥 基于Hock-Brown准则的层状岩体隧道开挖响应分析[期刊论文]-中南大学学报（自然科学版） 2009(6) 13.张雨霆.肖明.左双英 基于单元重构的岩土工程复杂地质断层建模方法[期刊论文]-岩石力学与工程学报 2009(9) 14.高盛翔.叶容章.唐伟励.徐清.贺金强 应用MATLAB接口程序建立复杂地质体FLAC3D模型[期刊论文]-煤田地质与勘探 2009(5) 15.刘晓明.罗周全.张保.李畅 上向分层回采采场稳定性数值模拟研究[期刊论文]-矿冶工程 2009(4) 16.袁海平.黄晖.苏成哲.周雄超 隐伏多空区对露天边坡稳定性的影响研究[期刊论文]-江西有色金属 2009(4) 17.蒋英礼.赵伯明.胡晓勇 软土地铁车站中柱在强震作用下的响应研究[期刊论文]-防灾减灾工程学报 2009(4) 18.赵伯明.蒋英礼.陈靖 软土地铁车站结构在三维强地震动作用下的响应分析[期刊论文]-中国铁道科学 2009(3) 19.王树仁.张海清.慎乃齐.曹海莹 下伏采空区桥隧工程变形及受力响应特征分析[期刊论文]-岩石力学与工程学报 2009(6) 20.艾显恒.谢建兵.朱和玲 地下矿山采场顶板厚度数值模拟研究[期刊论文]-现代矿业 2009(4) 21.柴红保.龙雪鸣 复杂地质体FLAC3D模型快速生成[期刊论文]-矿业工程研究 2009(2) 22.李鹏.林志.张亚兴 高地应力区断层带隧道施工力学研究[期刊论文]-重庆交通大学学报（自然科学版） 2009(1) 23.艾显恒.谢建兵.朱和玲 地下矿山采场顶板厚度数值模拟研究[期刊论文]-采矿技术 2009(2) 24.柴红保.万文 露天矿边坡与地下采空区相互作用的数值分析[期刊论文]-矿业工程研究 2009(3) 25.樊敬亮.金雪莲.梁玉红.曾钱帮 高山峡谷地区坝陵河大桥西岸边坡稳定性及工程效应[期刊论文]-解放军理工大学学报（自然科学版 ） 2009(5) 26.王树仁.张海清.慎乃齐 穿越采空区桥隧工程危害效应分析及对策[期刊论文]-解放军理工大学学报（自然科学版） 2009(5) 27.张玉灯 基于ADINA平台的FLAC3D建模[期刊论文]-西北水电 2008(4) 28.罗周全.刘晓明.吴亚斌.刘望平.张保 基于Surpac和Phase2耦合的采空区稳定性模拟分析[期刊论文]-辽宁工程技术大学学报（自然科 学版） 2008(4) 29.张玉灯 FLAC3D在岩质边坡稳定性分析中的应用[期刊论文]-路基工程 2008(6) 30.罗靖.邹荣华 桩基弦振效应的数值模拟与性能分析[期刊论文]-荆门职业技术学院学报 2008(6) 31.喻波.黄政宇 基于ANSYS的ABAQUS强度折减边坡稳定性分析[期刊论文]-公路工程 2008(2) 32.林杭.曹平.李江腾.江学良.何忠明 基于SURPAC的FLAC3D三维模型自动构建[期刊论文]-中国矿业大学学报 2008(3) 33.杨宇宏.许胜 基于Surfer和ANSYS复杂地质体的FLAC3D实体建模[期刊论文]-山西建筑 2011(5) 34.罗周全.吴亚斌.刘晓明.刘望平.杨彪 基于SURPAC的复杂地质体FLAC3D模型生成技术[期刊论文]-岩土力学 2008(5) 35.李荣伟.侯恩科.刘天林 基于三维数值模拟的露天煤矿复合边坡稳定性分析[期刊论文]-煤田地质与勘探 2008(2) 36.李根.赵娜 以ANSYS为平台的复杂模型到FLAC3D导入技术[期刊论文]-辽宁工程技术大学学报（自然科学版） 2008(z1) 37.邓红卫.朱和玲.周科平.肖雄 基于FLAC3D数值模拟的前后处理优化研究[期刊论文]-矿业研究与开发 2008(2) 38.李扬.梅林芳.周传波 露天转地下崩落法开采对高陡边坡影响的数值模拟[期刊论文]-矿冶工程 2008(3) 39.徐文杰.胡瑞林 基于逆向工程的三维复杂地质体精细建模及ADINA前处理在FLAC3D建模中的应用[期刊论文]-工程地质学报 2008(5) 40.崔芳鹏.胡瑞林.刘照连.张明.余文龙 基于Surfer平台的FLAC3D复杂三维地质建模研究[期刊论文]-工程地质学报 2008(5) 41.汪吉林.丁陈建.吴圣林 基于FLAC3D的复杂地貌三维地质建模[期刊论文]-地质力学学报 2008(2) 42.朱和玲.周科平.肖雄.胡建华 基于开采环境再造顶板最小安全厚度研究[期刊论文]-矿冶工程 2008(1) 43.罗周全.鹿浩.刘晓明.李畅.张保 基于CMS实测的采空区群稳定性数值模拟方法[期刊论文]-化工矿物与加工 2008(1) 44.张华.陆阳.程强 岩质边坡锚杆(索)框架梁加固的数值模拟[期刊论文]-公路交通科技 2008(1) 45.郑文棠.徐卫亚.童富果.石安池 复杂边坡三维地质可视化和数值模型构建[期刊论文]-岩石力学与工程学报 2007(8) 46.毕林.王李管.荆永滨 多介质复杂结构条件下岩体工程开挖力学响应三维数值模拟技术[期刊论文]-湖南科技大学学报（自然科学版） 2007(3) 47.胡义生.张玉峰 FLAC3D在公路边坡工程中的应用[期刊论文]-西部探矿工程 2007(10) 48.刘建华 连续介质力学行为分析的静态拉格朗日方法[期刊论文]-岩石力学与工程学报 2007(7) 49.Guoqing YIN.Wei JIN.Xianli YANG FEM simulation of formation of metamorphic core complex with ANSYS software[期刊论文 ]-世界地质（英文版） 2007(1) 50.董鹏.李长洪 基于ANSYS软件的FLAC3D实体建模[期刊论文]-山西建筑 2007(3) 51.卢敦华.何忠明.林杭 节理岩坡中注浆效应的FLAC3D数值分析[期刊论文]-煤田地质与勘探 2007(2) 52.邹力.彭雄志 浅谈FLAC-3D的应用原理、优缺点及改进措施[期刊论文]-四川建筑 2007(1) 53.单宏兰 金属矿山巷道围岩稳定性分析与监测点位置优化研究[学位论文]硕士 2007 54.徐鼎平.汪斌.江龙剑.曹作忠 太和铁矿西端帮边坡工程地质模型及其三维剖分方法[期刊论文]-矿业快报 2006(11) 55.吕淑然.刘红岩 矿体开采的FLAC3D数值模拟分析[期刊论文]-矿冶 2006(4) 56.朱昌星.阮怀宁.朱珍德.罗润林 锦屏深埋长大引水隧洞围岩蠕变特性仿真分析[期刊论文]-地下空间与工程学报 2006(6) 57.李小民 竹城公路边坡全长注浆锚杆支护参数优化[期刊论文]-湖南交通科技 2006(2) 58.李地元 高速公路连拱隧道围岩应力场和渗流场的耦合作用分析[学位论文]硕士 2006 59.林杭 竹城公路层状岩质边坡稳定性分析与防治措施研究[学位论文]硕士 2006 60.杨洪杰 土压盾构掘进对土体影响的试验研究与数值模拟[学位论文]博士后 2006 61.吕韬.杨建国.吴孟.何飞 建设用地地灾危险性评估的数值模拟研究[期刊论文]-地下空间与工程学报 2010(6) 本文链接：http://d.g.wanfangdata.com.cn/Periodical_yslxygcxb200506019.aspx