1.一种铁垫板下组合刚度的计算方法，其特征在于，所述计算方法包括如下步骤：S1.建立铁垫板板下组合刚度计算模型，所述铁垫板板下组合刚度计算模型基于铁垫板各位置处所承受的荷载大小和变形状态将铁垫板板下组合刚度计算模型划分为多个非均匀计算单元，每个所述非均匀计算单元的支承刚度与所述非均匀计算单元某一点处的位移值相关；其中，每个所述非均匀计算单元满足Winkler假定，相邻两个所述非均匀计算单元之间满足变形协调条件；S2.获取待计算铁垫板的工况参数以及荷载，将所述待计算铁垫板的工况参数代入所述铁垫板板下组合刚度计算模型中，采用中点刚度法，将荷载分成若干荷载步施加到所述铁垫板板下组合刚度计算模型的多个非均匀计算单元上，得到所有荷载步施加后的总形变；并通过所述总形变和所述荷载得到所述待计算铁垫板的板下组合刚度；所述变形协调条件为；y j+1 L ＝y j R ，θ j+1 L ＝θ j R ，M j+1 L ＝M j R ，Q j+1 L ＝Q j R ；其中，y j R ，θ j R ，M j R ，Q j R 分别表示第j个非均匀计算单元与第j+1个非均匀计算单元相邻位置的位移、转角、弯矩和剪力；y j+1 L ，θ j+1 L ，M j+1 L ，Q j+1 L 分别表示第j+1个非均匀计算单元与第j个非均匀计算单元相邻位置的位移、转角、弯矩、剪力。
2.根据权利要求1所述的计算方法，其特征在于，所述Winkler假定为：q(y j )＝k j y j ；其中，q(y j )为第j个非均匀计算单元板下垫板提供的反力；k j 为第j个非均匀计算单元对应地基系数；y j 为第j个非均匀计算单元的垂向位移。
3.根据权利要求1所述的计算方法，其特征在于，所述非均匀计算在无载条件下的变形曲线为 其中 为对应第j个非均匀计算单元的克雷洛夫函数， 并叠加修正项： 其中x 1j 、x 2j 为第j个非均匀计算单元上荷载作用的起始位置与终止位置，P j 为第j个非均匀计算单元所承受的荷载；对二者求导得到转角方程、弯矩方程和剪力方程为：
4.根据权利要求1所述的计算方法，其特征在于，所述铁垫板工况参数包括：铁垫板长度、铁垫板厚度、铁垫板宽度、螺栓紧固力作用范围和列车荷载作用宽度。
5.根据权利要求1所述的计算方法，其特征在于，所述步骤S2包括：当第i个荷载步加载时，根据加载步的初始位移y i，0 计算第i步对应的初始刚度k i，0 ；将初始刚度k i，0 代入Winker假定中得到在该荷载步一半荷载作用下的形变Δy i，Δ/2 ，再将Δy i，Δ/2 叠加至该荷载步的初始位移中获得该荷载步一半荷载作用下对应的形变y i，Δ/2 ，并根据y i，Δ/2 与刚度-位移曲线计算此时的中点刚度K i，Δ/2 ；将中点刚度k i，Δ/2 作为铁垫板在此荷载步下的刚度k i，Δ/2 ，计算铁垫板在此荷载步下的形变Δy i 并叠加至上一步的总形变中，最终得到所有荷载步下铁垫板下组合形变。
6.一种可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行实现上述权利要求1至5中任一项所述的方法。
7.一种铁垫板下组合刚度计算装置，其特征在于，包括至少一个处理器，以及与所述至少一个处理器通信连接的存储器；所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1至5中任一项所述的方法。
