One of the WRC bulletins will do the job, better if you can find some software.
As I recall, WRC-107 will calculate stresses in the main cylinder only, while WRC-297 does the nozzle as well. Can anybody confirm or deny?
WRC-107 does not use X,Y,Z coordinates, but uses vessel circumferential and longitudinal and nozzle torsional coordinates, but the X,Y,Z loads can be placed in the vessel coordinate system fairly easily. Not sure about 297.
Good luck,
Mike