Abstract: Most techniques to estimate the service life of coatings are experimental in nature, cost expensive and are computationally heavy. Present disclosure provides systems and methods that predict the combined effects of crack path propagation and zones of delamination, that form on coating material and its surface due to weathering. The system of the present disclosure implemented a combined Finite Element Method (FEM) and Monte Carlo based simulation approach to capture the effects of delamination and crack propagation, respectively. The crack paths are predicted using a probabilistic model, considering crack propagation, branching, and keeping a record of crack age. Stress distribution computations are performed using FEM to understand stress concentration zones and delamination behavior with time, which is methodically also combined with the time sequence of cracking as well.