An iterative eshelby method to calculate average stress fields and elastic moduli of composite and multi-phase materials