Ising model is a promising tool for simulation magnetic Barkhausen noise (MBN) signal for nondestructive evaluation. However, theoretical modeling of two-phase or multi-phase ferromagnetic materials are seldom reported by this model. In order to fill the blank, we extended Ising model combine with double Boltzmann transition function to simulated MBN signal of two-phase ferromagnetic materials. The volume fraction was considered as a vital parameter which effect the shape of MBN profiles directly. A continuous transition form was assumed and the proposed model was simulated by Monte Carlo algorithm. The second order Gaussian curve was used to fit the two-peak shape MBN profile. This simulation show that the first peak in lower field increased and the second peak in higher field gradually decreased with the first phase simulated volume fraction increases. In addition, the ratio feature about the peaks amplitude show a linear increased trend. The proposed model offer a feasible tool for predict MBN in two-phase materials.