说明：  基于weber能量法求解齿轮时变啮合刚度的程序，能够跑出刚度图，通过求解轮齿部分变形、基体变形及局部接触变形这三部分的变形，进而求得综合弹性变形，最终求出时变啮合刚度
(The program based on Weber energy method to solve the time-varying meshing stiffness of gears can run out of the stiffness diagram. By solving the deformation of three parts of gear teeth, including partial deformation, matrix deformation and local contact deformation, the comprehensive elastic deformation is obtained, and the time-varying meshing stiffness is finally obtained)