用于转子临界转速可靠性设计的蒙特卡洛传递矩阵法.pdfVIP

用于转子临界转速可靠性设计的 蒙特卡洛-传递矩阵法 李平 核工业理化工程研究院，天津300180 摘要本文将可靠性设计理论引入到旋转机械转子的临界转造设计申，将蒙特卡洛模拟法与传递 矩阵法蛄舍起来，给出一种用于旋转机械转子临界转速可靠性设计的蒙特卡洛．传递矩阵法，谊方 法具有应用简单，计算速度快等特点，并用一算饲说明了谊方法的设计过程． 关键词转于临界转速可靠性设计蒙特卡洛模拟传递矩阵 前言 随着可靠性设计在各个领域的应用．振动可靠性设计方法也逐步发展起来“。目前，已有基 于有限元方法的成熟的可靠性设计软件．如ANSYS、NESSUS等。但是．对于大型或复杂的结构， 用有限元法进行频率分析，是很耗时的，如果再进行蒙特卡洛模拟，计算时间就可能达到无法忍 受的程度。而在旋转机械转子临界转速的计算中，传递矩阵法是一种应用较广的方法，具有简单 快速的优点。因此．将其与蒙特卡洛模拟法结合，用于旋转机械转子临界转速的可靠性设计，可 以节约大量的计算时间，在设计初期具有较大的优势。 l转子临界转速的可靠性设计 转子临界转速设计的一个重要原则是要求转子的工作转速与临界转速保持一定的距离，以避 免发生共振。因此，临界转速的可靠性设计就是通过计算转子避开共振的可靠度来进行的，从丽 可以给出转子在工作状态下，不发生共振的概率。 1．1避开共振的可靠度 按传统的振动设计规范脚，当工作转速∞和临界转速地满足以下条件时 q∞时警(1+蠡)，0磊(0．3 (1a) W 峨∞时导(1-盈)．0盈0．3 (1b) 讲 设计是安全的。但在一般情况下，工作转速和临界转速并不是定值，而是服从一定分布的随机变 量，当其分布较宽时，就可能出现工作转速和临界转速接近的情况，发生共振(见图1)。图1就是 可靠性设计常用的应力-强度干涉模型，其中的阴影部分就是可能发生共振的区域。因此，在临界 转速可靠性设计中，避开共振的可靠度定义为 吐∞对R=P((嗥一％∞)o】，，ll=l+曩 (2a) 其中，”‘和“2是传统设计中的安全系数 81 1．2避开共振可靠度的一般表达式 现以弛油例，通过应力(工作转速)一强度(临界转速)分布干涉模型来推导避开共振可靠度的 一般表达式。 设临界转速峨的分布函数为只衄)，均值和标准差分别为‰。和S。。；工作转速∞的分布函数为 G(鲫，均值和标准差为Ⅳ。和＆。从概率论可知，／／I旃照县有相同的分布，只是均值和标准差变 为n，地和nl＆，为了推导方便，用斫E替m∞，且n-锄分布函数、均值和标准差仍为G(国)、p。和。 设临界转速珏的概率密度函数为。Ⅸ砖)，工作转速确概率密度函数为g(动，两者发生干涉，见图1。 如图2所示．转速值劬落于宽度为d∞的小区间内的概率等于该小区间所决定的单元面积 Aj．即 P∽一警卜≤卜警)J吲叫一^ ㈤ 觚) 刑∞I ，(雌 “∞ 掌(01) 图I干涉模型示意圈 圈2图1阴影部分放大图 临界转速纰大于工作转速曲的概率为 P(婢q)=J。，(q)8致=A C4) 设雠与硝盘立，则g佃，)dm与j。l，(啦)d啦是两个独立的随机事件，根据两种独立事件概率的乘法