Gauge Freedom in Orbital Mechanics
  In orbital and attitude dynamics the coordinates and the Euler angles are expressed as functions of the time and six constants called elements. Under disturbance, the constants are endowed with time dependence. The Lagrange constraint is then imposed to guarantee that the functional dependence of the perturbed velocity on the time and constants stays the same as in the undisturbed case. Constants obeying this condition are called osculating elements. The constants chosen to be canonical are called Delaunay elements, in the orbital case, or Andoyer elements, in the spin case. (As some Andoyer elements are time dependent even in the free-spin case, the role of constants is played by their initial values.) The Andoyer and Delaunay sets of elements share a feature not readily apparent: in certain cases the standard equations render them non-osculating. In orbital mechanics, elements furnished by the standard planetary equations are non-osculating when perturbations depend on velocities. To preserve osculation, the equations must be amended with extra terms that are not parts of the disturbing function. In the case of Delaunay parameterisation, these terms destroy canonicity. So under velocity-dependent disturbances, osculation and canonicity are incompatible. (Efroimsky and Goldreich 2003, 2004) Similarly, the Andoyer elements turn out to be non-osculating under angular-velocity-dependent perturbation. Amendment of only the Hamiltonian makes the equations render nonosculating elements. To make them osculating, more terms must enter the equations (and the equations will no longer be canonical). In practical calculations, is often convenient to deliberately deviate from osculation by substituting the Lagrange constraint with a condition that gives birth to a family of nonosculating elements. 
Michael Efroimsky