A function $f(x)$ is defined on the set of all real numbers. It satisfies the relation $f(x + y) = f(x) + f(y) + 2x^2y + 2xy^2$ for all $x, y \in \mathbb{R}$. If $f(1) = 1$, what is the value of $f(3)$?