东北大学学报:自然科学版  2020, Vol. 41 Issue (9): 1268-1273  
0

引用本文 [复制中英文]

张毅博, 孙志礼, 赵中强, 赵经武. 基于GPC的环肋耐压圆柱壳结构失稳概率分析[J]. 东北大学学报:自然科学版, 2020, 41(9): 1268-1273.
[复制中文]
ZHANG Yi-bo, SUN Zhi-li, ZHAO Zhong-qiang, ZHAO Jing-wu. Instability Probability Analysis of Ring Stiffened Pressure Cylindrical Shell Structures Based on GPC[J]. Journal of Northeastern University Nature Science, 2020, 41(9): 1268-1273. DOI: 10.12068/j.issn.1005-3026.2020.09.009.
[复制英文]

基金项目

国家自然科学基金资助项目(51775097)

作者简介

张毅博(1989-),男,河北石家庄人, 东北大学博士研究生;
孙志礼(1957-),男,山东巨野人,东北大学教授,博士生导师。

文章历史

收稿日期:2020-02-05
基于GPC的环肋耐压圆柱壳结构失稳概率分析
张毅博 1, 孙志礼 1, 赵中强 2, 赵经武 3     
1. 东北大学 机械工程与自动化学院, 辽宁 沈阳 110819;
2. 中车唐山机车车辆有限公司,河北 唐山 063000;
3. 中国人民解放军第93107部队, 辽宁 沈阳 110141
摘要:为评估具有小失效概率特性的深潜环肋耐压圆柱壳结构的失稳概率,提出了一种新的基于高斯过程分类和重要抽样的自适应分析方法.该方法通过引入马尔科夫链蒙特卡洛法和欧式距离,开发了一种新的考虑预测不确定性和取样均匀性的自适应试验设计策略,以便更高效地构造高斯过程分类器;采用核密度估计构造准最优重要抽样密度函数;基于失效概率估计的稳定性,提出了一种更精确的迭代停止准则.通过某一分段函数验证了所提分析方法的准确性及高效性.应用所提方法得到某深潜环肋耐压圆柱壳结构的失稳概率约为8.242×10-5.
关键词失稳概率    环肋耐压圆柱壳    高斯过程分类    小失效概率    自适应试验设计    
Instability Probability Analysis of Ring Stiffened Pressure Cylindrical Shell Structures Based on GPC
ZHANG Yi-bo 1, SUN Zhi-li 1, ZHAO Zhong-qiang 2, ZHAO Jing-wu 3     
1. School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819, China;
2. CRRC Tangshan Co., Ltd., Tangshan 063000, China;
3. Unit 93107 of the Chinese People′s Liberation Army, Shenyang 110141, China
Abstract: To evaluate the instability probability of deep submergence ring stiffened pressure cylindrical shell structures with small failure probability, an innovative adaptive analysis method based on Gaussian process classification (GPC) and importance sampling (IS) was proposed. By introducing the Markov Chain Monte Carlo (MCMC) and the Euclidean distance, a new adaptive strategy for design of experiments (DoE), considering the prediction uncertainty and the sampling uniformity, was developed to establish the Gaussian process classifier more efficiently. Furthermore, the quasi-optimal importance sampling density function was constructed by adopting the kernel density estimation (KDE). Based on the stability of failure probability estimation, a more accurate stopping criterion was also proposed. A piecewise function was utilized to verify the accuracy and efficiency of the proposed analysis method. The instability probability of a deep submergence ring stiffened pressure cylindrical shell structure obtained by the proposed method is about 8.242×10-5.
Key words: instability probability    ring stiffened pressure cylindrical shell    Gaussian process classification (GPC)    small failure probability    adaptive design of experiments    

随着对隐蔽性要求的不断提高,潜艇下潜深度逐渐增加,环肋耐压圆柱壳结构作为潜艇的基本结构,其稳定性问题随静水压力的增大而愈发突出[1].因此,开展深潜环肋耐压圆柱壳结构失稳概率分析具有重要意义和工程价值.

可靠性分析的是产品在特定条件下完成特定功能的能力.因此,本文采用可靠性方法评估环肋耐压圆柱壳结构的失稳概率.在可靠性领域,结构状态通常被分为两类:安全和失效.近年来,为提高可靠性分析效率,诸多学者提出了各种基于代理模型(如响应面[2]、Kriging模型[3]、支持向量机SVM[4]、人工神经网络ANN等)的分析方法.然而,现有可靠性分析方法主要针对具有连续响应的结构,而对于具有失稳(即屈曲)破坏的拱结构、桁架结构、薄壳结构等的分析研究甚少.尽管可以采用SVM和ANN处理结构状态二分类问题,然而SVM和ANN均需要大量已标注样本才能得到较为精确的分类结果.此外,SVM仍存在核函数、核参数以及损失函数较难选取等问题.

作为高斯过程的一个分支,高斯过程分类(GPC)不仅具有严格的统计理论基础,还能自适应获取超参数[5].然而,GPC执行推理的时间随训练样本数的增加呈立方增长.为此,学者们提出了各种自适应试验设计(DoE)策略,如Kapoor等[6]将最易分类错误的点(MEMP)作为新训练样本;Peng等[7]将最可能失效点(MPP)添加到DoE中.然而,现有DoE策略选择出来的训练样本容易聚集,即引起不必要的功能函数评估.此外,对于具有小失效概率结构的可靠性评估,尽管Yang等[8]提出了一种双循环且不依赖于MPP的基于Krging和重要抽样(IS)的可靠性分析方法,然而目前尚未发现相关基于GPC和IS的可靠性分析方法.

为此,本文提出了一种单循环的基于GPC和IS的自适应分析方法,既避免了双循环中优化算法易陷入局部最优解的缺陷,又保证了在分类面附近区域均匀取样.通过某一分段函数验证了所提方法的高效性与准确性,并将其用于评估某深潜环肋耐压圆柱壳结构的失稳概率.

1 基本理论 1.1 高斯过程分类

高斯过程分类通常使用类标签+1或-1处理二元分类问题,+1代表结构安全,-1代表结构失效.给定N个训练样本X(X=[x1, …, xN])及其类标签Y(Y=[y1, …, yN]),GPC通过建立XY间的映射关系(即潜函数f(x))来进行分类.

为对未知样本x*分类,需由式(1)计算x*相应潜变量f*的后验分布.

(1)

式中:f为训练样本X对应的潜向量(f=[f(x1), f(x2), …, f(xN)]T);p(f*|X, x*, f)为f*的条件先验分布;p(f|X, Y)为f的后验分布.

为求上述积分,GPC假定条件分布p(f|X)服从高斯分布,即

(2)

式中,K为协方差矩阵(Kij=k(xi, xj));k(·)为协方差函数,本文选用应用最广泛的平方指数协方差函数(也称高斯核函数)来计算K.

(3)

式中:σf用来控制局部相关程度;指数项用来表征xixj间的距离相关性,即如果xixj的距离相对于距离尺度l很小,则它们高相关,否则低相关.

由高斯过程可知,f*f的联合分布也服从高斯分布,即

(4)

式中:Kx*=[k(x*, x1), …, k(x*, xN)]Tkx*=k(x*, x*).

因此,条件先验分布p(f*|X, x*, f)为

(5)

由贝叶斯规则可得f的后验分布p(f|X, Y)为

(6)

式中:; p(Y|X)=∫p(f|X)p(Y|f)df.

将式(5)和式(6)代入式(1)中,则样本x*的类标签y*为+1的概率为

(7)

然而,式(1)和式(7)没有解析解,需采用数值近似法求解上述积分,本文选用Laplace近似法.Laplace采用高斯分布近似p(f*|X, Y, x*),即

(8)

式中:

因此,式(7)中的分类概率近似为

(9)

式中,σ(·)为响应函数,本文选择标准正态分布的累积分布函数作为响应函数.

显然,分类概率0.5为结构安全与否的界限.

1.2 基于重要抽样的可靠性分析理论

在可靠性分析中,失效概率定义为

(10)

式中:G(x)为结构的功能函数;fX(x)为x的联合概率密度函数;x=[x1, …, xM],M为影响结构性能的不确定因素(即随机变量)的个数.

尽管可以采用各种基于Monte Carlo采样的方法计算失效概率,但当失效概率很小时(如Pf < 10-4),需大量的Monte Carlo样本(>107)才能得到较为精确的结果.为此,需采用重要抽样等方差缩减技术来减少可靠性评估中随机抽样样本的数量.根据重要抽样,式(10)可改写为

(11)

式中:If(x)为失效指示函数,即时,If(x)=1,否则If(x)=0;h(x)为重要抽样密度函数;NIS为重要抽样样本数.

的方差和变异系数分别为

(12)
(13)
2 基于GPC和IS的自适应分析方法 2.1 所提自适应DoE策略

现有基于GPC的自适应DoE策略主要包含两种:基于最可能失效点(GPC+MPP)和基于最易分类错误点(GPC+MEMP).

在可靠性分析中,MPP是对失效概率影响最大的点,且MPP附近区域高斯过程分类器的拟合精度直接影响着失效概率估计的准确性.因此,文献[7]将MPP作为新训练样本,其可由式(14)近似获取.

(14)

由式(9)可知,分类概率越接近0.5,其类标签越容易分类错误.因此,文献[6]将满足式(15)的点,即最易分类错误的点MEMP作为新训练样本.

(15)

然而,上述选点策略存在一些缺陷:1)他们的候补样本集是基于Monte Carlo采样得到的,难以应用于小失效概率评估;2)选择出来的训练样本容易聚集,从而增加不必要的功能函数评估次数.为此,本文提出采用马尔科夫链蒙特卡洛法(MCMC)和欧式距离来避免上述缺陷.

由上述可知,高斯过程分类器在分类概率为0.5处的拟合精度对失效概率估计的准确性起决定性影响.因此,应将分类概率为0.5的样本作为候补样本.为此,本文采用MCMC法产生NCS个满足式(16)的样本并将其作为候补样本集.

(16)

式中,[ε]取为0.01.

为避免训练样本聚集,本文提出采用欧式距离来保证取样均匀性,其中欧式距离为

(17)

式中:xc为候补样本;xt为现有训练样本.

因此,本文将满足式(18)的点作为新训练样本,

(18)

式中,Φ, d分别为Ф(x)和d标准化后的值,标准化的主要目的是消除他们之间数量级差异的影响.Φ, d可由式(19)得到,

(19)

式中,dminxc距DoE中所有训练样本的最短距离.

2.2 基于核密度估计的重要抽样

重要抽样的核心是重要抽样密度函数h(x)的构造,理论上最优h(x)为

(20)

式(20)因包含真实失效概率而难以获取.为此,本文引入核密度估计(KDE)[8]来构造准最优.

基于构造的高斯过程分类器,通过MCMC产生Nk个失效样本(本文取Nk=2 000),则由核密度估计得到的准最优

(21)

式中:ω为带宽;K(·)为核概率密度函数,本文选取正态分布的概率密度函数作为K(·),其表达式为

(22)

带宽ω可由最小化平均积分平方误差获取,

(23)
2.3 停止准则

现有停止准则主要有两种,然而它们由于其固有缺陷而难以被工程实际所采纳.

1) 功能函数调用次数达到某一规定的阈值[9].显然,若规定的阈值过小,则易错估失效概率;反之,引起不必要的功能函数评估.

2) 前后两次迭代失效概率估计的相对误差小于某一阈值[7],即

(24)

显然,若选择的新训练样本对高斯过程分类器的精度改善甚微,则提前终止迭代且得到错误的失效概率估计.

为此,本文基于失效概率估计的稳定性,提出了一种新的、更加精确的停止准则,即

(25)

式中:为第i次迭代及后续Np-1次迭代得到的失效概率估计的平均值;(a∈{i, i+1, …, i+NP-1})为这NP次迭代中某次迭代得到的失效概率估计;ε为最大相对稳定性.只要ε小于阈值[δ],即失效概率趋于稳定,则停止迭代,本文取[δ]为0.03.此外,NP不宜过大或过小,本文令NP为10.

的相对误差定义为

(26)

式中,为由基于核密度估计的重要抽样(KDE-IS)得到的失效概率,本文将其作为精度比较的“真值”.

2.4 所提自适应分析方法

本文所提自适应分析方法的基本步骤如下:

步骤 1  通过Nataf变换将非独立随机变量转换为独立标准正态变量,采用拉丁超立方抽样随机生成N个初始训练样本SDoE=[x1, x2, …, xN]并调用真实功能函数获取类标签YDoE=[y1, y2, …, yN].

步骤 2  基于现有DoE,构造高斯过程分类器,采用MCMC法产生Nk个失效样本,并由式(21)~(23)构造准最优重要抽样密度函数.根据式(11)计算失效概率,如果不满足停止条件(25),则继续下一步,否则跳到步骤4.

步骤 3  基于当前高斯过程分类器,应用MCMC法生成NCS个满足式(16)的候补样本(本文令NCS=3 000).根据式(18)选择新训练样本,调用功能函数获取其对应的真实类标签,并将新训练样本添加到DoE中,返回步骤2.

步骤 4  由式(13)计算的变异系数,如果≤0.03,则停止迭代并输出相关结果,否则扩大重要抽样样本数并返回步骤2.

3 算例验证

以一分段函数为例,验证所提自适应分析方法的准确性及高效性,其表达式为[10-11]

(27)

式中,x1x2相互独立且均服从标准正态分布.

取12个拉丁超立方样本作为初始训练样本,且每种方法各运行10次以消除由不同初始训练样本及候补样本的随机性所产生的影响,其中GPC+MPP和GPC+MEMP以基于核密度估计的重要抽样(KDE-IS)生成的重要抽样样本作为候补样本.不同方法运行10次的平均结果如表 1所示.此外,图 1a1b还分别显示了由不同方法得到的失效概率估计和最大相对稳定性ε随迭代次数Nit的变化趋势.显然,由表 1图 1b可知,所提方法在满足失效概率分析精度要求的同时,还需要更少的功能函数评估次数.图 2比较了不同方法在调用真实功能函数28次时DoE中训练样本的分布及预测的高斯过程分类边界.显然,所提自适应分析方法选择的训练样本分布更均匀且能更好地收敛到真实分类边界,即需要调用真实功能函数次数更少.

表 1 不同方法运行10次的平均结果 Table 1 Average results of 10 runs of different methods
图 1 εNit的变化 Fig.1 Variations of and ε with Nit (a)— Nit的变化;(b)—εNit的变化.
图 2 不同方法下训练样本的分布及预测的高斯过程分类边界 Fig.2 Distribution of training samples and predicted Gaussian process classification boundary with different methods (a)—GPC+MPP;(b)—GPC+MEMP; (c)—所提方法.
4 环肋耐压圆柱壳失稳概率分析

某深潜环肋耐压圆柱壳结构采用矩形肋骨,其截面如图 3所示,其中,r为耐压壳体内径,t为耐压壳体厚度,l为肋骨间距,a为肋骨高度,b为肋骨宽度.各参数及环肋耐压圆柱壳所受外部压力p均服正态分布,相应分布参数如表 2所示[12].

图 3 矩形环肋圆柱壳结构截面形状 Fig.3 Section of cylindrical shell with rectangular ring ribs
表 2 各随机变量的分布特征 Table 2 Distribution of random variables

为评估该型环肋耐压圆柱壳结构的失稳概率,需采用ANSYS中的更符合工程实际和更具应用价值的非线性屈曲分析来分析其稳定性.由于环肋耐压圆柱壳结构的对称性,仅建立其1/4有限元模型,如图 4所示.模型材料为某高强度钢,其弹性模量为2.0×105 MPa,泊松比为0.3,屈服强度为785 MPa.圆柱壳采用壳单元shell181建模,矩形肋骨采用beam188单元建模.模型一端固定约束,另一端允许轴向位移,1/4截面处施加对称约束,壳结构外表面上施加压力载荷.

图 4 四分之一有限元模型 Fig.4 Quarter finite element model

由ANSYS中非线性屈曲分析得到的该型环肋耐压圆柱壳结构在稳定和失稳时的节点位移云图分别如图 5图 6所示.其中,图 6中该型环肋耐压圆柱壳结构的失稳类型为整体失稳.

图 5 稳定时节点位移云图 Fig.5 Cloud diagram of node displacement at stability
图 6 失稳时节点位移云图 Fig.6 Cloud diagram of node displacement at instability

通过建立该型环肋耐压圆柱壳结构参数化有限元模型及MATLAB和ANSYS的联合仿真,基于所提分析方法,评估其失稳概率.为此,在标准正态空间内通过拉丁超立方抽样随机选取30个初始训练样本并调用ANSYS仿真得到其相应最大节点位移;构造初始高斯过程分类器,根据2.1节所提自适应DoE策略选择新训练样本,调用ANSYS获取最大节点位移并加入到DoE中.经过238次循环仿真后,满足收敛条件,即迭代停止,所得该型深潜环肋耐压圆柱壳结构的失稳概率约为8.242×10-5.

5 结论

1) 通过某一分段函数表明,所提方法在满足失效概率分析精度要求的同时还极大地提高了分析效率,避免了不必要的功能函数评估.

2) 将所提自适应分析方法应用于评估某型深潜环肋耐压圆柱壳结构的失稳概率,结果表明:所提方法能够很好地适用于工程实际,这为预防潜艇失稳、缩短新型潜艇的研制周期及解决潜艇下潜的稳定性评估提供了一种切实可行的分析方法.

参考文献
[1]
季林帅.深潜耐压圆柱壳极限承载力研究[D].镇江: 江苏科技大学, 2015.
(Ji Lin-shuai.Research on the ultimate bearing capacity of deep submergence pressure cylinder shell[D].Zhenjiang: Jiangsu University of Science and Technology, 2015. )
[2]
Guimarães H, Matos J C, Henriques A A. An innovative adaptive sparse response surface method for structural reliability analysis[J]. Structural Safety, 2018, 73: 12-28. DOI:10.1016/j.strusafe.2018.02.001
[3]
Zhang Y B, Sun Z L, Yan Y T, et al. An efficient adaptive reliability analysis method based on Kriging and weighted average misclassification rate improvement[J]. IEEE Access, 2019, 7(1): 94954-94965.
[4]
Roy A, Manna R, Chakraborty S. Support vector regression based metamodeling for structural reliability analysis[J]. Probabilistic Engineering Mechanics, 2019, 55: 78-89. DOI:10.1016/j.probengmech.2018.11.001
[5]
Anh N T N, Salim B A, Lam P S. A scalable hierarchical Gaussian process classifier[J]. IEEE Transactions on Signal Processing, 2019, 67(11): 3042-3057. DOI:10.1109/TSP.2019.2911251
[6]
Kapoor A, Grauman K, Urtasun R, et al. Gaussian processes for object categorization[J]. International Journal of Computer Vision, 2009, 88(2): 169-188.
[7]
Peng L F, Su G S, Zhao W. Fast analysis of structural reliability using Gaussian process classification based dynamic response surface method[J]. Applied Mechanics and Materials, 2014, 501/502/503/504: 1067-1070.
[8]
Yang X F, Liu Y S, Mi C Y, et al. Active learning Kriging model combining with kernel-density-estimation-based importance sampling method for the estimation of low failure probability[J]. Journal of Mechanical Design, 2018, 140(5): 051402. DOI:10.1115/1.4039339
[9]
Li X, Gong C L, Gu L X, et al. A sequential surrogate method for reliability analysis based on radial basis function[J]. Structural Safety, 2018, 73: 42-53. DOI:10.1016/j.strusafe.2018.02.005
[10]
Guan X L, Melchers R E. Effect of response surface parameter variation on structural reliability estimates[J]. Structural Safety, 2001, 23(4): 429-444.
[11]
Yuan X K, Lu Z Z, Zhou C C, et al. A novel adaptive importance sampling algorithm based on Markov chain and low-discrepancy sequence[J]. Aerospace Science and Technology, 2013, 29(1): 253-261. DOI:10.1016/j.ast.2013.03.008
[12]
毛昭勇, 卫超, 樊瑜, 等. 基于多失效模式的水下航行器耐压壳体结构可靠性研究[J]. 机械强度, 2013, 35(6): 783-788.
(Mao Zhao-yong, Wei Chao, Fan Yu, et al. Study on reliability of the cylindrical pressure hull with multiple failure modes correlated for underwater vehicles[J]. Journal of Mechanical Strength, 2013, 35(6): 783-788.)