In this article we investigate transcendental entire solutions of any order of the following partial differential difference equations: \(\begin{aligned} f(z_1+c_1,z_2+c_2)^2+\left( \frac{\partial f(z_1,z_2)}{\partial z_1}\right) ^2=e^{2h(z_1,z_2)}\end{aligned}\) and \(\begin{aligned}{\left\{ \begin{array}{ll} f(z_1+c_1,z_2+c_2)^2+g(z_1+c_1,z_2+c_2)^2=e^{2h(z_1,z_2)}, \\ \frac{\partial g(z_1,z_2)}{\partial z_1}^2+\frac{\partial f(z_1,z_2)}{\partial z_1}^2=1-e^{2h(z_1,z_2)}, \end{array}\right. }\end{aligned}\) where \(h(z_1,z_2)\) is entire in \(\mathbb {C}^2\) and obtain two new important and interesting results, which are the generalization of some of the previous results. Some examples have been exhibited which show that our results are precise.