Coupling functions treatment in a Bi-Level optimization process