A very conservative method can be used when e is very small:
e = horizontal distance between cylinder shell center line and skirt cylinder center line (in.)
v = skirt vertival load #/in
t= head thickness in.
m= e*v ,in-#/in unit moment
bending stress in head fb= 6*m/t^2
add membrane stress to bending and check stress limits per Div 2.
For large vales of e , I recommend fea or elastic analysis.
Of course, this does not account for thermal.