We investigate the nonlinear growth stages of the bending instability in stellar disks with exponential radial density profiles. We found that the unstable modes are global (the wavelengths are larger than the disk scale lengths) and that the instability saturation level is much higher than that following from a linear criterion. The instability saturation time scales are of the order of one billion years or more. For this reason, the bending instability can play an important role in the secular heating of a stellar disk in the z direction. In an extensive series of numerical N-body simulations with a high spatial resolution, we were able to scan in detail the space of key parameters (the initial disk thickness z_0, the Toomre parameter Q, and the ratio of dark halo mass to disk mass M_h/M_d). We revealed three distinct mechanisms of disk heating in the z direction: bending instability of the entire disk, bending instability of the bar, and heating on vertical inhomogeneities in the distribution of stellar ma