(Key Laboratory of Geotechnical Mechanics and Engineering of Ministry of Water Resources，Changjiang River Scientific Research Institute，Wuhan，Hubei 430010，China)
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.
[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))
[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))
[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))
[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))
[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))
[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))
[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))
[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))
[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))
[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))
[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.
[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))
[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))
[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))
[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.
[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.
[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))
[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.
[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.
[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.
[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.
[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.
[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))
[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))
[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))
[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.
[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.
[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))
2019, Vol. 38 
