AbstractThe failure modes of complex systems are often highly dependent. Aiming at the reliability analysis of engineering systems with multiple failure modes, an efficient surrogate model for system reliability with multiple failure modes is proposed, namely, the Copula Adaptive Kriging surrogate model (CAKSM). First, the initial Kriging model is established based on the initial samples. In order to select more effective added samples to update the Kriging model, a new learning function, namely, Copula Unites with Expected Improvement Function (CUEIF), is established. The function can make a reasonable sample selection by considering the failure correlation. Then, based on the optimal Copula function, the reliability data of each failure mode are coupled and analyzed. In order to adapt to various computing scales, regularization is introduced to improve the prediction speed of the surrogate model, and the probability density function of each performance function is obtained by introducing nonparametric kernel density estimation. Finally, the effectiveness of the CAKSM method and the learning function CUEIF is verified through four examples.