I think there is a problem with beam line pressure load.
In the attached model the solution is all zeros, but nodes rotations should not be zero.
stiffness matrix is ok but load vector is all zeros but it shold be
0.
0.
-16666.666666667
0.
0.
16666.666666667
so solution shoud be
0.
0.
-0.0396825396825
0.
0.
0.0396825396825
Reaction forces are OK, so I think that there is a problem with nodal equivalent forces calculation, equivalent moments are missing.
Comments
The same is true for pressure on shells.
Would you say this is an important omission?
If there are only nodal loads it is not necesary to refine, and it would be fine if this can work with line pressure loads too.