The quantum phase transition (QPT) describes a sudden qualitative change of the macroscopic properties mapped from the eigenspectrum of a quantum many-body system. It has been studied intensively in quantum systems with the spin-boson model, but it has barely been explored for systems in coupled spin-boson models. In this paper, we study the QPT with coupled spin-boson models consisting of coupled two-level atoms embedded in three-dimensional anisotropic photonic crystals. The dynamics of the system is derived exactly by means of the Laplace transform method, which has been proven to be equivalent to the dissipationless non-Markovian dynamics. Drawing on methods for analyzing the ground state, we obtain the phase diagrams through two exact critical equations and two QPTs are found: one QPT is that from the phase without one bound state to the phase with one bound state and another is that from one phase with the bound state having one eigenvalue to another phase where the bound state has two eigenvalues. Our analytical results also suggest a way of control to overcome the effect of decoherence by engineering the spectrum of the reservoirs to approach the non-Markovian regime and to form the bound state of the whole system for quantum devices and quantum statistics.