The thermodynamic bending behavior of multilayered functionally graded carbon nanotube reinforced composite (FG-CNTRC) plates has received limited attention in the existing literature, despite their promising applications in advanced aerospace, mechanical, and nano-engineering structures. To address this gap, the present study develops, for the first time, an analytical model to investigate the thermodynamic bending response of multilayered FG-CNTRC plates subjected to thermal loading. The formulation is based on the higher-order shear deformation theory (HSDT), while the governing equations and boundary conditions are derived using Hamilton's principle. An analytical solution for simply supported plates is obtained by employing Navier's method. The material properties are assumed to be temperature-dependent and are evaluated according to the Touloukian principle. In addition, the volume fraction of carbon nanotubes (CNTs) is considered to vary through the thickness according to several distribution patterns, namely uniform distribution (UD) and four functionally graded distributions (FG-X, FG-O, FG-V, and FG-A). A comprehensive parametric investigation is performed to examine the effects of CNT distribution patterns, CNT volume fraction, temperature difference, and number of layers on the bending response of the plates. The numerical results reveal that the FG-X distribution exhibits the highest stiffness and consequently the lowest deflection among all the considered distributions. Furthermore, increasing the CNT volume fraction and the number of layers enhances the structural stiffness, whereas increasing the temperature difference leads to a significant reduction in stiffness and an increase in deflection. The proposed model provides new insights into the thermodynamic bending behavior of multilayered FG-CNTRC plates and may serve as an efficient analytical tool for the design and optimization of advanced nanocomposite structures.