Combining the symplectic variations theory, the homogeneous control equation and isoparametric element homogeneous formulations for piezothermoelastic hybrid laminates problems were deduced. Firstly, based on the generalized Hamilton variation principle, the non-homogeneous Hamilton canonical equation for piezothermoelastic bodies was derived. Then the symplectic relationship of variations in the thermal equilibrium formulations and gradient equations was considered, and the non-homogeneous canonical equation was transformed to homogeneous control equation for solving independently the coupling problem of piezothermoelastic bodies by the incensement of dimensions of the canonical equation. For the convenience of deriving Hamilton isoparametric element formulations with four nodes, one can consider the temperature gradient equation as constitutive relation and reconstruct new variation principle. The homogeneous equation simplifies greatly the solution programs which are often performed to solve non-homogeneous equation and second order differential equation on the thermal equilibrium and gradient relationship.