Title | ||
---|---|---|
Rapid non-linear finite element analysis of continuous and discontinuous Galerkin methods in MATLAB. |
Abstract | ||
---|---|---|
MATLAB is adept at the development of concise Finite Element (FE) routines, however it is commonly perceived to be too inefficient for high fidelity analysis. This paper aims to challenge this preconception by presenting two optimised FE codes for both continuous Galerkin (CG) and discontinuous Galerkin (DG) methods. Although this has previously been achieved for linear-elastic problems, no such optimised MATLAB script currently exists, which includes the effects of material non-linearity. To incorporate these elasto-plastic effects, the externally applied load is split into a discrete number of loadsteps. Equilibrium is determined at each loadstep between the externally applied load and the arising internal forces using the Newton–Raphson method. The optimisation of the scripts is primarily achieved using vectorised blocking algorithms, which minimise RAM-to-cache overheads and maximise cache reuse. |
Year | DOI | Venue |
---|---|---|
2019 | 10.1016/j.camwa.2019.03.012 | Computers & Mathematics with Applications |
Keywords | Field | DocType |
Elasto-plasticity,Finite element analysis,Discontinuous Galerkin,MATLAB code vectorisation | Discontinuous Galerkin method,Bottleneck,Mathematical optimization,Nonlinear system,MATLAB,Galerkin method,Direct stiffness method,Finite element method,Computational science,Solver,Mathematics | Journal |
Volume | Issue | ISSN |
78 | 9 | 0898-1221 |
Citations | PageRank | References |
0 | 0.34 | 0 |
Authors | ||
4 |
Name | Order | Citations | PageRank |
---|---|---|---|
S O'Sullivan | 1 | 0 | 0.34 |
R.E. Bird | 2 | 0 | 0.34 |
William M. Coombs | 3 | 0 | 2.37 |
Stefano Giani | 4 | 36 | 9.55 |