Modelling of ground motion in the vicinity of massive structures

A two-dimensional elastic Chebyshev spectral element method (SPEM) is used to model the seismic wavefield within a massive structure and in its vicinity. We consider 2-D models where a linear elastic structure, with quadrangular cross-section, resting on an elastic homogeneous half-space, is impinged upon by the waves generated by a surface impulse at some distance. The scattering of Rayleigh waves and the response of the structure are extensively analysed in a parametric way, varying size, mechanical parameters and shape of the load. Some of the models considered are representative of embankments and earth dams. The simulation shows that some models resonate, storing part of the incoming energy. With realistic parameters, the lowest resonance frequency is due to pure shear deformation and is controlled by the shear velocity and height of the load. Flexural modes are excited only at higher frequencies. The acceleration at the top of the structure may be five/seven times higher than at the base, depending on the mass of the structure. The gradual release of trapped energy produces a ground roll lasting several seconds after the wave front has passed. The ground-roll amplitude depends on the sturcture's mass and can be as large as 30% of the peak acceleration. Outside resonance conditions, the ground motion is almost unaffected by the presence of the artefact; the horizontal motion on top of it is nearly twice the motion at ground level. Similar results should be expected when the incident field is an upcoming shear wave. A qualitative discussion shows that the presence of anelastic attenuation in the embankment does not significantly alter the preceding conclusions, unless it is of very low values (e.g. Q < 15).The modelling results that we discuss indicate that the soil-structure interaction may substantially alter the ‘free-field’ ground motion. From a practical point of view, the main conclusions are: (1) careful analysis is necessary when interpreting seismic records collected in the vicinity of large artefacts; (2) seismic hazard at a site may depend on the presence of man-made structures such as embankments, dams, tall and massive buildings.     

A numerical method for the study of shear band propagation in soft rocks

This paper investigates the possibility of interpreting progressive shear failure in hard soils and soft rocks as the result of shear propagation of a pre‐existing natural defect. This is done through the application of the principles of fracture mechanics, a slip‐weakening model (SWM) being used to simulate the non‐linear zone at the tips of the discontinuity. A numerical implementation of the SWM in a computation method based on the boundary element technique of the displacement discontinuity method (DDM) is presented. The crack and the non‐linear zone at the advancing tip are represented through a set of elements, where the displacement discontinuity (DD) in the tangential direction is determined on the basis of a friction law. A residual friction angle is assumed on the crack elements. Shear resistance decreases on elements in the non‐linear zone from a peak value at the tip, which is characteristic of intact material, to the residual value. The simulation of a uniaxial compressive test in plane strain conditions is carried out to exemplify the numerical methodology. The results emphasize the role played by the critical DD on the mechanical behaviour of the specimen. A validation of the model is shown through the back analysis of some experimental observations. The results of this back analysis show that a non‐linear fracture mechanics approach seems very promising to simulate experimental results, in particular with regards to the shear band evolution pattern. Copyright © 2009 John Wiley & Sons, Ltd.     

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

Fundamental physics and absolute positioning metrology with the MAGIA lunar orbiter

MAGIA is a mission approved by the Italian Space Agency (ASI) for Phase A study. Using a single large-diameter laser retroreflector, a large laser retroreflector array and an atomic clock onboard MAGIA we propose to perform several fundamental physics and absolute positioning metrology experiments: VESPUCCI, an improved test of the gravitational redshift in the Earth?CMoon system predicted by General Relativity; MoonLIGHT-P, a precursor test of a second generation Lunar Laser Ranging (LLR) payload for precision gravity and lunar science measurements under development for NASA, ASI and robotic missions of the proposed International Lunar Network (ILN); Selenocenter (the center of mass of the Moon), the determination of the position of the Moon center of mass with respect to the International Terrestrial Reference Frame/System (ITRF/ITRS); this will be compared to the one from Apollo and Lunokhod retroreflectors on the surface; MapRef, the absolute referencing of MAGIA??s lunar altimetry, gravity and geochemical maps with respect to the ITRF/ITRS. The absolute positioning of MAGIA will be achieved thanks to: (1) the laboratory characterization of the retroreflector performance at INFN-LNF; (2) the precision tracking by the International Laser Ranging Service (ILRS), which gives two fundamental contributions to the ITRF/ITRS, i.e. the metrological definition of the geocenter (the Earth center of mass) and of the scale of length; (3) the radio science and accelerometer payloads; (4) support by the ASI Space Geodesy Center in Matera, Italy. Future ILN geodetic nodes equipped with MoonLIGHT and the Apollo/Lunokhod retroreflectors will become the first realization of the International Moon Reference Frame (IMRF), the lunar analog of the ITRF.     

Modelling mechanical behaviour of limestone under reservoir conditions

High porosity and low permeability limestone has presented pore collapse. As fluid is withdrawn from these reservoirs, the effective stresses acting on the rock increase. If the strength of the rock is overcome, pore collapse may occur, leading to irreversible compaction of porous media with permeability and porosity reduction. It impacts on fluid withdrawal. Most of reservoirs have been discovered in weak formations, which are susceptible to this phenomenon. This work presents a study on the mechanical behaviour of a porous limestone from a reservoir located in Campos Basin, offshore Brazil. An experimental program was undergone in order to define its elastic plastic behaviour. The tests reproduced the loading path conditions expected in a reservoir under production. Parameters of the cap model were fitted to these tests and numerical simulations were run. The numerical simulations presented a good agreement with the experimental tests. Copyright © 2006 John Wiley & Sons, Ltd.     

A cosmological study of the star formation history in the solar neighbourhood

We use a cosmological galactic evolutionary approach to model the Milky Way. A detailed treatment of the mass aggregation and dynamical history of the growing dark halo is included, together with a self-consistent physical treatment for the star formation processes within the growing galactic disc. This allows us to calculate the temporal evolution of star and gas surface densities at all galactic radii, in particular, the star formation history (SFH) at the solar radius. A large range of cosmological mass aggregation histories (MAHs) is capable of producing a galaxy with the present-day properties of the Milky Way. The resulting SFHs for the solar neighbourhood bracket the available observational data for this feature, the most probable MAH yielding the optimal comparison with these observations. We also find that the rotation curve for our Galaxy implies the presence of a constant density core in its dark-matter halo.     

Stress–dilatancy based modelling of granular materials and extensions to soils with crushable grains

Stress–dilatancy relations have played a crucial role in the understanding of the mechanical behaviour of soils and in the development of realistic constitutive models for their response. Recent investigations on the mechanical behaviour of materials with crushable grains have called into question the validity of classical relations such as those used in critical state soil mechanics. In this paper, a method to construct thermodynamically consistent (isotropic, three‐invariant) elasto‐plastic models based on a given stress–dilatancy relation is discussed. Extensions to cover the case of granular materials with crushable grains are also presented, based on the interpretation of some classical model parameters (e.g. the stress ratio at critical state) as internal variables that evolve according to suitable hardening laws. Copyright © 2004 John Wiley & Sons, Ltd.     

Correlation of Plio-Pleistocene deposits of the Lower Rhine Basin (North-West Germany) and the Valle Ricca pits (Central Italy)

Biostratigraphic and palaeomagnetic research has been carried out on selected profiles in North-Western Germany and Central Italy in order to correlate Plio-Pleistocene sections.Around the Réunion subchronozone, vegetation in the Valle Ricca was dominated by a mountain forest type pointing to cool climatic conditions. In the Lower Rhine Basin vegetation had a slightly warmer character, but was also dominated strongly by a mountain type of forest vegetation, characteristic of cooler summers and higher precipitation than at present.Two warmer periods have been determined in the Central Italian profiles; the youngest around the Olduvai reversal, should correspond to the youngest Tiglian beds of the Lower Rhine Basin.Deposits older than the Réunion subzone have not been found in the Valle Ricca. In the Lower Rhine Basin, on the other hand, the Gauss-Matuyama boundary is located between the uppermost Pliocene, the Reuverian C and the Praetiglian.     

Current confinement in solar coronal loops

In this paper solar coronal loops are regarded as regions of localized current flows. The main purpose is to investigate the consequences of current confinement rather than to produce a model. The physical and observational basis for this assumption are presented as well as the connection with previous studies on loop structure. A proper choice of the current profile allows us to treat quantitatively the equilibrium structure of the loops and their MHD and resistive stability properties. Regions of absolute stability against ideal kink modes are found. Explicit growth rates for the tearing-mode instability are computed. The possible relevance of other resistive effects is also discussed and the crucial importance of the small-scale geometry of the magnetic field outlined.