Large space structures would exhibit orbit-attitude-structure-thermal coupled effects in complicated space environment. Thus, a new thermal-structure coupled modelling method is proposed using gradient-deficient absolute nodal coordinate formulation beam elements. Compared to the previous methods, both axial and circumferential heat conduction are considered. The cubic interpolation in temperature field instead of the linear interpolation is adopted to obtain a unified-mesh temperature-displacement description. The gravitational force and gravity gradient are modelled, and the effects of attitude motion and structural deformation on solar radiation intensity are considered. In addition, three kinds of quasi-static thermal-structure coupled models are proposed based on Euler-Bernoulli beam theory. Based on the quasi-static models, the theoretical formulas are derived to effectively predict the thermally induced vibrations of the beam in space. Four numerical examples are studied to validate the proposed models, including the heat conduction examples with three boundary conditions and a thermal-structure coupled cantilever beam example. Based on the proposed formulation, the orbit-attitude-structure-thermal coupled dynamics of a multibody system consisted of a rigid body and a large flexible appendage is investigated considering the gravity gradient, solar radiation, Earth's shadow, and self-shadow.