Obtaining the Puiseux series of multiple imaginary (characteristic) roots (MIRs) is a fundamental issue in the stability analysis of time-delay systems. However, to the best of the authors' knowledge, this issue has not been fully investigated up to date. This note focuses on the Puiseux series expansion of MIRs of linear time-invariant systems including commensurate delays. For an MIR of any multiplicity, we propose an algorithm for defining the structure of the Puiseux series, as well as the explicit computation of the corresponding coefficients. By using the proposed method, we can find all the Puiseux series corresponding to all the root loci.