Variational interface element model for 2D and 3D hydraulic fracturing simulations