Abstract—Recently, there are many studies of the nonlinear structural analysis. This paper presents the method to predict nonlinear response structures. A numerical approach is proposed for the nonlinear analysis structures problems using displacement method. The formulation is based on the rope stiffness matrix of the beam element, which is presented taking into account the nonlinearity of the materials and second-order effects due to node displacements. The tangent stiffness matrix of the element is obtained by the virtual work principal. In this numerical nonlinear study, the equilibrium of the structure is solved using iterative methods. The Simpson integration scheme is used to integrate the surface of the section. The system of equations is solved in the absolute coordinate. The method is implemented in Fortan 90. The numerical examples with bending fracture modes are presented. Considering the load of the structure, verification of the reliability of the developed software is done by comparing numerical result with experimental and analytical results.