This paper investigates numerically the post-bifurcation evolution of ellipsoidal or spherical balloons subject to internal pressure, where the primary -shaped curves of pressure vs. principal stretch and the corresponding bifurcated branch, i.e. pear-shaped deformation, are captured quantitatively. We quantify and discuss the range of pear-shaped bifurcation intervals of ellipsoids and the associated critical points. For rugby-shaped ellipsoidal balloons, we find that there exists a threshold for the aspect ratio of the major and minor axes that leads to pear-shaped bifurcation, which is detected by the finite element method. When the rugby shape becomes sufficiently slender, the nonlinear response of the ellipsoidal balloons is well approximated by the deformation of a tube with localized bulging instead of pear-shaped configuration. We obtain the nonlinear evolution of the localized bulging of rugby-like ellipsoids numerically. Furthermore, we examine the influence of various aspect ratios for rugby-shaped balloons on the localized bulging response. We find that pumpkinshaped ellipsoids can always bifurcate into pear shape. Lastly, we provide a unified phase diagram on instability mode selection of various aspect ratios of ellipsoidal balloons, with diverse representative deformed configurations.