非线性功能函数串联桁架结构体系非概率可靠性分析方法

1.本发明属于串联桁架结构体系可靠性分析技术领域，具体涉及非线性功能函数串联桁架结构体系非概率可靠性分析方法。


背景技术：

2.桁架结构以其受力简单、重量轻、便于制造和组装等优点广泛应用与工程结构和机械机构，如桥梁、厂房、塔吊、桁架梁和桁架臂等结构中。受材料属性、制造条件、服役环境、运行工况等因素的影响，机构中往往存在着与弹性模量、强度、刚度、尺寸和载荷等参数相关的不确定参数。结构中的不确定参数的存在会使结构在服役过程中存在一定的安全隐患，为了确保结构的安全可靠，须对结构进行可靠性分析。传统可靠性分析采用随机模型度量结构中的不确定参数，进而发展了概率可靠性分析方法。然而构建随机模型需大量的不确定参数样本信息，而对于试验难度大和试验成本高的结构，结构中不确定参数的样本信息通常不足。鉴于此，只需少量样本信息确定不确定参数变化区间而不需明确其具体分布的凸集模型被用来度量结构中的不确定参数，并发展了相应的非概率可靠性分析方法。现有的串联结构体系可靠性分析的解析方法，如非概率区间估计法和非概率逐步等效平面法等，均是针对线性功能函数结构体系，针对非线性结构功能函数的结构体系解析可靠性分析方法尚不足。考虑到实际结构中结构功能函数大多是关于不确定参数的非线性函数，因此，为了兼顾结构体系可靠性分析过程中的精度与效率，发展针对非线性结构功能函数的结构体系解析可靠性分析方法成为结构可靠性分析中亟需解决的问题。


技术实现要素：

3.本发明针对现有的用于计算非线性功能函数结构体系非概率可靠度或失效度的方法不足的问题，提出非线性功能函数串联桁架结构体系非概率可靠性分析方法，在采用凸集模型度量结构中不确定参数的背景下，基于o.ditlevsen界限法的基本思想，通过在单个失效模式的设计点对其进行泰勒二阶展开并求解一阶失效域体积和在任意两个失效模式的联合验算点对两个失效模式进行泰勒一阶展开并求解两个失效模式共失效域体积，发展用于计算非线性功能函数串联结构体系非概率失效度的区间估计法，通过该方法计算非线性功能函数串联结构体系的失效度不仅能够提高可靠性分析的精度，还能解决数值模拟方法效率低的问题，其应用前景广泛且便于推广使用。
4.为了解决上述技术问题，本发明通过以下技术方案予以实现：
5.非线性功能函数串联桁架结构体系非概率可靠性分析方法，包括以下步骤：
6.步骤1、构建用于度量桁结构中不确定参数的椭球模型；
7.步骤2、获取具有单位半径等价球体模型并计算球体体积；
8.步骤3、获取桁架结构体系失效模式的结构功能函数并转换至标准参数空间；
9.步骤4：获取结构功能函数的二阶近似功能函数并计算单个结构功能函数对应的一阶失效域体积；
10.步骤5、在两两功能函数的联合验算点进行一阶泰勒展开并计算两个失效模式共失效域体积；
11.步骤6、采用o.ditlevsen界限法计算串联结构体系失效域体积并求解最短区间；
12.步骤7、计算串联桁架结构体系的非概率失效度。
13.进一步地，步骤1中，所述构建用于度量桁结构中不确定参数的椭球模型具体为：
14.根据结构中不确定参数的数字特征构建度量不确定参数的椭球模型，如下：
[0015][0016]
其中，x为不确定参数向量且x＝(x1,x2,
…
,xn)
t
，n为不确定参数向量x的维数，x中第i个元素xi的边缘区间为为不确定参数xi的下界，为不确定参数xi的上界，为不确定参数的均值向量且其第i个元素为ω
x
为用于确定椭球模型尺寸和方向的特征矩阵且ρ
ij
为第i个不确定参数xi和第j个不确定参数xj的相关系数，j为正整数且取值范围为1～n，当i＝j时，ρ
ij
＝1，为第i个不确定参数边缘区间的区间半径且rn为n维实数域。
[0017]
进一步地，步骤2中，所述获取具有单位半径等价球体模型并计算球体体积，具体包括：
[0018]
步骤201、不确定参数向量的归一化处理：
[0019]
根据公式获取不确定参数向量x的归一化不确定参数向量u，其中，u＝(u1,u2,
…
,un)
t
，ui为第i个不确定参数xi对应的归一化不确定参数；
[0020]
步骤202、构建不确定参数的多维归一化等价椭球模型：
[0021]
根据步骤201中对不确定参数的归一化处理获取归一化等价椭球模型其中ωu为椭球模型ω
x
在归一化空间u中的等价特征矩阵且ωu＝diag(xr)ω
x
diag(xr)，diag(xr)为以xr中元素为对角元素的n维对角矩阵；
[0022]
步骤203、对归一化等价椭球模型的特征矩阵ωu进行choleskey分解，即其中，l0为choleskey分解得到的下三角矩阵；
[0023]
步骤204、将归一化等价椭球模型转化到标准参数空间获取具有单位半径的等价球体模型e
δ
＝{δ|δ
t
δ≤1,δ∈rn}，其中，
为标准参数空间中的不确定参数，δi为对应于不确定参数xi的标准不确定参数；
[0024]
步骤205、计算单位半径等价球体模型的体积：
[0025]
根据公式计算多维等价单位球体模型e
δ
的体积vn，其中，γ(
·
)为gamma函数。
[0026]
进一步地，步骤3中，所述获取桁架结构体系失效模式的结构功能函数并转换至标准参数空间，具体包括：
[0027]
步骤301、确定桁架结构体系各失效模式的功能函数：
[0028]
根据桁架结构失效准则确定桁架结构体系各失效模式的功能函数g
l
(x)，其中l为桁架结构体系失效模式的编号且l＝1,2,
…
,m，m为桁架结构体系失效模式总个数；
[0029]
步骤302、获取标准参数空间中的结构功能函数：
[0030]
按照不确定参数向量x与标准参数空间δ空间中的标准化向量δ之间的关系对桁架结构体系各失效模式功能函数g
l
(x)进行变量代换，得到标准参数空间δ空间中相应失效模式的结构功能函数g
l
(δ)。
[0031]
进一步地，步骤4中，所述获取结构功能函数的二阶近似功能函数并计算单个结构功能函数对应的一阶失效域体积，具体包括：
[0032]
步骤401、判断串联桁架结构体系各失效模式的功能函数g
l
(δ)是否为不确定参数向量δ的线性函数：
[0033]
对功能函数g
l
(δ)进行求导处理，当功能函数g
l
(δ)关于不确定参数向量δ的一阶偏导数均为常数，表明功能函数g
l
(δ)是关于不确定参数向量δ的线性函数，执行步骤402；当功能函数g
l
(δ)关于不确定参数向量δ的一阶偏导数存在非常数，表明功能函数g
l
(δ)是关于不确定参数向量δ的非线性函数，执行步骤403；
[0034]
步骤402、功能函数g
l
(δ)是不确定参数向量δ的线性函数，g
l
(δ)可记为：其中，a
l0
表示功能函数g
l
(δ)的常数项，a
l
表示功能函数g
l
(δ)中不确定参数向量δ的系数向量且a
l
＝(a
l1
,a
l2
,
…
,a
ln
)
t
，假设该结构功能函数的非概率可靠性指标为β
l
，则相应的一阶失效体积可表示为
[0035]
步骤403、结构功能函数g
l
(δ)是不确定参数向量δ的非线性函数，根据优化问题
求得标准参数空间中第l个结构功能函数设计点并在该设计点对结构功能函数进行泰勒二阶展开，得第l个结构功能函数的二阶近似结构功能函数为其中，表示结构功能函数g
l
(δ)在设计点处的梯度向量且(δ)在设计点处的梯度向量且表示海森矩阵且考虑到设计点在失效面上，因此有同时令α
lδ
表示与结构功能函数g
l
(δ)在设计点处的梯度向量同方向的单位向量，且令则第l个结构功能函数的非概率可靠性指标可表示为令则可得通过正交变换δ＝h
l
y，其中，h
l
＝[h
l1
,h
l2
,
…
,h
l(n-1)
,α
lδ
]，可通过施密特正交化构建，可得y变量空间的极限状态方程为g
ly
(y)＝g
ly
(y1,y2,
…
,yn)，通过上述变换即可将失效面近似成一个以的方向为对称轴的抛物面，其功能函数可表示为其中是一个n-1维实对称方阵且有n-1个特征向量及相应的特征值κ
li
,i＝1,2,
…
,n-1，然后根据公式求得近似曲面在设计点处的平均曲率，其中其中分别为第l个结构功能函数对应的矩阵m
l
和q
l
的对角元素，可得
考虑到y空间中椭球域可表示为y1+y2+
…
+yn＝1，将代入y1+y2+
…
+yn＝1可得可得该区域的体积为最后单个失效模式的二阶失效体积可表示为
[0036]
进一步地，步骤5中，所述在两两功能函数的联合验算点进行一阶泰勒展开并计算两个失效模式共失效域体积，具体包括：
[0037]
取标准参数空间δ空间中的m个结构功能函数中的任意两个记为g
l
(δ)和gk(δ)，通过求解优化问题求得两个结构功能函数的联合验算点并在该验算点对两个结构功能函数进行泰勒一阶展开，可得两个结构功能函数的非概率可靠性指标和相关系数可分别表示为结构功能函数的非概率可靠性指标和相关系数可分别表示为可得两个失效模式的共失效域体积
其中γ
lk
＝arccos(ρ
lk
)。
[0038]
进一步地，步骤6中，所述采用o.ditlevsen界限法计算串联结构体系失效域体积并求解最短区间，具体包括：
[0039]
步骤601、串联结构体系非概率可靠性分析的o.ditlevsen界限法：
[0040]
根据o.ditlevsen界限法的基本思想，分别采用一阶失效域体积和二阶失效域体积代替o.ditlevsen界限法中的一阶失效概率和二阶失效概率，即可得到计算串联结构体系失效域体积的非概率窄界限法公式然后针对非概率窄界限法的上下界表达式寻找合理失效模式排序方式从而确定窄界限法中的最大下界和最小上界；
[0041]
步骤602、o.ditlevsen界限法的最大下界：
[0042]
非概率窄界限法中下界如果对任意j都满足则该下界可表示为此时失效模式的排序对计算串联结构体系失效域体积的下界没有影响；而当非概率窄界限法下界中存在j使得这时失效模式的排列顺序对下界的求解就会产生影响，根据图论理论，最大下界的求解在于再相关权重图中选择一些点使得最大，具体步骤为：(1)选择最大权值点并设置为序号a；(2)在其余点中选择满足max(v
x-v
ax
)的x点并将其设置为序号b；(3)再在其余点中选择满足max(v
x-v
ax-v
bx
)的x点并将其设置为序号c；(4)重复上述过程直到所有剩余点满足可得o.ditlevsen界限法最大下界
[0043]
步骤603、o.ditlevsen界限法的最小上界：非概率窄界限法的上界为
求最小上界实际是求的最大值根据图论理论，可通过寻找最大树的运算法则求解，具体步骤为：(1)将两个失效模式共失效域体积大小作为联结图中线的权重值，假设联结图中第l
t
条线对应的权重值为v(l
t
)，v(l
t
)与v
jk
一一对应；(2)取权值最大的两条线作为最大树的分支；(3)选择另一条线，使之不与已选择的线构成封闭路且权重值为除已选中线外的最大者；(4)重复步骤(3)的过程直到联结图中的所有点都被以上所选择的线连接起来，即可得到最大树，相应线条权重值之和即为将其代入公式即可得最小上界
[0044]
进一步地，所述不确定参数向量x包括惯性矩参数向量、荷载参数向量和弹性模量参数向量。
[0045]
进一步地，不确定参数向量维数n满足n≥2。
[0046]
进一步地，当不确定参数向量维数n＝2时，步骤205、步骤4和步骤5所求为面积；当n＝3时，步骤205、步骤4和步骤5所求为体积；当n＞3时，步骤205、步骤4和步骤5所求为多维体积；
[0047]
步骤402和步骤5所述一阶失效体积公式和两个失效模式共失效域体积公式中β取值范围为(0,1)。
[0048]
与现有技术相比，本发明至少具有以下有益效果：
[0049]
1.可用于非线性功函数串联结构体系的非概率可靠性分析。弥补了现有方法竟能求解线性功能函数的串联结构体系可靠度的不足。
[0050]
2.提出在来联合验算点对两个结构功能函数进行泰勒一阶展开，解决了两个非线性功能函数的失效模式共失效域体积难以求解的问题。
[0051]
3.采用基于图论理论的排序方法求解o.ditlevsen界限法的最大下界和最小上界，该方法精度高、效果好。
[0052]
4.所采用的o.ditlevsen界限法具有足够的求解精度和较之于数值模拟法更高的求解效率。
[0053]
综上所述，本发明通过采用o.ditlevsen界限法计算非线性功能函数串联结构体系的非概率失效度，使求解复杂的非线性功能函数串联结构体系可靠度问题成为可能。该方法极大简化了桁架结构非概率可靠性分析的复杂性。在满足工程实际问题精度要求的前提下实现了高效求解，其适用面广且应用前景广泛。本发明通过在设计点对非线性结构功能函数进行二阶近似，提高了对单个失效模式失效域体积的近似精度，然后在对两个失效模式影响最大的联合验算点进行一阶近似从而求解两个失效模式共失效域体积，在提高串联桁架结构体系非概率可靠度求解精度的同时还能保持较高的求解效率。
[0054]
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
[0055]
为了更清楚地说明本发明具体实施方式中的技术方案，下面将对具体实施方式描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0056]
图1为本发明的方法流程框图。
[0057]
图2为本实施例中平面串联二十九杆静定桁架结构体系。
[0058]
图3为本实施例中以失效域体积作为权重值的权重图。
具体实施方式
[0059]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0060]
如图1和图2所示，本发明的非线性功能函数串联桁架结构体系非概率可靠性分析方法，包括以下步骤：
[0061]
步骤1、构建用于度量结构中不确定参数的椭球模型：根据结构中不确定参数的数字特征构建度量不确定参数的椭球模型为其中，x为不确定参数向量且x＝(x1,x2,
…
,xn)
t
，n为不确定参数向量x的维数，x中第i个元素xi的边缘区间为为不确定参数xi的下界，为不确定参数xi的上界，为不确定参数的均值向量且其第i个元素为ω
x
为用于确定椭球模型尺寸和方向的特征矩阵且ρ
ij
为第i个不确定参数xi和第j个不确定参数xj的相关系数，j为正整数且取值范围为1～n，当i＝j时，ρ
ij
＝1，x
ir
为第i个不确定参数边缘区间的区间半径且rn为n维实数域；
[0062]
本实施例中，以图2所示的平面二十九杆静定桁架结构体系为例，在平面二十九杆静定桁架结构体系中，按杆件截面积大小将所有杆件分为四类，即第1类杆件：上弦杆和侧杆，第2类杆件：斜腹杆，第3类杆件：竖直腹杆，第4类杆件：下弦杆。结构系统中的不确定参数包括对应于四类杆件的惯性矩x1,x2,x3,x4，作用在结构体系上的荷载x5和结构体系中材料的弹性模量x6，即结构中不确定参数向量x的维数为n＝6，相应的不确定参数向量为x＝
(x1,x2,x3,x4,x5,x6)
t
。其中，第1类杆件惯性矩x1的边缘区间为第2类杆件惯性矩x2的边缘区间为第3类杆件惯性矩x3的边缘区间为第4类杆件惯性矩x4的边缘区间为作用在结构体系上的荷载x5的边缘区间为材料的弹性模量x6的边缘区间为四个惯性矩之间的相关系数为ρ
12
＝ρ
13
＝ρ
14
＝ρ
23
＝ρ
24
＝ρ
34
＝0.7，惯性矩、载荷和材料弹性模量之间的相关系数为0。不确定参数的均值为边缘区间的区间半径为其中为第i个不确定参数的上界，为第i个不确定参数的下界，故其边缘区间的均值向量为xc＝(6774100mm4,966800mm4,1554300mm4,3754200mm4,68.2kn,210gpa)
t
，边缘区间的区间半径向量为xr＝(203223mm4,29004mm4,46629mm4,112626mm4,2.046kn,0.63gpa)
t
，用于确定椭球模型方向和形状的特征矩阵为则用于度量该桁架结构中不确定参数的椭球模型为
[0063]
步骤2、获取具有单位半径等价球体模型并计算球体体积，过程如下：
[0064]
步骤201、不确定参数向量的归一化处理：根据公式获取不确定参数向量x的归一化不确定参数向量u，其中，u＝(u1,u2,
…
,un)
t
，ui为第i个不确定参数xi对应的归一化不确定参数；
[0065]
步骤202、构建不确定参数的多维归一化等价椭球模型：根据步骤201中对不确定参数的归一化处理获取归一化等价椭球模型其中ωu为椭球模型ω
x
在归一化空间u中的等价特征矩阵且ωu＝diag(xr)ω
x
diag(xr)，diag(xr)为以xr中元素为对角元素的n维对角矩阵；
[0066]
本实施例中，不确定参数向量x的归一化不确定参数向量u为
根据公式ωu＝diag(xr)ω
x
diag(xr)可得归一化空间中等价椭球模型的特征矩阵则可得归一化空间中等价椭球模型
[0067]
步骤203、对归一化等价椭球模型的特征矩阵ωu进行choleskey分解，即其中，l0为choleskey分解得到的下三角矩阵；
[0068]
步骤204、将归一化等价椭球模型转化到标准参数空间获取具有单位半径的等价球体模型e
δ
＝{δ|δ
t
δ≤1,δ∈rn}，其中，为标准参数空间中的不确定参数，δi为对应于不确定参数xi的标准不确定参数；
[0069]
对归一化空间中等价椭球模型的特征矩阵进行choleskey分解，可得下三角矩阵为根据公式即可将原空间中的不确定参数向量转化到标准参数空间δ空间，得具有单位半径的等价球体模型e
δ
＝{δ|δ
t
δ≤1,δ∈rn}；
[0070]
步骤205、计算单位半径等价球体模型的体积：根据公式计算多维等价单位球体模型e
δ
的体积vn，其中，γ(
·
)为gamma函数；
[0071]
本实施例中，不确定参数维数为n＝6，根据公式可得单位半径等
价球体模型的体积为
[0072]
步骤3、获取桁架结构体系失效模式的结构功能函数并转换至标准参数空间，过程如下：
[0073]
步骤301、确定串联桁架结构体系各失效模式的功能函数：根据桁架结构失效准则确定串联桁架结构体系各失效模式的功能函数g
l
(x)，其中，l为结构体系失效模式的编号且l＝1,2,
…
,m，m为结构体系失效模式的总个数；
[0074]
本实施例中，根据受压杆件满足压杆稳定性条件和桁架下弦中点最大竖直位移不超过13mm确定了该二十九杆静定桁架结构的七个主要失效模式，即失效模式总个数m＝7，七个失效模式对应的结构功能函数为：
[0075][0076][0077][0078][0079][0080][0081][0082]
步骤302、获取标准参数空间中的结构功能函数：按照不确定参数向量x与标准参数空间δ空间中的标准化向量δ之间的关系对串联桁架结构体系各失效模式的功能函数g
l
(x)进行变量代换，得到标准参数空间δ空间中相应失效模式的结构功能函数g
l
(δ)；
[0083]
本实施例中，根据不确定参数向量x与标准参数空间δ空间中的标准化向量δ之间的关系：可得
将其代入原变量空间的七个失效模式对应的结构功能函数中，得标准参数空间δ空间中相应失效模式的结构功能函数g
l
(δ)：
[0084]
g1(δ)＝2.1622e-8(6.2179δ6+2072.6)(126500δ1+38574δ2+59759δ3+142266δ4+6.7741e6)-8.8546δ
5-295.15
[0085]
g2(δ)＝6.25e-8(6.2179δ6+2072.6)*(126500δ1+38574δ2+59759δ3+142266δ4+6.7741e6)-8.928δ
5-297.6
[0086]
g3(δ)＝3.3058e-8(6.2179δ6+2072.6)(33299δ3+32640δ4+1.5543e6)-3.069δ
5-102.3
[0087]
g4(δ)＝6.25e-8(6.2179δ6+2072.6)(126500δ1+38574δ2+59759δ3+142266δ4+6.7741e6)-11.16δ
5-372；
[0088]
g5(δ)＝2.1622e-8(6.2179δ6+2072.6)(18876δ2+8528.9δ3+20303δ4+966800)-1.2649δ
5-42.165
[0089]
g6(δ)＝6.25e-8(6.2179δ6+2072.6)(126500δ1+38574δ2+59759δ3+142266δ4+6.7741e6)-11.16δ
5-372
[0090][0091]
步骤4、获取结构功能函数的二阶近似功能函数并计算单个结构功能函数对应的一阶失效域体积，过程如下：
[0092]
步骤401、判断串联桁架结构体系各失效模式的功能函数g
l
(δ)是否为不确定参数向量δ的线性函数：对功能函数g
l
(δ)进行求导处理，当功能函数g
l
(δ)关于不确定参数向量δ的一阶偏导数均为常数，表明功能函数g
l
(δ)是关于不确定参数向量δ的线性函数，执行步骤402；当功能函数g
l
(δ)关于不确定参数向量δ的一阶偏导数存在非常数，表明功能函数g
l
(δ)是关于不确定参数向量δ的非线性函数，执行步骤403；
[0093]
步骤402、功能函数g
l
(δ)是不确定参数向量δ的线性函数，g
l
(δ)可记为：其中，a
l0
表示功能函数g
l
(δ)的常数项，a
l
表示功能函数g
l
(δ)中不确定参数向量δ的系数向量且a
l
＝(a
l1
,a
l2
,
…
,a
ln
)
t
，假设该结构功能函数的非概率可靠性指标为β
l
，则相应的一阶失效体积可表示为：
[0094][0095]
步骤403、结构功能函数g
l
(δ)是不确定参数向量δ的非线性函数。根据优化问题求得标准参数空间中第l个结构功能函数设计点并在该设计点对结构功能函数进行泰勒二阶展开，得第l个结构功能函数的二阶近似结构功能函数为其中表示结构功能函数g
l
(δ)在设计点处的梯度向量且(δ)在设计点处的梯度向量且表示海森矩阵且考虑到设计点在失效面上，因此有同时令α
lδ
表示与结构功能函数g
l
(δ)在设计点处的梯度向量同方向的单位向量，且令则第l个结构功能函数的非概率可靠性指标可表示为令则可得通过正交变换δ＝h
l
y，其中且可通过施密特正交化构建，可得y变量空间的极限状态方程为g
ly
(y)＝g
ly
(y1,y2,
…
,yn)，通过上述变换即可将失效面近似成一个以的方向为对称轴的抛物面，其功能函数可表示为
其中是一个n-1维实对称方阵且有n-1个特征向量及相应的特征值κ
li
,i＝1,2,
…
,n-1，然后根据公式求得近似曲面在设计点处得平均曲率，其中其中分别为第l个结构功能函数对应的矩阵m
l
和q
l
的对角元素，可得考虑到y空间中椭球域可表示为y1+y2+
…
+yn＝1，将代入y1+y2+
…
+yn＝1可得可得该区域的多维体积为最后单个失效模式的二阶失效面积可表示为
[0096]
本实施例中，通过求导可知，七个结构功能函数对标准参数空间中的不确定参数向量δ的一阶偏导均存在非常数，表明七个结构功能函数均是关于不确定参数的非线性函数，故执行步骤403，七个结构功能函数的非概率可靠性指标和设计点分别为：
[0097][0098][0099]
根据表中数据可知失效模式g2(δ)、g4(δ)和g6(δ)的非概率可靠性指标分别为20.9222、17.7184和17.7184，非概率可靠性指标大于1表明该失效模式的失效域与凸集模型无交集，既该失效模式为无效的，因此，以下仅对非概率可靠性指标小于1的四个失效模式g1(δ)、g3(δ)、g5(δ)和g7(δ)进行分析，根据步骤403计算得四个结构功能函数对应的单个失效模式的二阶失效面积a1为：
[0100] g1(δ)g3(δ)g5(δ)g7(δ)v
cap
0.1358590.0003303420.1658340.106155v
pan
0.000015360.00000011810.000029780.0000003188凹凸性凸凸凹凸v10.135874360.00033046010.165804220.1061553188
[0101]
步骤5、在两两功能函数的联合验算点进行一阶泰勒展开并计算两个失效模式共失效域体积，过程如下：
[0102]
取标准参数空间δ空间中的m个结构功能函数中的任意两个记为g
l
(δ)和gk(δ)。通过求解优化问题求得两个结构功能函数的联合验算点并在该验算点对两个结构功能函数进行泰勒一阶展开，可得两个结构功能函数的非概率可靠性指标和相关系数可分别表示为结构功能函数的非概率可靠性指标和相关系数可分别表示为可得两个失效模式的共失效域体积其中γ
lk
＝arccos(ρ
lk
)；
[0103]
本实施例中，通过优化问题求得的四个结构功能函数g1(δ)、g3(δ)、g5(δ)和g7(δ)两两之间的联合验算点为：
[0104][0105]
在联合验算点对结构功能函数进行泰勒一阶展开，然后求得两个失效模式的共失效域体积为：
[0106][0107]
步骤6、采用o.ditlevsen界限法计算串联结构体系失效域体积并求解最短区间，过程如下：
[0108]
步骤601、串联结构体系非概率可靠性分析的o.ditlevsen界限法：根据o.ditlevsen界限法的基本思想，分别采用一阶失效域体积和二阶失效域体积代替o.ditlevsen界限法中的一阶失效概率和二阶失效概率，即可得到计算串联结构体系失效域体积的非概率窄界限法公式然后针对非概率窄界限法的上下界表达式寻找合理失效模式排序方式从而确定窄界限法中的最大下界和最小上界；
[0109]
步骤602、o.ditlevsen界限法的最大下界：非概率窄界限法中下界如果对任意j都满足则该下界可表示为此时失效模式的排序对计算串联结构体系失效域体积的下界没有影响。而当非概率窄界限法下界中存在j使得这时失效模式的排列顺序对下界的求解就会产生影响，根据图论理论，最大下界的求解在于再相关权重图中选择一些点使得最大，具体步骤为：(1)选择最大权值点并设置为序号a；(2)在其余点中选择满足max(v
x-v
ax
)的x点并将其设置为序号b；(3)再在其余点中选择满足max(v
x-v
ax-v
bx
)的x点并将其设置为序号c；(4)重复上述过程直到所有剩余点满
足可得o.ditlevsen界限法最大下界
[0110]
本实施例中，根据步骤4中求得的单个失效模式的失效域体积和步骤5中求得的两个失效模式的共失效域体积，绘制出以失效域体积作为权重值的权重图如图3所示，根据步骤602求得o.ditlevsen界限法最大下界对应失效模式排序为g5(δ)
→
g7(δ)
→
g1(δ)
→
g3(δ)，对应最大下界为
[0111]
步骤603、o.ditlevsen界限法的最小上界：非概率窄界限法的上界为求最小上界实际是求的最大值根据图论理论，可通过寻找最大树的运算法则求解，具体步骤为：(1)将两个失效模式共失效域体积大小作为联结图中线的权重值，假设联结图中第l
t
条线对应的权重值为v(l
t
)，v(l
t
)与v
jk
一一对应；(2)取权值最大的两条线作为最大树的分支；(3)选择另一条线，使之不与已选择的线构成封闭路且权重值为除已选中线外的最大者；(4)重复步骤(3)的过程直到联结图中的所有点都被以上所选择的线连接起来，即可得到最大树，相应线条权重值之和即为将其代入公式即可得最小上界
[0112]
本实施例中，基于步骤602中的权重图，然后根据步骤603求得o.ditlevsen界限法最小上界对应失效模式排序为g1(δ)
→
g5(δ)
→
g3(δ)
→
g7(δ)，对应最大下界为
[0113]
步骤7、计算串联桁架结构体系的非概率失效度：考虑到结构体系的非概率可靠度rs和非概率失效度fs满足rs+fs＝1，为了方便起见，本发明采用非概率失效度作为结构体系可靠性的度量，根据非概率失效度的定义，即非概率失效度定义为凸集模型陷入结构体系失效域部分的体积与整个椭球模型体积之比，即串联桁架结构体系的非概率失效度区间为
[0114]
本实施例中，根据步骤602获得的串联结构体系失效域体积的最大下界步骤603获得的串联结构系体失效域体积的最小上界以及步骤步骤205获得的单位半径等价球体模型的体积即可通过给出该串联桁架结构体系的非概率失效度区间为[0.0546884,0.0599835]。
[0115]
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明
的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。