Numerical modelling of wrinkled hyperelastic membranes with topologically complex internal boundary conditions