Dislocation-based finite element method for homogenized limit domain characterization of structured metamaterials

Abstract

PurposeHeterogeneous and micro-structured materials have been the object of multiscale and homogenization techniques aimed at recognizing the elastic properties of the equivalent continuum. The proposed investigation deals with the mechanical characterization of the heterogeneous material structured metamaterials through analyzing the ultimate strength using the limit analysis of the Representative Volume Element (RVE). To get the desired material strength, a novel finite element formulation based on the derivation of self-equilibrated solutions through the finite elements devoted to calculating the lower bound theorem has been implemented together with the limit analysis in Melàn’s formulation.Design/methodology/approachThe finite element formulation is based on discrete mapping of Volterra dislocations in the structure using isoparametric representation. Using standard finite element techniques, the linear operator V, which relates the self-equilibrated internal solicitation to displacement-like nodal parameters, has been built through finite element discretization of displacement and strain.FindingsThe proposed work presented an elastic homogenization of the mechanical properties of an elementary cell with a geometry known in the literature, the isotropic truss. The matrix of elastic constants was calculated by subjecting the RVE to numerical load tests, simulated with a commercial FEM calculation code. This step showed the dependence of the isotropy properties, verified with Zener theory, on the density of the RVE. The isotropy condition of the material is only achieved for certain section ratios between body-centered cubic (BCC) and face-centered cubic (FCC), neglecting flexural effects at the nodes. The density that satisfies Zener’s conditions represents the isotropic geomatics of the isotropic truss.Originality/valueFor the isotropic case, the VFEM procedure was used to evaluate the isotropy of the limit domain and was compared with the Mises–Schleicher limit domain. The evaluation of residual ductility and dissipation energy allowed a measurement parameter for the limit anisotropy to be defined. The novelty of the proposal consisted in the formulation of both the linearized and the nonlinear limit locus of the material; hence, it furnished the starting point for further limit analysis of the structures whose elementary volume has been described through the proposed approach.