Research has been carried out to study the effects of new tunnelling on an existing adjacent tunnel to ensure the safety and serviceability of tunnels. Prior studies on twin-tunnel interaction have mostly centred on simplifying perpendicularly crossing tunnelling in a single-layered soil stratum. New tunnel excavation beneath an existing tunnel at different skew angles in two-layered strata can lead to different patterns of stress redistribution and adverse impacts on the existing tunnel. In this paper, results of three-dimensional centrifuge and numerical modelling carried out to study the twin-tunnel interaction with varying advancing orientations and layered soils will be reported. The influence of new tunnel excavation on an existing tunnel was simulated in-flight by controlling both the tunnel weight and volume losses. An advanced hypoplastic constitutive model that can capture stress-, path, and strain-dependency of soil behaviour is utilised for numerical back-analyses and parametric studies. Cases investigated include twin-tunnel interaction at three different skew angles (30 degrees, 60 degrees, 90 degrees) in a uniform sand layer and at skew angle of 90 degrees in two-layered sand with different relative densities and thicknesses. Distinct load redistribution patterns will be presented to explain deformation mechanisms of the existing tunnel at different tunnel advancing skew angles to highlight the effects of tunnelling orientation. The results of perpendicularly crossing tunnelling in twolayered sand will also be reported and compared to reveal the influence of layered soil. The findings and new insights can help engineers better estimate advancing tunnelling effects on existing tunnels and enhance the safety of tunnel construction.

Energy shallow foundations represent an innovative technology that can simultaneously support structural loads and harvest geothermal energy. During geothermal operations, the underlying soils are subjected to structural loads and temperature fluctuations. Despite the potential, knowledge regarding the thermo-hydro-mechanical behavior of the multilayered soils beneath the energy foundations remains scarce. This study proposed an analytical approach to investigate the thermo-hydro-mechanical response of soft fine-grained soils beneath energy shallow foundations. The analysis focused on the evolutions of the temperature, pore water pressure, and vertical displacement of the underlying soils. The results indicate that the generation and development of the thermally induced excess pore pressure are controlled by thermal transfer processes and soil hydraulic properties. Furthermore, the mechanical load-induced ground settlement decreases upon heating and increases upon cooling, primarily due to the development of thermally induced pore pressure and the thermal volume changes of the soil skeleton. Under the considered conditions, ignoring the thermally induced mechanical effects could result in a settlement prediction error of nearly 120%. Therefore, the thermo-hydro-mechanical interactions within the soils should be appropriately considered in the analysis and prediction of the displacement behavior of the energy foundations.

Suction bucket jackets have been used as foundations for offshore wind turbines in intermediate water depths where layered soil stratigraphies are often encountered. Although suction installation in layered soils has been studied, experimental data on the in-service response is scarce. During installation in stratigraphies containing a low permeability layer underlain by a high permeability layer, suction is transferred to the underlying layer when the pressure at the lid invert is sufficient to uplift the low permeability plug. This suction-transfer mechanism also affects the in-service response, albeit the load-sharing mechanism is not well understood. This paper presents data from centrifuge tests of suction buckets subjected to constant amplitude and varying amplitude cyclic vertical loading in two stratigraphies-a sand with an overlying clay layer and in a sand with a sandwiched clay layer. These experiments show that tensile stresses exceeding the vented tensile resistance can be withstood without significant uplift of the bucket in both stratigraphies, even under a zero mean stress. Plug uplift was shown to have an important effect on the amount of stress transferred to the skirts, with the load-sharing mechanism depending on the stratigraphy. Additionally, the load-sharing mechanism and the bucket in-service resistance was shown to depend on the effectiveness of the clay in sealing the soil plug within the bucket, with a more effective seal resulting in higher tensile resistance and therefore better performance. A limiting loading condition was not identified in the sand with a sandwiched clay layer, with the data indicating that the suction pressure to cause plug uplift during cyclic loading may be much higher than during suction installation.

In this study, a theoretical approach is presented for analyzing how rectangular barrettes respond laterally in layered transversely isotropic soil deposits. To do this analysis, a modified Vlasov model is used. In this study, the barrette and the soil around it are treated as a continuum system. The deformation of the barrette is analyzed using the Timoshenko beam theory. By multiplying the barrette's displacement with a pair of decay functions, the horizontal soil displacement can be quantified. The equations that govern the barrette and soil are derived based on the principle of minimum energy, along with the appropriate boundary conditions. These equations are then solved using an iterative method. The accuracy of the results is confirmed by comparing the barrette response to two previously published results. Additionally, the impact of the shape of the rectangular cross and the anisotropy of the soil on the static responses of a barrette are explored. The results show that the ratio Esh/Esv between the horizontal modulus and vertical modulus for the transversely isotropic soil has significant influences for the response of barrette. An increase of Esh/Esv from 0.5 to 3.0 can lead to a reduction of around 75%, 54%, 30%, 40% for the maximums of lateral displacement, rotation, moment, and soil reaction, respectively.

Behind a retaining wall, the mean effective stress mainly decreases during an excavation phase following an unloading stress path. The volumetric strains generated by purely elastic soils are systematically dilative which induces aberrant ground uplifts. The introduction of plasticity along with a nonlinear elastic domain turns out to be essential for a realistic prediction of ground movements. In this paper, a numerical analysis is carried out using a finite element code considering an advanced soil constitutive model called Generalized Hardening Soil which has been recently developed. This model contains the exact same set of features as the Hardening Soil Small Strain model but with the possibility to activate each of its plastic and nonlinear elastic mechanisms independently. The role of these mechanisms are investigated to assess their impact on the shape and the amplitude of the ground movements. Numerical results demonstrated that plasticity triggers the main contractive volumetric strains leading to settlements. Nevertheless it cannot fully compensate the elastic uplifts due to unloading. The insertion of strain dependent stiffness is essential as well as the stress dependency. A back analysis of the historical excavation of the Taipei National Enterprise Center permitted to validate these findings.

This paper investigates the consolidation behavior of multi-layered viscoelastic soils considering groundwater. First, the fractional Merchant viscoelastic model is introduced to describe the behavior of multi-layered viscoelastic soils considering groundwater. Later, the governing equations are extended to a viscoelastic medium by virtue of the elastic-viscoelastic corresponding principle in the Laplace-Hankel domain. According to the extended precise integration method, the soil layer is divided into a series of layer units. Then the relationship between general stress vector and general displacement vector on the top and bottom planes is established. Every two adjacent layer units are combined into one layer in each computational iteration. The solutions in the Laplace-Hankel domain are obtained by considering the boundary conditions, and numerical inversion is performed to obtain the solutions in the physical domain. The practicability of the present method is assessed by comparing the numerical results with those in the existing literature and done by ABAQUS. Finally, the effects of groundwater table, properties of the soils above groundwater table, load depth, viscoelastic parameters, and soil stratification are investigated.

Prefabricated vertical drains (PVDs) combined with vacuum and/or surcharge loading have been widely adopted to improve the strength of soft soils. Precise consolidation analysis is the theoretical basis for the design of preloading method with PVD. Current consolidation theories for layered soils with PVD seldom consider the influence of large strain, nonlinear creep, and self-weight loading simultaneously. This paper, thus, presents a finite strain elastic visco-plastic consolidation model, called RCS-EVP, for radial consolidation of layered soils with PVD. RCS-EVP is developed based on the piecewise-linear method. It takes into account nonlinear creep with limit creep strain, variable boundary conditions, anisotropy of soil hydraulic conductivity, and variable compressibility and hydraulic conductivity during the consolidation under self-weight, time-dependent surcharge and/or vacuum loading. The performance of RCS-EVP is evaluated by comparing with the results from finite element simulations and a laboratory physical model test. The variations of settlement and pore pressure of a soft soil ground improved by vacuum preloading with PVD are estimated using RCS-EVP. The results indicate that RCS-EVP provides good estimates of long-term consolidation of layered soils with PVD under both laboratory and in-situ conditions.

The mechanical response of energy pile groups in layered cross-anisotropic soils under vertical loadings is studied with the aid of the coupled finite element method- boundary element method (FEM-BEM). The single energy pile is simulated based on the finite element theory, which then is extended to energy pile groups. The global flexibility matrix for soils is obtained by considering the coupling effects of vertical and thermal loadings. The coupled FEM-BEM equation for the interaction between energy pile groups and soils is derived based on the displacement compatibility condition at the pile-soil interface. According to the displacement coordination condition and force balance in the rigid cap, the displacement of the cap and axial forces of pile groups can be solved. The presented theory is validated by comparing the calculated results with numerical simulations and field test results in existing literature. Finally, effects of the thermal loading, pile-soil stiffness ratio, pile spacing, cross-anisotropy of Young's modulus and the stratification are discussed.

Geo-materials naturally display a certain degree of anisotropy due to various effects such as deposition. Besides, they are often two-phase materials with a solid skeleton and voids filled with water, and commonly known as poroelastic materials. In the past, despite numerous studies investigating the vibrations of strip foundations, dynamic impedance functions for multiple strip footings bonded to the surface of a multi-layered transversely isotropic poroelastic half-plane have never been reported in the literature. They are first presented in this paper. All strip foundations are assumed to be rigid, fully permeable, and subjected to three types of time-harmonic loadings. The dynamic interaction problem is investigated by using an exact stiffness matrix method and a discretization technique. The flexibility equations are established by enforcing the appropriate rigid body displacement boundary conditions at each footing-layered soil interface. Numerical results for dynamic impedance functions of two-strip system are presented to illustrate the influence of various effects on dynamic responses of multiple rigid strip foundations.