sympy physics mechanics

[Peter Stahlecker]

I calculated a very simple pendulum system, using one non - holonomic constraint. All worked as per expectation. Then I tried to calculate the reactive forces at the point of attachment, using KM.auxiliary_eqs. (KM being the Kanes object).

I replaced the accelerations of the generalised coordinates by the respective entries of rhs.

(rhs = mass_matrix)**(-1) * forcing )

However, now the results did not meet expectations at all.

Did I do something wrong, or is it just not possible to do what I tried?

Thanks a lot for any comments / suggestions!

[Peter Stahlecker]

Forget my question above: I DID make a stupid mistake!

All works fine!