Stabilizer-free polygonal and polyhedral virtual elements

A variationally consistent numerical approach based on the Virtual Element Method (VEM) is presented for the analysis of 2D elastoplasticity problems. The mixed Hu-Washizu functional of elasticity is extended to incorporate the energy contributions specific to the finite-step elastoplastic problem. It is demonstrated how the governing equations of the discretized elastoplastic problem - including the loading-unloading conditions - emerge naturally as the stationarity conditions of the VEM-discretized functional. Spurious hourglass modes are prevented by formulating a self-stabilized version of Virtual Elements (VEs) that exploits the possibility offered by the mixed approach to define strain and displacement approximations of the same order. The insensitivity of VEs to element distortion and the possibility to use polygonal elements with any shape and number of edges is tested with the analysis of several benchmarks from the literature. It is shown how accurate solutions can be obtained also in the case of non-convex quadrilateral or pentagonal elements. Additionally, the role of internal moment degrees of freedom in preventing elastoplastic locking at the plastic failure limit is elucidated.

By introducing the shape functions of virtual node polygonal (VP) elements into the standard extended finite element method (XFEM), a conforming elemental mesh can be created for the cracking process. Moreover, an adaptively refined meshing with the quadtree structure only at a growing crack tip is proposed without inserting hanging nodes into the transition region. A novel dynamic crack growth method termed as VP-XFEM is thus formulated in the framework of fracture mechanics. To verify the newly proposed VP-XFEM, both quasi-static and dynamic cracked problems are investigated in terms of computational accuracy, convergence, and efficiency. The research results show that the present VP-XFEM can achieve good agreement in stress intensity factor and crack growth path with the exact solutions or experiments. Furthermore, better accuracy, convergence, and efficiency of different models can be acquired, in contrast to standard XFEM and mesh-free methods. Therefore, VP-XFEM provides a suitable alternative to XFEM for engineering applications.

The present work deals with the formulation of a virtual element method for two dimensional structural problems. The contribution is split in two parts: in part I, the elastic problem is discussed, while in part II (Artioli et al. in Comput Mech, 2017) the method is extended to material nonlinearity, considering different inelastic responses of the material. In particular, in part I a standardized procedure for the construction of all the terms required for the implementation of the method in a computer code is explained. The procedure is initially illustrated for the simplest case of quadrilateral virtual elements with linear approximation of displacement variables on the boundary of the element. Then, the case of polygonal elements with quadratic and, even, higher order interpolation is considered. The construction of the method is detailed, deriving the approximation of the consistent term, the required stabilization term and the loading term for all the considered virtual elements. A wide numerical investigation is performed to assess the performances of the developed virtual elements, considering different number of edges describing the elements and different order of approximations of the unknown field. Numerical results are also compared with the one recovered using the classical finite element method.

The present work deals with the formulation of a virtual element method for two dimensional structural problems. The contribution is split in two parts: in part I, the elastic problem is discussed, while in part II (Artioli et al. in Comput Mech, 2017) the method is extended to material nonlinearity, considering different inelastic responses of the material. In particular, in part I a standardized procedure for the construction of all the terms required for the implementation of the method in a computer code is explained. The procedure is initially illustrated for the simplest case of quadrilateral virtual elements with linear approximation of displacement variables on the boundary of the element. Then, the case of polygonal elements with quadratic and, even, higher order interpolation is considered. The construction of the method is detailed, deriving the approximation of the consistent term, the required stabilization term and the loading term for all the considered virtual elements. A wide numerical investigation is performed to assess the performances of the developed virtual elements, considering different number of edges describing the elements and different order of approximations of the unknown field. Numerical results are also compared with the one recovered using the classical finite element method.

Phase-field modeling of brittle fracture using an efficient virtual element scheme

The geometric shape of an element plays a key role in computational methods. Triangular and quadrilateral shaped elements are utilized by standard finite element methods. The pioneering work of Wachspress laid the foundation for polygonal interpolants which introduced polygonal elements. Tessellations may be considered as the next stage of element shape evolution. In this work, we investigate the topology optimization of tessellations as a means to coalesce art and engineering. We mainly focus on M.C. Escher's tessellations using recognizable figures. To solve the state equation, we utilize a Mimetic Finite Difference inspired approach, known as the Virtual Element Method. In this approach, the stiffness matrix is constructed in such a way that the displacement patch test is passed exactly in order to ensure optimum numerical convergence rates. Prior to exploring the artistic aspects of topology optimization designs, numerical verification studies such as the displacement patch test and shear loaded cantilever beam bending problem are conducted to demonstrate the accuracy of the present approach in two-dimensions.