Reference is made to OS-V: 0085 which uses nodal pressure to compute a pressure vessel (modelled as quarter segment).
I'm using different parameters to study a thick-walled cylinder. However, results show a kind of "hourglassing" effect with nodal pressure. It is not real hourglassing because it is only localised to the inner radius, elements are already fully integrated (second order) and the effect does not disappear with mesh refinement. See the figure below:

After checking the reference manual, it seems to be possible to try "PLOAD2" as alternative to radial nodal pressure. The figure below shows that PLOAD2 pressure applies to the element edges instead of the nodes:

But now I'm stuck, because these loads are not included in the analysis. Did I miss something? Could CQPSTN elements be incompatible? Or is there a better way to define edge pressure on a 2D model? See the model in attachment.