Abstract: Based on the single-particle theory, we inspect the effects of harmonic wiggler field components on the outputs of free electron lasers (including the components of fundamental and harmonic frequecies), and derive the complete simulation equations for millimeter wave FELs. Meanwhile under the condition that the harmonic wiggler fields are much smaller than the fundamental one, we describe the effects simulated. Using the improved 3D WAGFEL code.