Analysis of undrained cavity expansion in elasto‐plastic soils with non‐linear elasticity

A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd.     

Upper bound limit analysis using linear finite elements and non‐linear programming

A new method for computing rigorous upper bounds on the limit loads for one‐, two‐ and three‐dimensional continua is described. The formulation is based on linear finite elements, permits kinematically admissible velocity discontinuities at all interelement boundaries, and furnishes a kinematically admissible velocity field by solving a non‐linear programming problem. In the latter, the objective function corresponds to the dissipated power (which is minimized) and the unknowns are subject to linear equality constraints as well as linear and non‐linear inequality constraints. Provided the yield surface is convex, the optimization problem generated by the upper bound method is also convex and can be solved efficiently by applying a two‐stage, quasi‐Newton scheme to the corresponding Kuhn–Tucker optimality conditions. A key advantage of this strategy is that its iteration count is largely independent of the mesh size. Since the formulation permits non‐linear constraints on the unknowns, no linearization of the yield surface is necessary and the modelling of three‐dimensional geometries presents no special difficulties. The utility of the proposed upper bound method is illustrated by applying it to a number of two‐ and three‐dimensional boundary value problems. For a variety of two‐dimensional cases, the new scheme is up to two orders of magnitude faster than an equivalent linear programming scheme which uses yield surface linearization. Copyright © 2001 John Wiley & Sons, Ltd.     

Interpretation of torsional shear results for nonlinear stress–strain relationship

Torsional tests are frequently used in material testing to determine shear stress–strain relationships. In the tests, torque is applied to a solid or hollow cylindrical specimen, whose twist angle is recorded as the response. While the raw test results are readily available in terms of torque–twist sets, converting them into intrinsic stress–strain relationships is not always straightforward, owing to the nonlinear stress–strain behavior and the non‐uniform radial distribution of the stress and strain fields. It follows that the approaches currently used to interpret torsional data are essentially empirical based. The present paper addresses this data interpretation issue from a more rigorous perspective. The paper includes two parts: the first part concerns hollow cylindrical specimens under torsion, where an analytical‐based recursive procedure is developed; and the second part deals with solid torsional specimens, for which an analytical solution given by Nadai is directly employed. Examples pertaining to these two cases are examined and the resulting shear stress–strain curves are compared with those based on a number of popular empirical methods. The proposed analytical approaches can serve as a benchmark for checking the accuracy and validity of those empirical propositions. It is worth noting that, for brittle materials exhibiting softening behavior, only the analytical method gives the correct interpretation. Copyright © 2007 John Wiley & Sons, Ltd.     

Compressibility and permeability of sand–kaolin mixtures. Experiments versus non‐linear homogenization schemes

This study deals with the behaviour of mixtures of sand and saturated kaolin paste considered as composite materials made of permeable and deformable (with non‐linear behaviour) matrix (the kaolin paste) with rigid and impervious inclusions (the sand grains). Oedometric and permeability tests highlight the key role of the state of the clay paste, and show the existence of a threshold of sand grain concentration above which a structuring effect influences both compressibility and permeability. At the light of these experiments two homogenization schemes (with simplifying assumptions to make the problem manageable) are considered to model these two parameters. Qualitative and quantitative comparisons with experimental data point out their respective domain of interest and limitations: a tangent homogenization scheme is shown to be sufficient to describe the macroscopic properties for dilute sand concentration; above the concentration threshold, the structuring effect is captured by the new homogenization scheme developed in this paper. Copyright © 2010 John Wiley & Sons, Ltd.     

Study of bond stress–slip relationship and radial dilation in prestressed concrete

This paper presents two test procedures for evaluating the bond stress–slip and the slip–radial dilation relationships when the prestressing force is transmitted by releasing the steel (wire or strand) in precast prestressed elements. The bond stress–slip relationship is obtained with short length specimens, to guarantee uniform bond stress, for three depths of the wire indentation (shallow, medium and deep). An analytical model for bond stress–slip relationship is proposed and compared with the experimental results. The model is also compared with the experimental results of other researchers. Since numerical models for studying bond‐splitting problems in prestressed concrete require experimental data about dilatancy angle (radial dilation), a test procedure is proposed to evaluate these parameters. The obtained values of the radial dilation are compared with the prior estimated by numerical modelling and good agreement is reached. Copyright © 2012 John Wiley & Sons, Ltd.     

Response evaluation for horizontally loaded fixed‐head pile groups using 3‐D non‐linear analysis

The response of laterally loaded pile foundations may be significantly important in the design of structures for such loads. A static horizontal pile load test is able to provide a load–deflection curve for a single free‐head pile, which significantly differs from that of a free‐ or fixed‐head pile group, depending on the particular group configuration. The aim of this paper is to evaluate the influence of the interaction between the piles of a group fixed in a rigid pile cap on both the lateral load capacity and the stiffness of the group. For this purpose, a parametric three‐dimensional non‐linear numerical analysis was carried out for different arrangements of pile groups. The response of the pile groups is compared to that of the single pile. The influence of the number of piles, the spacing and the deflection level to the group response is discussed. Furthermore, the contribution of the piles constituting the group to the total group resistance is examined. Finally, a relationship is proposed allowing a reasonable prediction of the response of fixed‐head pile groups at least for similar soil profile conditions. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.     

A two‐step Taylor–Galerkin algorithm applied to shock wave propagation in soils

This paper presents a two‐step Taylor–Galerkin algorithm formulated in terms of velocities and stresses which can be applied to solid dynamics problems requiring good resolution of small wavelengths, such as propagation of shocks. The proposed model is both fast, as it uses simple linear elements (triangles in 2D and tetrahedra in 3D), and accurate. It avoids locking and mesh alignment problems, and therefore can be applied to localized failure computations without the limitations exhibited by the classical displacement formulations. Copyright © 2003 John Wiley & Sons, Ltd.     

Desiccation shrinkage of non‐clayey soils: a numerical study

A mesoscale model of desiccation of soil based on the evolution of the pore system idealized as bimodal is numerically examined. A simplified evolution of the model reveals a series of characteristics that qualitatively agree with the observed macroscopic experimental findings. The principal mechanism is deemed to be driven by the surface evaporation and water outflow generating a pore pressure gradient resulting in the shrinkage mainly of the largest pores. The amount of shrinkage is a function of (negative) pore pressure and is controlled by the compressibility of the solid matrix. The numerical model includes also the ensuing partial saturation stage initiated by the air entry simulated as a scenario with a moving phase interface inside the pore. The proposed model can be extended beyond the two‐mode porosity soils, to include the multi‐modal porosity, or its statistical representation.Copyright © 2012 John Wiley & Sons, Ltd.     

Incrementally non‐linear plasticity applied to rock joint modelling

Rock joint constitutive modelling is discussed through two new rock joint constitutive relations and a discrete numerical model. Regarding the constitutive relations, we emphasise the number of ‘tensorial zones’, that is, domains of constitutive incremental linearity; they involve four zones for the first (called ‘quadrilinear’) and an infinite number for the second one (called ‘incrementally nonlinear’). Using these formulations, a large class of loading paths can be considered. Hardening through shearing and relations between the normal and tangential directions of the joint (e.g., dilatancy) can be described. Their predictive abilities are checked. Plastic features are included even if the relations are defined outside the elasto‐plastic formalism. These relations obey, hence, the physical evidence as the plastic limit criterion and flow rule. The flow rule is nonassociated, and the corresponding link with the nonsymmetry of the constitutive matrix is examined. Comparisons between the two relations and the discrete numerical model, that is, a direct numerical simulation, which is fundamentally different, also are discussed within the context of infilled rock joints. Copyright © 2011 John Wiley & Sons, Ltd.     

Inextensible reinforcement on non‐linear elasto‐plastic subgrade under oblique pull

In this paper, a rational analysis of pullout resistance of inextensible sheet reinforcement subjected to oblique end force has been presented considering a non‐linear (hyperbolic), elasto‐plastic, normal stress–displacement relationship of the subgrade. Under an oblique pull, high normal stresses develop on stronger subgrades, thus mobilizing high shearing resistance at the reinforcement–soil interface. The higher the bearing resistance of the subgrade, the higher the horizontal component of pullout force and the lower the end displacement of the reinforcement. On the other hand, the end displacement at pullout can become very high for weaker subgrades especially at high values of the angle of obliquity. Also, the pullout capacity under oblique loading for weaker subgrades may approach or even fall below the axial pullout capacity at high values of the angle of obliquity. These adverse pullout responses owing to a low value of bearing resistance of subgrade are magnified when the subgrade stiffness is also small. On weaker subgrades, improvement in angle of interface shear is not advisable as this leads to further reduction in the pullout force and increase in the end displacement. Results are compared with back analysis of published test data on model reinforced soil walls. The comparison suggests that the present model leads to a more rational and better prediction of the pullout failure. Copyright © 2008 John Wiley & Sons, Ltd.     

Consolidation analysis for partly saturated clay by using an elastic–plastic effective stress–strain model

The theory of consolidation is extended to partly saturated clay soils, and formulated for finite element analyses. This formulation couples the effects of both stress and flow. It takes account of variations of this permeability of the soil and compressibility of the pore fluid with changes in void ratio, and the non-linear stress–strain behaviour of soil. The Cam Clay model is revised to model the stress–strain behaviour of compacted soils. The compressibility of pore fluid is derived using Boyle's Law and Henry's Law, taking into account the effect of surface tension. An empirical equation is developed for permeability of pore fluid. An example of settlement of a footing on partly saturated soil is described and discussed.     

Numerical modelling of competitive components transport with non‐linear adsorption

A characteristic finite element (CFE) algorithm for modelling contaminant transport problems coupled with non‐linear competitive adsorption is presented. An alternative algorithm, termed as the transport‐equilibrium Petrov–Galerkin (TEPG) methods in this paper, is also presented for comparison. Through analyses of a number of examples with Peclet number ranging from zero to infinity, it is shown that the CFE algorithm is very competitive with the middle–point TEPG method in terms of accuracy, stability and efficiency. The fully explicit and fully implicit TEPG methods are found to be less appropriate for transport problems coupled with non‐linear equilibrium equations. Copyright © 2000 John Wiley & Sons, Ltd.     

ACMEG‐TS: A constitutive model for unsaturated soils under non‐isothermal conditions

This paper introduces an unconventional constitutive model for soils, which deals with a unified thermo‐mechanical modelling for unsaturated soils. The relevant temperature and suction effects are studied in light of elasto‐plasticity. A generalized effective stress framework is adopted, which includes a number of intrinsic thermo‐hydro‐mechanical connections, to represent the stress state in the soil. Two coupled constitutive aspects are used to fully describe the non‐isothermal behaviour. The mechanical constitutive part is built on the concepts of bounding surface theory and multi‐mechanism plasticity, whereas water retention characteristics are described using elasto‐plasticity to reproduce the hysteretic response and the effect of temperature and dry density on retention properties. The theoretical formulation is supported by comparisons with experimental results on two compacted clays. Copyright © 2008 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.     

Capturing the complete stress–strain behaviour of jointed rock using a numerical approach

This paper presents the results of a series of numerical experiments using the synthetic rock mass (SRM) approach to quantify the behaviour of jointed rock masses. Field data from a massive sulphide rock mass, at the Brunswick mine, were used to develop a discrete fracture network (DFN). The constructed DFN model was subsequently subjected to random sampling whereby 40 cubic samples, of height to width ratio of two, and of varying widths (0.05 to 10 m) were isolated. The discrete fracture samples were linked to 3D bonded particle models to generate representative SRM models for each sample size. This approach simulated the jointed rock mass as an assembly of fractures embedded into the rock matrix. The SRM samples were submitted to uniaxial loading, and the complete stress–strain behaviour of each specimen was recorded. This approach provided a way to determine the complex constitutive behaviour of large‐scale rock mass samples. This is often difficult or not possible to achieve in the laboratory. The numerical experiments suggested that higher post‐peak modulus values were obtained for smaller samples and lower values for larger sample sizes. Furthermore, the observed deviation of the recorded post‐peak modulus values decreased with sample size. The ratio of residual strength of rock mass samples per uniaxial compressive strength intact increases moderately with sample size. Consequently, for the investigated massive sulphide rock mass, the pre‐peak and post‐peak representative elemental volume size was found to be the same (7 × 7 × 14 m). Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of localization failure with strain‐gradient‐enhanced damage mechanics

Strain gradient implies an important characteristic in localized damage deformation, which can be observed in the softening state of brittle materials, and strain gradients constitute the basic behaviours of localization failure area of the materials. The most important point in strain gradient is its damaging function including an internal length scale, which can be used to express the scale effects of mechanical responses of brittle rock mass. By extending the strain gradient theory and introducing an intrinsic material length scale into the constitutive law, the authors develop an isotropic damage model as well as a micro‐crack‐based anisotropic damage model for rock‐like materials in this paper. The proposed models were used to simulate the damage localization under uniaxial tension and plain strain compression, respectively. The simulated results well illustrated the potential of these models in dealing with the well‐known mesh‐sensitivity problem in FEM. In the computation, elements with C₁ continuity have been implemented to incorporate the proposed models for failure localization. When regular rectangle elements are encountered, the coupling between finite difference method (FDM) and conventional finite element method (FEM) is used to avoid large modification to the existing FEM code, and to obtain relatively higher efficiency and reasonably good accuracy. Application of the anisotropic model to the 3D‐non‐linear FEM analysis of Ertan arch dam has been conducted and the results of its numerical simulation coincide well with those from the failure behaviours obtained by Ertan geophysical model test. In this paper, new applications of gradient theories and models for a feasible approach to simulate localized damage in brittle materials are presented. Copyright © 2002 John Wiley & Sons, Ltd.