Thermo-mechanical coupling response of a layered isotropic medium around a cylindrical heat source

Based on the governing equations of the thermo-elastic problem, the analytical layer-elements of a finite layer and an underlying half-space are obtained using the Laplace-Hankel transform and the characteristic value method. The cylindrical heat source is divided into several micro cylindrical heat sources, which can be approximately simulated by plane heat sources. Then, the global stiffness matrix for the problem is assembled and solved in the transformed domain, and a Laplace-Hankel transform inversion is taken to obtain the real solution. Finally, the influence of heat source types, division numbers, embedded depths and layered properties on the thermo-mechanical coupling response is investigated.     

Consolidation around a heat source buried deep in a porous thermoelastic medium with anisotropic flow properties

A solution is derived for the heat flow and consolidation which occur when a heat source is buried deep in a porous thermoelastic soil having anisotropic flow properties. This solution is used to examine the pore pressure generation and dissipation near both point and cylindrical heat sources. An increase in temperature will tend to generate an increase in excess pore pressure. However, the pore water will tend to flow from regions of high excess pore pressure to regions of low excess pore pressure, and so consolidation will occur, and temperature-generated excess pore pressures will tend to dissipate. Many natural soils exhibit horizontal layering and so have a higher horizontal than vertical permeability. It is shown that in soils the excess pore pressure generated by a heat source is significantly less than that in an isotropic soil having an equal vertical permeability.     

Thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source

We develop a new numerical model based on a precise integration method to investigate the coupled thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source. To obtain the relational matrices of the extended precise integration method, we first convert the governing equations of the problem into ordinary differential matrix equations through the Laplace–Hankel transform. Then, the cylindrical heat source is divided into a series of plane heat sources, and the plane temperature load term is added to the state vector between layer elements. By combining the layer elements, we build a layered transversely isotropic numerical model containing a cylindrical heat source in the transformed domain. Finally, we solve the model in the transformed domain and obtain the solution of the problem in the real domain through the Laplace–Hankel transform inversion. The accuracy of this method is verified by comparing the solutions with the results of the analytical method and the finite element method. Then, we study the influence of the anisotropy of thermal parameters, the embedded depth, the length/radius ratio, the type of heat source and the stratification of the medium on the thermo-mechanical coupled performance.

     

Consolidation around a point heat source

When a heat source such as a canister of radioactive waste is buried in a saturated soil the temperature changes that occur will cause the pore water to expand a greater amount than the voids of the soil. The temperature change will thus usually be accompanied by an increase in pore pressure. If the soil is sufficiently permeable these pore pressures will dissipate. This paper develops an analytic solution for the fundamental problem of a point heat source buried deep in a saturated soil.     

地埋点热源作用下层状半空间的热-力耦合响应

采用解析层元法对存在地下点热源的岩土工程问题进行解答。首先从热弹性力学三维问题的基本控制方程出发,利用拉普拉斯-傅里叶积分变换推导出其在变换域内单层介质及下卧半空间的解析层元;然后结合有限单元法原理组装得到总刚度矩阵,结合边界条件,得到其在变换域内的解答,最后应用拉普拉斯-傅里叶积分逆变换技术,得到物理域内的解。编制了相应的计算程序,算例分析表明：该结果与已有文献吻合较好,该方法求解层状半空间的热-力耦合响应问题具有较好的适用性和较高的精度;层状岩土介质体系中,热扩散系数对温度及地表隆起的变化过程影响显著,但对其初始值和最终稳定值影响不明显;分层特性对岩土介质温度分布及地表位移变化过程均有显著影响。     

The time–deflection behaviour of a rigid under-reamed anchor in a deep clay layer

The behaviour of an under-reamed anchor in a consolidating soil is examined by approximating the anchor as an impermeable circular plate embedded in a deep soil layer. Hankel and Laplace transforms are applied to the equations governing the consolidation process, and this greatly simplifies the equations, allowing a solution to be obtained for the transformed variables. Numerical inversion of both the Hankel and Laplace transforms is used to obtain the solution at any time. A particular feature of the solution method is that the unknown contact stress between the anchor and the soil and the unknown flows in the plane of the anchor are approximated as a series of simple functions with unknown coefficients. By determining the coefficients of the terms in the series, the complete solution may be found. Computations have been carried out using the method proposed, and results are presented for the time–settlement behaviour of an impermeable anchor. These results are compared with some published and some recomputed finite element data, and this highlights some of the difficulties encountered in using such numerical techniques.     

Thermomechanical effects on adsorbed water in clays around a heat source

Numerical results of the simulations of the saturated clay response to a line heat source with a constant power output are presented and compared to the analytical solutions. Mixture theory of two interacting thermoelastic constituents is employed for the modelling. In this theory, mass transfer between the adsorbed water and bulk water is included to simulate the degeneration of adsorbed water into bulk water occurring at elevated temperatures. A finite difference method with non-uniform mesh is adopted to determine both the near- and far-field behaviour. The degeneration of adsorbed water into bulk water is more pronounced close to the heat source. This may produce a lower excess pore water pressure, a higher permeability increase, and accelerate the water flow in the vicinity of the heat source.     

Consolidation around a spherical heat source with a decaying power output

When a heat source such as a cannister of radioactive waste is placed in a saturated soil the temperature of both the soil and the pore water will rise. The coefficient of expansion of the pore water is greater than that of the soil skeleton and so there will be a tendency for a differential volume change and a consequent increase in pore pressure. This will in turn lead to a reduction in effective stress and so if the temperature is too great cracking of the surrounding soil may occur.If the soil is sufficiently permeable the temperature generated pore pressures will dissipate during consolidation and so the severity of these effects will be reduced. A further reduction will also occur because the power output of the radioactive waste will decrease with time due to radioactive decay.This paper develops an analytic solution for the temperature and pore pressure field around a rigid impermeable spherical heat source which takes account of the reduction of power output due to radioactive decay.     

Post-installation pore-pressure changes around spudcan and long-term spudcan behaviour in soft clay

This paper presents a dual-stage Eulerian–Lagrangian analysis for modelling the entire process of spudcan installation in soft clay, followed by consolidation and working load operation. The analysis consists of three components, namely undrained effective stress Eulerian analysis of spudcan installation, mesh-to-mesh variable mapping and coupled-flow Lagrangian analysis for the post-installation spudcan working behaviour. The results show good agreement with centrifuge model data but also highlight the importance of replicating the hysteretic behaviour of the soil. The findings also show that while a wished-in-place approach was able to model the long-term bearing response of the spudcan, rotational stiffness was over-estimated. This is due to the fact that, while the wished-in-place analysis was able to model the hardening of the soil ahead of the spudcan, it was unable to model the softening of back-flowed soil behind spudcan. The latter influences the spudcan fixity significantly, but not bearing response. Although the analyses were conducted using ABAQUS, they can, in principle, be conducted using other codes.     

The behaviour of layered soil or rock containing a decaying heat source

Solutions are presented for the behaviour of layered soil or rock deposits which contain a heat source. Such a problem arises when high level nuclear waste is placed in deep underground depositaries, as the waste continues to generate heat for many years after placement. This heating of the surrounding soil or rock may lead to expansion and cracking with subsequent contamination of ground water. Results are presented for heat soureces with different decay rates and for heat sources in layers of material with different coefficients of expansion. An example using realistic data for rock is also given. The solution method involves applying Fourier or Hankel transforms to the field quantities and this reduces the two-dimensional or axisymmetric problem to one involving a single spatial dimension. In cases where the soil or rock is horizontally layered, the method has great advantages over other numerical methods such as finite element or finite difference techniques, since little computer storage and data preparation time is required. Solution of the time-dependent problem is carried out by applying Laplace transforms to the field variables, obtaining solutions and then using numerical means to invert the transformed solutions. This enables easy solution of problems involving time-dependent (i.e. decaying) heat sources.     

增强地源热泵竖直埋管地下换热器换热性能的研究

地下换热器换热性能直接影响到地源热泵系统的初投资和运行成本.地下换热器可分为单U型、双U型和1 2型三类.采用非稳态法对膨润土、水泥 砂两类回填材料进行了室内实验;在相同的边界条件下,采用扶正器和定位器施工可以保证下管,并能降低钻孔内热阻.因地制宜地选择埋管形式,采用稳定、高效的回填材料和减小热短路是提高地下换热器换热性能的必要途径.     

非饱和上海软土水力和力学特性耦合弹塑性模拟

目前大多数非饱和土的弹塑性本构模型用非饱和击实土的试验结果进行验证,但现场其他类型的土,如沉积土经常有在非饱和状态下外部环境变化的情况。现有的非饱和土弹塑性模型是否适用于沉积土一类的现场土是需要研究的课题。进行非饱和上海第③层土的吸力控制排水排气三轴剪切试验,使用文中提出的能统一考虑非饱和土水力性状和力学性状的弹塑性本构模型,预测上述三轴试验结果,并与试验数据进行比较。比较结果显示,建立的本构模型能够很好地预测非饱和上海软土的水力和力学性质,说明该模型不仅可以适用击实土的预测,还能够很好地适用于其他类型非饱和土的水力和力学性质的模拟。     

Scoping analyses of the coupled thermal-hydrological-mechanical behaviour of the rock mass around a nuclear fuel waste repository

Scoping calculations were performed in order to assess the influence of radiogenic heat on the performance of the rock mass around a nuclear fuel waste repository. The full coupling between the thermal, mechanical and hydrological processes involved was considered by using the finite element code, FRACON, developed through an extension of Biot's classical theory of soil consolidation. By considering the full T---H---M coupling, several important safety features, which would otherwise be omitted in uncoupled analyses, were detected in the present study. In particular, it was shown that the heat-induced pore pressure increase around the repository has the potential to significantly increase the rate of groundwater flow, and affect the structural integrity of the rock mass.     

基于耦合欧拉-拉格朗日法的锚板在黏土中的极限承载力数值分析

当结构在土体中运动时,往往导致土体发生较大的变形,此类问题采用大变形数值分析方法更为恰当。耦合欧拉-拉格朗日（Coupled Eulerian-Lagrangian, 简称CEL）法是大变形数值分析方法中的一种,在分析大变形问题时具有很强的适用性,但在国内尚未开展CEL法分析锚板承载力的研究。以方形锚板在均质土及线性土中的拔出试验为原型,基于CEL法建立数值模型,对锚板的极限承载力及破坏机制进行研究,并通过用户自定义子程序,实现了线性土的强度分布随锚板拔出而变化。计算结果表明,土体杨氏模量越大,锚板的极限承载力越大;随着位移增大,锚板的抗拔力先增大,后降低;当埋深小于临界埋深时,土体发生整体破坏;当埋深大于等于临界埋深时,土体发生局部破坏。数值计算反映的规律与试验结果基本吻合,体现了CEL法模拟锚板在海床中大位移响应的出色能力。     

严寒地区地源热泵地埋管周围土壤冻结影响因素的实验研究

为了研究影响地埋管周围土壤冻结的因素,文章建立了模拟严寒地区热失衡状态下地源热泵冬季运行情况的实验装置。实验研究了进口流体温度、土壤体积含水率和土壤初始温度对地埋管周围土壤冻结的影响。结果表明:进口流体温度变化改善了埋管周围的温度场,有利于埋管换热器换热;随着土壤初始温度升高,减小了土壤物性参数的变化,增加了地埋管换热器的换热量;土壤体积含水率的增加,强化了周围土壤与地埋管换热器的换热作用,有利于提高地埋管的出口流体温度和取热值。     

Upper bound solution of a laterally loaded rigid monopile in normally consolidated clay

A centrifuge test is performed to investigate the ultimate lateral capacity of a free-head rigid monopile foundation in normally consolidated (NC) clay. Based on the test result, an upper-bound velocity field is constructed to illustrate the evolution of the soil resisting mechanisms with the progressive rotation and the effect on reaction forces during loading. Simultaneously finite element analyses are also employed to model the centrifuge test and verify the conclusion of the upper bound analysis. Finally, an empirical expression is presented for predicting the upper-bound collapse loads of monopiles with different eccentricity/length ratios and length/diameter ratios in NC clay.     

Stress analysis of a deep rigid axially-loaded cylindrical anchor in an elastic medium

The elastostatic problem of an infinite elastic medium containing an axially-loaded rigid cylindrical inclusion is investigated. This problem is of interest in connection with the geotechnical study of the time-independent, load-deflection characteristics of deep rigid anchors embedded in cohesive soil or rock media. The problem is formulated by means of Hankel integral transforms and reduced to a system of four coupled sing ular integral equations, where the unknown quantities are the normal and the shear stresses acting on the entire surface of the anchor. Numerical solutions are investigated for various Poisson's ratios and several values of the aspect ratio of radius to length of the cylindrical anchor.     

多向荷载作用下层状地基刚性桩筏基础分析

提出一种多向荷载作用下层状地基中刚性桩筏基础的计算方法。基于剪切位移法,采用传递矩阵形式分析了竖向荷载下桩顶面-桩顶面相互作用;引入修正桩侧地基模量,采用有限差分法分析了水平荷载下桩顶面-桩顶面相互作用;基于层状弹性半空间理论,分析了多向荷载下桩顶面-土表面、土表面-桩顶面、土表面-土表面的相互作用关系。建立了桩土体系柔度矩阵,得到了多向荷载下层状地基中刚性桩筏基础的受力和变形的关系以及桩的内力和变形沿桩身分布规律。通过与有限元对比,验证了该方法的合理性和修正地基模量的优越性,并对多向荷载作用下的桩筏基础进行了计算分析,计算结果表明,水平力将会引起桩筏基础的倾斜。     

New method for non-linear analysis of rigid piles in clay

Using the numerical technique of solving nonlinear simultaneous equations, a new method of analysis is developed for design of eccentric and/or inclined loaded rigid piles in clay. The method takes account of the non-linear strain-stress characteristics of the soil around the piles and can be used for round, rectangular and other symmetric rigid piles.

The computed results using this method were found to be in good agreement with the measured values of some field tests and laboratory model tests.