Thermal-electrical-chemical-mechanical coupling controls the behavior of many trans port and electrochemical reactions processes in physical, chemical and biological sys tems. Hence, advanced understanding of the coupled behavior is crucial and attracting a large research interest. However, most o f the existing coupling theories are limited to the partial coupling or particular process. In this paper, on the basis of irreversible thermo dynamics, a variational principle for the thermal electrical chemical mechanical fully coupling problems is proposed. The complete fully coupling governing equations, includ ing the heat conduction, mass diffusion, electrochemical reactions and electrostatic potential, are derived from the variational principle. Here, the piezoelectricity, conduc tivity, and electrochemical reactions are taken into account. Both the constitutive rela tions and evolving equations are fully coupled. This theory can be used to deal with coupling problems in solids, including conductors, semiconductors, piezoelectric and nonpiezoelectric dielectrics. As an application of this work, a developed boundaiy value problem is solved numerically in a mixed ion-electronic conductor (M1EC). Numerical results show that the coupling between electric field, diffusion, and chemical reactions influence the defect distribution, electrostatic potential and mechanical stress.