It seems to me that many different momenta $\dot{\bf p}_j $ can satisfy d'Alembert's principle:
$$\tag{1} \sum_{j=1}^N ( {\bf F}_j^{(a)} - \dot{\bf p}_j ) \cdot \delta {\bf r}_j~=~0 $$
in a constrained system.
For example, take two particles connected by a rod, with an applied non-zero force of the same size on each particle, along the line of the rod and in opposite directions (i.e. forces trying to pull the particles apart). In this case, any pair of $\dot{\bf p}_1, \dot{\bf p}_2$ of the same size and in opposite directions along the line of the rod seems to satisfy the d'Alembert's principle (taking into account the dependencies between $\delta {\bf r}_1$ and $\delta {\bf r}_2$ imposed by the rod), even though some of these do not satisfy the constraint.
Am I misunderstanding something, or does d'Alembert's principle not provide any meaningful information in this example?