Abstract:This paper presents a novel hybrid-stress finite element method for solving elastic fields around elliptical inclusions. A super polygonal-sided element containing an inclusion is constructed to reflect elastic behavior around the inclusion in terms of the variational principle of modified Hellinger-Reissner functionals. Displacement and stress fields in the element are expressed as complex series with Laurent series and Faber series by using the Muskhelishvili complex potential method and the conformal transformation technique. The super element is in conjunction with standard 4-node hybrid-stress elements to establish a new hybrid-stress finite element method. Benchmark example in the paper shows present method is rather general and efficient, as well as yields numerical results with higher accuracy and fewer elements. As an extended application of the method, an infinite isotropic plate containing two inclusions is analyzed under remote tension load. Effects of two inclusion span and elastic stiffness ratio on stress concentration coefficient are also discussed.