FINITE ELEMENT MODELLING FOR LAND DISPLACEMENTS DUE TO PUMPING

Equations of equilibrium (force balance) and flow in multidimensions were coupled in this paper to describe land displacements due to pressure decline in aquifers. A Galerkin finite element model based on these equations was developed. The saturated/unsaturated behaviour and the isotropic/anisotropic properties of permeability and elasticity were considered when the model was formulated. This model was verified by comparing its simulation results with those of known analytical solutions for simplified cases. The simulation of displacements due to pressure decline in unsaturated media was also performed. Those results demonstrated that the choice of boundary ranges for an aquifer with infinite domain may significantly affect the estimated horizontal and vertical displacements. To obtain a good estimation of land displacements, the boundary ranges should be carefully chosen. The displacements occurring in unconfined aquifers are caused by the drop of the water table and the change in body force in the dewatering zone. Simulation results also indicated that the change in body force should be considered once an unconfined aquifer has been pumped. Otherwise, the horizontal and vertical displacements in unconfined aquifers would be overestimated and underestimated, respectively. The behaviour of land displacements due to pumping was shown to be affected by changes in the total stresses in aquifers.     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

Multiscale approach to geo‐composite cellular structures subjected to rock impacts

Geo‐composite cellular structures are an efficient technological solution for various applications in civil engineering. This type of structure is particularly well adapted to resisting rockfalls and can act as a defensive structure. However, the design of such structures is for the most part empirically based; this lack of research‐based design stagnates optimization and advanced development. In this paper, the mechanical behaviour of a geo‐composite cellular structure is investigated using a multi‐scale approach, from the individual cell made up of an assembly of rocky particles contained in a wire netting cage to the entire structure composed of a regular array of cells. Based on discrete modelling of both the cell and structure scales, a computational tool has been developed for design purposes. Copyright © 2007 John Wiley & Sons, Ltd.     

Importance of footing–soil separation on dynamic stiffness of piled embedded footings

Different phenomena such as soil consolidation, erosion, and scour beneath an embedded footing supported on piles may lead to loss of contact between soil and the pile cap underside. The importance of this separation on the dynamic stiffness and damping of the foundation is assessed in this work. To this end, a numerical parametric analysis in the frequency domain is performed using a rigorous three‐dimensional elastodynamic boundary element–finite element coupling scheme. Dimensionless plots relating dynamic stiffness functions computed with and without separation effects are presented for different pile–soil configurations. Vertical, horizontal and rocking modes of oscillation are analyzed for a wide range of dimensionless frequencies. It is shown that the importance of separation is negligible for frequencies below those for which dynamic pile group effects start to become apparent. Redistribution of stiffness contributions between piles and footing is also addressed. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of frictional forces acting on sidewalls of buried box culverts

The effects of frictional forces acting on the sidewalls of buried box culverts are presented as determined with finite element method (FEM) and detailed soil modelling. The possibility of reducing earth pressure on deeply buried concrete box culverts by the imperfect trench installation (ITI) method has been contemplated during the last several decades. There have been limited research results published primarily regarding the qualitative aspect of load reduction in ITIs. It was found during the course of this study that significant frictional forces develop along the sidewalls of box culverts and adjacent sidefills in ITIs. Current American Association of State Highway and Transportation Officials provisions do not consider these frictional forces, but they cannot be neglected in ITIs, as their effect is dominant. An optimum geometry for the soft zone in ITIs is presented to maximize earth load reductions. The soil–structure interaction at the box culvert–soil interface was found to have a significant effect on total earth pressure acting on the bottom slab. Predictor equations for earth load reduction rates were formulated for ITIs incorporating the optimum soft zone geometry based on the FEM. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.     

White dwarf + red dwarf binary systems and the period gap of cataclysmic variables

Cataclysmic variables (CV) and pre-cataclysmic binaries (PCB) are discussed. The main difference between them (accretion or its absence) is shown to be a consequence of the evolutionary process and the properties of their progenitors. Both types of system have a bimodal distribution of their periods, but the extrema are in counterphase. Luminosity-effective temperature diagrams for the primary components are used to show that both systems have approximately the same age, which conflicts with the notion of PCBs as precursors of CVs. Calculations of the maximum distance between components for which the system remains stable show that CVs have passed this limit, while PCBs maintain their stability during this evolution. It is suggested that, after ejecting a common shell, future CVs immediately become semi-detached systems. It this is so, then there must be cataclysmic variables which are the central stars of planetary nebulae. __________ Translated from Astrofizika, Vol. 50, No. 3, pp. 427–439 (August 2007).     

Two‐dimensional physically based finite element runoff model for small agricultural watersheds: I. Model development

Soil erosion by water is the root cause of ecological degradation in the Shiwalik foothills of Northern India. Simulation of runoff and its component processes is a pre‐requisite to develop the management strategies to tackle the problem, successfully. A two‐dimensional physically based distributed numerical model, ROMO2D has been developed to simulate runoff from small agricultural watersheds on an event basis. The model employs the 2‐D Richards equation with sink term to simulate infiltration and soil moisture dynamics in the vadoze zone under variable rainfall conditions, and 2‐D Saint‐Venant equations under the kinematic wave approximation along with Manning's equation as the stage‐discharge equation for runoff routing. The various flow‐governing equations have been solved numerically by employing a Galerkin finite element method for spatial discretization using quadrilateral elements and finite difference techniques for temporal solutions. The ROMO2D computer program has been developed as a class‐based program, coded in C + + in such a way that with minor modifications, the model can be used to simulate runoff on a continuous basis. The model writes output for a runoff hydrograph of each storm. Model development is described in this paper and the results of model testing and field application are to be presented in a subsequent paper. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of an improved semi-Lagrangian procedure to fully-nonlinear simulation of sloshing in non-wall-sided tanks

The semi-Lagrangian procedure is widely used for updating the fully-nonlinear free surface in the time domain. However, this procedure is only available to cases when the body surface is vertical near the waterline. Present study introduces an improved semi-Lagrangian procedure which removes this ‘vertical-wall’ limitation. Coupling with the boundary element method, the improved semi-Lagrangian procedure is applied to the simulation of fully-nonlinear sloshing waves in non-wall-sided tanks. From the result comparison with the open source CFD software OpenFOAM, it is confirmed that this numerical scheme could guarantee a sufficient accuracy. Further series studies on 2D and 3D fully-nonlinear sloshing waves in wedged tanks are performed. Featured phenomena are observed which are distinct from those in wall-sided tanks.     

基于多尺度有限元的海上升压站脚靴式桩套筒连接结构分析研究

脚靴式海上升压站导管架的桩套筒与主导管的连接结构是关键荷载传递部位。整体梁系有限元模型无法反应局部板型结构的失效状态;截断板系有限元模型则受制于边界条件的合理性,结果不够精确。文中基于某工程海上升压站,建立不同精细化范围的多尺度有限元模型及不同边界条件的截断模型,对结构屈服失效及屈曲失效复核的同时论证了局部网格对计算精细程度的影响及截断模型对边界条件的敏感性;此外,研究了轭板板厚及应力比对屈曲失效的敏感程度,为今后海上升压站脚靴结构不同阶段的设计、优化提供参考和借鉴。