Development of computational modeling tools has revolutionized studies of magmatic processes over the last four decades. Their refinement from binary mixing equations to thermodynamically controlled geochemical assimilation models has provided more comprehensive and detailed modeling constraints of an array of magmatic systems. One of the questions that has not yet been vigorously studied using thermodynamic constraints is the origin of massif-type anorthosites. The parental melts to these intrusions are hypothesized to be either mantle-derived high-Al basaltic melts that undergo crustal contamination or monzodioritic melts derived directly from lower crust. On the other hand, many studies suggest that the monzodioritic rocks do not represent parental melts but instead represent crystal remnants of residual liquids left after crystal fractionation of parental melts. Regardless of the source or composition, magmas that produce massif-type anorthosites have been suggested to have undergone polybaric (~1000–100 MPa) fractional crystallization while ascending through the lithosphere. We conducted lower crustal melting, assimilation-fractional crystallization, and isobaric and polybaric fractional crystallization major element modeling using two thermodynamically constrained modeling tools, the Magma Chamber Simulator (MCS) and rhyolite-MELTS, to test the suitability of these tools and to study the petrogenesis of massif-type anorthosites. Comparison of our models with a large suite of whole-rock data suggests that the massif-type anorthosite parental melts were high-Al basalts that were produced when hot mantle-derived partial melts assimilated lower crustal material at Moho levels. These contaminated basaltic parental magmas then experienced polybaric fractional crystallization at different crustal levels (~40 to 5 km) producing residual melts that crystallized as monzodioritic rocks. Model outcomes also support the suggestion that the cumulates produced during polybaric fractional crystallization likely underwent density separation, thus producing the plagioclase-rich anorthositic rocks. The modeled processes are linked to a four-stage model that describes the key petrogenetic processes that generate massif-type anorthosites. The presented framework enables further detailed thermodynamic and geochemical modeling of individual anorthosite intrusions using MCS and involving trace element and isotope constrains. • Thermodynamic constraints on the origin of massif-type anorthosites using MCS. • Models provide indications for the composition and source of parental magmas. • The parental magmas were generated by AFC processes of mantle-derived melts. • The models are linked to general magmatic processes during anorthosite formation. • The general models set the foundation for detailed modeling of individual intrusions.