The restoration of electricity distribution networks determines the switching operations needed to isolate a faulted area, aiming to minimize operating costs related to the disconnection of loads and to the number of switching operations. Additionally, the final topology after the system restoration must ensure that the operational constraints are respected. In other words, the limits of nodal voltage and branch current magnitudes, the feeders capacities and the radial operation of the distribution networks must be respected. In this context, this work aims to develop modern convex models to solve the problem of restoration in electricity distribution networks. These models bring robustness, generality, and flexibility. Moreover, they are efficiently solved using commercial computational solvers, which makes them interesting for utilities.