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

, можно остаться в конфигурационном пространстве, вспомнить про метрику Якоби, геодезические которой есть траектории движения. Тогда уравнение связи должно определять геодезические подмногообразия чтобы интегрируемость точно сохранилась.
В общем, если уравнение связи инвариантно относительно потоков, порожденных первыми интегралами, то интегрируемость сохраняется.
В противном случае, видимо, нет. Есть, наверное, и другие варианты рассмотрения.
Это я рассмотрел случаи, когда связь не мешает движению первоначальной системы.
Однако, интегрируемость сохраняется, например, когда сначала точка свободно движется в

,
А после наложения связи

движется по единичной сфере, и импульсы

не находятся в инволюции с функцией

.
Так что всякие варианты здесь возможны.