基于MATLAB的可靠度四阶矩方法及其优化设计.pdfVIP

学兔兔 第1期 (总第164期) 机 械 工 程 与 自 动 化 No．1 2011年 2月 MECHANICAL ENGINEERING AUT0MAT10N Feb． 文章编号：1672—6413(2011)O1—0027—03 基于MATLAB的可靠度四阶矩方法及其优化设计 刘彦明，张学良，杨 波 (太原科技大学，山西 太原 030024) 摘要：介绍了一种基于MATLAB计算机械可靠度的方法，当数据不充分但有足够资料写出样本四阶矩时可以 考虑此方法。并用实例给出了在计算可靠度的基础上进行机械可靠性优化设计的方法。 关键词：MATLAB；四阶矩；可靠性；优化设计 中图分类号：TP301 文献标识码：A 0 引言 法的可行性。 随着现代科学技术的发展，机械产品的结构日益 1 可靠性设计 复杂，性能参数越来越高，工作环境条件更加严酷，往 由应力分布和强度分布的干涉理论可知，可靠度 往因一个零件的失效而造成灾难性后果。有人预言，今 是 “强度大于应力的整个概率”，表示为： 后在激烈的国际市场竞争中，只有可靠性高的产品及 R(f)一P( S)一P I(S—S)0)I—P(S／s 其企业才能幸存下来。结构可靠性的定义为 ]：在规定 1) 。 的时间内和规定的条件下，结构机构完成规定功能的 其中：5为零件的强度；s为零件受到的应力。 能力。结构机构可靠度为：在规定的时间内和规定的 可靠性设计的一个目标就是计算可靠度： 条件下，结构机构完成规定功能的概率。在目前的机 r R— I／ (X)dx 。 ……………………… (1) 械设计中，往往将可靠性设计与可靠性优化设计相结 #( 0 合，使零件朝着小型化和高可靠度方向发展。日本和 其中： (X)是基本随机变量向量 一( 。，z ，-z 一， 美国在可靠性方面应用的比较广泛，从目前来看，现 )的联合概率密度，这些变量各自代表着载荷和零部 在可靠性设计_2 理论和方法已经发展到一定的水平， 件的特性；g( )为状态函数，可表示机械结构的两种 在实际工程中取得了一定的效果。但是在这种设计中 状态，g( )≥0为失败状态，g(-z)0为安全状态。将 往往将设计变量和设计参数按正态分布来计算，在其 g( )在X处展开成二阶Taylor级数 ，有： 基础上发展了一次二阶矩方法(first order and second moment，F0SM)和改进的一次二阶矩方法(advanced g—g c 十 c 一 +吉 · FOSM，AFOSM)。由于工程实际问题的复杂性和统计 数据的相对缺乏，单单按正态分布来设计会与实际问 ( 一 )( 一 ，)。 ……………………………… (2) 题不符，给计算带来大的误差，甚至有时得不到理想 这里，取二阶近似均值、一阶近似协方差、一阶近 的设计结果。对于随机变量服从非正态分布的可靠性 似三阶矩和一阶近似四阶矩，如果各随机变量相互独 优化设计情况，Rackwitz和Fiessler提出了一种等价 立，则前四阶矩为： 正态变量算法__4]，但这种方法需要知道随机变量的概 率分布函数(CDF)或者概率密度函数(PDF)。在实际 g =EE