Cracks in solids can be detected by ultrasound if they are open. However, their detection is not easy when they are closed with a closure stress, and it is a fundamental problem in ultrasonic testing. Subharmonics with half the input frequency is potentially useful in the detection and evaluation of such cracks. We developed analytical and numerical theories accounting for the crack parameters, such as closure stress and crack surface conditions, for the first time and proved their validity by comparison with experiments on a well-defined fatigue crack in aluminum alloy. Based on these theories, we propose a novel method to estimate size of partially closed cracks, which solves the fundamental problem in ultrasonic testing.