Problem with cylindrical roller bearings

Hello everyone,

I am creating a numerical model of a cylindrical roller bearings using Hypermesh.
I started with one cylindrical roller with diameter 5mm. The contact area I mesh with 2D-automesh, to get the mesh finer at the contact.
Afterwards I use 'tetra mesh'-panel to create a mesh to all solids.
I have refine the contact to 0.025mm, but I do not get the analytically results of Hertzian pressure (more than 30% deviation).
The analysis is non-linear-quasi-static, the contact SLIDE, NLPARAM NINC=5 and the loadstep is set to NLSTAT.

If someone could give me tips for mesh of cylindrical rollers, I would be infinitely grateful.

And sometimes comes the ERROR # 941 NON-CONVERGENT NONLINEAR ITERATIONS FOR SUBCASE 1.

What can I do to check what I need to change in my model?

Regards,
Martin