Да, это интересная задача. Ведь это поможет понять, почему Луна движется подобным образом? Почему притормаживает и ускоряется?

9000 км - это слишком мало. Для 940 млн км земной орбиты, даже для 400 тыс км расстояния до Луны, 9000 - это ничто. Вклад этих вихляний в момент импульса Луны, полагаю, слишком мал, хотя, разумеется, должен быть. Согласен, однако, что надо учесть все факторы и все посчитать.

А вы не согласны, что осевое вращение Луны может быть неравномерным?

Иван Евгеньевич! Мне кажется, я понял, куда девается импульс. Я размышляю так:

Итак Луна движется относительно Земли не равномерно, "внутри" земной орбиты медленнее, а снаружи "быстрее". При этом она продолжает оставаться повернутой к Земле одной и той же стороной. Это значит, что на внутренней стороне орбиты Земли она должна вращаться медленнее, а на внешней - быстрее. Соответственно, если мы добавим к орбитальному импульсу осевой, он должен первый уравновесить. Как Вы считаете?

"Вихляние" Земли, я думаю, на импульс Луны оказывает мало влияния. Но с учетом вышесказанного, надо допустить, что угловая скорость Земли тоже неравномерна, хотя этого "практически" незаметно.

Тимофей Александрович, посмотрите Вот эту книгу, может найдете там полезное?

http://eqworld.ipmnet.ru/ru/library/mechanics/celestial.htm

Чеботарев Г.А. Аналитические и численные методы небесной механики. М.-Л.: Наука, 1965 (djvu)

ДА, конечно, на бумажке так и написано, переписал с ошибкой.

Итак, момент импульса гелиоцентрической луны зависит от времени, т.е. не постоянный.

А что если теперь учесть вращение Луны вокруг своей оси? Не поможет оно скомпенсировать переменность? Впрочем, вряд ли, оно ведь постоянное... А вот движение относительно орбиты Земли, вверх-вниз относительно плоскости эклиптики?

Или Вы хотите сказать, что есть две эклиптики - на эпоху J2000, и текущая, и их пересечением определяется линия узлов Земли?

Насколько я понимаю, ДВУ Земли служит отправной точкой для вычисления прецессии. Как и наклонение. ТВР постоянно смещается из-за прецессии, чтобы знать долготу ТВР нужна какая-то точка отсчета, за таковую и принимают ДВУ. Я только не понимаю физический смысл ДВУ Земли, ведь тут плоскость орбиты и плоскость эклиптики совпадают. Откуда же берется у орбиты Земли линия узлов?

Выходит ДВУ определяют по отношению к эклиптике эпохи J2000? 

Берут текущий наклон оси (obliquity) $\epsilon_0$ и прецессию $P_1$ за 1000 лет, и фиксируют эпоху J2000. 

Для вычисления вековых смещений (secular variations) вводят две переменных $sin \pi sin \Pi$ и $sin \pi cos \Pi$, где $\pi$ - наклонение эклиптики, а $\Pi$ - ДВУ.

Из этих данных вычисляют среднюю эклиптику на дату относительно фиксированной эклиптики J2000.

Далее считают производные относительно масс планет и относительно констант $\epsilon_0$ и $P_1$ и вычисляют погрешности.

Как-то так.

Недопонимание присутствует. Так как они определяют? Есть это в статье? Я не читал просто.

Поскольку ТВР смещается вдоль эклиптики в направлении орбитального движения Земли, а ДВУ измеряется от ТВР в ту же сторону, соответственно она должна уменьшаться.

Вас не затруднит разобраться, как они вычисляют долготу восходящего узла Земли? http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?db_key=AST&bibcode=1994A%26A...282..663S&letter=.&classic=YES&defaultprint=YES&whole_paper=YES&page=663&epage=663&send=Send+PDF&filetype=.pdf