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.

A nonlinear parallel-bonded stress corrosion (NPSC) model is proposed to simulate the fatigue characteristics of artificial rock (concrete) during cyclic loading. Numerical simulations of fatigue tests replicate the main mechanical features of concrete specimens subjected to cyclic loading observed in the laboratory. A nonlinear reduction speed of the bond diameter between two bonded particles represents the damage rate induced by the fatigue load. The damage rate is proportional to the maximum cyclic load level when the minimum cyclic load level is fixed. Compared with laboratory data, a logarithmic function of bond diameter in the NPSC model resulted in the best fit to simulate the fatigue behaviour of concrete. The simulation includes acoustic emission (AE) monitoring during fatigue tests. The axial strain of the assembly is governed by the evolution of bond breakages. The sum of released bond strain energy is documented as value proportional to cumulative AE energy. The simulation results show very similar evolution compared with laboratory data, which verifies the effectiveness of AE energy simulation.