Στην παρούσα διπλωματική εργασία μελετάται η συμπεριφορά του ανθρώπινου δέρματος κατά την εφαρμογή σημειακής (κατακόρυφης) δύναμης αλλά και μη σημειακής (επιφανειακής). Για τη μελέτη κατασκευάζεται εικονικό μοντέλο τμήματος του δέρματος με χρήση πεπερασμένων στοιχείων. Στόχος είναι η εισαγωγή του μαθηματικού μοντέλου σε ρομποτικό σύστημα το οποίο θα κάνει ένεση προκαλώντας το λιγότερο δυνατό πόνο ή μασάζ με τη μεγαλύτερη δυνατή αποτελεσματικότητα. Στη διπλωματική δεν προγραμματίζεται κάποιο ρομπότ που να επιτελεί το έργο αυτό, αλλά αναπτύσσεται ο κώδικας σε matlab ο οποίος μετασχηματισμένος θα μπορούσε να χρησιμοποιηθεί από ένα ρομποτικό σύστημα. Αρχικά, αναλύεται η μέθοδος των πεπερασμένων στοιχείων ώστε να επιλεγεί το κατάλληλο μοντέλο για το πρόβλημά μας. Με την παραδοχή ότι υπάρχει αξονική συμμετρία επιλέγεται να χρησιμοποιηθεί μια δομή δύο διαστάσεων αποτελούμενη από δύο στρώματα διαφορετικής ελαστικότητας (υποθετικά, δέρμα και υποδόριος ιστός). Για τη μοντελοποίηση αυτή αναλύονται οι διαδικασίες που πρέπει να ακολουθηθούν και οι παραδοχές που πρέπει να ληφθούν υπ’ όψιν. Επίσης, αναλύεται με μαθηματικό τρόπο η προσέγγιση του προβλήματος παραμόρφωσης σε δύο διαστάσεις. Στη φάση αυτή, ανάπτυξης και επιλογής του βέλτιστου μοντέλου, η προσομοίωση γίνεται εξ' ολοκλήρου σε matlab. Στη συνέχεια, αφού ορισθεί το μοντέλο σε δύο διαστάσεις πρέπει να προσδιοριστούν οι παράμετροι για το σχετικό πάχος των δύο στρωμάτων και τη σχετική τους ελαστικότητα ώστε το μοντέλο να ανταποκρίνεται στην πραγματικότητα. Για να επιτευχθεί η συγκριτική αυτή αξιολόγηση, χρησιμοποιούνται δύο ομάδες πραγματικών πειραματικών τιμών που περιγράφονται σε σχετικές μελέτες διαθέσιμες στη διεθνή βιβλιογραφία. Δοκιμάζονται συγκριτικά διάφορα ζεύγη παραμέτρων με στόχο να γίνει μια ανάλυση ευαισθησίας του μοντέλου, και να εντοπιστεί αυτός ο συνδυασμός τιμών για τον οποίο υπάρχει η ελάχιστη απόκλιση από τα πειραματικά δεδομένα. Τέλος, αξίζει να αναφέρουμε πως στο δεύτερο τμήμα της εργασίας όπου μοντελοποιείται η συμπεριφορά του δέρματος σε μη σημειακή καταπόνηση χρησιμοποιούνται τα ίδια πειραματικά δεδομένα καθώς και δεδομένα από άλλες προσπάθειες μοντελοποίησης. Ο κώδικας που χρησιμοποιείται είναι κατά ένα μεγάλο ποσοστό ίδιος με αυτόν του 1ου μέρους. / The scope of this thesis is the study of the behavior of the human skin surface at the application of a vertical point as well as a non-point force on it. A model of a skin part is developed for the needs of this study, with use of finite elements method. The aim is to introduce the mathematical model in a robotic system which will inject causing the least possible pain or massage with maximum efficiency. In this thesis the robot is not being programmed, but it is been developed the code in matlab which by an appropriate transformation it could be used by the robot. First of all, the finite element method is been analyzed in order to select an appropriate model for our problem. Assuming axial symmetry is chosen to use a two-dimensional structure consisting of two layers of different elasticity (supposedly, skin and subcutaneous tissue). For the modeling the procedures we have to follow and the standards we should take into account are being analyzed. Furthermore we analyze the approach of the problem of deformation in two dimensions mathematically. At this stage, construction and selection of the optimal model, the simulation is entirely in matlab application. Then, after having built the model in two dimensions, the parameters for the relative thickness of both layers and their relative flexibility to the model to reflect reality must be defined. To achieve comparative evaluation, two groups of actual experimental values described in related studies available in the literature. Several pairs of parameters are being tested and compared in order to do a sensitivity analysis of the model, and found this combination of values for which there is little deviation from the experimental data. Finally it is worth mentioning that the second part of the work where the behavior of the skin in non-spot stress is being modeled there are used the same experimental data and data from other modeling efforts. The code used is a large proportion same as the one used in the 1st part.