When ultrasound propagates to a crack, transmitted and reflected waves are generated. These waves have useful information for the detection of the crack lying in a structure. When a crack is under residual stress, crack surfaces will contact each other and a closed crack is formed. For closed cracks, the fundamental component of the reflected and transmitted waves will be weak, and as such it is not easy to detect them. In this case, higher harmonic components will be useful. In this paper, nonlinear characteristic of a closed crack is modeled by a continuum material having a tensile-compressive unsymmetry, and the amplitude of the second harmonic wave was obtained by spectrum analysis. Variation of the second harmonic component depending on the nonlinearity of the inclusion was investigated. Two-dimensional plane strain model is considered, and finite element software ABAQUS/Explicit is used.