结构机构失效概率灵敏度分解方法、计算方法及应用与流程

本公开涉及可靠性及稳健性分析技术领域，尤其涉及一种结构机构失效概率灵敏度分解方法、计算方法及应用。

背景技术

失效概率是在规定条件下机构不能完成既定任务事件的概率，是衡量机构可靠性的重要指标。灵敏度分析是衡量变量不确定性对感兴趣输出量(如输出方差、分布、失效概率等)的贡献程度，为改进输出提供依据。通过失效概率灵敏度的分析可以得到决定结构失效的输入变量分布参数的相对重要程度，从而为结构可靠性分析、预测和优化提供指导。

然而，工程结构中输入变量之间通常存在着相关性，影响结构失效概率的各个输入变量不服从独立分布。输入变量相关性与独立性共同影响失效概率灵敏度。通过现有技术计算得到的失效概率灵敏度为输入变量相关性与独立性共同作用的灵敏度，其并不能用来衡量输入变量相关性与独立性在失效概率灵敏度分析中作用大小。

需要说明的是，在上述背景技术部分公开的信息仅用于加强对本公开的背景的理解，因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。

技术实现要素：

本公开的目的在于提供一种结构机构失效概率灵敏度分解方法、计算方法及应用，该方法通过引入copula函数来描述输入变量的联合概率分布，将失效概率灵敏度分解为独立灵敏度和相关灵敏度，进而分离出变量相关性与独立性分别对失效概率灵敏度计算的影响，研究得到的技术方法为设计人员提供了衡量输入变量相关性与独立性在失效概率灵敏度分析中作用大小的理论依据。

根据本发明的一个方面，提供一种基于copula函数的结构机构失效概率灵敏度分解方法，包括：

将多维随机变量联合分布函数表示为copula函数形式，从而建立所述多维随机变量联合分布函数与边缘概率密度函数的关系式：

其中，fx(x)为所述联合概率密度函数，c(u1,...,un)为所述copula函数，为所述边缘概率密度函数；

基于所述关系式将所述结构机构的失效概率灵敏度分解为独立灵敏度与相关灵敏度之和。

在本发明的一种示例性实施例中，所述基于所述关系式将所述结构机构的失效概率灵敏度分解为独立灵敏度与相关灵敏度之和包括：

将所述联合概率密度函数对数的偏导数表示为：

其中，为边缘分布的核函数为copula密度函数的核函数

计算所述结构机构的失效概率：

其中，为独立灵敏度为相关灵敏度if(x)为失效概率指示函数。

在本发明的一种示例性实施例中，所述copula函数包括为双变量copula函数或者多变量copula函数。

在本发明的一种示例性实施例中，所述双变量copula函数为u1u2(1+θ(1-u1)(1-u2))、φg[φ^-1(u1),φ^-1(u2)；θ]、u1u2(1-θ(1-u1)(1-u2))^-1中的一种。

在本发明的一种示例性实施例中，所述多变量copula函数为fgm函数族和clayton函数族中的一种。

在本发明的一种示例性实施例中，所述边缘概率密度函数为

中的一种。

根据本发明的一方面，提供一种结构机构失效概率灵敏度计算方法，应用上述的基于copula函数的结构机构失效概率灵敏度分解方法，包括：

根据所述结构机构的边缘概率密度函数随机抽取多个输入变量样本向量xj＝(x1,j,...,xn,j)(j＝1,...,n)，基于所述变量样本向量xj＝(x1,j,...,xn,j)(j＝1,...,n)根据所述结构机构的功能函数计算所述结构机构的失效概率指示函数样本if(xj)。

计算对应于所述输入变量样本向量的累积分布函数值向量uj＝(u1,j,...,un,j)，并根据copula密度函数计算密度函数样本值c(uj)；

根据边缘输入变量样本计算核函数样本并根据所述累积分布函数样本向量uj得到的copula核函数样本cki(uj)。

根据公式计算独立灵敏度

根据公式计算相关灵敏度

根据公式计算所述结构机构失效概率灵敏度

根据本发明的一方面，提供一种y型节点管失效概率灵敏度计算方法，其特征在于，应用上述的结构机构失效概率灵敏度计算方法。

本发明提出一种结构机构失效概率灵敏度分解方法，该方法通过引入copula函数来描述输入变量的联合概率分布，将失效概率灵敏度分解为独立灵敏度和相关灵敏度，进而分离出变量相关性与独立性分别对失效概率灵敏度计算的影响。一方面，该方法得到的技术方法为设计人员提供了衡量输入变量相关性与独立性在失效概率灵敏度分析中作用大小的理论依据；另一方面，该方法计算过程简单，可实施性强。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本发明。

附图说明

此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本公开基于copula函数的结构机构失效概率灵敏度分解方法的流程图；

图2为本公开结构机构失效概率灵敏度计算方法一种示例性实施例的流程图；

图3为本公开y型节点管失效概率灵敏度计算方法一种示例性实施例中y型节点管的受力图。

具体实施方式

现在将参考附图更全面地描述示例实施方式。然而，示例实施方式能够以多种形式实施，且不应被理解为限于在此阐述的范例；相反，提供这些实施方式使得本公开将更加全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中，提供许多具体细节从而给出对本公开的实施方式的充分理解。然而，本领域技术人员将意识到，可以实践本公开的技术方案而省略所述特定细节中的一个或更多，或者可以采用其它的方法、组元、装置、步骤等。在其它情况下，不详细示出或描述公知技术方案以避免喧宾夺主而使得本公开的各方面变得模糊。

此外，附图仅为本公开的示意性图解，并非一定是按比例绘制。图中相同的附图标记表示相同或类似的部分，因而将省略对它们的重复描述。附图中所示的一些方框图是功能实体，不一定必须与物理或逻辑上独立的实体相对应。可以采用软件形式来实现这些功能实体，或在一个或多个硬件模块或集成电路中实现这些功能实体，或在不同网络和/或处理器装置和/或微控制器装置中实现这些功能实体。

为方便理解本实施例中失效概率灵敏度的分解过程，本实施例性实施例首先对失效概率及失效概率灵敏度进行定义。

考虑结构机构的功能函数为y＝g(x)，其中x＝(x1,x2,...,xn)为相关输入变量，其不确定性可以通过联合概率密度函数fx(x)来描述。结构失效概率pf可以表示为积分形式：

其中，设存在n个随机输入变量影响结构机构可靠性，并将影响该结构机构可靠性的n维随机输入变量表述为x＝(x1,x2,...,xn)∈rⁿ，rⁿ为n维变量空间，n为大于等于1的整数。if(x)为失效概率指示函数(当x∈f时if(x)＝1，当时if(x)＝0)，e[·]为求均值函数。

为衡量分布参数μi的微小扰动对失效概率的影响，定义基于偏导数的失效概率灵敏度：

式中μi只作为输入变量xi分布参数的符号表示，实际分析过程中需根据分布类型确定分布参数。式(2)将失效概率对分布参数的偏导数转化为联合概率密度函数对数对分布参数偏导数的积分，避免了直接利用差分法对失效概率灵敏度进行计算所带来的计算精度和计算效率问题。

基于此，本示例性实施例首先提供一种基于copula函数的结构机构失效概率灵敏度分解方法，如图1所示，为本公开基于copula函数的结构机构失效概率灵敏度分解方法的流程图。包括：

步骤s1：将多维随机变量联合分布函数表示为copula函数形式，从而建立所述多维随机变量联合分布函数与边缘概率密度函数的关系式：

其中，fx(x)为所述联合概率密度函数，c(u1,...,un)为所述copula函数，为所述边缘概率密度函数；

步骤s2：基于所述关系式将所述结构机构的失效概率灵敏度分解为独立灵敏度与相关灵敏度之和。

本发明提出一种结构机构失效概率灵敏度分解方法，该方法通过引入copula函数来描述输入变量的联合概率分布，将失效概率灵敏度分解为独立灵敏度和相关灵敏度，进而分离出变量相关性与独立性分别对失效概率灵敏度计算的影响。一方面，该方法得到的技术方法为设计人员提供了衡量输入变量相关性与独立性在失效概率灵敏度分析中作用大小的理论依据；另一方面，该方法计算过程简单，可实施性强。

本示例性实施例中，将多维随机变量联合分布函数表示为copula函数形式，从而建立所述多维随机变量联合分布函数与边缘概率密度函数的关系式的具体过程如下：

根据sklar理论，多维随机变量的联合分布函数可以表示为copula函数的形式：

其中，为各随机变量的边缘分布函数值。

根据式(3)，多维随机变量的联合概率密度函数可以表示为

其中，可以通过积累统计数据得到输入变量的边缘概率密度函数

需要说明的是，copula函数的相关性可以通过kendall的秩相关系数τ来衡量，它不仅能够描述变量之间的线性相关性，也可以对变量之间的非线性相关性进行描述。因此，该方法可以适用于输入变量为非线性关系的失效概率灵敏度分析。

本示例性实施例中，所述基于所述关系式将所述结构机构的失效概率灵敏度分解为独立灵敏度与相关灵敏度之和可以包括：

根据式(3)将所述联合概率密度函数对数的偏导数表示为：

其中，为边缘分布的核函数为copula密度函数的核函数

将式(5)带入式(2)计算所述结构机构的失效概率：

其中，为独立灵敏度为相关灵敏度if(x)为失效概率指示函数。

需要说明的是，由式(5)可知，联合概率密度函数对数的偏导数可以分解为两部分，第一部分为边缘概率密度函数对数的偏导数，本示例性实施例将其定义为边缘分布的核函数第二部分为copula密度函数对数的偏导数，本示例性实施例类似地定义为copula密度函数的核函数

由式(6)可知，针对变量相关的情况，分布参数μi对失效概率的灵敏度可以分解为两部分，第一部分e(if(x)·ki)表示分布参数μi的微小扰动通过边缘概率密度函数(即变量独立部分)对失效概率的影响，称为独立灵敏度第二部分e(if(x)·cki)表示分布参数μi的微小扰动通过copula密度函数(即变量相关部分)对失效概率的影响，称为相关灵敏度如果则表示分布参数μi对失效概率的影响主要是通过输入变量xi的边缘统计特性进行传递；如果则表示分布参数μi对失效概率的影响主要是通过输入变量xi与其他变量之间的相关性统计特征进行传递。因此，式(6)可以很清晰地反映出分布参数μi对失效概率的灵敏度与变量之间相关性的关系。

上述式(6)中包含有边缘分布的核函数以下介绍针几种常用的边缘概率密度函数分布形式以及其对应的边缘分布的核函数如表1所示：

表1常用的几种分布类型及其核函数

由表1可以看到，在给定分布参数μi的条件下，对于常用分布类型的核函数其数值只和边缘输入变量xi的值有关，所以可以将核函数表示为输入变量的函数，记为因此，在利用抽样方法进行失效概率计算时，只需利用计算过程中得到的输入变量xi的样本值，便可以根据表1中的核函数解析表达式计算得到核函数进而结合已有的功能函数计算结果得到独立灵敏度该过程将在下面结构机构失效概率灵敏度计算方法中得以体现。需要说明的是，在其他实施例中，边缘概率密度函数以及其对应的边缘分布的核函数还有更多的表现形式，这些都属于本公开的保护范围。

上述式(6)中还包含有copula密度函数的核函数以下介绍针几种常用的copula函数以及其对应的copula密度函数的核函数

本示例性实施例中，copula函数可以为双变量copula函数。如果输入变量xi之间存在两两相关的情况，则可以利用双变量copula函数描述两个变量之间的相关性，以下介绍几种常用的copula函数以及其对应的copula核函数，如表2所示。

表2五种常见的双变量copula函数及对应的核函数

需要说明的是，在其他实施例中，双变量copula函数及对应的核函数还有更多的表现形式，这些都属于本公开的保护范围。

上述双变量copula函数只能用于描述两个变量之间的相关性，实际中还存在着多个变量相关的情况，本示例性实施例中，copula函数还可以为多变量copula函数，多变量copula函数可以用于描述多个输入变量的相关性。

本示例性实施例中，多变量copula函数可以为fgm函数族和clayton函数族，以下分别对fgm函数族和clayton函数族进行介绍。

fgm型copula函数族的一般表达式可以表示为：

其中，θ为维度空间{1,2,...,d}的子集，并且至少包含两个元素。θs为任意子集s∈θ的连接参数，并且当s＝{i1,i2,...,ik}，

不难证明，任意fgm型copula函数是绝对连续的，并且可以给出其密度函数cd(u1,u2,...,ud)如下:

根据式(5)中核函数的定义，fgm型copula函数族的核函数可以表示为

clayton型copula函数族是archimedean型copula的一种，其一般表达式可以表示为：

其中，参数且θ≠0。当θ＝0时，表示变量相互独立。

相应地，clayton型copula的密度函数可以表示为：

根据copula核函数的定义，fgm型copula函数族的核函数可以表示为

需要说明的是，在其他示例性实施例中，双变量copula函数还有更多的表现形式，这些都属于本公开的保护范围。

通过上述内容对五种双变量copula函数以及两类典型的多变量copula函数族的核函数进行了推导。可以看到，copula核函数表达式中包含着边缘累积分布函数对分布参数的偏导数根据累积分布函数的积分表达式，可以作如下推导：

其中，为边缘分布的核函数。因此，对于偏导数可以直接利用已有样本值和得到的核函数进行计算，进而将其代入选定的copula核函数计算式中得到其计算结果，该过程将在下面结构机构失效概率灵敏度计算方法中得以体现。

本示例性实施例还提供一种结构机构失效概率灵敏度计算方法，应用上述的基于copula函数的结构机构失效概率灵敏度分解方法，如图2所示，为本公开结构机构失效概率灵敏度计算方法一种示例性实施例的流程图。该方法包括：

步骤s1：根据所述结构机构的边缘概率密度函数随机抽取多个输入变量样本向量xj＝(x1,j,...,xn,j)(j＝1,...,n)，基于所述变量样本向量xj＝(x1,j,...,xn,j)(j＝1,...,n)根据所述结构机构的功能函数计算所述结构机构的失效概率指示函数样本if(xj)。

步骤s2：计算对应于所述输入变量样本向量的累积分布函数值向量uj＝(u1,j,...,un,j)，并根据copula密度函数计算密度函数样本值c(uj)；

步骤s3：根据边缘输入变量样本计算核函数样本并根据所述累积分布函数样本向量uj得到的copula核函数样本cki(uj)。

步骤s4：根据公式：

计算独立灵敏度

步骤s5：根据公式：

计算相关灵敏度

步骤s6：根据公式计算所述结构机构失效概率灵敏度

需要说明的是，步骤s1中，xj表示第j次抽取的样本向量，(x1,j,...,xn,j)分别表示第j次抽取的n个输入变量样本。步骤s2中，u1,j分别表示第j次抽取的输入变量样本对应的累计分布函数值。uj表示n个累计分布函数值组合而成的积分布函数值向量。步骤s3中，表示(x1,j,...,xn,j)中第j次抽取的第i个输入变量样本。步骤s4中，公式(7)可以通过公式推导得出。步骤s5中，公式(8)可以通过公式推导得出。

本示例性实施例还提供一种y型节点管失效概率灵敏度计算方法，其特征在于，应用上述的结构机构失效概率灵敏度计算方法。

如图3所示，为本公开y型节点管失效概率灵敏度计算方法一种示例性实施例中y型节点管的受力图。该y型节点管剖面a-a的载荷包括轴向力f、面内弯矩mi、截面扭矩m0。d为主管外径，d为侧管外直径，b为侧管内半径，θ为主管与侧管夹角。

根据该管道的受力情况可以建立极限状态方程如下式：

g＝0.8-2×10^-5f-(2×10^-4|mi|)^1.2-(2×10^-5|m0|)^2.1

其中，f、mi和m0为正态分布随机变量，其分布参数见表3。考虑轴向力f、面内弯矩mi、截面扭矩m0具有相关性，并可以通过clayton型copula函数建立其联合分布。考虑copula函数的kendall秩相关系数τ＝0.5时，y型节点管的失效概率及其灵敏度分析结果列于表4中。

表3y型节点管中随机变量的分布参数

表4y型节点管的失效概率灵敏度计算结果

表4中给出了输入变量独立和相关时的失效概率计算结果，可以看到：对于y型节点管，输入变量之间存在相关性时将导致结构的失效概率的增加。对于失效概率灵敏度，面内弯矩mi的分布参数(μ2，σ2)对失效概率的影响最大，并且其影响主要是由独立灵敏度引起的；而对于截面扭矩m0的分布参数(μ3，σ3)，尽管其失效概率灵敏度最小，但其独立灵敏度和相关灵敏度都比较大。因此在y型节点管结构设计中，针对截面扭矩m0的分布参数对失效概率的影响，可以单独考虑该输入变量的独立统计信息；而针对截面扭矩m0的分布参数对失效概率的影响，则必须综合考虑截面扭矩m0的统计信息以及与其他变量的相关性统计信息。

本领域技术人员在考虑说明书及实践这里发明的发明后，将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未发明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本发明的真正范围和精神由所附的权利要求指出。

上述所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中，如有可能，各实施例中所讨论的特征是可互换的。在上面的描述中，提供许多具体细节从而给出对本发明的实施方式的充分理解。然而，本领域技术人员将意识到，可以实践本发明的技术方案而没有特定细节中的一个或更多，或者可以采用其它的方法、组件、材料等。在其它情况下，不详细示出或描述公知结构、材料或者操作以避免模糊本发明的各方面。