In the present study, we propose an evaluation of a non-linear homogenization model applied to hyperelastic composites having random microstructure. This modelling approach consists in a 3D implementation of the secondorder method introduced by Ponte Castañeda & Tiberio [1]. We first recall the basic principles of this method. Then, we investigate a composite made up of an hyperelastic matrix reinforced by spherical deformable or rigid particles.Computational issues of the micromechanical model are discussed and some obtained results allowto demonstrate the reinforcement effect of the particles. In order to provide a rigorous evaluation of the methodology, finite elements computations, on an unit cell, are performed and compared to the predictions of the model. Finally, a confrontation with experimental results is provided.