Compositional dependence of the magnetic properties such as the saturation magnetization and the magnetic anisotropy energy were analysed using the three magnetic subnetwork mean field model and these results were compared to the radial distribution functions obtained from X-ray diffraction spectra and the hyperfine field distribution from Mossbauer spectroscopies. The ratio of the ferromagnetic Fe-subnetwork to the antiferromagnetic Fe-subnetwork and the spin value of Fe varied with the composition and all of Fe-subnetworks are coupled antiparallel to the Tb-subnetwork. Also, it was found that the hyperfine field at Fe nuclear is affected significantly by the conduction electron polarization due to the magnetic mean field near the center Fe atom as well as by the core electron polarization.