Ce travail de recherche a pour but d'étudier l'analyse de flexion thermomécanique des plaques fonctionnellement graduées épaisses en proposant une nouvelle théorie quasi-3D de déformation de cisaillement. Le modèle mathématique utilisé propose seulement 5 variables comme dans le cas de la théorie de la déformation au cisaillement du premier ordre (FSDT). Contrairement aux théories d'ordres élevés (HSDT) cette théorie présente un nouveau champ de déplacement qui comprend des variables intégrales indéterminées. Les propriétés mécaniques de la plaque fonctionnellement graduée sont supposées changer dans le sens de l'épaisseur selon une loi de puissance (P-FGM). Les équations gouvernantes pour l'étude de flexion thermomécanique sont obtenues par le principe du travail virtuel et résolues par une méthode de type Navier. Des résultats intéressants sont déterminés et comparés aux résultats trouvés par les théories HSDT et 2D. Les influences de l'épaisseur de la plaque et de l'index de loi de puissance sur la flèche et les contraintes des plaques FGM épaisses seront représentées.