The Eshelby-type problems of an arbitrary-shape polygonal or polyhedral inclusion embedded in an infinite homogeneous isotropic elastic material are analytically solved using a simplified strain gradient elasticity theory (SSGET) that contains a material length scale parameter. The Eshelby tensors for a plane strain inclusion with an arbitrary polygonal cross section and for an arbitrary-shape polyhedral inclusion are analytically derived in general forms in terms of three potential functions. These potential functions, as area integrals over the polygonal cross section and volume integrals over the polyhedral inclusion, are evaluated. For the polygonal inclusion problem, the three area integrals are first transformed to three line integrals using the Green's theorem, which are then evaluated analytically by direct integration. In the polyhedral inclusion case, each of the three volume integrals is first transformed to a surface integral by applying the divergence theorem, which is then transformed to a contour (line) integral based on Stokes' theorem and using an inverse approach. In addition, the Eshelby tensor for an anti-plane strain inclusion with an arbitrary polygonal cross section embedded in an infinite homogeneous isotropic elastic material is analytically solved. Each of the newly derived Eshelby tensors is separated into a classical part and a gradient part. The latter includes the material length scale parameter additionally, thereby enabling the interpretation of the inclusion size effect. For homogenization applications, the area or volume average of each newly derived Eshelby tensor over the polygonal cross section or the polyhedral inclusion domain is also provided in a general form. To illustrate the newly obtained Eshelby tensors and their area or volume averages, different types of polygonal and polyhedral inclusions are quantitatively studied by directly using the general formulas derived. The numerical results show that the components of the each SSGET-based Eshelby tensor for all inclusion shapes considered vary with both the position and the inclusion size. It is also observed that the components of each averaged Eshelby tensor based on the SSGET change with the inclusion size.