本发明公开了一种功能梯度变厚度叶片模型的颤振分析方法，具体包括以下步骤：构建功能梯度变厚度叶片的模型，提取NURBS的基本参数；通过几何细化步骤，得到新的控制点和节点向量信息，并以此构建颤振分析的插值基函数；用上述的基函数来描述变厚度叶片的颤振位移变量，通过一阶活塞理论，精化板理论和相应边界条件，建立其能量泛函；计算单个参数单元的刚度矩阵和质量矩阵；以此为基础循环形成整体刚度和质量矩阵，得到变厚度叶片的颤振特性结果。本发明利用等几何方法，精化板理论和一阶活塞理论的结合，利用较少的控制点描述复杂几何形状，在保证变厚度叶片几何精确性的同时，节省了计算成本，提高了分析的效率。