Within the framework of the displacement‐based virtual element method (VEM), namely, for plane elasticity, an important topic is the development of optimal techniques for the evaluation of the stress field. In fact, in the classical VEM formulation, the same projection operator used to approximate the strain field (and then evaluate the stiffness matrix) is employed to recover, via constitutive law, the stress field. Considering a first‐order formulation, strains are locally mapped onto constant functions, and stresses are piecewise constant. However, the virtual displacements might engender more complex strain fields for polygons, which are not triangles. This leads to an undesirable loss of information with respect to the underlying virtual stress field. The recovery by compatibility in patches, originally proposed for finite element schemes, is here extended to VEM, aiming at mitigating such an effect. Stresses are recovered by minimizing the complementary energy of patches of elements over an assumed set of equilibrated stress modes. The procedure is simple, efficient, and can be readily implemented in existing codes. Numerical tests confirm the good performance of the proposed technique in terms of accuracy and indicate an increase of convergence rate with respect to the classical approach in many cases.
An equilibrium-based stress recovery procedure for the VEM / Artioli, E; de Miranda, S; Lovadina, C; Patruno, L. - In: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING. - ISSN 0029-5981. - 117:8(2019), pp. 885-900. [10.1002/nme.5983]
An equilibrium-based stress recovery procedure for the VEM
Artioli E;
2019
Abstract
Within the framework of the displacement‐based virtual element method (VEM), namely, for plane elasticity, an important topic is the development of optimal techniques for the evaluation of the stress field. In fact, in the classical VEM formulation, the same projection operator used to approximate the strain field (and then evaluate the stiffness matrix) is employed to recover, via constitutive law, the stress field. Considering a first‐order formulation, strains are locally mapped onto constant functions, and stresses are piecewise constant. However, the virtual displacements might engender more complex strain fields for polygons, which are not triangles. This leads to an undesirable loss of information with respect to the underlying virtual stress field. The recovery by compatibility in patches, originally proposed for finite element schemes, is here extended to VEM, aiming at mitigating such an effect. Stresses are recovered by minimizing the complementary energy of patches of elements over an assumed set of equilibrated stress modes. The procedure is simple, efficient, and can be readily implemented in existing codes. Numerical tests confirm the good performance of the proposed technique in terms of accuracy and indicate an increase of convergence rate with respect to the classical approach in many cases.File | Dimensione | Formato | |
---|---|---|---|
Artioli_et_al-2018-International_Journal_for_Numerical_Methods_in_Engineering.pdf
Accesso riservato
Dimensione
1.94 MB
Formato
Adobe PDF
|
1.94 MB | Adobe PDF | Visualizza/Apri Richiedi una copia |
Pubblicazioni consigliate
I metadati presenti in IRIS UNIMORE sono rilasciati con licenza Creative Commons CC0 1.0 Universal, mentre i file delle pubblicazioni sono rilasciati con licenza Attribuzione 4.0 Internazionale (CC BY 4.0), salvo diversa indicazione.
In caso di violazione di copyright, contattare Supporto Iris