结构可靠性优化设计的证据理论和微分演化方法

2014-09-27 13:58唐和生苏瑜薛松涛邓立新

湖南大学学报·自然科学版 2014年4期

唐和生+苏瑜+薛松涛+邓立新

文章编号：16742974(2014)04003306

收稿日期：20130417

基金项目：国家自然科学基金资助项目（51178337，50708076）；科技部国家重点实验室基础研究项目（SLDRCE11B01）；同济大学土木工程学院光华基金资助项目

作者简介：唐和生（1973-），男，安徽安庆人，同济大学副教授，博士

通讯联系人，E-mail:thstj@tongji.edu.cn

摘要：采用证据理论作为传统概率的替代方法处理不精确的数据信息，提出了基于证据理论的可靠性优化设计方法.该方法针对给定的失效概率许用值Pf，通过计算证据理论的不确定测度Pl(F)，以Pl(F)<Pf作为可靠性约束条件的替代模型进行可靠性优化设计.为了降低基于证据理论不确定量化分析的计算成本，引入微分演化优化算法计算区间边界值.以典型桁架结构形状优化问题为例，验证了所提出方法的准确性和有效性.



关键词：证据理论；微分演化；可靠性优化设计；形状优化



中图分类号：TU318；TU323.4 文献标识码：A

AMethodofReliabilityDesignOptimizationUsingEvidence

TheoryandDifferentialEvolution



TANGHesheng1,2,SUYu1,XUESongtao1,DENGLixin1

（1.ResearchInstituteofStructuralEngineeringandDisasterReduction,TongjiUniv,Shanghai200092,China;

2.StateKeyLaboratoryofDisasterPreventioninCivilEngineering,TongjiUniv,Shanghai200092,China）

Abstract:Anewmethodofreliabilitydesignoptimizationusingtheevidencetheorywasproposed.Evidencetheorywaspresentedasanalternativetotheclassicalprobabilitytheorytohandletheimprecisedatasituation.TheplausibilitymeasurePl(F)basedonevidencetheory,withPl(F)



Keywords:evidencetheory;differentialevolution;reliabilitydesignoptimization;shapeoptimization



基于可靠性的优化设计是不确定性结构优化设计的有效途径.传统的可靠性优化设计常采用概率模型，但概率方法需要足够的统计数据信息来拟合其概率特征，而实际工程中这些数据通常是无法准确获知的，所以传统概率方法面临巨大的挑战.

近年来，国内外很多学者致力于发展非概率的不确定建模手段，并在其基础上提出非概率可靠性优化设计方法，其中由Dempster[1]和Shafer[2]提出的证据理论具有较强的不确定处理能力，已经成为不确定信息表达和量化的有力工具，在多目标识别、信息融合、多属性决策等领域获得了广泛应用[3-6].而基于证据理论的工程可靠性优化分析刚刚起步，并且主要应用于机械与航空领域.Mourelatos等[7]将证据理论应用于内压容器可靠性优化设计中，研究了基于证据理论的失效概率或可靠度指标的计算问题.Bae等[8]运用证据理论解决机械工程中的不确定问题，实现了飞机机翼结构的可靠性优化设计.郭慧昕等[9-11]提出了证据理论和区间分析相结合的可靠度优化设计方法，将此应用于内压容器和气门弹簧的可靠性研究中.尽管已经取得了一些进展，但是证据理论仍然很少应用于实际工程的可靠性优化问题中，计算成本是导致该问题的主要原因.同时，基于证据理论的可靠性优化设计在土木工程中的应用还是一个新课题.

为此,本文采用证据理论处理不确定情况，引入微分演化算法降低证据理论在可靠性优化中的计算成本，提出了证据理论和微分演化算法相结合的可靠性优化设计方法，并将该方法应用于桁架结构的形状优化问题中，来验证本文所提方法的有效性.

1证据理论的基本原理

证据理论是由Dempster和Shafer提出的，又称为DS理论.它是建立在辨识框架上的一种不确定理论，设Θ为辨识框架，它表示关于命题互不相容的所有可能答案的有限集合，类似于概率论中有限的样本空间，幂集2Θ定义为辨识框架Θ中所有子集的集合，证据理论是对幂集元素进行基本概率赋值[12].定义函数m：2Θ→[0,1]，AΘ.当满足：

m(Φ)=0，∑AΘm(A)=1，（1）

称m为框架Θ上的基本信任分配函数(BBA)，m(A)为A的基本信任度，表示证据对A的信任程度，若m(A)>0，则称A为焦元.

若m为框架Θ上的基本信任分配函数，则称由

Bel(A)=∑BAm(B),（2）

Pl(A)=1-Bel()=∑B∩A≠m(B)（3）

所定义的函数Bel：2Θ→[0,1]为Θ上的信任函数，函数Pl：2Θ→[0,1]为Θ上的似然函数，Bel(A)表示对A为真的信任程度，Pl(A)表示对A为非假的信任程度，也称为命题A的似然度，两者之间的关系如图1所示，Bel(A)和Pl(A)提供了概率P(A)的上限和下限.由此可见，以概率论为基础的传统可靠性问题只是证据理论的一个特例.



图1对命题A的不确定描述

Fig.1UncertaintyrepresentationofpropositionA

对于认识不够透彻的不确定参数，可能会有多个专家提出不同的理论或不同的数据来组成多源证据，证据理论可以通过合成规则综合考虑各种证据源的影响.经典的DS合成规则为：假定m1和m2是同一辨识框架Θ上的2个基本信任分配函数，焦元分别为Ai和Bj，则新的基本信任分配函数m为：

m(A)=∑Ai∩Bj=Am1(Ai)m2(Bj)1-K,A≠.（4）

式中：K=∑Ai∩Bj=m1(Ai)m(Bj)，表示证据冲突性的大小.此即为2个证据合成的Dempster法则.当证据冲突比较大时，应选用其他的合成方法[13].

2证据理论和微分演化算法相结合的可靠

性优化设计方法

2.1基于证据理论的可靠性优化设计模型

一般来说，基于概率理论的可靠性优化问题可表述为：

minF(xN,d),

s.t.P{gj(x,d)>g0}

P{gj(x,d)R,

j=1,2,…,k,

dL≤d≤dU,

xL≤xN≤xU．(5)

式中：x=[x1,x2,…,xn]T为n维不确定性量，分别服从一定的概率分布，[xL,xU]为不确定量x的名义值xN的取值区间；d=[d1,d2,…,dm]T为m维确定性量，[dL,dU]为d的取值区间；F和g分别为目标函数和约束函数；g(x,d)>g0表示结构发生失效，g(x,d)<g0表示结构安全，g0为许可的响应值；k为约束条件的数目；P{}表示真实概率；Pf和R分别为结构设计给定的失效概率和可靠度许用值.

当问题中不确定量的认识信息较少或不完整时，上述优化模型中的约束条件不能采用概率理论来建立，此时，可以利用证据理论的不确定建模手段解决这一问题.如图1所示，证据理论用似然函数和信任函数来进行不确定性度量，可以证明[Bel,Pl]是真实概率的区间估计，真实的失效概率或可靠度夹逼在该区间内：

Pl{gj(x,d)>g0}

P{gj(x,d)>g0}

Bel{gj(x,d)R→

P{gj(x,d)R．（6）

由此，将Pl{gj(x,d)>g0}<pf或Bel{gj(x,d)<g0}>R作为传统可靠性约束条件的替代模型，建立基于证据理论的广义可靠度优化设计数学模型：

minF(xN,d)

Pl{gj(x,d)>g0}

g0}>R，

j=1,2,…,k，

dL≤d≤dU，

xL≤xN≤xU.（7）

模型中Pl{gj(x,d)>g0}或Bel{gj(x,d)<g0}是进行非概率可靠性分析的前提，由于结构失效占整个设计空间的比例较小，故以Pl{gj(x,d)>g0}<pf作为约束条件可以提高优化效率，下面就Pl{gj(x,d)>g0}的计算进行详细阐述.

2.2基于微分演化方法的失效似然度计算

对于结构分析中出现的不确定量，证据理论将其表述为区间数.在计算结构失效似然度时，首先根据不确定量的可能取值范围，将其划分为有限个互不相容的基本区间作为辨识框架.以图2为例，不确定量x1的辨识框架X1={x11,x12,x13}，幂集2X1={Ф,{x11},{x12},{x13},{x11,x12},{x11,x13},{x12,x13},X1}，基于证据（专家意见或实测数据）分析，对幂集中的焦元进行基本信任度赋值，得到x1的基本信任分配函数m.



图2不确定参数基本区间

Fig.2Basicintervalsofuncertaintyparameter



然后，对每个不确定变量的焦元进行笛卡尔运算，得到联合焦元区间，以二维不确定参数为例：

C=x1×x2={ck=［x1m,x2n］:x1m∈X1,

x2n∈X2}.（8）

式中:x1m，x2n和ck分别为X1，X2和C的焦元区间.考虑到x1和x2的独立性，二维联合焦元的基本信任度m(ck)=m(x1m)m(x2n).

由于x1m和x2n都是区间，焦元ck在集合上为一矩形，显然，对于n维问题，联合辨识框架中的焦元为n维“超立方体”.令y(x,d)表示结构极限状态函数，结构失效域F为：

F={x:y=g0－g(x,d)<0,x=

［x1,…,xn］∈ck}.（9）

在联合BBA和失效域F的基础上，可根据式(2)和式(3)计算结构失效测度Bel(F)和Pl(F)为：

Bel(F)=∑ckFm(Xc),

Pl(F)=∑ck∩F≠m(Xc)．(10)

可见，在计算失效测度时需要确定联合焦元ck是否满足ckF或ck∩F≠,图3描述了ck对Bel(F)和Pl(F)的贡献．从图3可以看出：1）若ymax>0，ymin>0，则ck∩F=，ck对Pl(F)和Bel(F)没有贡献，即ck不参与Pl(F)和Bel(F)的计算.2）若ymax<0，ymin<0，则ck∩F=ck，ck对Pl(F)和Bel(F)的计算都有贡献.3）若ymax>0,ymin<0，ck仅对Pl(F)有贡献.



图3焦元区间对失效似然度的贡献

Fig.3Focalelementcontribution

toplausibilityofthefailureregion



因此，为准确判断，需要求解y(x,d)在ck对应的“超立方体”域上的极值，即

［ymin,ymax］=［mincky(x,d),mincky(x,d)］.（11）

求解式（11）中区间极值的主要方法有采样法和优化方法，采样法的精度很大程度上取决于采样点数目，计算代价巨大.优化方法会极大降低计算量，但由于不确定量x的焦元区间数目多，而且结构响应并不是简单的显式而是通过有限元分析得到的，故利用传统的优化算法求解复杂多维非凸的极限状态函数y(x,d)在ck上的极值显得非常困难.

近年来仿生智能优化算法被广泛引入到结构优化中，例如模拟退火法(SA)[14]、遗传算法(GA)[15]、微分演化法(DE)[16]等，其中DE是一种新颖的启发式智能算法，采用变异、交叉和选择3项基本操作，通过若干代种群演化操作不断舍弃劣质解，保留优质解，最终获取近似全局最优解.研究表明，微分演化算法在求解非凸、多峰、非线性优化问题中表现出较强的稳健性，同时具有收敛较快的优点[17].因此，本文采用DE提高y(x,d)区间极值的求解速度，如图4所示，从而减少优化设计的计算成本.

图4区间函数极值求解

Fig.4Intervaloptimizationforcomputingbounds

根据以上描述，基于证据理论的可靠性优化设计是利用微分演化算法在满足可靠性约束Pl{g(x,d)>g0}<pf的条件下，寻求结构最优设计变量和最优目标值.同时，基于证据理论的结构约束失效似然度又是通过微分演化算法来提高其计算效率的，可靠性约束条件的计算流程如图5所示.



图5可靠性约束计算流程

Fig.5Flowchartofcalculationofreliabilityconstraint

3算例分析

为了便于比较，取文献[18]中的25杆桁架形状优化进行讨论，结构形式见图6，弹性模量名义值E=68950MPa，作用于桁架上的荷载名义值列于表1，容许拉压应力[σ]=±275.6MPa，各节点三向允许的最大位移为8.89mm.其他参数见文献[18].



图625杆空间桁架结构

Fig.625barspacetrussstructure



表125杆桁架节点荷载名义值

Tab.1Normalvalueofjointloadfor25bartruss

节点号

Fx/kN

Fy/kN

Fz/kN

4.448

-44.48

2.224

2.668

0

文献[18]对该桁架进行了确定性优化，本文在此基础上考虑不确定情况，将外荷载和弹性模量视为不确定的，假定其不确定信息（焦元区间及基本信任度）如表2所示，在2种不确定因素存在的情况下，进行既满足可靠度约束条件又使结构总质量最小的最优设计，该不确定优化问题的数学模型为：

findd=［A1,A2,…,A8,X4,Y4,…,Y8］,

minF(d)=∑8i=1ρAiLi+λM,

s.t.Pl{gi(x,d)<0}

g1(x,d)=278.6-maxσk(x,d),

g2(x,d)=8.89-maxujl(x,d),

x=［F1x,F1y,…,F6z,E］．(12）

式中:d为尺寸和形状设计变量;x为不确定参数;Pl{gi(x,d)<0}

考虑结构允许的失效概率Pf为0.05和0.1二种情况，采用本文所提方法对25杆桁架进行可靠性形状优化.Pf=0.05情况的评价函数收敛曲线和最终形状分别见图7和图8.图9给出最优设计时位移约束函数g2(x,d)的信任度和似然度累计分布曲线.为了与文献[18]的确定性形状优化结果相比较，根据图9，表3详细列出了应力和位移约束失效似然度.

表2外荷载和弹性模量的不确定信息

Tab.2Theuncertaininformationofloadandelasticmodulus

F1x/kN

F1y/kN

F1z/kN

F2y/kN

区间

BPA

区间

BPA

区间

BPA

区间

BPA

[3.54.6]

0.15

[-48.9-44.5]

0.2

[-48.9-44.5]

0.2

[-48.9-44.5]

0.2

[4.24.6]

0.65

[-44.5-40.0]

0.5

[-44.5-40.0]

0.5

[-44.5-40.0]

0.5

[4.65.3]

0.2

[-40.0-33.4]

0.3

[-40.0-33.4]

0.3

[-40.0-33.4]

0.3

F2z/kN

F3x/kN

F6x/kN

E/103MPa

[-48.9-44.5]

0.2

[1.82.3]

0.15

[2.12.8]

0.15

[6065]

0.1

[-44.5-40.0]

0.5

[2.12.3]

0.65

[2.52.8]

0.65

[6570]

0.5

[-40.0-33.4]

0.3

[2.32.7]

0.2

[2.83.2]

0.2

[7080]

0.4

迭代次数

图725杆桁架形状优化的收敛曲线

Fig.7Shapeoptimizationconvergence

historyof25bartruss





图825杆桁架的形状优化结果

Fig.8Optimumshapeof25bartruss

由图7可知，该算法具有很高的计算效率，25杆桁架形状的优化计算在迭代大约100次后已经收敛.从图9和表3可以看出，在Pf=0.1和Pf=0.05两种情况下，位移约束失效概率［P(g2<0)］的不确定区间

g2(x,d)(a)Pf=0.05

g2(x,d)(b)Pf=0.1

图9位移约束的信任度和似然度累积分布

Fig.9Cumulativebeliefandplausibility

distributionfordisplacementconstraint



[Bel,Pl]分别为[0,0.066]和[0,0.036]，失效似然度即概率上界均小于相应的失效概率许用值，满足设计可靠度的要求.由表3可知，由于考虑不确定性的存在，基于DS可靠性优化结果总质量要比确定性优化结果有所增加，但是从失效似然度来看，前者的可靠性（93.4%，96.4%）要明显高于后者（5%）.由此可见，对不确定量的认识信息较少，无法采用概率理论时，证据理论以区间测度[Bel,Pl]代替传统概率单值来描述这种认知不确定显得更为合理.由于基于DS的可靠性优化将使结构具有良好的鲁棒性，有效避免由于错误估计而造成优化结果的偏差.

表325杆空间桁架形状优化结果比较

Tab.3Comparisonofoptimaldesignsforthe25bartruss

设计变量

/mm2

文献[18]

结果/mm2

本文不确定分析

结果/mm2

Pf=0.05

Pf=0.1

A1

64.5

A2

64.5

A3

645

774.2

A4

64.5

A5

64.5

A6

64.5

A7

64.5

A8

580.6

838.7

709.6

X4

949.9

984.6

961.1

Y4

1406.6

1770.5

1375.6

Z4

3283.9

2599.1

3283.7

X8

1315.9

1322.5

1468.0

Y8

3544.8

3400.5

3533.2

总质量/kg

Pl(g1<0)

Pl(g2<0)

53.1

0

0.95

66.5

0

0.036

62.2

0

0.066

4结论

可靠性优化设计中，由于不确定信息较少无法构造精确概率分布时，证据理论代替传统的概率理论进行不确定信息描述是一种理想的选择.该方法用不确定区间测度[Bel,Pl]代替不可知的真实概率来处理不完备的数据信息，以Pl(F)<Pf作为概率可靠性约束条件的替代模型，本文同时引入微分演化算法提高求解联合焦元内反应区间极值的计算效率以及可靠性优化的寻优速度，降低证据理论在可靠性优化中的计算成本.

本文以25杆桁架形状优化为例，在考虑荷载和弹性模量均为不确定的情况下，基于DS进行可靠性优化设计得到了很好的结果.分析结果表明，本文所提方法在实际工程中具有一定的应用前景.

参考文献1］ DEMPSTER A P. Upper and lower probabilities induced by a multiplicand mapping[J]. Annals of Mathematical Statistics, 1967,38:325-339.［2］ SHAFER G. A mathematical theory of evidence [M]. Princeton: Princeton University Press, 1976:9-50.［3］ MCCLEAN S, SCOTNEY B. Using evidence theory for the integration of distributed database[J]. International Journal of Intelligent Systems, 1997,12(10):763-776.［4］ FRANCOIS J, GRANDVALET Y, DENOEUX T, et al. Addendum to resample and combine: an approach to improving uncertainty representation in evidential pattern classification[J]. Information Fusion, 2003(4):75-85.［5］张晓明,王航宇,黄达.基于DS证据理论的多平台协同数据融合[J].计算机工程,2007,33(11):242-243.ZHANG Xiaoming, WANG Hangyu, HUANG Da. Cooperative multiplatform data fusion based on DS evidence theory[J]. Computer Engineering, 2007,33(11):242-243.(In Chinese)［6］郝军,吴炜,杨栈.改进DS算法在船舶汽轮机故障诊断中的应用研究[J].船海工程,2007,36(3):49-51.HAO Jun, WU Wei, YANG Zhan. Application of the improved DS algorithm in fault diagnosis of ships turbinemachinery[J]. Ship & Ocean Engineering, 2007,36(3):49-51.(In Chinese)［7］ MOURELATOS Z P, ZHOU J. A design optimization method using evidence theory[J]． Journal of Mechanical Design, ASME, 2006,128(7):901-908.［8］ BAE H R, GRANDHI R V, CANFIELD R A. Reliabilitybased design optimization under imprecise uncertainty[C]//46th Structures, Structural Dynamics and Materials Conference, AIAA/ASME/ASCE/AHS/ASC, 2005:3301-3313.［9］郭慧昕,刘德顺,胡冠昱,等.证据理论和区间分析相结合的可靠性优化设计方法[J].机械工程学报,2008,44(12)：36-41.GUO Huixin, LIU Deshun, HU Guanyu，et al. Method of reliability design optimization using evidence theory and interval analysis[J]. Chinese Journal of Mechanical Engineering, 2008,44(12):36-41.(In Chinese)［10］郭慧昕,夏立龙,戴娟.基于证据理论的结构失效概率计算方法[J].应用基础与工程科学学报,2008, 16(3):457-464.GUO Huixi, XIA Lilong, DAI Juan. Method for calculating the failure probability of engineering structure using evidence theory[J]. Journal of Basic Science and Engineering, 2008,16(3):457-464.(In Chinese)

［11］任丕顺．基于证据理论的气门弹簧可靠性设计优化[J]．装备制造技术，2012（5）:49-52.REN Pishun. The reliability optimal design of the valve spring based on evidence theory[J]. Equipment Manufacturing Technology, 2010（5）:49-51.(In Chinese)［12］段新生.证据理论与决策、人工智能[M].北京:中国人民大学出版社,1993:13-32.DUAN Xinsheng. Evidence theory and decision & artificial intelligence [M]. Beijing: China Renmin University Press, 1993:13-32. (In Chinese)［13］SENTZ K, FERSON S. Combination of evidence in dempstershafer theory, SAND20020835 report［R］. Sandia: Sandia National Laboratories, 2002.［14］BALLING R. Optimal steel frame design by simulated annealing[J]. Journal of Structural Engineering, 1991,117(6):1780-1795.［15］GOLDBERG D, KUO C. Genetic algorithms in pipeline optimization[J]. Journal of Computing in Civil Engineering, 1987,1(2):128-141.［16］唐和生,范德伟,王兆亮,等.桁架尺寸优化微分演化算法[J].湖南大学学报：自然科学版,2011,38(11):13-18.TANG Hesheng, FAN Dewei, WANG Zhaoliang, et al. Differential evolution algorithm to size the optimization of truss structures[J]. Journal of Hunan University: Natural Sciences, 2011,38(11):13-18.(In Chinese)［17］STORN R, PRICE K. Differential evolution—a simple and efficient adaptive scheme for global optimization over continuous spaces[J]. Journal of Global Optimization, 1997,11(4):341-359.［18］GHOLIZADEH S, BARATI H. A comparative study of three metaheuristics for optimum design of trusses[J]. International Journal of Optimization in Civil Engineering, 2012,3:423-441.endprint

湖南大学学报·自然科学版2014年4期

湖南大学学报·自然科学版的其它文章: 长周期随机地震作用下超大跨斜拉桥的行波效应分析; 隔层桁架式立体停车结构地震倒塌易损性分析; 竖向新型连接装配式剪力墙抗震性能试验研究; 土海底沉管隧道体系三维地震响应分析; CFRP中等约束钢筋混凝土方柱反复受压本构模型; 铝合金板式节点受弯承载力试验研究