Your psglength constraints have the form (sum of products of variables) = constant. The left side is a quadratic function, and is not positive semidefinite. CPLEX only handles quadratic constraints of the form (positive semidefinite express <= RHS) or what are called second order cone constraints. Your constraint is neither, and in fact makes the feasible region non-convex.
The CPLEX user manual has a section on what sorts of quadratic constraints are allowed.
Paul
#CPLEXOptimizers#DecisionOptimization