This article deals with an algorithm for numerical modelling of hyperelastic heterogeneous materials undergoing large deformations. The microstructure of these materials is changing (deforming) during
a loading process, the changes in the microstructure depend on macroscopic deformations. To compute macroscopic responses, we must know material stiffness parameters and stresses in the heterogeneous structure. These effective parameters are obtained by solving microscopic problems. The number of microproblems is enormous, because in each iteration step (due to geometrical and material nonlinearities) it is needed to evaluate the effective material
parameters in each macroscopic quadrature point. To reduce a computational time a parallel algorithm is presented.