Equivalent linearization method is invented for nonlinear random vibration analysis. Gaussian mixture model is fitted to nonlinear dynamic response via optimization. GM-ELM can capture the non-Gaussianity of the nonlinear response effectively. GM-ELM can compute crossing rate, first-passage probability and mean peak response. Numerical examples show GM-ELM is highly accurate compared to the conventional ELM.