Судя по всему, это более чистый вопрос о пиомо, чем тот, который имеет какое-либо отношение к CPLEX. Самое близкое, что я могу найти, это example4.py (в pyomo github). Этот пример pyomo имеет только квадратичную целевую функцию, но, похоже, будет просто изменить метод constraint_rule
, чтобы включить квадратичное ограничение.