Non‐linear consolidation analysis of soil with variable compressibility and permeability under cyclic loadings

In this paper, a simple semi‐analytical method has been developed to solve the one‐dimensional non‐linear consolidation problems by considering the changes of compressibility and permeability of the soil layer, subjected to complicated time‐dependent cyclic loadings at the ground surface. The solution presented here takes into account e ～ lg k_v and e ～ lg σ′ linear responses. With c_k the slope of the e ～ lg k_v line and c_c as the slope of the e ～ lg σ′ line, the identified parameter c_c/c_k is found to control the rate of consolidation. Using the solutions obtained, some diagrams are prepared and the relevant behaviours of one‐dimensional non‐linear consolidation of saturated soft soil under cyclic loadings are discussed. The method in this paper does not require any special data; conventional oedometer data can be used. Therefore, the method is particularly efficient and convenient for engineering practice. Copyright © 2006 John Wiley & Sons, Ltd.     

Wave propagation in nonlinear one‐dimensional soil model

The objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of shear modulus as a function of strain. In this paper the meaning of the material stiffness obtained from impact and harmonic excitation tests on a surface slab is discussed. A one‐dimensional discrete model with the nonlinear material stiffness is used for this purpose. When a static load is applied followed by an impact excitation, if the amplitude of the impact is very small, the measured wave velocity using the cross‐correlation indicates the wave velocity calculated from the tangent modulus corresponding to the state of stress caused by the applied static load. The duration of the impact affects the magnitude of the displacement and the particle velocity but has very little effect on the estimation of the wave velocity for the magnitudes considered herein. When a harmonic excitation is applied, the cross‐correlation of the time histories at different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress–strain loop under steady‐state condition. Copyright © 2008 John Wiley & Sons, Ltd.     

Rigid cylinder in a transversely isotropic half‐space under lateral loads

A boundary integral equation method is presented for a rigid cylindrical pipe‐pile of finite length embedded in a transversely isotropic half‐space under lateral loads. In the framework of three‐dimensional elastostatics, the complicated soil‐structure interaction problem is shown to be reducible to three coupled Fredholm integral equations. Through an analysis of the associated Cauchy singular kernels, the intrinsic singular characteristics of the radial, angular, and vertical interfacial load transfers are rendered explicit. By means of a complicated numerical procedure, detailed results on the three‐dimensional load–transfer process are provided for benchmark comparison and practical applications. Copyright © 2011 John Wiley & Sons, Ltd.     

A modulus‐multiplier approach for non‐linear analysis of laterally loaded pile groups

A modulus‐multiplier approach, which applies a reduction factor to the modulus of single pile p–y curves to account for the group effect, is presented for analysing the response of each individual pile in a laterally loaded pile group with any geometric arrangement based on non‐linear pile–soil–pile interaction. The pile–soil–pile interaction is conducted using a 3D non‐linear finite element approach. The interaction effect between piles under various loading directions is investigated in this paper. Group effects can be neglected at a pile spacing of 9 times the pile diameter for piles along the direction of the lateral load and at a pile spacing of 6 times the pile diameter for piles normal to the direction of loading. The modulus multipliers for a pair of piles are developed as a function of pile spacing for departure angle of 0, 90, and 180sup>/sup> with respect to the loading direction. The procedure proposed for computing the response of any individual pile within a pile group is verified using two well‐documented full‐scale pile load tests. Copyright © 2006 John Wiley & Sons, Ltd.     

Analysis of the compression of structured soils using the disturbed state concept

The aim of this note is to quantify the influence of soil structure on the compression behaviour of natural soils using the disturbed state concept (DSC). The behaviour of the fully adjusted state is chosen to be that of the corresponding soil in a reconstituted condition so that the disturbance function is a direct measure of the effects of soil structure. A new DSC compression model is proposed. This model is able to describe the compression behaviour of structured soils under loading, swelling and reloading. Special versions of the proposed model are also described for situations (a) where the compression behaviour of the corresponding reconstituted soils is linear in the e–ln p′ space and (b) where the compression is one‐dimensional. The ability of the proposed model and its various versions to describe the compression behaviour of structured soils has been verified. Copyright © 2000 John Wiley & sons, Ltd.     

Indentation of a free‐falling lance penetrometer into a poroelastic seabed

A solution is developed for the build‐up, steady and post‐arrest dissipative pore fluid pressure fields that develop around a blunt penetrometer that self‐embeds from freefall into the seabed. Arrest from freefall considers deceleration under undrained conditions in a purely cohesive soil, with constant shear strength with depth. The resulting decelerating velocity field is controlled by soil strength, geometric bearing capacity factors, and inertial components. At low impact velocities the embedment process is controlled by soil strength, and at high velocities by inertia. With the deceleration defined, a solution is evaluated for a point normal dislocation penetrating in a poroelastic medium with a prescribed decelerating velocity. Dynamic steady pressures, P_D, develop relative to the penetrating tip geometry with their distribution conditioned by the non‐dimensional penetration rate, U_D, incorporating impacting penetration rate, consolidation coefficient and penetrometer radius, and the non‐dimensional strength, N_D, additionally incorporating undrained shear strength of the sediment. Pore pressures develop to a steady peak magnitude at the penetrometer tip, and drop as P_D=1/x_D with distance x_D behind the tip and along the shaft. Peak induced pressure magnitudes may be correlated with sediment permeabilities, post‐arrest dissipation rates may be correlated with consolidation coefficients, and depths of penetration may be correlated with shear strengths. Together, these records enable strength and transport parameters to be recovered from lance penetrometer data. Penetrometer data recorded off La Palma in the Canary Islands (J. Volcanol. Geotherm. Res. 2000; 101 :253) are used to recover permeabilities and consolidation coefficients from peak pressure and dissipation response, respectively. Copyright © 2004 John Wiley & Sons, Ltd.     

Analytical theory for one‐dimensional consolidation of clayey soils exhibiting rheological characteristics under time‐dependent loading

This paper presents analytical solutions to the one‐dimensional consolidation problem taking into consideration the rheological properties of clayey soil under variable loadings. A four‐element rheological model is introduced, and different loading types are involved, i.e. constant loading, one‐step loading, triangular loading, rectangular loading, and isosceles–trapezoidal cyclic loading. The differential equations governing consolidation are solved by the Laplace transform. Based on the solutions obtained, the influences of the rheological parameters and loading conditions on the consolidation process are investigated. It has been shown that the consolidation behavior is mainly governed by four dimensionless parameters, a₁, a₂, b, and T_v0. Load shape has a great influence on the rate of consolidation. A decrease either in the modulus of the spring in the Kelvin body or in the viscosity coefficient of independent dashpot will slow down the rate of consolidation. An increase in the viscosity coefficient of the dashpot in the Kelvin body will make the rate of consolidation increase at an early stage but decrease at a later stage. For isosceles–trapezoidal cyclic loading, the consolidation rate in each cycle reaches a maximum at the end of the constant loading phase and the minimum at the end of this cycle. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical analysis of a one‐dimensional infiltration problem in unsaturated soil by a seepage–deformation coupled method

A multiphase coupled elasto‐viscoplastic finite element analysis formulation, based on the theory of porous media, is used to describe the rainfall infiltration process into a one‐dimensional soil column. Using this framework, we have numerically analyzed the generation of pore water pressure and deformations when rainfall is applied to the soil. A parametric study, including rainfall intensity, soil–water characteristic curves, and permeability, is carried out to observe their influence on the changes in pore water pressure and volumetric strain. From the numerical results, it is shown that the generation of pore water pressure and volumetric strain is mainly controlled by material parameters α and n′ that describe the soil–water characteristic curve. A comparison with the laboratory results shows that the proposed method can describe very well the characteristics observed during the experiments of one‐dimensional water infiltration into a layered unsaturated soil column. Copyright © 2010 John Wiley & Sons, Ltd.     

Inverse analysis for geomaterial parameter identification using Pareto multiobjective optimization

An inverse analysis method that combines the back propagation neural network (BPNN) and vector evaluated genetic algorithm (VEGA) was proposed to identify mechanical geomaterial parameters for a more accurate prediction of deformation. The BPNN is used to replace the time‐consuming numerical calculations, thus enhancing the efficiency of the inverse analysis. The VEGA is used to find the Pareto‐optimal solutions to multiobjective functions. Unlike traditional back‐analysis methods which are based on only 1 type of field measurement and a single objective function, this proposed method can consider multiple field observations simultaneously. The proposed method was applied to the Shapingba foundation pit excavation located in Chongqing city, China. Two types of measurements are considered in the method simultaneously: the displacements in the x‐direction (north orientation) and those in the y‐direction (east orientation). Five deformation modulus parameters for artificial backfill soil, silty clay, siltstone, sandstone, and mudstone were selected as the inversion parameters. Compared with the weighted sum approach, the proposed method was demonstrated as an efficient multi‐objective optimization tool for back calculating undetermined parameters. After performing a forward‐calculation using the optimized parameters obtained by the inverse analysis, the predicted results were well consistent with the practical deformation in magnitude and trend.     

橡胶砂弹性动力学参数的弯曲-伸缩元试验研究

橡胶砂作为一种廉价环保的土工材料,应用前景十分广阔,关于其力学特性的研究具有重要意义。采用弯曲–伸缩元法对橡胶砂进行了剪切波和压缩波联合测试,分别采用离散频率扫描法和初达波法确定了剪切波和压缩波的传播时间,进而得到其弹性动力学参数：初始剪切模量、侧限模量和泊松比,分析了橡胶含量和围压对橡胶砂弹性动力学参数的影响。结果表明：在同样围压下,随着橡胶含量的增加,橡胶砂初始剪切模量和侧限模量逐渐减小,初始泊松比逐渐增大;在同样橡胶含量下,随着围压的增大,橡胶砂初始剪切模量和侧限模量逐渐增大,初始泊松比逐渐减小。最后,在此基础上进行了两种因素耦合效应的分析以及相关力学机制的探讨。     

Three‐dimensional shakedown solutions for anisotropic cohesive‐frictional materials under moving surface loads

Previous work on three‐dimensional shakedown analysis of cohesive‐frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three‐dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower‐bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross‐anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only. Copyright © 2013 John Wiley & Sons, Ltd.     

变荷载下不排水桩复合地基桩间土和下卧层土固结分析

基于悬浮不排水桩复合地基超静孔隙水压力解答,建立了变荷载条件下复合地基桩间土和下卧层土平均固结度的解析解。通过与有限元解答的对比,证明了解析解的正确性和有效性。利用解析解对桩间土和下卧层土的固结性状进行了分析。研究发现,桩间土的固结速率远大于下卧层土,并随置换率和桩土模量比的增加而增大。在工程设计采用的置换率范围内,桩的贯入比较小时,下卧层土的固结速率受置换率和桩土模量比的影响很小;桩的贯入比较大时,下卧层土的固结速率随置换率和桩土模量比的增加而增大。桩间土和下卧层土的固结存在临界置换率和临界桩土模量比,超过临界值后,桩间土和下卧层土的固结速率不随桩土模量比和置换率的增加而增大。桩间土和下卧层土的固结速率随加荷速率的增加而增大。     

Modulus of subgrade reaction that varies with magnitude of displacement of cohesionless soil

Modulus of subgrade reaction is one of the required design parameters in any structural analyses of shallow foundations. However, the constant values of modulus of subgrade reaction that are determined from either from literature studies or the results of plate load tests, regardless of magnitude of soil’s displacement under design loads, have been used in structural designs of foundations. In this study, the results of 43 published full-scale field plate load tests in cohesionles soils were gathered to expose any variation in the values of modulus of subgrade reaction as soil’s displacement increases. Extensive finite element (FE) analyses were carried out while the results of FE analyses were compared with actual measured results of field load tests. The results of analyses indicated that both the modulus of subgrade reaction decrease with increase in magnitude of displacement of soils and internal forces of a design of structural frame is higher with the values of modulus of subgrade reaction that are sensitive to soil’s displacement are implemented into analyses. Therefore, structural dimensions of any structure with the constant values of modulus of subgrade reaction would not be a precise engineering solution.     

Consolidation in spatially random unsaturated soils based on coupled flow‐deformation simulation

This paper integrates random field simulation of soil spatial variability with numerical modeling of coupled flow and deformation to investigate consolidation in spatially random unsaturated soil. The spatial variability of soil properties is simulated using the covariance matrix decomposition method. The random soil properties are imported into an interactive multiphysics software COMSOL to solve the governing partial differential equations. The effects of the spatial variability of Young's modulus and saturated permeability together with unsaturated hydraulic parameters on the dissipation of excess pore water pressure and settlement are investigated using an example of consolidation in a saturated‐unsaturated soil column because of loading. It is found that the surface settlement and the pore water pressure profile during the process of consolidation are significantly affected by the spatially varying Young's modulus. The mean value of the settlement of the spatially random soil is more than 100% greater than that of the deterministic case, and the surface settlement is subject to large uncertainty, which implies that consolidation settlement is difficult to predict accurately based on the conventional deterministic approach. The uncertainty of the settlement increases with the scale of fluctuation because of the averaging effect of spatial variability. The effects of spatial variability of saturated permeability k_sat and air entry parameters are much less significant than that of elastic modulus. The spatial variability of air entry value parameters affects the uncertainties of settlement and excess pore pressure mostly in the unsaturated zone. Copyright © 2016 John Wiley & Sons, Ltd.     

Three‐dimensional analysis of the seismic behaviour of micropiles used in the reinforcement of saturated soil

This paper includes a numerical study of the behaviour of micropiles used for the reinforcement of saturated soil. Analysis is carried out using the (u–p) formulation (displacement for the solid phase and pore‐pressure for the fluid phase) implemented in a three‐dimensional finite element program. The soil behaviour is described by means of a cyclic elastoplastic constitutive relation which was developed within the framework of the bounding surface concept. The paper is composed of three parts. The first one is concerned with a presentation of the numerical model; the second includes analysis of the seismic behaviour of a single micropile; the last part deals with the group effect under seismic loading. Copyright © 2001 John Wiley & Sons, Ltd.     

A new approach for vertical impedance in radially inhomogeneous soil layer

In the past studies on pile vibrations, the soil around the pile is mainly regarded as homogeneous or multi‐layered piecewise homogeneous. However, under most engineering conditions, the surrounding soil becomes seriously disturbed due to construction effects. This may strengthen or weaken the shear modulus of the soil resulting in the soil becoming radially inhomogeneous. As a consequence of this, El Naggar extended Novak's plane‐strain model to account for the radial inhomogeneity by the use of multiple springs connected in series. Rather than using this approach, this paper proposes a new model which is thought to be theoretically more rigorous and one which may be described as complex stiffness transfer model. It is shown that the solution developed in this study agrees well with the more limited solutions of Novak and Dotson and Veletsos under several special conditions. Finally, the scope of application has been enlarged as a result of the generalizations made in the present model. Copyright © 2011 John Wiley & Sons, Ltd.     

A coupled analysis of cavity and pore volume changes for pulse tests conducted in soft clay deposits

This paper presents a new interpretation method for pulse tests, a field permeability test that allows for rapid measurements of hydraulic conductivity k in aquitards. This new method applies to soft clay deposits. To initiate a pulse test, a known volume of water is injected into the sand filter of a monitoring well isolated by a packer. The resulting pressure increase yields an outward movement of the sand filter cavity wall. After presenting the usual interpretation methods and their limits, this paper proposes a new interpretation method based on a coupled analysis of the pressure and displacement fields in the soil using the Biot–Darcy formulation. A series of analytical and numerical non‐dimensional velocity graphs, normalized plots of the mean hydraulic head difference versus its rate of change, are given. For a linear elastic material, these type curves show relatively small variations with the sand filter aspect ratio and the Poisson ratio of the tested clay. The type curves are also found to be independent of the clay compressibility (m_v) and k, an important result. A series of pulse tests conducted in a soft marine clay deposit near Montreal, Canada, are interpreted with the proposed method. The hydraulic conductivity values calculated from these tests are closely correlated with independent estimates obtained using long‐duration variable‐head tests. Compared with previous interpretation methods, the proposed method allows soil volume changes to be reconciled with cavity expansion phenomena and the range of type curves available for the interpretation of test data to be constrained. Copyright © 2013 John Wiley & Sons, Ltd.     

Physics‐based simulation of multiple interacting crack growth in brittle rocks driven by thermal cooling

Crack growth in hot brittle rocks, driven by thermal cooling, was simulated using a coupled two‐dimensional discrete element and heat transport model that explicitly includes the random initiation and subsequent propagation of interacting cracks. The model clearly predicts that a quasi‐hierarchical array of subparallel cracks, oriented along the direction of the temperature gradient, is formed under small to moderately large thermally generated strain load conditions. The simulation results also demonstrate that, after an initial transient, thermal cracks propagate in a stable fashion with a velocity that scales with ~ t^{− 1/2}. However, under large thermal strain loads, a more complicated geometry composed of cracks that curve and coalesce develops during the later stages of crack growth. Copyright © 2016 John Wiley & Sons, Ltd.     

小应变条件下隧道开挖引起的横向沉降槽分析

利用配备有局部高精度位移传感器STDTTS+UNSAT型号的高级应力路径三轴测试系统,对小应变条件下的土体进行了试验。采用三维数值分析方法和应用小应变模型及摩尔-库仑模型,分别对隧道施工引起横向沉降槽的影响展开了研究。结果表明,在小应变条件下土体表现出高模量以及模量随应变增加而衰减的特性,且当剪切应变应用对数坐标（lg）表示时,剪切模量随剪切应变增加而衰减的规律可以表示为带有2个拐点的反S型曲线,同时,利用式（1）对其进行很好地模拟;两类模型计算结果的规律性是相似的,采用小应变模型所得到的最大下沉位移符合工程实际,且横向沉降槽的宽度也能符合现场实测,即能够适当地解决采用一般模型所得到的最大下沉位移合理,而横向沉降槽的宽度过宽,或者横向沉降槽宽度合理,而最大下沉位移过小,两者间存在对立的问题。通过计算结果对比进一步验证了考虑小应变条件下土体特性的必要性以及二次开发所得到的模型的合理性与可适用性。