A new analytical model for thermoelastic responses of a multi-layered composite plate with imperfect interfaces is developed. The composite plate contains an arbitrary number of layers of dissimilar materials and is subjected to general mechanical loads (both distributed internally and applied on edges for each layer) and temperature changes, which can vary from layer to layer and along two in-plane directions. Each layer is regarded as a Kirchhoff plate, and each imperfect interface is described using a spring-layer interface model, which can capture discontinuities in the displacement and stress fields across the interface. Unlike existing models, the governing equations and boundary conditions are simultaneously derived for each layer by using a variational procedure based on the first and second laws of thermodynamics, which are then combined to obtain the global equilibrium equations and boundary conditions for the multi-layered composite plate. A general analytical solution is developed for a symmetrically loaded composite square plate with an arbitrary number of layers and imperfect interfaces by using a new approach that first determines the interfacial normal and shear stress components on one interface. Closed-form solutions for two- and three-layer composite square plates are obtained as examples by directly applying the general analytical solution. Numerical results for two-, three- and five-layer composite plates under different loading and boundary conditions predicted by the current model are provided, which compare well with those obtained from finite element simulations using COMSOL, thereby validating the newly developed analytical model.