This paper describes a higher-order global–local theory for thermal/mechanical response of moderately thick laminated composites with general lamination configurations. In-plane displacement fields are constructed by superimposing the third-order local displacement field to the global cubic displacement field. To eliminate layer-dependent variables, interlaminar shear stress compatibility conditions have been employed, so that the number of variables involved in the proposed model is independent of the number of layers of laminates. Imposing shear stress free condition at the top and the bottom surfaces, derivatives of transverse displacement are eliminated from the displacement field, so that C 0 interpolation functions are only required for the finite element implementation. To assess the proposed model, the quadratic six-node C 0 triangular element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate. Comparing to various existing laminated plate models, it is found that simple C 0 finite elements with non-zero normal strain could produce more accurate displacement and stresses for thick multilayer composite plates subjected to thermal and mechanical loads. Finally, it is remarked that the proposed model is quite robust, such that the finite element results are not sensitive to the mesh configuration and can rapidly converge to 3-D elasticity solutions using regular or irregular meshes.