拉格朗日乘子法求解结构可靠度.pdfVIP

貉牵气肯 固体力学学报 ACTAMECHANICA SOLIDASINICA 拉格朗日乘子法求解结构可靠度+ 杨绿峰 赵艳林 韦树英 (广四大学土木建筑上程学院，南宁，530004) 摘要在结构可靠度分析中．由于极限状态方程通常是多变量、非线性的，必须借助迭代 法，通过烦琐重复的计算近似确定设计验算点和可靠指标该文提出求解设计验算点和结构可 靠度指标的一种方法利用极值理论，在识别拉格朗13乘子后，建立了一组方程，用于求解设计 验算点，进而求解结构可靠指标和可靠度．实例表明，这种方法可以简化结构可靠度的计算过 程，同是提高计算精度 关键词结构可靠度，可靠指标，拉格朗日乘子，设计验算点 O引言 工程结构中普遍存在不确定性，包括模糊性’11和随机性．以可靠度理论为基础的建筑 结构设计规范充分考虑了结构的随机性．但是，可靠性分析中的极限状态方程通常是多变量 非线性的，因而在极限状态面上寻求设计验算点并不是一件容易的工作．人们通常采用迭代 法，通过烦琐的重复计算来确定验算点和结构可靠指标，进而求解可靠度．一次二阶矩法、Jc 法”J、实用分析法”o及映射变换法”o都采用了类似的迭代方法，甚至在利用随机有限元 法”’计算结构可靠度时，迭代方法仍是人们倚重的手段． 本文利用极值理论和拉格朗日乘子法，在识别拉氏乘子后建立了一组方程．求解这组方 程可以确定设计验算点．并据此可以求解结构可靠指标和 结构可靠度，由于在计算设计验算点过程中．没有涉及到可 靠指标．从而大大简化了计算过程，提高了计算精度． 1拉氏乘子及其识别 如图1所示，结构可靠指标的几何含义，是标准正态空 间坐标系Oz，x：…‰中原点0到极限状态面(线)的最短距 离OP’．点P’称为设计验算点．在结构可靠度分析中，可靠 指标口的计算常转化为求坐标原点0到设计验算点P’的距 图1坐标系 离，即计算直线OP’的长度．由于极限状态面通常是曲面， 所以在确定点P‘的位置和OP’的长度时，通常采用试算、迭代方法．由于这种方法同时涉 及到可靠指标口的求解和验算点P+(zj，z；，x；)的确定，过程较繁，且耗时较多． 本文提出的方法将P’的确定与口的计算分二步进行．首先根据极值理论，确定P+点在 极限状态面是的位置(#j，*；，…．z：)；然后计算OP’的长度求得可靠指标届在n维笛卡尔 空间坐标系0以z：…z。中，点P’到原点0的距离的平方为 西十百千人才工程专项墓金资助． 专辑 杨绿峰等： 拉格朗日乘子法求解结构可靠度 ·233· d2=lOP’I (1) 2=∑(xj)2 由于点P’位于极限状态面上，所以应满足极限状态方程 G(XI*，z；，…，z：)=G(P’)=0 (2) 引人拉格朗日乘子A，可以根据式(1)、(2)构造新的函数 (3) D(P‘)=D(xi，z；，…，z：)=∑(xj)2+2,G(P+) 点P‘应满足D(P+)的驻值条件 (4a) 掣：2xj+AG。(P·)：o，江1，2，…，。 塑(￡：!：c(P．1：0 (4b) dA 式中