hertzian contact pressure

Hello all,

I have applied a moving ellyptical hertzian pressure load on a surface using a DLOAD which works, but I also need the traction load which is the multiplication of this pressure with a friction coeffcient.

Does anybody have any clue how to do it? Directly through abaqus or inside the same DLOAD?

Kind regards,

Omid