A basal slip model for Lagrangian finite element simulations of 3D landslides          下载免费PDF全文

Massimiliano Cremonesi  Francesco Ferri  Umberto Perego 《国际地质力学数值与分析法杂志》2017,41(1):30-53

A Lagrangian numerical approach for the simulation of rapid landslide runouts is presented and discussed. The simulation approach is based on the so‐called Particle Finite Element Method. The moving soil mass is assumed to obey a rigid‐viscoplastic, non‐dilatant Drucker–Prager constitutive law, which is cast in the form of a regularized, pressure‐sensitive Bingham model. Unlike in classical formulations of computational fluid mechanics, where no‐slip boundary conditions are assumed, basal slip boundary conditions are introduced to account for the specific nature of the landslide‐basal surface interface. The basal slip conditions are formulated in the form of modified Navier boundary conditions, with a pressure‐sensitive threshold. A special mixed Eulerian–Lagrangian formulation is used for the elements on the basal interface to accommodate the new slip conditions into the Particle Finite Element Method framework. To avoid inconsistencies in the presence of complex shapes of the basal surface, the no‐flux condition through the basal surface is relaxed using a penalty approach. The proposed model is validated by simulating both laboratory tests and a real large‐scale problem, and the critical role of the basal slip is elucidated. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Surface flow boundary conditions in modeling land subsidence due to fluid withdrawal     

Baú D  Ferronato M  Gambolati G  Teatini P 《Ground water》2004,42(4):516-525

Land subsidence due to subsurface fluid (water, gas, oil) withdrawal is often predicted by either finite element or finite difference numerical models based on coupled poroelastic theory, where the soil is represented as a semi-infinite medium bounded by the traction-free (ground) surface. One of the variables playing a most important role on the final outcome is the flow condition used on the traction-free boundary, which may be assumed as either permeable or impermeable. Although occasionally justified, the assumption of no-flow surface seems to be in general rather unrealistic. A permeable boundary where the fluid pressure is fixed to the external atmospheric pressure appears to be more appropriate. This paper addresses the response, in terms of land subsidence, obtained with a coupled poroelastic finite element model that simulates a distributed pumping from a horizontal aquifer confined between two relatively impervious layers, and takes either a permeable boundary surface, i.e., constant hydraulic potential, or an impermeable boundary, i.e., a zero Neumann flow condition. The analysis reveals that land subsidence is rather sensitive to the flow condition implemented on the traction-free boundary. In general, the no-flow condition leads to an overestimate of the predicted ground surface settlement, which could even be 1 order of magnitude larger than that obtained with the permeable boundary.  相似文献   

The Vallo di Diano and Auletta extensional basins in the southern Apennines (Italy): a simple model for a complex setting     

Laura Amicucci  Massimiliano R. Barchi  Paola Montone  Nicola Rubiliani 《地学学报》2008,20(6):475-482

The aim of this paper was to provide a significant case‐history concerning the evolution of a segmented system of extensional faults and related basins, investigated by a set of seismic reflection profiles. We investigated two kinematically linked semi‐grabens, developed at the hangingwall of two opposite‐dipping normal faults: the Vallo di Diano and Auletta basins, located in the southern Apennines, one of the most seismically active regions of the Italian peninsula. Our interpretation suggests that the Pliocene–Quaternary tectonic history consists of a single extensional phase, where the major NW–SE trending normal faults, generating the basins, and the strike‐slip and transtensional faults connecting adjacent normal fault segments, simultaneously act, being part of the same extensional system. We also conclude that major normal faults, bordering the basins, should be considered as potential seismogenic sources in the seismic hazard evaluation.  相似文献   

From present rivers to old terrigenous sediments: the evolution of the drainage system in the eastern Southern Alps     

G. Monegato  C. Stefani  M. Zattin 《地学学报》2010,22(3):218-226

  相似文献   

Development of an incised valley‐fill at an evolving rift margin: Pleistocene eustasy and tectonics on the southern side of the Gulf of Corinth,Greece     

Katarina Gobo  Massimiliano Ghinassi  Wojciech Nemec  Eivind Sjursen 《Sedimentology》2014,61(4):1086-1119

This study from the southern margin of the Gulf of Corinth documents a Late Pleistocene incised valley‐fill succession that differs from the existing facies models, because it comprises gravelly shoal‐water and Gilbert‐type deltaic deposits, shows strong wave influence and lacks evidence of tidal activity. The valley‐fill is at least 140 m thick, formed in 50 to 100 ka between the interglacials Marine Isotope Stage 9a and Marine Isotope Stage 7c. The relative sea‐level rise left its record both inside and outside the incised valley, and the age of the valley‐fill is estimated from a U/Th date of coral‐bearing deposits directly outside the palaeovalley outlet. Tectonic up‐warping due to formation of a valley‐parallel structural relay ramp contributed to the valley segmentation and limited the landward extent of marine invasions. The valley segment upstream of the ramp crest was filled with a gravelly alluvium, whereas the downstream segment accumulated fluvio‐deltaic deposits. The consecutive deltaic systems nucleated in the ramp‐crest zone, forming a bathymetric gradient that promoted the ultimate growth of thick Gilbert‐type delta. The case study contributes to the spectrum of conceptual models for incised valley‐fill architecture. Four key models are discussed with reference to the rates of sediment supply and accommodation development, and it is pointed out that not only similarity, but also all departures of particular field cases from these end‐member models may provide valuable information on the system formative conditions. The Akrata incised valley‐fill represents conditions of high sediment supply and a rapid, but stepwise development of accommodation that resulted from the spatiotemporal evolution of normal faulting at the rift margin and overprinted glacioeustatic signals. This study adds to an understanding of valley‐fill architecture and provides new insights into the Pleistocene tectonics and palaeogeography of the Corinth Rift margin.  相似文献   

A numerical study of the size of the homoclinic tangle of hyperbolic tori and its correlation with Arnold diffusion in Hamiltonian systems     

Elena Lega  Massimiliano Guzzo  Claude Froeschlé 《Celestial Mechanics and Dynamical Astronomy》2010,107(1-2):129-144

Using a three degrees of freedom quasi-integrable Hamiltonian as a model problem, we numerically compute the unstable manifolds of the hyperbolic manifolds of the phase space related to single resonances. We measure an exponential dependence of the splitting of these manifolds through many orders of magnitude of the perturbing parameter. This is an indirect numerical verification of the exponential decay of the normal form, as predicted by the Nekhoroshev theorem. We also detect different transitions in the topology of these manifolds related to the local rational approximations of the frequencies. The variation of the size of the homoclinic tangle as well as the topological transitions turn out to be correlated to the speed of Arnold diffusion.  相似文献   

Extended analytical formulas for the perturbed Keplerian motion under a constant control acceleration     

Federico Zuiani  Massimiliano Vasile 《Celestial Mechanics and Dynamical Astronomy》2015,121(3):275-300

  相似文献   

Pulsar sensitivity studies of the GLAST large area telescope     

Massimiliano Razzano  Alice K. Harding 《Astrophysics and Space Science》2007,309(1-4):241-244

In this contribution we present our preliminary investigation on pulsar sensitivity of the Large Area Telescope, the main instrument aboard the GLAST mission. In particular we concentrated our attention to pulsars located at low galactic latitudes. We created a set of simulated pulsars having different fluxes in an array of galactic coordinates separated by a distance greater than the LAT Point Spread Function in order to avoid confusion between adjacent sources. Galactic gamma-ray sky background as used during the second LAT Data Challenge (DC2) is also included. We then run an automatic routine for testing periodicity for all the pulsars considering an opportune timing solution. In this way we can obtain a map of the sensitivity of the periodic searches for different fluxes and for various Galactic latitudes. Some assumptions have been made by simulating the pulsar sources, but this study is a first step toward an estimate for pulsed emission sensitivity of the GLAST LAT. The pulsed flux sensitivity profile we generate could also be input to a population synthesis code of Galactic pulsars in order to obtain more accurate predictions of the number of expected pulsar detections by GLAST. On behalf of the GLAST LAT Collaboration.  相似文献   

The nekhoroshev theorem and the asteroid belt dynamical system     

Alessandro Morbidelli  Massimiliano Guzzo 《Celestial Mechanics and Dynamical Astronomy》1996,65(1-2):107-136

The present paper reviews the Nekhoroshev theorem from the point of view of physicists and astronomers. We point out that Nekhoroshev result is strictly connected with the existence of a specific structure of the phase space, the existence of which can be checked with several numerical tools. This is true also for a degenerate system such as the one describing the motion of an asteroid in the so called main belt. The main difference is that in some parts of the belt, the Nekhoroshev result cannot apply a priori. Mean motion resonances of order smaller than the logarithm of the mass of Jupiter and first order secular resonances must be excluded. In the remaining parts, conversely, the Nekhoroshev theorem can be proved, provided someparameters, such as the masses, the eccentricities and the inclinations of the planets are small enough. At the light of this result, a massive campaign of numerical integrations of real and fictitious asteroids should allow to understand which is the real dynamical structure of the asteroid belt.  相似文献   

Numerical performance of projection methods in finite element consolidation models     

Giuseppe Gambolati  Giorgio Pini  Massimiliano Ferronato 《国际地质力学数值与分析法杂志》2001,25(14):1429-1447

Projection, or conjugate gradient like, methods are becoming increasingly popular for the efficient solution of large sparse sets of unsymmetric indefinite equations arising from the numerical integration of (initial) boundary value problems. One such problem is soil consolidation coupling a flow and a structural model, typically solved by finite elements (FE) in space and a marching scheme in time (e.g. the Crank–Nicolson scheme). The attraction of a projection method stems from a number of factors, including the ease of implementation, the requirement of limited core memory and the low computational cost if a cheap and effective matrix preconditioner is available. In the present paper, biconjugate gradient stabilized (Bi‐ CGSTAB) is used to solve FE consolidation equations in 2‐D and 3‐D settings with variable time integration steps. Three different nodal orderings are selected along with the preconditioner ILUT based on incomplete triangular factorization and variable fill‐in. The overall cost of the solver is made up of the preconditioning cost plus the cost to converge which is in turn related to the number of iterations and the elementary operations required by each iteration. The results show that nodal ordering affects the perfor mance of Bi‐CGSTAB. For normally conditioned consolidation problems Bi‐CGSTAB with the best ILUT preconditioner may converge in a number of iterations up to two order of magnitude smaller than the size of the FE model and proves an accurate, cost‐effective and robust alternative to direct methods. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献