 |
| Titel |
3-D FEM modelling of geological structures caused by geometrical instabilities and contrasts in rock strength |
| VerfasserIn |
Marina von Tscharner, Stefan Schmalholz |
| Konferenz |
EGU General Assembly 2014
|
| Medientyp |
Artikel
|
| Sprache |
Englisch
|
| Digitales Dokument |
PDF |
| Erschienen |
In: GRA - Volume 16 (2014) |
| Datensatznummer |
250090665
|
| Publikation (Nr.) |
 EGU/EGU2014-4918.pdf |
|
|
|
|
|
| Zusammenfassung |
| Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth’s crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These structures occur from the centimeter to the kilometer scale and the related deformation processes control the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D deformation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. We study the 3-D evolution of geometrical instabilities with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FEM code combines a numerical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FEM mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D deformation for a visco-elasto-plastic rheology. Stresses are limited by a yield stress using a visco-plastic formulation and the viscous rheology is described by a power-law flow law. The 3-D FEM code is applied to model 3-D power-law folding and power-law Rayleigh-Taylor instabilities (diapirs) with different re-meshing scenarios. The results are tested with the analytical solution for small amplitudes and with 2-D numerical results for large amplitudes. Thereby,the small initial geometrical perturbations for folding and necking are exactly followed by the FEM mesh. In order to test and measure the numerical properties for an Eulerian mesh we use the analytical solution for a two-dimensional viscous inclusion in pure shear. The FEM code is further applied to study the evolution of lithospheric-scale 3-D slab detachment and its lateral propagation. |
|
|
|
|
|
|