The model order reduction of second-order dynamical multi-input and multi-output (MIMO) systems with proportional damping arising in the numerical simulation of mechanical structures is discussed. Based on finite element modeling the systems describing the mechanical structures are large and sparse, either undamped or proportionally damped. This work concentrates on a new model reduction algorithm for such second order MIMO systems which automatically generates a reduced system approximating the transfer function in the lower range of frequencies. The method is based on the rational global Arnoldi method. It determines the expansion points iteratively. The reduced order and the number of moments matched per expansion point are determined adaptively using a heuristic based on some error estimation. Numerical examples comparing our results to modal reduction and reduction via the rational block Arnoldi method are presented. Key words: Model Order Reduction, Simulation, Krylov Subspace,...