Structures suffer from damages in their lifetime due to time-dependant effects, such as fatigue, creep and shrinkage, which can be expressed by concrete strains. These processes could be seen in the context of strain estimation of pre-stressed structures in two phases by using numerical methods. Their aim is checking and validating existing code procedures in determination of deformations of pre-tensioned girders by solving a system of nonlinear equations with strains as unknown parameters. This paper presents an approach based on the Newton-Raphson method for obtaining the stresses and strains in middle span section of pre-stressed girders according the equilibrium state.