This article introduces a new semi-analytical nonlinear finite element formulation for thin cylinders according to a continuum-based approach. A comparison between the continuum-based approach and the classical approach for the buckling behavior of isotropic and orthotropic perfect cylinders validates the results. The classical approach is defined according to thin shell theories based on the von Karman approximation. A mathematical modeling for geometry imperfection of cylinders is derived according to a continuum-based approach whose results are compared with the results of the classical approach for imperfect cylinders. The influence of neglecting some nonlinear terms in the classical approach for perfect and imperfect cylinders on the buckling path is investigated. In the buckling analysis, two methods, i.e. the perturbation and load disturbance methods, which undertake to switch to the post-buckling path, are compared to each other.