Efficient Implementation of Adaptive P1-FEM in Matlab

Funken, Stefan; Praetorius, Dirk; Wissgott, Philipp
October 2011
Computational Methods in Applied Mathematics;2011, Vol. 11 Issue 4, p460
Academic Journal
We provide a MATLAB package p1afem for an adaptive P1-finite element method (AFEM). This includes functions for the assembly of the data, different error estimators, and an indicator-based adaptive mesh-refining algorithm. Throughout, the focus is on an efficient realization by use of MATLAB built-in functions and vectorization. Numerical experiments underline the efficiency of the code which is observed to be of almost linear complexity with respect to the runtime. Although the scope of this paper is on AFEM, the general ideas can be understood as a guideline for writing efficient MATLAB code.


Related Articles

  • Numerical examples of variable Three-node beam elements based on Positional FEM. Lvzhou Ma; Yuqin Yan; Xunlin Diao; Jian Liu // Advanced Materials Research;2014, Vol. 988, p471 

    Based on MATLAB platform, program TBEP (Three-node Beam Element based on Positional FEM) is compiled. Pure bending of Euler beam and cubic diamond-shaped rigid frame are calculated. The calculation results show that positional FEM uses fewer elements and gains higher calculation precision and...

  • Dual-primal FETI algorithms for edge finite-element approximations in 3D. Toselli, Andrea // IMA Journal of Numerical Analysis;Jan2006, Vol. 26 Issue 1, p96 

    A family of dual-primal finite-element tearing and interconnecting methods for edge-element approximations in 3D is proposed and analysed. The key part of this work relies on the observation that for these finite-element spaces there is a strong coupling between degrees of freedom associated...

  • Porównanie adekwatnoÅ›ci wybranych metod numerycznych do opisu pola przepÅ‚ywu w wentylatorze poprzecznym. Stacharska-Targosz, Jolanta; Chmielowiec, Monika // Systems: Journal of Transdisciplinary Systems Science;2006, Vol. 11 Issue 1, p268 

    The complexity of flow structure and the interactions of design parameters are the main reason of lack of an acceptable design procedure and accurate theoretical flow field model describing the real flow conditions in the cross flow fans. There were some theoretical approaches created the...

  • A Generalized Hard Thresholding Pursuit Algorithm. Li, Haifeng; Fu, Yuli; Zhang, Qiheng; Rong, Rong // Circuits, Systems & Signal Processing;Apr2014, Vol. 33 Issue 4, p1313 

    Compressed sensing ensures the accurate reconstruction of sparse signals from far fewer samples than required in the classical Shannon-Nyquist theorem. In this paper, a generalized hard thresholding pursuit (GHTP) algorithm is presented that can recover unknown vectors without the sparsity level...

  • Computing intersections between non-compatible curves and finite elements. Durand, Raul; Farias, Márcio; Pedroso, Dorival // Computational Mechanics;Sep2015, Vol. 56 Issue 3, p463 

    This paper presents a method to find all intersections between curved lines such as structural line elements and finite element meshes with intentions to generate smaller, non-compatible, line cells (e.g. bar elements) between crossings. The intersection finding algorithm works for two and...

  • An alternative contact/impact identification algorithm for 2d structural problems. Greco, M.; Coda, H. B.; Venturini, W. S. // Computational Mechanics;Oct2004, Vol. 34 Issue 5, p410 

    This paper presents a new algorithm to identify the occurrence of impact (contact) among structures. Structures are modeled here by finite elements and the resulting formulation is employed to solve two-dimensional frame impact problems by the Lagrange Multiplier technique. The proposed...

  • A Filled Function Algorithm for Multiobjective Optimization. LIUYANG YUAN; ZHONGPING WAN; JIAWEI CHEN // WSEAS Transactions on Mathematics;Aug2012, Vol. 11 Issue 8, p659 

    In this paper a filled function algorithm is applied to compute one of the nonisolated Pareto optimal points of an unconstrained multiobjective optimization problem. Firstly, the original problem is converted into an equivalent global optimization problem. Subsequently, a novel filled function...

  • Finite Element Analysis of Polycrystalline Deformation with the Rate-dependent Crystal Plasticity. Yoon, J. H.; Huh, H.; Lee, Y. S. // AIP Conference Proceedings;2007, Vol. 908 Issue 1, p1325 

    Constitutive models for the crystal plasticity have the common objective which relates the behavior of microscopic single crystals in the crystallographic texture to the macroscopic continuum response. This paper presents the texture analysis of polycrystalline materials using the rate-dependent...

  • Optimisation of Shape Parameters and Process Manufacturing for an Automotive Safety Part. Gildemyn, Eric; Dal Santo, Philippe; Potiron, Alain; Saïdane, Delphine // AIP Conference Proceedings;2007, Vol. 908 Issue 1, p1219 

    In recent years, the weight and the cost of automotive vehicles have considerably increased due to the importance devoted to safety systems. It is therefore necessary to reduce the weight and the production cost of components by improving their shape and manufacturing process. This work deals...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics