Development of semi-implicit midpoint and Romberg stress integration algorithms for single hardening soil constitutive models
Autor: | Divyanshu Kumar Lal, Arghya Das |
---|---|
Rok vydání: | 2020 |
Předmět: |
Computer science
Numerical analysis Constitutive equation 0211 other engineering and technologies General Engineering 02 engineering and technology 01 natural sciences Midpoint Finite element method Computer Science Applications 010101 applied mathematics Computational Theory and Mathematics Romberg's method Projection method Boundary value problem 0101 mathematics Algorithm Software Cutting-plane method 021101 geological & geomatics engineering |
Zdroj: | Engineering Computations. 37:3477-3503 |
ISSN: | 0264-4401 |
DOI: | 10.1108/ec-08-2019-0358 |
Popis: | Purpose Semi-implicit type cutting plane method (CPM) and fully implicit type closest point projection method (CPPM) are the two most widely used frameworks for numerical stress integration. CPM is simple, easy to implement and accurate up to first order. CPPM is unconditionally stable and accurate up to second order though the formulation is complex. Therefore, this study aims to develop a less complex and accurate stress integration method for complex constitutive models. Design/methodology/approach Two integration techniques are formulated using the midpoint and Romberg method by modifying CPM. The algorithms are implemented for three different classes of soil constitutive model. The efficiency of the algorithms is judged via stress point analysis and solving a boundary value problem. Findings Stress point analysis indicates that the proposed algorithms are stable even with a large step size. In addition, numerical analysis for solving boundary value problem demonstrates a significant reduction in central processing unit (CPU) time with the use of the semi-implicit-type midpoint algorithm. Originality/value Traditionally, midpoint and Romberg algorithms are formulated from explicit integration techniques, whereas the present study uses a semi-implicit approach to enhance stability. In addition, the proposed stress integration algorithms provide an efficient means to solve boundary value problems pertaining to geotechnical engineering. |
Databáze: | OpenAIRE |
Externí odkaz: |