CN105956336B - 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法 - Google Patents

一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法 Download PDF

Info

Publication number
CN105956336B
CN105956336B CN201610407978.9A CN201610407978A CN105956336B CN 105956336 B CN105956336 B CN 105956336B CN 201610407978 A CN201610407978 A CN 201610407978A CN 105956336 B CN105956336 B CN 105956336B
Authority
CN
China
Prior art keywords
ready
vector
function
built structure
design
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610407978.9A
Other languages
English (en)
Other versions
CN105956336A (zh
Inventor
王磊
王晓军
王睿星
耿新宇
樊维超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201610407978.9A priority Critical patent/CN105956336B/zh
Publication of CN105956336A publication Critical patent/CN105956336A/zh
Application granted granted Critical
Publication of CN105956336B publication Critical patent/CN105956336B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法。该方法综合考虑实际工程结构中的多源不确定性情况,建立了一种有效评估静、动态混合不确定性参数对结构全寿命期内安全性能的量化指标,并以此为设计的约束条件,完成了针对拟建结构的非概率可靠性优化。首先,利用区间数和区间过程合理表征静、动态混合不确定性的本征规律;其次,结合区间数学运算法则,构建基于拟建结构响应历程的动态功能函数,并借助离散化手段,定义和解析功能函数的数字特征;利用首次穿越理论,定义并给出结构混合可靠性指标的显式表达式;最后,建立基于混合可靠性的拟建结构非概率可靠性优化模型,并结合智能蚁群算法,实现优化历程的快速迭代与稳健求解。

Description

一种面向静动态混合不确定性的拟建结构非概率可靠性优化 设计方法
技术领域
本发明涉及多源不确定性条件下拟建结构的综合优化设计技术领域,特别涉及考虑静态、动态混合不确定性作用下结构安全态势的合理评估与优化设计,为指导工程技术人员制定大型复杂拟建结构的设计方案具有可借鉴的理论与应用价值。
背景技术
拟建结构(如航空航天飞行器、桥梁、建筑等)在人类生活、生产、国防建设及社会进步等发展方面扮演着极其重要的角色,是人类生存必需的物质基础。然而,考虑到结构长期使用过程中(飞行器服役期10-20年,桥梁、建筑服役期20-50年),在内部与外部、人为与自然因素作用下,随着时间的推移将发生不可逆转的材料老化和结构损伤。这种过程的累积将导致结构性能急剧劣化、耐久性远低于设计预期;加之服役环境下载荷的多变性,各种因素的叠加严重影响结构使用周期内的安全性能。针对工程上拟建结构静动力特性的可靠性设计可以有效寻找到确保结构在未来整个生命周期内安全的综合性能最佳方案,相关研究是目前工程领域所关注的核心。
近几十年来,随着工程结构的系统化进程,资源的损耗速度已大大超过人们的预估。如何在结构设计之初,尤其是针对大型结构,在确保其安全的基础上,合理利用有限资源实现配置优化,是当前理论研究和工程应用领域中一个亟待解决的热点问题。随着计算机技术的飞速发展,基于可靠性的优化设计成为目前国内外都在积极探索和研究的重要课题。然而,由于实际工程结构中载荷不确定性、材料分散性等体现出的静、动态特性,基于混合不确定性的可靠性度量指标计算及基于可靠性的优化求解过程相对复杂;此外,受到小样本条件的制约,现有的基于概率统计理论的可靠性优化方法无法应用于真实工程结构设计方案的制定。因此,静、动态多源不确定性条件下针对拟建结构的非概率可靠性优化设计方法目前相对较少,理论基础薄弱。
本发明将重点探究材料分散性、动态载荷不确定性综合作用下拟建结构的优化问题,利用所建立的混合可靠性指标作为约束条件,确保结构安全的同时,实现轻量化的最终目的。本发明所涉及的方法将为拟建结构的精细化方案设计提供技术保障。
发明内容
本发明要解决的技术问题是:克服现有技术的不足,提供一种面向静、动态混合不确定性的拟建结构非概率可靠性优化设计方法,充分考虑实际工程问题中普遍存在的静态不确定性(如材料分散性等)和动态不确定性(如载荷多变性等),以提出的混合可靠性优化设计模型为拟建结构方案更新的指导策略,所得到的设计结果更符合真实情况,工程适用性更强。
本发明采用的技术方案为:一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法,该方法实现步骤如下:
第一步:考虑存在于拟建结构中的多源不确定性参数,利用区间数向量X={Xi|i=1,2,...,n1}T表征静态不确定性信息,n1为区间数向量包含的维度,Xi为区间数向量中的第i个元素,静态不确定性信息通常指材料参数的分散性,区间数向量X的本征规律由均值向量Xc和半径向量Xr来体现,即:
X∈[Xc-Xr,Xc+Xr]
其具体表达式为:
其中,Xi 是区间数变量Xi的下界和上界。
利用区间过程向量Y(t)={Yj(t)|j=1,2,...,n2}T表征动态不确定性信息,n2为区间过程向量包含的维度,Yj(t)为区间过程向量中的第j个元素,动态不确定性信息通常指动力载荷的波动性。区间过程向量Y(t)的边界规律由均值过程向量Yc(t)和半径过程向量Yr(t)来体现,即:
Y(t)∈[Yc(t)-Yr(t),Yc(t)+Yr(t)]
其具体表达式为:
此外,区间过程向量Y(t)任意不同时刻t1和t2下的相关性规律可由自相关系数向量ρY(t1,t2)表征,即:
其中,Yj(t)是区间过程Yj(t)的下界和上界,表示任意不同时刻t1和t2下区间过程Yj(t)的自相关系数函数,为对应的自协方差函数。
第二步:利用第一步提出的n1维静态不确定性区间数向量X和n2维动态不确定性区间过程向量Y(t),构建n1+n2维混合不确定性条件下拟建结构的线性动态化功能函数:
G(t,X,Y(t),d)=f(t,d)+a(d)X+b(t,d)Y(t)
其中,d表示可设计的结构形状与尺寸向量,如工程梁、板及三维实体结构,f(t,d)为确定性函数,a(d)={ai(d),i=1,2,...,n1}表示时不变系数向量,b(t,d)={bj(t,d),j=1,2,...,n2}表示时变系数向量;基于区间数学运算法则,利用第一步中定义的均值向量Xc和半径向量Xr,均值过程向量Yc(t)和半径过程向量Yr(t)以及自相关系数向量ρY(t1,t2)可进一步得到描述功能函数特征的均值函数Gc(t,X,Y(t),d)、半径函数Gr(t,X,Y(t),d)和自相关系数函数ρG(t1,t2)的数学表达式如下:
以及
其中,CovG(t1,t2)表示功能函数G(t,X,Y(t),d)的协方差函数。
第三步:基于首次穿越理论并结合第二步构建的线性动态化功能函数,利用时间离散化方法,定义如下穿越事件Ek发生的可能度指标:
PI{Ek}=Pos{G(kΔt,X,Y(kΔt),d)>0∩G((k+1)Δt,X,Y((k+1)Δt),d)≤0}
其中,PI{·}表示事件发生的可能度,G(kΔt,X,Y(kΔt),d)>0表示拟建结构在kΔt时刻安全,即功能函数大于零,G((k+1)Δt,X,Y((k+1)Δt),d)≤0表示拟建结构在(k+1)Δt时刻失效,即功能函数小于等于零,符号“∩”表示事件的交运算,k为计数指标,Δt表示时间增量。这里,事件Ek即表示为拟建结构在时间段[kΔt,(k+1)Δt]内发生了一次穿越,Δt通常为一微小量,其取值设定为完整生命周期T的1/1000。
第四步:遍历所有时间段内的穿越可能度PI{Ek},计算基于静、动态不确定性的拟建结构混合可靠度计算指标:
其中,分别表示整个生命周期T内拟建结构的混合可靠度和混合失效度,Pos(0)表示结构在构建之初即发生失效的可能度,求解上式即可实现拟建结构动力安全态势的有效评估;
第五步:以第四步定义的拟建结构混合可靠度和常规确定性约束作为限制条件,以结构的轻量化功能作为优化目标,以结构的形状与尺寸参数作为优化的设计变量,综合考虑静、动态不确定性(如材料分散性与载荷不确定性)的综合效应,构建出如下面向拟建结构减重优化的混合可靠性设计模型:
find:d
min M(T,Xc,Yc(t),d)
s.t.gi*(t,Xc,Yc(t),d)≤0,i*=1,2,...,l1
其中,M(T,Xc,Yc(t),d)为结构的名义质量函数,表示第i*个确定性约束条件,表示第j*个混合可靠性约束条件,和PI{Ek,j*}分别表示第j*个可靠性约束下结构初始失效的可能度和时间段[kΔt,(k+1)Δt]内结构发生穿越失效的可能度,为可靠性设计的许用值,i*和j*为计数指标,l1和l2分别表示确定性约束与可靠性约束的个数,Ωd为设计变量的可行域。上述优化模型通过嵌入智能蚁群求解算法,完成设计变量的更新并实现整个迭代历程的收敛。
第六步:迭代过程中,如果当前设计不满足第五步混合可靠性设计模型中拟建结构混合可靠度和常规确定性约束,或者尽管满足上述限制条件,但当前设计并不是最优设计,即相较于迭代过程中的前一组设计,目标函数的相对变化百分比大于预设值ξ时,决定设计变量取值的种群重置更新,将已经完成迭代次数的值增加1,并返回第三步,否则,进行第七步;
第七步:如果相较于迭代过程中的前一组设计,目标函数的相对变化百分比小于预设值ξ时,终止计算,将得到的全局最优设计方案中的变量参数作为最终的拟建结构设计方案。
本发明与现有技术相比的优点在于:本发明提供了同时考虑静态、动态混合不确定性效应下的拟建结构非概率可靠性优化设计的新思路,弥补和完善了传统基于单源不确定性的结构静态可靠性设计以及基于概率统计模型的结构时变可靠性设计理论方法的局限性。所构建的拟建结构优化设计模型,一方面综合评估了含多源不确定性结构全寿命周期内的安全性问题,另一方面,以建立的混合可靠性指标作为约束条件,实现了结构安全前提下减重效果的最优,为拟建结构精细化设计方案的规划奠定了一定的理论基础。
附图说明
图1是本发明针对静、动态混合不确定性下拟建结构的非概率可靠性优化设计流程图;
图2是本发明建立的线性动态功能函数在微小时间段内的几何可行域示意图;
图3是本发明定义的穿越事件发生可能度的计算方法示意图;
图4是本发明实施例中拟建的类X-37B机翼结构模型示意图;
图5是本发明实施例中结构可靠性优化设计目标函数迭代历程示意图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
如图1所示,本发明提出了一种面向静、动态混合不确定性的拟建结构非概率可靠性优化设计方法,包括以下步骤:
(1)考虑存在于拟建结构中的多源不确定性参数,利用区间数向量X={Xi|i=1,2,...,n1}T表征静态不确定性信息,n1为区间数向量包含的维度,Xi为区间数向量中的第i个元素,静态不确定性信息通常指材料参数的分散性,区间数向量X的本征规律由均值向量Xc和半径向量Xr来体现,即:
X∈[Xc-Xr,Xc+Xr]
其具体表达式为:
其中,Xi 是区间数变量Xi的下界和上界。
利用区间过程向量Y(t)={Yj(t)|j=1,2,...,n2}T表征动态不确定性信息,n2为区间过程向量包含的维度,Yj(t)为区间过程向量中的第j个元素,动态不确定性信息通常指动力载荷的波动性。区间过程向量Y(t)的边界规律由均值过程向量Yc(t)和半径过程向量Yr(t)来体现,即:
Y(t)∈[Yc(t)-Yr(t),Yc(t)+Yr(t)]
其具体表达式为:
此外,区间过程向量Y(t)任意不同时刻t1和t2下的相关性规律由自相关系数向量ρY(t1,t2)表征,即:
其中,Yj(t)是区间过程Yj(t)的下界和上界,表示任意不同时刻t1和t2下区间过程Yj(t)的自相关系数函数,为对应的自协方差函数。
(2)利用第一步提出的n1维静态不确定性区间数向量X和n2维动态不确定性区间过程向量Y(t),构建n1+n2维混合不确定性条件下拟建结构的线性动态化功能函数:
G(t,X,Y(t),d)=f(t,d)+a(d)X+b(t,d)Y(t)
其中,d表示可设计的结构形状与尺寸向量,如工程梁、板及三维实体结构,f(t,d)为确定性函数,a(d)={ai(d),i=1,2,...,n1}表示时不变系数向量,b(t,d)={bj(t,d),j=1,2,...,n2}表示时变系数向量;基于区间数学运算法则,利用第一步中定义的均值向量Xc和半径向量Xr,均值过程向量Yc(t)和半径过程向量Yr(t)以及自相关系数向量ρY(t1,t2)可进一步得到描述功能函数特征的均值函数Gc(t,X,Y(t),d)、半径函数Gr(t,X,Y(t),d)和自相关系数函数ρG(t1,t2)的数学表达式如下:
以及
其中,CovG(t1,t2)表示功能函数G(t,X,Y(t),d)的协方差函数,具体表达式如下:
显然,基于静、动态混合不确定性的线性功能函数的数字特征可以通过上面公式直接获得。一旦G(t,X,Y(t),d)的数学表达式确知后,功能函数的不确定性变换特征亦可完全确定。
(3)基于首次穿越理论并结合第二步构建的线性动态化功能函数,利用时间离散化方法,定义如下穿越事件Ek发生的可能度指标:
PI{Ek}=Pos{G(kΔt,X,Y(kΔt),d)>0∩G((k+1)Δt,X,Y((k+1)Δt),d)≤0}
其中,PI{·}表示事件发生的可能度,G(kΔt,X,Y(kΔt),d)>0表示拟建结构在kΔt时刻安全,即功能函数大于零,G((k+1)Δt,X,Y((k+1)Δt),d)≤0表示拟建结构在(k+1)Δt时刻失效,即功能函数小于等于零,符号“∩”表示事件的交运算,k为计数指标,Δt表示时间增量。这里,事件Ek即表示为拟建结构在时间段[kΔt,(k+1)Δt]内发生了一次穿越,Δt通常为一微小量,其取值设定为完整生命周期T的1/1000。
上式中,为了求解PI{Ek}的显式表达式,需要首先构建表征G(kΔt,X,Y(kΔt),d)和G((k+1)Δt,X,Y((k+1)Δt),d)的几何可行域。这里引入标准化思想,于是有:
G(kΔt,X,Y(kΔt),d)
=Gc(kΔt,X,Y(kΔt),d)+Gr(kΔt,X,Y(kΔt),d)·U1
G((k+1)Δt,X,Y((k+1)Δt),d)
=Gc((k+1)Δt,X,Y((k+1)Δt),d)+Gr((k+1)Δt,X,Y((k+1)Δt),d)·U2
其中,U1和U2表示标准化区间变量。结合第二步中功能函数关于不同时刻的相关性定义,进而构建出对应的偏转矩形域来合理表征静、动态混合不确定性下拟建结构功能函数的变化范围(如图2所示),且有:
上图中L1和L2分别表示偏转矩形域的边长(L1≤L2)。利用坐标变换,事件,G(kΔt,X,Y(kΔt),d)>0和G((k+1)Δt,X,Y((k+1)Δt),d)≤0可以进一步表示为:
上面不等式所表述的约束边界在几何上与图2所形成的矩形域相干涉,于是可得到表征穿越事件发生的几何区域(如图3所示)。基于体积比思想,穿越可能度指标PI{Ek}可被定义为干涉区域面积与总可行域面积之比,即:
(4)遍历所有时间段内的穿越可能度PI{Ek},计算基于静、动态不确定性的拟建结构混合可靠度计算指标:
其中,分别表示整个生命周期T内拟建结构的混合可靠度和混合失效度,Pos(0)表示结构在构建之初即发生失效的可能度,求解上式即可实现拟建结构动力安全态势的有效评估。
(5)以第四步定义的拟建结构混合可靠度和常规确定性约束作为限制条件,以结构的轻量化功能作为优化目标,以结构的形状与尺寸参数作为优化的设计变量,综合考虑静、动态不确定性(如材料分散性与载荷不确定性)的综合效应,构建出如下面向拟建结构减重优化的混合可靠性设计模型:
find:d
min M(T,Xc,Yc(t),d)
s.t.gi*(t,Xc,Yc(t),d)≤0,i*=1,2,...,l1
其中,M(T,Xc,Yc(t),d)为结构的名义质量函数,表示第i*个确定性约束条件,表示第j*个混合可靠性约束条件,分别表示第j*个可靠性约束下结构初始失效的可能度和时间段[kΔt,(k+1)Δt]内结构发生穿越失效的可能度,为可靠性设计的许用值,i*和j*为计数指标,l1和l2分别表示确定性约束与可靠性约束的个数,Ωd为设计变量的可行域。上述优化模型通过嵌入智能蚁群求解算法,完成设计变量的更新并实现整个迭代历程的收敛,其核心公式表示为:
其中,τuv(T0+N)和τuv(T0)分别表示时刻T0+N和T0下蚁群间从设计点u到v进行信息交换与传递的信息素,W为种群总数,反映的是蚁群个体w从设计点u到v的信息素损失,ΨN为权重系数。
第六步:迭代过程中,如果当前设计不满足第五步混合可靠性设计模型中拟建结构混合可靠度和常规确定性约束,或者尽管满足上述限制条件,但当前设计并不是最优设计,即相较于迭代过程中的前一组设计,目标函数的相对变化百分比大于预设值ξ时,决定设计变量取值的种群重置更新,将已经完成迭代次数的值增加1,并返回第三步,否则,进行第七步;
第七步:如果相较于迭代过程中的前一组设计,目标函数的相对变化百分比小于预设值ξ时,终止计算,将得到的全局最优设计方案中的变量参数作为最终的拟建结构设计方案。
实施例:
为了更充分地了解该发明的特点及其对工程实际的适用性,本发明针对如图4所示的拟建的类X-37B高超声速飞行器机翼结构进行含静、动态混合不确定性的非概率可靠性优化设计。该机翼由铝合金蒙皮、钛合金的梁和肋以及复合材料蜂窝三明治结构组合而成,与机身相连接一侧的机翼面固支约束,高马赫下结构的气动载荷分布作用于完整蒙皮表面。
本实施例中,考虑大气压力参数patm(t)为区间过程,合金结构的材料特性参数(包括钛合金的弹性模量ETi、泊松比μTi以及铝合金的弹性模量EAl、泊松比μAl)为区间数变量,具体参数信息见表1。此外,分别针对强度准则和刚度准则构建如下动态功能函数:
其中,表示拟建机翼结构的最大节点应力,对应为最大节点位移。于是,以结构有效质量Mt为优化目标,以蒙皮厚度Ts以及梁和肋的宽度Wr和Wb作为设计变量的可靠性优化模型可表示为:
find:Ts,Wr,Wb
min Mt(Ts,Wr,Wb)
其中,分别表示静、动态混合不确定性下的强度可靠性指标和刚度可靠性指标,约束条件上述优化问题利用智能蚁群算法迭代求解直至收敛,结果见图5。从结果中可以看出,相较于初始设计,可靠性设计的减重比达到17.4%;此外,由于采用了智能寻优策略,求解效率显著提升,有效迭代步数仅为49步即实现收敛,这对于处理大型复杂拟建结构的设计问题具有很好的借鉴作用。
表1
综上所述,本发明提出了一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法。该方法通过区间数向量和区间过程向量分别表征拟建结构的静态、动态不确定性参量,利用时间离散策略和标准化方法,定义动态功能函数的数字特征和几何表达;进而,基于时变可靠性理论中的首次穿越思想,建立面向多源不确定性的混合可靠性度量指标;最后,以可靠性水平作为约束条件,以减重作为寻优目标,以拟建结构的形状和尺寸作为设计变量,以智能蚁群算法进行高效求解,完成了面向静、动态混合不确定性的拟建结构非概率优化设计过程。
本发明未详细阐述部分属于本领域技术人员的公知技术。
以上仅是本发明的具体步骤,对本发明的保护范围不构成任何限制;其可扩展应用于含多源不确定性的复合材料拟建结构可靠性设计领域,凡采用等同变换或者等效替换而形成的技术方案,均落在本发明权利保护范围之内。

Claims (6)

1.一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法,其特征在于实现步骤如下:
第一步:考虑存在于拟建结构中的多源不确定性参数,利用区间数向量X={Xi|i=1,2,...,n1}T表征静态不确定性信息,n1为区间数向量包含的维度,Xi为区间数向量中的第i个元素,静态不确定性信息指材料参数的分散性,区间数向量X的本征规律由均值向量Xc和半径向量Xr来体现,即:
X∈[Xc-Xr,Xc+Xr]
利用区间过程向量Y(t)={Yj(t)|j=1,2,...,n2}T表征动态不确定性信息,n2为区间过程向量包含的维度,Yj(t)为区间过程向量中的第j个元素,动态不确定性信息指动力载荷的波动性;区间过程向量Y(t)的边界规律由均值过程向量Yc(t)和半径过程向量Yr(t)来体现,即:
Y(t)∈[Yc(t)-Yr(t),Yc(t)+Yr(t)]
此外,区间过程向量Y(t)任意不同时刻t1和t2下的相关性规律由自相关系数向量ρY(t1,t2)表征;
第二步:利用第一步提出的n1维静态不确定性区间数向量X和n2维动态不确定性区间过程向量Y(t),构建n1+n2维混合不确定性条件下拟建结构的线性动态化功能函数:
G(t,X,Y(t),d)=f(t,d)+a(d)X+b(t,d)Y(t)
其中,d表示可设计的结构形状与尺寸向量,包括工程梁、板及三维实体结构;f(t,d)为确定性函数,a(d)={ai(d),i=1,2,...,n1}表示时不变系数向量,b(t,d)={bj(t,d),j=1,2,...,n2}表示时变系数向量;基于区间数学运算法则,利用第一步中定义的均值向量Xc和半径向量Xr,均值过程向量Yc(t)和半径过程向量Yr(t)以及自相关系数向量ρY(t1,t2)进一步得到描述功能函数特征的均值函数Gc(t,X,Y(t),d)、半径函数Gr(t,X,Y(t),d)和自相关系数函数ρG(t1,t2)的数学表达式;
第三步:基于首次穿越理论并结合第二步构建的线性动态化功能函数,利用时间离散化方法,确定如下穿越事件Ek发生的可能度指标:
PI{Ek}=Pos{G(kΔt,X,Y(kΔt),d)>0∩G((k+1)Δt,X,Y((k+1)Δt),d)≤0}其中,PI{Ek}表示事件发生的可能度,Ek表示为拟建结构在时间段[kΔt,(k+1)Δt]内发生了一次穿越,G(kΔt,X,Y(kΔt),d)>0表示拟建结构在kΔt时刻安全,即功能函数大于零,G((k+1)Δt,X,Y((k+1)Δt),d)≤0表示拟建结构在(k+1)Δt时刻失效,即功能函数小于等于零,符号“∩”表示事件的交运算,k为计数指标,Δt表示时间增量;
第四步:遍历所有时间段内的穿越可能度PI{Ek},计算基于静态不确定性、动态不确定性的拟建结构混合可靠度计算公式:
其中,分别表示整个生命周期T内拟建结构的混合可靠度和混合失效度,Pos(0)表示结构在构建之初即发生失效的可能度,求解上式即可实现拟建结构动力安全态势的有效评估;
第五步:以第四步的拟建结构混合可靠度和常规确定性约束作为限制条件,以结构的轻量化功能作为优化目标,以结构的形状与尺寸参数作为优化的设计变量,综合考虑静态不确定性、动态不确定性的综合效应,构建出面向拟建结构减重优化的混合可靠性设计模型,并通过嵌入智能蚁群迭代算法,完成设计变量的更新并实现整个优化历程的收敛;
第六步:迭代过程中,如果当前设计不满足第五步混合可靠性设计模型中拟建结构混合可靠度和常规确定性约束,或者尽管满足上述限制条件,但当前设计并不是最优设计,即相较于迭代过程中的前一组设计,目标函数的相对变化百分比大于预设值ξ时,设计变量取值的种群重置更新,将已经完成迭代次数的值增加1,并返回第三步,否则,进行第七步;
第七步:如果相较于迭代过程中的前一组设计,目标函数的相对变化百分比小于预设值ξ时,终止计算,将得到的全局最优设计方案中的变量参数作为最终的拟建结构设计方案。
2.根据权利要求1所述的一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法,其特征在于:所述第一步中区间数向量X的均值向量Xc和半径向量Xr表示为:
其中,Xi 是区间数变量Xi的下界和上界,区间过程向量Y(t)的均值过程向量Yc(t)、半径过程向量Yr(t)以及任意不同时刻t1和t2下的自相关系数向量ρY(t1,t2)表示为:
以及
其中,Yj(t)是区间过程Yj(t)的下界和上界,表示任意不同时刻t1和t2下区间过程Yj(t)的自相关系数函数,为对应的自协方差函数。
3.根据权利要求1所述的一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法,其特征在于:所述第二步中描述线性动态化功能函数特征的均值函数Gc(t,X,Y(t),d)和半径函数Gr(t,X,Y(t),d)分别表示为:
此外,自相关系数函数ρG(t1,t2)的数学表达式为:
其中,CovG(t1,t2)表示功能函数G(t,X,Y(t),d)的协方差函数。
4.根据权利要求1所述的一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法,其特征在于:所述第三步中Δt为微小量,取值设定为完整生命周期T的1/1000。
5.根据权利要求1所述的一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法,其特征在于:所述第五步中构建出面向拟建结构减重优化的混合可靠性设计模型描述为:
find:d
min M(T,Xc,Yc(t),d)
s.t.gi*(t,Xc,Yc(t),d)≤0,i*=1,2,...,l1
j*=1,2,...,l2
d∈Ωd
其中,M(T,Xc,Yc(t),d)为结构的名义质量函数,表示第i*个确定性约束条件,表示第j*个混合可靠性约束条件,分别表示第j*个可靠性约束下结构初始失效的可能度和时间段[kΔt,(k+1)Δt]内结构发生穿越失效的可能度,为可靠性设计的许用值,i*和j*为计数指标,l1和l2分别表示确定性约束与可靠性约束的个数,Ωd为设计变量的可行域。
6.根据权利要求1所述的一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法,其特征在于:所述第六步中目标函数的相对变化百分比的预设值ξ设定为1%。
CN201610407978.9A 2016-06-12 2016-06-12 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法 Expired - Fee Related CN105956336B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610407978.9A CN105956336B (zh) 2016-06-12 2016-06-12 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610407978.9A CN105956336B (zh) 2016-06-12 2016-06-12 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法

Publications (2)

Publication Number Publication Date
CN105956336A CN105956336A (zh) 2016-09-21
CN105956336B true CN105956336B (zh) 2019-01-18

Family

ID=56908802

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610407978.9A Expired - Fee Related CN105956336B (zh) 2016-06-12 2016-06-12 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法

Country Status (1)

Country Link
CN (1) CN105956336B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106407586B (zh) * 2016-09-28 2019-08-13 北京航空航天大学 一种基于不确定优化的飞行器结构功能一体化设计平台系统
CN106650141B (zh) * 2016-12-29 2018-12-21 北京航空航天大学 一种预测周期性材料性能的不确定性分析方法
CN106709215B (zh) * 2017-02-21 2019-02-22 北京航空航天大学 一种基于级数展开的连续体结构非概率可靠性拓扑优化方法
CN107451375A (zh) * 2017-08-18 2017-12-08 电子科技大学 考虑混合不确定性的结构可靠性优化设计方法
CN107729647B (zh) * 2017-10-13 2020-09-25 电子科技大学 一种将时变可靠性转换成静态可靠性的分析方法
CN107958310B (zh) * 2017-12-07 2019-03-19 北京航空航天大学 一种基于区间时变可靠性的在役结构最优维护设计方法
CN108710713B (zh) * 2018-02-07 2021-07-23 北京航空航天大学 一种求解含区间参数结构的静力学响应的贝叶斯配点方法
CN108776433B (zh) * 2018-07-20 2021-02-02 北京航空航天大学 一种静态与动态信息混合的容错控制方法
CN111967140B (zh) * 2020-07-20 2023-07-07 北京航空航天大学 一种考虑混合不确定性的性能退化实验建模与分析方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103324782A (zh) * 2013-05-29 2013-09-25 北京航空航天大学 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法
CN105653794A (zh) * 2015-12-30 2016-06-08 北京航空航天大学 一种含初始缺陷矩形板结构的时变可靠性设计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103324782A (zh) * 2013-05-29 2013-09-25 北京航空航天大学 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法
CN105653794A (zh) * 2015-12-30 2016-06-08 北京航空航天大学 一种含初始缺陷矩形板结构的时变可靠性设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Hybrid reliability analysis of structures with multi-source uncertainties;Lei Wang等;《Acta Mechanica》;20140228;第413-430页
基于非概率凸模型可靠性的结构优化设计;王晓军等;《北京航空航天大学学报》;20120531;第630-635页

Also Published As

Publication number Publication date
CN105956336A (zh) 2016-09-21

Similar Documents

Publication Publication Date Title
CN105956336B (zh) 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法
CN103366065B (zh) 一种基于区间可靠性的飞行器热防护系统尺寸优化设计方法
Zhao et al. Extended VIKOR method based on cross-entropy for interval-valued intuitionistic fuzzy multiple criteria group decision making
CN105930562A (zh) 一种非概率条件下的结构性能优化设计方法
CN105144206B (zh) 多目标优化方法和设备
CN102999678A (zh) 一种基于非线性多目标区间鲁棒优化的汽车降噪方法
CN107958310B (zh) 一种基于区间时变可靠性的在役结构最优维护设计方法
CN105893716A (zh) 一种基于分形理论的结构断裂非概率可靠性分析方法
CN105939026B (zh) 基于混合Laplace分布的风电功率波动量概率分布模型建立方法
Lin et al. Uncertainty quantification in state estimation using the probabilistic collocation method
Shao et al. A comprehensive G-Lz-ADC effectiveness evaluation model for the single communication satellite system in the context of poor information
Zhang et al. Spherical fuzzy hamacher power aggregation operators based on entropy for multiple attribute group decision making
CN109446557A (zh) 一种基于概率密度演化的随机气动弹性系统稳定性分析方法
CN103617337B (zh) 一种飞机维修间隔辅助分析方法
CN106484957B (zh) 一种重复使用运载器再入飞行制导控制性能评估系统
Yang et al. MAGDM in hesitant interval-valued Pythagorean linguistic Z-number based on combined score function and entropy
Nikbay et al. Flutter based aeroelastic optimization of an aircraft wing with analytical approach
Wijayasekara et al. Fuzzy linguistic knowledge based behavior extraction for building energy management systems
CN107656901B (zh) 一种配电网严格安全边界和非严格安全边界的求解方法
Gao et al. Distributed learning control for high-speed trains subject to operation safety constraints
Tang et al. Nonprobabilistic reliability analysis for an inside flap of an aircraft
Xing et al. Balanced Evaluation for Flight Ground Support Process Based on Approximate Ideal Point Sort Method
Nielsen Design of Aerospace Laminates for Multi-Axis Loading and Damage Tolerance
Wang et al. Application of gray fuzzy analytic hierarchy process in air traffic safety of tower control
Wang et al. Fast calculation of probabilistic TTC with static voltage stability constraint

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190118