温度场及温度应力模拟的高阶非连续变形分析方法

doi:10.13722/j.cnki.jrme.2019.0626

摘要
图/表
参考文献
相关文章 (15)

全文: PDF (330 KB) HTML (1 KB)
输出: BibTeX | EndNote (RIS)

摘要对于某些岩体工程问题，如核废料地下存储，温度应力是必须考虑的一个重要影响因素，因为它可能会导致岩石发生破裂而带来危害性后果。传统非连续变形分析方法(DDA)，作为岩体稳定分析的主流数值方法之一，尚无法模拟岩体内部的温度场和温度应力分布。针对这一情况，通过在DDA块体外覆盖一等边三角形来构造块体上的高阶位移函数和高阶温度场函数，辅之以强制位移连续和温度场连续的接缝单元，建立可用于模拟温度场及温度应力分布的高阶非连续变形分析方法。该方法建立起DDA与基于单位分解的有限元法(FEM)之间的联系，可直接借鉴FEM成熟的温度场及温度应力数值分析理论。最后，将所建立的非连续变形分析方法应用到解析算例的温度场和温度应力数值求解中。结果表明，模拟结果与解析解非常接近，具有较高的精度。证实了所提方法的正确性及有效性。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	徐栋栋
	邬爱清
	卢波
	汪斌
	蒋昱州
	李聪
	刘小红

关键词 ：数值分析, 温度场, 温度应力, 非连续变形分析, 有限元法, 接缝单元

Abstract：For some rock engineering problems，such as underground storage of nuclear waste，thermal stress is an important factor to be considered，because it may lead to rock fracture and bring harmful consequences. Traditional discontinuous deformation analysis method (DDA)，as one of the main numerical methods for stability analysis of rock mass，cannot simulate the distribution of temperature field and thermal stress in rock mass. In order to solve this problem，a high-order discontinuous deformation analysis method for simulating the distribution of temperature field and thermal stress is established by constructing higher-order displacement function and higher-order temperature field function on DDA block covered by an equilateral triangle outside it，supplemented by joint elements used to force the continuities of displacement and temperature between blocks. This method has established the relationship between DDA and partition of unity based finite element method (FEM). The mature theory of FEM in the simulation of temperature field and thermal stress can be directly used for reference. Finally，the established discontinuous deformation analysis method is applied to the numerical solution of temperature field and thermal stress in the analytical examples. The results show that the simulated results are very close to the analytical solutions and have very high accuracy. The correctness and effectiveness of the proposed method have been verified.

Key words： numerical analysis temperature field thermal stress discontinuous deformation analysis finite element method joint element

引用本文:

徐栋栋，邬爱清，卢波，汪斌，蒋昱州，李聪，刘小红. 温度场及温度应力模拟的高阶非连续变形分析方法[J]. 岩石力学与工程学报, 2019, 38(S2): 3595-3602.
XU Dongdong，WU Aiqing，LU Bo，WNAG Bin，JIANG Yuzhou，LI Cong，LIU Xiaohong. High-order discontinuous deformation analysis method for simulation of temperature field and thermal stress. , 2019, 38(S2): 3595-3602.

链接本文:

http://rockmech.whrsm.ac.cn/CN/10.13722/j.cnki.jrme.2019.0626 或 http://rockmech.whrsm.ac.cn/CN/Y2019/V38/IS2/3595

[2]	黄润秋，王贤能，唐胜传. 热应力的交变作用对边坡危岩体形成的影响[J]. 自然科学进展，1999，9(8)：716–722.(HUANG Runqiu，WANG Xianneng，TANG Shengchuan. Thermal stress effect on slope rock mass[J]. Nature Science Process，1999，9(8)：716–722.(in Chinese))
[9]	林绍忠，明峥嵘，祁勇峰. 用数值流形法分析温度场及温度应力[J]. 长江科学院院报，2007，24(5)：72–75.(LIN Shaozhong，MING Zhengrong，QI Yongfeng. Thermal field and thermal stress analysis based on numerical manifold method[J]. Journal of Yangtze River Scientific Research Institute，2007，24(5)：72–75.(in Chinese))
[12]	谭育新，张慧华，胡国栋. 二维稳态热传导问题的正六边形流形元研究[J]. 应用数学和力学，2017，38(5)：594–604.(TAN Yuxin，ZHANG Huihua，HU Guodong. 2D steady heat conduction analysis with the regular hexagon numerical manifold method[J]. Applied Mathematics and Mechanics，2017，38(5)：594–604.(in Chinese))
[19]	JIAO Y Y，HUANG G H，ZHAO Z Y，et al. An improved three-dimensional spherical DDA model for simulating rock failure[J]. Science China Technological Sciences，2015，58(9)：1 533–1 541.
[1]	唐世斌，唐春安，朱万成，等. 热应力作用下的岩石破裂过程分析[J]. 岩石力学与工程学报，2006，25(10)：2 071–2 078.(TANG Shibin，TANG Chunan，ZHU Wancheng，et al. Numerical investigation on rock failure process induced by thermal stress[J]. Chinese Journal of Rock Mechanics and Engineering，2006，25(10)：2 071–2 078.(in Chinese))
[11]	ZHANG H H，MA G W，REN F. Implementation of the numerical manifold method for thermo-mechanical fracture of planar solids[J]. Engineering Analysis with Boundary Elements，2014，44 ：45–54.
[21]	NING Y，YANG J，AN X，et al. Modelling rock fracturing and blast-induced rock mass failure via advanced discrimination within the discontinuous deformation analysis framework[J]. Computers and Geotechnics，2011，38(1)：40–49.
[22]	FU X，SHENG Q，ZHANG Y，et al. Boundary setting method for the seismic dynamic response analysis of engineering rock mass structures using the discontinuous deformation analysis method[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2015，39(15)：1 693–1 712.
[29]	杨帆. 数值流形法的接缝单元及其应用研究[硕士学位论文][D]. 武汉：长江科学院，2010.(YANG Fan. Research on the joint elements of numerical manifold method and its application[M. S. Thesis][D]. Wuhan：Changjiang River Scientific Research Institute，2010.(in Chinese))
[7]	王峰. 改进的无网格法求解温度场和温度应力及其在水工结构分析中的应用[博士学位论文][D]. 大连：大连理工大学，2015. (WANG Feng. Improved meshless method for the solution of temperature field and thermal stress and its application to hydraulic structure analysis[Ph. D. Thesis][D]. Dalian：Dalian University of Technology，2015.(in Chinese))
[17]	徐栋栋，邬爱清. 基于显式时间积分的非连续变形分析方法[J]. 岩石力学与工程学报，2016，35(增2)：3 526–3 534.(XU Dongdong，WU Aiqing. Discontinuous deformation analysis method based on the explicit time integration scheme[J]. Chinese Journal of Rock Mechanics and Engineering，2016，35(Supp.2)：3 526–3 534.(in Chinese))
[27]	JIANG Q H，YEUNG M R，SUN N. Modelling groundwater pressure and thermal loading in three-dimensional discontinuous deformation analysis[M]. [S. l.]：[s. n.]，2004：471–476.
[5]	李连崇，唐春安，唐世斌，等. 含损伤演化的TM耦合数值模型及其应用研究[J]. 力学学报，2006，38(4)：514–521.(LI Lianchong，TANG Chunan，TANG Shibin，et al. The damage coupled thermal- mechanical(TM)model for rock failure process and applications[J]. Chinese Journal of Theoretical and Applied Mechanics，2006，38(4)：514–521.(in Chinese))
[15]	SHI G H. Discontinuous deformation analysis： a new numerical model for the statics and dynamics of deformable block structures[J]. Engineering computations，1992，9(2)：157–168.
[25]	王书法，朱维申，郑少河. 节理岩体热应力问题的非连续变形分析方法[J]. 岩石力学与工程学报，2000，19(5)：643–646.(WANG Shufa，ZHU Weishen，ZHENG Shaohe. DDA method of thermal stress problem for jointed rock mass[J]. Chinese Journal of Rock Mechanics and Engineering，2000，19(5)：643–646.(in Chinese))
[13]	李腊梅，冯春. 一种非连续介质中热传导过程的数值模拟方法[J]. 工程力学，2016，33(1)：25–31.(LI Lamei，FENG Chun. A numerical simulation method for heat conduction in discontinuous media[J]. Engineering Mechanics，2016，33(1)：25–31.(in Chinese))
[6]	张杨，强晟. 非稳定温度场的三维p型有限元方法研究[J]. 岩土力学，2009，30(2)：487–491.(ZHANG Yang，QIANG Sheng. Research on 3D p-version hierarchical FEM for unsteady temperature field[J]. Rock and Soil Mechanics，2009，30(2)：487–491.(in Chinese))
[26]	张勇慧. 关于DDA中推导热应力问题公式的讨论[J]. 岩石力学与工程学报，2001，20(2)：277–278.(ZHNAG Yonghui. Discussion of the thermal stress problem in DDA method[J]. Chinese Journal of Rock Mechanics and Engineering，2001，20(2)：277–278.(in Chinese))
[3]	王贤能，黄润秋，黄国明. 边坡岩体浅层破坏的热应力效应研究[J]. 工程地质学报，1997，5(3)：262–268.(WANG Xianneng，HUANG Runqiu，HUANG Guoming. Study on the thermal stress effect on large failure of shallow of slope rock mass[J]. Journal of Engineering Geology，1997，5(3)：262–268.(in Chinese))
[8]	王润英. 计算三维温度场及温度应力的径向点插值无网格法[J]. 三峡大学学报：自然科学版，2009，31(2)：18–22.(WANG Runying. Radial point interpolation meshfree method for calculating 3D temperature field and thermal stress[J]. Journal of China Three Gorges University：Natural Sciences，2009，31(2)：18–22.(in Chinese))
[10]	李义. 基于独立覆盖数值流形法的大体积混凝土温度场仿真计算[硕士学位论文][D]. 武汉：长江科学院，2015.(LI Yi. Simulation of mass concrete temperature field based on the numerical manifold method with independent cover[M. S. Thesis][D]. Wuhan：Changjiang River Scientific Research Institute，2015.(in Chinese))
[16]	YANG Y，XU D，ZHENG H. Explicit discontinuous deformation analysis method with lumped mass matrix for highly discrete block system[J]. International Journal of Geomechanics，2018，18(9)：04018098.
[18]	JIANG Q H，YEUNG M R. A model of point-to-face contact for three- dimensional discontinuous deformation analysis[J]. Rock Mechanics and Rock Engineering，2004，37(2)：95–116.
[20]	MORGAN W E，ARAL M M. An implicitly coupled hydro- geomechanical model for hydraulic fracture simulation with the discontinuous deformation analysis[J]. International Journal of Rock Mechanics and Mining Sciences，2015，73 ：82–94.
[23]	JIAO Y Y，ZHANG X L，ZHAO J. Two-dimensional DDA contact constitutive model for simulating rock fragmentation[J]. Journal of Engineering Mechanics，2011，138(2)：199–209.
[4]	朱伯芳. 大体积混凝土温度应力与温度控制[M]. 北京：中国电力出版社，1999：122–141.(ZHU Bofang. Thermal stresses and temperature control of mass concrete[M]. Beijing：China Electric Power Press，1999：122–141.(in Chinese))
[14]	严成增. FDEM-TM方法模拟岩石热破裂[J]. 岩土工程学报，2018，40(7)：1 198–1 204.(YAN Chengzeng. Simulating thermal cracking of rock using FDEM-TM method[J]. Chinese Journal of Geotechnical Engineering，2018，40(7)：1 198–1 204.(in Chinese))
[24]	徐栋栋，邬爱清，李聪，等. 破裂全过程模拟的改进非连续变形分析方法[J]. 岩土力学，2019，40(3)：1 169–1 178.(XU Dongdong，WU Aiqing，LI Cong，et al. An improved discontinuous deformation analysis method for simulation of whole fracturing process[J]. Rock and Soil Mechanics，2019，40(3)：1 169–1 178.(in Chinese))
[28]	JIAO Y Y，ZHANG X L，ZHANG H Q，et al. A coupled thermo-mechanical discontinuum model for simulating rock cracking induced by temperature stresses[J]. Computers and Geotechnics，2015，67 ：142–149.