混合不确定性下的结构可靠性分析方法

2018-10-18肖宁聪周成宁张成林刘志亮

电子科技大学学报 2018年5期

肖宁聪，周成宁，张成林，刘志亮

(1.湖南大学汽车车身先进设计制造国家重点实验室长沙 410082；2.电子科技大学机械电子工程学院成都 611731；3.上海航天设备制造总厂上海闵行区 200245)

在工程中，通常会出现各种各样的不确定性问题，如外部载荷、材料的属性、加工尺寸、近似模型误差等，不确定性是影响结构可靠性最为关键的因素之一[1-2]。不确定性通常可以分为两大类：随机不确定性和认知不确定性[3]。随机不确定性是由于事物的固有波动性引起的，在数学中通常可以用随机变量进行建模。认知不确定性由于数据不足、信息量少等引起的，可用区间理论、模糊数学、可能性理论等进行描述。

机械零部件通常有多种失效模式(如疲劳、磨损、腐蚀等)，由于各失效模式通常有共享变量及相同的因素影响，因而有较强的非线性相关性。多失效模式可靠性方法一直是待解决的难点问题，相关学者相继提出了一系列方法，如二阶窄方法[4]、互补交集模型[5]、鞍点近似法[6]等。然而，大多数现有多失效模式可靠性方法只能刻画线性相关性，而对非线性问题不能有效解决。不仅如此，现有方法大多集中在随机不确定性下，混合不确定性下的多模式可靠性方法研究鲜有报道。近年来，Copula函数被用于对变量及失效模式间复杂非线性相关问题进行建模[7]。文献[8]提出了基于鞍点近似及Copula函数的结构系统可靠性分析方法；文献[9]提出了基于高斯Copula函数的可靠性设计优化方法；文献[10]用Copula函数对相关性进行建模，详细分析了假定失效模式间相互独立在可靠性分析中存在的误差。虽然基于Copula函数多失效模式可靠性分析有了相应的研究，但是绝大部分建立在随机不确定下，混合不确定下基于Copula函数的结构系统可靠性分析方法鲜有报道。针对工程中普遍存在的混合不确定性问题，分别用随机变量和区间变量对随机和认知不确定性进行建模，为了提高计算效率，用先进的仿真方法对随机变量进行采样从而确定性能函数的极值响应。为了最大限度地减少人为主观信息等带来的误差，利用最大熵方法对极值响应进行分布近似，通过Copula函数最终近似确定系统失效概率的最大和最小值。最后给出一个工程算例验证本文方法的精度和有效性。

1 确定混合不确定性下性能函数响应的极值

式中，Zi为第i个性能函数的输出响应，

设随机变量Xi的分布函数为则Xi的Nr个样本可表示为：

基于随机矢量X的Nr个样本，则第i个性能函数在区间变量上的极值响应(最大和最小值)分别表示为：

现有很多优化方法可用来求解式(3)和式(4)的优化模型，如遗传算法、序列规划算法等。由式(3)和式(4)可得系统Nr个最小和最大值响应分别为：

2 基于最大熵理论近似极值响应的概率密度函数

最大熵原理为在现有的数据和样本情况下，确定出熵最大的一种概率分布作为现有数据和样本量的概率分布函数。由于最大熵在近似随机变量的概率密度函数时不需要任何额外的人为主观信息及假设，因此是一种精度较高的方法[11]。利用最大熵原理近似随机变量Z的分布可表示为[12]：

满足以下约束：

式中，p(Z)为随机变量Z的概率密度函数；Sz为积分域；为第i阶原点矩。

根据拉格朗日乘子法，p(Z)可表示为：

式中，λi为待定的拉格朗日乘子，可用牛顿迭代优化算法求得。

研究表明，在一般运用中仅用前四阶矩即可满足精度要求，由式(5)和式(6)可知，第k个失效模式最小和最大值响应的前四阶矩可表示为：

由式(10)～式(12)，则第k个失效模式最小和最大值响应的最大熵密度函数为：

由式(13)、式(14)可知第k个失效模式的失效概率最小和最大值分别可计算为：

式(15)、式(16)中的积分没有解析解，可借助辛普森或梯形积分算法求解。

3 利用Copula函数确定系统失效概率的最大和最小值

3.1 Copula函数

式中，θ为Copula函数的参数。若全为连续的边缘分布函数，则Copula函数唯一确定。为了方便和简单起见，仅考虑二元Copula函数的情况，多元Copula函数情况可在此基础上类推，常见二元Copula函数如下[14]：

Clayton Copula函数：

Gumbel Copula函数：

Gauss Copula函数：

3.2 混合不确定性下Copula函数可靠性建模

由式(17)可知，具有两个失效模式串联系统的失效概率可表示为：

需指出的是，当系统存在区间变量时候，系统失效概率是区间而非精确值，此时系统失效概率的最大和最小值可分别近似表示为下述优化模型：

式中，θU和θL分别为Copula函数参数的上下界，为第k个失效模式的失效概率最小和最大值。

为了确定Copula函数参数的上下界，假定Copula函数的结构不变，首先在区间矢量Y的取值范围内随机产生由最大熵方法可得出两个失效模式分别在下的最大熵密度函数、分布函数和失效概率对所得的最大熵分布函数进行采样，并根据分步估计和最大似然估计法，可估计出下的Copula函数参数样本的对数似然函数可表示为：

式中，F1和F2分别表示失效模式1、2的最大熵分布函数；下的最大熵概率密度函数；Nh为获取的样本量。

在满足一定精度要求的条件下，Uθ和Lθ可分别近似表示为：

4 算例分析

某空心压杆，如图1所示，具有稳定性及强度两种失效模式，其性能函数分别为：

式中，F、E、S、d、t、l分别表示轴向载荷、材料弹性模量、材料屈服极限、截面中径、壁厚和杆长，单位分别为kN、GPa、MPa、mm、mm、mm。为了说明本文方法的有效性，壁厚假定为区间变量，其上下界分别为2.42 mm和2.38 mm。变量的相关分布信息如表1所示。

图1 空心杆

表1 变量的分布信息

为了求解各失效模式失效概率的最大和最小值，首先用拉丁方采样方法对各随机变量进行采样，样本量为5 000。由式(3)、式(4)可得样本在区间变量上的最大和最小值，并分别进行最大熵密度函数估计，如g1、g2响应最大值的最大熵密度函数分别如图2和图3所示。由式(15)、式(16)可得失效模式g1、失效概率的最大和最小值。需特别指出的是，变量t为区间变量，只代表其可在区间范围内任意取值，而并不知道其任何其他信息，因此不能等同于其在所取值范围内均匀分布。另外，选择Copula函数对失效模式相关性进行建模，由式(22)、式(23)可得参数θ的最大和最小值估计分别为1.92、1.55。

图2 g1响应最大值的密度函数

为了说明本文方法的精度和有效性，采用蒙特卡罗方法进行对比验证，样本量为610。当系统中存在区间变量时，为了确定系统失效概率最大和最小值，可把区间变量平均分成若干小区间，分别计算区间变量在小区间端点时系统的失效概率并确定失效概率的极值，所得系统失效概率最大和最小值如表2所示。从表2可知，用本文方法所得的结果与蒙特卡罗所得的结果较为接近。然而，蒙特卡罗方法一般情况下需要大量的样本，因此计算效率低，特别是在混合不确定性下，对于区间变量任意取值，都需要大量的样本进行仿真，因此计算效率极低，在工程中难以适用。