If you're willing to make some approximations, here's a way to derive it.
Remember the relationship between the traction vector $T$ and the stress tensor $\sigma$:
$$T_{j} = n_{i}\sigma_{ij}$$
where $n$ is the unit normal vector of a differential area element being analyzed.
If we take the cross-sectional area $A_0$ of the object in question to not change under application of force (only true for very small deformations), and that the tractions are approximately uniform over the surface of application, we can multiply by this constant to get a physically relevant force on the left-hand side:
$$F_{j} = A_0n_{i}\sigma_{ij}$$
Inserting the linear stress-strain relationship involving the stiffness tensor $C$, we find:
$$F_{j} = A_0n_{i}C_{ijkl}\epsilon_{kl}$$
If we also approximate that the displacements are approximately linear, that is, that $u_{i,j} \approx \frac{x_i}{(x_0)_j}$ where $x_i$ is a displacement from rest and $(x_0)_j$ is a vector representing rest length in the $j$th coordinate, we can recall the definition of the infinitesimal strain tensor $\epsilon_{ij} = \frac{1}{2}(u_{i,j} + u_{j,i})$ to get a relationship between forces and displacements:
$$F_{j} = \frac{1}{2}A_0 n_i C_{ijkl}\left(\frac{x_k}{(x_0)_l} + \frac{x_l}{(x_0)_k}\right)$$
This looks hairy, but the linear relationship you're looking for between forces and displacements is in there! To make it look nicer, let's define a tensor of the form:
$$\boxed{k_{jk} = \frac{A_0 n_i C_{ijkl}}{\left(x_0\right)_{l}}}$$
Noting the symmetries of the stiffness matrix, we find that we can rewrite the relationship as:
$$F_i = k_{ij}x_j$$
where the previously defined second-order tensor $k_{ij}$ is precisely the tensor you're looking for!