An analytical method based on the nonlinear layered finite element method is developed to simulate an overall load-deflection behavior of strengthened beams. The developed model distinguishes itself by its capability to trace residual flexural behavior of a beam after the fracture of brittle strengthening materials at peak load. The model. which uses a rather advanced numerical technique for iterative convergence to equilibrium, can be regarded as superior to the two models based on load control and displacement control The model predictions were compared with the experimental results and it was observed that there was good agreement between them.