2d plane stress gear pair

Sorry, I should have told you that I previously checked this post: Zero reaction moment when reading from ROT NODE in 2D

However, I do not want to control my simulation applying a rotation instead of a torque. Indeed, I reached to make it work if I apply a rotation to the pinion… but however the results is very strange since the teeth do not deform as expected…

In my case i want to obtain load distribution information along the line of action.

Any work around or any additional suggestions?
Thank you very much!