This paper considers the optimal model reduction problem of matrix second-order linear systems in the sense of Hilbert-Schmidt-Hankel norm,with the reduced order systems preserving the structure of the original systems.The expressions of the error function and its gradient are derived.Two numerical examples are given to illustrate the presented model reduction technique.