We consider the possibility that the Universe is made of a dark fluid described by a quadratic equation of state P = Kρ2, where ρ is the rest-mass density and AT is a constant. The energy density ϵ = pc2 + Kρ2 is the sum of two terms, a rest-mass term ρc2 that mimics "dark matter" (P = 0) and an internal energy term 11 = Kρ2 = P that mimics a "stifffluid" (P = ϵ) in which the speed of sound is equal to the speed of light. In the early universe, the internal energy dominates and the dark fluid behaves as a stifffluid (P ∼ ϵ, ϵ α(a-6). In the late universe, the rest-mass energy dominates and the dark fluid behaves as pressureless dark matter (P ≃0, ϵ α a-3). We provide a simple analytical solution of the Friedmann equations for a universe undergoing a stiffmatter era, a dark matter era, and a dark energy era due to the cosmological constant. This analytical solution generalizes the Einstein-de Sitter solution describing the dark matter era, and the ACDM model describing the dark matter era and the dark energy era. Historically, the possibility of a primordial stiffmatter era first appeared in the cosmological model of Zel'dovich where the primordial universe is assumed to be made of a cold gas of baryons. A primordial stiffmatter era also occurs in recent cosmological models where dark matter is made of relativistic self-gravitating Bose-Einstein condensates (BECs). When the internal energy of the dark fluid mimicking stiffmatter is positive, the primordial universe is singular like in the standard big bang theory. It expands from an initial state with a vanishing scale factor and an infinite density. We consider the possibility that the internal energy of the dark fluid is negative (while, of course, its total energy density is positive), so that it mimics anti-stiffmatter. This happens, for example, when the BECs have an attractive self-interaction with a negative scattering length. In that case, the primordial universe is nonsingular and bouncing like in loop quantum cosmology. At t = 0, the scale factor is finite and the energy density is equal to zero. The universe first has a phantom behavior where the energy density increases with the scale factor, then a normal behavior where the energy density decreases with the scale factor. For the sake of generality, we consider a cosmological constant of arbitrary sign. When the cosmological constant is positive, the Universe asymptotically reaches a de Sitter regime where the scale factor increases exponentially rapidly with time. This can account for the accelerating expansion of the Universe that we observe at present. When the cosmological constant is negative (anti-de Sitter), the evolution of the Universe is cyclic. Therefore, depending on the sign of the internal energy of the dark fluid and on the sign of the cosmological constant, we obtain analytical solutions of the Friedmann equations describing singular and nonsingular expanding, bouncing, or cyclic universes. [ABSTRACT FROM AUTHOR]