Why do we get the same differential equations from both principles? Surely there is a fundamental connection between them? When written out, the two seem to have nothing in common.
$$\sum _i ( \mathbf F _i - \dot{\mathbf p}_i) \cdot \delta \mathbf r _i = 0$$
$$S[q(t)] = \int ^{t_2} _{t_1} \mathcal L (q,\dot{q},t)dt$$
After playing with d'Alembert's principle we find that we can rewrite the whole thing as $$\sum _i \left[ \frac{d}{dt} \left( \frac{\partial T}{\partial \dot{q}_i} \right) - \frac{\partial T}{\partial q_i}-Q_i\right] \delta q_i$$
This can further be rewritten under certain conditions so we get the exact form of the E-L equation.
It seems to me that both ways of arriving at the result are fundamentally different.
A function must obey the E-L equations in order to minimize the action over a path, but when we look at the virtual work, it appears that they come from the fact that (quoting Goldstein) "particles in the system will be in equilibrium under a force equal to the actual force plus a 'reversed effective force'."
I think I understand the principle of stationary action, I can see how it leads to the E-L equations, but d'Alembert's Principle seems so arbitrary, I can't see any motivation for it.