The proliferation of distributed generation on distribution feeders triggers a large number of integration and planning studies. Further, the complexity of distribution feeder models, short simulation time steps, and long simulation horizons rapidly render studies computational burdensome. To mend this issue, we propose a methodology for reducing the number of nodes, loads, generators, line, and transformers of p-phase distribution feeders with unbalanced loads and generation, nonsymmetric wire impedance, mutual coupling, shunt capacitance, and changes in voltage and phase. The methodology is derived on a constant power load assumption and employs a Gaussian elimination inversion technique to design the reduced feeder. Compared to previous work by the authors, the inversion reduction takes half the time and voltage errors after reduction are reduced by an order of magnitude. Using a snapshot simulation the reduction is tested on six additional publicly available feeders with a maximum voltage error 0.0075 p.u. regardless of feeder size or complexity, and typical errors on the order of 1×10 −4 p.u. For a day long QSTS simulation on the UCSD A feeder, errors are shown to increase with changes in loading when a large number of buses removed, but shows less variation for less than 85% of buses removed.