A linear time-invariant dissipative Hamiltonian (DH) systemẋ = (J − R)Qx, with a skew-Hermitian J, an Hermitian positive semi-definite R, and an Hermitian positive definite Q, is always Lyapunov stable and under weak further conditions even asymptotically stable. In various applications there is uncertainty on the system matrices J, R, Q, and it is desirable to know whether the system remains asymptotically stable uniformly against all possible uncertainties within a given perturbation set. Such robust stability considerations motivate the concept of stability radius for DH systems, i.e., what is the maximal perturbation permissible to the coefficients J, R, Q, while preserving the asymptotic stability. We consider two stability radii, the unstructured one where J, R, Q are subject to unstructured perturbation, and the structured one where the perturbations preserve the DH structure. We employ characterizations for these radii that have been derived recently in [SIAM J. Matrix Anal. Appl., 37, pp. 2016 ] and propose new algorithms to compute these stability radii for large scale problems by tailoring subspace frameworks that are interpolatory and guaranteed to converge at a super-linear rate in theory. At every iteration, they first solve a reduced problem and then expand the subspaces in order to attain certain Hermite interpolation properties between the full and reduced problems. The reduced problems are solved by means of the adaptations of existing level-set algorithms for H∞-norm computation in the unstructured case, while, for the structured radii, we benefit from algorithms that approximate the objective eigenvalue function with a piece-wise quadratic global underestimator. The performance of the new approaches is illustrated with several examples including a system that arises from a finite-element modeling of an industrial disk brake.