Most models currently used for complex phases in the calculation of phase diagrams (Calphad) method are based on the compound energy formalism. The way this formalism is presently used, however, is prone to poor extrapolation behavior in higher-order systems, especially when treating phases with complex crystal structures. In this paper, a partition of the Gibbs energy into effective bond energies, without changing its configurational entropy expression, is proposed, thereby remarkably improving the extrapolation behavior. The proposed model allows the use of as many sublattices as there are occupied Wyckoff sites and has great potential for reducing the number of necessary parameters, thus allowing shorter computational time. Examples for face centered cubic (fcc) ordering and the σ phase are given.