Simulation algorithms for hierarchical Archimedean copulas beyond the completely monotone case