We present a method that simplifies the refinement of force constants in normal mode calculations and makes the results more reliable. The method avoids normal mode crossings by constraining the force constants during refinement. It was tested with pyrrole, imidazole, benzene, pyridine, pyrimidine, aniline and adenine using a Urey-Bradley force field. The global error of the frequency fit for these molecules was 0.61%. The method reproduced with fewer parameters the accuracy of similar calculations of the single ring aromatic compounds. It improved the accuracy and isotopic shifts of previous empirical calculations of adenine by 40%. The C-C and C-N stretchings differed by less than 7% from the values of force constant-bond length empirical relations.