Εξ’ αιτίας της ραγδαίας αύξησης στη χρήση σύνθετων υλικών, στην παρούσα εργασία ασχοληθήκαμε με σύνθετο υλικό μήτρας εποξικής ρητίνης και ενίσχυσης κόκκους αλουμινίου με μέση διάμετρο 63 μm. Σκοπός της έρευνας αυτής είναι κατά πρώτον να μελετήσει τη συμπεριφορά του συνθέτου σε στατική κάμψη αναλόγως της κατά βάρος περιεκτικότητά του σε κόκκους και να συγκρίνει τα πειραματικά αποτελέσματα με αντίστοιχες θεωρητικές προβλέψεις και κατά δεύτερον να μελετήσει την επίδραση βλάβης λόγω ύπαρξης εγκοπής ή θερμικής κόπωσης στην υποβάθμιση του σύνθετου υλικού και να προβλέψει την παραμένουσα μηχανική συμπεριφορά χρησιμοποιώντας το μοντέλο R.P.M. Τρία διαφορετικά μοντέλα πρόβλεψης, ανεπτυγμένα απ’ το CMG (Composite Material Group) του Πανεπιστημίου Πατρών, θα εφαρμοστούν για τις προβλέψεις και δύο απ’ αυτά θα συγκριθούν και με άλλα μοντέλα ανεπτυγμένα από άλλους ερευνητές. Κατ’ αρχήν, για τον υπολογισμό της μέγιστης τάσης εφαρμόστηκε το μοντέλο του τμηματικού κόκκου. Σύμφωνα με το μοντέλο αυτό, ο κόκκος χωρίζεται σε άπειρο αριθμό ινών και σύμφωνα με την θεωρία του Cox μπορεί να υπολογιστεί η μέση τάση που ασκείται σε κάθε ίνα του κόκκου. Στη συνέχεια, για τον υπολογισμό του μέτρου ελαστικότητας συναρτήσει της κατά βάρους περιεκτικότητας σε κόκκους αλουμινίου εφαρμόστηκε το μοντέλο της μεσοφάσης. Το μοντέλο αυτό λαμβάνει υπόψιν την ύπαρξη μιας τρίτης φάσης ανάμεσα στις δύο κύριες (μήτρα-ενίσχυση) και η οποία αποτελεί μια σημαντική παράμετρο επίδρασης στη συμπεριφορά οποιουδήποτε σύνθετου υλικού. Τέλος, στις περιπτώσεις της ύπαρξης εγκοπής και της θερμικής κόπωσης εφαρμόστηκε το μοντέλο R.P.M. Το μοντέλο αυτό μπορεί να εφαρμοστεί για την περιγραφή της υποβάθμισης των ιδιοτήτων των υλικών μετά από βλάβη και στη συγκεκριμένη εργασία εφαρμόστηκε για να περιγράψει τις μεταβολές στην μέγιστη τάση του υλικού και στο μέτρο ελαστικότητάς του συναρτήσει του μήκους εγκοπής ή του αριθμού των θερμικών κύκλων. / -

Στα πλαίσια της παρούσας διατριβής αναπτύχθηκε ένα τρισδιάστατο παραµετρικό µοντέλο Προοδευτικής Βλάβης (ΠΒ) για την προσοµοίωση της µηχανικής συµπεριφοράς των σύνθετων υλικών σε εφελκυσµό και κόπωση. Το µοντέλο αναπτύχθηκε µε τέτοιο τρόπο, ώστε να µπορεί να εφαρµοστεί µε µικρές τροποποιήσεις σε οποιοδήποτε κατασκευαστικό στοιχείο από σύνθετα υλικά. Αποτελείται από τις συνιστώσες της τρισδιάστατης ανάλυσης τάσεων, η οποία πραγµατοποιείται µε την µέθοδο των πεπερασµένων στοιχείων, της ανάλυσης αστοχίας, η οποία πραγµατοποιείται µε πολυωνυµικά κριτήρια αστοχίας και της υποβάθµισης των ιδιοτήτων του σύνθετου υλικού, η οποία πραγµατοποιείται µε χρήση κανόνων υποβάθµισης. Οι συνιστώσες αυτές, λειτουργούν µε βάση έναν επαναληπτικό αλγόριθµο, ο οποίος σταµατά όταν επαληθευτεί ένα προκαθορισµένο κριτήριο τελικής αστοχίας. Στο σηµείο αυτό, καθορίζεται η αποµένουσα αντοχή του κατασκευαστικού στοιχείου. Όλες οι συνιστώσες του µοντέλου προγραµµατίστηκαν στον κώδικα πεπερασµένων στοιχείων ANSYS, δηµιουργώντας µια εύχρηστη µακρο-ρουτίνα. Το µοντέλο ΠΒ έχει την δυνατότητα να εκτιµά την έναρξη και διάδοση βλάβης συναρτήσει του φορτίου, την δυσκαµψία και την αποµένουσα αντοχή του κατασκευαστικού στοιχείου στο οποίο εφαρµόζεται. Εφαρµογή του µοντέλου ΠΒ έγινε στο πρόβληµα των µηχανικών συνδέσεων πολύστρωτων πλακών υπό εφελκυστικά φορτία, το οποίο δεν έχει αντιµετωπιστεί επιτυχώς θεωρητικά µέχρι σήµερα, κυρίως λόγω της αδυναµίας υπολογισµού του τρισδιάστατου πολύπλοκου τασικού πεδίου, που αναπτύσσεται γύρω από τον ήλο και των πολύπλοκων µηχανισµών αστοχίας που αναπτύσσονται στο σύνθετο υλικό. Η µηχανική σύνδεση που µοντελοποιήθηκε, αποτελούνταν από µια πολύστρωτη πλάκα, µια πλάκα από αλουµίνιο και έναν ήλο µε προεξέχον κεφάλι. Για τον υπολογισµό των τάσεων αναπτύχθηκε ένα τρισδιάστατο µοντέλο πεπερασµένων στοιχείων της σύνδεσης στον 102 ΚΕΦΑΛΑΙΟ 10: Περίληψη-Συµπεράσµατα • • • κώδικα ANSYS. Ιδιαίτερη έµφαση δόθηκε στην λεπτοµερή µοντελοποίηση της περιοχής γύρω από τον ήλο. Όλες οι πιθανές επαφές µεταξύ των µελών της σύνδεσης µοντελοποιήθηκαν. Η σύσφιξη της σύνδεσης προσοµοιώθηκε ως ένα αρχικό-ξεχωριστό βήµα φόρτισης. Το µοντέλο πεπερασµένων στοιχείων επαληθεύτηκε εκτενώς, µέσω συγκρίσεων διαφορετικών τασικών κατανοµών µε αντίστοιχες αριθµητικές και αναλυτικές λύσεις. Η εφαρµογή του µοντέλου ΠΒ στην µηχανική σύνδεση, αποτελεί συγχρόνως και έλεγχο των δυνατοτήτων του. Τα αποτελέσµατα της εφαρµογής είναι: η εκτίµηση της έναρξης και διάδοσης βλάβης στην πολύστρωτη πλάκα συναρτήσει του φορτίου, ο µακροσκοπικός µηχανισµός αστοχίας, η δυσκαµψία και η αποµένουσα αντοχή της σύνδεσης. Για την εκτίµηση της επίδρασης της γεωµετρίας της σύνδεσης στην διάδοση βλάβης και στον µακροσκοπικό µηχανισµό αστοχίας της σύνδεσης, πραγµατοποιήθηκε µια παραµετρική ανάλυση, στην οποία µεταβάλλονταν οι γεωµετρικοί λόγοι e/d και w/d που καθορίζουν την θέση του ήλου. Παρατηρήθηκε ότι: στους µικρούς λόγους w/d (ο ήλος είναι κοντά στο πλαϊνό άκρο της πλάκας) ο µακροσκοπικός µηχανισµός αστοχίας ήταν ‘εφελκυσµός’, ενώ στους µεγάλους ‘εισχώρηση’, στους µικρούς λόγους e/d (ο ήλος είναι κοντά στο κάτω άκρο της πλάκας) ο µηχανισµός ήταν ‘διάτµηση’ και για µεγάλους ‘εισχώρηση’, ενώ σε δύο γεωµετρίες, παρατηρήθηκαν µικτές αστοχίες (‘διάτµηση’-‘εισχώρηση’ και ‘εφελκυσµός’-‘εισχώρηση’). Οι µηχανισµοί αστοχίας της ‘διάτµησης’ και του ‘εφελκυσµού’ οδήγησαν σε µικρότερες αντοχές των συνδέσεων σε σχέση µε την ‘εισχώρηση’. Όλες οι εκτιµήσεις που αφορούν την έναρξη και διάδοση βλάβης, καθώς και τον µακροσκοπικό µηχανισµό αστοχίας των συνδέσεων, επαληθεύτηκαν από αντίστοιχα αριθµητικά και πειραµατικά αποτελέσµατα που αντλήθηκαν από την βιβλιογραφία. Η καταλληλότητα των κριτηρίων αστοχίας και των κανόνων υποβάθµισης ιδιοτήτων, καθώς και η επίδραση τους στα αποτελέσµατα του µοντέλου, µελετήθηκαν για πρώτη φορά σε ένα µοντέλο ΠΒ. Η µελέτη πραγµατοποιήθηκε, µέσω σύγκρισης της εκτίµησης της δυσκαµψίας και αποµένουσας αντοχής συνδέσεων διαφορετικής γεωµετρίας και αλληλουχίας στρώσεων, µε αντίστοιχα πειραµατικά αποτελέσµατα από την βιβλιογραφία. Στα πλαίσια αυτής της µελέτης, προτάθηκαν δύο νέες οµάδες κριτηρίων αστοχίας και κανόνων υποβάθµισης η καταλληλότητα των οποίων, συγκρίθηκε µε αυτήν των οµάδων που χρησιµοποιήθηκαν αρχικά στο µοντέλο. Τα αποτελέσµατα της µελέτης συνοψίζονται ως εξής: 103 ΚΕΦΑΛΑΙΟ 10: Περίληψη-Συµπεράσµατα • • • • • η εκτίµηση της δυσκαµψίας της σύνδεσης, ήταν ικανοποιητική σε όλες τις περιπτώσεις και δεν επηρεάστηκε από τον εκάστοτε συνδυασµό των κριτηρίων αστοχίας και κανόνων υποβάθµισης, η εκτίµηση της αποµένουσας αντοχής της σύνδεσης, επηρεάστηκε σε µεγάλο βαθµό από τον εκάστοτε συνδυασµό των κριτηρίων αστοχίας και κανόνων υποβάθµισης. Συγκεκριµένα, η ανάλυση µε τον συνδυασµό των αρχικών οµάδων κριτηρίων αστοχίας και κανόνων υποβάθµισης, έδωσε αρκετά συντηρητικές προβλέψεις σε όλες τις γεωµετρίες της σύνδεσης (απόκλιση από τα πειράµατα 35-40%), ενώ η ανάλυση µε τον συνδυασµό των νέων οµάδων, έδωσε πολύ ικανοποιητικές εκτιµήσεις της αποµένουσας αντοχής που εµφάνισαν απόκλιση από τα πειράµατα η οποία κυµαινόταν από 1.6 µέχρι 6%. Στην συνέχεια, το µοντέλο ΠΒ επεκτάθηκε-τροποποιήθηκε, ώστε να δύναται να εφαρµοστεί σε περιπτώσεις φόρτισης κόπωσης (µοντέλο Προοδευτικής Βλάβης Κόπωσης (ΠΒΚ)). Η βασική διαφορά του µοντέλου ΠΒΚ µε το µοντέλο ΠΒ έγκειται στην προσοµοίωση της κυκλικής φόρτισης και την προσθήκη στην υποβάθµιση του σύνθετου υλικό λόγω τοπικών αστοχιών, της υποβάθµισης λόγω κόπωσης, η οποία οφείλεται στην φύση της κυκλικής φόρτισης. Για την προσοµοίωση της κυκλικής φόρτισης αναπτύχθηκε στην παρούσα εργασία µια νέα µεθοδολογία, η οποία: η εφαρµογή της οποίας απαιτεί την εκτέλεση µικρού αριθµού βηµάτων φόρτισης, λαµβάνει υπόψη τον εκάστοτε λόγο τάσεων R της φόρτισης, και µπορεί να εφαρµοστεί για κυκλική φόρτιση οποιουδήποτε λόγου τάσεων. Η υποβάθµιση λόγω κόπωσης εφαρµόζεται σε επίπεδο δυσκαµψίας και αντοχής και είναι συνάρτηση των κύκλων φόρτισης. Για την προσοµοίωση της στην παρούσα εργασία, αναπτύχθηκε µια µεθοδολογία, η οποία απαιτεί για να εφαρµοστεί µικρό αριθµό πειραµάτων. Το µοντέλο ΠΒΚ είναι επίσης παραµετρικό, αφού µπορεί να εφαρµοστεί σε οποιοδήποτε κατασκευαστικό στοιχείο υποβαλλόµενο σε κόπωση σταθερού εύρους οποιουδήποτε λόγου τάσεων, απαιτώντας δεδοµένα από µικρό αριθµό πειραµάτων, µέσα από τα οποία θα γίνει η προσοµοίωση της υποβάθµισης λόγω κόπωσης του συγκεκριµένου υλικού. Το µοντέλο ΠΒΚ έχει την δυνατότητα να εκτιµά την διάδοση βλάβης κόπωσης συναρτήσει του αριθµού των κύκλων και την διάρκεια ζωής της κατασκευαστικού στοιχείου στο οποίο εφαρµόζεται. Εφαρµογή του µοντέλου ΠΒΚ έγινε σε πολυµερή πολύστρωτες πλάκες από δύο διαφορετικά σύνθετα υλικά (Fiberdux-HTA/6376 και APC-2), οι οποίες υποβάλλονται σε κόπωση εφελκυσµού-θλίψης (R=-1). Για την προσοµοίωση της υποβάθµισης λόγω κόπωσης και την επαλήθευση των τελικών αποτελεσµάτων του µοντέλου, 104 ΚΕΦΑΛΑΙΟ 10: Περίληψη-Συµπεράσµατα χρησιµοποιήθηκαν πειραµατικά αποτελέσµατα που υπήρχαν στο εργαστήριο Τεχνολογίας & Αντοχής Υλικών. Παρουσιάστηκε η εκτίµηση του µοντέλου για την διάδοση της διαστρωµατικής αποκόλλησης και της αστοχίας των ινών σε µια περίπτωση χαµηλής και µια υψηλής µέγιστης τάσης κόπωσης, αντίστοιχα. Η διάδοση αυτών των µηχανισµών αστοχίας επιλέχτηκε διότι ο βαθµός συσσώρευσης τους χρησιµοποιήθηκε από το µοντέλο ως κριτήριο τελικής αστοχίας. Επίσης, η εκτίµηση της διάδοσης της διαστρωµατικής αποκόλλησης, για διαφορετικά ποσοστά της διάρκειας ζωής µιας Fiberdux-HTA/6376 πολύστρωτης πλάκας, συγκρίθηκε µε αντίστοιχά πειραµατικά αποτελέσµατα (C-scan γραφήµατα). Η σύγκριση έδειξε καλή συµφωνία, όσον αφορά την διαστρωµατική αποκόλληση που αναπτύχθηκε στα ελεύθερα άκρα των πλακών. Και στις δύο περιπτώσεις, η διαστρωµατική αποκόλληση αναπτύχθηκε στα ελεύθερα άκρα λόγω υψηλών διαστρωµατικών τάσεων στους πρώτους κύκλους φόρτισης και διαδόθηκε προς το εσωτερικό της πλάκας. Ωστόσο, η διάδοση της διαστρωµατικής αποκόλλησης που αναπτύχθηκε στις εσωτερικές περιοχές αρχικής βλάβης, δεν ήταν δυνατόν να προβλεφτεί, εξαιτίας του ότι στο µοντέλο ΠΒΚ δεν προσοµοιώθηκε η αρχική βλάβη. Τέλος, µε βάση την διάδοση βλάβης και την χρήση κατάλληλων κριτηρίων τελικής αστοχίας, εκτιµήθηκε, για διαφορετικά επίπεδα µέγιστης τάσης, η διάρκεια ζωής σε κόπωση εφελκυσµού-θλίψης (R=-1) των δύο πολύστρωτων πλακών και κατασκευάστηκαν οι αντίστοιχες καµπύλες S-N. Η σύγκριση µεταξύ των εκτιµηθέντων και πειραµατικών καµπυλών S-N, έδειξε µια αρκετά ικανοποιητική συσχέτιση, η οποία ήταν όµως διαφορετική για τα δύο υλικά. Η διαφοροποίηση αυτή οφείλεται στο συνδυασµό της πειραµατικής διασποράς µε το σφάλµα που επεισέρχεται στο µοντέλο από τον τρόπο προσοµοίωσης της υποβάθµισης της δυσκαµψίας και αντοχής λόγω κόπωσης (προσαρµογή των πειραµατικών σηµείων σε κάθε επίπεδο τάσης). / In this thesis a three-dimensional progressive damage model was developed to simulate the mechanical behaviour of composite structural components under tesnile and fatigue loading. First the model was used to simulate the mechanical behaviour of composite bolted joints under tensile loading. The model was able to simulate the initation and propagation of damage with increased load and to assess the macroscopic failure mode, stiffness and residual strenght of the joint. The numerical predictions were verified with experimental results. The model was then used to simulate the mechanical behaviour of CFRP laminates in tension-compression fatigue. The model was able to simulate the progression of fatigue damage and to evaluate the fatigue life of CFRP plates. The different model components were implemented into the ANSYS FE CODE using a parametric macro-routine.

Τα δομικά στοιχεία των μηχανολογικών κατασκευών υπόκεινται σε επαναλαμβανόμενες κυκλικές καταπονήσεις, από τις οποίες δημιουργούνται και διαδίδονται ρωγμές. Οι καταπονήσεις αυτές, οι οποίες προκαλούν κόπωση στις κατασκευές, μπορεί να είναι είτε καθαρά μηχανικές είτε θερμικές ή να προκύπτουν σα συνδυασμός θερμικής και μηχανικής φόρτισης. Τυπικές περιπτώσεις θερμικών και θερμομηχανικών φορτίσεων εμφανίζονται σε κατασκευές, όπως σωλήνες κυκλωμάτων ψύξης, πιεστικά δοχεία, συνιστώσες ηλεκτρικών κυκλωμάτων, θάλαμοι μηχανών εσωτερικής καύσης και πτερύγια στροβιλοκινητήρων. Η κυκλική μεταβολή του θερμικού φορτίου στις προαναφερθείσες περιπτώσεις, συνιστά συνθήκες θερμικής κόπωσης. Επίσης, λόγω της σχετικά υψηλής συχνότητας του φορτίου η θερμοκρασία παρουσιάζει έντονη μεταβολή στο χώρο και στο χρόνο. Ο προσδιορισμός της διάρκειας ζωής ενός δομικού στοιχείου κατά τη φάση του σχεδιασμού μπορεί να γίνει με τη βοήθεια πειραματικών διαδικασιών. Τα πειράματα όμως κόπωσης είναι δαπανηρά και χρονοβόρα και προφανώς απαιτούνται περισσότερες από μια πειραματικές δοκιμές. Οπότε, είναι εύλογο να υπάρχουν υπολογιστικά εργαλεία που να δίνουν τη δυνατότητα στο μηχανικό να εκτιμήσει την διάρκεια ζωής ή τη σοβαρότητα της βλάβης ενός εξαρτήματος. Τα περισσότερα υπολογιστικά μοντέλα αναφέρονται σε καθαρά μηχανικές καταπονήσεις. Έτσι υπάρχει πρόσφορο έδαφος για την ανάπτυξη υπολογιστικών εργαλείων για την ανάλυση προβλημάτων θερμικής και θερμομηχανικής κόπωσης. Τέτοιου είδους εργαλεία θα πρέπει να λαμβάνουν υπόψη το κλείσιμο των ρωγμών, που συμβαίνει λόγω των θερμικών παραμορφώσεων, διότι είναι δυνατόν να επηρεάζεται τοπικά το θερμοκρασιακό πεδίο. Επομένως, χρειάζεται επαναληπτική διαδικασία για τον προσδιορισμό του θερμικού και τασικού πεδίου που αλληλεπιδρούν. Είναι προφανές ότι η ανάλυση της θερμικής κόπωσης εξελίσσεται σε συνθέτη διαδικασία, που θα πρέπει να συμπεριλαμβάνει τον υπολογισμό της κατανομής της θερμοκρασίας, την τοπική επίδραση του άκρου της ρωγμής στο τασικό πεδίο καθώς και την επαφή των επιφανειών της ρωγμής. Η μέθοδος των συνοριακών στοιχείων είναι ικανή να αντιμετωπίζει τέτοιου είδους τοπικές επιδράσεις. Η παρούσα διατριβή επικεντρώνεται στην ανάπτυξη υπολογιστικού εργαλείου βασισμένου στα συνοριακά στοιχεία, για την πρόβλεψη της διάδοσης ρωγμών και την εκτίμηση της διάρκειας ζωής, εξαρτημάτων υπό θερμική και θερμομηχανική κόπωση. Έμφαση δίνεται σε περιπτώσεις που το θερμικό φορτίο προκαλεί κλείσιμο της ρωγμής και σε περιπτώσεις διεπιφανειακών ρωγμών, όπου το θερμοκρασιακό πεδίο επηρεάζεται από την θερμική αντίσταση ανάμεσα στις επιφάνειες της ρωγμής. Στο πρώτο κεφάλαιο γίνεται βιβλιογραφική ανασκόπηση σε εργασίες που εστιάζουν σε φαινόμενα κόπωσης και διάδοσης ρωγμών, καθώς και στην ανάπτυξη υπολογιστικών μοντέλων για την πρόβλεψη της διάδοσης ρωγμών. Επιπλέον, προσδιορίζεται λεπτομερώς το αντικείμενο της παρούσας διατριβής και εξηγείται η συνεισφορά της και τα καινοτόμα σημεία της. Στο δεύτερο κεφάλαιο περιγράφεται η ιδιόμορφη συμπεριφορά του άκρου της ρωγμής, δίνονται οι διατυπώσεις των μεγεθών θραύσης που χρησιμοποιούνται στην ανάλυση της κόπωσης και αναφέρονται τρόποι με τους οποίους μελετάται η διάδοση ρωγμών. Στο τρίτο κεφάλαιο περιγράφονται λεπτομερώς οι ολοκληρωτικές συνοριακές διατυπώσεις για την επίλυση προβλημάτων θερμοελαστικότητας. Στο τέταρτο κεφάλαιο περιγράφονται οι υπολογιστικές διαδικασίες που ακολουθούνται στην παρούσα εργασία για τον προσδιορισμό του πεδίου θερμοκρασιών και μετατοπίσεων, καθώς και ο τρόπος που προσομοιώνεται η διάδοση ρωγμής. Στο πέμπτο κεφάλαιο παρατίθενται τα αποτελέσματα που προέκυψαν από τις αναλύσεις για διάφορες περιπτώσεις, ενώ στο έκτο κεφάλαιο εξάγονται συμπεράσματα και διατυπώνονται προτάσεις για μελλοντική έρευνα. / The prediction of fatigue life is essential for the integrity and reliability of a structure when designing engineering components that undergo cyclic loading. In most cases, the mechanical cyclic loads are taken into account in order to evaluate the life and damage tolerance of structures with existing cracks. However, there exists a category of structures that experience severe thermal cycling that acts alongside the mechanical loads. Such structures include cooling system pipes, pressure vessels, pistons and combustion chambers of internal combustion engines, gas turbine blades and components of electrical circuits. Interfacial crack growth is of paramount importance when designing components that are protected by thermal barrier coatings in order to increase their endurance and efficiency. These types of structures are exposed to very intense thermo-mechanical cycling, which gradually causes delamination and eventually leads to spallation of the coating Numerical simulations, via the finite element method, are a common trend, when analysing the endurance of coated components. However, important aspects such as the heat exchange between the contacting faces and friction are not taken into account in fracture assessments of these components. The boundary element method is very attractive for crack-growth analyses because only the boundary is meshed, rather than the whole domain of the problem. In the present thesis, the boundary integral equations of uncoupled, time-dependent thermo-elasticity are employed to account for the time-varying nature of the thermal load. Our study discusses the influence of crack closure on quasi-static, sub-critical crack extension in the presence of thermo-mechanical cyclic loading. Appropriate thermal and mechanical boundary conditions are imposed on the numerical model to account for the contact state. The validity of the code to compute the temperature distribution under thermal cycling is checked through analytical solutions. Afterwards, a pure mode-I and mixed mode fracture problems in homogeneous material are analysed and the results are compared to other boundary element solutions. The singularity resulting from tractions and heat flux around the crack tip is effectively captured by singular quarter-point elements, while the fracture magnitudes can be computed using appropriate traction formulas. In these problems, the fatigue life is evaluated in terms of load cycle when the crack closure is considered. The number of cycles required for an existing crack to grow a certain length can be empirically predicted using the Paris’ law. The crack extension angle is evaluated by means of the maximum circumferential stress. The results are discussed, clearly indicating the impact of crack closure on fatigue life evaluation. The main conclusion is that crack closure should be incorporated into the analysis whenever the contact effect is inevitable. Otherwise, the fatigue life may be underestimated, leading to a conservative design. Finally, the sub-domain boundary element procedure is applied to interfacial cracks where the crack closure is more pronounced. Specifically, a case of a thermal barrier coating system is investigated. The thermal resistance between the contacting crack faces is incorporated into the procedure and it is assumed to be dependent on the contact pressure. If crack closure due to thermal distortion takes place, then the displacement and traction field may affect the heat flux between the crack faces, and the thermal and mechanical parts of the problem will need to be solved repeatedly until thermo-mechanical convergence is achieved. The results suggest that there are significant effects on the behaviour of stably growing cracks and the evaluation of failure capacity, emanating from crack closure, the amount of thermal resistance and the phase angle between the mechanical and thermal loads.