A variational asymptotic composite beam model has been developed for thermoelastic analysis. Composite beams, including sandwich structure and laminates, under different boundary conditions are examined. Previously developed beam model, which is based on variational-asymptotic method, is extended to incorporate temperature-dependent materials experiencing large temperature changes. The recovery relations have been derived so that the temperatures, heat fluxes, stresses, and strains can be recovered over the cross-section. The present theory is implemented into the computer program VABS (Variational Asymptotic Beam Sectional analysis). Numerical results are compared with the 3D analysis for the purpose of demonstrating advantages of the present theory and use of VABS.