This work details a methodology for parameterization of the kinetics and thermodynamics of the thermal decomposition of polymers blended with reactive additives. This methodology employs Thermogravimetric Analysis, Differential Scanning Calorimetry, Microscale Combustion Calorimetry, and inverse numerical modeling of these experiments. Blends of glass-fiber-reinforced polybutylene terephthalate (PBT) with aluminum diethyl phosphinate and melamine polyphosphate were used to demonstrate this methodology. These additives represent a potent solution for imparting flame retardancy to PBT. The resulting lumped-species reaction model consisted of a set of first- and second-order (two-component) reactions that defined the rate of gaseous pyrolyzate production. The heats of reaction, heat capacities of the condensed-phase reactants and products, and heats of combustion of the gaseous products were also determined. The model was shown to reproduce all aforementioned experiments with a high degree of detail. The model also captured changes in the material behavior with changes in the additive concentrations. Second-order reactions between the material constituents were found to be necessary to reproduce these changes successfully. The development of such models is an essential milestone toward the intelligent design of flame retardant materials and solid fuels.