matlab用jc法计算可靠度,基于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可知,在考虑截尾分布的情况下可靠指标均有一定程度的提高,其主要是因为截去了由于水位过高导致重力坝的可靠指标降低的不安全部分。当然,这也只有在保证满足截尾分布的要求下,如不产生漫坝的情况,才是真正的计算结果。因此在采用截尾分布处理时,首先需要了解截尾分布的实用性和有效性。
五、结论
本文充分利用MATLAB软件编制程序的语言紧凑,使用方便灵活,以及利用丰富的库函数避开了繁杂的子程序编程任务,成功实现了截尾分布下JC法求解可靠度。
通过数学算例验证了截尾分布下JC法的可行性和实用性,并且通过重力坝的应用可知,对于一些实际工程中截去不可能出现的情况时,其可靠指标的求解更符合实际情况。但是,对于一些难以确定截尾点数值的变量,如摩擦系数、凝聚力等,还需更加深入的研究。
参考文献:
[1] 赵国藩.工程结构可靠性理论与应用[M].大连:大连理工大学出版社,1996.
[2] 吴世伟,李宏.变量分布截尾下可靠度计算的JC法[J].河海大学学报,1988,16(6):75-80.
[3] 许福友,陈艾荣.基于截尾概率分布的结构可靠性分析[J].工程力学,2006,23(11):52-57.
[4] 林明财.截尾JC法在拱座稳定分析中的应用研究[J].水利科技,2007,3:40-42,46.
[5] 张志涌.精通MATLAB6.5版教程[M].北京:北京航空航天大学出版社,2003.
[6] 王法银,程心恕.JC法求重力坝可靠度的电算程序[J].农田水利与小水电,1995,3:22-27.
[7] 吴世伟.结构可靠度分析[M].北京:人民交通出版社,1990.
matlab用jc法计算可靠度,基于MATLAB的截尾分布下JC法计算可靠度相关推荐
- 基于matlab 的燃油喷雾图像处理方法,基于MATLAB的燃油喷雾图像处理方法.doc
基于MATLAB的燃油喷雾图像处理方法 基于MATLAB的燃油喷雾图像处理方法 摘要:提出了基于MATLAB的燃油喷雾图像处理方法.通过对喷雾图像的采集, 分割和滤波处理, 将多幅图像转换融合为一幅 ...
- matlab x(n)16点DFT,[基于MATLAB的数字信号处理实例分析.doc
[基于MATLAB的数字信号处理实例分析 湖北文理学院理工学院 学生结业论文 课程名称:MATLAB教程 结业论文名称:基于MATLAB的数字信号处理实例分析 专业名称:通信工程 班级:1011 学号 ...
- matlab各种潮流计算代码,基于MATLAB的潮流计算源程序代码
<基于MATLAB的潮流计算源程序代码>由会员分享,可在线阅读,更多相关<基于MATLAB的潮流计算源程序代码(7页珍藏版)>请在人人文库网上搜索. 1.*电力系统直角坐标系下 ...
- matlab计算表面粗糙度,基于MATLAB图像处理的铸造表面粗糙度测量方法
铸造表面粗糙度严重影响铸件的外观质量和使用性能,表面粗糙度还影响机械设备配合性质的稳定性和动力的消耗,而且表面越粗糙,越容易产生噪声.振动等[1-2].所以,控制和提高铸件表面粗糙度,不仅能够有效提高 ...
- MATLAB计算杨氏模量,基于MATLAB处理动态图像的杨氏模量测量方法.doc
7A版优质实用文档 PAGE PAGE 16 7A版优质实用文档 基于MATLAB处理动态图像的杨氏模量测 量方法 黄浩钊1,陈玥琦2,李家卉3 第三军医大学学员13营,重庆,400038:2.第三军 ...
- matlab滤波器设计双陷滤波,基于MATLAB 双线性变换法IIR 滤波器的设计
摘要:滤波是信号处理的基础,因此滤波器的设计也就成为数字信号处理的基本问题之一.数字滤波是对数字信号处理不可或缺的环节,其好处远远大于模拟滤波器.本文应用MATLAB,实现了用双线性变换法设计IIR滤 ...
- matlab为什么要升维数,基于MATLAB和升维投影法的手写字符识别输入系统
256 现代交际·2019年6期 作者简介:王思达,长春师范大学学生,研究方向:计算机科学与技术:刘勇,长春师范大学学生,研究方向:计算机科学与技术:袁汇灵, 长春师范大学学生,研究方向:计算机科学与 ...
- matlab计算正交变换,基于Matlab的LFM信号的正交变换和脉冲压缩
基于 Matlab的 LFM 信号的正交变换和脉冲压缩 付银娟 (西安邮电学院 陕西 西安 710121) 摘 要 :正交变换和脉冲压缩是雷达信号处理中常用的两个基本技术.介绍了正交变换和脉冲压缩的基 ...
- 自控matlab设计,自动控制原理课程设计--基于MATLAB软件的自动控制系统仿真
自动控制原理课程设计--基于MATLAB软件的自动控制系统仿真 < 自动控制原理 >课程设计说 明 书日期: 2012.5.21-2012.6.1 自动控制原理 课程设计评阅书题目 基于 ...
最新文章
- 宁波大学计算机专业复试,2016年宁波大学信息科学与工程学院计算机专业考研复试题库. (1)...
- 《资本说》如何帮助企业从“0”到“1”,看王波玩转投资攻略
- 并非所有图像都值16x16个词--- 清华华为提出一种自适应序列长度的动态ViT
- 安全狗:云时代的服务器安全守护者
- Python多模块文件共享变量
- 如何配置jenkins 与代理服务器吗?
- 082_html5Web存储
- 解读设计模式----单例模式(Singleton Pattern)
- java字符串怎么拼接字符串_Java中String使用+ 拼接字符串的原理是什么?
- 通过binlog恢复mysql数据库
- 国家码信道顺从表(含数据库和配置)
- c语言循环教案,C语言教学(七-上)for循环
- Pytorch深度学习笔记(02)--损失函数、交叉熵、过拟合与欠拟合
- SDN多控制器资源池的设计与实现
- arm linux打印机驱动编译方法,嵌入式linux驱动打印机HP系列
- ATMEL芯片系列介绍2
- 养蚕日记软件测试,养蚕的过程日记300字
- iOS视图成像理论及性能优化
- Java SE 基础部分经典100道笔试题
- vue qrcodejs2生成二维码实现手机APP扫码进行web网页登录
热门文章
- Android 框架学习2:源码分析 EventBus 3.0 如何实现事件总线
- 第12章 智能优化算法(ACO)
- CTF-综合测试(高难度)【超详细】
- php 公众号 群发,微信公众号【服务号】群发策略调整,并开通高级群发接口
- [PTA]练习5-3 数字金字塔
- php laravel 教程,Laravel框架学习之新手教程
- 基于华视身份证读卡器读取身份证信息的Android demo
- Tampermonkey扩展下载
- 服务器U盘安装虚拟化,用Proxmox ISO镜像制作引导U盘
- ffmpeg AV转换常用命令