We apply the Wigner function formalism from quantum optics via two approaches, Wootters' discrete Wigner function and the generalized Wigner function, to detect quantum phase transitions in critical spin-12 systems. We develop a general formula relating the phase space techniques and the thermodynamical quantities of spin models, which we apply to single, bipartite and multi-partite systems governed by the XY and the XXZ models. Our approach allows us to introduce a novel way to represent, detect, and distinguish first-, second- and infinite-order quantum phase transitions. Furthermore, we show that the factorization phenomenon of the XY model is only directly detectable by quantities based on the square root of the bipartite reduced density matrix. We establish that phase space techniques provide a simple, experimentally promising tool in the study of many-body systems and we discuss their relation with measures of quantum correlations and quantum coherence.