Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of deformation response and chemo‐mechanical coupling in chalk

A constitutive relation is derived for describing the mechanical response of chalk. The approach is based on a phenomenological framework which employs chemo‐plasticity. The properties of the material are assumed to be affected by the physico‐chemical processes that occur through the interaction between the skeleton and the pore fluid. The underlying mechanism is discussed by invoking a micromechanical analysis. The performance of the framework is illustrated by examining the evolution of mechanical characteristics in the presence of different pore fluids. Copyright © 2006 John Wiley & Sons, Ltd.     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.     

Impact evaluation on groundwater of the infiltration of an urban wastewater geopurification system placed over a detrital aquifer in a semiarid region (Spain)

Over a period of 4 years and 4 months, the geopurification installations at Dehesas de Guadix (Granada, Spain) were monitored to determine the impact on soil and groundwater of the controlled discharge of urban wastewater, and also to identify the best indicators of the entry of the recharged water into the aquifer. The installations are located in an area where the climate is Mediterranean sub-arid, with an average precipitation of less than 287 mm/year, and a rate of evapotranspiration that is almost three times greater. The system was controlled by determining the balance of majority nutrients and boron in the soil and in the groundwater, both at the points affected directly by the wastewater discharge and at others. The quantity of mass discharged was relatively large (COD 14,656 g/m², NO₃ 85 g/m², NO₂ 4 g/m², NH₄ 2,425 g/m², PO₄ 1,143 g/m², K 1,531 g/m², B 63 g/m²). It was observed that the elimination of nutrients within the soil (COD 97.5%, PO₄ 94.4%, K 59.17%, N _total 18.8%, B 12.69%) was very efficient except for the nitrogen, which nevertheless did not reach the groundwater, as it was eliminated at deep levels of the unsaturated zone. Only 12.69% of the boron was removed, and appreciable, increasing amounts of this element did reach the groundwater. Unexpectedly, none of the majority nutrients behaved as a reliable indicator of the impact on groundwater; despite this, the boron and the bicarbonate did clearly reflect the arrival of the recharged water, and are proposed as the best indicators.     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

基于神经网络的平底结构砰击压力预报

对利用神经网络预报平底结构入水砰击压力的方法进行了探讨。首先利用仿真软件计算各种情况下平底结构入水所产生的砰击压力，以此形成训练神经网络的数据集。其次利用数据集对三层反馈式网络进行了训练，讨论了不同隐含层节点数对该非线性系统的拟合能力，并且对梯度下降法、动量修正法和基于优化的LM算法的有效性和精度进行了比较，最后得出了适合平底结构入水砰击系统的网络结构。     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.     

Influence of large-diameter pipe pile driving on surrounding marine soils and adjacent submarine pipelines

Abstract

With the large-scale development and utilization of ocean resources and space, it is inevitable to encounter existing submarine facilities in pile driving areas, which necessitates a safety assessment. In this article, by referring to a wharf renovation project as a reference, the surrounding soil response and buried pipe deformation during pile driving in a near-shore submarine environment are investigated by three-dimensional (3D) numerical models that consider the pore water effect. Numerical studies are carried out in two different series: one is a case of a single pile focusing on the effect of the minimum plane distance of the pile–pipe, and the other is a case of double piles focusing on the effect of the pile spacing.     

A higher-order σ-coordinate non-hydrostatic model for nonlinear surface waves

A higher-order non-hydrostatic model in a σ-coordinate system is developed. The model uses an implicit finite difference scheme on a staggered grid to simultaneously solve the unsteady Navier-Stokes equations (NSE) with the free-surface boundary conditions. An integral method is applied to resolve the top-layer non-hydrostatic pressure, allowing for accurately resolving free-surface wave propagation. In contrast to the previous work, a higher-order spatial discretization is utilized to approximate the large horizontal pressure gradient due to steep surface waves or rapidly varying topographies. An efficient direct solver is developed to solve the resulting block hepta-diagonal matrix system. Accuracy of the new model is validated by linear and nonlinear standing waves and progressive waves. The model is then used to examine freak (extreme) waves. Features of downshifting focusing location and wave asymmetry characteristics are predicted on the temporal and spatial domains of a freak wave.