A solution method for linear and geometrically nonlinear MDOF systems with random properties subject to random excitation