In the present study, the experimental and finite element (FE) analyses are first carried out to investigate the deboning behavior of metal‐composite adhesive joints under modes of I and mode II loading. To conduct an FE on the debonding propagation, cohesive zone method (CZM), as well as maximum nominal stress and energy criteria, is applied. In the reliability analysis, to achieve the probability of debonding growth (PODG), limit state functions are formulated based on the energy release rate. To that end, the first‐order reliability method (FORM), the second‐order reliability method (SORM), and the Monte Carlo simulation (MCS) are used to calculate the PODG. The effect of initial debonding length on the PODG in all mentioned modes is investigated. Results obtained from reliability analysis reveal that the random variables including the initial debonding length, width, and thickness are the most sensitive variables to ascertain the PODG.