In my opinion the method I outlined is always conservative: as I pointed out it is normally used for flanged joints.
Of course you are right when you point the finger on a possible effect of non uniformity around the circumference: after all on large girth flanges there are dozens of bolts, whilst for a support ring to have 4 support points only is a common situation.
However I think that this effect may be accounted for quite easily.
First one has to note that the increased membrane stress near the support would be classified as a local stress under Div.2 coverage, and thus limited to 1.5
S. Also this is a longitudinal stress, so its value is unlikely to impact on shell thickness, that, I assume, is normally determined by circumferential stress under pressure.
As a consequence you might conservatively estimate the width of shell subject to the longitudinal support stress as follows: take the effective supported width of ring (e.g. width of the supporting beam), add the height of the ring, this will give a conservative width of shell onto which to distribute the support load. If you can pass this check with 1.5
S as the allowable (also considering buckling effects under compression, if you have this), then you should be on the conservative side.
prex
Online tools for structural design