How Force boundary condition works for tet and hex mesh?
Solution
The force on a surface is distributed to the underlying mesh in a manner that yields a uniform pressure across the mesh. The forces are distributed to the nodes using the element shape functions.
The linear shape functions distribute that in a manner that results in each node on that face getting an equal percentage of that faces load (1/3 for trias, 1/4 for quads). With parabolic shape functions, the forces end up being 0 at the corner nodes so you only see values on the midnodes.
This is the standard one used across all commercial FEA codes
Notes