Nutritional rickets is a worldwide public health problem; however, the current diagnostic methods retain shortcomings for accurate diagnosis of nutritional rickets. To identify urinary biomarkers associated with nutritional rickets and establish a noninvasive diagnosis method, urinary metabonomics analysis by ultra-performance liquid chromatography/quadrupole time-of-flight tandem mass spectrometry and multivariate statistical analysis were employed to investigate the metabolic alterations associated with nutritional rickets in 200 children with or without nutritional rickets. The pathophysiological changes and pathogenesis of nutritional rickets were illustrated by the identified biomarkers. By urinary metabolic profiling, 31 biomarkers of nutritional rickets were identified and five candidate biomarkers for clinical diagnosis were screened and identified by quantitative analysis and receiver operating curve analysis. Urinary levels of five candidate biomarkers were measured using mass spectrometry or commercial kits. In the validation step, the combination of phosphate and sebacic acid was able to give a noninvasive and accurate diagnostic with high sensitivity (94.0%) and specificity (71.2%). Furthermore, on the basis of the pathway analysis of biomarkers, our urinary metabonomics analysis gives new insight into the pathogenesis and pathophysiology of nutritional rickets.