In this paper, numerical methods of fatigue life prediction for elastomers subjected to multidirectional, variable amplitude loadings are presented. Because experiments and numerical methods use different stress measures in large deformation, transformation between nominal stress and the second Piola–Kirchhoff stress is performed before fatigue life calculation. In order to incorporate the Mullins effect, the material properties of elastomers are calculated after an initial transition period. An efficient interpolation scheme using load stress/strain curves under unidirectional loading is proposed based on the fatigue characteristic of elastomers. A rainflow counting method with multi‐stress components is developed for variable amplitude loadings, and the critical plane method is applied to find the plane with the maximum damage parameter. Fatigue life predictions using the proposed numerical method are validated against experimental results. As a practical example, the fatigue life of a rubber engine mount is predicted using the proposed numerical method.