CN113111450B

CN113111450B - 一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法

Info

Publication number: CN113111450B
Application number: CN202110260039.7A
Authority: CN
Inventors: 马洪顺; 房心强; 汪兵兵; 季娟; 李浩峰; 侯臣; 吴静奇
Original assignee: Bosch Huayu Steering Systems Co Ltd
Current assignee: Bosch Huayu Steering Systems Co Ltd
Priority date: 2021-03-10
Filing date: 2021-03-10
Publication date: 2022-07-19
Anticipated expiration: 2041-03-10
Also published as: CN113111450A

Abstract

本发明涉及机械零件强度分析技术领域，具体地说是一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，本发明同现有技术相比，将温度影响纳入疲劳强度分析中，建立纳入温度因子的疲劳破坏寿命计算方法，可准确预测塑料蜗轮在不同温度下的疲劳破坏寿命，为塑料蜗轮设计、疲劳强度校核和制造使用提供了有力的理论和实验支撑，可有效避免塑料蜗轮的早期疲劳失效，减少塑料蜗轮疲劳断裂发生的概率，提高塑料蜗轮运行的安全可靠性。

Description

一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法

技术领域

本发明涉及机械零件强度分析技术领域，具体地说是一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法。

背景技术

塑料蜗轮摩擦系数低，耐磨性好，传动效率高，可在无润滑或少润滑条件下运行。同时，塑料蜗轮具有较好的弹性，能吸震防冲，运行时噪音低，传动平稳。但由于塑料蜗轮硬度低，弹性模量小，疲劳强度不高，因此限制了其使用范围。

现有的塑料蜗轮疲劳强度设计方法一般均采用传统的疲劳S-N曲线就行校核评估，其未能充分纳入温度的影响。但由于塑料蜗轮散热效果较差，在实际运行过程中，温度可能高达100℃以上，较高的使用温度会使塑料达到玻璃化转变温度以上而发生软化，使得塑料蜗轮疲劳强度急剧降低，导致塑料齿轮早期高温疲劳破坏。

因此，需要设计一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，以确保塑料蜗轮在高温运行时的安全稳定性。

发明内容

本发明的目的是克服现有技术的不足，提供了一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，以确保塑料蜗轮在高温运行时的安全稳定性。

为了达到上述目的，本发明提供一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，包括如下步骤：步骤1，控制塑料蜗轮疲劳实验箱的环境温度T，并针对不同环境温度T分别施加设定扭矩P，得到塑料蜗轮材料在不同环境温度下施加不同扭矩所对应的疲劳破坏寿命N；步骤2，建立疲劳破坏寿命S-N公式：N＝N0×(P/P0)^-a(T)×10^{-b(T)×(T/T0-1)}，a_(T)为塑料蜗轮材料温度相关的疲劳破坏寿命指数，b_(T)为塑料蜗轮材料温度相关的疲劳破坏寿命指数；步骤3，确定额定工作扭矩P0和额定环境温度T0；步骤4，建立额定工作扭矩P0和额定环境温度T0所对应的额定疲劳破坏寿命N0计算公式：Log^N0＝a_(T)×Log^P/P0+b_(T)×(T/T0-1)+Log^N；步骤5，计算不同环境温度T下的对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N；步骤6，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的平均值；步骤7，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的偏差平方和；步骤8，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的协方差；步骤9，通过线性回归模型最小二乘法拟合计算塑料蜗轮材料的疲劳破坏寿命系数b；步骤10，通过线性回归模型最小二乘法拟合计算塑料蜗轮材料的疲劳破坏寿命系数a；步骤11，将疲劳破坏寿命系数b、疲劳破坏寿命系数a代入步骤4公式，计算得额定疲劳破坏寿命N0；步骤12，将额定疲劳破坏寿命N0代入步骤2公式，获得另一额定工作温度下的疲劳破坏寿命计算公式；步骤13，使用步骤12获得的公式拟合不同环境温度T条件下塑料蜗轮的疲劳破坏寿命曲线，得到该环境温度T下塑料蜗轮的疲劳破坏寿命指数a_(T)，塑料蜗轮的疲劳破坏寿命指数b_(T)；步骤14；根据步骤13所获得的疲劳破坏寿命指数a_(T)，获得环境温度T相关的疲劳破坏寿命指数a_(T)计算公式为a_(T)＝-5.224-0.009×T；步骤15；根据步骤13所获得的疲劳破坏寿命指数b_(T)，获得环境温度T相关的疲劳破坏寿命指数b_(T)计算公式为b_(T)＝0.06+0.028×T-1.933×T²；步骤16，将步骤14、步骤15的计算公式代入步骤12中另一额定工作温度下的疲劳破坏寿命计算公式，获得疲劳破坏寿命预测公式。

所述的塑料蜗轮疲劳实验箱为带温控箱的塑料蜗轮疲劳实验台架。

所述的额定工作扭矩P0和额定环境温度T0根据塑料蜗轮实际工作条件确定。

所述的Log^N0＝a_(T)×Log^P/P0+b_(T)×(T/T0-1)+Log^N由N＝N0×(P/P0)^-a(T)×10^{-b(T)×(T/T0-1)}双边取对数并移项获得。

所述的步骤6，对数归一化扭矩Log^P/P0平均值计算公式为：E_(X)＝(X₁+X₂+X₃+……X_n)/n；归一化温度T/T0-1⁰平均值计算公式为：E_(Y)＝(Y₁+Y₂+Y₃+……Y_n)/n；对数疲劳破坏寿命Log^N平均值计算公式为：E_(Z)＝(Z₁+Z₂+Z₃+……Z_n)/n。

所述的步骤7，对数归一化扭矩Log^P/P0偏差平方和计算公式为：D_(X)＝S[X_i(i＝1,2,3…n)-E_(X)]²；归一化温度T/T0-1⁰偏差平方和计算公式为：D_(Y)＝S[Y_i(i＝1,2,3…n)-E_(Y)]²；对数疲劳破坏寿命Log^N偏差平方和计算公式为：D_(Z)＝S[Z_i(i＝1,2,3…n)-E_(Z)]²。

所述的步骤8，对数归一化扭矩Log^P/P0协方差计算公式为：Cov(X,Y)＝E[X_{i(i＝1,2,3……n)}Y_{i(i＝1,2,3……n)}]-E[X_{i(i＝1,2,3……n)}]×E[Y_{i(i＝1,2,3……n)}]；归一化温度T/T0-1⁰协方差计算公式为：Cov(X,Z)＝E[X_{i(i＝1,2,3……n)}Z_{i(i＝1,2,3……n)}]-E[X_{i(i＝1,2,3……n)}]×E[Z_{i(i＝1,2,3……n)}]；对数疲劳破坏寿命Log^N协方差计算公式为：Cov(Y,Z)＝E[Y_{i(i＝1,2,3……n)}Z_{i(i＝1,2,3……n)}]-E[Y_{i(i＝1,2,3……n)}]×E[Z_{i(i＝1,2,3……n)}]。

所述的步骤9，疲劳破坏寿命系数b的计算公式为：b＝[Cov(X,Y)×Cov(X,Z)-Cov(Y,Z)×D_(X)]/[D_(X)D_(Y)-Cov(X,Y)]。

所述的步骤10，疲劳破坏寿命系数a的计算公式为：a＝[Cov(X,Y)×Cov(Y,Z)-Cov(X,Z)×D_(Y)]/[D_(X)D_(Y)-Cov(X,Y)]。

本发明同现有技术相比，将温度影响纳入疲劳强度分析中，建立纳入温度因子的疲劳破坏寿命计算方法，可准确预测塑料蜗轮在不同温度下的疲劳破坏寿命，为塑料蜗轮设计、疲劳强度校核和制造使用提供了有力的理论和实验支撑，可有效避免塑料蜗轮的早期疲劳失效，减少塑料蜗轮疲劳断裂发生的概率，提高塑料蜗轮运行的安全可靠性。

附图说明

图1为塑料蜗轮与蜗杆啮合示意图。

图2为不同扭矩与温度下蜗轮疲劳破坏寿命实验数据。

图3为温度T＝40℃时塑料蜗轮在不同扭矩下疲劳破坏寿命实验数据与拟合曲线对比结果。

图4为温度T＝60℃时塑料蜗轮在不同扭矩下疲劳破坏寿命实验数据与拟合曲线对比结果。

图5为温度T＝100℃时塑料蜗轮在不同扭矩下疲劳破坏寿命实验数据与拟合曲线对比结果。

图6为塑料蜗轮材料疲劳破坏寿命系数a与温度T的关系。

图7为塑料蜗轮材料疲劳破坏寿命系数b与温度T的关系。

图8为塑料蜗轮在运行温度T＝80℃时纳入温度因子的疲劳破坏寿命预测曲线与疲劳实验数据对比。

具体实施方式

现结合附图对本发明做进一步描述。

本发明提供一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，包括如下步骤：步骤1，使用带温控箱的塑料蜗轮疲劳实验台架，控制塑料蜗轮疲劳实验箱的环境温度T＝40℃、60℃和100℃，并针对不同环境温度T分别施加设定扭矩P＝60Nm、70Nm、80Nm、90Nm、100Nm、110Nm和120Nm，得到塑料蜗轮材料在不同环境温度下施加不同扭矩所对应的疲劳破坏寿命N；塑料蜗轮与蜗杆啮合示意图参见图1，实验数据结果见图2。

步骤2，建立疲劳破坏寿命S-N公式：N＝N0×(P/P0)^-a(T)×10^{-b(T)×(T/T0-1)}。式中，P0为额定工作扭矩，N0为额定工作扭矩P0和额定环境温度T0所对应的额定疲劳破坏寿命，T0为额定环境温度，a_(T)为塑料蜗轮材料温度相关的疲劳破坏寿命指数，b_(T)为塑料蜗轮材料温度相关的疲劳破坏寿命指数。

步骤3，确定额定工作扭矩P0和额定环境温度T0。额定工作扭矩P0和额定环境温度T0根据塑料蜗轮实际工作条件确定，本例中，P0＝80Nm，T0＝70℃。

步骤4，建立额定工作扭矩P0和额定环境温度T0所对应的额定疲劳破坏寿命N0计算公式：Log^N0＝a_(T)×Log^P/P0+b_(T)×(T/T0-1)+Log^N。

Log^N0＝a_(T)×Log^P/P0+b_(T)×(T/T0-1)+Log^N由N＝N0×(P/P0)^-a(T)×10^{-b(T)×(T/T0-1)}双边取对数并移项获得。

步骤5，计算不同环境温度T下的对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N。

以环境温度T＝40℃时为例，塑料蜗轮的对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N结果如下表。

步骤6，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的平均值。

对数归一化扭矩Log^P/P0平均值计算公式为：E_(X)＝(X₁+X₂+X₃+……X_n)/n；归一化温度T/T0-1⁰平均值计算公式为：E_(Y)＝(Y₁+Y₂+Y₃+……Y_n)/n；对数疲劳破坏寿命Log^N平均值计算公式为：E_(Z)＝(Z₁+Z₂+Z₃+……Z_n)/n。

式中，E_(X)、E_(Y)和E_(Z)表示平均值，X_{i(i＝1,2,3……n)}、Y_{i(i＝1,2,3……n)}和Z_{i(i＝1,2,3……n)}为变量，其中X、Y、Z分别表示为Log^P/P0、T/T0-1和Log^N。

本例中，E_(X)＝0.04，E_(Y)＝-0.05，E_(Z)＝4.40。

步骤7，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的偏差平方和。

对数归一化扭矩Log^P/P0偏差平方和计算公式为：D_(X)＝S[X_i(i＝1,2,3…n)-E_(X)]²；归一化温度T/T0-1⁰偏差平方和计算公式为：D_(Y)＝S[Y_i(i＝1,2,3…n)-E_(Y)]²；对数疲劳破坏寿命Log^N偏差平方和计算公式为：D_(Z)＝S[Z_i(i＝1,2,3…n)-E_(Z)]²。

式中，D_(X)、D_(Y)和D_(Z)为X_{i(i＝1,2,3……n)}、Y_{i(i＝1,2,3……n)}和Z_{i(i＝1,2,3……n)}偏差平方和，E_(X)、E_(Y)和E_(Z)为X_{i(i＝1,2,3……n)}、Y_{i(i＝1,2,3……n)}和Z_{i(i＝1,2,3……n)}的平均值，其中X、Y、Z分别为Log^P/P0、T/T0-1和Log^N。

本例中，D_(X)＝0.42，D_(Y)＝5.33，D_(Z)＝19.98。

步骤8，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的协方差。

对数归一化扭矩Log^P/P0协方差计算公式为：Cov(X,Y)＝E[X_{i(i＝1,2,3……n)}Y_{i(i＝1,2,3……n)}]-E[X_{i(i＝1,2,3……n)}]×E[Y_{i(i＝1,2,3……n)}]；归一化温度T/T0-1⁰协方差计算公式为：Cov(X,Z)＝E[X_{i(i＝1,2,3……n)}Z_{i(i＝1,2,3……n)}]-E[X_{i(i＝1,2,3……n)}]×E[Z_{i(i＝1,2,3……n)}]；对数疲劳破坏寿命Log^N协方差计算公式为：Cov(Y,Z)＝E[Y_{i(i＝1,2,3……n)}Z_{i(i＝1,2,3……n)}]-E[Y_{i(i＝1,2,3……n)}]×E[Z_{i(i＝1,2,3……n)}]。

式中，Cov(X,Y)、Cov(X,Z)和Cov(Y,Z)分别为X～Y、X～Z和Y～Z的协方差，E(XY)、E(XZ)和E(YZ)分别为XY、XZ和YZ的平均值，其中X、Y、Z分别为Log^P/P0、T/T0-1和Log^N。

本例中，Cov(X,Y)＝1.11E-16，Cov(X,Z)＝-2.45，Cov(Y,Z)＝-4.94。

步骤9，通过线性回归模型最小二乘法拟合计算塑料蜗轮材料的疲劳破坏寿命系数b，b＝[Cov(X,Y)×Cov(X,Z)-Cov(Y,Z)×D_(X)]/[D_(X)D_(Y)-Cov(X,Y)]。

式中，Cov(X,Y)、Cov(X,Z)和Cov(Y,Z)分别为X～Y、X～Z和Y～Z的协方差，D_(X)和D_(Y)为X_{i(i＝1,2,3……n)}和Y_{i(i＝1,2,3……n)}偏差平方和，其中X、Y、Z分别表示为Log^P/P0、T/T0-1和Log^N。

本例中，根据步骤7和步骤8的计算结果，可得到塑料蜗轮材料的疲劳破坏寿命系数b＝0.927。

步骤10，通过线性回归模型最小二乘法拟合计算塑料蜗轮材料的疲劳破坏寿命系数a，a＝[Cov(X,Y)×Cov(Y,Z)-Cov(X,Z)×D_(Y)]/[D_(X)D_(Y)-Cov(X,Y)]。

式中Cov(X,Y)、Cov(X,Z)和Cov(Y,Z)分别为X～Y、X～Z和Y～Z的协方差，D_(X)和D_(Y)为X_{i(i＝1,2,3……n)}和Y_{i(i＝1,2,3……n)}偏差平方和，其中X、Y、Z分别为Log^P/P0、T/T0-1和Log^N。

本例中，根据步骤7和步骤8的计算结果，可得到塑料蜗轮材料的疲劳破坏寿命系数a＝5.827。

步骤11，将疲劳破坏寿命系数b、疲劳破坏寿命系数a代入步骤4公式，计算得额定疲劳破坏寿命N0。

本例中，得到额定工作扭矩P0＝80Nm及额定环境温度T0＝80℃时，所对应的额定疲劳破坏寿命N0＝38911。

步骤12，将额定疲劳破坏寿命N0代入步骤2公式，获得另一额定工作温度下的疲劳破坏寿命计算公式。

本例中，得到以额定工作扭矩P0＝80Nm及额定环境温度T0＝70℃时的疲劳破坏寿命计算公式为：N＝38911×(P/80)^-a(T)×10^{-b(T)×(T/70-1)}。

步骤13，使用步骤12获得的公式拟合不同环境温度T条件下塑料蜗轮的疲劳破坏寿命曲线，得到该环境温度T下塑料蜗轮的疲劳破坏寿命指数a_(T)，塑料蜗轮的疲劳破坏寿命指数b_(T)。

本例中，在环境温度T＝40℃条件下，拟合曲线与实验数据对比结果见图3，得到环境温度T＝40℃时，塑料蜗轮的疲劳破坏寿命指数a_(T)为-5.56，塑料蜗轮的疲劳破坏寿命指数b_(T)为0.87。

在环境温度T＝60℃条件下，拟合曲线与实验数据对比结果见图4，得到环境温度T＝60℃时，塑料蜗轮的疲劳破坏寿命指数a_(T)为-5.82，塑料蜗轮的疲劳破坏寿命指数b_(T)为1.04。

在环境温度T＝100℃条件下，拟合曲线与实验数据对比结果见图5，得到环境温度T＝100℃时，塑料蜗轮的疲劳破坏寿命指数a_(T)为-6.11，塑料蜗轮的疲劳破坏寿命指数b_(T)为0.92。

步骤14；根据步骤13所获得的疲劳破坏寿命指数a_(T)，获得环境温度T相关的疲劳破坏寿命指数a_(T)计算公式为a_(T)＝-5.224-0.009×T。

步骤15；根据步骤13所获得的疲劳破坏寿命指数b_(T)，获得环境温度T相关的疲劳破坏寿命指数b_(T)计算公式为b_(T)＝0.06+0.028×T-1.933×T²。

步骤16，将步骤14、步骤15的计算公式代入步骤12中另一额定工作温度下的疲劳破坏寿命计算公式，获得疲劳破坏寿命预测公式：N＝38911×(P/80)^{5.224+0.009×T}×10^{-(0.06+0.028×T}-^{1.93×T^2)×(T/70-1)}。

通过步骤16计算塑料蜗轮环境温度T＝80℃的疲劳破坏寿命，并与实验数据对比结果见图8，结果显示本发明计算得到的疲劳破坏寿命预测曲线与实验数据吻合较好，验证了本发明的可行行和有效性。

Claims

1.一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，其特征在于：包括如下步骤：步骤1，控制塑料蜗轮疲劳实验箱的环境温度T，并针对不同环境温度T分别施加设定扭矩P，得到塑料蜗轮材料在不同环境温度下施加不同扭矩所对应的疲劳破坏寿命N；步骤2，建立疲劳破坏寿命S-N公式：N=N0×(P/P0)^-a(T) ×10^{-b(T)× (T/T0-1)} ，a_(T)为塑料蜗轮材料温度相关的疲劳破坏寿命指数，b_(T)为塑料蜗轮材料温度相关的疲劳破坏寿命指数；步骤3，确定额定工作扭矩P0和额定环境温度T0；步骤4，建立额定工作扭矩P0和额定环境温度T0所对应的额定疲劳破坏寿命N0计算公式：Log^N0= a_(T)×Log^P/P0+b_(T)×(T/T0-1)+ Log^N；步骤5，计算不同环境温度T下的对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N；步骤6，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的平均值；步骤7，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的偏差平方和；步骤8，分别计算对数归一化扭矩Log^P/P0、归一化温度T/T0-1以及对数疲劳破坏寿命Log^N的协方差；步骤9，通过线性回归模型最小二乘法拟合计算塑料蜗轮材料的疲劳破坏寿命系数b；步骤10，通过线性回归模型最小二乘法拟合计算塑料蜗轮材料的疲劳破坏寿命系数a；步骤11，将疲劳破坏寿命系数b、疲劳破坏寿命系数a代入步骤4公式，计算得额定a疲劳破坏寿命N0；步骤12，将额定疲劳破坏寿命N0代入步骤2公式，获得另一额定工作温度下的疲劳破坏寿命计算公式；步骤13，使用步骤12获得的公式拟合不同环境温度T条件下塑料蜗轮的疲劳破坏寿命曲线，得到该环境温度T下塑料蜗轮的疲劳破坏寿命指数a_(T)，塑料蜗轮的疲劳破坏寿命指数b_(T)；步骤14；根据步骤13所获得的疲劳破坏寿命指数a_(T)，获得环境温度T相关的疲劳破坏寿命指数a_(T)计算公式为a_(T)=-5.224 -0.009×T；步骤15；根据步骤13所获得的疲劳破坏寿命指数b_(T)，获得环境温度T相关的疲劳破坏寿命指数b_(T)计算公式为b_(T)=0.06+0.028×T-1.933×T²；步骤16，将步骤14、步骤15的计算公式代入步骤12中另一额定工作温度下的疲劳破坏寿命计算公式，获得疲劳破坏寿命预测公式；所述的步骤9，疲劳破坏寿命系数b的计算公式为：b=[Cov (X,Y) × Cov (X,Z)- Cov (Y,Z)×D_(X)]/[D_(X)D_(Y)- Cov (X,Y)]，Cov(X, Y)、Cov(X, Z)和Cov(Y, Z)分别为X~Y、X~Z和Y~Z的协方差，D_(X)和D_(Y)为X_{i(i=1,2,3……n)}和Y_{i(i=1,2,3……n)}偏差平方和，其中X、Y、Z分别表示为Log^P/P0、T/T0-1和Log^N；步骤10，疲劳破坏寿命系数a的计算公式为：a=[ Cov (X, Y) × Cov (Y,Z)- Cov (X, Z) ×D_(Y)]/[D_(X)D_(Y)- Cov (X, Y)], Cov(X, Y)、Cov(X, Z)和Cov(Y, Z)分别为X~Y、X~Z和Y~Z的协方差，D_(X)和D_(Y)为X_{i(i=1,2,3……n)}和Y_{i(i=1,2,3……n)}偏差平方和，其中X、Y、Z分别为Log^P/P0、T/T0-1和Log^N。

2.根据权利要求1所述的一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，其特征在于：所述的塑料蜗轮疲劳实验箱为带温控箱的塑料蜗轮疲劳实验台架。

3.根据权利要求1所述的一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，其特征在于：所述的额定工作扭矩P0和额定环境温度T0根据塑料蜗轮实际工作条件确定。

4.根据权利要求1所述的一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，其特征在于：所述的Log^N0= a_(T)×Log^P/P0+b_(T)×(T/T0-1)+ Log^N由N=N0×(P/P0)^-a(T) ×10^{-b(T)× (T/T0-1)}双边取对数并移项获得。

5.根据权利要求1所述的一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，其特征在于：所述的步骤6，对数归一化扭矩Log^P/P0平均值计算公式为：E_(X) =(X₁+X₂+X₃+……X_n)/n；归一化温度T/T0-1⁰平均值计算公式为：E_(Y) =(Y₁+Y₂+Y₃+……Y_n)/n；对数疲劳破坏寿命Log^N平均值计算公式为：E_(Z) =(Z₁+Z₂+Z₃+……Z_n)/n。

6.根据权利要求1所述的一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，其特征在于：所述的步骤7，对数归一化扭矩Log^P/P0偏差平方和计算公式为：D_(X)=S[X_i(i=1,2,3…n)-E_(X)]²；归一化温度T/T0-1⁰偏差平方和计算公式为：D_(Y)=S[Y_i(i=1,2,3…n)-E_(Y)]²；对数疲劳破坏寿命Log^N偏差平方和计算公式为：D_(Z)=S[Z_i(i=1,2,3…n)- E_(Z)]²。

7.根据权利要求1所述的一种纳入温度因子的塑料蜗轮疲劳破坏寿命预测方法，其特征在于：所述的步骤8，对数归一化扭矩Log^P/P0协方差计算公式为：Cov(X, Y)=E[X_{i(i=1,2,3……n)}Y_{i(i=1,2,3……n)}]-E[X_{i(i=1,2,3……n)}]×E[Y_{i(i=1,2,3……n)}]；归一化温度T/T0-1⁰协方差计算公式为：Cov(X, Z)=E[X_{i(i=1,2,3……n)}Z_{i(i=1,2,3……n)}]-E[X_{i(i=1,2,3……n)}]×E[Z_{i(i=1,2,3……n)}]；对数疲劳破坏寿命Log^N协方差计算公式为：Cov(Y, Z)=E[Y_{i(i=1,2,3……n)}Z_{i(i=1,2,3……n)}]-E[Y_{i(i=1,2,3……n)}]×E[Z_{i(i=1,2,3……n)}]。