The engineering application of cylindrical heat source usually influences the characteristics of the surrounding soil and induces significant coupling responses of temperature, seepage and stress fields. This paper presents a numerical investigation into the thermo-hydraulic-mechanical (THM) coupling behavior of layered soils surrounding a cylindrical heat source by combining the analytical layer element method (ALEM) and the finite element method (FEM). The FEM is utilized to analyze the soil-cylindrical heat source interaction and the ALEM is used to obtain the fundamental solution of the layered soil. Then the THM coupling solution between the heat source and the soil is obtained by the simplified FEM-ALEM coupling method. Detailed validation against experimental result, and verifications against semi-analytical solution and numerical results by Comsol Multiphysics are performed to confirm the robustness of the present method, followed by extensive parametric studies to examine the effects of decaying half-life and type of the time-varying heat source, and the soil layered characteristics on the behavior of soils. It’s evident that the present method is accessible and can be implemented via common programming language, e.g. Fortran.

圆柱形热源的工程应用通常会影响周围土壤的特性，并引起温度、渗流和应力场的显着耦合响应。本文通过结合分析层元法 (ALEM) 和有限元法 (FEM) 对圆柱形热源周围的层状土壤的热力水力 (THM) 耦合行为进行了数值研究。有限元法用于分析土壤-圆柱形热源相互作用，而ALEM用于获得层状土壤的基本解。然后通过简化的FEM-ALEM耦合方法得到热源与土壤的THM耦合解。针对实验结果的详细验证，Comsol Multiphysics 对半解析解和数值结果进行了验证，以确认本方法的稳健性，随后进行了广泛的参数研究，以检查衰减半衰期和时变热源类型的影响，以及土壤层状特征对土壤行为的影响。很明显，本方法是可访问的，并且可以通过通用编程语言（例如 Fortran）来实现。