First, since $t$ does not enter into the given relation we suppress the time-dependence of $f$. Additionally, while the problem presumably assumes ordinary 3D space I'll consider $f$ as a function $\mathbb{R}^n\to \mathbb{R}$ (i.e., $n$-dimensional space). Now let $\mathbf{x}=\mathbf{r}_{CA}-\mathbf{r}_{BA}=\mathbf{r}_{CB}$, $\mathbf{y}=\mathbf{r}_{BA}$ and rearrange the relation as $$f(\mathbf{x})+f(\mathbf{y})=f(\mathbf{x}+\mathbf{y})$$
Since the points $A,B,C$ are arbitrary, the above relation should likewise be valid for arbitrary real $\mathbf{x},\mathbf{y}$. It is easy to see that the $f(\mathbf{x}):=\mathbf{c}\cdot \mathbf{x}=\sum_{i=1}^n c_i x_i$ is a solution for any $\mathbf{c}\in\mathbb{R}^n$ which I take to be the desired "Hubble law". Is this the only solution?
This is a multidimensional version of Cauchy's (additive) functional equation. Via the MathOverflow question here, the multidimensional version (assuming continuous $f$) appears in Kuczma's An Introduction to the Theory of Functional Equations and Inequalities. Quoting from the answer to the aforementioned Q&A: "[Theorem 5.52] shows that the only continuous solutions are precisely the functions of the form $\mathbb{R}^n\to \mathbb{R}:(x_1,\ldots,x_n)\mapsto \sum_{i=1}^n c_i x_i$ with $c_1,\ldots,c_n\in\mathbb{R}$." The proof is a straightforward extension of the one-variable version; this has appeared on MSE many times (see for instance Martin Sleziak's summary here) so I won't consider this further.
What if we drop the condition of continuity? In that case there do exist nontrivial solutions to Cauchy's additive functional equation. However, these solutions are well-known to be quite pathological; in particular, such functions are known to not be measurable. For a summary of basic results on Cauchy's functional equation in the 1D case, see Martin Sleziak's summary as referenced above; for more definitive results in the multidimensional case, Kuczma seems the appropriate reference.