On the equivalence between the cell-based smoothed finite element method and the virtual element method

We revisit the cell-based smoothed finite element method (SFEM) for quadrilateral elements and extend it to arbitrary polygons and polyhedrons in 2D and 3D, respectively. We highlight the similarity between the SFEM and the virtual element method (VEM). Based on the VEM, we propose a new stabilization approach to the SFEM when applied to arbitrary polygons and polyhedrons. The accuracy and the convergence properties of the SFEM are studied with a few benchmark problems in 2D and 3D linear elasticity. Later, the SFEM is combined with the scaled boundary finite element method to problems involving singularity within the framework of the linear elastic fracture mechanics in 2D

CVEM-BEM Coupling with Decoupled Orders for 2D Exterior Poisson Problems

For the solution of 2D exterior Dirichlet Poisson problems, we propose the coupling of a Curved Virtual Element Method (CVEM) with a Boundary Element Method (BEM), by using decoupled approximation orders. We provide optimal convergence error estimates, in the energy and in the weaker L-2-norm, in which the CVEM and BEM contributions to the error are separated. This allows for taking advantage of the high order flexibility of the CVEM to retrieve an accurate discrete solution by using a low order BEM. The numerical results confirm the a priori estimates and show the effectiveness of the proposed approach

Virtual Element based formulations for computational materials micro-mechanics and homogenization

In this thesis, a computational framework for microstructural modelling of transverse behaviour of heterogeneous materials is presented. The context of this research is part of the broad and active field of Computational Micromechanics, which has emerged as an effective tool both to understand the influence of complex microstructure on the macro-mechanical response of engineering materials and to tailor-design innovative materials for specific applications through a proper modification of their microstructure. While the classical continuum approximation does not account for microstructural details within the material, computational micromechanics allows detailed modelling of a heterogeneous material's internal structural arrangement by treating each constituent as a continuum. Such an approach requires modelling a certain material microstructure by considering most of the microstructure's morphological features. The most common numerical technique used in computational micromechanics analysis is the Finite Element Method (FEM). Its use has been driven by the development of mesh generation programs, which lead to the quasi-automatic discretisation of the artificial microstructure domain and the possibility of implementing appropriate constitutive equations for the different phases and their interfaces. In FEM's applications to computational micromechanics, the phase arrangements are discretised using continuum elements. The mesh is created so that element boundaries and, wherever required, special interface elements are located at all interfaces between material's constituents. This approach can be effective in modelling many microstructures, and it is readily available in commercial codes. However, the need to accurately resolve the kinematic and stress fields related to complex material behaviours may lead to very large models that may need prohibitive processing time despite the increasing modern computers' performance. When rather complex microstructure's morphologies are considered, the quasi-automatic discretisation process stated before might fail to generate high-quality meshes. Time-consuming mesh regularisation techniques, both automatic and operator-driven, may be needed to obtain accurate numeric results. Indeed, the preparation of high-quality meshes is today one of the steps requiring more attention, and time, from the analyst. In this respect, the development of computational techniques to deal with complex and evolving geometries and meshes with accuracy, effectiveness, and robustness attracts relevant interest. The computational framework presented in this thesis is based on the Virtual Element Method (VEM), a recently developed numerical technique that has proven to provide robust numerical results even with highly-distorted mesh. These peculiar features have been exploited to analyse two-dimensional representations of heterogeneous materials' microstructures. Ad-hoc polygonal multi-domain meshing strategies have been developed and tested to exploit the discretisation freedom that VEM allows. To further simplify the preprocessing stage of the analysis and reduce the total computational cost, a novel hybrid formulation for analysing multi-domain problems has been developed by combining the Virtual Element Method with the well-known Boundary Element Method (BEM). The hybrid approach has been used to study both composite material's transverse behaviour in the presence of inclusions with complex geometries and damage and crack propagation in the matrix phase. Numerical results are presented that demonstrate the potential of the developed framework

A virtual element method for the solution of 2D time-harmonic elastic wave equations via scalar potentials

In this paper, we propose and analyse a numerical method to solve 2D Dirichlet timeharmonic elastic wave equations. The procedure is based on the decoupling of the elastic vector field into scalar Pressure (P-) and Shear (S-) waves via a suitable Helmholtz– Hodge decomposition. For the approximation of the two scalar potentials we apply a virtual element method associated with different mesh sizes and degrees of accuracy. We provide for the stability of the method and a convergence error estimate in the L 2 -norm for the displacement field, in which the contributions to the error associated with the P- and S- waves are separated. In contrast to standard approaches that are directly applied to the vector formulation, this procedure allows for keeping track of the two different wave numbers, that depend on the P- and S- speeds of propagation and, therefore, for using a high-order method for the approximation of the wave associated with the higher wave number. Some numerical tests, validating the theoretical results and showing the good performance of the proposed approach, are presented

A survey of Trefftz methods for the Helmholtz equation

Trefftz methods are finite element-type schemes whose test and trial functions are (locally) solutions of the targeted differential equation. They are particularly popular for time-harmonic wave problems, as their trial spaces contain oscillating basis functions and may achieve better approximation properties than classical piecewise-polynomial spaces. We review the construction and properties of several Trefftz variational formulations developed for the Helmholtz equation, including least squares, discontinuous Galerkin, ultra weak variational formulation, variational theory of complex rays and wave based methods. The most common discrete Trefftz spaces used for this equation employ generalised harmonic polynomials (circular and spherical waves), plane and evanescent waves, fundamental solutions and multipoles as basis functions; we describe theoretical and computational aspects of these spaces, focusing in particular on their approximation properties. One of the most promising, but not yet well developed, features of Trefftz methods is the use of adaptivity in the choice of the propagation directions for the basis functions. The main difficulties encountered in the implementation are the assembly and the ill-conditioning of linear systems, we briefly survey some strategies that have been proposed to cope with these problems.Comment: 41 pages, 2 figures, to appear as a chapter in Springer Lecture Notes in Computational Science and Engineering. Differences from v1: added a few sentences in Sections 2.1, 2.2.2 and 2.3.1; inserted small correction