The existing finite element software mainly uses solid and shell elements of the preset crack to simulate the bending effect of the crack structure but not the beam element. Therefore, in this paper, in the simplified computational model based on the Euler-Bernoulli theory, the Dirac-delta function is introduced into the beam element"s curvature function to represent the crack, the displacement distribution function are obtained. Then, the cracked beam element is proposed by deriving the strain matrix and the element stiffness matrix. Based on the theory of co-rotational coordinate method (C.R formulation), the geometric nonlinear crack beam element considering crack structure is further developed, and the specific algorithm and process of the element are given. Using the finite element solution framework developed by our research group, the crack beam element and the co-rotational crack beam element are realized in the program. Through typical examples, the effects of crack location, crack depth and external load on the bending performance of cracked beams are compared and analyzed. The results show that the element can accurately describe the deflection change of the crack structure, and has reliable accuracy and computational efficiency in linear and geometric nonlinear analysis compared with the existing elements for simulating the deformation of the crack structure.