Počet záznamů: 1
Efficient and flexible MATLAB implementation of 2D and 3D elastoplastic problems
- 1.
SYSNO ASEP 0504439 Druh ASEP J - Článek v odborném periodiku Zařazení RIV J - Článek v odborném periodiku Poddruh J Článek ve WOS Název Efficient and flexible MATLAB implementation of 2D and 3D elastoplastic problems Tvůrce(i) Čermák, Martin (UGN-S)
Sysala, Stanislav (UGN-S) RID, ORCID
Valdman, Jan (UTIA-B) RID, ORCIDCelkový počet autorů 3 Zdroj.dok. Applied Mathematics and Computation. - : Elsevier - ISSN 0096-3003
Roč. 355, August 2019 (2019), s. 595-614Poč.str. 20 s. Forma vydání Online - E Jazyk dok. eng - angličtina Země vyd. US - Spojené státy americké Klíč. slova MATLAB code vectorization ; elastoplasticity ; finite element method ; tangential stiffness matrix ; semismooth Newton method Vědní obor RIV BA - Obecná matematika Obor OECD Applied mathematics Vědní obor RIV – spolupráce Ústav teorie informace a automatizace - Obecná matematika CEP LQ1602 GA MŠMT - Ministerstvo školství, mládeže a tělovýchovy Způsob publikování Omezený přístup Institucionální podpora UGN-S - RVO:68145535 ; UTIA-B - RVO:67985556 UT WOS 000464930500044 EID SCOPUS 85063371842 DOI 10.1016/j.amc.2019.02.054 Anotace Fully vectorized MATLAB implementation of various elastoplastic problems formulated in terms of displacement is considered. It is based on implicit time discretization, the finite element method and the semismooth Newton method. Each Newton iteration represents a linear system of equations with a tangent stiffness matrix. We propose a decomposition of this matrix consisting of three large sparse matrices representing the elastic stiffness operator, the strain-displacement operator, and the derivative of the stress-strain operator. The first two matrices are fixed and assembled once and only the third matrix needs to be updated in each iteration. Assembly times of the tangent stiffness matrices are linearly proportional to the number of plastic integration points in practical computations and never exceed the assembly time of the elastic stiffness matrix. MATLAB codes are available for download and provide complete finite element implementations in both 2D and 3D assuming von Mises and Drucker–Prager yield criteria. One can also choose several finite elements and numerical quadrature rules. Pracoviště Ústav geoniky Kontakt Lucie Gurková, lucie.gurkova@ugn.cas.cz, Tel.: 596 979 354 Rok sběru 2020 Elektronická adresa https://www.sciencedirect.com/science/article/pii/S0096300319301584
Počet záznamů: 1