The main objective of this thesis consists on the development and application of a computational program, ANLOG (Non Linear Analysis of Geotechnical Constructions), to analyse geotechnical problems coupling flow and deformation, and considering constitutive non-linear relations for the soil modelling.At the global level, the solution of the non linear system of equations, generated by the FEM formulation of this problem, is done using the Newton-Raphson Modified method with automatic load and time increments. On the other hand, at the local level, the constitutive equations are integrated using the explicit algorithm of stress integration with sub steps. Non-linear models, elastic (hyperbolic) and elastoplastic (modified Camclay, Lade and Lade-Kim) are utilised to represent the stress-strain-strength behaviour of the soil. The procedures of turn on gravity and the one proposed by Brown and Booker (1985) are used to simulate numerically, respectively, embankments and excavations.Initial analyses of excavations and embankments are performed in which it can be observed that the response is a function of time and depends strongly on the permeability coefficient and the loading velocity. In particular, drained and undrained situations are special cases of a coupling analysis. With respect to the use of the automatic load and time increment, it is shown that its use is fundamental to define the best increment size to adopt in the incremental iterative solution. About the stress integration algorithm utilised, it is shown that the error in the consistency condition decreases exponentially with the number sub steps. Finally, the ANLOG program is used to analyse the experimental excavation of Camboinhas (Silva, 1979). Comparisons between numerical and field displacements present a reasonable agreement. However, the pore pressure given by the numerical simulation varied between 10 percent to 30 percent in relation to its initial value, whereas in the field this variation was smaller than 3 percent.