Based on the CR full nonlinear calculation model and the Timoshenko beam theory，and considering the beam element shear effect，the beam-column theory（stability function format） is adopted to establish a fine and efficient plane beam element analysis theory，which considers the axial force，bending moment，shear and their coupling effect of the element. The beam-column theory and the CR coordinate method are used to improve the fast solving method for the geometric nonlinear problem of plane beam element under the high precision target. On the basis of the two order beam-column theory，the concrete form of the angular displacement equation for the shear effect is derived，and the exact stability function expression of the stiffness coefficient is obtained. The power series of the above stability function is transformed and the corresponding truncation criterion is put forward on the premise of ensuring the numerical accuracy and the stability of the calculation process. Subsequently，the formula for calculating the node force of the element considering the effect of shear deformation on the bowstring effect is derived，and the program is compiled and compared with classic examples. Moreover，it is analyzed with the actual engineering. The results show that，for the general multi-span continuous beam structure，the calculation results of the classical calculation can be obtained by the calculation theory proposed in this paper. For the class of large displacement，such as the elbow rigid frame structure，the full process deformation curve of the load-displacement considering the shear effect can be obtained by this theory，and the first and two critical points of leap buckling are also consistent with the classical literature. For the long-span cable-stayed bridge with different main girder sections，reasonable overall deformation characteristics can be obtained by using this theory，and the efficiency of analysis can also be significantly improved.