Abstract: Considering nonlinear aerodynamic, the aero-elastic equation in supersonic flow is constructed by applying the third piston theory based on Von Karman's nonlinear panel theory. Then the equation is scattered in modal space by Galerkin method. The influence of thermal buckling amplitude and flutter amplitude by nonlinear aerodynamic is analyzed. The stability of heated panels and the boundary of flutter is determined using Lyapunov indirect method. The nonlinear flutter responses of the heated panel are numerically simulated in time domain. The results show that there are complex dynamic phenomenons such as convergent motion, limit cycle oscillation and chaotic motion.