Optimal control simulation is applied to numerically design nonresonant laser pulses that maximize the degrees of three-dimensional (3D) alignment of SO2 using the lowest-order induced-dipole interaction. In our trials, combinations of more than two mutually orthogonal, linearly polarized subpulses are always obtained as the optimal solutions. Each subpulse in the optimal pulses impulsively excites the rotational wave packet. The optimal pulses effectively cooperate with the rotational dynamics up to only a few partial revival timings owing to the rotational dephasing that determines the effective control periods. The control mechanisms are interpreted in terms of the time derivatives of the expectation values of the squares of the direction cosines, that of the rotational energy, and the interplay between them. We find a special and important role of the last subpulses as they align the molecular axes using the interaction through the two smallest polarizability components, while the other subpulses excite the rotational wave packet mainly through the largest polarizability component. The control pulses composed of the specified number of subpulses are also numerically optimized by actively utilizing the instantaneous penalty to systematically show the superiority of the use of more than two subpulses over that of two subpulses, the latter of which leads to the saturation of the degree of 3D alignment as a function of total fluence.