Development Of Excess Pore Pressure Research Articles

Undrained cyclic responses of sandy soils improved by soybean-induced calcium carbonate precipitation

ABSTRACT This study investigates the behaviour of sandy soils under cyclic conditions, emphasising the efficacy of soybean-induced calcium carbonate precipitation (SICP) treatment in mitigating liquefaction. Analysis of shear wave velocity and undrained cyclic responses was conducted, with a focus on understanding deformation response and development of excess pore pressure. Results revealed CaCO3 content as a definitive indicator for improvements in soil stiffness, with its distribution influenced by soil’s relative density and number of SICP injections. Enhanced cementation efficacy was observed in denser or deeper layers of sand column, attributed to prolonged SICP retention. A marked correlation between shear wave velocity and CaCO3 content highlighted SICP’s role in enhancing the liquefaction resistance of loose sandy soils. Dynamic triaxial testing showed the nuanced strain behaviours across varying relative densities and cyclic stress ratio (CSR) values, emphasising the strengthening influence of CaCO3 formations in the soil matrix. The empirical equation served as a testament to the potential of SICP in soil stabilisation endeavours. The intervention of SICP treatment presents a potent tool for bolstering soil stability, especially in regions susceptible to seismic activities. The insights garnered advocate for further research aimed at harnessing and optimising this soil improvement technique on a global scale.

Numerical analysis of seismic response of a circular tunnel-rectangular underpass system in liquefiable soil

Earthquake-induced liquefaction can lead to uplift displacements of underground structures. While the behaviour of a single underground structure has been extensively studied, the coupled response of a shallow buried circular tunnel adjacent to an underpass with a rectangular section in liquefiable ground remains unclear. In this study, the finite difference method based on the fully coupled effective stress is employed to analyse the seismic response of two adjacent underground structures. Firstly, comprehensive parameters of Hostun sand for P2Psand constitutive model was calibrated. The accuracy of the numerical method is verified through a centrifuge test. Subsequently, comparisons are made between the seismic response of the tunnel-soil-underpass interaction and other cases: free field, tunnel, and underpass. The tunnel-soil-underpass interaction in liquefiable ground is analysed, with a special emphasis on the development of excess pore pressure, acceleration response of the soil, displacement and deformation of the structures. A section deformation coefficient was proposed to depict the tunnel deformation. Finally, horizontal spacing, buried depth of tunnel and underpass are discussed to clarify the interaction of two structures. The optimal relative position of two structures is suggested, providing valuable insights for the spatial design of underground structures in liquefiable ground.

Numerical And Physical Modelling Of The Pore Pressure Development Around A Monopile Foundation

Offshore wind has favoured the use of monopile foundation due to its simplicity in design, construction and industrial scalability. The stability of the monopile foundations can be affected not only by the direct action of wave loads but also by the response of the surrounding seabed. Numerical and physical modelling can be used to simulate the wave-structure-seabed interaction and accurately predict the wave-induced seabed response around the monopile foundation. Within this context, a 3D coupled numerical model is developed to investigate the excess pore pressure development around a monopile foundation and the accompanying changes in the effective stress of the seabed soil. In addition to the coupled hydrodynamic-geotechnical analyses, physical model tests have been performed at the Coastal & Ocean Basin (COB) in Ostend (BE) within the SOILTWIN project in December 2023. These tests provide insight into the soil behaviour around the monopile foundation based on pore pressure measurements. The comparison of the numerical results with experimental data is essential for an improved calibration of the numerical model as well as for a better understanding of the soil response under various wave loading conditions. The experimental setup and the initial findings will be discussed during the conference.

Stochastic analysis of the seismic performance of sheet-pile walls supporting heterogeneous liquefiable ground

This paper presents the findings of a stochastic analysis conducted to evaluate the effects of the spatial soil variability and temporal base motion variability on the seismic response of sheet pile walls supporting liquefiable ground. A random finite element analysis is used to model the centrifuge experiment performed during the 2020 phase of the Liquefaction Experiments and Analysis Project (LEAP). The relative density of the main soil deposit in the LEAP-2020 experiments ranged from 50 to 75 %. In this study the spatial variability of each centrifuge model was determined from the inflight cone penetration tests (CPT) performed before the shaking phase in each experiment. The vertical spatial correlation length was directly determined from the CPT measurements, while the mean and coefficient of variation of the soil relative density were obtained from an empirical relationship correlating the cone tip resistance with the relative density. The variability observed in the achieved base motions for the LEAP-2020 centrifuge experiments was modeled by generating synthetic base motion time histories that matched the variability observed in the response spectra of the achieved base motions. The results obtained from this analysis shed light on the effects of the base motion variability and soil density on the excess pore pressure development and the triggering of soil liquefaction. It also provides an understanding on the sensitivity of the wall lateral displacement to such variabilities. The observed variability in the simulation response is compared with the variability observed in the centrifuge experiments. This comparison is used to evaluate the capability of the current numerical simulation tools in modeling the variability in the wall response.

Fluid effects in model granular flows

Pore fluid plays a crucial role in many granular flows, especially those in geophysical settings. However, the transition in behaviour between dry flows and fully saturated flows and the underlying physics that relate to this are poorly understood. In this paper, we report the results of small-scale flume experiments using monodisperse granular particles with varying water content and volume in which the basal pore pressure, total pressure, flow height and velocity profile were measured at a section. We compare the results with theoretical profiles for granular flow and with flow regimes based on dimensional analysis. The runout and the centre of mass were also calculated from the deposit surface profiles. As the initial water content by mass was increased from zero to around 10%, we first observed a drop in mobility by approximately 50%, as surface tension caused cohesive behaviour due to matric suction. As the water content was further increased up to 45%, the mobility also increased dramatically, with increased flow velocity up to 50%, increased runout distance up to 240% and reduced travel angle by up to 10° compared to the dry case. These effects can be directly related to the basal pore pressure, with both negative pressures and positive pore pressures being measured relative to atmospheric during the unsteady flow. We find that the initial flow volume plays a role in the development of relative pore pressure, such that, at a fixed relative water content, larger flows exhibit greater positive pore pressures, greater velocities and greater relative runout distances. This aligns with many other granular experiments and field observations. Our findings suggest that the fundamental role of the pore fluid is to reduce frictional contact forces between grains thus increasing flow velocity and bulk mobility. While this can occur by the development of excess pore pressure, it can also occur where the positive pore pressure is not in excess of hydrostatic, as shown here, since buoyancy and lubrication alone will reduce frictional forces.Graphical abstract

Physical modelling of subsea structure removal via thermally induced pore pressure

Given the requirements of environmental agencies, the decommissioning of subsea structures, including removal, recycling, restructuring or re-adaptation, is inherent to deepwater oil exploration activities. During removal by lifting, the possible development of negative excess pore pressure (i.e., suction) in the soil beneath the structure significantly hampers the recovery of subsea infrastructure. This study presents a method for potentially reducing this suction by heating the fine soil under the foundation during uplift to produce excess thermally-induced positive pore pressure, which may reduce or eliminate suction. Reduced-scale physical modelling studies of mudmat/skirted like foundation structures showed that the negative pore pressure generated during pullout was in certain amount counterbalanced by the thermally-induced positive pore pressure, facilitating their removal. However, this process depends on the magnitude of the sustained loading with respect to the pullout resistance. Upward displacements take place just after heating, without further load increments and the displacement rate is dependent on the level of sustained load.

Modeling heterogeneity and permeability evolution in a compaction band using a phase-field approach

Compaction bands are tabular zones of localized compressive deformation associated with porosity and permeability reduction. Depending on their orientation, compaction bands can act as barriers to fluid flow, and can be detrimental to fluid production in oil and gas reservoirs, as well as in CO2 sequestration. The process of permeability reduction and the development of excess pore pressures during compaction band formation in a heterogeneous rock mass are not fully understood. Furthermore, few studies have modeled compaction band formation considering coupled hydromechanical processes. In this study, we propose a coupled hydromechanical, phase-field approach for capturing the formation and propagation of compaction bands in heterogeneous porous media. Breakage mechanics is adopted to characterize the free energy function in the intact and damaged material. The resulting phase-field variable provides a measure of the degree of grain crushing. Permeability reduction in the zone of compaction localization is modeled using the Kozeny–Carman equation accounting for microstructural evolution. Numerical simulations demonstrate the ability of the model to capture compaction band formation, porosity reduction, and permeability evolution under drained and undrained conditions. The results highlight the role of effective confining pressure, drainage conditions, and material parameters on the styles of compaction bands that form.

Modelling seasonal landslide motion: Does it only depend on fluctuations in normal effective stress?

AbstractLandslide motion is often simulated with interface‐like laws able to capture changes in frictional strength caused by the growth of the pore water pressure and the consequent reduction of the effective stress normal to the plane of sliding. Here it is argued that, although often neglected, the evolution of all the 3D stress components within the basal shear zone of landslides also contributes to changes in frictional strength and must be accounted for to predict changes in seasonal velocity. For this purpose, an augmented sliding‐consolidation model is proposed which allows for the computation of excess pore pressure development and downslope sliding with any constitutive law with 3D stress evolution. Simulations of idealised infinite slope models subjected to hydrologic forcing are used to study the role of in‐situ stress conditions and stress rate multiaxiality. Specifically, a Drucker‐Prager perfectly plastic model is used to replicate frictional failure and shear deformation at the base of landslides. The model reveals that conditions amenable to the shearing of a frictional interface are met only after numerous rainfall cycles, that is, when multiaxial stress rates are suppressed. In this case, the landslide is predicted to move through a seasonal ratcheting controlled only by the effective stress component normal to the plane of sliding. By contrast, in newly formed landslides, the multiaxial stress evolution is found to produce further regimes of motion, from plastic shakedown to cyclic failure, neither of which can be captured by interface‐like frictional laws. Notably, the model suggests that a transition across these regimes can emerge in response to an aggravation of the magnitude of forcing, implying that (i) fluctuations in climate may alter the seasonal trends of motion observed today; (ii) our ability to quantify landslide‐induced risks is impaired unless proper geomechanical models are used to examine their long‐term dynamics.

Numerical simulation of pore pressure development beneath suction anchor under undrained condition during uplift

Suction anchors are increasingly used in the deep-sea anchoring infrastructure, and the accurate evaluation of the uplift behavior is of great importance. In this paper, the development of pore pressure under the suction anchor lid during the undrained pullout process was investigated. The numerical model based on the solid-fluid two-phase theory and by updated Lagrangian method was employed. Laboratory experimental studies with three different aspect ratios were carried out to verify the effectiveness and accuracy of numerical method. Parametric studies were performed to analyze the effect of key soil properties on the accumulation of pore pressure. The results show that the numerical method in terms of capturing the failure pattern and simulating pore pressure evolution was significantly efficient. The development of excess pore pressure increased faster at the beginning, followed by a slower increase and then tended to level off. The maximum excess pore pressure appeared just beneath the suction anchor, and decreased with the increasing soil depth. Parametric analysis showed that a greater negative pore pressure needed a larger shear modulus, MCC swelling index and critical stress ratio, or a smaller porosity, permeability and MCC compression index. The softening coefficient had little effect on the pore pressure generation.

Variation in axial load distribution of piles in liquefiable slope by centrifuge test

During seismic loading, the axial pile response is significantly affected, due to the development and dissipation of excess pore pressure in the foundation soil. Positive skin friction is reduced by liquefaction during shaking, while negative skin friction is caused by reconsolidation after the end of the shaking. In particular, the negative skin friction causes a large drag load, resulting in pile compression and settlement. Thus, this study investigated the axial pile response under seismic loading by performing dynamic centrifuge tests on liquefied sand with the slopes of 15° and 27°. The centrifuge tests were conducted to simulate the end-bearing single and 2 × 2 group piles embedded in a liquefied slope. During the centrifuge spinning, the obtained drag load due to the negative skin friction agreed with the value calculated with the beta method for shaft resistance. The model was shaken with sinusoidal base motion with peak amplitudes of 0.1 and 0.2 g, and it was found that the liquefaction caused a loss of spin-induced drag load. The dissipation of excess pore pressure resulted in reconsolidation settlement, leading to large drag loads acting on the pile. In addition, the negative skin friction was found to be larger than the liquefied residual strength recommended in the current practice codes. As a result, it is seen that the consideration of the strength reduction is not necessary for the calculation of negative skin friction induced by reconsolidation, since the liquefied ground regains skin friction after the dissipation.

Insights from in-situ pore pressure monitoring around a wind turbine monopile

Offshore wind turbine monopiles transfer large bending moments into the soil. Due to the large diameter of the foundation, a significant volume of soil is mobilised to absorb the incoming forces from wind, waves and current. If excess pore pressures are generated during load cycles, the long drainage path length could lead to accumulation of pore pressures and loss of strength and stiffness in the soil. During monopile design, undrained behaviour of the soil during cyclic loading is often assumed to conservatively assess the lateral pile resistance during storms. In-situ monitoring of pore pressures in the soil surrounding a monopile foundation allows the assumption of undrained behaviour to be checked for several episodes during the operational life of the foundation. A 3-year period of monitoring data was analysed for a North Sea monopile which included normal operation, storm episodes and rotor stops. Pore pressure and earth pressure sensors were positioned on the pile at two different depths in the most intensely loaded layer of surface sand. The data provides insight into the amount of excess pore pressure accumulation with storm loading not leading to significant accumulation of excess pore pressure. The only episodes showing a significant development of excess pore pressures were rotor stops with suction pressures being developed on the unloaded side of the monopile. The observed response suggests that significant pile-soil stiffness reductions in the top layer due to cyclic loading are not expected for the monitored location. The study shows the value of monitoring pore pressures in the soil mass surrounding offshore foundations for checking design assumptions.

Effective stress analysis of dynamically installed anchor installation in rate-dependent clay

Accurate prediction of the final embedment depth and installation-induced excess pore pressure of dynamically installed anchors (DIAs) are prerequisites for the design of subsequent anchor set-up during a waiting period and capacity under operational loadings. This paper presents a dynamic large deformation finite element approach within the effective stress framework and investigates the installation behaviour of DIAs in clay. An elasto-viscoplastic model is adopted to capture strain rate dependency of the undrained shear strength. Development of excess pore pressure for various impact velocities, anchor geometries, viscosity parameters are quantified. An expression is proposed for estimating the final embedment depth of DIAs in clay.

Natural rubber latex treatment of sand: A novel remediation technique for soil liquefaction

We present a novel method to improve liquefaction resistance in sandy soil by treating it with natural rubber latex (NRL) using a pressurised permeation technique. Sand samples treated with NRL solutions of 10%, 15%, 20%, and 30% concentrations were tested. The liquefaction resistance of the untreated and NRL-treated sand samples was evaluated by performing strain-controlled cyclic triaxial tests with 0.15%, 0.3%, and 0.5% axial strain amplitudes. We find that NRL treatment reduces excess pore pressure development in the sand, thereby increasing liquefaction resistance. Further, this study also explores the effect of NRL treatment on stiffness deterioration in the sand under cyclic loading. NRL-treated sand samples showed resistance to progressive stiffness degradation even after many cycles due to interparticle bonding developed by virtue of the solidified NRL fibres. From a pilot study on the behaviour of treated sand under stress-controlled loading, it was seen that irrespective of the pore pressure development, NRL effectively resist the strain development in sand. Hence, we propose the utilisation of NRL as an environment-friendly soil stabiliser for mitigating liquefaction.

Centrifuge modeling studies on effects of composition on liquefaction of mine waste rock

Liquefaction behavior is investigated in four centrifuge modeling tests performed on loose mine waste rock composed of different ratios of gravel, sand, and fines. These specimens have similar shear wave velocity of about 117–128 m/s. Each centrifuge model was subjected to a 0.27-g sinusoidal base shaking amplitude. It was found that soil liquefaction was observed for all models with the exception of a model with a mixture of 70/30% gravel/sand. Dilative behavior was observed from the acceleration time histories for all models and was dependent on the excess pore pressure development and effective stress. Pore pressure dissipation during shaking was found in the models without fines. Degradation and distortion of shear stress-strain loops were also observed. Larger settlement and volumetric strain were observed after the tests for the models with fines. Based on the results from this study, an increase in sand and fines content causes an increase in liquefaction potential of mine waste rock.

Experimental investigation of transient bending moment of piles during seismic liquefaction

The purpose of this paper is to study the behavior of pile-supported structures in liquefiable soils, specifically when the soil surrounding the pile transits from no-liquefaction to full-liquefaction. A series of shaking table tests were performed on four pile-supported structures subjected to different input motions and, as a result, with different times to reach full-liquefaction. The bending moment of the piles during the transient phase is compared with those predicted in pre- and post-liquefaction stages. The experimental results showed that the maximum bending moment may occur during the transient phase (i.e. during the development of excess pore pressure before the soil is fully liquefied). Arguably, the observed amplification in bending moment is caused by the the tuning effect between the predominant frequency of the input motion and the frequency of the pile-supported structure, which is progressively decreasing during the liquefaction process. Results are presented using a non-dimensional framework whose parameters are derived from the governing mechanics. A new parameter TAF (Transient Amplification Factor) is defined to predict the design bending moment during the transient phase. It is shown that the transient bending moment can be obtained from the newly introduced parameter TAF, and the values of maximum bending moments in the pre- and post-liquefaction stages. It is found that TAF is a function of two easily obtainable parameters: (a) time taken to reach full liquefaction, this can be obtained through site response analysis; (b) elongation of natural period of vibration, expressed as ratio of the time period of the structure at full liquefaction to the time period at zero-liquefaction. Finally, practical implications of the main findings are discussed.

Numerical Assessment of the Loading Factors Affecting Liquefaction-Induced Failure

This paper presents parametric studies that assess the role of loading factors (i.e., number of cycles, frequency, and amplitude) on liquefaction-induced failure by performing numerical simulations. Most of the existing literature considers the effects of the soil properties on the development of excess pore pressure with few research endeavours focusing on the effects of the input motion itself. Numerical simulations are performed herein, via the advanced software platform OpenSees, to generate several finite element models that consider non-linear development of pore pressure inside the soil. Several sinusoidal inputs were considered to study the effects of the various loading factors and compare the responses. The main findings arise from evaluating the effects of several input motion parameters (number of cycles, frequency, and amplitude) on soil liquefaction through numerical simulations. This research study, based on state-of-the-art knowledge, may be applied to assess future seismic events and to update or propose new code provisions for soil liquefaction.

Strength Characteristics of Clay–Rubber Waste Mixtures in Low-Frequency Cyclic Triaxial Tests

This paper presents the results of consolidated and undrained (CU) triaxial cyclic tests related to the influence of tire waste addition on the strength characteristics of two different soils from Southern Poland: unswelling kaolin and swelling red clay. The test procedure included the normally consolidated remolded specimens prepared from pure red clay (RC) and kaolin (K) and their mixtures with two different fractions of shredded rubber powder (P) and granulate (G) in 5%, 10%, and 25% mass proportions. All samples were subjected to low-frequency cyclic loading carried out with a constant stress amplitude. Analysis of the results includes consideration of the effect of rubber additive and number of load cycles on the development of excess pore pressure and axial strain during the cyclic load operation and on the maximum stress deviator value. A general decrease in the shear strength due to the cyclic load operation was observed, and various effects of shear strength depended on the mixture content and size of the rubber waste particles. In general, the use of soil–rubber mixtures, especially for expansive soils and powder, should be treated with caution for cyclic loading.

Centrifuge modelling on monotonic and cyclic lateral behaviour of monopiles in kaolin clay

This paper details a programme of centrifuge tests (100g) on the monotonic and cyclic lateral behaviour of large-diameter semi-rigid monopiles in kaolin clay. Piles of different diameters were used, with a maximum prototype value of 5·9 m. The pile–soil interaction was investigated in terms of the overall load–displacement behaviour, distribution of bending moment and pile deflection, p–y curves and excess pore pressure response. For monotonic loading, procedures were developed to fit the monotonic p–y curves obtained from different pile and soil properties using a hyperbolic tangent model. In addition, an explicit model was derived to assess the development of excess pore pressure within the surficial soil, which was found to approach the vertical effective soil stress at large displacement. For cyclic loading, the accumulation of lateral pile displacement could be accompanied by either increase or decrease in soil resistance, which strongly depends on the loading rate. Strong correlation between the development of peak soil resistance and a dimensionless velocity group was obtained, and was supported by the measured excess pore pressure. Under undrained conditions, the peak soil resistance decreases almost linearly with the accumulation of excess pore pressure.

Groundwater erosion of coastal gullies along the Canterbury coast (New Zealand): a rapid and episodic process controlled by rainfall intensity and substrate variability

Abstract. Gully formation has been associated to groundwater seepage in unconsolidated sand- to gravel-sized sediments. Our understanding of gully evolution by groundwater seepage mostly relies on experiments and numerical simulations, and these rarely take into consideration contrasts in lithology and permeability. In addition, process-based observations and detailed instrumental analyses are rare. As a result, we have a poor understanding of the temporal scale of gully formation by groundwater seepage and the influence of geological heterogeneity on their formation. This is particularly the case for coastal gullies, where the role of groundwater in their formation and evolution has rarely been assessed. We address these knowledge gaps along the Canterbury coast of the South Island (New Zealand) by integrating field observations, luminescence dating, multi-temporal unoccupied aerial vehicle and satellite data, time domain electromagnetic data and slope stability modelling. We show that gully formation is a key process shaping the sandy gravel cliffs of the Canterbury coastline. It is an episodic process associated to groundwater flow that occurs once every 227 d on average, when rainfall intensities exceed 40 mm d−1. The majority of the gullies in a study area southeast (SE) of Ashburton have undergone erosion, predominantly by elongation, during the last 11 years, with the most recent episode occurring 3 years ago. Gullies longer than 200 m are relict features formed by higher groundwater flow and surface erosion > 2 ka ago. Gullies can form at rates of up to 30 m d−1 via two processes, namely the formation of alcoves and tunnels by groundwater seepage, followed by retrogressive slope failure due to undermining and a decrease in shear strength driven by excess pore pressure development. The location of gullies is determined by the occurrence of hydraulically conductive zones, such as relict braided river channels and possibly tunnels, and of sand lenses exposed across sandy gravel cliffs. We also show that the gully planform shape is generally geometrically similar at consecutive stages of evolution. These outcomes will facilitate the reconstruction and prediction of a prevalent erosive process and overlooked geohazard along the Canterbury coastline.

SANISAND-MSf: a sand plasticity model with memory surface and semifluidised state

A new constitutive model for sand is formulated by incorporating two new constitutive ingredients into the platform of a reference critical state compatible bounding surface plasticity model with kinematic hardening, in order to address primarily the undrained cyclic response. The first ingredient is a memory surface for more precisely controlling stiffness affecting the plastic deviatoric and volumetric strains and ensuing excess pore pressure development in the pre-liquefaction stage. The second ingredient is the concept of a semifluidised state and the related formulation of stiffness and dilatancy degradation, aiming at modelling large shear strain development in the post-liquefaction stage. In parallel, a modified flow rule aimed at providing a better description of non-proportional monotonic and cyclic loading is introduced. With a single set of constants, for which a detailed calibration procedure is provided, this new model successfully simulates undrained cyclic torsional and triaxial tests with different cyclic stress ratios, separately for the pre- and post-liquefaction stages, as well as liquefaction strength curves based on [Formula: see text] and shear strain criteria for initial liquefaction. The successful reproduction of the sand element response under undrained cyclic shearing contributes to future applications in realistic and thorough seismic site response analysis.