When a crack exists under a residual stress, for example in welds, the crack can be closed and it shows non symmetric behavior for tension and compression. Ultrasonic detection method for those nonlinear cracks has been developed recently. The method uses the higher order harmonics generating at the crack surface. In this study, parameter study was carried out for the analysis of the harmonics generation at a nonlinear contact interface as a preliminary study for general 3-dimensional cracks. One-dimensional problem with simple bilinear behavior for the contacting surface was considered. The amplitude of second harmonic to the fundamental wave was obtained for various stiffness ratios, incident frequencies, and the contacting layer thicknesses.