In this work, the stress distribution along cylindrical composite pressure vessels with different dome geometries is investigated. The dome contours are generated through an integral method based on shell stresses. Here, the influence of each dome contour on the stress distribution at the interface of the dome-cylinder is evaluated. At first, the integral formulation for dome curve generation is presented and solved for the different dome contours. An analytical approach for the calculation of the secondary stresses in a cylindrical pressure vessel is introduced. For the analysis, three different cases were investigated: (i) a polymer liner; (ii) a single layer of carbon-epoxy composite wrapped on a polymer liner; and (iii) multilayer carbon-epoxy pressure vessel. Accounting for nonlinear geometry is seen to have an effect on the stress distribution on the pressure vessel, also on the isotropic liner. Significant secondary stresses were observed at the dome-cylinder interface and they reach a maximum at a specific distance from the interface. A discussion on the trend in these stresses is presented. The numerical results are compared with the experimental results of the multilayer pressure vessel. It is observed that the secondary stresses present in the vicinity of the dome-cylinder interface has a significant effect on the failure mechanism, especially for thick walled cylindrical composite pressure vessel. It is critical that these secondary stresses are directly accounted for in the initial design phase.