[...] that the angular frequency should match with the damped frequency so that the force is in phase with the motion.
That's not really true.
You're also conflating frequency and phase: two sines can have the same frequency yet different phase. I think you're a little puzzled as to why the natural oscillator frequency isn't the same as the driving force frequency. That I will address here.
The equation of motion of the driven, damped oscillator is:
$$m\ddot{x}+c\dot{x}+kx=F_0\cos(\omega t+\varphi_d)\tag{1}$$
In the underdamped case the solution is:
$$x(t)=A_he^{-\gamma t}\sin(\omega't+\varphi_h)+A\cos(\omega t-\varphi)$$
The first part of the RHS is the transient solution, the second part of the RHS is the steady state solution. $\omega'$ is the damped angular velocity (see derivation here).
It's clear that for $t\to +\infty$, the transient solution decays to zero, i.e. $A_he^{-\gamma t}\sin(\omega't+\varphi_h)=0$.
So in steady state the driving force and the motion have the same frequency: that is the frequency of the driving force ($\omega=2\pi f$)
This is the consequence of how differential equations (DE) like $(1)$ are solved. First, a particular solution is found for the homogeneous DE:
$$m\ddot{x}+c\dot{x}+kx=0$$
That particular solution is the transient solution.
Then a second solution is 'guessed', that is the steady state solution.
Superposition then demands the actual solution of $(1)$ is the sum (linear combination) of both solutions.