›› 2010, Vol. 31 ›› Issue (2): 638-644.

• 数值分析 • 上一篇    下一篇

不连续面在双重介质热-水-力三维耦合分析中的有限元数值实现

赵延林1, 2, 3,王卫军1, 2,曹 平3,王 军3,赵阳升4   

  1. 1.湖南科技大学 煤矿安全开采技术湖南省重点实验室,湖南 湘潭 411201;2.湖南科技大学 能源与安全工程学院,湖南 湘潭 411201; 3.中南大学资源与安全工程学院,长沙 410083;4.太原理工大学 采矿工艺研究所,太原 030024
  • 收稿日期:2008-01-05 出版日期:2010-02-10 发布日期:2010-03-24
  • 作者简介:赵延林,男,1973年生,博士研究生,主要从事岩体多场耦合理论和数值仿真研究和教学。
  • 基金资助:

    国家重点基础研究发展计划(973)项目(No. 2007CB209400);国家自然科学基金资助项目(No. 50774093);煤矿安全开采技术湖南省重点实验室开放基金项目(No. 2008001)

Numerical implementation of discontinuities in dual media 3D model for thermo-hydro-mechanical coupling

ZHAO Yan-lin 1, 2, 3, WANG Wei-jun 1, 2, CAO Ping3, WANG-jun3, ZHAO Yang-sheng4   

  1. 1. Hunan Provincial Key Laboratory of Safe Mining Techniques of Coal Mines, Hunan University of Science and Technology, Xiangtan 411201, China; 2. School of Energy and Safety Engineering, Hunan University of Science and Technology, Xiangtan 411201, China; 3. School of Resources and Safety Engineering, Central South University, Changsha 410083, China; 4. Institute of Mining Technology, Taiyuan University of Technology, Taiyuan 030024, China
  • Received:2008-01-05 Online:2010-02-10 Published:2010-03-24

摘要:

在建立双重介质热-水-力耦合微分控制方程的基础上,提出了裂隙岩体热-水-力耦合的三维力学模型,对不同介质分别建立以节点位移、水压力和温度为求解量的三维有限元格式,开发了双重介质热-水-力耦合分析的的三维有限元计算程序,在有限元数值分析中不连续面应力计算采用等厚度空间8节点节理单元进行离散,而不连续面渗流和热能计算时采用平面4节点等参单元进行离散,这样保证了不同介质之间的水量、热量交换和两类模型接触处节点水头、温度和位移相等。通过高温岩体地热开发算例,揭示了在热-水-力耦合作用下不连续面处于低应力区,其张开度随运行时间的延长呈非线性增加,非稳定渗流阶段不连续面显著地控制着渗流场的整体分布,它的水头远高于拟连续岩体介质的水头,而进入稳定渗流阶段不连续面的控渗作用不明显,由于高温岩体地热开发系统中存在大规模的热量补给,不连续面对岩体温度场分布的影响并不显著。

关键词: 不连续面, 双重介质, 热-水-力耦合, 节理单元, 有限元

Abstract:

Based on the differential control equations of the coupled thermo-hydro-mechanical processes of dual media, a three dimensional coupled model of dual media to research coupled thermo-hydro-mechanical(THM) processes of fractured rock masses, is proposed. In coupled model, the values of nodal displacement, water pressure and temperature in different media were solved. A three dimensional FEM code. for analysis of coupled thermo-hydro-mechanical processes of dual media is developed tentatively, in FEM program, 3D jointed elements with 8 nodes were used to simulate the crack elements during the stress analysis of discontinuities; and 2D isoparametric elements with 4 nodes were used to simulate discontinuities during seepage and heat energy analysis of discontinuities. So the different media can exchange water and heat flux, water head, temperature and displacement of nodes which are in the interface between two models were equal. Through the numerical computation for heat extraction in HDR; it is shown: that in the thermo-hydro-mechanical coupling process, discontinuous is in low stress state; the crack width increases nonlinearly with the heat extraction; in transient fluid flow stage, the water head in discontinuities was higher than equivalently continuous rock medinm; in steady fluid flow phase, the dominant effect of discontinuities on fluid flow would unconspicuous. Because of massive heat replenishment, the compact of discontinuities on temperature distribution was indistinctive.

Key words: discontinuities, dual media, thermo-hydro-mechanical coupling, jointed element, finite elements

中图分类号: 

  • TU 452
[1] 孙锐, 杨峰, 阳军生, 赵乙丁, 郑响凑, 罗静静, 姚捷, . 基于二阶锥规划与高阶单元的 自适应上限有限元研究[J]. 岩土力学, 2020, 41(2): 687-694.
[2] 叶观宝, 郑文强, 张 振, . 大面积填土场地中摩擦型桩负摩阻力分布特性研究[J]. 岩土力学, 2019, 40(S1): 440-448.
[3] 赵密, 欧阳文龙, 黄景琦, 杜修力, 赵旭, . P波作用下跨断层隧道轴线地震响应分析[J]. 岩土力学, 2019, 40(9): 3645-3655.
[4] 张海廷, 杨林青, 郭芳, . 基于SBFEM的层状地基埋置管道动力 响应求解与分析[J]. 岩土力学, 2019, 40(7): 2713-2722.
[5] 曹洪, 胡瑶, 骆冠勇. 滤管两端均不在含水层层面的承压不 完整井近似计算方法研究[J]. 岩土力学, 2019, 40(7): 2774-2780.
[6] 王翔南, 李全明, 于玉贞, 喻葭临, 吕禾, . 基于扩展有限元法对土体滑坡破坏过程的模拟[J]. 岩土力学, 2019, 40(6): 2435-2442.
[7] 周小文, 程 力, 周 密, 王 齐, . 离心机中球形贯入仪贯入黏土特性[J]. 岩土力学, 2019, 40(5): 1713-1720.
[8] 梅慧浩, 冷伍明, 聂如松, 刘文劼, 伍晓伟, . 重载铁路路基面动应力峰值随机分布特征研究[J]. 岩土力学, 2019, 40(4): 1603-1613.
[9] 吴顺川, 马 骏, 程 业, 成子桥, 李建宇, . 平台巴西圆盘研究综述及三维启裂点研究[J]. 岩土力学, 2019, 40(4): 1239-1247.
[10] 邱 敏, 袁 青, 李长俊, 肖超超, . 基于孔穴扩张理论的黏土不排水抗剪强度 计算方法对比研究[J]. 岩土力学, 2019, 40(3): 1059-1066.
[11] 徐栋栋, 邬爱清, 李 聪, 汪 斌, 蒋昱州, 曾 平, 杨永涛, . 破裂全过程模拟的改进非连续变形分析方法[J]. 岩土力学, 2019, 40(3): 1169-1178.
[12] 郑安兴, 罗先启, 陈振华, . 基于扩展有限元法的岩体水力劈裂耦合模型[J]. 岩土力学, 2019, 40(2): 799-808.
[13] 王冬勇, 陈曦, 于玉贞, 吕彦楠, . 基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析[J]. 岩土力学, 2019, 40(12): 4890-4896.
[14] 刘庆彬, 潘懋, 刘洁, 郭艳军, 张小双, 姚健鹏, 李芳玉, . 基于ParaView的Abaqus有限元输出结果的 可视化与虚拟现实[J]. 岩土力学, 2019, 40(12): 4916-4924.
[15] 刘锋涛, 张绍发, 戴北冰, 张澄博, 林凯荣, . 边坡稳定分析刚体有限元上限法的锥规划模型[J]. 岩土力学, 2019, 40(10): 4084-4091.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 张文杰,陈云敏. 垃圾填埋场抽水试验及降水方案设计[J]. , 2010, 31(1): 211 -215 .
[2] 宫伟力,安里千,赵海燕,毛灵涛. 基于图像描述的煤岩裂隙CT图像多尺度特征[J]. , 2010, 31(2): 371 -376 .
[3] 万 智,董 辉,刘宝琛. 基于正交设计下SVM滑坡变形时序回归预测的超参数选择[J]. , 2010, 31(2): 503 -508 .
[4] 孙曦源,栾茂田,唐小微. 饱和软黏土地基中桶形基础水平承载力研究[J]. , 2010, 31(2): 667 -672 .
[5] 王明年,郭 军,罗禄森,喻 渝,杨建民,谭忠盛. 高速铁路大断面黄土隧道深浅埋分界深度研究[J]. , 2010, 31(4): 1157 -1162 .
[6] 胡勇刚,罗 强,张 良,黄 晶,陈亚美. 基于离心模型试验的水泥土搅拌法加固斜坡软弱土地基变形特性分析[J]. , 2010, 31(7): 2207 -2213 .
[7] 谭峰屹,姜志全,李仲秋,颜惠和. 附加质量法在昆明新机场填料压实密度检测中的应用研究[J]. , 2010, 31(7): 2214 -2218 .
[8] 柴 波,殷坤龙,肖拥军. 巴东新城区库岸斜坡软弱带特征[J]. , 2010, 31(8): 2501 -2506 .
[9] 杨召亮,孙冠华,郑 宏. 基于潘氏极大值原理的边坡稳定性的整体分析法[J]. , 2011, 32(2): 559 -563 .
[10] 王光进,杨春和,张 超,马洪岭,孔祥云,侯克鹏. 超高排土场的粒径分级及其边坡稳定性分析研究[J]. , 2011, 32(3): 905 -913 .