The free Student license has a node or element limit of 32,000. This model easily fits within that limit, it has 8,775 nodes.

You can use a minimum of 1 N-R Residual Force plot. It doesn't cost much to have more than one. On some problems, the location of the maximum NR residual bounces around the model so 3 to 6 is a good number to see if the location is consistent so you can be confident that you are refining the mesh at the right location.

You must have a solid body with a cylindrical face on it. If you need more help, become a member on studentcommunity.ansys.com and create a New Discussion.

Ok...Two things...1 - My hug stress level is a singularity due to the way I need to apply preload in Nastran. 2 - The author has shown X-axis directional displacement instead of Total displacement. My Nastran outcomes align nicely now. Fields looks similar, however displacement calues are a bit different in X direction, range = -0.16 to +0.0187 compared to author's -0.1 to 0.023

When I model a bolt preload for a bolted joint, I generally don't look at the stress in the bolt because I idealize the bolt leaving a sharp interior corner that creates a singularity. I use a separate Excel spreadsheet for bolt stress calculations. The Ansys model is used to evaluate the stress in the two parts that are being loaded after being bolted together.

One substep might be too few and cause the solver to do a bisection on the stepsize, wasting time. It might solve with 3 substeps, but each one might need more iterations.