“…Computational tools based on integral equations and on the employment of the boundary element method (BEM) [15,16] for their solution have been already successfully developed for the analysis of polycrystalline materials, both in 2D [17] and 3D [18,19], and have been applied to either computational homogenization [20], multiscale materials modelling [21,22,23] or micro-cracking analysis [24,25,26], also considering piezo-electric polycrystals [27], high-cycle and low-cycle fatigue [28,29] and hydrogen assisted cracking [30]. Such formulations, often built on Voronoi tessellations [31], which provide a reasonable approximation of the microstructural morphology [32,33], are expressed uniquely in terms of displacements and tractions of points belonging to the boundary of the crystals in the aggregate, thus providing a reduction in the number of degrees of freedom needed for analysing a give microstructure.…”