Volume 24 Issue 6
Dec.  2018
Turn off MathJax
Article Contents
YANG Shaohua, LI Zhonghai, 2018. A NUMERICAL CALCULATION APPROACH BASED ON FEM FOR LONG-TERM DEFORMATION OF LITHOSPHERE. Journal of Geomechanics, 24 (6): 768-775. DOI: 10.12090/j.issn.1006-6616.2018.24.06.079
Citation: YANG Shaohua, LI Zhonghai, 2018. A NUMERICAL CALCULATION APPROACH BASED ON FEM FOR LONG-TERM DEFORMATION OF LITHOSPHERE. Journal of Geomechanics, 24 (6): 768-775. DOI: 10.12090/j.issn.1006-6616.2018.24.06.079

A NUMERICAL CALCULATION APPROACH BASED ON FEM FOR LONG-TERM DEFORMATION OF LITHOSPHERE

doi: 10.12090/j.issn.1006-6616.2018.24.06.079
More Information
  • Received: 2018-08-13
  • Revised: 2018-08-31
  • Published: 2018-12-01
  • Finite element method (FEM), which is very important in almost all fields of solid earth sciences, has been widely used in numerical experiments in solid earth sciences due to its flexibility and precision, ranging from short-term seismic activity to long-term lithospheric deformation, mantle convection and even planetary evolution. However, with the deepening of the research, some specific geological problems bring challenges to the finite element calculation, especially the numerical calculation of large-scale lithosphere deformation, such as the evolution of subduction zone and stress localization caused by plastic rheology in shear zone. Based on explicit FEM, an attempt is made to numerical calculation in the process of large deformation of lithosphere considering visco-elastic-plasticity in this study. Marker-In-Cell (MIC) approach was used in processing each material migration. On the basis of describing the basic principle and process, the core modules of visco-elastic deformation, elastic-plastic deformation, large deformation and heat transfer were tested. These four modules are the key to simulate the long-term deformation of the lithosphere. According to the comparison between the test results and the previous simulation results, the tested core modules basically meet the test requirements. It can be predicted that the existing basic algorithms can meet the needs of studying the large deformation of the lithosphere, and further specific research work will explore this kind of problem.

     

  • Full-text Translaiton by iFLYTEK

    The full translation of the current issue may be delayed. If you encounter a 404 page, please try again later.
  • loading
  • [1]
    Hughes T J R. The finite element method:linear static and dynamic finite element analysis[M]. New York:Dover Publication, 2000.
    [2]
    Ismail-Zadeh A, Tackley P. Computational methods for geodynamics[M]. Cambridge:Cambridge University Press, 2010.
    [3]
    Tucker G E, Hancock G R. Modelling landscape evolution[J]. Earth Surface Processes and Landforms, 2010, 35(1):28~50. doi: 10.1002/esp.v35:1
    [4]
    Yang S H, Shi Y L. Three-dimensional numerical simulation of glacial trough forming process[J]. Science China Earth Sciences, 2015, 58(9):1656~1668. doi: 10.1007/s11430-015-5120-8
    [5]
    Beaumont C, Jamieson R A, Nguyen M H, et al. Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation[J]. Nature, 2001, 414(6865):738~742. doi: 10.1038/414738a
    [6]
    Van Keken P E, Kiefer B, Peacock S M. High-resolution models of subduction zones:implications for mineral dehydration reactions and the transport of water into the deep mantle[J]. Geochemistry, Geophysics, Geosystems, 2002, 3(10):1~20. doi: 10.1029/2001GC000256/pdf
    [7]
    Leng W, Zhong S J. Viscous heating, adiabatic heating and energetic consistency in compressible mantle convection[J]. Geophysical Journal International, 2008, 173(2):693~702. doi: 10.1111/gji.2008.173.issue-2
    [8]
    杨少华, 许志琴, 李忠海, 等.欧罗巴星陨石坑对冰层厚度的制约[J].地球物理学报, 2017, 60(3):935~940. http://www.cnki.com.cn/Article/CJFDTotal-DQWX201703009.htm

    YANG Shaohua, XU Zhiqin, LI Zhonghai, et al. Constraint of impact craters on ice thickness on the Europa[J]. Chinese Journal of Geophysics, 2017, 60(3):935~940(in Chinese with English abstract) http://www.cnki.com.cn/Article/CJFDTotal-DQWX201703009.htm
    [9]
    Zhang S Q, O'Neill C. The early geodynamic evolution of Mars-type planets[J]. Icarus, 2016, 265:187~208. doi: 10.1016/j.icarus.2015.10.019
    [10]
    He Y, Puckett E G, Billen M I. A discontinuous Galerkin method with a bound preserving limiter for the advection of non-diffusive fields in solid Earth geodynamics[J]. Physics of the Earth and Planetary Interiors, 2017, 263:23~37. doi: 10.1016/j.pepi.2016.12.001
    [11]
    Enright D, Losasso F, Fedkiw R. A fast and accurate semi-Lagrangian particle level set method[J]. Computers & Structures, 2005, 83(6~7):479~490. http://www.ams.org/mathscinet-getitem?mr=2143508
    [12]
    Weinberg R F, Schmeling H. Polydiapirs:multiwavelength gravity structures[J]. Journal of Structural Geology, 1992, 14(4):425~436. doi: 10.1016/0191-8141(92)90103-4
    [13]
    Gerya T V. Introduction to numerical geodynamic modelling[M]. Cambridge, UK:Cambridge University Press, 2010.
    [14]
    Ranalli G. Rheology of the Earth[M]. UK:Chapman and Hall, 1995.
    [15]
    Kronbichler M, Heister T, Bangerth W. High accuracy mantle convection simulation through modern numerical methods[J]. Geophysical Journal International, 2012, 191(1):12~29. doi: 10.1111/gji.2012.191.issue-1
    [16]
    Kaus B J P. Factors that control the angle of shear bands in geodynamic numerical models of brittle deformation[J]. Tectonophysics, 2010, 484(1~4):36~47. doi: 10.1016/j.tecto.2009.08.042
    [17]
    Moresi L, Dufour F, Mühlhaus H B. Mantle convection modeling with viscoelastic/brittle lithosphere:numerical methodology and plate tectonic modeling[J]. Pure and Applied Geophysics, 2002, 159(10):2335~2356. doi: 10.1007/s00024-002-8738-3
    [18]
    Kaus B J P, Popov A A, Baumann T S, et al. Forward and inverse modelling of lithospheric deformation on geological timescales[A]. NIC Symposium 2016[C]. NIC Series, 2016, 298~307.
    [19]
    Gerya T V, Yuen D A. Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems[J]. Physics of the Earth and Planetary Interiors, 2007, 163(1~4):83~105. doi: 10.1016/j.pepi.2007.04.015
    [20]
    Jaeger J C, Cook N G W, Zimmerman R W. Fundamentals of rock mechanics[M]. Malden:Blackwell Publishingl, 2007.
    [21]
    王勖成.有限单元法[M].北京:清华大学出版社, 2003.

    WANG Xucheng. Finite element method[M]. Beijng:Tsinghua University Press, 2003. (in Chinese)
    [22]
    夏志皋.塑性力学[M].上海:同济大学出版社, 1991.

    XIA Zhigao, Plasticity mechnics[M]. Shanghai:Tongji University Press, 1991. (in Chinese)
    [23]
    Shewchuk J R. Delaunay refinement algorithms for triangular mesh generation[J]. Computational Geometry, 2002, 22(1~3):21~74. doi: 10.1016/S0925-7721(01)00047-5
    [24]
    Burden R L, Faires J D. Numerical analysis[M]. Boston:Brooks/Cole, 2011.
    [25]
    李庆扬, 王能超, 易大义.数值分析[M]. 5版.北京:清华大学出版社, 2008.

    LI Qingyang, WANG Nengchao, YI Dayi. Numerical analysis[M]. 5th ed. Beijing:Tsinghua University Press, 2008. (in Chinese)
    [26]
    Van Keken P E, King S D, Schmeling H, et al. A comparison of methods for the modeling of thermochemical convection[J]. Journal of Geophysical Research:Solid Earth, 1997, 102(B10):22477~22495. doi: 10.1029/97JB01353
    [27]
    Thieulot C. ELEFANT:a user-friendly multipurpose geodynamics code[J]. Solid Earth Discussion, 2014, 6(2):1949~2096. doi: 10.5194/sed-6-1949-2014
    [28]
    Turcotte D L, Schubert G. Geodynamics[M]. Cambridge, UK:Cambridge University Press, 2002.
    [29]
    le Pourhiet L, May D A, Huile L, et al. A genetic link between transform and hyper-extended margins[J]. Earth and Planetary Science Letters, 2017, 465:184~192. doi: 10.1016/j.epsl.2017.02.043
    [30]
    Ruh J B, Gerya T, Burg J P. 3D effects of strain vs. velocity weakening on deformation patterns in accretionary wedges[J]. Tectonophysics, 2014, 615~616:122~141. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=a43bd9c299bbc09c6d683626a1ea36c3
  • 加载中

Catalog

    Figures(5)

    Article Metrics

    Article views (236) PDF downloads(21) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return