Digital twin of design and manufacturing process of reinforced elastomeric matrices are crucial for evaluating the mechanical behavior of structural component such as tires and gaskets. This study focuses on estimating the effective behavior of an elastomeric matrix reinforced with mica spherical particles using full-field homogenization. Mean as well as the statistical variation of local field quantities i.e. norm of deviatoric Cauchy stress, and Eulerian strain in both of the phases are investigated. Effect of Poisson’s ratio, volume fraction and loading conditions are also assessed for the same. The hyperelastic matrix is modeled with the Saint Venant-Kirchhoff (SVK) and Neo-Hookean (NH) models to consider the effect of geometric and material nonlinearity respectively. The inhomogeneities are treated as SVK. A representative random artificial microstructure is generated using the random sequential adsorption algorithm. The representative volume element is simulated in finite element(FE) framework for strain based loading conditions by applying periodic boundary conditions. Hybrid elements are employed in FE simulations to prevent numerical instabilities occurring due to incompressible nature of hyperelastic materials. It is observed that stress fluctuates strongly compared to the strain in the particulate phase. In nearly incompressible scenarios the Poisson’s ratio of the matrix phase is very sensitive in estimating effective behavior.