In this paper, a suboptimal model reduction algorithm is presented for models with time delay. This method is an extension to the present authors’ method for models without a time delay. An objective function is defined to minimize the weighted integral squared value between the outputs of the reduced order model and the original model when they are subjected to the same input signal. Examples are given to show the advantages of using the reduced-order modelling technique provided.