For monogenic and a special case of non-holonomic constraints where we have$$ \sum_{k} a_{l k} d q_{k}+a_{t t} d t=0 \tag{2-20} $$ we use lagrange multipliers and hamiltons principle to reach the below equation:
Goldstein, 2 edition
$$\frac{d}{d t} \frac{\partial L}{\partial \dot{q}_{k}}-\frac{\partial L}{\partial q_{k}}=\sum_{l} \lambda_{l} a_{l k}, \quad k=1,2, \ldots, n.\tag{2-29}$$
And he says further
What is the physical significance of the $\lambda_{i}$ 's? Suppose one removed the constraints on the system, but instead applied external forces $Q_{k}^{\prime}$ in such a manner as to keep the motion of the system unchanged. The equations of motion would likewise remain the same. Clearly these extra applied forces must be equal to the forces of constraint, for they are the forces applied to the system so as to satisfy the condition of constraint. Under the influence of these forces $Q_{k}^{\prime}$, the equations of motion are $$ \frac{d}{d t} \frac{\partial L}{\partial \dot{q}_{k}}-\frac{\partial L}{\partial q_{k}}=Q_{k}^{\prime}.\tag{2-31} $$ But these must be identical with Eqs. (2-29). Hence we can identify $\sum \lambda_{1} a_{l k}$ with $Q_{k}^{\prime}$, the generalized forces of constraint.
Question: How does it follow that "Under the influence of these forces $Q_{k}^{\prime}$, the equations of motion are eq. (2-31)?"
Since we've a monogenic system, potential maybe a function of velocities $V=V(r,\dot r,t)$ therefore replacing constraints by equivalent forces and using d'Alemberts principle we might get $$\frac{d}{d t}\left(\frac{\partial T}{\partial \dot{q}_{j}}\right)-\frac{\partial T}{\partial q_{j}}=Q_{j}$$ and not eq. (2-31).
Can anyone shed some light please.