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, using 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. It provides an alternative to the classical random vibration approach. It is also an alternative when the basic MCS becomes very inefficient.