In this paper, we present a numerical method for solving nonlinear Hammerstein fractional integral equations. The method approximates the solution by Picard iteration together with a numerical integration designed for weakly singular integrals. Error analysis of the proposed method is also investigated. Numerical examples approve its efficiency in terms of accuracy and computational cost