For non-destructive evaluation (NDE) of fatigue damage accumulation in composites, this research proposed a combined linear and a nonlinear ultrasonic guided wave method. For the linear Lamb waves approach, a laser-generation based imaging system (LGBI) is utilized to measure the phase velocities of guided waves in composites. The elastic moduli of the specimen are then obtained by inverting the measured phase velocities using genetic algorithms (GAs). The variation of the above two parameters (phase velocity and elastic moduli), together with the guided wave amplitudes, are then observed during the fatigue process. Nonlinear second harmonics in composites are studied theoretically and numerically. A third-order strain energy function of transversely isotropic materials is expressed by five invariants of the Green-Lagrange strain tensor. Results enable intelligent selection of primary modes for cumulative second harmonics generation. Meanwhile, finite element simulations are conducted to characterize second harmonics in light of the theory.