1.一种考虑非均匀分布齿根裂纹故障的仿真方法，其特征在于，包括：步骤1：将齿轮沿齿宽方向分解为多个独立的切片齿，基于每个所述切片齿的角位移建立切片齿动态传递误差计算公式；

以及，基于变截面悬臂梁理论建立裂纹参数作用下的切片齿时变啮合刚度计算公式；

基于所述切片齿动态传递误差计算公式与齿根裂纹影响下的切片齿时变啮合刚度计算公式求解切片齿的动态啮合力；并根据所述切片齿的动态啮合力建立基于切片齿的齿轮传动系统动力学方程；

其中，所述切片齿动态传递误差计算公式为：δi＝‑Rbgβgy‑Rbpβpy‑[(Zg‑Liβgx)‑(Zp‑Liβpx)]cosαr‑[(Xp+Liβpz)‑(Xg+Liβgz)]sinαr‑erri其中，δi为第i个切片齿的动态传递误差，Rbg\Rbp为大\小齿轮基圆半径，单位为：m；βg\βp为大\小齿轮角位移，单位为：rad，下标x\y\z分别表示绕该方向的转动；Zg\Zp为大\小齿轮质心垂向位移，单位为：m；Xg\Xp为大/小齿轮质心纵向位移，单位为：m；erri表示第i个切片齿装配制造误差和齿廓修形；Li表示切片齿质心与整个齿轮质心的距离，单位为：m；

步骤2：对当前待仿真齿轮的使用状况进行评估，根据评估结果设定当前待仿真齿轮的裂纹参数，采集当前待仿真齿轮的设计参数，并将所述设计参数与裂纹参数代入所述基于切片齿的齿轮传动系统动力学方程中，求解得到当前待仿真齿轮的动力学响应。

2.根据权利要求1所述的考虑非均匀分布齿根裂纹故障的仿真方法，其特征在于，所述设计参数包括：齿轮模数、齿轮齿数、齿轮压力角、齿宽、弹性模量、泊松比、齿顶高系数、顶隙系数；

所述裂纹参数，包括：裂纹长度、裂纹深度、裂纹交角。

3.根据权利要求1所述的考虑非均匀分布齿根裂纹故障的仿真方法，其特征在于，所述根据多个所述切片齿的动态啮合力建立齿轮系统动力学方程，包括：步骤101，基于所述切片齿的动态啮合力以及每个所述切片齿质心到整个齿轮质心的距离求解每个切片齿的动态啮合力相对于齿轮质心产生的非平衡力矩；

步骤102，基于多个所述切片齿的动态啮合力相对于齿轮质心产生的非平衡力矩，齿轮啮合力和轴承支承力来建立所述基于切片齿的齿轮传动系统动力学方程。

4.根据权利要求1所述的考虑非均匀分布齿根裂纹故障的仿真方法，其特征在于，所述基于不同裂纹深度切片齿相对于齿轮质心产生的侧滚和摇头振动的动力学方程为：其中，Nz(k,1)表示左侧轴承内外圈垂向接触力；Nz(k,2)表示右侧轴承内外圈垂向接触力；Nx(k,2)表示右侧轴承内外圈纵向接触力；Nx(k,1)表示左侧轴承内外圈纵向接触力；Jx\Jz为侧滚和摇头的转动惯量；ns为切片齿数量；k＝1\k＝2分别表示大\小齿轮； 为非均匀啮合力产生的侧滚非平衡力矩； 为非均匀啮合力产生的摇头非平衡力矩。

5.根据权利要求1‑4任一所述的考虑非均匀分布齿根裂纹故障的仿真方法，其特征在于，所述基于变截面悬臂梁理论建立裂纹参数作用下的切片齿时变啮合刚度计算公式，包括：

利用势能原理分别计算裂纹参数作用下的切片齿弯曲刚度、切片齿剪切刚度和切片齿轴向压缩刚度；

基于变截面悬臂梁模型求解裂纹参数作用下的切片齿轮体变形在啮合线上的等效刚度；

根据赫兹接触理论得到的切片齿赫兹接触刚度；

基于所述切片齿弯曲刚度、切片齿剪切刚度和切片齿轴向压缩刚度、齿根裂纹影响下的切片齿轮体变形在啮合线上的等效刚度、切片齿的赫兹接触刚度得到所述切片齿时变啮合刚度计算公式。

6.根据权利要求1所述的考虑非均匀分布齿根裂纹故障的仿真方法，其特征在于，所述步骤2中，采用数值积分法求解得到当前待仿真齿轮的动力学响应。