Dear Professor Taylor,
Thank you very much for your help.
The problem has been solved by using BLENd command for both of mesh generation and contact surfaces. Of course, SNODes are defined including the center of the disk. I have tried to fix the center of disk and apply a tangential displacement through contact friction between the lower disk and a upper brick block. The disk rotate well. But when I change the disk with a ring (disk with a hole in the middle), FEAP seem not to set up the SNOD center nodes as a rotating axis.
Thank you for your time again.
Best regards,
Vinh Hiep Nguyen