Efficient and flexible MATLAB implementation of 2D and 3D elastoplastic problems

Kód výsledku v IS VaVaI

<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F68145535%3A_____%2F19%3A00504439" target="_blank" >RIV/68145535:_____/19:00504439 - isvavai.cz</a>

Nalezeny alternativní kódy

RIV/67985556:_____/19:00504439 RIV/60076658:12310/19:43899348 RIV/61989100:27120/19:10241990 RIV/61989100:27730/19:10241990

Výsledek na webu

<a href="https://www.sciencedirect.com/science/article/pii/S0096300319301584" target="_blank" >https://www.sciencedirect.com/science/article/pii/S0096300319301584</a>

DOI - Digital Object Identifier

<a href="http://dx.doi.org/10.1016/j.amc.2019.02.054" target="_blank" >10.1016/j.amc.2019.02.054</a>

Jazyk výsledku

angličtina

Název v původním jazyce

Efficient and flexible MATLAB implementation of 2D and 3D elastoplastic problems

Popis výsledku v původním jazyce

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.

Název v anglickém jazyce

Efficient and flexible MATLAB implementation of 2D and 3D elastoplastic problems

Popis výsledku anglicky

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.

Druh

J_imp - Článek v periodiku v databázi Web of Science

CEP obor

—

OECD FORD obor

10102 - Applied mathematics

Projekt

Výsledek vznikl pri realizaci vícero projektů. Více informací v záložce Projekty.

Návaznosti

P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)

Rok uplatnění

2019

Kód důvěrnosti údajů

S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů

Název periodika

Applied Mathematics and Computation

ISSN

0096-3003

e-ISSN

—

Svazek periodika

355

Číslo periodika v rámci svazku

August 2019

Stát vydavatele periodika

US - Spojené státy americké

Počet stran výsledku

20

Strana od-do

595-614

Kód UT WoS článku

000464930500044

EID výsledku v databázi Scopus

2-s2.0-85063371842