一种基于混沌多项式展开的可靠性全局灵敏度分析方法
【专利摘要】本发明公开了一种基于混沌多项式展开的可靠性全局灵敏度分析方法，包括如下步骤：S1.根据待分析模型的输入参数和输出参数，确定关键输入参数与关键输出参数；S2.构建用于表征所述关键输入参数与关键输出参数之间函数关系的满足预设条件的混沌多项式展开；S3.计算所述混沌多项式展开的可靠度数值；S4.根据所述可靠度数值，构建可靠性混沌多项式展开；S5.基于Sobol’指标，解析计算可靠性全局灵敏度指标。本发明具有可大幅减少全局灵敏度分析的计算量，分析结果精确度高，更加符合工程实际，可体现不同输入参数间对全局可靠性交叉影响等优点。
【专利说明】
一种基于混沌多项式展开的可靠性全局灵敏度分析方法
技术领域
[0001] 本发明涉及一种可靠性灵敏度分析方法，尤其涉及一种基于混沌多项式展开的可 靠性全局灵敏度分析方法。
【背景技术】
[0002] 灵敏度分析研究的是模型输出受输入参数变化的影响。由于其预测性和诊断性， 通常将其作为建模及模型分析的首要条件。可靠性灵敏度分析将可靠性作为研究对象，主 要分析模型输入变量分布参数的变化引起失效概率变化的程度，借助可靠性灵敏度分析可 以找到对可靠性影响较大/小的因素，从而为可靠性建模分析、参数识别、可靠性优化设计 等工作提供支持。
[0003] 传统的可靠性灵敏度分析，计算的仅是在输入参数的均值点或者标准差的微小变 动对可靠性的影响，是一种典型的局部灵敏度分析方法。具有下列局限性：（1)无法探索输 入参数的整个取值空间对可靠性的影响，从而无法找到输入参数的最佳变化区域；（2)在某 参数的概率分布变化范围内，各个点处的偏导数是不同的，甚至差距甚大，仅选择定义域内 某个点处的偏导数来作为灵敏度判据是不恰当的；（3)在对某一参数进行可靠性灵敏度指 标计算的时候，需假定其余参数为定值，无法考虑输入参数同时变化的情形，从而不能研究 各输入参数的交叉作用对可靠性的影响，无法找到影响可靠性的风险因子，该风险因子的 不确定性会掩盖其他参数对可靠性的影响，导致其他参数难以识别。
【发明内容】
[0004] 本发明要解决的技术问题就在于:针对现有技术存在的技术问题，本发明提供一 种可大幅减少全局灵敏度分析的计算量，分析结果精确度高，更加符合工程实际，可体现不 同输入参数间对全局可靠性交叉影响的基于混沌多项式展开(PCE)的可靠性全局灵敏度分 析方法。
[0005] 为解决上述技术问题，本发明提出的技术方案为：一种基于混沌多项式展开的可 靠性全局灵敏度分析方法，包括如下步骤：
[0006] S1.根据待分析模型的输入参数和输出参数，确定关键输入参数与关键输出参数；
[0007] S2.构建用于表征所述关键输入参数与关键输出参数之间函数关系的满足预设条 件的混沌多项式展开；
[0008] S3.计算所述混沌多项式展开的可靠度数值；
[0009] S4.根据所述可靠度数值，构建可靠性混沌多项式展开；
[0010] S5.基于Sobol'指标，解析计算可靠性全局灵敏度指标。
[0011] 作为本发明的进一步改进，所述步骤S1中确定所述关键输入参数与关键输出参数 的步骤包括：
[0012] SI. 1.计算所述输入参数和输出参数的概率分布、随机变量特性和随机过程特性；
[0013] S1.2.判断是否满足任意一个预设的判定准则，是则判定所述输入参数为关键输 入参数，所述输出参数为关键输出参数；
[0014] 所述预设的判定准则包括：
[0015] T1.所述概率分布是否满足预设的分布条件；
[0016] T2.所述随机变量特性是否满足预设的特性阈值；
[0017] T3.所述随机过程特性是否满足预设的过程特性阈值。
[0018] 作为本发明的进一步改进，所述步骤S2的具体步骤包括：
[0019] S2.1.根据所述关键输入参数的概率密度函数分布，结合Wiener-Askey方案，确定 关键输入参数的基底类型，所述基底类型为标准随机变量的函数；
[0020] S2.2.将所述关键输入参数与关键输出参数表征为所述基底类型的混沌多项式展 开，所述混沌多项式展开的阶数为n，n多2，n的初始值为2,如式(1)所示：
[0021]
[0022] 式（1)中，y(η)为混纯多项式展开，η为混纯多项式展开的阶，AA,.，％，a均为 混沌多项式展开的系数，么)，Γ成，《)，Γ,私，'…尤）均为所确定的基底类型， 式，4，· ?均为标准随机变量；
[0023] S2.3.从所述关键输入参数中随机采样，生成第一关键输入参数集，将所述第一关 键输入参数集中的关键输入参数表示为所述标准随机变量的转换函数，将n+1阶混沌多项 式展开所确定的基底类型的根作为随机配点输入至所述转换函数，计算得到第一输入参 数;所述第一关键输入参数集中元素个数至少为所述η阶混沌多项式展开中系数个数的2 倍；
[0024] S2.4.将所述第一输入参数输入预设的第一蒙特卡洛仿真(MCS)模型，计算得到与 所述第一输入参数对应的第一输出参数，生成由第一输入参数与第一输出参数构成的第一 样本数据；
[0025] S2.5.计算所述第一样本数据的条件数，并判断所述条件数是否小于预设的条件 数阈值，是则跳转到步骤S2.6;否则跳转到步骤S2.3;
[0026] S2.6.根据所述第一样本数据，通过改进的概率配点法和回归分析法计算如式（1) 所示的η阶和n+1阶混沌多项式展开的系数，计算η阶混沌多项式展开的值、n+1阶混沌多项 式展开的值，并判断所述η阶混沌多项式展开的值、n+1阶混沌多项式展开的值和所述第一 输出参数值之间的误差是否小于预设的误差阈值，是则确定所述混沌多项式展开的阶为n， 否则，将所述混沌多项式展开的阶数加1，跳转至步骤S2.2。
[0027] 作为本发明的进一步改进，所述步骤S3包括如下步骤：
[0028] S3.1.对所述关键输入参数采用内外表直积法进行水平组合，得到2倍于由所述步 骤S2.6中所确定的η阶混沌多项式展开中包含的未知系数个数的组合数；
[0029] S3.2.计算所述组合数中各关键输入参数的水平偏离中心值的大小Ah，对于每一 个组合数，将所述如式(1)所示的混沌多项式展开变形为如式(2)所示的形式；
[0030]
[0031 ] 式⑵中，y (η)为混纯多项式展开，η为混纯多项式展开的阶，％，*3/,,均为混纯 多项式展开的系数，Γι(?, +Δ~),Γ2((< +Δ气)，(戋 +ΔΑ,))，Γη((裘 +Δ/\),钱.+Δ/?;>??成 ^ 均为变形后的基底类型，么，彖，…，'均为标准随机变量，，Δ/??:，…均为各关键输入参 数的水平偏离中心值的大小；
[0032] 将所述式(2)展开合并，变形为如式(3)所示形式：
[0033]
[0034] 式(3)中，y (η)为混纯多项式展开，η为混纯多项式展开的阶，:<V_…均为 混纯多项式展开的系数，r 1 (4 )，L (A，' )，r?(A，，· ?,)均为所确定的基底类型， A，4，…，4均为标准随机变量；
[0035] S3.3.对于每一个变形的混纯多项式展开，构建基于该混纯多项式展开的MCS仿真 模型，分别进行MCS仿真，计算每一个变形的混沌多项式展开的可靠度数值。
[0036] 作为本发明的进一步改进，所述步骤S3.1的具体步骤包括：
[0037] S3.1.1.从所述关键输入参数中随机采样，生成第二关键输入参数集；
[0038] S3.1.2.确定所述第二关键输入参数集中各关键输入参数的均值的离散水平数；
[0039] S3.1.3.根据预设的误差影响比例，计算第二关键输入参数集中各关键输入参数 的均值对于离散水平数的误差影响值，对所述误差影响值进行正交组合，生成2倍于由所述 步骤S2.6中所确定的η阶混沌多项式展开中包含的未知系数个数的组合数。
[0040] 作为本发明的进一步改进，所述步骤S4包括如下步骤：
[0041] S4.1.根据步骤S2.6中所确定的混沌多项式展开的阶数，将如式（1)所示的混沌多 项式展开表示为如式(4)所示的缩减形式，
[0042]
[0043] 式⑷中，R(l)为可靠度估计值，是混沌多项式展开的系数，!= + 为P阶混沌多项式展开的基底类型，N。为混沌多项式展开所包含系数的总个数；
[0044] S4.2.根据所述可靠度数值，通过加归分析求解方法，计算如式(4)所示的混沌多 项式展开的系数，得到可靠性混沌多项式展开。
[0045] 作为本发明的进一步改进，所述步骤S5的具体步骤包括：
[0046] S5.1.将式(4)按照Sobol'形式展开，并按照变量的个数进行分组，得到分组后的 混沌多项式展开；
[0047] S5.2.计算分组后的混沌多项式展开中各组分的方差，以及式(4)的总方差；
[0048] S5.3.根据Sobol'计算公式，计算混沌多项式展开各阶的可靠性全局灵敏度指标 和单个输入变量的总的可靠性全局灵敏度指标。
[0049] 与现有技术相比，本发明的优点在于：
[0050] 1、本发明采用混沌多项式展开来建立可靠性与输入参数之间的函数关系，仅需要 少量仿真作为基础，便可解析计算Sobol'指标，能够在保证精确性的前提下，大幅减少全局 灵敏度分析的计算量。
[0051] 2、本发明摒弃传统基于参数定义域内某一点来作为可靠性灵敏度大小的评判标 准，本发明考虑不确定性的存在，输入参数的变化区域通常为一个不确定性区间，从整个变 化区间出发去评价可靠性灵敏度，更符合工程实际。
[0052] 3、本发明不仅可以计算某输入参数对可靠性的影响，还可以计算输入参数间的交 叉作用对可靠性的影响，从而可以识别风险因子，为参数识别工作提供支持。
[0053] 4、本发明可以用于校准影响可靠性的输入参数的最佳变化区域，从而可以指导参 数的容差控制等工作。
【附图说明】
[0054]图1为本发明具体实施例流程示意图。
[0055] 图2为本发明具体实施例中组合数示意图。
【具体实施方式】
[0056] 以下结合说明书附图和具体优选的实施例对本发明作进一步描述，但并不因此而 限制本发明的保护范围。
[0057]如图1所示，本实施例一种基于混沌多项式展开的可靠性全局灵敏度分析方法，包 括如下步骤:S1.根据待分析模型的输入参数和输出参数，确定关键输入参数与关键输出参 数;S2.构建用于表征关键输入参数与关键输出参数之间函数关系的满足预设条件的混沌 多项式展开;S3.计算混沌多项式展开的可靠度数值;S4.根据可靠度数值，构建可靠性混沌 多项式展开;S5.基于Sobol'指标，解析计算可靠性全局灵敏度指标。
[0058] 在本实施例中，步骤S1中确定关键输入参数与关键输出参数的步骤为:SI. 1.计算 输入参数和输出参数的概率分布、随机变量特性和随机过程特性;S1.2.判断是否满足任意 一个预设的判定准则，是则判定输入参数为关键输入参数，输出参数为关键输出参数;在本 实施例中，预设的判定准则为:T1.概率分布是否满足预设的分布条件;T2.随机变量特性是 否满足预设的特性阈值;T3.随机过程特性是否满足预设的过程特性阈值。
[0059] 在本实施例中，步骤S2的具体步骤为:S2.1.根据关键输入参数的概率密度函数分 布，结合Wiener-Askey方案，确定关键输入参数的基底类型，基底类型为标准随机变量的函 数;S2.2.将关键输入参数与关键输出参数表征为基底类型的混沌多项式展开，混沌多项式 展开的阶数为n，n多2，n的初始值为2,如式(1)所示：
[0060]
[0061] 式（1)中，y(n)为混沌多项式展开，η为混沌多项式展开的阶，％，土，％，％...4均为 混沌多项式展开的系数，),r 2g4),Γ,成4，.··4")均为所确定的基底类型， A，&,···，&均为标准随机变量;S2.3.从关键输入参数中随机采样，生成第一关键输入参数 集，将第一关键输入参数集中的关键输入参数表示为标准随机变量的转换函数，将n+1阶混 沌多项式展开所确定的基底类型的根作为随机配点输入至转换函数，计算得到第一输入参 数;第一关键输入参数集中元素个数至少为η阶混沌多项式展开中系数个数的2倍;在初始 情况下，混沌多项式展开的阶数为η = 2,则第一关键输入参数集中元素的个数至少为3阶混 沌多项式展开中系数个数的2倍;S2.4.将第一输入参数输入预设的第一蒙特卡洛仿真模 型，计算得到与第一输入参数对应的第一输出参数，生成由第一输入参数与第一输出参数 构成的第一样本数据;S2.5.计算第一样本数据的条件数，并判断条件数是否小于预设的条 件数阈值，是则跳转到步骤S2.6;否则跳转到步骤S2.3; S2.6.根据第一样本数据，通过改进 的概率配点法（EPCM，Extension of probabilistic collocation method)和回归分析法 计算如式（1)所示的n阶和n+1阶混沌多项式展开的系数，计算n阶混沌多项式展开的值、n+1 阶混沌多项式展开的值，并判断η阶混沌多项式展开的值、n+1阶混沌多项式展开的值和第 一输出参数值之间的误差是否小于预设的误差阈值，是则确定混沌多项式展开的阶为n，否 贝1J，将混沌多项式展开的阶数加1，跳转至步骤S2.2。
[0062] 在本实施例中，步骤S3的具体步骤为:S3.1.对关键输入参数采用内外表直积法进 行水平组合，得到2倍于由步骤S2.6中所确定的η阶混沌多项式展开中包含的未知系数个数 的组合数;S3.2.计算组合数中各关键输入参数的水平偏离中心值的大小△ h，对于每一个 组合数，将如式(1)所示的混沌多项式展开变形为如式(2)所示的形式；
[0063]
[0064] 式(2)中，y(n)为混沌多项式展开，η为混沌多项式展开的阶A 均为混沌 多项式展开的系数，r.i 鐵 +Δ4 +Δ久)，(式，(C: +十从^^ 均为变形后的基底类型，4,我,· · ·4"均为标准随机变量，，機2，…，Δ/\均为各关键输入参 数的水平偏离中心值的大小;将式(2)展开合并，变形为如式(3)所示形式：
[0065]
[0066] 式（3)中，y(n)为混沌多项式展开，η为混沌多项式展开的阶，心，气，气*2，~ 2: λ均为 混沌多项式展开的系数，),Γ2(式,4)均为所确定的基底类型， 均为标准随机变量;S3.3.对于每一个变形的混沌多项式展开，构建基于该混沌 多项式展开的MCS仿真模型，分别进行MCS仿真，计算每一个变形的混沌多项式展开的可靠 度数值。
[0067]在本实施例中，步骤S3.1的具体步骤为:S3.1.1.从关键输入参数中随机采样，生 成第二关键输入参数集;S3.1.2.确定第二关键输入参数集中各关键输入参数的均值的离 散水平数;S3.1.3.根据预设的误差影响比例，计算第二关键输入参数集中各关键输入参数 的均值对于离散水平数的误差影响值，对误差影响值进行正交组合，生成2倍于由步骤S2.6 中所确定的η阶混沌多项式展开中包含的未知系数个数的组合数。
[0068]在本实施例中，假定所研究的模型具有三个设计变量，三个设计变量分别所对应 的标准随机变量为:L，ξ3，每个标准随机变量具有关于零点对称的两个离散水平数，分 另IJ为土a，土b，土c，选定的正交表类型为L 4(23)，预设的误差影响比例为±10%，则对离散水 平数的误差影响值采用内外表直积法，可得到16组组合数，如图2所示。
[0069] 在本实施例中，步骤S4的具体步骤为:S4.1.根据步骤S2.6中所确定的混沌多项式 展开的阶数，将如式(1)所示的混沌多项式展开表示为如式(4)所示的缩减形式，
[0070]
[0071 ]式(4)中，R")为可靠度估计值，尤是混纯多项式展开的系数為(&Π 二< ◎ =「,4炙 为P阶混沌多项式展开的基底类型，N。为混沌多项式展开所包含系数的总个数;S4.2.根据 可靠度数值，通过加归分析求解方法，计算如式(4)所示的混沌多项式展开的系数，得到可 靠性混沌多项式展开(RPCE)，即RPCE的最终形式。
[0072] 在本实施例中，步骤S5的具体步骤为:S5.1.将式(4)按照Sobol'形式展开，并按照 变量的个数进行分组，得到分组后的混沌多项式展开;S5.2.计算分组后的混沌多项式展开 中各组分的方差，以及式(4)的总方差;S5.3.根据Sobol'计算公式，计算混沌多项式展开各 阶的可靠性全局灵敏度指标和单个输入变量的总的可靠性全局灵敏度指标。
[0073] 在本实施例中，步骤S5.1中按照变量的个数进行分组的具体步骤为:对被加项不 展开，将单独的各个变量的各次项分为一组，将具有两个变量的被加项的各次项分为一组， 依此类推，得到分组后的混沌多项式展开，其数学表达式为： 「00741
[0075] 式(5)中，(6=(?，…，βη)是一个整数序列，且满足队彡0，
Φ为PCE的阶数， a的定义为
，?为混沌多项式展开的系数，％! (L，···")为混沌多项式展开所选择的基底类型；
[0076] 计算分组后的混沌多项式展开中各组分的方差，如式(6)所示，
[0077]
[0078] 式(6)中，￡\,..4为输入参数组合{h，…，is}的方差，β= (ft，…，βη)是一个整数序列，且 满足私多0，
Φ为PCE的阶数，的定义为：
&为混沌多项式展开的系数，为混沌多项式展开所选择基底类型的数学期望。
[0079] 混沌多项式展开的总方差如式(7)所示， 「00801
[0081]式(7)中，Dr为混沌多项式展开的方差，Var[RU)]为可靠度估计值RU)的方差，；；: 为混沌多项式展开的系数，幻]为混沌多项式展开所选择基底类型的数学期望;Nc为 混沌多项式展开所包含系数的总个数。
[0082]根据如式(8)所示各阶Sobol'计算公式，
[0083] Sh^^Dh..jD (8)
[0084] 式⑶中，5；.1为输入参数组合{^，…山}的全局灵敏度指标，A,为输入参数组 合，…，is}的方差，D为混沌多项式展开的总方差。
[0085] 将式(6)和式(7)代入式(8)，即可计算混沌多项式展开各阶的可靠性全局灵敏度 指标，如式(9)所示，
[0086]
[0087] 式(8)中，表示输入参数组合{ii,…，is}的全局灵敏度，β= (βι,…，βη)是一 个整数序列，且满足β :彡〇，
，P为P C E的阶数，气...4的定义为：
为混沌多项式展开的系数，五为混沌多项 式展开所选择基底类型的数学期望。
[0088] 单个输入变量的总的可靠性全局灵敏度指标，如式(10)所示，
[0089]
[0090] 在式（10)中，Sid..#为单个输入变量的总的可靠性全局灵敏度，上，…，jt为一个给 定的整数序列的定义为 ={(/ρ·.·Λ),α，···,乂）ch,·.·,/,)}，
[0091] 在本实施例中，采用混沌多项式展开来建立可靠性与输入参数之间的函数关系， 仅需要少量仿真作为基础，便可解析计算Sobol'指标，能够在保证精确性的前提下，大幅减 少全局灵敏度分析的计算量;摒弃传统基于参数定义域内某一点来作为可靠性灵敏度大小 的评判标准，本发明考虑不确定性的存在，输入参数的变化区域通常为一个不确定性区间， 从整个变化区间出发去评价可靠性灵敏度，更符合工程实际。同时，本实施例不仅可以计算 某输入参数对可靠性的影响，还可以计算输入参数间的交叉作用对可靠性的影响，从而可 以识别风险因子，为参数识别工作提供支持;可以用于校准影响可靠性的输入参数的最佳 变化区域，从而可以指导参数的容差控制等工作。
[0092] 上述只是本发明的较佳实施例，并非对本发明作任何形式上的限制。虽然本发明 已以较佳实施例揭露如上，然而并非用以限定本发明。因此，凡是未脱离本发明技术方案的 内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在 本发明技术方案保护的范围内。
【主权项】
1. 一种基于混浊多项式展开的可靠性全局灵敏度分析方法，其特征在于，包括如下步 骤：51. 根据待分析模型的输入参数和输出参数，确定关键输入参数与关键输出参数；52. 构建用于表征所述关键输入参数与关键输出参数之间函数关系的满足预设条件的 混浊多项式展开；53. 计算所述混浊多项式展开的可靠度数值；54. 根据所述可靠度数值，构建可靠性混浊多项式展开；55. 基于Sobol'指标，解析计算可靠性全局灵敏度指标。2. 根据权利要求1所述的基于混浊多项式展开的可靠性全局灵敏度分析方法，其特征 在于，所述步骤S1中确定所述关键输入参数与关键输出参数的步骤包括： 51.1. 计算所述输入参数和输出参数的概率分布、随机变量特性和随机过程特性； 51.2 .判断是否满足任意一个预设的判定准则，是则判定所述输入参数为关键输入参 数，所述输出参数为关键输出参数； 所述预设的判定准则包括： T1.所述概率分布是否满足预设的分布条件； T2.所述随机变量特性是否满足预设的特性阔值； T3.所述随机过程特性是否满足预设的过程特性阔值。3. 根据权利要求2所述的基于混浊多项式展开的可靠性全局灵敏度分析方法，其特征 在于:所述步骤S2的具体步骤包括： 52.1. 根据所述关键输入参数的概率密度函数分布，结合Wiener-Askey方案，确定关键 输入参数的基底类型，所述基底类型为标准随机变量的函数； 52.2 .将所述关键输入参数与关键输出参数表征为所述基底类型的混浊多项式展开， 所述混浊多项式展开的阶数为n，n>2，n的初始值为2,如式(1)所示：式(1)中，y(n)为混浊多项式展开，η为混浊多项式展开的阶，"《，和％，V·诚为混浊多 项式展开的系数，Γ,(妾,)，Γ2抬，名ΧΓ,,指,，4，''·，每}均为所确定的基底类型，妾,，4,…為均 为标准随机变量； 52.3. 从所述关键输入参数中随机采样，生成第一关键输入参数集，将所述第一关键输 入参数集中的关键输入参数表示为所述标准随机变量的转换函数，将η+1阶混浊多项式展 开所确定的基底类型的根作为随机配点输入至所述转换函数，计算得到第一输入参数;所 述第一关键输入参数集中元素个数至少为所述η阶混浊多项式展开中系数个数的2倍； 52.4. 将所述第一输入参数输入预设的第一蒙特卡洛仿真模型，计算得到与所述第一 输入参数对应的第一输出参数，生成由第一输入参数与第一输出参数构成的第一样本数 据； 52.5. 计算所述第一样本数据的条件数，并判断所述条件数是否小于预设的条件数阔 值，是则跳转到步骤S2.6;否则跳转到步骤S2.3; 52.6. 根据所述第一样本数据，通过改进的概率配点法和回归分析法计算如式（1)所示 的η阶和n+1阶混浊多项式展开的系数，计算η阶混浊多项式展开的值、n+1阶混浊多项式展 开的值，并判断所述η阶混浊多项式展开的值、n+1阶混浊多项式展开的值和所述第一输出 参数值之间的误差是否小于预设的误差阔值，是则确定所述混浊多项式展开的阶为n，否 贝1J，将所述混浊多项式展开的阶数加1，跳转至步骤S2.2。4. 根据权利要求3所述的基于混浊多项式展开的可靠性全局灵敏度分析方法，其特征 在于，所述步骤S3包括如下步骤： 53.1. 对所述关键输入参数采用内外表直积法进行水平组合，得到2倍于由所述步骤 S2.6中所确定的η阶混浊多项式展开中包含的未知系数个数的组合数； 53.2. 计算所述组合数中各关键输入参数的水平偏离中屯、值的大小Δ h，对于每一个组 合数，将所述如式(1)所示的混浊多项式展开变形为如式(2)所示的形式；式(2)中，y(n)为混浊多项式展开，η为混浊多项式展开的阶，"nA,，％，％.、,均为混浊多项 式展开的系数，均为变形后的基底类型，《,，妾;，…，焉均为标准随机变量，ΔΑ,,，ΔΑ,:，…，Δ^。均为各关键输入参 数的水平偏离中屯、值的大小； 将所述式(2)展开合并，变形为如式(3)所示形式：式(3)中，y(n)为混浊多项式展开，η为混浊多项式展开的阶，0。，"。，％，％-蝴为混浊多 项式展开的系数，。(每)，。(写，妾;)，Γ,,(、;，、;，···，妾均为所确定的基底类型，《,,，《:，···，《。均 为标准随机变量； 53.3. 对于每一个变形的混浊多项式展开，构建基于该混浊多项式展开的MCS仿真模 型，分别进行MCS仿真，计算每一个变形的混浊多项式展开的可靠度数值。5. 根据权利要求4所述的基于混浊多项式展开的可靠性全局灵敏度分析方法，其特征 在于，所述步骤S3.1的具体步骤包括： 53.1.1. 从所述关键输入参数中随机采样，生成第二关键输入参数集； 53.1.2. 确定所述第二关键输入参数集中各关键输入参数的均值的离散水平数； 53.1.3. 根据预设的误差影响比例，计算第二关键输入参数集中各关键输入参数的均 值对于离散水平数的误差影响值，对所述误差影响值进行正交组合，生成2倍于由所述步骤 S2.6中所确定的η阶混浊多项式展开中包含的未知系数个数的组合数。6. 根据权利要求5所述的基于混浊多项式展开的可靠性全局灵敏度分析方法，其特征 在于，所述步骤S4包括如下步骤： S4.1.根据步骤S2.6中所确定的混浊多项式展开的阶数，将如式（1)所示的混浊多项式 展开表示为如式(4)所示的缩减形式，(4； 式(4)中，R似为可靠度估计值，/,是混浊多项式展开的系！为P阶混浊多项式展开的基底类型，N。为混浊多项式展开所包含系数的总个数； S4.2.根据所述可靠度数值，通过加归分析求解方法，计算如式(4)所示的混浊多项式 展开的系数，得到可靠性混浊多项式展开。7.根据权利要求6所述的基于混浊多项式展开的可靠性全局灵敏度分析方法，其特征 在于，所述步骤S5的具体步骤包括： 55.1. 将式(4)按照Sobol'形式展开，并按照变量的个数进行分组，得到分组后的混浊 多项式展开； 55.2. 计算分组后的混浊多项式展开中各组分的方差，W及式(4)的总方差； 55.3. 根据Sobol'计算公式，计算混浊多项式展开各阶的可靠性全局灵敏度指标和单 个输入变量的总的可靠性全局灵敏度指标。
【文档编号】G06F17/50GK106096138SQ201610414171
【公开日】2016年11月9日
【申请日】2016年6月14日 公开号201610414171.8, CN 106096138 A, CN 106096138A, CN 201610414171, CN-A-106096138, CN106096138 A, CN106096138A, CN201610414171, CN201610414171.8
【发明人】杜绍华, 周桂法, 汪旭, 陈旭鸿, 匡芬, 潘宇雄, 袁莹莹
【申请人】中车株洲电力机车研究所有限公司