0504439 20240111141019.020190507235959.9 85063371842 000464930500044 10.1016/j.amc.2019.02.054 20 s. E cav_un_epca*02561600096-3003355595-614Elsevier MATLAB code vectorization elastoplasticity finite element method tangential stiffness matrix semismooth Newton method cav_un_auth*0296439 Čermák Martin UGN-S Oddělení aplikované matematiky a informatiky & Oddělení IT4Innovations Department of applied mathematics and computer science and Department IT4Innovations Applied Mathematics and Computer Science & IT4Innovations CZ K Ústav geoniky AV ČR, v. v. i. cav_un_auth*0221817 Sysala Stanislav UGN-S Oddělení aplikované matematiky a informatiky & Oddělení IT4Innovations Department of applied mathematics and computer science and Department IT4Innovations Applied Mathematics and Computer Science & IT4Innovations CZ Ústav geoniky AV ČR, v. v. i. cav_un_auth*0292941 Valdman Jan UTIA-B Matematická teorie rozhodování Department of Decision Making Theory MTR MTR Department of Decision Making Theory Ústav teorie informace a automatizace AV ČR, v. v. i. textový soubor https://www.sciencedirect.com/science/article/pii/S0096300319301584 cav_un_auth*0333099 LQ1602 GA MŠk 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. WOS BA 10000 10100 10102 2020 3 UTIA-B 10000 10100 10102 4 X hod 4xh 20231122143959.3 UTIA-B BA RVO:68145535 RVO:67985556 http://hdl.handle.net/11104/0296072 1 Applied Mathematics and Computer Science & IT4Innovations UGN-S 10102 MATHEMATICS, APPLIED S 3+4 Article Chemistry Physical|Physics Atomic Molecular Chemical C MATHEMATICS.APPLIED 2.709 1.442 6.4 0.04193 0.667 30975 182 0.969 34.25 3.79 7682 Q1 3.102 862 Q1 49.057 97.5 Q3 Q4 2.4 Q1 97.51 2019 Soubory v repozitáři: UGN_0504439.pdf 85063371842 SCOPUS PUBMED 000464930500044 WOS textový soubor cav_un_epca*0256160 Applied Mathematics and Computation 0096-3003 1873-5649 Roč. 355 August 2019 2019 595 614 Elsevier