“…The use of first derivatives of the deflection in kernels of integrals related to the GNL effect, instead of the second derivatives, and the fact that there was no need of relating the derivatives of in-plane forces were the main features of the alternative formulation. This study employs the DRM to replace the domain integral by equivalent boundary integrals in the Soares and Palermo formulation [23]. The numerical implementation employed quadratic shape functions to approximate displacements, plate rotations, distributed shears and moments in the boundary elements as implemented in the study [23], as well as the inverse iteration and Rayleigh quotient to compute the lowest eigenvalue with the corresponding eigenvector.…”