Littérature scientifique sur le sujet « 3D Runout Modeling »

This study presents a field-scale simulation of the Hongshiyan landslide in China. It uses an advanced numerical approach (material point method (MPM)) and a constitutive model (the Drucker–Prager model + μ(I) rheological relation) for the three-dimensional (3D) simulation. The performance of the developed MPM model is validated with laboratory-scale experimental data on granular collapse before being applied to field-scale analyses. ArcGIS data are used to create a 3D MPM model of the soil body with complicated geometry. Although the developed model can describe the multiple phases of granular flow, it focuses on the runout behavior of the landslide in this work. The landslide is assumed to have occurred suddenly due to an earthquake, and global sudden failure rather than progressive failure is modeled. The MPM simulation results match reasonably well with the measured post-earthquake topography (e.g., deposit height of about 120 m and stretch length of about 900 m in the river) and landslide duration of about 1 min. The velocity of the sliding mass increases rapidly during flow, especially in the first 20 s. The velocity profiles along the depth direction at different locations of the sliding body exhibit an exponential distribution similar to that of a Bagnold-type profile, indicating that the sliding body is fully mobilized. The rate-dependent dissipation parameter β used in the model significantly influences the runout behavior (e.g., flow speed, velocity distribution, and deposit shape).

Tsunami waves were observed along the Bulgarian Black Sea coastline on 7 May 2007. The maximum rise and fall of the sea level were 1.2 and 2.0 m, respectively, with wave oscillations between 4 and 8 min. At first, submarine landsliding and then later on atmospheric disturbance were suggested as the cause of the tsunami. Numerical modeling, assuming a landslide displacing 30–60 million m3 material on the slope with a thickness range of more than 20–40 m, could reproduce the main characteristics of the recorded tsunami. In this early model, the landslide initiated on the shelf at a water depth of 100 m with a runout of approximately 20 km into 1000 m water depth. Subsequent and recent numerical modeling suggested that the failure may have initiated on the slope, anywhere between 200 and 1500 m seafloor depth. The runout of the transported sediments in these latest model was at 1850 m water depth. Just a few years after the tsunami, OMV and its joint venture partners, TOTAL and Repsol, acquired modern deepwater data sets in the same area where the submarine landsliding was assumed to occur. These data sets included multibeam swath bathymetry area and 3D reflection seismic data. These data sets offer a possibility to establish the presence of speculative submarine landslide responsible for the tsunami, with its geometry and nature. Our results provide direct evidence for the occurrence of large nonseismic catastrophic sediment failures along the Bulgarian coast. In this study, we illustrate Quaternary submarine landslides on 3D seismic reflection data immediately below the one responsible for the 2007 event; we also briefly point out the potential interpretation pitfall related to sediment waves and mass transport complexes.

Modeling texture of milled surfaces using analytic methods requires explicit knowledge of a large number of variables some of which change during machining. These include dynamically changing tool runout, deflection, workpiece material properties, displacement of the workpiece within its fixture and others. Due to the complexity of all factors combined, an alternative approach is presented utilizing the ability of neural networks and fractals to implicitly account for these combined conditions. In the initial model, predicted surface points are first connected using splines to model 3D surface maps. Results are presented over varying several cutting parameters. Then, replacing splines, an improved fractal method is presented that determines fractal characteristics of milled surfaces to model more representative surface profiles on a small scale. The fractal character of surfaces as manifested by the fractal dimension provides evidence of chaos in milling.

Geological disasters, especially landslides, frequently occur in Enshi County, Hubei Province, China. On 21 July 2020, a large-scale landslide occurred in Enshi due to continuous rainfall. The landslide mass blocked the Qingjiang River, formed a dammed lake and caused great damage to surrounding roads and village buildings. In this study, the geomechanical properties of the landslide mass were obtained through field surveys. A three-dimensional topography model of the slope was established using the particle flow code (PFC) and the numerical parameters of the model were calibrated. A 3D discrete element model (DEM) was used to simulate the propagation of Shaziba landslide, and the dynamic behavior of the landslide was divided into five stages. The simulation results show that the landslide movement lasted approximately 1000 s. The maximum average velocity of the landslide reached up to 7.5 m/s and the average runout distance was about 1000 m. The simulated morphology of the landslide deposits was in good agreement with the field data. In addition, the influence of effective modulus on the calculation results was analyzed. The results indicate that the propagation behavior of a landslide and the morphology of landslide deposits are closely related to the effective modulus in the contact model of the PFC3D.

Rockfalls present a significant hazard to human activities; therefore, their identification and knowledge about potential spatial impacts are important in planning protection measures to reduce rockfall risk. Remote sensing with unmanned aerial vehicles (UAVs) has allowed for the accurate observation of slopes that are susceptible to rockfall activity via various methods and sensors with which it is possible to digitally collect information about the rockfall activity and spatial distributions. In this work, a three-dimensional (3D) reconstruction of rock deposits (width, length, and height) and their volumes are addressed, and the results are used in a rockfall trajectory simulation. Due to the availability of different sensors on the UAV, the aim was also to observe the possible differences in the dimension estimations between photogrammetric and LiDAR (light detection and ranging) point clouds, besides the most traditional method where rock deposit dimensions are measured on the field using a measuring tape. The motivation for reconstructing rock dimensions and volumes was solely for obtaining input parameters into a rockfall model. In order to study the differences between rock-measuring methods, rock dimensions were used as input parameters in a rockfall model, and additionally, modeling results such as propagation probability, maximum kinetic energies, and maximum passing heights were compared. The results show that there are no statistically significant differences between the measurement method with respect to rock dimensions and volumes and when modeling the propagation probability and maximum passing heights. On the other hand, large differences are present with maximum kinetic energies where LiDAR point cloud measurements achieved statistically significantly different results from the other two measurements. With this approach, an automated collection and measurement process of rock deposits is possible without the need for exposure to a risk of rockfall during fieldwork.

In five-axis high speed milling of freeform surface with ball-end cutters, unwanted machining results are usually introduced by some error effects. Hence precise modeling and simulation of milled sculptured surfaces topography and roughness is the key to obtain optimal process parameters, satisfactory surface quality and high machining efficiency. In this paper, a predictive model for sculptured surface topography and roughness of ball-end milling is developed. Firstly, a mathematical model including both the relative motion of the cutter-workpiece couple and some influential factors on machined surface quality such as the tool runout, tool deflection and tool wear is proposed, and subsequently the analytical form of the tool swept envelope is derived by means of homogeneous coordinate transformation. Then the minimal z-values of the corresponding points lied in discrete cutting edges model and Z-map workpiece model are used to update the workpiece surface topography and to calculate 3D surface roughness. Finally, the simulation algorithm is realized with Matlab software. A series of machining tests on 3Cr2MoNi steel are conducted to validate the model, and the machined surface topography is found in good accordance with the simulation result.

The use of dynamic computational methods has become indispensable for addressing problems related to rockfall hazard. Although a number of models with various degrees of complexity are available, model parameters are rarely calibrated against observations from rockfall experiments. A major difficulty lies in reproducing the apparent randomness of the impact process related to both ground and block irregularities. Calibration of rigorous methods capable of explicitly modeling trajectories and impact physics of irregular blocks is difficult, as parameter spaces become too vast and the quality of model input and observation data are insufficient. The model presented here returns to the simple “lumped-mass” approach and simulates the characteristic randomness of rockfall impact as a stochastic process. Despite similarities to existing approaches, the model presented here incorporates several novel concepts: (i) ground roughness and particle roughness are represented as a random change of slope angle at impact; (ii) lateral deviations of rebound direction from the trajectory plane at impact are similarly accounted for by perturbing the ground orientation laterally, thus inducing scatter of run-out directions; and (iii) a hyperbolic relationship connects restitution factors to impact deformation energy. With these features, the model is capable of realistically accounting for the influence of particle mass on dynamic behaviour. The model only requires four input parameters, rendering it flexible for calibration against observed datasets. In this study, we calibrate the model against observations from the rockfall test site at Vaujany in France. The model is able to reproduce observed distributions of velocity, jump heights, and runout at observation points. In addition, the spatial distribution of the trajectories and landing points has been successfully simulated. Different parameter sets have been used for different ground materials such as an avalanche channel, a forest road, and a talus cone. Further calibration of the new model against a range of field datasets is essential. This study is part of an extensive calibration program that is still in progress at this first presentation of the method, and focuses on fine-tuning the details of the stochastic process implemented both in two-dimensional (2D) and three-dimensional (3D) versions of the model.

Abstract The modeling of soil erosion processes is affected by several factors that reflect the physical-geographic conditions of the study site together with the land use linkage. The soil parameters are significant in the modeling of erosion and also runoff processes. The correct determination of a soil's parameters becomes a crucial part of the model's calibration. This paper deals with a sensitivity analysis of seven soil input parameters to the physically-based Erosion 3D model. The results show the variable influence of each soil parameter. The Erosion 3D model is very sensitive to initial soil moisture, bulk density, and erodibility.

A highly modular and scale-consistent Terrestrial Systems Modeling Platform (TerrSysMP) is presented. The modeling platform consists of an atmospheric model (Consortium for Small-Scale Modeling; COSMO), a land surface model (the NCAR Community Land Model, version 3.5; CLM3.5), and a 3D variably saturated groundwater flow model (ParFlow). An external coupler (Ocean Atmosphere Sea Ice Soil, version 3.0; OASIS3) with multiple executable approaches is employed to couple the three independently developed component models, which intrinsically allows for a separation of temporal–spatial modeling scales and the coupling frequencies between the component models. Idealized TerrSysMP simulations are presented, which focus on the interaction of key hydrologic processes, like runoff production (excess rainfall and saturation) at different hydrological modeling scales and the drawdown of the water table through groundwater pumping, with processes in the atmospheric boundary layer. The results show a strong linkage between integrated surface–groundwater dynamics, biogeophysical processes, and boundary layer evolution. The use of the mosaic approach for the hydrological component model (to resolve subgrid-scale topography) impacts simulated runoff production, soil moisture redistribution, and boundary layer evolution, which demonstrates the importance of hydrological modeling scales and thus the advantages of the coupling approach used in this study. Real data simulations were carried out with TerrSysMP over the Rur catchment in Germany. The inclusion of the integrated surface–groundwater flow model results in systematic patterns in the root zone soil moisture, which influence exchange flux distributions and the ensuing atmospheric boundary layer development. In a first comparison to observations, the 3D model compared to the 1D model shows slightly improved predictions of surface fluxes and a strong sensitivity to the initial soil moisture content.

Abstract. Rockfall presents an ongoing challenge to the safe operation of transportation infrastructure, creating hazardous conditions which can result in damage to roads and railways, as well as loss of life. Rockfall risk assessment frameworks often involve the determination of rockfall runout in an attempt to understand the likelihood that rockfall debris will reach an element at risk. Rockfall modelling programs which simulate the trajectory of rockfall material are one method commonly used to assess potential runout. This study aims to demonstrate the effectiveness of a rockfall simulation prototype which uses the Unity 3D game engine. The technique is capable of simulating rockfall events comprised of many mobile fragments, a limitation of many industry standard rockfall modelling programs. Five fragmental rockfalls were simulated using the technique, with slope and rockfall geometries constructed from high-resolution terrestrial laser scans. Simulated change detection was produced for each of the events and compared to the actual change detection results for each rockfall as a basis for testing model performance. In each case the simulated change detection results aligned well with the actual observed change in terms of location and magnitude. An example of how the technique could be used to support the design of rockfall catchment ditches is shown. Suggestions are made for future development of the simulation technique with a focus on better informing simulated rockfall fragment size and the timing of fragmentation.

Gli eventi sismici sono riconosciuti come una delle maggiori cause per l’innesco di frane (Keefer, 1984). Le frane sismo-indotte sono documentate sin dal IV secolo (Seed, 1968). È stata condotta un’analisi sulla distribuzione spaziale delle frane sismo-indotte nell’area circostante la sorgente sismogenetica per meglio comprendere il loro innesco in aree sismiche e per delimitare la massima distanza alla quale un sisma con data magnitudo possa indurre frane. Tuttavia, quando si applicano tali approcci a eventi storici si pone un problema legato al sottocampionamento delle frane più piccole, che possono essere obliterate dall'erosione e dall'evoluzione del paesaggio. Per questo motivo è importante caratterizzare accuratamente la distribuzione delle frane, in termini di dimensione, in funzione della distanza dalla sorgente sismica. Sono stati analizzati sei terremoti in tutto il mondo che hanno innescato un significativo numero di frane (Finisterre 1993, Northridge 1994, Niigata 2004, Wenchuan 2008, Iwate 2008 and Tohoku 2011) per meglio comprendere le relazioni esistenti tra la distribuzione spaziale delle frane, l’accelerazione di picco al suolo (PGA), la distanza dalla sorgente, il relief e le litologie presenti nell’area. Si è osservata una forte relazione tra la PGA e la dimensione delle frane, mentre una la relazione tra la loro dimensione e la distanza dalla sorgente non è altrettanto chiara, ciò legato all’interazione tra diversi fattori quali ad esempio il relief e la litologia. Sono state realizzate e analizzate le curve magnitudo-frequenza (MFC) per differenti distanze dall’area sorgente attraverso varie metodologie: stimatore di massima verosimiglianza per distribuzioni di tipo potenza cumulate (Clauset et al, 2009), stimatore di massima verosimiglianza per distribuzioni di tipo potenza non cumulate, regressione ai minimi quadrati per funzioni di tipo potenza non cumulate in scala logaritmica e stimatore di massima verosimiglianza per la distribuzione Double Pareto. Dalle analisi si è potuto osservare un decrescere della densità spaziale delle frane con la distanza, ma un basso impatto della dimensione delle frane. Inoltre la funzione Double Pareto è stata scelta come miglior strumento per il fittaggio dei dati (Valagussa et al, 2014). Allo scopo di definire il rischio legato alle frane sismo-indotte è stata sviluppata una metodologia per la zonazione probabilistica quantitativa del rischio da frane da crollo (Valagussa et al, 2014). Il metodo è stato applicato e dimostrato nell’area del Friuli (Apli orientali) colpita da un terremoto di magnitudo 6.4 nel 1976. Quattro inventari sono stati realizzati sia tramite attività di terreno che da dati storici. La metodologia si basa sul vettore di rischio tridimensionale (RHVmod) le cui componenti includo l’energia cinetica, l’altezza di volo e la frequenza annua. I primi due valori sono calcolati per ogni cella del versante per mezzo del programma Hy-STONE. La frequenza annua è invece determinata moltiplicando la frequenza d’innesco annua per il numero di transiti simulati in ogni cella. La frequenza d’innesco annua è calcolata combinando l’area instabile, calcolata per 10 differenti scenari con differente frequenza annua di occorrenza sulla base di caratteristiche morfometriche e sismiche, e la curva magnitudo-frequenza relativa dei blocchi identificati da attività di terreno. Una serie di analisi discriminanti sono state condotte per determinare le variabili che controllano l’area in frana, sulla base degli inventari redatti e di DEMs a differenti risoluzioni (1 e 10m). L’analisi ha dimostrato il ruolo rilevante della curvatura nella definizione dell’area instabile. Per verificare la validità della mappa di PGA utilizzata nelle analisi, una nuova mappa è stata redatta sulla base delle Precarious Balanced Rocks identificate sul terreno.
Earthquakes have been recognized as a major cause of landsliding (Keefer, 1984), and landslides triggered by earthquakes have been documented since the IV century (Seed, 1968). The spatial distribution of earthquake-induced landslides around the seismogenetic source has been analysed to better understand the triggering of landslides in seismic areas and to forecast the maximum distance at which an earthquake, with a certain magnitude, can trigger landslides. However, when applying such approaches to old earthquakes one should be concerned about the undersampling of smaller landslides, which can be cancelled, by erosion and landscape evolution. For this reason, it is important to characterize carefully the size distribution of landslides as a function of distance from the earthquake source. I analysed six earthquakes in the world that triggered significant amount of landslides (Finisterre 1993, Northridge 1994, Niigata 2004, Wenchuan 2008, Iwate 2008 and Tohoku 2011) to better understand the relation between the spatial distribution of the landslides, the peak ground acceleration (PGA), the distance from the sources, the relief and the lithologies of the area. I observed a strong relationship between landslides size and PGA, while the relationship between the distance from the source and the landslide size distribution is not clear, due to the interaction of different factors such as relief and lithology. I also developed magnitude frequency curves (MFC) for different distances from the source area by using different methods, such as: the maximum likelihood estimator of cumulative power-law distribution (Clauset et al, 2009); the maximum likelihood estimator of non-cumulative power-law function; the least square regression of non-cumulative log power-law function and the maximum likelihood estimator of Double Pareto distribution. I observed a decrease of the spatial density of landslides with distance, with a small effect of the size of these landslides. I also identify the Double Pareto function as the best tool for the fitting of the data (Valagussa et al., 2014a). In order to define the hazard due to earthquake-induced landslides, I developed a methodology for quantitative probabilistic hazard zonation for rockfalls (Valagussa et al., 2014b). I applied and demonstrated the method in the area of Friuli (Eastern Italian Alps) that was affected by the 1976 Mw 6.5 earthquake. Four rockfall datasets have been prepared from both historical data and field surveys. The methodology relies on a three-dimensional hazard vector (RHVmod), whose components include the rockfall kinetic energy, the fly height, and the annual frequency. The values of the first two components are calculated for each location along the slope using the 3D rockfall runout simulator Hy-STONE. The rockfall annual frequency is assessed by multiplying the annual onset frequency by the simulated transit frequency. The annual onset frequency is calculated 2 through a procedure that combines the extent of unstable areas, calculated for 10 different seismichazard scenarios with different annual frequencies of occurrence, and the magnitude relativefrequency relationship of blocks as derived from the collected field data. For each annual frequency of occurrence, the unstable area is calculated as a function of morphometric and earthquake characteristics. A series of discriminant-analysis models, using the rockfall datasets and DEMs of different resolution (1 and 10 m), identified the controlling variables and verified the model robustness. In contrast with previously published research, I show that the slope curvature plays a relevant role in the computation of the unstable area. To ensure the validity of the peak ground acceleration used as seismic parameter in the discriminant function, I also try to define a map of PGA based on the precarious balanced rocks surveyed on the field.

Modelling hydrologic processes is important for understanding how the water cycle works in different environments. Cities which undergo constant changes are subject to flood hazards resulting from severe rainfall. This paper aims to simulate severe rainfall, visualize the results, incorporating both spatial and temporal dimensions, and to make future recommendations for further studies on flood visualization. Visualizing the results from a rainfall simulation using GIS provides urban planners and others the means to view the dynamics of the surface runoff. At the same time, it makes accessible advanced querying and analytical tools. A hydrological model for the study area in Gävle, Sweden was used to simulate a 100-year rainfall. Through FME, the data was reduced, time-stamped and combined to a shapefile. Both 2D software, ArcGIS, and 3D software, ArcScene, were used for creating an animated flood visualization. This study shows that although 2D tested better by a group of planners and water professionals, the 3D was still considered more intuitive. The heightened sense of realism from 3D outweighs its drawbacks, and further studies are required to test different methods of 3D visualization.

The aim of this research project is to improve the ability to forecast the propagation of shallow and rapid landslides through the combination of remote sensing and advanced methods of numerical modelling techniques. The activities carried out during the PhD have been done in order to achieve three objectives: - to define the procedures and methods of numerical analysis of the propagation of shallow rapid landslides, through the use of existing codes that allow to model the post- failure behaviour dynamic; - to use remote sensing data derived from LiDAR images (aerial or terrestrial) or radar monitoring data (in particular GB-InSAR), to obtain information about the geometric and geomorphological features of the landslides and to find new future sources areas that can be used with numerical models in order to obtain new risk scenarios; - to use abovementioned the new procedure on case studies of high-risk landslide. Two numerical modelling DAN-W and DAN-3D (Hungr, 1995 McDougall, 2006) are used, based on Lagrangian numerical method to solve the equations of St. Venant. Consider the mass as an "equivalent fluid", governed by an internal frictional rheology and basal rheology that must be chosen with back analysis by trial-and-error procedure. These parameters are then used to forecast analysis using new source areas with new volumes obtained by LiDAR and GB-InSAR monitoring data analysis. The latter, in particular, are processed through a numerical code developed in the MATLAB language. The MATLAB code uses the cumulative displacement maps relative to a selected time interval, to calculate possible sources areas according to a threshold and finally analysed the frequency of movement occurrence. This procedure has been applied at four case studies: gravity induced-Pyroclastic Density Currents (gi-PDC) that take place on the side of Stromboli volcanic island; Mount Rotolon debris flow; Gessi-Mazzalasino debris slide and Santa Trada rock/debris slide.

Carbide re-hobbing is a variation of the gear hobbing process. It is typically used for finishing fully-hardened gear blanks that have been semi-finished, generally by a previous hobbing operation. This paper will discuss a new approach to modeling the carbide re-hobbing process with the goal of improving part quality for a typical pinion. Prior modeling approaches have been based on analytical chip calculation methods. Such approaches, however, limit the geometry of the tool and candidate workpiece to such profiles as would be implemented in the model initially. This new modeling approach involves the use of CAD/CAM/CAE tools to simulate the hobbing process in a virtual 3D environment. As such, the models may now take into account the specific tool geometries, workpiece geometries, setup errors and various cutting conditions with much greater ease. The results of the simulation in predicting cutting forces, part deflection and the resulting profile deviations will be presented. Further, the effect of tool setup error, in particular both synchronous and asynchronous runout, on part quality will be examined in simulation. The simulation results reveal that each type of runout provides a unique signature of profile deviation error for the left and right flanks. The relationship between these setup errors and resulting profile errors will be examined in detail and compared with data from controlled machining tests.

This paper introduces an approach for modeling and representation of geometric tolerances on any 3D solid model using the Objected Oriented Programming (OOP) paradigm. The modeling scheme is supported by a comprehensive validation engine, which certifies the tolerance type against the 3D geometry context both syntactically and semantically. The major objective of this work is to develop a methodology for interfacing tolerance modeling with boundary representation (B-Rep) based 3D solid model geometry. We will demonstrate that the OOP paradigm is very efficient and flexible for tolerance model representation, which is required within the interactive design process. Six categories of tolerance classes have been developed for size, form, orientation, position, runout and profile, which extend a general tolerance class through inheritance. An instance of the general tolerance class will be initialized when picking a feature or a group of features to tolerance, depending upon feature(s’) characteristics and attributes. To apply a tolerance object the system obtains the 3D geometric data from the solid model using the feature extraction paradigm. When the required tolerance type is selected for modeling, an instance from the specified tolerance type class will be initialized through inheritance from the general feature tolerance class and gathers the necessary information / tolerance data. An intelligent validation engine that supports the modeler is introduced. The engine validates any selected tolerancing activity in two stages. First, it ensures that the selected feature or group of features is suitable for the selected tolerance type. Second, it ensures that the data specified does not lead to over/under-dimensioning. The paper also discusses a prototype system implemented to test the modeler and the validation engine. The results have been very encouraging while testing the system on a number of engineering models.

Abstract Modeling texture of milled surfaces using analytic methods requires explicit knowledge of a large number of variables some of which change during machining. These include dynamically changing tool runout, deflection, work-piece material properties, displacement of the workpiece within its fixture and others. Due to the complexity of all factors combined, an alternative approach is presented utilizing the ability of neural networks and fractals to implicitly account for these combined conditions. In the initial model, predicted surface points are first connected using splines to reconstruct 3D surface maps. Results are presented over varying several cutting parameters. Then, replacing splines, an improved fractal method is presented that determines fractal characteristics of milled surfaces to reconstruct more representative surface maps on a small scale. The fractal character of self-similarity within surfaces as manifested by the fractal dimension provides evidence of chaos in milling.