Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

Liquefaction potential evaluations by energy-based method and stress-based method for various ground motions

An energy-based liquefaction potential evaluation method (EBM) previously developed was applied to a uniform sand model shaken by seismic motions recorded at different sites during different magnitude earthquakes. It was also applied to actual liquefaction case histories in Urayasu city during the 2011 M9.0 Tohoku earthquake and in Tanno-cho during the 2003 M8.0 Tokachi-oki earthquake. In all these evaluations, the results were compared with those by the currently used stress-based method (SBM) under exactly the same seismic and geotechnical conditions. It was found that EBM yields similar results with SBM for several ground motions of recent earthquakes but has easier applicability without considering associated parameters. In Urayasu city, the two methods yielded nearly consistent results by using an appropriate coefficient in SBM for the M9.0 earthquake, though both overestimated the actual liquefaction performance, probably because effects of plasticity and aging on in situ liquefaction strength were not taken into account. In Tanno-cho, EBM could evaluate actual liquefaction performance due to a small-acceleration motion during a far-field large magnitude earthquake while SBM could not.     

The multi-component Hekla Ö Tephra,Iceland: a complex widespread mid-Holocene tephra layer

Large Plinian eruptions from Hekla volcano, Iceland, produce compositionally zoned tephra used as key markers in tephrochronology. However, spatial variations in chemical composition of a tephra layer may complicate its identification. An example is the 5950–6180 cal a bp Hekla Ö tephra layer, which shows compositional spread from rhyolite, dacite and andesite to basalt. In soil sections north of Hekla, the SiO₂ content of the tephra glass reaches 76 wt% in the lowest unit of the Hekla Ö deposit and decreases to 62–63 wt% in the uppermost unit. Intermingled within the whole deposit are basalt tephra grains having 46–47 wt% SiO₂. The composition of the basalt glass includes primitive basalt and a more evolved basalt (MgO >6 and <6 wt%, respectively). Together with literature data, the Hekla Ö tephra and the so-called T-Tephra/Hekla-T are most likely from contemporaneous eruptions of different vents on the Hekla volcanic system, forming a single important marker tephra (Hekla ÖT) deposited over 80% of Iceland. Identification is complicated by its spatial compositional heterogeneity, such as systematic decrease in SiO₂ content from the east to the west of Hekla volcano. Consequently, an individual tephra layer from a large explosive eruption can have different composition at different locations. © 2020 John Wiley & Sons, Ltd.     

Discrete element model for crack propagation in brittle materials

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

A cavity expansion–based solution for interpretation of CPTu data in soils under partially drained conditions

A cavity expansion–based solution is proposed in this paper for the interpretation of CPTu data under a partially drained condition. Variations of the normalized cone tip resistance, cone factor, and undrained-drained resistance ratio are examined with different initial specific volume and overconsolidation ratio, based on the exact solutions of both undrained and drained cavity expansion in CASM, which is a unified state parameter model for clay and sand. A drainage index is proposed to represent the partially drained condition, and the critical state after expansion and stress paths of cavity expansion are therefore predicted by estimating a virtual plastic region and assuming a drainage-index–based mapping technique. The stress paths and distributions of stresses and specific volume are investigated for different values of drainage index, which are also related to the penetration velocity with comparisons of experimental data and numerical results. The subsequent consolidation after penetration is thus predicted with the assumption of constant deviatoric stress during dissipation of the excess pore pressure. Both spherical and cylindrical consolidations are compared for dissipation around the cone tip and the probe shaft, respectively. The effects of overconsolidation ratio on the stress paths and the distributions of excess pore pressure and specific volume are then thoroughly investigated. The proposed solution and the findings would contribute to the interpretation of CPTu tests under a random drained condition, as well as the analysis of pile installation and the subsequent consolidation.     

Material selection for the gas containment system of a compressed natural gas carrier fleet

In this work, the possible exploitation of fiber-reinforced composites in the context of maritime transportation of compressed natural gas (CNG) is investigated. In addition to a more conventional steel configuration, two different fiber materials, carbon and glass, are considered as construction materials for pressure vessels (PVs) to be stored on board ships, with thickness optimized by FEM analysis.The considered scenario is represented by the transportation of CNG from an offshore well to a terminal on shore. Fleets of ships carrying CNG in pressure vessels manufactured with the investigated materials are generated by means of a ship synthesis model (SSM) software and compared on the basis of technical and economical indicators.The choice of the construction material influences considerably the weight of the PVs, which represent a major item of total ship weight and reflects directly on the general transport performances in terms of resistance, seakeeping and reliability in the service. On the other hand, capital as well as operating expenditures are considerably affected by the choice. When exploring the design space, the ship synthesis model is able, at a preliminary stage of the design, to account for the various technical and economical aspects, their implications and relationships. Results are presented of computations carried out in a specific case, identified by the annual gas production and other characteristics of the well terminal and a cruising route for the ships. The comparison is carried out on the basis of the cost per transported unit of gas and of the percentage of success in the transportation process. The computations show that the choice of the PV material has a key influence on the results in terms of optimal number, dimensions and speed of the ships.     

Modeling the 2D behavior of dry‐stone retaining walls by a fully discrete element method

The dry‐stone retaining walls (DSRW) have been tipped as a promising solution for sustainable development. However, before recently, their behavior is relatively obscure. In this study, discrete element method (DEM) approach was applied to simulate the plane strain failure of these walls. A commercial DEM package (PFC2D™) was used throughout this study. The authors used a fully discrete approach; thus, both the wall and the backfill were modeled as discrete elements. The methodology for obtaining the micromechanical parameters was discussed in detail; this includes the three mechanical sub‐systems of DSRWs: wall, backfill and interface. The models were loaded progressively until failure, and then the results were compared with the full‐scale experimental results where the walls were loaded, respectively, with hydrostatic load and backfill. Despite its complexity and its intensive calculation time, DEM model can then be used to validate a more simplified approach. Copyright © 2015 John Wiley & Sons, Ltd.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.