基于乘子法的结构可靠性分析.docxVIP

基于乘子法的结构可靠性分析 0 不确定结构可靠性分析方法 在结构设计中，材料不均匀、安装误差和其他问题的存在将给结构性能带来各种影响。传统结构设计没有考虑施工过程中的不确定因素，而实际施工中载荷和材料性质等不断变化，这些不确定因素会使失效概率不断增大。为了确保产品的结构可靠性，有必要对不确定结构可靠度分析方法进行研究。目前，大量研究者为了解决实际工程问题，发展和提出了许多先进的计算方法 基于此，本文提出一种基于乘子法的结构可靠分析方法，该方法稳定有效地解决不确定结构可靠度指标求解问题，数值算例和工程实例验证了所提出的方法在解决可靠性分析问题时具有较好的稳定性和有效性。 1 结构可靠度计算方法 在最优化理论中，罚函数法是求解约束优化问题的一个重要方法，把约束优化问题转化成一个或一系列的无约束优化问题，通过无约束问题来求解约束极值问题。通常采用的方法是在原目标函数加上一个“惩罚项”来迫使迭代点逼近可行域。通常情况下，必须使用罚参数与约束条件来构造惩罚项。当前迭代点不是可行点时，就要构建惩罚项，并且使得惩罚随着不可行点到可行域距离的增大而变大；可行点处不实施惩罚。惩罚项强制迭代点逼近可行域，最终落入可行域。外罚函数法结构简单，可以直接调用无约束优化算法的通用程序,但是也具有局限性，罚参数σ 对于等式约束问题，数学模型可以表示为： 其中，h(x)=[h 其中，乘子向量λ=（λ 对于任意的x∈D，有 由式（4）可知，如果乘子向量λ 使用外罚函数法求解，式（5）的增广目标函数为 可以得出，当σ0适当大时，x 选取一个 而在原问题的KT点（x 为使{x 当h(x 在结构可靠度分析中，结构可靠度指标β的几何意义是标准正态坐标系中原点到极限状态曲面的最短距离，垂足即为验算点，验算点到标准坐标系的原点的距离即为结构可靠度指标。 假设有n个独立的影响结构可靠度的非标准正态分布的随机变量，该结构的极限状态函数可表示为 设计验算点为x 通过优化处理可靠度指标才能确定验算点，先设β为极限状态曲面上点P∑X 基于上述讨论，给定功能方程及其分布概型、随机变量的基本参数，即可通过迭代求解并得到可靠度指标，迭代计算步骤如下： 2）求解子问题。以x 3）检验终止条件。若‖h(x 4）更新罚参数。若‖h(x 5）更新乘子向量，令λ 6）令k:=k+1，转步骤2）。 2 可靠性计算示例 2.1 可靠度检验 极限状态方程为 x 根据表2和表3可知，与可靠度指标β越大、失效概率P 从图1～图7中可看出，本文方法和AFOSM法在求解可靠度指标时都能得到稳定和有效的解，但本文算法很快收敛，且得到精度较高的解。总之，本文所提出的方法对于不确定结构可靠度优化设计问题具有较好的适用性。 2.2 基于b柱前端最大允许路径b的可靠度分析方法 在交通事故中，汽车发生侧面碰撞时缓冲空间小，车体变形大。侧围结构侵入量、侵入速度和侵入形态等对乘员安全性具有较大影响，B柱最大侵入量是衡量汽车侧碰中耐撞性的重要指标。可移动障壁以50 km/h的速度从侧面撞向汽车侧身，侧碰有限元模型如图8所示，车辆侧碰设计变量如图9所示，以汽车关键设计尺寸中左侧车体框架厚度x 借助实验设计，运用拉丁超立方设计方法从设计空间中选28个样本点，调用有限元模型分析，构建最大侵入量b的响应面函数。B柱下端最大允许侵入量b 列出以下3种参数取值，求解可靠度指标。 1）取第1组中各随机变量的均值为区间的中间值，BPA数值为标准差，具体的随机变量分布参数取值情况见表6。 根据本文方法，经过3次迭代，求得可靠度指标β=2.1421，失效概率P 根据表7和表8可知，可靠度指标β越大，失效概率P 2）取第2组中各随机变量的均值为区间的中间值，BPA数值为标准差，具体的随机变量分布参数取值情况见表10。 根据本文方法，经过3次迭代，求得可靠度指标β=0.1298，失效概率P 根据表11和表12可知，计算数据和可靠度指标β越大，失效概率P 3）取第3组中各随机变量的均值为区间的中间值，BPA数值为标准差，具体的随机变量分布参数取值情况见表14。 根据本文方法，经过3次迭代，求得可靠度指标β=0.4721，失效概率P 根据表15和表16可知，计算数据和可靠度指标β越大、失效概率P 经过上述分析，汽车侧碰在该设计变量和侵入量前提下是比较安全的，比较了AFOSM法和本文方法，发现本文方法收敛较快并能得到较高精度的解。总之，本文所提出的方法在求解汽车侧碰结构可靠度优化设计问题时，能快速收敛且得到稳定和有效的解，具有较好的适用性。 3 考虑稳定性问题的外函数算法 本文针对具有一定非线性程度功能函数的可靠性分析问题，提出了一种基于乘子法的结构可靠性分析方法。该方法充分利用外函数求解非线性优化问题的优点，从乘子法的KT条件出发，快速求解追踪到最优点。