•We determine tunnel stability in soil with a soft upper layer and hard lower layer.•The effects of multiple soil layers and different soft soil parameters on an underground tunnel in China are ...studied.•We present the upper-bound finite element method with a mesh adaptive strategy.•Two main reasons for tunnel instability and large surface settlement are presented.•Two construction proposals are presented to ensure normal working conditions and safety of engineering.
This paper presents a limit analysis of a tunnel in non-homogeneous soil using the upper-bound limit analysis in combination with finite element technology and a mesh adaptive strategy by investigating the critical loads and failure mechanisms of a tunnel cross-section and tunnel face. A case study of an underground tunnel in China is conducted, with a focus on the effects of multiple soil layers and different soft soil parameters. This work demonstrates that the variation in the soil layering plays an important role in the initiation of tunnel collapse, possibly shifting the failure region and even changing the failure mechanism. When the soil changes from a strongly weathered conglomerate to a soft soil, not increasing the low tunnel pressure quickly enough and using the wrong excavation mode are the two main reasons for tunnel instability and large surface settlement.
Piles are one-dimensional (1D) beam–column elements used to support tall buildings, onshore and offshore structures, bridges, wind turbines, etc. In practice, piles are commonly designed using the ...Euler–Bernoulli Theory. However, this theory is only considered appropriate for long, slender piles because it neglects the contribution of the shear deformation on the pile response. The effect of the shear deformation might be of importance in short, rigid and large-diameter piles such as monopiles, and it should be considered in the analysis. In geotechnical engineering, studies addressing this effect and its interaction with the surrounding soil are scarce and very limited. The aim of this work is to conduct a parametric study to investigate the influence of the shear deformation, semi-rigid connections and pile/soil stiffness ratio on the head deflection, head rotation and first-mode of critical load of an embedded pile. For this purpose, the well-known commercial FE software SAP2000 is employed. Four examples are presented to show the effect of the shear deformation on reinforced concrete (RC) and hollow circular piles embedded in either homogeneous or non-homogeneous soil. For lateral deformation analysis, it is found that deflections and rotations at the pile head considerably increase when the shear deformation is taken into account and they are of most significance at low values of pile/soil stiffness ratios. For buckling analysis, Timoshenko’s theory shows that the critical load notably increases as the stiffness of the surrounding soil increases and, regardless of the soil/pile stiffness ratio, it rapidly decreases as the shear stiffness of the pile decreases.
•The effects of shear deformation on the lateral response of piles are investigated.•The effects of shear deformation on the critical load of piles are examined.•The influence of the soil inhomogeneity in conjunction with the effect of shear deformation is examined.•Ranges of pile/soil stiffness ratios where shear deformation are of significance are provided.•Results from Timoshenko’s and Euler–Bernoulli’s beam theories are discussed.
The dynamic response of a three-dimensional multi-layered pavement system subjected to a moving load is calculated by an analytical solution method. In this method, the pavement structure is assumed ...to be a transverse isotropic viscoelastic multi-layered structure with partially bonded interface, and the subgrade is considered to be a transverse isotropic half-space medium with its modulus varying with depth. The load is considered as a three-dimensional elliptical moving load with a constant speed. This study is divided into three key steps. Firstly, the ordinary differential equations of pavement and subgrade structures are obtained via the double Fourier transform, and the general solutions of dynamic response are obtained according to the theory of ordinary differential equations and the Frobenius method. Secondly, combined with the boundary conditions, the method of wave vector matrix, and the inverse Fourier integral transform, the dynamic responses of pavement structure in time domain are calculated. The accuracy of the analytical solution is verified by comparing against the solution from a finite element method. Lastly, the sensitivity analysis is conducted to assess the influence of four parameters on the dynamic response, namely, the horizontal load, the condition of interlayer contact, the transverse isotropy, and the non-uniform distribution of subgrade moduli. The calculation methods and the corresponding results can further develop a theory of the multi-layered system.
•Propose a fast method for calculating the dynamic response of pavement structure.•Consider the non-uniform distribution of subgrade modulus along depth.•Comprehensively investigate the effects of load and interlayer contact conditions.
•Beam-column elements with generalized end-boundary conditions in Pasternak soil.•Differential Transformation Method (DTM) applied to solve the governing equation.•Flexibility to incorporate ...semirigid connections and lateral restraints.•Comparison with other methods to highlight accuracy and implementation.
A simplified analytical approach to analyze soil-structure interaction of beam-column elements (i.e., beams, columns, piles, etc.) with generalized end-boundary conditions on a homogeneous or non-homogeneous Pasternak foundation is developed. The Differential Transformation Method (DTM) is employed to solve the equation that governs the mechanical response of the proposed beam-column element. Using the same set of expressions, the proposed formulation can be used to conduct both static and stability analyses. It also provides the flexibility to incorporate the effects of semirigid connections, lateral restraints, and axial and lateral load at the ends of the element. The effects of any slenderness ratio and element/soil stiffness ratio are captured with the proposed formulation. Five examples are examined to illustrate the simplicity of the method to conduct static and stability analyses. The proposed method is compared with more rigorous, but cumbersome, methods to highlight the accuracy and relative ease of implementation.
A simplified method to study the response of embedded piles in a non-homogeneous Pasternak elastic soil is developed. The governing differential equation (DE) derived from equilibrium is solved using ...the Differential Transformation Method (DTM), providing a compact solution for static and stability analyses that is easy to implement. The solution derived with the DTM overcomes the limitation of conventional approaches used to solve the proposed structural model, where the mathematical derivation and solution become cumbersome and tedious to implement. The proposed formulation is a continuation of the work recently presented by the authors for frictionless beam-column elements, and now incorporates the magnitude and distribution of friction along the element. The applied axial load is carried by both end-bearing resistance and skin friction, so that purely end-bearing, partially frictional, and purely frictional piles can be studied. In practice, the influence of side friction on pile stability is typically neglected, but can have a substantial influence on the buckling capacity of embedded elements---a practical benefit that has received little attention. This work provides an analytical solution for the governing DE with the flexibility to account for different distributions and magnitudes of skin friction, different combinations of geotechnical end- and side-resistance, and generalized end-boundary conditions. The model is applied to elucidate the effects of intermediate end-boundary conditions, as well as the distribution and magnitude of soil-pile interface friction on pile stability.
Heat transfer through the soil is important in shallow geothermal applications, in plant growth as well as in the surface energy balance. In this work, an analytical solution for a ...conduction–convection equation that models heat transfer is presented for non-homogeneous soil. The main assumption underlying the solution is that the thermal diffusivity of soil is piecewise-constant. A method to determine the depth-dependent thermal diffusivity and the water flux density through the porous medium based on temperature measurements is also developed and applied to field data from a site in Ioannina, Greece. The thermal diffusivity was found to increase in the layer below the surface and decrease for larger depths, while convection through the porous medium was found to be present in wet conditions and to account for about 10% of the heat flux in terms of the annual variability of temperature. Finally, the capability of the novel method in capturing the spatio-temporal variability of soil temperature is compared to three commonly used algorithms: the amplitude, the phase, and the conduction–convection algorithm. The novel method is able to reduce the root-mean-square error for the predicted variability of temperature at all depths by an order of magnitude compared to the other three algorithms.
•A new model for wave-induced seabed response considering 3D non-homogeneous soil properties.•Effects of the horizontally non-homogeneous soil property on the local seabed response.•Wave-induced ...dynamic response and liquefaction potential of the 3D non-homogeneous seabed around a mono-pile.
Marine seabed is often composed of non-homogeneous soil in various directions, where its response to dynamic wave loading shows much difference from the case of homogeneous soil. Most existing numerical models for wave-induced non-homogeneous seabed response have been limited to vertically layered seabed. In this study, a new model is developed to simulate wave-induced seabed response considering non-homogeneous soil properties in three dimensions. The present model (1) directly solves the original fully dynamic governing equations for the overall equilibrium of soil, the equilibrium of pore fluid flow and the mass balance of porous seabed, (2) allows non-uniform distribution of eight soil parameters in the three-dimensional space, and (3) retains all the spatial derivative terms of soil parameters to fully represent the three-dimensional non-homogeneous seabed behavior. The model well reproduces the existing experimental data, numerical solution and field observation. Two preliminary model tests are performed to investigate the effects of horizontally non-homogeneous soil property on the local seabed response, and the three-dimensional non-homogeneous seabed response around a mono-pile. It is found that the presence of coarser soil will reduce the liquefaction depth of the adjacent finer soil. This impact is more pronounced for larger wave period, wave height, soil permeability, saturation degree, and smaller water depth. The seabed response around a mono-pile can be considerably modified if the soil sorting due to long-time lateral rocking of mono-pile is taken into account, which decreases the liquefaction depth in the vicinity of the mono-pile.
It is important to ensure that a grounding system is designed with a low magnitude of earth resistance, so the protection system can divert the large fault current to earth effectively. The ...performances and protection level of a grounding system need to be acknowledged as the condition of soil structure changes with different soil characteristics. At present, there is a lack of systematic guide or standards for grounding grid designs that consider non-uniform soil and its impact on the grounding systems. By computing the grid safety threshold parameters consisting of the grid impedance, step, and touch voltages, a comparison has been made between uniform soil and two-layer soil models. Where the competence and level of safety of the grounding systems depend on the soil attributes, the significant impact of various soil conditions is seen. The evaluations on performance and safety assessment in two-layer soil conditions hold the novelty and originality as there is no such comparison and discussion have been made to date. These comparisons would help in forecasting the behavior and safety of the grounding system in various soil environments, which would provide engineers with additional expertise to design an effective and secure grounding system. This research would contribute to the existing body of knowledge by differentiating and predicting the performance of a grounding system when the characteristics of the soil differ significantly from uniform soil as most of the standards and guidelines only consider uniform soil while designing a grounding system, owing to its complexity at the site.