结构可靠度分析与设计的编程实践.docx

结构可靠度分析与设计的编程实践
摘 要：基于《荷载与结构设计方法》中讲授的结构可靠性分析与设计的基本原理，对课程中所给出的例题利用Matlab软件编制了相应的计算机程序。通过此次编程实践，加深了自己对结构可靠性分析与设计的认识了理解。
1.引言
不确定性是自然界中普遍存在的一种客观现象，工程设计中的不确定性有多种不同的形式，人们认识最早、目前得到广泛应用的是随机性，人们用概率的方法研究。结构可靠度方法则是结构可靠性设计方法中的重要一项。结构的可靠度是结构在规定时间内，规定条件下结构能够完成预定功能的概率。从简单到复杂或精确程度的不同，先后提出的可靠度计算方法有MVFOSM（一次二阶矩法）、AFOSM（改进的一次二阶矩法）、RF、MCS（蒙特卡洛数值模拟法）等方法。本文根据《荷载与结构设计方法》中讲授的结构可靠度分析与设计的基本原理，对课程中所给出的例题利用Matlab软件编制了相应的计算机程序。
2.结构可靠度分析的基本原理
当用概率描述结构的可靠性时，就需要根据结构中基本随机变形变量或综合随机变量的概率分布进行计算。在实际工程中，结构的功能函数往往是由多个随机变量组成的非线性函数，而且这些随机变量并不都服从正态分布或对数正态分布，因此不能直接采用相应的公式计算可靠指标，而需要作出某些近似简化后进行计算。下面本文将介绍分析结构可靠度的几种常用方法：
2.1.AFOSM法（改进的一次二阶矩法）
1974年Hasofer和Lind科学地对可靠指标进行了定义，引入了验算点的概念，使得一阶二次矩模式有了进一步的发展，由于分析中要迭代求解验算点，验算点是可靠分析中的一个关键点，所以人们又将这种方法称为验算点法。
随机变量服从正态分布的情形。、
2.1.1.功能函数为线性函数

式中：(i=0,1,2,…n)为常数。
为进一步在标准正太坐标系中研究可靠指标的几何意义，按下式将随机变量变换为标准正态随机变量

则结构功能函数可由表示为：

从而功能函数的平均值和标准差为：

可靠度指标为：

由上式计算的可靠指标与结构失效概率的对应关系：

2.1.2.功能函数为非线性函数
为计算可靠指标，将非线性功能函数在验算点处泰勒级数展开，则结构功能函数的一次展开式为：

则Z的平均值和标准差可分别近似表示为

由此可求得可靠指标：

2.2.RF法
当随机变量不服从正态分布时，则需要将其进行当量正态化将随机变量等效为正态随机变量，当量正态化的条件是，在验算点处使非正态随机变量的概率分布函数值与当量正态随机变量的概率分布函数值相等，的概率密度函数值与当量正态随机变量的概率密度函数值相等，即

则：

当量正态化后：

上述三式构成非正态随机变量情况下可靠指标的迭代计算公式。
3.结构可靠性设计的基本原理
结构设计首先要保证结构的安全性，其次是保证结构的适用性、耐久性和偶然作用下的整体稳定性。为实现这些目的，就要对结构进行合理的设计。
3.1结构可靠度校准
如果按照结构设计规范规定的设计表达式进行可靠度计算，按规范设计的结构或结构构件的可靠度指标，在这种情况下反映的是一种结构或一类结构构件的可靠指标，代表了设计规范的可靠水平，这一过程称为结构可靠度校准。
结构可靠度校准的步骤：
(1).假设一个荷载效应比；
(2).确定构件抗力的特征值：

(3).确定基本变量的均值和标准差：
均值：
标准差：
(4). 确定极限状态方程：

(5).利用RF法求解可靠度指标；
(6).变化荷载效应比，计算不同可靠指标的均值
3.2分项系数的校准
如果给定了结构或构件的设计表达式且已知变量的概率分布和统计参数，即可求得按设计表达式设计的结构或构件的可靠指标，可靠指标至于作用效应标准值的比值有关，而与作用效应标准值本身的大小无关。这样可以通过先设定作用分项系数、抗力分项系数和作用组合系数，然后计算相应的可靠指标，并与目标可靠指标进行比较，从而确定一组使计算的可靠指标与目标可靠指标最接近的系数值，即确定的、、和使下式最小

下面是用迭代算法求分项系数的过程：
（1）确定极限状态方程和设计方程的表达式。
①确定基本变量的概率分布类型和近似参数
②最多有两个未知的均值需要求解。一个是抗力的均值，另一个是相应的荷载效应均值，荷载效应比率用于建立荷载与抗力均值之间的联系
③第一次迭代，利用极限状态方程 获得两个未知变量均值的关系。
（2）利用上述均值获得初始的设计点。
（3）对于非正态分布对应的设计点取值，利用当量正态化公式确定等效的正态均值和正态标准差。
（4）计算n个变量的方向余弦

（5）计算n个变量设计点的坐标值

（6）通过求解极限状态方程确定两个未知均值新的关系

（7）重复3-6步直到收敛
（8）当收敛条件满足，计算分项系数

4.例题的程序运行结果展示
4.1Homework2.1的结果展示
The example 5.4

The problem 5.3
迭代次数
1
2

3.1537
3.1537

-0.9997
-0.9997

0.0232
0.0232
4.2.Homework2.2的结果展示
The example 3.5

The example 5.9

The example 5.10

The example 5.11

The problem 5.4

4.3Homework2.3的结果展示
The example 5.11
The problem 5.4
4.4 Homework3的结果展示
The example 6.3

The example 6.7

5.结论与讨论
结构可靠度设计与分析方法，从简单到复杂或精确程度的不同，先后提出的可靠度计算方法有MVFOSM、AFOSM、RF、MCS等方法。通过此次编程实践，加深了自己对结构可靠度分析与设计的认识了理解。
6.参考文献
1. 中华人民共和国国家标准. 工程结构可靠性设计统一标准（GB50153—2008），北京：中国建筑工业出版社，2009.
2. 中华人民共和国国家标准. 建筑结构可靠度设计统一标准（GB50068—2001），北京：中国建筑工业出版社，2001.
3. 张明著. 结构可靠度分析 — 方法与程序. 北京：科学出版社，2009年.
4. 李国强等编著. 工程结构荷载与可靠度设计原理（第三版）. 北京：中国建筑工业出版社，2005.
5. 贡金鑫，魏巍巍. 工程结构可靠性设计原理. 北京：机械工业出版社，2007.
6. Nowak, A. S., and K. R. Collins (2000). Reliability of Structures. McGraw-Hill, New York, NY. 张川导读，重庆：重庆大学出版社，2005年3月

附录
附件1：Homework 2
%AFOSM homework 2_1_1
% limit state g=EI-310.5w
% E ,I and w are all normal variables

while

end
'2_1_1'

*******************************************************************************

*******************************************************************************
%RF 2_2_1
%limit state g=R-G-Q
%R lognorm G normal Q extrem I

while
%R服从对数正态分布

%Q服从极值Ⅰ型分布

%G服从正态分布

%正态化后极限状态方程变为g=R1-G1-Q1

end
'2_2_1'

*******************************************************************************

%RF 2_2_21
%limit state g=R-Q
%R lognorm Q extrem I

while
%R服从对数正态分布

%Q服从极值Ⅰ型分布

%正态化后极限状态方程变为g=R1-Q1

end
'2_2_21'

*******************************************************************************
%RF 2_2_22
%limit state g=R-Q
%R lognorm Q normal

'mu_R=''deta_R='
'mu_Q=''deta_Q='

while
%R服从对数正态分布

%Q服从正态分布

%正态化后极限状态方程变为g=R1-Q1

end
'2_2_22'

*******************************************************************************
%RF 2_2_31
%g=Z*Fy-M
% Z are normal variables; Fy are lognormal variables; M are extrem type I variables

while
%Z服从正态分布

%Fy服从对数正态分布

%M服从极值Ⅰ型分布

%正态化后极限状态方程变为g=Z1*Fy1*M1
'(mu_Z1+alpha_Z1*x*sigma_Z1)*(mu_Fy1+alpha_Fy1*x*sigma_Fy1)-(mu_M1+alpha_M1*x*sigma_M1)=0'

end
'2_2_31'

*******************************************************************************
%RF 2_2_32
%limit state g=3Y-X
%Y lognormal; X extrem type I

while
%Y服从对数正态分布

%X服从极值Ⅰ型分布

%正态化后极限状态方程变为g=3Y1-X1
'3*(mu_Y1+alpha_Y1*x*sigma_Y1)-(mu_M1+alpha_X1*x*sigma_X1)=0'

end
'2_2_32'

*******************************************************************************
%2_31
Monte Carlo Simulation
%judgement:Fy*Z-M<0 (limit state g=Fy*Z-M)
%Fy are lognormal variables Z are normal variables M are extreme i variables

*******************************************************************************%2_32
Monte Carlo Simulation homework 2_32
%judgement3Y-X<0 (limit state g=3Y-X)
%Y are lognormal variables X are extreme i variables

*******************************************************************************

附录2 Homework 3
%homework 3_1
%reliability-based partial safety factor design
%R are lognorm variables Q are extreme I viariables G are normal variables

while

%R服从对数正态分布

%Q服从极值Ⅰ型分布

%G服从正态分布

%正态化后极限状态方程变为g=R1-G1-Q1

end

*******************************************************************************%RF3_2
%g=R-Q
%R、Q值未给出，不妨令R=Q=1

while

end