The junction between cylindrical part and end dome is the most critical place of the filament wound composite pressure vessels because the shear forces and bending moments occur in this area. The strength analyses of junction between cylindrical part and different types of known end domes based on influence coefficients method were performed for filament wound composite pressure vessels manufactured using glass-epoxy and carbon-epoxy composites. Four known end domes were considered that is, a spherical shell, a geodesic-isotensoid shell, a shell based on minimizing of the Tsai-Hill’s strength criterion and shell designed in authors’ previous work (Hoffman 1 shell). Moreover, the study is also dealing with design of optimal end dome shape which takes account to the junction area. Finding of optimal end dome shape was based on data-driven evolutionary algorithms. Using input dataset from the analytical simulation surrogate models were created using Evolutionary Deep Neural Nets (EvoDN2) algorithm. The multicriterial optimization procedure was performed using Constrained Reference Vector Algorithm (cRVEA) optimization algorithm. Results indicates that the optimal end dome shape highly depend on material properties.