360 likes | 587 Vues
Influence of buried craters on wrinkle ridges on Mars. Erwan Mazarico 12.521 Final Project. MOLA. Outline. Motivation Wrinkle Ridges on Mars Observations of Buried craters Modeling with Abaqus Preliminary Results Conclusion and Future Work. Motivation. Wrinkle ridges
E N D
Influence of buried craters on wrinkleridges on Mars Erwan Mazarico 12.521 Final Project MOLA
Outline • Motivation • Wrinkle Ridges on Mars • Observations of Buried craters • Modeling with Abaqus • Preliminary Results • Conclusion and Future Work
Motivation • Wrinkle ridges • Globally observed • Compressional features • On volcanic plains, some “ghost craters” are observed Hypothesis = buried craters that influences the distribution of stresses at the surface of the volcanic layer
Wrinkle Ridges • Topographic observations • Maximum relief of several 100’s meters • Spacing ~20-35km Solis Planum Lunae Planum Watters(2004)
Wrinkle formation • Several scenarios • Single blind listric thrust fault • Thermal stresses • Multiple folds/faults (thrust, back-thrust) Montési and Zuber (2003)
Single blind listric thrust fault Coulomb Watters (2004)
Single blind listric thrust fault • Listric fault does not break the surface • Quite sensitive to fault geometry and slip • Thin-skinned formation • Consistent with buried craters localizing the thrust faults Coulomb
Thermal Stresses • Thermal stresses due to global cooling are large • Compressional mechanism to form wrinkle ridges far from obvious loading/uplift zones • Shallow expression because elastic tensional stresses at greater depth • Explains ‘globality’
Deep thrust faults • Large-scale topographic gradient on top of ridges • Large, deep (4-15km) thrust faults(<30°) • Near surface: fault stops; spay into smaller forward- or back-thrusts; accommodated by folding Golombek et al. (2001) Montési and Zuber (2003) Montési and Zuber (2003)
Ghost Crater • Formation Hypothesis Cartoon initial length
Ghost Crater • Formation Hypothesis Cartoon initial length initial length
Modeling with Abaqus • High-end FEM code with GUI, automatic meshing, etc. • Constructing the model and putting the BC • Different geometries tested • Crater modeled as cylinder, smooth cylinder, part of a sphere • Different displacements applied • Pushing on the X+ face, or both X-/X+ faces
Using Abaqus Model Material • = 0.6 • = 0.25 • E = 20GPa • = 2700 kg.m-3 176km x 176km x 25km
Using Abaqus Mesh Tetrahedral elements
Using Abaqus Boundary Conditions
Using Abaqus Gravity Load
Using Abaqus Impose displacement
Abaqus Rheology Test Xtotal=19km=12.5% Elastic ( = 0.25 ) Elasto-Plastic Associated V≈9% V≈30.5% !! V≈2% V≈0% Elasto-Plastic Non-Associated Elastic Incompressible( = 0.5 )
Abaqus Rheology Test • Elastic is definitely not the rheology to use in our case • Elasto-Plastic Associated (EPA) causes some concern from the test … • Elasto-Plastic Non-Associated (EPNA) often crashes due to non-uniqueness (negative eigenvalues), or too high local plastic strains ►most runs are done in EPA
Abaqus Results • Cylindrical crater • Cylindrical crater with smoothed edges • Spherical crater
Cylindrical crater von Mises stresses
Cylindrical crater Vertical displacement
Cylindrical Crater Vertical displacement
Cylindrical Crater Vertical displacement
Cylindrical Crater AC Yield
Smoothed Cylindrical Crater Vertical displacement
Smoothed Cylindrical Crater Plastic Strain
Spherical Crater Vertical Displacement
Spherical Crater Total Strain
Spherical Crater Concern AC Yield
Problems encountered • Strange behavior of Elasto-Plasticity Associated material … (V, AC Yield) • Could not produce real wrinkles • Volcanic plain not superimposed (attempted but not successful)
(very) preliminary conclusions • Localization of uplift at the crater rims seems to happen but needs futher work to be conclusive, especially after adding a layer above • Hard to form wrinkles and to use Abaqus properly !
Future Work • Model crater with top layer (2-8km thick) with appropriate contact interaction • Use Abaqus output stresses to produce model of faults • Use 3d-def to get surface displacements • See conditions to have visible wrinkles following the buried craters • Count ghost craters from MOLA and see if consistent