Um modelo de dano baseado em energia com método de solução explícita no tempo utilizando o método dos elementos finitos é proposto para a simulação de delaminação induzida por fadiga de alto ciclo em compósito laminado de carbono/epoxi. Dois modelos utilizando zona coesiva para unir a mecânica do dano e mecânica da fratura são propostos, sendo que um é dependente de uma estimativa de comprimento da zona coesiva e outro é completamente independente deste comprimento. Ambos os modelos necessitam de parâmetros determinados experimentalmente, mais precisamente os coeficientes da equação que descreve a taxa de crescimento de trinca, comumente conhecida como equação de Paris. É mostrado que os modelos constitutivos apresentam resposta ligeiramente diferente entre si dependendo da condição de carregamento, evidenciando a utilização do comprimento da zona coesiva, mas quando utilizados para modelar ensaios estruturais reproduzem com qualidade a propagação para delaminação em modo puro e misto sem a necessidade de ajuste de dados.