We construct finite element spaces of symmetric stress tensors that are exactly divergence-free. Moreover, their basis functions can be chosen so that they have small supports. These properties are highly desired in a number of important applications. Approximation properties of finite element spaces of divergence-free tensor functions are derived from properties of C1 finite elements.