验算点法在工程结构可靠度编程中的运用.docVIP

验算点法在工程结构可靠度编程中的运用 （北京航空航天大学 交通科学与工程学院，北京，100191） 摘 要：用验算点法对可靠度编程中的关键问题进行了探讨并给出了解决方法。这些问题主要包括随机变量服从正态分布的情形且功能函数为非线性和随机变量不服从正态分布时的当量正态化方法。针对这两种情况运用了两个编程算例来说明验算点法在可靠度分析中的运用。 关键字：验算点法、可靠指标、正态分布函数、非正态随机变量、当量正态化 前 言 结构可靠度为结构在规定的时间内和规定的条件下完成预定功能的概率。结构可靠性理论的研究，起源于对结构设计、施工和使用过程中存在的不确定性的认识，以及结构设计风险决策理论中计算结构失效概率的需要。结构可靠度计算方法有一次二阶矩方法、二次二阶矩方法、蒙特卡洛方法及其他方法。一次二阶矩方法又分为中心点法和验算点法，其中验算点法是目前可靠度分析中最常用的方法。由于这两种方法都是将非线性功能函数作为一次泰勒级数展开，并使用了随机变量的平均值（一阶矩）和方差（二阶矩），故称为一次二阶矩方法。利用验算点法计算结构的可靠指标时，需要预先知道验算点的坐标值，而对于非线性结构功能函数和非正态随机变量的情形，验算点坐标值是不能预先求得的，因此一般需要迭代求解。 随机变量服从正态分布的情形 功能函数为线性函数 功能函数随机变量是一个正态随机变量，其概率密度函数和U的密度曲线如图1示 图1 一个随机变量时的可靠指标（左图为正态随机变量，右图为标准正态随机变量） 假定存在n个相互独立的随机变量，，其均值为，标准差为结构功能函数为： （1） 其中为常数 将随机变量变换为标准正态随机变量 （2） 则由（1）表示的功能函数表示成 从而功能函数的平均值和标准差表示为 按照严格的可靠度指标定义 （3） 可靠度指标和结构失效概率存在精确的对应关系 对极限状态方程 两端同时除以得到： （4） 与公式（3）比较，有 （5） 令 公式（5）可以写成： （6） 公式（6）表示的是一法线式的直线方程，为法线与坐标轴夹角余弦 图2 可靠度指标的几何意义及验算点 验算点在空间（标准正态空间）表示为： 在空间表示为： 两者之间的关系为： 根据几何关系有： 在空间，验算点坐标值： 通常表示为： 功能函数为线性函数 假定随机变量服从正态分布，但结构功能函数不再是线性函数，显然，这时精确求解的平均值和标准差是非常困然的。同结构功能函数为非线性的情形一样，如果将可靠指标定义为标准正态坐标系中坐标原点到极限状态曲面的距离，垂足为验算点，则不管结构极限状态方程的数学表达形式如何，只要具有相同的力学或物理含义，在标准正态坐标系中，所表示的都将是同一个曲面，曲面上与坐标原点距离最近的点也只有一个。因而，所得到的可靠指标是唯一的，不像中心点法那样，随结构极限状态方程数学表达式的形式而变。 图3验算点取法 如果验算点已知 可以在该点一次项展开： 其均值和标准差为： 所以可靠度指标： 实际上验算点不可知，需要补充条件： 对比表达式得到： 编程算例 假定结构功能函数为。随机变量的平均值和标准差分别为，，，，均服从正态分布。用验算点法就算可靠指标和失效概率，允许误差取 2.3.1 算法分析 假定验算点，一般取 所以， 计算 功能函数对的一阶偏导数为， 所以 计算 重新计算验算点 则， 若，为规定的允许误差，则停止迭代，所求即为要求的可靠 指标；否则，取，转2）继续迭代，当验算点误差小于，结束。 2.3.2 源程序 #include iostream #include cmath using namespace std; //example 3-4 double Beta(double x1, double x2); double CosOx1(double x1, double x2); double CosOx2(double x1, double x2); double X1 = 38.0; double X2 = 54.0; double e = 1.0e-3; double Beta(double x1, double x2) { return (54.0 * x1 + 38.0 * x2 - x1 * x2 - 1000.0) / sqrt(pow(3.8 * x2, 2) + pow(5.4 * x1, 2)); } double CosOx1(double x1, double x2) { retur