1.一种行星滚柱丝杠副的不确定性分析方法，其特征在于，包括以下步骤：S1、通过结构参数、材料性能及边界条件求解行星滚柱丝杠副螺纹牙载荷分布和接触特性，考虑结构参数、材料性能及边界条件中存在的不确定性，基于S‑N曲线和应力‑强度干涉理论构建行星滚柱丝杠副螺纹牙接触疲劳失效模式下的极限状态函数；

步骤S1中，通过结构参数、材料性能及边界条件求解行星滚柱丝杠副螺纹牙载荷分布和接触特性，包括：采用迭代算法求解行星滚柱丝杠副螺纹牙载荷分布模型，如下：其中，PS,PR和PN分别为丝杠、滚柱和螺母的螺距，τ为滚柱参与接触的螺纹牙数，z为滚柱个数，螺纹从丝杠固定端到自由端依次编号，i＝1,2,…,τ，F为轴向外载荷，FSRi为滚柱与丝杠接触侧第i对螺纹牙上承受的轴向载荷，FNRi为滚柱与螺母接触侧第i对螺纹牙上承受的轴向载荷； 为丝杠轴段刚度， 为滚柱轴段刚度， 为螺母轴段刚度， 为丝杠螺纹牙刚度， 为滚柱螺纹牙刚度， 为螺母螺纹牙刚度；

求解丝杠与滚柱的接触刚度，如下：

其中， 为丝杠与滚柱的接触刚度； 为丝杠与滚柱第i对螺纹牙的弹性接触变形；

求解螺母与滚柱螺纹副接触刚度，如下：

其中， 为螺母与滚柱螺纹副接触刚度， 为螺母与滚柱第i对螺纹牙的弹性接触变形；

S2、引入确定性的低偏差点集Halton序列，采用试验设计方法分析不确定性因素对行星滚柱丝杠副螺纹牙承载能力及接触性能影响的敏感性和贡献程度，将重要度排序靠前的参数作为随机变量来降低极限状态功能函数的维数，得到简化后的极限状态函数；

S3、引入期望风险学习函数构建行星滚柱丝杠副主动学习代理模型，采用拟蒙特卡洛方法产生抽样样本点，调用所建立的主动学习代理模型计算每组样本点对应的简化后的极限状态函数，对行星滚柱丝杠副进行不确定性分析，得到行星滚柱丝杠副在任一工况下的失效概率及可靠性灵敏度。

2.如权利要求1所述的行星滚柱丝杠副的不确定性分析方法，其特征在于，丝杠与滚柱第i对螺纹牙的弹性接触变形 的计算公式如下：其中，K(e)和L(e)是第一类和第二类完全椭圆积分，a为接触椭圆的长半轴，b为接触椭圆的短半轴， 为接触椭圆的偏心率，ke＝b/a，为等效弹性模量，ES和vS为丝杠的弹性模量和泊松比，ER和vR为滚柱的弹性模量和泊松比，θ和λ为接触角和滚柱的螺旋升角。

3.如权利要求1所述的行星滚柱丝杠副的不确定性分析方法，其特征在于，基于S‑N曲线和应力‑强度干涉理论构建行星滚柱丝杠副螺纹牙接触疲劳失效模式下的极限状态函数，如下：其中，x＝(x1,x2,…xi…,xn)为影响极限状态函数的随机变量，xi为第i个随机变量，n为随机变量个数，σHlim为接触疲劳极限，σSRi为丝杠与滚柱第i对螺纹牙上的接触应力， 为滚柱与丝杠接触侧上所有螺纹牙之间的最大接触应力，σNRi为螺母与滚柱第i对螺纹牙上的接触应力， 为滚柱与螺母接触侧上所有螺纹牙之间的最大接触应力；FSRi为滚柱与丝杠接触侧第i对螺纹牙上承受的轴向载荷，FNRi为滚柱与螺母接触侧第i对螺纹牙上承受的轴向载荷；a为接触椭圆的长半轴，b为接触椭圆的短半轴，θ和λ为接触角和滚柱的螺旋升角。

4.如权利要求3所述的行星滚柱丝杠副的不确定性分析方法，其特征在于，步骤S2包括：S21、对于行星滚柱丝杠副的每个随机变量xi在[0,1]的区间内生成Halton序列，由Halton序列得到输入参数矩阵为X＝(x1,x2,…,xi,…,xn)，其中xi＝(xi1,xi2,…,xij,…,TxiN) 为xi的样本向量；

S22、响应量包括：滚柱两接触侧的最大载荷分布系数：最大接触应力 和螺纹

副上的局部接触应力σSRi,σNRi，通过计算得到响应矩阵为：Y＝(y1,y2,…,yk,…,ym)，其中ykT＝(yk1,yk2,…,ykj,…,ykN) 为第k个响应yk的样本向量；

S23、将矩阵X和Y中的数据在[‑1,1]范围内进行归一化，再用多项式响应面法将这些数据进行拟合：S24、将多项式系数 以百分比的形式来表示，可以反映每个输入参数对第k个响应量的影响程度；

S25、选择重要度排序靠前的参数来降低随机变量的维数，得到简化的行星滚柱丝杠副的极限状态函数。

5.如权利要求4所述的行星滚柱丝杠副的不确定性分析方法，其特征在于，步骤S3中包括：S31、在不确定域中随机生成N0＝20个样本，计算简化后的极限状态功能函数值，引入期望风险学习函数构建行星滚柱丝杠副主动学习代理模型，首先构造初始的代理模型，如下：T ‑1

u＝FR r‑f(x)

T

其中，Y＝[g(x1),g(x2),…,g(xn)] 为具有n个样本点的真实响应函数，f(x)＝[f1(x),T Tf2(x),…fi(x),…,fp(x)]为回归多项式基函数向量，β＝[β1,β2,…βi,…,βp]为回归系数T向量，p为回归多项式的个数，r为待测点与样本点之间的相关向量函数，F＝[f(x1) ,f(x2T T T) ,…,f(xn) ]为n×p阶扩展矩阵，R为n×n阶对称正定相关矩阵， 为近似极限状态函数g(x)的估计值， 为 的方差，g(x)计算如下：2

其中，z(x)为服从N(0,σ)正态分布的随机过程，则任意两个样本点w和x之间的协方差为：2

cov[z(w),z(x)]＝σR(θ,w,x)

T

其中，θ＝[θ1,θ2,…,θn] 为相关性参数，R(θ,w,x)为表示变量θ,w,x相关性的函数且可以利用连续可微的高斯相关性函数来描述为：2

其中，回归系数β及方差σ可以用广义最小二乘回归方法表示为：T ‑1 ‑1 T ‑1

β＝(FR F) FR Y

5

S32、在不确定性域中生成Nc＝10个候选样本；

S33、根据初始代理模型计算Nc个样本点的估计值 和期望风险的估计值取ERF值最大的点为训练点x*， 计算方法如下：

其中，sign(x)为符号函数，当x＞0时，sign(x)＝1，当x＜0时，sign(x)＝‑1，Φ()为标准正态分布的概率密度函数，Ψ()为标准正态分布的累积分布函数；

‑5

S34、设置阈值或收敛条件为ε＝10 ；如果max(ERF)>ε，增加x*作为新的训练点，更新行星滚柱丝杠副主动学习代理模型；

S35、重复步骤S33～S35，直到max(ERF)≤ε；

S36、利用更新后的星滚柱丝杠副主动学习代理模型，采用拟蒙特卡洛方法产生抽样样本点，调用所建立的主动学习代理模型计算每组样本点对应的简化后的极限状态函数，对行星滚柱丝杠副进行不确定性分析，得到行星滚柱丝杠副在任一工况下的失效概率及可靠性灵敏度。

6.如权利要求5所述的行星滚柱丝杠副的不确定性分析方法，其特征在于，所述采用拟蒙特卡洛方法产生抽样样本点，调用所建立的主动学习代理模型计算每组样本点对应的简化后的极限状态函数，对行星滚柱丝杠副进行不确定性分析，得到行星滚柱丝杠副在任一工况下的失效概率及可靠性灵敏度，包括：S361、对于行星滚柱丝杠副的每个随机变量xi在[0,1]的区间内生成Halton序列，各变量相互独立且均为正态分布，随机变量的累积概率分布函数F(xij)可表示为：其中，xij为第i个随机变量的第j个样本值，f(xi)为随机变量xi的概率密度函数，N为抽样样本数， 为第xij对应的Halton数，其中Halton序列的具体生成过程为：若q为任意质数，则任意自然数j有唯一的q进制表达式：

2 K

j＝j0+j1q+j2q+···+jKq

K K+1

q≤j＜q ；

ji∈{0,1,···,q‑1}；i＝0,1,···,K其中，K表示lnj/lnq的整数部分，以q为底的基逆函数定义为：对于任意自然数j＞0，均满足 若前n个质数为q1,q2,…,qn，则n维Halton序列可以表示为：

S362、获得随机变量xij的样本值

S363、获得许用接触应力的样本点

T

S364、计算每组随机向量的样本xj＝(x1j,x2j,…,xnj) 作用下的载荷分布和局部接触特性，得到滚柱与丝杠或螺母接触侧上的最大接触应力，然后进一步计算极限状态函数g(xj)，如果g(xj)≤0，状态指示函数IF(xj)＝1，否则IF(xj)＝0；

S365、失效概率 的估计及其变异系数 计算如下：

其中，Nf为失效样本数；

S366、可靠性灵敏度的估计值如下：

其中，fx(x)为随机变量的联合概率密度函数， 为基本随机变量xi的第k个分布参数，mi为第i个随机变量xi的分布参数的数量，正态分布随机变量具有均值 和标准差 两个分布参数，可靠性灵敏度可以计算为：

7.一种计算机可读存储介质，其特征在于，所述存储介质包括存储的程序，其中，所述程序执行如权利要求1‑6任意一项所述的行星滚柱丝杠副的不确定性分析方法。

8.一种电子设备，其特征在于，包括：一个或多个处理器，存储器以及一个或多个程序，其中，所述一个或多个程序被存储在所述存储器中，并且被配置为由所述一个或多个处理器执行，所述一个或多个程序包括用于执行上述如权利要求1‑6任意一项所述的行星滚柱丝杠副的不确定性分析方法。

9.一种行星滚柱丝杠副的不确定性分析系统，其特征在于，包括以下模块：极限状态函数构建模块，用于通过结构参数、材料性能及边界条件求解行星滚柱丝杠副螺纹牙载荷分布和接触特性，考虑结构参数、材料性能及边界条件中存在的不确定性，基于S‑N曲线和应力‑强度干涉理论构建行星滚柱丝杠副螺纹牙接触疲劳失效模式下的极限状态函数；

通过结构参数、材料性能及边界条件求解行星滚柱丝杠副螺纹牙载荷分布和接触特性，包括：采用迭代算法求解行星滚柱丝杠副螺纹牙载荷分布模型，如下：其中，PS,PR和PN分别为丝杠、滚柱和螺母的螺距，τ为滚柱参与接触的螺纹牙数，z为滚柱个数，螺纹从丝杠固定端到自由端依次编号，i＝1,2,…,τ，F为轴向外载荷，FSRi为滚柱与丝杠接触侧第i对螺纹牙上承受的轴向载荷，FNRi为滚柱与螺母接触侧第i对螺纹牙上承受的轴向载荷； 为丝杠轴段刚度， 为滚柱轴段刚度， 为螺母轴段刚度， 为丝杠螺纹牙刚度， 为滚柱螺纹牙刚度， 为螺母螺纹牙刚度；

求解丝杠与滚柱的接触刚度，如下：

其中， 为丝杠与滚柱的接触刚度； 为丝杠与滚柱第i对螺纹牙的弹性接触变形；

求解螺母与滚柱螺纹副接触刚度，如下：

其中， 为螺母与滚柱螺纹副接触刚度， 为螺母与滚柱第i对螺纹牙的弹性接触变形；

极限状态函数简化模块，用于引入确定性的低偏差点集Halton序列，采用试验设计方法分析不确定性因素对行星滚柱丝杠副螺纹牙承载能力及接触性能影响的敏感性和贡献程度，将重要度排序靠前的参数作为随机变量来降低极限状态功能函数的维数，得到简化后的极限状态函数；

不确定性分析模块，用于引入期望风险学习函数构建行星滚柱丝杠副主动学习代理模型，采用拟蒙特卡洛方法产生抽样样本点，调用所建立的主动学习代理模型计算每组样本点对应的简化后的极限状态函数，对行星滚柱丝杠副进行不确定性分析，得到行星滚柱丝杠副在任一工况下的失效概率及可靠性灵敏度。