An angular analysis of the rare decay $B^0_s\rightarrow\phi\mu^+\mu^-$ is presented, using proton-proton collision data collected by the LHCb experiment at centre-of-mass energies of $7$, $8$ and $13~\rm{TeV}$, corresponding to an integrated luminosity of $8.4~\rm{fb}^{-1}$. The observables describing the angular distributions of the decay $B^0_s\rightarrow\phi\mu^+\mu^-$ are determined in regions of $q^2$, the square of the dimuon invariant mass. The results are consistent with Standard Model predictions.