一种基于乘子法的结构可靠性分析方法

徐柳， 王锬， 杜义贤， 刘晋玮， 黄文超， 王林军

（三峡大学a.水电机械设备设计与维护湖北省重点实验室；b.机械与动力学院，湖北 宜昌443002）

0 引 言

在工程设计中，材料不均匀、安装误差等问题的存在，往往会对结构性能造成各种影响甚至破坏。传统结构设计没有考虑施工过程中的不确定因素，而实际施工中载荷和材料性质等不断变化，这些不确定因素会使失效概率不断增大。为了确保产品的结构可靠性，有必要对不确定结构可靠度分析方法进行研究。目前，大量研究者为了解决实际工程问题，发展和提出了许多先进的计算方法[1-6]。徐强等[7]提出了大坝体系可靠度改进计算方法；黄逸群等[8]对实际工程问题中隧道型钢喷混凝土初期支护进行可靠度计算；高东川等[9]使用一次二阶矩理论中心点法预测了支护后的软岩巷道围岩的稳定性；肖宇峰[10]提出了一种基于子网同构判定的高效计算方法；李怀龙等[11]计算了轨道板横向宽轨枕单元的轨下截面和板中截面的可靠度；李奎等[12]研究了深埋隧道素混凝土衬砌可靠度计算模型，并进行了相关可靠性分析计算；王鹏等[13]提出了基于应力-强度模型的DTECS -2 设备可靠度计算方法；Hamed Fazlollahtabara等[14]提出了一个用于可靠度计算的集成马尔可夫和反向传播神经网络。通过分析发现，乘子法在处理实际工程结构可靠性分析问题时具有较好的效果。

基于此，本文提出一种基于乘子法的结构可靠分析方法，该方法稳定有效地解决不确定结构可靠度指标求解问题，数值算例和工程实例验证了所提出的方法在解决可靠性分析问题时具有较好的稳定性和有效性。

1 乘子法理论介绍

在最优化理论中，罚函数法是求解约束优化问题的一个重要方法，把约束优化问题转化成一个或一系列的无约束优化问题，通过无约束问题来求解约束极值问题。通常采用的方法是在原目标函数加上一个“惩罚项”来迫使迭代点逼近可行域。通常情况下，必须使用罚参数与约束条件来构造惩罚项。当前迭代点不是可行点时，就要构建惩罚项，并且使得惩罚随着不可行点到可行域距离的增大而变大；可行点处不实施惩罚。惩罚项强制迭代点逼近可行域，最终落入可行域。外罚函数法结构简单，可以直接调用无约束优化算法的通用程序, 但是也具有局限性，罚参数σk→+∞,会使得增广目标函数变得“越来越病态”。为了改善这种局限，在乘子法中引入拉格朗日函数，并加上适当的罚函数。Powell和Hestenes在1969年针对等式约束优化问题同时独立提出了乘子法，后来Rockfellar在1973年将乘子法推广到求解不等式约束优化问题。从原问题的拉格朗日函数出发,加上适当的罚函数,可以将原问题转化为一系列的无约束优化子问题。

对于等式约束问题，数学模型可以表示为：

其中，h（x）=[h1（x）,h2（x）,…,hl（x）]T，可行域D={x∈Rn|h（x）=0}，拉格朗日函数为

2 可靠度算例

2.1 数值算例

极限状态方程为

x1与x2相互独立，且均服从正态分布。变量分布参数取值如表1所示。根据本文方法，经过4次迭代，求得可靠度指标β=2.3302，迭代计算数据见表2。使用改进一次二阶矩（AFOSM）法检验，经过11次迭代，求得可靠度指标β=2.3302，迭代计算数据见表3。两种方法的可靠指标计算结果及失效概率见表4。

表1 各随机变量分布参数取值情况

表2 乘子法迭代计算数据

表3 AFOSM法迭代计算数据

表4 两种方法的可靠指标计算结果及失效概率

根据表2和表3可知，与可靠度指标β越大、失效概率Pf越小的一般规律相吻合。随着迭代过程的进行，可靠度与失效概率均逐渐趋于稳定，本文方法只需经过4次迭代即可得出结果，而AFOSM法则需要经过11次迭代才能得出结果。根据表4，使用本文方法和AFOSM法，最终都收敛于点（11.1855,1.6549），可靠度指标β=2.3302，失效概率Pf=0.0099。两种方法的可靠度指标和迭代次数之间的关系如图1所示。本文方法的目标函数收敛曲线如图2所示，目标函数等值线如图3所示，目标函数等值线局部放大图如图4所示。AFOSM法的目标函数曲线图如图5所示,目标函数等值线如图6所示，目标函数等值线局部放大图如图7所示。

图1 迭代次数和可靠度之间的关系

图2 乘子法目标函数曲线图

图3 乘子法目标函数等值线

图4 乘子法目标函数等值线局部放大图

从图1～图7中可看出，本文方法和AFOSM法在求解可靠度指标时都能得到稳定和有效的解，但本文算法很快收敛，且得到精度较高的解。总之，本文所提出的方法对于不确定结构可靠度优化设计问题具有较好的适用性。

图5 AFOSM法目标函数曲线图

图6 AFOSM 法目标函数等值线

图7 AFOSM法目标函数等值线局部放大图

2.2 汽车侧碰分析

在交通事故中，汽车发生侧面碰撞时缓冲空间小，车体变形大。侧围结构侵入量、侵入速度和侵入形态等对乘员安全性具有较大影响，B柱最大侵入量是衡量汽车侧碰中耐撞性的重要指标。可移动障壁以50 km/h的速度从侧面撞向汽车侧身，侧碰有限元模型如图8所示，车辆侧碰设计变量如图9所示，以汽车关键设计尺寸中左侧车体框架厚度x1、左前车门防撞梁厚度x2、左侧A柱内板厚度x3和外板厚度x4为设计变量，边缘BPA如表5所示，x1、x2、x3、x4相互独立[15]。

图8 汽车侧碰有限元模型

图9 车辆侧碰设计变量

借助实验设计，运用拉丁超立方设计方法从设计空间中选28个样本点，调用有限元模型分析，构建最大侵入量b的响应面函数。B柱下端最大允许侵入量bmax=290 mm，可建立如下极限状态方程：

列出以下3种参数取值，求解可靠度指标。

3 结 论

本文针对具有一定非线性程度功能函数的可靠性分析问题，提出了一种基于乘子法的结构可靠性分析方法。该方法充分利用外函数求解非线性优化问题的优点，从乘子法的KT条件出发，快速求解追踪到最优点。数值算例和工程算例的计算结果表明，相对于AFOSM方法，本文所提出的方法对求解非线性程度较高的极限状态方程具有较好的收敛速率。同时，本文方法将来可应用于一些复杂工程问题结构的可靠度反问题中，还可以应用在串并联可靠性分析求解问题。