A numerical algorithm for ground reaction curves of circular openings in strain-softening rock masses satisfying the generalized Hoek-Brown criterion