<abstract> <p>In this paper, in order to provides a powerful new tool for quantitative and qualitative analysis of dynamics properties in flexible mechanical multibody systems, the symmetry theory and numerical algorithms for preserving structure in modern analytical mechanics is introduced into flexible multibody dynamics. First, taking the hub-beam systems as an example, the original nonlinear partial differential-integral equations of the system dynamics model are discretized into the finite-dimensional Lagrange equations by using the assumed modal method. Second, the group analysis theory is introduced and the criterion equations and the corresponding conserved quantities of Noether symmetries are given according to the invariance principle, which provide an effective way for analytic integral theory of dynamic equations. Finally, a conserved quantity-preserving numerical algorithm is constructed by coordinates incremental discrete gradient, which makes full use of the invariance of conserved quantity to eliminate the error consumption for a long time. The simulation results show that the deeper mechanical laws and motion characteristics of flexible mechanical multibody systems dynamics can be obtained with the help of symmetries and conserved quantities, which can provide reference for more precise dynamic optimization design and advanced control of systems.</p> </abstract>