As there are geometrical nonlinearity and concrete shrinkage and creep in long-span flexible structures such as concrete cable-stayed bridges, so, based on differential method, the geometrical nonlinearity balance equation for large rotation displacement - small strain analysis was deduced under co-rotational coordinate system, which considered the initial strain.Combining with finite element formula for equivalent nodal force of shrinkage and creep by using initial strain method, then, through building total and incremental relationships derived from differential equations of nodal displacements and forces, a practical algorithm for geometrical nonlinearity analysis of the plane beam considering concrete shrinkage and creep under global coordinate was obtained, and its calculation flowchart was also given.The geometrical nonlinearity analysis of the tower of long-span concrete cable-stayed bridge with hybrid girder considering concrete creep was performed.The results demonstrate that the algorithm developed can solve the above mentioned problems and has some engineering application value.