The nonlinear evolution of the thermal Weibel instability is studied by using three-dimensional particle-in-cell simulations. After a fast saturation due to a reduction in the temperature anisotropy, the instability evolves to a quasistationary state which includes a single mode long wavelength helical magnetic field and a finite degree of temperature anisotropy. The nonlinear stability of this state is explained by periodic variations of the temperature anisotropy axis. At long time scales the magnetic field, wave number, and temperature anisotropy slowly evolve to the decreasing magnitudes.