Le giunzioni incollate guadagnano sempre più mercato, nel campo delle costruzioni in genere, dove è necessario un alleggerimento delle strutture. Nel caso di geometrie di giunto semplici il dimensionamento avviene attraverso relazioni analitiche che restituiscono il valore massimo delle
tensioni, il quale deve essere inferiore al limite di utilizzo dell’adesivo stesso. Quando le geometrie sono complesse l’approccio analitico diventa impraticabile, di conseguenza si provvede a verificare la correttezza della soluzione mediante analisi agli elementi finiti (EF). L’introduzione del modello di zona coesiva nell'analisi EF permette di simulare il danneggiamento ed il cedimento del giunto in condizioni quasi-statiche e impulsive. In questo articolo si vuole implementare il modello per la simulazione della propagazione a fatica di difetti, utilizzando il software agli elementi finiti ABAQUS assieme a subroutine esterne interagenti con il modello EF stesso. Un punto focale dell'implementazione sarà il calcolo automatico del tasso di rilascio di energia G in modo indipendente dalla geometria del difetto stesso. I parametri del modello ricavati da prove di tenacità a frattura e propagazione di difetti a fatica in modo I, saranno utilizzati come riferimento per la convalida dell'implementazione.