The bin pack problem denotes the process of assigning a set of n items into a minimal number of bins of capacity c. It can be simply formulated as an ILP as per the below description:
My question is : how to express an additional constraint in the ILP formulation such that the load on the used bins is as much as possible balanced or fair across the used bins.
Thank you.

You could use minimizing the sum of absolute deviations approach.
Minimizing the sum of absolute deviations
Let deviations be represented by $$\epsilon_i = Y_i - \sum_{j} X_{ji}b_i,$$ where $i$ is the $i^{th}$ observation, $\epsilon_i$ gives the deviation, $Y_i$ is an observation. To minimize the deviation, the problem is formulated in a basic form as: $$\ Min \sum_{i} |\epsilon_i|$$ as the objective function, and linear constraints are $$\ s.t. \epsilon_i + \sum_{j}X_{ji}b_i=Y_i \text{ } \forall i, \ \epsilon_i, b_i \in \mathbb{R} $$
See primary source: https://optimization.mccormick.northwestern.edu/index.php/Optimization_with_absolute_values for details on how to linearize this.