Contributions of the pre-ionized H2 (PI-H2) and ionized H + 2 subsystems of the two-electron H2 system to its high-order harmonic generation in 8-cycle sin 2 -like ultrafast intense laser pulses are calculated and analyzed based on the solution of the time-dependent Schrödinger equation (TDSE) for the one-dimensional two-electronic H2 system with fixed nuclei. The laser pulses have λ = 390 & 532 nm wavelengths and I = 1 × 10 14 , 5×10 14 , 1×10 15 & 5×10 15 Wcm −2 intensities. It is found that at the two lower intensities, the PI-H2 subsystem dominantly produces the HHG spectra. While, at the two higher intensities, both PI-H2 and ionized H + 2 subsystems contribute comparably to the HHG spectra. In the H + 2 subsystem, the symmetry of the populations of H + 2 (I) and H + 2 (II) regions (left and right regions of H + 2 subsystem) is broken by increasing the laser intensity. Complex patterns and even harmonics also appear at these two higher intensities. For instance, at 1 × 10 15 Wcm −2 intensity and λ = 532 nm wavelength, the even harmonics are appeared near cut-off region. Interestingly, at 5 × 10 15 Wcm −2 intensity and λ = 390 nm wavelength, the even harmonics replaced by the odd harmonics with red shift. At λ = 390 & 532 nm wavelengths and I = 1×10 15 intensity, the two-electron cutoffs corresponding to nonsequential double-recombination (NSDR) with maximum return kinetic energy of 4.70Up are detected. The HHG spectra of the whole H2 system obtained with and without nuclear dynamics treated classically are approximately similar. However, at 1 × 10 15 Wcm −2 intensity and λ = 532 nm wavelength, if we take into account nuclear dynamics, the even harmonics which are appeared near cutoff region, replaced by the odd harmonics with blue shift. *