基于IDEPSO
转向架构架作为轨道车辆结构设计中的大型复杂部件，通常依靠确定性的有限元分析预判结构在设计阶段是否满足标准要求。然而，无论是静强度分析，还是疲劳强度分析，分析结果均较为理想化，与试验数据偏差较大。其原因在于分析中采用的确定性有限元模型忽略了构架在设计与制造过程中的不确定性因素，例如：材料属性、设计公差、加工偏差和试验载荷等[1]。因此，将可靠性理论与有限元分析相结合，利用可靠度评判构架在设计阶段的静态性能，可以更好地反映设计的合理性。
今年前3季度我国的中高密度纤维板进口完成14万t，比上年同期增长7.65%。主要进口来源国有：新西兰、德国、比利时等。进口金额达到1.07亿美元，比上年同期增长8.86%。
针对结构可靠性分析方法如何更好地应用于工程实际，部分学者进行了深入研究。卢耀辉等[2]考虑了影响构架最大主应力的随机变量，采用概率设计方法建立参数化模型，采用Monte Carlo 模拟（MCS）方法计算转向架构架的可靠度及灵敏度。滑林等[3]分别从区间变量随机化和随机变量区间化的角度，对船体进行可靠性分析，提高了可靠性分析的准确性和实用性。杨瑞刚等[4]鉴于起重机金属结构不确定性设计变量的随机性、模糊性和非概率性，提出一种多源不确定性混合可靠性分析方法，使可靠度结果更加符合工程实际。黄洪钟等[5]为解决现有方法进行数控机床不确定性问题可靠性分析时的局限性，基于不精确概率理论提出一种混合不确定性量化方法。此外，为解决大型复杂结构存在的物理模型难以获取、可靠度求解难度大等问题，姜潮等[6]在考虑不确定性因素的基础上，采用径向基函数代理模型表征工程机械臂参数与其固有频率之间的函数关系，采用概率-区间混合可靠性分析方法对失效概率进行估计；Zhi 等[7]建立转向架构架的多级响应面代理模型，运用Monte Carlo 模拟计算转向架构架的可靠度。为提高复杂结构可靠性分析的精度，Echard 等[8]提出一种结合Kriging模型和MCS的主动学习可靠性分析方法（AK-MCS）；Huang 等[9]采用基于Kriging模型和子集模拟（Subset Simulation，SS）的主动学习小失效概率评估方法（AK-SS），对某盾构隧道的可靠性进行分析。然而，上述研究仅能反映结构在单工况下的可靠性水平，未考虑结构的多工况和多随机响应特性。
Li等[10]提出改进的SS方法，解决了实际工程中结构多随机响应的可靠性问题。该方法虽然能够对多随机响应的失效概率进行快速求解，但是计算结果受每层样本个数和初始条件概率等参数随机性的影响，并且仅适用于随机变量相同的性能函数。Xiao等[11]提出一种基于主动Kriging模型的混合变量结构可靠性分析方法，解决了复杂结构中子结构的失效概率问题，但是未考虑各工况之间的相关性，并且不便于多个独立随机变量结构的工程应用。
本文以某型转向架构架为研究对象，采用改进差分进化粒子群算法（IDEPSO）获取最优的每层样本个数和初始条件概率，以提高SS 方法的求解精度，进而采用几何分析确定联合失效概率，计算多工况下构架的可靠度及最优失效次序，并与传统方法进行对比，验证基于IDEPSO-SS可靠性分析的准确性，为设计阶段转向架构架在多种组合工况下的总体性能分析方法由传统的确定性分析向不确定性分析转变提供理论参考。
五山堰河村的乡村建设和乡村旅游发展道路可以概括为“茶产业+文化+旅游”，以茶坛为代表的茶文化的乡村文化空间表征和空间再现是堰河成功的关键。旅游学术研究表明，旅游发展对地方文化的影响是双面的，一方面会促进地方文化自觉，另一方面又会导致地方文化的变迁。五山茶坛和堰河茶文化旅游的例子可以说明在文化空间生产理论的指导下，乡土文化会因旅游的发展而得到复兴，并为了更好的适应旅游发展的需要而进行创造性传承。
1 基于IDEPSO-SS的多工况可靠性分析
1.1 基于SS方法的失效概率求解
SS方法的基本思想是将小失效概率的求解过程，通过马尔可夫链蒙特卡罗（Markov Chain Monte Carlo，MCMC）方法生成条件样本，进而转换为一系列较大条件失效概率乘积的模拟过程，解决了MCS方法在求解工程小失效概率时效率较低的问题。
若自适应引入的中间事件满足嵌套关系，L1
L2
…
Lm=L，则目标失效概率PL为
式中：L为目标失效事件；Li为第i层的中间失效事件；P(L1)为首次失效的概率；P(Li|Li-1)为Li对应的条件概率。
2.1 补种和间苗 当幼苗出土后，应及时进行查苗补苗，以保证全苗间苗宜早，一般以子叶展开时间苗效果最好。如果苗期虫害严重时，可适当晚定苗。
在采用SS 方法求解之前，首先需要设置所有条件概率P(Li|Li-1)的初始值P0，进而确定各中间失效事件。若第i层模拟中使用的样本个数为N，则第i层条件概率的估计值为
式中：
为第i层模拟得到的失效概率估计值；ILi+1为示性函数；{ xk|k=1，2，…，N }为独立同分布样本集合。
将式（2）带入式（1）得
下行横缝缝长均值为2.750m，上行为2.267m，即下行横缝缝长大于上行；下行横缝缝宽均值为2.148mm，上行为2.034mm，即下行横缝缝宽大于上行。下行纵缝缝长均值为4.100m，上行为2.141m，即下行纵缝缝长大于上行；下行纵缝缝宽均值为2.006 mm，上行为0.667 mm，即下行纵缝缝宽大于上行。
式中：P0(i-1)为第i层的初始条件概率；NL(i)为第i层落入失效域的样本个数。
1.2 基于IDEPSO算法的失效概率优化
通过以上分析可知，每层样本个数和初始条件概率对失效概率估计值的影响较大。传统依据经验确定的每层样本个数和初始条件概率很难得到结构的最优失效概率。基于此，采用IDEPSO算法优化高阶响应面代理模型，获得最优每层样本个数和初始条件概率，从而计算单工况最优失效概率。
假设1个D维的搜索空间中包含了B个粒子，则第i（i=1，2，…，m）个粒子可表示为1个D维向量：Xi=(xi1 xi2 …xiD)，个体最优位置为Z=(Zi1 Zi2 …ZiD)，群体最优位置为G=(Gi1 Gi2 …GiD)。当2个最优位置确定后，可根据式（4）更新粒子的速度和位置。
式中：vid′为迭代过程中的粒子速度；t为迭代次数；Zid′为迭代过程中的个体最优位置；Xid′为迭代过程中的粒子位置；ZGd′为迭代过程中的群体最优位置；ω为惯性权重；c1和c2为学习因子；r1和r2为[0，1]内的均匀随机数。
PSO算法虽然能够以较大的概率收敛于全局最优解，并且避免复杂的遗传操作，但是由于最优粒子的速度和位置更新后可能为零，没有自主进化机制，因此极易陷入局部最优解，出现早熟停滞现象。
为验证IDEPSO算法的求解精度，基于图3所示各工况下高阶响应面代理模型计算得到的失效概率，并将其与MCS的差值定义为适应度函数。将表10中的数据代入IDEPSO算法中，对适应度函数进行寻优，并对正交试验设计结果进行分析，得到各工况下最优控制参数水平组合分别为：J=0.6，P=0.4，M=50，ωmax=1.3，ωmin=0.3；J=0.7，P=0.7，M=30，ωmax=1.5，ωmin=0.3；J=0.7，P=0.7，M=50，ωmax=1.2，ωmin=0.4；J=0.7，P=0.3，M=30，ωmax=1.1，ωmin=0.3。
其中，
式中：H i为变异操作中任一变异个体集；J为缩放因子，J ∈(0，1)；Xbest为当前代中最优个体；Xr1，Xr2和Xr3为当前代中随机选取的个体，其中r3为[0，1]内的均匀随机数；ζ为（0，1）之间均匀分布的随机数；f (xi)为第i个粒子的适应度；xmin和xmax为粒子中的最小值和最大值。
1.3 多工况结构失效概率计算
对于转向架构架而言，在设计阶段需要满足的工况较多，仅对单一工况进行可靠性分析得到的设计结果局限性较大。为全面分析结构设计的合理性，基于最优失效概率，采用二阶窄边界法（Ditlevsen法）和最优准则计算多工况下结构的最优失效次序和可靠度。
二阶窄边界法[15]由于考虑了各工况失效的相关性，被广泛应用于多工况下的可靠度计算。多工况失效概率Pf的表达式为
式中：Pij为工况i 与工况j 同时失效的概率，即联合失效概率。
攻坚三类重点对象。重点片区、重点经济薄弱村和重点低收入群体，是我省脱贫难度大的三类对象，更大力度集聚资源、集中力量攻坚的“三个重点”。对重点片区，明确按照整体帮扶规划，抓好排定重点项目和关键工程建设，改善交通出行、农田水利等基础设施，补齐发展短板。对123个省定重点经济薄弱村和“病残孤老灾”等重点低收入群体，要严格落实《关于加强农村重点低收入群体和重点经济薄弱村脱贫攻坚工作的意见》，强化有针对性的帮扶政策措施和工作举措。对相关政策支持和工作安排，尽量体现向“三个重点”倾斜。
基因漂流在自然界中极为普遍，对生物进化有重要意义。随着转基因作物的广泛种植，其对相关物种的生存状况影响引人思考。
由式（7）可知，Pij的计算精度直接影响多工况下失效概率的界宽。虽然采用数值积分和MCS方法能够较为精确地计算Pij值，但是不利于工程应用。
为了兼顾联合失效概率的求解精度和效率，采用基于几何分析的联合概率计算方法确定Pij[16]。Pij的表达式为
为解决这一问题，采用改进差分进化算法（IDE）优化PSO算法。研究表明，采用不同的变异策略使得算法的搜索能力不同，文献[12]采用的变异策略具有较强的全局搜索能力，而文献[13-14]采用的变异策略具有较强的局部搜索能力。为增强PSO算法的综合搜索能力，避免不同变异策略出现寻优早熟现象，提出IDE算法的混合变异策略为
式中：φ为标准正态分布的分布函数；Qn(n=i，j)为几何关系确定的2个工况间的失效概率，其求解表达式相同；βn为工况i 和工况j的最优失效概率对应的可靠性指标；ρ为工况失效时的相关系数；β0为中间参数；θn为各工况可靠性指标间的夹角。
将最优失效概率代入式（8）计算得到联合失效概率Pij，并代入式（7）得到多工况下失效概率的界宽估计值。虽然该值在计算精度上有所提高，但是受多工况下结构失效次序的影响较大。
1.4 最优失效次序的确定
由式（7）可知，多工况下失效概率界宽的最大下界即为最优可靠度（可靠度的最小上界）。为保证计算结果不受工况失效次序的影响，若对任意事件i不等式成立，则式（7）的下界可改写为
由式（9）可知：相比采用式（7）求解的多工况失效概率，采用式（9）求解时，工况的失效次序对Pij的值没有影响，其最大值即为所需的最大下界。为此，给出采用最优准则法求解式（9）的最大值，具体步骤如下：
（1）结合第1.1 和1.2节计算各工况的最优失效概率，根据1.3节计算各工况的联合失效概率Pij；
（2）选择失效概率中的最大值，并将其排序为λ；在余下的失效概率中选择工况序号ξ，使得pξ-pλξ的值最大，并将其排序为λ′；
在保证弹簧原长l0不变的情况下，可以通过改变O点的位置(即改变l1的长度)来改变连杆AB的初始角度θ0，同时弹簧与连杆之间的初始角度α0也随之改变。初始角度可由以下关系式求得：
导师要引导团队成员及时总结经验教训，把实践中获得的感性认识及时升华为理性认识，用于指导后期其他类型的实践。同样，这些经验总结对于其他团队的实践，以及其他类型的实践项目，同样具有借鉴意义。
（3）在除最大失效概率外的其他失效概率中选择工况序号ξ，使得pξ-pλξ-pλ′ξ 最大，并将其排序为λ′′；
（4）重复上述步骤，直至剩余的失效概率满足不等式pξ-pλξ
虽说我国已经实施了很多的发展策略，包括：丝路基金、亚投行等，但就实际情况而言，并未完全覆盖多样化的融资需求。在保障风险控制的前提下，我国金融机构应当积极探索全新的融资模式，不断吸引各项资本的参加，为“一带一路”项目建设奠定基础。与此同时，强化PPP模式的应用，联合民间资本与社会资本，为“一带一路”建设提供资金支持，政府部门应当强化各类形式合作模式的鼓励。
综上所述，可得基于IDEPSO-SS的多工况下转向架构架可靠性分析流程如图1所示。
图1 多工况下转向架构架可靠性分析流程
2 基于SS的构架失效概率响应面代理模型
2.1 构架静强度响应面函数
考虑不确定性因素对构架静强度的影响，建立不确定性参数关于静强度的响应面函数，是对构架进行可靠性分析的基础和前提。以某型转向架构架为对象进行分析，该构架由横梁、纵梁、牵引座等组成，分别采用壳单元和实体单元建模。考虑一系悬挂对转向架构架的影响，建立COMBIN14 单元模拟弹簧，单元刚度与一系悬挂相同，有限元模型如图2所示。
根据标准UIC 615-4《移动动力装置-转向架和走行装置-动力转向架构架强度试验方法》和EN 13749《铁路应用-轮对和转向架-转向架结构要求的规定方法》分别计算构架主体超常载荷、牵引拉杆纵向冲击载荷、电机振动冲击载荷、电机短路载荷和制动座载荷。依据设计需求对载荷进行组合，限于篇幅，仅选取4种组合工况为例进行分析，具体见表1和表2。
干旱指标应反映所关注的干旱类型，包括需水、供水、面对干旱表现出来的脆弱性和干旱潜在的影响。干旱指标的表现取决于数据的可获取性和数据质量，在应用于干旱前应该经过测试，在干旱后应该进行评价。干旱指标的选取必须考虑时间和空间的敏感性和特殊性，对不同时期和不同地区采用不同的干旱指标，还要综合考虑各类因素对干旱的影响，尽可能反映水循环及伴生过程中各组成部分之间复杂的相互关系及其影响，建立多指标的干旱指标体系。同时，还要确定与干旱水平在统计上具有一致性的干旱指标阈值，界定干旱的发展和消退。
图2 转向架构架的有限元模型
表1 超常工况载荷
组合工况1 4空气弹簧垂向力/kN左166.9 166.9右166.9 166.9空气弹簧横向力/kN 16横向止挡力/kN 108.6扭曲位移/mm 22
表2 特殊超常工况载荷
组合工况2 3空气弹簧垂向力/kN左166.9 166.9右166.9 166.9纵向力/kN 336 153齿轮箱吊杆反力/kN前123.6后-15.2
转向架构架结构静强度采用当量应力表示，在超常载荷工况下，母材当量应力应小于相应材料安全系数S=1.0时的许用应力，焊缝当量应力应小于相应材料安全系数S=1.1时的许用应力[17]。转向架构架组成材料及其力学性能见表3。
表3 构架组成材料的力学性能
?
根据构架的有限元模型及其组合工况计算静强度，基于分析结果确定对最大应力影响较大的设计参数，其统计特征见表4。表中：d1-d5为对各工况下对静强度影响较大的尺寸参数；F1-F5为对各工况下对静强度影响较大的载荷。
表4 设计参数的统计特征
设计参数限值下限均值上限结构尺寸/mm d1 11.50 12.00 12.50 d2 9.50 10.00 10.50 d3 7.50 8.00 8.50 d4 15.50 16.00 16.50 d5 15.50 16.00 16.50载荷/kN F1 97.74 108.60 119.46 F2 302.40 336.00 369.60 F3 111.24 123.60 135.96 F4 158.55 169.90 170.25 F5 158.55 169.90 170.25弹性模量E/GPa 195.70 206.00 216.30泊松比μ 0.25 0.30 0.35
基于表4提供的不确定性设计参数，建立各工况对应的构架静强度响应面函数，其表达式S1-S4分别为
2.2 静强度响应面函数的精度检验
为检验各工况下构建静强度响应面函数的拟合精度，对式（10）—式（13）进行方差分析（Analysis of Variance，ANOVA），其中工况1和工况4的部分分析结果分别见表5和表6。
表5 工况1下构架静强度响应面函数ANOVA结果
?
表6 工况4下构架静强度响应面函数ANOVA结果
?
通过方差分析可知，以上列举的2个工况下构架的响应面函数拟合精度较高，可以用于表征不确定性设计参数与输出的响应函数间的关系；且各设计参数的显著性较为明显，验证了不确定性设计参数选择的合理性。
2.3 构架失效概率响应面代理模型的建立
根据文献[18]确定每层样本个数和初始条件概率的取值范围，对构架进行D-最优试验设计。基于式（10）—式（13）、表3和试验设计，采用SS方法求解对应的失效概率，其中工况1和工况4下的试验设计结果分别见表7和表8。
表7 工况1下构架D-最优试验设计及响应值
?
表8 工况4下构架D-最优试验设计及响应值
?
采用麦夸特法对各工况下试验数据进行插值拟合，得到每层样本个数和初始条件概率关于失效概率的高阶响应面代理模型，如图3所示。为验证构架失效概率响应面模型的拟合精度，得到试验值与预测值的残差分布如图4所示。图4中：样本点表示由试验设计得到的试验值，灰色平面表示响应面的预测标准面。
由图3可知：每层样本个数和初始条件概率的不确定性对失效概率的影响较大，虽然与SS 方法的基本表达式相吻合，但是在采用SS方法过程中，依据经验确定每层样本个数和初始条件概率取值范围，由此计算得到的失效概率未必最优，降低了可靠性分析的精度。因此，寻求最优的每层样本个数和初始条件概率取值范围，对确定最优失效概率、保证可靠性分析的准确性具有重要的意义。
图3 不同工况下P0和N关于PL的响应面代理模型
图4 不同工况下响应面代理模型的残差分布
由图4可知：各工况残差的数量级均为10-4，说明响应面代理模型的拟合精度较高，可用于后续的最优失效概率的求解。
3 基于IDEPSO-SS的单工况下构架可靠性分析
3.1 IDEPSO算法最优控制参数的确定
在IDEPSO算法中，缩放因子J、变异概率P、种群规模M 及惯性权重的最大值ωmax和最小值ωmin等控制参数对其优化的效率和精度具有显著的影响[19]。为此，采用正交试验设计获取IDEPSO算法最优的控制参数组合。各控制参数及其水平值见表9，正交试验表L25(55)见表10。
2．2．2 发病时间分布 480例现患哮喘儿童发病好发于冬季和秋季，发病人数分别占42．71%（205例）和21．55%（102例），不定季节发作者占9．58%（46例），常年发作者占1．04%（5例），换季节发作者占14．58%（70例）。好发时间以午夜43．13%（207例）和清晨29．38%（141例）为主，睡前和午后发病分别为17．08%（82例）、1．04%（5例），无规律者为9．28%（45例）。
表9 不同水平下IDEPSO算法控制参数值
水平值1 2 3 4 5 J 0.5 0.6 0.7 0.8 0.9 P 0.3 0.4 0.5 0.6 0.7 M 30 50 70 80 100 ωmax 1.1 1.2 1.3 1.4 1.5 ωmin 0.3 0.4 0.5 0.6 0.7
表10 各控制参数正交试验设计值
试验号1 2…2 4 25 J 0.5 0.5…0.9 0.9 P 0.3 0.4…0.6 0.7 M 30 50…70 80 ωmax 1.1 1.2…1.2 1.3 ωmin 0.3 0.4…0.3 0.4
其中，
3.2 构架单工况下最优失效概率计算
采用各工况下最优控制参数再次对适应度函数进行优化，得到各工况下最优失效概率及每层样本个数和初始条件概率。不同优化算法下各工况的迭代曲线如图5所示。图中：GA为遗传算法；GAPSO为遗传算法优化的粒子群算法。
由图5可知：各工况下每层样本个数和初始条件概率的最优值分别为（805，0.25）（803，0.26）（801，0.11）和（770，0.182），可见由于适应度函数不同，得到的每层样本个数和初始条件概率最优值也不同，若按照经验对其取值得到的失效概率很难保证最优；与GAPSO，PSO 和GA算法比，IDEPSO算法在收敛速度和精度上最优，尤其在收敛速度上具有明显优势，同时避免了优化过程陷入局部最优解的情况。
图5 不同优化算法下各工况的迭代曲线
各工况下基于最优失效概率的构架累积分布函数如图6所示。IDEPSO-SS算法与MCS算法的结构失效概率（模拟精度）和调用函数次数（模拟效率）对比分别见表11和表12。
不耐贮运也与翠香吸钙能力差有关。钙增加耐贮藏性是已有的明确结论。因此，补钙可能对增强翠香耐贮运性或有一定作用。
图6 各工况下基于最优失效概率的构架累积分布函数
表11 IDEPSO-SS与MCS精度对比
方法MCS IDEPSO-SS误差/%各工况下失效概率/10-4工况1 105.00 104.81 0.181工况2 39.00 38.69 0.795工况3 34.00 33.60 1.176工况4 57.00 56.98 0.035
表12 IDEPSO-SS与MCS效率对比
方法MCS IDEPSO-SS效率提高/%各工况下调用函数次数工况1 9 543 353 96.3工况2 9 210 378 95.9工况3 9 439 364 96.1工况4 9 347 357 96.2
由表11和表12可知：与采用MCS算法相比，采用IDEPSO-SS算法得到的最优失效概率误差较小，精度较高，且计算效率显著提高，可以较为准确地表示构架在设计参数波动下的可靠性水平。
4 多工况下构架可靠性分析及最优失效次序确定
多工况结构可靠性分析的前提是获得各工况的联合失效概率。为分析相关系数对失效概率的影响，将图6得到的最优失效概率代入式（7）和式（8），计算得到相关系数ρ≥0.6时各工况下的联合失效概率，计算结果见表13。
Excel作为电子表格软件拥有强大的计算、分析功能。由于实测样本数据量很大，因此采用Excel 2007软件对这些数据进行统计、计算和分析，得出物种丰富度指数、物种多样性指数和物种均匀度指数等如下指标数据值(表1)。
根据表11和表13，采用Ditlevsen 方法计算多工况下构架的结构可靠度，并利用最优准则法确定最大下界及结构最优失效次序，计算结果如图7所示。图中：阴影部分为相关系数取0.6～0.9时多工况下构架失效概率界宽。
表13 各相关系数下工况间联合失效概率
相关系数0.60 0.65 0.70 0.75 0.80 0.85 0.90工况间联合失效概率/10-4 P12 12.0 14.0 17.0 20.0 24.0 28.0 33.0 P13 11.0 14.0 16.0 19.0 23.0 27.0 31.0 P14 16.0 19.0 23.0 27.0 32.0 38.0 44.0 P23 5.2 6.4 7.9 9.7 12.0 15.0 18.0 P24 6.9 8.5 10.0 13.0 15.0 19.0 23.0 P34 6.6 8.1 9.9 12.0 15.0 18.0 22.0
图7 多工况下构架可靠度
由图7可知：随着相关系数的增加，失效概率呈现下降趋势，表明考虑各工况间的相关性后，各工况类似于并联系统，导致失效概率降低；优化得到的失效概率界宽较窄，说明最优失效次序下的可靠度计算结果较为保守，降低了在设计阶段按照工况顺序依次失效导致过高估计结构可靠度的风险。对比图6可知，单工况的结构可靠度均高于多工况结构可靠度，表明仅进行单工况可靠性分析，计算结果过高，设计偏于危险。
相关系数ρ=0.6时最优准则法的求解过程见表14。表中：*代表构架各工况的最优失效次序。
表14 相关系数为0.60时最优准则法求解过程
ξ 1 2 3 4 Pf Pξ 10.481×10-3*3.869×10-3 3.360×10-3 5.698×10-3 18×10-3 Pξ-Pλξ 2.669×10-3 2.260×10-3 4.098×10-3*Pξ-Pλξ-Pλ′ξ 1.980×10-3*1.600×10-3 Pξ-Pλξ-Pλ′ξ-Pλ″ξ 1.460×10-3*
由表14可知，采用最优准则法得到的最优失效次序为工况1、工况4、工况2 和工况3，基于该次序计算得到的最优可靠度减少了Ditlevsen 方法的界宽，而且偏于安全。
5 结 论
（1）采用IDEPSO算法优化SS的每层样本个数和初始条件概率，提出基于改进SS 方法即IDEPSO-SS的单工况下最优失效概率计算方法，并且将其与Ditlevsen 方法、最优准则结合，提出多工况下最优可靠度及最优失效次序计算方法。
（2）采用IDEPSO-SS算法对转向架构架进行单工况下的可靠性分析，得到各工况下的最优失效概率。通过与传统方法对比，表明采用IDEPSOSS算法在保证计算精度的前提下，计算效率明显提高。
（3）构架的多工况下可靠性分析实例表明，本文方法较单工况下的可靠性分析，结果不仅更贴近工程实际，而且更加安全；获得的最优失效次序有利于获得结构在多种组合工况下的最优可靠度，为提高轨道车辆关键零部件的设计水平，保证可靠性分析的准确性提供了方法支撑。
参考文献
[1]智鹏鹏，李永华，陈秉智.考虑参数不确定性的转向架构架结构强度分析[J].中国机械工程，2019，30（1）：22-29.
（ZHI Pengpeng，LI Yonghua，CHEN Bingzhi.Structural Strength Analysis of Bogie Frames Considering Parameter Uncertainty[J].China Mechanical Engineering，2019，30（1）：22-29.in Chinese）
[2]卢耀辉，曾京，邬平波，等.铁道车辆转向架构架可靠性参数灵敏度分析[J].中国铁道科学，2010，31（1）：111-115.
（LU Yaohui，ZENG Jing，WU Pingbo，et al.Parametric Sensitivity Analysis for the Bogie Frame Reliability of Railway Vehicle[J].China Railway Science，2010，31（1）：111-115.in Chinese）
[3]滑林，吴梵，牟金磊.随机-非概率模型下的船体结构屈服强度可靠性分析[J].国防科技大学学报，2018，40（6）：177-182.
（HUA Lin，WU Fan，MU Jinlei.Reliability Analysis on Yielding Strength of Hull Structure Based on Random Non-Probabilistic Model[J].Journal of National University of Defense Technology，2018，40（6）：177-182.in Chinese）
[4]杨瑞刚，孟令军，胡立春.桥式起重机主梁结构多源不确定性混合可靠性分析[J].安全与环境学报，2018，18（4）：1246-1251.
（YANG Ruigang，MENG Lingjun，HU Lichun.Hybrid Reliability Analysis for the Multiple Source Uncertainty of the Dominant Girder Structure of the Bridge-Frame Crane[J].Journal of Safety and Environment，2018，18（4）：1246-1251.in Chinese）
[5]黄洪钟，刘征，米金华，等.混合不确定性下机床主轴可靠性建模与分析[J].中国科学：物理学、力学、天文学，2018，48（1）：42-53.
（HUANG Hongzhong，LIU Zheng，MI Jinhua，et al.Reliability Modeling and Analysis of Heavy-Duty CNC Machine Tool Spindle under Hybrid Uncertainty[J].Scientia Sinica：Physica，Mechanica& Astronomica，2018，48（1）：42-53.in Chinese）
[6]曾光，姜潮，倪冰雨.工程机械柔性机械臂振动的混合可靠性分析[J].中国机械工程，2017，28（12）：1400-1405，1412.
（ZENG Guang，JIANG Chao，NI Bingyu.Hybrid Reliability Analysis for Flexible Manipulator Vibrations in Engineering Machinery[J].China Mechanical Engineering，2017，28（12）：1400-1405，1412.in Chinese）
[7]ZHI P P，LI Y H，CHEN B Z，et al.Fuzzy Optimization Design-Based Multi-Level Response Surface of Bogie Frame[J].International Journal of Structural Integrity，2019，10（2）：134-148.
[8]ECHARD B，GAYTON N，LEMAIRE M.AK-MCS：An Active Learning Reliability Method Combining Kriging and Monte Carlo Simulation[J].Structural Safety，2011，33（2）：145-154.
[9]HUANG X X，CHEN J Q，ZHU H P.Assessing Small Failure Probabilities by AK-SS：An Active Learning Method Combining Kriging and Subset Simulation[J].Structural Safety，2016，59：86-95.
[10]LI H S，MA Y Z，CAO Z J.A Generalized Subset Simulation Approach for Estimating Small Failure Probabilities of Multiple Stochastic Responses[J].Computers&Structures，2015，153：239-251.
[11]XIAO N C，YUAN K，ZHOU C N.Adaptive Kriging-Based Efficient Reliability Method for Structural Systems with Multiple Failure Modes and Mixed Variables[J].Computer Methods in Applied Mechanics and Engineering，2020，359：112649.
[12]林健，黄家才，汪木兰，等.DEPSO算法在永磁直线电机结构优化设计中的应用[J].中国机械工程，2012，23（24）：2908-2912.
（LIN Jian，HUANG Jiacai，WANG Mulan，et al.Application of DEPSO Algorithm to Permanent Magnet Linear Motor Structure Optimization Design[J].China Mechanical Engineering，2012，23（24）：2908-2912.in Chinese）
[13]王军祥，董建华，陈四利.基于DEPSO 混合智能算法的岩土体应力-渗流-损伤耦合模型多参数反演研究[J].应用基础与工程科学学报，2018，26（4）：872-887.
（WANG Junxiang，DONG Jianhua，CHEN Sili.Multi-Parameters Inversion of Stress-Seepage-Damage Coupling Model Based on DEPSO Intelligent Algorithm[J].Journal of Basic Science and Engineering，2018，26（4）：872-887.in Chinese）
[14]陈璟华，邱明晋，唐俊杰，等.基于改进差分进化和粒子群混合算法的电力系统最优潮流计算[J].广东工业大学学报，2017，34（5）：22-28.
（CHEN Jinghua，QIU Mingjin，TANG Junjie，et al.A Hybrid Algorithm Based on Improved Differential Evolution and Particle Swarm Optimization for Power System Optimal Power Flow Calculation[J].Journal of Guangdong University of Technology，2017，34（5）：22-28.in Chinese）
[15]DITLEVSEN O.Narrow Reliability Bounds for Structural Systems[J].Journal of Structural Mechanics，1979，7（4）：453-472.
[16]姚继涛，赵国藩，浦聿修.二维标准正态联合概率的计算[J].建筑结构学报，1996，17（4）：10-19.
（YAO Jitao，ZHAO Guofan，PU Yuxiu.Methods for Computing Two-Dimension Joint Probability of Standard Normal Random Vectors[J].Journal of Building Structures，1996，17（4）：10-19.in Chinese）
[17]FENG Y S.A Method for Computing Structural System Reliability with High Accuracy[J].Computers & Structures，1989，33（1）：1-5.
[18]AU S K，BECK J L.Estimation of Small Failure Probabilities in High Dimensions by Subset Simulation[J].Probabilistic Engineering Mechanics，2001，16（4）：263-277.
[19]李永华，盛自强，宫琦.基于GAPSO-RBFNN的动车组电机吊架多目标稳健优化设计[J].中国铁道科学，2020，41（3）：103-110.
（LI Yonghua，SHENG Ziqiang，GONG Qi.GAPSO-RBFNN-Based Multi-Objective Robust Optimal Design for Motor Hanger of EMU[J].China Railway Science，2020，41（3）：103-110.in Chinese）