Přihlášení pro studenty

Přihlášení pro zaměstnance

Publikace detail

Vectfem: a generalized MATLAB-based vectorized algorithm for the computation of global matrix/force for finite elements of any type and approximation order in linear elasticity
Rok: 2024
Druh publikace: článek v odborném periodiku
Název zdroje: Zeitschrift fur Angewandte Mathematik und Physik
Strana od-do: 1-37
Tituly:
Jazyk Název Abstrakt Klíčová slova
cze Vectfem: zobecněný vektorizovaný algoritmus založený na MATLABu pro výpočet globální matice/síly pro konečné prvky libovolného typu a aproximace v lineární elasticitě V tomto článku představujeme nový vektorizovaný algoritmus založený na MATLABu pro efektivní sériový výpočet globální matice/síly vyplývající z metody konečných prvků (MKP) pro sítě libovolného typu a řádu aproximace v lineární elasticitě. metoda konečných prvků; algoritmus; vektorizace; pružnost; prostor; MATLAB
eng Vectfem: a generalized MATLAB-based vectorized algorithm for the computation of global matrix/force for finite elements of any type and approximation order in linear elasticity In this paper, we introduce a new vectorized MATLAB-based algorithm for efficient serial computation of global matrix/force arising from finite element method (FEM) for meshes of any type and approximation order in linear elasticity. Because for-loops in MATLAB are very slow, we propose a modified process that takes advantage of vectorization and sparse assembly to achieve good performance while using the same memory as the standard algorithm. For this purpose, by using good programming practices, the implementation of this scheme is succinctly described and can be integrated into any MATLAB package dealing with FEM. Specifically, attention is paid to the calculation of the triplet (row index, column index, matrix components) as well as the assembly of the global stiffness matrix, mass matrix and force vector. Additionally, an extension of the proposed approach for Mindlin plate theory and functionally graded materials is outlined. Finally, the accuracy of this strategy is verified on selected numerical tests after comparing the obtained results with those of ABAQUS. In terms of performance, the study conducted on a set of meshes considering the standard algorithm and two other well-known MATLAB vectorized algorithms revealed that: (i) for a 2D beam problem meshed with P1-triangle elements, a speedup of about 8 and 15 is achieved with sparse and fsparse, respectively. (ii) for a 3D plate problem meshed with P1-tetrahedral elements, a speedup of about 4 and 8 is achieved with sparse and fsparse, respectively. When compared to ABAQUS performance, the proposed scheme results in a computational time that is about five times smaller. finite element method; algorithm; vectorization; elasticity; sparse; MATLAB

Katedry

Pracoviště

Jak nás najdete?

zobrazit na mapy.cz

Služby

Často hledáte

© 2023 Univerzita Pardubice, Studentská 95, 532 10 Pardubice 2
Telefon: 466 036 111, 466 036 112, 466 036 113
Správce webu, RSS
ID datové schránky: f5vj9hu
Prohlášení o přístupnosti
CC BY-NC-ND 4.0 CZ