One open question remains regarding the theory of the generalized variational principle, that is, why the stress-strain relation still be derived from the generalized variational principle while the Lagrangian multiplier method is applied in vain? This study shows that the generalized variational principle can only be understood and implemented correctly within the framework of thermodynamics. As long as the functional has one of the combination $A(\epsilon_{ij})-\sigma_{ij}\epsilon_{ij}$ or $B(\sigma_{ij})-\sigma_{ij}\epsilon_{ij}$, its corresponding variational principle will produce the stress-strain relation without the need to introduce extra constraints by the Lagrangian multiplier method. It is proved herein that the Hu-Washizu functional $\Pi_{HW}[u_i,\epsilon_{ij},\sigma_{ij}]$ and Hu-Washizu variational principle comprise a real three-field functional.