A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Developing constitutive models from EPR‐based self‐learning finite element analysis

A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.     

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.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Morphology and Luminosity Segregations in Galaxy Groups

We study morphology and luminosity segregation of galaxies in groups. We analyze the two catalogs of (∼2×400) groups which have been identified in the Nearby Optical Galaxy sample, by means of hierarchical and percolation `friends-of-friends' methods. We find that earlier-type (brighter) galaxies are more clustered and lie closer to the group centers, both in position and in velocity, than later-type (fainter) galaxies. Spatial segregations are stronger than kinematical segregations. These effects are generally detected at the ≳3-sigma level, with the exception of morphological segregation in velocity, which is the weakest effect. Our main results are confirmed by the analysis of statistically more reliable groups (with at least five members), and are strengthened by the detection of segregation in both hierarchical and percolation catalogs. Luminosity segregation is shown to be independent of morphology segregation. Our conclusions agree with a continuum of segregation properties of galaxies in systems, from low-mass groups to massive clusters. This revised version was published online in July 2006 with corrections to the Cover Date.     

Size and shape evolution of pores in a viscoplastic matrix under compression

The frequent use of soils and earth materials for hydraulic capping and for geo‐environmental waste containment motivated our interest in detailed modelling of changes in size and shape of macro‐pores to establish links between soil mechanical behaviour and concurrent changes in hydraulic and transport properties. The objective of this study was to use finite element analysis (FEA) to test and extend previous analytical solutions proposed by the authors describing deformation of a single macro‐pore embedded in linear viscoplastic soil material subjected to anisotropic remote stress. The FEA enables to consider more complex pore geometries and provides a detailed picture of matrix yield behaviour to explain shortcomings of approximate analytical solutions. Finite element and analytical calculations agreed very well for linear viscous as well as for viscoplastic materials, only limited for the case of isotropic remote stress due to the simplifications of the analytical model related to patterns and onset of matrix‐yielding behaviour. FEA calculations were compared with experimental data obtained from a compaction experiment in which pore deformation within a uniform modelling clay sample was monitored using CAT scanning. FEA predictions based on independently measured material properties and initial pore geometry provided an excellent match with experimentally determined evolution of pore size and shape hence lending credence to the potential use of FEA for more complex pore geometries and eventually connect macro‐pore deformation with hydraulic properties. Copyright © 2006 John Wiley & Sons, Ltd.     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

A photometric study of the polar ring galaxy UGC 5600

We analyze our BVR _c photometry for UGC 5600, a candidate polar ring galaxy, obtained with the 6-m telescope. We have confirmed the existence of an inner polar ring and show that the outer ring-shaped structure represents spiral arms; i.e., UGC 5600 belongs to the rare class of gas-rich spiral galaxies with inner polar rings.     

Effect of element size on the static finite element analysis of steep slopes

The accuracy of the computed stress distribution near the free surface of vertical slopes was evaluated in this study as a function of the element size, including aspect ratio. To accomplish this objective, a parametric study was carried out comparing stresses computed using the finite element method (FEM) to those obtained from a physical model composed of photoelastic material. The results of the study indicate a reasonable agreement between a gelatin model and the FEM model for shear stresses, and an overall good agreement between the two models for the principal stresses. For stresses along the top of the slope, the height of the element tends to be more important than width or aspect ratio, at least for aspect ratios up to 4. In all cases, the greatest difference between the two models occurs in the vicinity of the slope. Specifically, if H is defined as the slope height, an element height of H/10 appears to be adequate for the study of stresses deep within the slope, such as for typical embankment analyses. However, for cases where tensile stresses in the vicinity of the slope face which are critical, such as for the stability analysis of steep slopes, element heights as small as H/32, or higher‐order elements, are necessary. Copyright © 2001 John Wiley & Sons, Ltd.