Nodes of solid elements don’t have rotational degrees of freedom and some workarounds are required to apply boundary conditions or loads to rotational DOFs. In PrePoMax you can use rigid body constraint and apply BCs/loads to its reference point. Like here: PrePoMax (CalculiX FEA) - Tutorial 2 - Torsion of an elliptical bar - YouTube
Local cylindrical coordinate system could also help but it’s not yet supported in PrePoMax and would require keyword edits.
Thank you very much for your help.
So, if my understanding is correct, your recommendation is to make a kinematic coupling of journal surface to an origin (bearing center) through Rigid Body & Reference point. And then, assuming the shaft direction is Z, I guess you recommend to set the reference point U1, U2, U4, U5 fixed, and U3 and U6 uncostrained. Is my understanding correct, or?