Hello, i am new with static simulations so sorry for this basic question- i made very simple model of shaft and applied force on one area. I wanted to compare results by hand and by software. I made somewhere with my calculatons mistake and i am not able to get result like from software. Could you please correct my hand calculation? or could you send me link for some similar examples (where each step is explained)?

Thanks a lot, please follow my shares excel file with attached picture.

You have not got the restraints right in ansys. You have issues related to Poisson's ratio when you fix the end of the shaft that are not accounted for in the hand calc. If you use a qtr symmetry model and and only use a sliding restraint on the three faces you will get the average stress (the hand calc) to show up as the maximum (or pretty close) as the induced stress risers will be minimized at the BC. If you check the average stress in the Ansys model it will be very close to the right number.