смысл метода лагранжевых множителей состоит в следующем. Пусть для простоты есть конечномерное пространство

, и мы хотим найти экстремум функции

на подпространстве

, заданном условиями

. Для этого нам надо занулить вариацию функции
вдоль подпространства

. При этом, вообще говоря, вариации в трансверсальных направлениях не будут равны нулю. Вместо того, чтобы разрешать явно условия на

и варьировать по независимым координатам, можно добавить к функции множители Лагранжа. Тогда уравнение будет выглядеть как

, а также уравнения связей. Смысл суммы с лямбдами в вариации тот, что она параметризует (ко)нормальное расслоение к

, т.е. вариация

в трансверсальных направлениях к

не обязана быть равной нулю, как мы и хотели.
Соответственно, в вашем лагранжиане связь никуда не девается. Поле лямбда нельзя просто так исключить (а подставлять его в лагранжиан бессмысленно, т.к. оно умножается на нечто, что равно нулю). Если вы его выкинули, то будьте добры не забывать его уравнения движения, т.е. связь, и допускать вариации

только вдоль нее.