An analytical study of elastoviscoplastic fluid flow in tubes of non-circular cross section is presented. The constitutive structure of the fluid is described by a linear frame invariant combination of the Phan-Thien−Tanner model of viscoelastic fluids and the Bingham model of plastic fluids. Non-circular tube cross sections are modeled by the shape factor method a one-to-one mapping of the circular base contour into a wide spectrum family of arbitrary tube contours. Field variables are expanded into asymptotic series in terms of the elasticity measure, the Weissenberg number We, coupled with an asymptotic expansion in terms of the geometrical mapping parameter ε leading to a set of hierarchical momentum balance equations which are solved successively up to and including the third order in We when the secondary field appears for the first time. The computational algorithm developed is applied to the study of the non-rectilinear flow in tubes with triangular and square cross sections. We find that the presence of the yield stress dampens the intensity of the purely viscoelastic vortices, the higher the yield stress the lower the intensity of the vortices in the cross-section, and the further away the vortices are from the center of the cross section as compared to the purely viscoelastic vortices. The results also evidence that viscoelasticity increases the axial flow for given viscoplastic conditions and pressure drop, and consequently increases the rate of flow, a phenomenon that may find applications in optimizing material transportation.