Application of Orthogonal Functions to Equivalent Linearization Method for MDOF Duffing-Van der Pol Systems under Nonstationary Random Excitations