Advanced air defense methods are essential to address the growing complexity of aerial threats. The increasing number of targets necessitates better defensive coordination, and a promising strategy involves the use of interceptors together to protect a specific area. This task fundamentally depends on accurately predicting their kinematic envelopes, or reachable sets. This paper presents a novel approach to determine the boundaries of reachable sets for aerodynamic interceptors, accounting for energy loss from drag, energy gain from thrust, variable acceleration limits, and autopilot dynamics. The devised numerical method approximates reachable sets for nonlinear problems using a constrained model predictive programming concept. Results demonstrate that explicitly accounting for input constraints, such as acceleration limits, significantly impacts the shape and area of the reachable boundaries. Furthermore, a sensitivity analysis was conducted to demonstrate the impact of parameter variations on the reachable set. Revealing the reachable set’s sensitivity to variations in thrust and drag coefficients, this analysis serves as a framework for considering parameter uncertainty and enables the evaluation of these effects prior to embedding the reachability boundaries into an offline database for guidance applications. The resulting boundaries, representing minimum and maximum ranges for various initial parameters, can be stored offline, allowing interceptors to estimate their own or allied platforms’ kinematic capabilities for cooperative strategies.