Στην παρούσα διπλωματική εργασία αναπτύχθηκε τρισδιάστατο αριθμητικό μοντέλο με βάση την μέθοδο των πεπερασμένων στοιχείων για την προσομοίωση της μηχανικής συμπεριφοράς συνδέσεων με κόλλα πολύστρωτων πλακών. Το μοντέλο αναπτύχθηκε χρησιμοποιώντας το εμπορικό πακέτο πεπερασμένων στοιχείων ANSYS. Για την προσομοίωση της συμπεριφοράς της κόλλας χρησιμοποιήθηκαν δύο προσεγγίσεις: η μοντελοποίηση της ζώνης συνοχής και η μοντελοποίηση της βλάβης του συνεχούς μέσου. Οι δύο αυτές προσεγγίσεις συγκρίθηκαν τόσο ως προς την αξιοπιστία τους, η οποία καθορίζεται από την σύγκριση με πειραματικά αποτελέσματα, όσο και ως προς την ευκολία εφαρμογής τους, η οποία καθορίζεται από τα δεδομένα που απαιτούνται και τον υπολογιστικό χρόνο. Η σύγκριση των δύο μεθοδολογιών έγινε στην βάση της εφαρμογής τους για την προσομοίωση της μηχανικής συμπεριφοράς σε φόρτιση τύπου Ι σύνδεσης με κόλλα μεταξύ δύο ψευδοισότροπων CFRP πολύστρωτων πλακών. Το συγκεκριμένο πρόβλημα επελέγη διότι υπήρχαν διαθέσιμα πειραματικά αποτελέσματα προς σύγκριση στο Εργαστήριο. Οι πολύστρωτες πλάκες μοντελοποιήθηκαν χρησιμοποιώντας το στρωματικό στοιχείο του ANSYS SOLID185. Στο στοιχείο αυτό κάθε στρώση μοντελοποιείται ξεχωριστά ως ορθότροπο υλικό. Η εφαρμογή της μοντελοποίησης της ζώνης συνοχής έγινε μέσω της χρήσης του στοιχείου του ANSYS INTER205. Για την εφαρμογή της μοντελοποίησης της βλάβης του συνεχούς μέσου αναπτύχθηκε μακρο-ρουτίνα χρησιμοποιώντας την γλώσσα προγραμματισμού του κώδικα ANSYS. Τα αριθμητικά αποτελέσματα έδειξαν ότι και οι δύο μεθοδολογίες προσομοιώνουν με ικανοποιητική ακρίβεια την καμπύλη δύναμης-μετατόπισης της σύνδεσης. Σχετικά με την ευκολία εφαρμογής των δύο μεθόδων, η σύγκριση έδειξε ότι η μέθοδος της μοντελοποίησης της ζώνης συνοχής υπερτερεί έναντι της μεθόδου μοντελοποίησης της βλάβης του συνεχούς μέσου διότι απαιτεί μικρότερο αριθμό δεδομένων, μειονεκτεί όμως ως προς τον απαιτούμενο υπολογιστικό χρόνο. Και οι δύο μέθοδοι κρίνονται κατάλληλες για χρήση στην αριθμητική σχεδίαση συνδέσεων με κόλλα. / This work is based on the development of three-dimensional numerical model based on the finite element method to simulate the mechanical behavior of adhesive bonded joints in composite materials. The model was developed in finite element procedures under the framework of the commercial software ANSYS. To simulate the behavior of the adhesive used two approaches: the cohesive zone modeling (CZM) and the continuum damage modeling (CDM). These two approaches are compared both in terms of reliability, which is determined by comparison with experimental results, and applicability, which is determined by the parameters required and the computational time. The comparison between the two methodologies was the basis of their application to simulate the mechanical behavior under mode-I fracture behavior of adhesively bonded joints between two CFRP plates. This problem was chosen because there were experimental results to compare in the laboratory. The sandwich plates are modeled using the stromal element of ANSYS SOLID185. This item each layer separately modeled as orthotropic material. The adhesive is modeled using the interface element of ANSYS INTER205. For the purpose of modeling the failure of continuous medium developed macro routine using the programming language code ANSYS. The numerical results showed that both methodologies simulate sufficient precision the curve force-displacement of the connection. About applicability of the two methods, the comparison showed that the process of cohesive zone modeling outweighs the process of continuum damage modeling because it requires less number of parameters, but falls to the computational time. Both methods are suitable for use in numerical design of adhesively bonded joints.