The invention provides a CPU+GPU heterogeneous parallel computing based natural frequency characteristic analysis method for a turbomachinery blade. The method comprises the steps of firstly, establishing a finite element model according to a three-dimensional model and material parameters of a to-be-analyzed turbomachinery blade, and performing pre-stress analysis on the blade; secondly, performing blade grid data preprocessing, simultaneously calculating a rigidity matrix and a mass matrix of each unit in a CPU and a GPU, and assembling the rigidity matrix and the mass matrix into total rigidity and mass matrixes; setting constraint conditions of the blade and a rim, wherein the setting process comprises constraint of rigid and elastic displacement of a boundary, contact coupling of a blade root and the rim or contact coupling of connecting pieces and correction of the total rigidity matrix; thirdly, extracting generalized eigenvalues and eigenvectors of the total rigidity and mass matrixes by using a CPU+GPU heterogeneous parallel algorithm; fourthly, converting the eigenvalues and eigenvectors into a frequency and a vibration mode of the blade, and performing output; and finally, judging a vibration type of a natural vibration mode, and drawing frequency curve distribution, a vibration safety graph or a Campbell chart of the blade according to the vibration type.