In the present paper, the development of "special hole-elements", to enable an efficient and accurate analysis of stress concentration around through-thickness holes in angle-ply laminates, is presented. In these "hole-elements", the development of which is based on a modified complementary energy principle, the analytical asymptotic solutions for the stress-state near the hole are embedded. The fully 3-D stress-state in the laminate is accounted for, and the interlaminar traction reciprocity is satisfied a priori, while the interelement reciprocity condition is satisfied a posteriori through a Lagrange multiplier method. In addition, a simple method of estimation of stress concentration factors is also given. Results obtained from the present "special-hole-element" procedure are compared with the solutions by the present simple estimation method as well as other available solutions in literature.