A novel reliability evaluation method for large nonlinear engineering systems excited by dynamic loading applied in time domain is presented. For this class of problems, the performance functions are expected to be function of time and implicit in nature. Available first-or second-order reliability method (FORM/SORM) will be challenging to estimate reliability of such systems. Because of its inefficiency, the classical Monte Carlo simulation (MCS) method also cannot be used for large nonlinear dynamic systems. In the proposed approach, only tens instead of hundreds or thousands of deterministic evaluations at intelligently selected points are used to extract the reliability information. A hybrid approach, consisting of the stochastic finite element method (SFEM) developed by the author and his research team using FORM, response surface method (RSM), an interpolation scheme, and advanced factorial schemes, is proposed. The method is clarified with the help of several numerical examples. (C) 2016 Faculty of Engineering, Ain Shams University. Production and hosting by Elsevier B.V.