Αντικείμενο του παρόντος άρθρου αποτελεί η ανάλυση της συμπεριφοράς πασσάλου, υπό αξονική φόρτιση μεγάλου εύρους, η οποία μπορεί να οδηγήσει σε απώλεια φέρουσας ικανότητας. Συγκεκριμένα, εξετάζεται ο συντελεστής στατικής δυσκαμψίας μεμονωμένου κατακόρυφου πασσάλου εμπεδωμένου σε ομοιογενές ή πολυστρωματικό έδαφος τυχαίας γεωμετρίας και μηχανικών ιδιοτήτων. Για την επίλυση του προβλήματος αναπτύσσεται αναλυτική λύση κλειστού τύπου βασισμένη στη θεωρία Winkler. Στο εν λόγω μοντέλο η προσομοίωση της μηχανικής συμπεριφοράς του εδάφους γίνεται μέσω μή-γραμμικών ελατηρίων “t-z” τοποθετημένων κατά μήκος του άξονα και στη βάση του πασσάλου, σε συνδυασμό με συναρτήσεις σχήματος, οι οποίες περιγράφουν αξιόπιστα την μεταβολή της κατακόρυφης μετακίνηση του πασσάλου με το βάθος. Με επιλογή κατάλληλης συνάρτησης σχήματος και καμπυλών “t-z”, και μετά από επαναληπτική διαδικασία εφαρμογής, επιτυγχάνεται ικανοποιητική ακρίβεια στην τιμή της δυσκαμψίας για κατακόρυφη μετακίνηση στην κεφαλή του πασσάλου. Σε αντίθεση με τις κλασικές αριθμητικές λύσεις, η προτεινόμενη μέθοδος δεν απαιτεί διακριτοποίηση του πασσάλου σε πεπερασμένα στοιχεία (και στη συνέχεια επίλυση ενός συστήματος γραμμικών εξισώσεων μεγάλης τάξης), παρά μόνο σε "κελιά" με στόχο την ολοκλήρωση με το βάθος. Έτσι, τα παραγόμενα αποτελέσματα είναι διαχείρισιμα ακόμα και μέσω απλού φύλλου εργασίας σε Excel ή και υπολογιστή τσέπης. Η μέθοδος προγραμματίστηκε σε περιβάλλον Visual Basic 2010, κυρίως λόγω της δυνατότητας γραφικής παρουσίασης των αποτελεσμάτων και τη σύγκρισή τους με αντίστοιχα αποτελέσματα από άλλες μεθόδους. Τα αποτελέσματα κρίνονται ως ιδιαίτερα ενθαρρυντικά, καθώς συγκλίνουν ικανοποιητικά σε αυτά αυστηρότερων μεθόδων, χωρίς ανάγκη περίπλοκης αριθμητικής ανάλυσης η οποία να ξεφεύγει από τις γνώσεις και δυνατότητες του Γεωτεχνικού Μηχανικού. / In the present work the behavior of a pile submitted to large range lateral loading is analyzed, which may lead to failure of both the surrounding soil and the pile itself either at the head or in depth. Namely, we examine the static stiffness coefficients for displacement and rotation of a flexible pile, vertically embedded in a homogeneous or multilayer soil of random geometry and mechanical properties. To solve the problem, a simple analytical method is developed, based on Euler–Bernoulli classic beam model, incremented with non linear Winkler Springs. The non-linear behaviour of the pile is described in a cross-sectional plane through moment-curvature diagram. The model is used in combination with the principle of work and suitable shape functions, which describe reliably the elastic line of the pile when the lateral load is gradually increasing. By iterative implementation of the method, realistic predictions are achieved in the stiffness coefficients in swaying, rocking and cross-swaying-rocking. The number of iterations is relatively small if the stress level of the system is not significantly increased compared with the previous load step. Unlike classic numerical solutions, the proposed method does not require discretization of the pile into finite elements (resulting to solve a system of linear equations), but only in "cells", to integrate with depth. In this way, results can be generated throughout a simple worksheet or even a calculator. The method was implemented in a Visual Basic 2010 environment, mainly for reasons of graphical presentation and comparison of the results to other coming from relevant methods. The results of the aforementioned method are considered satisfactory, as they converge fairly well with those coming from more rigorous methods based on complicated numerical analyses. The results of the herein proposed method are also compared to experimental in situ results relatively successfully.