nI ran the same model and got the result I expected, 990N, compared to your 970N, but the actual reaction force is 1000N. (I am not sure why we get such differences in rounding)nMy guess is that other software you used are either beam analysis packages, and not full FE packages.nTo understand this result, you need to understand that the FE method does not solve for force, and force is preserved only on average, across the domain. To check this out for yourself, under output controls, turn nodal forces on, then resolve. Plot the user defined result, enfovectors. You will see the 1000N reaction force, but also that each element has it's own reaction force, and at nodes shared by two elements, there are two results, both different.n