In this work, Monte Carlo simulation was used to model the dynamic hysteresis behavior of magnetic/nonmagnetic multilayers in the framework of Heisenberg model using the spin-flip algorithm. The purpose is to investigate the hysteresis properties while varying the direction of external applied magnetic field and temperature taking account interface composition of magnetic/non-magnetic atoms. From the results, with thermal variation and interface composition, the simulated hysteresis loops change significantly due to the competition between magnetic anisotropy in magnetic films and the presence of isolated clusters at the interface.