IDEAS home Printed from https://ideas.repec.org/a/eee/apmaco/v424y2022ics0096300322001345.html
   My bibliography  Save this article

Fast MATLAB evaluation of nonlinear energies using FEM in 2D and 3D: Nodal elements

Author

Listed:
  • Moskovka, Alexej
  • Valdman, Jan

Abstract

Nonlinear energy functionals appearing in the calculus of variations can be discretized by the finite element (FE) method and formulated as a sum of energy contributions from local elements. A fast evaluation of energy functionals containing the first order gradient terms is a central part of this contribution. We describe a vectorized implementation using the simplest linear nodal (P1) elements in which all energy contributions are evaluated all at once without the loop over triangular or tetrahedral elements. Furthermore, in connection to the first-order optimization methods, the discrete gradient of energy functional is assembled in a way that the gradient components are evaluated over all degrees of freedom all at once. The key ingredient is the vectorization of exact or approximate energy gradients over nodal patches. It leads to a time-efficient implementation at higher memory-cost. Provided codes in MATLAB related to 2D/3D hyperelasticity and 2D p-Laplacian problem are available for download and structured in a way it can be easily extended to other types of vector or scalar forms of energies.

Suggested Citation

  • Moskovka, Alexej & Valdman, Jan, 2022. "Fast MATLAB evaluation of nonlinear energies using FEM in 2D and 3D: Nodal elements," Applied Mathematics and Computation, Elsevier, vol. 424(C).
  • Handle: RePEc:eee:apmaco:v:424:y:2022:i:c:s0096300322001345
    DOI: 10.1016/j.amc.2022.127048
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0096300322001345
    Download Restriction: Full text for ScienceDirect subscribers only

    File URL: https://libkey.io/10.1016/j.amc.2022.127048?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    As the access to this document is restricted, you may want to search for a different version of it.

    References listed on IDEAS

    as
    1. Čermák, M. & Sysala, S. & Valdman, J., 2019. "Efficient and flexible MATLAB implementation of 2D and 3D elastoplastic problems," Applied Mathematics and Computation, Elsevier, vol. 355(C), pages 595-614.
    2. Anjam, I. & Valdman, J., 2015. "Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements," Applied Mathematics and Computation, Elsevier, vol. 267(C), pages 252-263.
    Full references (including those not matched with items on IDEAS)

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Miroslav Frost & Jan Valdman, 2022. "Vectorized MATLAB Implementation of the Incremental Minimization Principle for Rate-Independent Dissipative Solids Using FEM: A Constitutive Model of Shape Memory Alloys," Mathematics, MDPI, vol. 10(23), pages 1-17, November.
    2. Innerberger, Michael & Praetorius, Dirk, 2023. "MooAFEM: An object oriented Matlab code for higher-order adaptive FEM for (nonlinear) elliptic PDEs," Applied Mathematics and Computation, Elsevier, vol. 442(C).

    Most related items

    These are the items that most often cite the same works as this one and are cited by the same works as this one.
    1. Dang-Bao Tran & Jaroslav Navrátil & Martin Čermák, 2021. "An Efficiency Method for Assessment of Shear Stress in Prismatic Beams with Arbitrary Cross-Sections," Sustainability, MDPI, vol. 13(2), pages 1-17, January.
    2. Yan-Ran Li & Xin-Hui Shao & Shi-Yu Li, 2021. "New Preconditioned Iteration Method Solving the Special Linear System from the PDE-Constrained Optimal Control Problem," Mathematics, MDPI, vol. 9(5), pages 1-13, March.
    3. Miroslav Frost & Jan Valdman, 2022. "Vectorized MATLAB Implementation of the Incremental Minimization Principle for Rate-Independent Dissipative Solids Using FEM: A Constitutive Model of Shape Memory Alloys," Mathematics, MDPI, vol. 10(23), pages 1-17, November.
    4. Vasileios E. Kontosakos, 2020. "Fast Quadratic Programming for Mean-Variance Portfolio Optimisation," SN Operations Research Forum, Springer, vol. 1(3), pages 1-15, September.
    5. Lukáš Pospíšil & Martin Čermák & David Horák & Jakub Kružík, 2020. "Non-Monotone Projected Gradient Method in Linear Elasticity Contact Problems with Given Friction," Sustainability, MDPI, vol. 12(20), pages 1-11, October.
    6. Innerberger, Michael & Praetorius, Dirk, 2023. "MooAFEM: An object oriented Matlab code for higher-order adaptive FEM for (nonlinear) elliptic PDEs," Applied Mathematics and Computation, Elsevier, vol. 442(C).
    7. Voet, Yannis, 2023. "On the fast assemblage of finite element matrices with application to nonlinear heat transfer problems," Applied Mathematics and Computation, Elsevier, vol. 436(C).
    8. Meijun Zhou & Jiayu Qin & Zenan Huo & Fabio Giampaolo & Gang Mei, 2022. "epSFEM: A Julia-Based Software Package of Parallel Incremental Smoothed Finite Element Method (S-FEM) for Elastic-Plastic Problems," Mathematics, MDPI, vol. 10(12), pages 1-25, June.

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:apmaco:v:424:y:2022:i:c:s0096300322001345. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If CitEc recognized a bibliographic reference but did not link an item in RePEc to it, you can help with this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Catherine Liu (email available below). General contact details of provider: https://www.journals.elsevier.com/applied-mathematics-and-computation .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.