Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Modelling the poroelastoplastic behaviour of soils subjected to internal erosion by suffusion

Granular soils subjected to seepage flow may suffer suffusion, ie, a selective internal erosion. Extending the classical approach of poromechanics, we deduce a new form of the Clausius-Duhem inequality accounting for dissipation due to suffusion, and we deduce restrictions on the constitutive laws of the soil. We suggest (a) a possible coupling between the seepage forces and the suffusion kinetics and (b) an extension of an existing elastoplastic model for the skeleton mechanical behaviour. Numerical integrations of the elastoplastic model are carried out under drained axisymmetric triaxial and oedometric conditions. As a result, we prove that the extended model is able to qualitatively reproduce the suffusion induced strains and the strength reduction experimentally observed. Predictions on the oedometric behaviour of suffusive soils are also provided.     

低轨道人造卫星(CHAMP、GRACE、GOCE)与高精度地球重力场——卫星重力大地测量的最新发展及其对地球科学的重大影响

评述了卫星重力大地测量的最新发展及其对地球科学的重大影响。为了更好地理解地球内部物理构造与海洋动力学，以及大陆，冰川和海洋的相互作用，改善现有地球重力场模型（包括精度和空间解析度）是非常重要的。IUGG等国际组织对此已经强调了很多年。最近，由德国的GFZ（GeoForschungsZentrum)，美国的NASA（National Aeronautics and Space Adminitration)以及欧洲宇航局ESA（European Space Agency)开发研制了最先进的地球监测技术－SST（Satellite-to-Sateilite Tracking)。其主要特点是利用现有的GPS连续追踪新发射低轨道卫星，并由低轨道卫星对地球重力场作精密观测。已经发射和即将发射的卫星有3颗：GHAMP（Challenging Mini-Satellite Payload for Geophysical Research an Application)已经于2000年发射；GRACE（Gravity Recovery and Climate Experimert)定于2002年发射；GOCE（Gravity Field and Steady-state Ocean Cirulation Explorer)计划2004年发射，它们可以统称为重力卫星。载有SST技术的人造卫星的主要目的是获得具有前所未有的高精度和高空间解析度的全球重力场和大地水准面模型，加强人们对地球内部构造的理解并为海洋和气象研究提供更好地参考。上述3个重力卫星工作在有明显区别的不同波谱内，它们有不同的科学应用，仅有一小部分重合。所以，就应用而言它们是完全互补的。它们在地球科学中的应用将是广泛的，特别对于固体地球物理学，海洋学以及大地测量学等领域，它们将会带来革命性的变化，其意义不亚于GPS。     

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.     

亚毫米波射电天文用超导SIS接收机的研制(英文)

为了研制亚毫米波射电天文用超导SIS(超导 -绝缘体 -超导 )接收机 ,我们重点开展了如下研究 ,1 )Nb超导SIS结在其能隙频率附近的量子混频特性 ,及其结合高能隙超导薄膜 (NbTiN)和高电导率金属薄膜 (Al)分布结阵在 780 - 950GHz频率区间的量子混频特性 ;2 )亚毫米波超导混频器嵌入阻抗的数值和实验表征 ;3)高电流密度小面积Nb超导SIS结的制备和特性表征 ;4)一个 60 0 - 72 0GHz超导SIS混频器的研制和特性表征。本文详细介绍了相关的数值分析和实验测量结果。     

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.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

PROPAGATION OF 30—60 DAY LOW FREQUENCY OSCILLATIONS AND THEIR INFLUENCE UPON THE SUBTROPICAL WESTERLIES JET STREAM DURING NORTHERN HEMISPHERE WINTER

Based on daily ECMWF gridpoint data of two winters during 1981—1983 including an ENSOyear,propagation of low frequency oscillations(LFO)during Northern Hemisphere winters andtheir influences upon 30—60 day oscillations of the subtropical jet stream are studied with the sta-tistical methods as complex empirical orthogonal function(CEOF)and so on.Results show that inthe winter of a normal year(1981—1982),30—60 day oscillations in the subtropical zone aremainly in the northern and southern flanks of exit region of jet stream.In the ENSO year(1982—1983),they are mainly in the vicinity of entrance and exit regions of jet stream.Intraseasonalchanges of subtropical jet stream manifested themselves as latitudinal fluctuation or longitudinalprogression or regression of about 40 day period.There are marked differences between propagat-ing passages of low frequency modes responsible for changes of subtropical jet stream in the normalyear(1981—1982)and in the ENSO year(1982—1983).Changes of oscillation amplitude showobvious phases.In general,the one in late winter is stronger than that in early winter,strongestone occurs in February.     

Evolution of 2:1 layered silicates in low-grade metamorphosed Liassic shales of Central Switzerland

White mica from the Liassic black shales and slates in Central Switzerland was analysed by transmission electron microscopy (TEM) and electron microprobe to determine its textural and compositional evolution during very low-grade prograde metamorphism. Samples were studied from the diagenetic zone, anchizone and epizone (T ≈100°–450 °C). Phyllosilicate minerals analysed include illite/smectite (I/S), phengite, muscovite, brammallite, paragonite, margarite and glauconite. Textural evolution primarily is towards larger, more defect-free grains with compositions that approach those of their respective end-members. The smectite-to-illite transformation reduced the amounts of the exchange components SiK_?1Al_?1, MgSiAl_?2, and Fe³⁺Al_?1. These trends continue to a lesser degree in the anchizone and epizone. Correlations between the proportion of smectite in I/S and the composition of I/S indicate that smectite layers may contain a high layer charge. Illite in I/S bears a compositional resemblance to macrocrystalline phengite in some samples, but is different in others. Paragonite first appears in the upper diagenetic zone or lower anchizone as an interlayer-deficient brammallite, and it may be mixed with muscovite on the nanometre scale. Owing to the small calculated structure factor for paragonite-muscovite superstructures, conventional X-ray powder diffraction cannot distinguish between mixed-layer structures and a homogeneous compositionally intermediate solid solutions. However, indirect TEM evidence shows that irregularly shaped domains of Na- and K-rich mica exist below 10 nm. Subsequent coarsening of domains at higher grades produced discrete paragonite grains at the margins of muscovite crystals or in laths parallel to the basal plane of the host muscovite. Margarite appears in the epizone and follows a textural evolution similar to paragonite in that mixtures of margarite, paragonite, and muscovite may initially occur on the nanometre scale. However, no evidence of interlayer-poor margarite has been found.