We develop an effective Lagrangian for QCD in terms of order parameters of the gluon condensate. (A useful analogy to the Ginzburg-Landau theory of superconductivity may be made.) In this work the order parameters are identified with averages of the low-momentum components of the vector potential. The theory is formulated such that only a single order parameter, which describes the number of gluon pairs in the condensate (of all colors), is relevant. This order parameter becomes spatially dependent in the presence of quarks and one obtains a dynamical model for the formation of hadrons as nontopological soliton solutions of the effective Lagrangian of our model.