基于MATLAB截尾分布下JC法计算可靠度.docVIP

基于MATLAB截尾分布下JC法计算可靠度摘 要： 在水工结构可靠度分析中，随机变量的分布形式常因几何尺寸、物理环境等条件限制，传统JC法已经不适用，因此需要对部分变量进行截尾分布处理。在此借助MATLAB丰富的函数资源，编制出截尾分布处理后的改进JC法求解程序。算例研究表明该方法的可行性，以及程序的简易、实用性。 关键词： 可靠度； 截尾分布； JC法； MATLAB 中图分类号： TV314 文献标识码： A 文章编号： 1009-8631（2011）01-0053-02 一、引言 随着结构可靠性理论研究和工程设计方法的发展，近似概率设计方法已进入实用阶段，因此掌握结构可靠指标的求解方法是非常重要。其中，JC法[1]以其通俗易懂、计算速度快、适合编制计算程序和便于一般工程技术人员掌握而得到了广泛应用，且已经成为国际安全度联合委员会（JCSS）推荐采用。 然而，在水工结构的可靠度分析中，常由于随机变量变异系数过大，传统JC法求解验算点容易出现负值，而实际工程中的材料抗拉强度、抗压强度、摩擦系数、凝聚力等取值一般不小于零。另外，由于随机变量分布受到几何尺寸、物理环境等条件限制，如上游水位其最大值不可能大于坝高等，此时，一般JC法便不适用。因此，在结构可靠度分析时必须对其分布形式进行截尾分布处理[2-4]。 MATLAB已经成为了国际上最流行的科学与工程计算的软件工具，且已经成为了一种具有广泛应用前景的全新的计算机软件编程语言，具有简洁紧凑，函数丰富；程序可移植性好；图形功能强大等特点。它在科学计算、信号处理、自动控制和绘图等领域得到了广泛地应用[5]。 本文以MATLAB7.6为开发环境，编制了截尾分布下JC法的计算程序。将该程序应用于数学算例求解验证了该程序的简易、实用性，最后将该程序移植于重力坝可靠度分析中，取得了较好的效果。 二、JC法 设假设结构的极限状态方程为： Z=g（x）=0 （1） 随机变量x的均值为μx，标准差为σx，设n维度的x*=（x*1，x*2，…，x*n）T为极限状态面上一点，称 以上公式推导均是建立在随机变量服从正态分布的基础上。然而，如果随机变量不符合正态分布，如对数分布，极值I型分布等，则可以根据国际安全度联合委员会（JCSS）推荐使用的当量正态化法，事先将非正态随机变量进行当量正态化处理。 假设非正态变量xi的均值为μx ，标准差为σx ，概率密度函数为fx （xi），累积分布函数为Fx （xi），与其对应的当量正态化变量为xi′，其均值为μx ′，标准差为σx ′，概率密度函数为fx ′（xi′），累积分布函数为Fx ′（xi′），“当量正态化”需要满足的条件(参见图1)为： 三、截尾分布下的JC法 下面针对不同的截尾方式列出相关的计算公式，其推导的具体过程见文献（2）。 （一）左截尾处理时 当x*≥xp时，令x*=0.99xp，并分别用（15）和（16）计算fxi′（x*）和Fxi′（x*）。 四、截尾分布下JC法的程序编制及算例研究 （一）数学算例 对下面的极限状态函数进行可靠度分析： g（x）=3x1-x2 其中，x1~LN（30，0.12），x2服从极值I型分布，且均值为55，方差为0.1。 在无截尾的情况下采用编制的程序求得β=2.9274。对于x1进行左截尾，对x2进行右截尾，选择不同截尾点xp的计算结果见表1。其中对于x1p=22和x2p=80时的求解程序如下所示。 由表1的计算结果与文献（2）的计算结果很接近，且求解程序非常简洁，仅需要数十行即可，而FORTRAN编制需要几百行[6]，而且还需要编制不同的子程序。 （二）工程应用 采用文献[7]中的实例进行计算。混凝土重力坝断面图见图2。考虑下列随机变量：混凝土的抗拉强度Rt、抗压强度Rc、上游水位H、混凝土容重γ、坝顶附加荷载Q、建基面摩擦系数f和粘聚力C，各随机变量情况见表2。重力坝的失效模式为抗拉抗压强度失效、沿建基面滑动和整体倾覆三种。鉴于该坝为中小型重力坝，则只考虑坝体强度不够和沿建基面滑动两种失效模式。由材料力学方法[7]计算可得： 在抗滑稳定条件下，功能函数为： Z1=1000fγ-75fH+fQ+40C-5H2 在强度条件下，坝趾抗压功能函数为： Z2=Rc+6.25γ-0.147H-0.00625H3+0.0256Q 坝踵的抗拉功能函数为： Z3=Rt+43.75γ-3.903H-0.00625H3+0.0756Q 对以上三种情况，仅考虑上游水位右截尾的方式即H≥53时，用本文程序计算结果见表3。 由表3可知，在