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

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

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

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

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

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