The rigid-flexible coupling dynamic characteristic of a rotating hub-beam system in a temperature field was investigated. Considering geometric nonlinearity and thermal effect, based on the exact nonlinear strain-displacement relationship, the rigid-flexible coupling dynamic equations of the hub-beam system were established using virtual work principle and finite element method. Since the nonlinear stiffness matrices are related to the high order terms of deformation, each element of the nonlinear stiffness matrices is expressed as the product of generalized coordinate vectors and constant matrices. Simulation results show that, by using such formulation, duplicated integration can be avoided and computational efficiency can be significantly improved. Furthermore, the effect of geometric nonlinearity on rigid-flexible coupling dynamics in case of temperature increase was investigated. Spectrum analysis method was used to investigate the change of natural frequencies with the moment of inertia of the hub and the temperature.