收藏 分享(赏)

高等转子动力学.doc

上传人:wspkg9802 文档编号:6258532 上传时间:2019-04-03 格式:DOC 页数:21 大小:1.18MB
下载 相关 举报
高等转子动力学.doc_第1页
第1页 / 共21页
高等转子动力学.doc_第2页
第2页 / 共21页
高等转子动力学.doc_第3页
第3页 / 共21页
高等转子动力学.doc_第4页
第4页 / 共21页
高等转子动力学.doc_第5页
第5页 / 共21页
点击查看更多>>
资源描述

1、高等转子动力学理论、技术与应用闻邦春 顾家柳 夏松波 王正 1999.8 机械工业出版社读 书 笔 记Areo Chan目录 第 1 章 转子动力学的计算和分析方法 .1油膜力 .2各向同性转子系统计算分析 .3转子临界转速和振型计算 .4转子系统的动不平衡响应 .5各向异性支承转子系统的计算分析 .5传递矩阵法求转子系统的瞬态响应 .7传递矩阵法阻抗耦合法 .8传递矩阵分振型综合法 .9传递矩阵直接积分法 .10关于盘、轴单元的传递矩阵及数值积分差分格式的讨论 .11 第 2 章 轴承的动力特性 .15固定瓦径向滑动轴承的油膜刚度和阻尼系数 .15油膜刚度和阻尼系数的表达式 .16一般径向滑

2、动轴承的油膜刚度系数和油膜阻尼系数的计算 .17有限元法 .17 第 3 章 转子系统的动力稳定性 .19转子动力学读书笔记1 第 1 章 转子动力学的计算和分析方法主要内容:转子弯曲振动的临界转速,不平衡响应,稳定性,各种激励下的瞬态响应计算,扭转振动的固有频率和响应。MzCKF系数为系统的质量、阻尼和刚度矩阵,z 为广义坐标矢量,F 是广义外力转子回转效应陀螺矩阵结构动力学 ()()GzSz其中 G 为反对成阵,陀螺矩阵, K 是刚度矩阵对称部分, S 为不对称部分。复杂多转子、多支承的自由振动微分方成形式: 00iMZCZ求解转子动力学的两种方法:传递矩阵法和有限元法转子本体刚度建模,周

3、断弯曲刚度的简化比较复杂。轴承建模(油膜力,弹性阻尼支承):在轴径围绕静态平衡位置小幅涡动时: 0xxyyRxKCyK,C 分别 4 个元素,刚度系数和阻尼系数,八个系数为油膜动力特性系数。滚动轴承一般模化为弹性支承,阻尼很小若轴承座弹性不能忽略,可简化为一个有质量、阻尼和弹簧组成的单自由度系统。当刚度和转子相似或更低,可以作为整体考虑两种模型:质量沿轴线分布的分布参数模型和质量简化到 n 个节点上的有限自由度的集总参数模型(模态的截断) 15.34Nr为要求计算的固有频率(或临界转速的)最高阶次r根据清华大学出版社转子动力学集总化结点 i 的运动微分方程 LRibjibj xxij j yy

4、i iixxFtQmyyKC0()RLp xxdi yyi ii iiI tMI转子动力学读书笔记211111122LRxi i iiyi i iiiiiiiiiQxxyyxy 1111 112323RLxi i iiy i i iiir iiri iiiMxxyyl 其中 , ,321iiEJl26iiEJl3iiEJl支承 j 的运动微分方程 ibjibjb bbjj j jjj jjxxxmyyy KCKC和 分别是作用在结点 i 上的外力和外力矩。FL油膜力轴径上的线性油膜力: xbbbyRxymCK1122xy tttDE式中:, , , ,112bCAa1,2btbttEAbmI1

5、2taCK,2bbbttmCK21bbttD 21124114bbt t t t tb bb bt xxxxyyyytt CCmt转子动力学读书笔记3,21124b bb bbt t tbbtxxyyttxy mmECC各向同性转子系统计算分析沿轴线一个平面内的弯曲振动模态临界转速和相应振型右手坐标系 ,原点在转子左端, 沿轴线,从 端部看,转子逆时针旋转。oxyzozz转子第 i 个截面:状态矢量为 ,由界面径向位移 、挠角 ,弯矩 和剪力 的幅值组成。iZixiiMiQTi iZXAMQ1ii(1)带弹性支承的刚性薄片圆盘:第 j 个支承的总刚度2()bsjKmK 是刚度,m 是质量,b

6、表示轴承座, 为涡动角速度,I 为转动惯量(p 极轴,d 直径)sjKRiLiiQRiipdiI ixi圆盘惯性力: 2imX惯性力矩: ()dipiiIA由达朗伯原理(理论力学) 22()RLiiiisjiiidipiiQKMI则传递矩阵 221001()i dipiisjDImK 转子动力学读书笔记4无质量等截面弹性轴段 ,其中2321(1)601ilEJlB26tEJkGAl转子临界转速和振型计算Prohl 传递矩阵法随着试算频率的提高,运算精度会降低尤其是大型轴系。原因是传递矩阵中大部分的项含有 ,若 很大,最后将导致计算两个大数的差值,导致计算精度下降。表现2N为轴尾部幅值急剧增加。

7、书中给定模型两端自由,有 ,0QM把状态向量 r 个元素分成零值和非零值部分,有 , , 有,TiifQ,TiieXA,TiiZfe12i iuffee其中元素 u 来自于21iBD引入 Riccati 变换 :iifse121122i iieusf u 所以 11i iis初始条件: 10,0fe末端条件: (系统频率方程式)1211 1, 0NNNsf利用频率扫描法求解比例解 ,各截面状态矢量的比例解,即各截面位移的比例解,就是对应此临界转速的振1AX型。 1210NsA转子动力学读书笔记51211112NNNssAXXNNeX转子系统的动不平衡响应在薄圆盘上加不平衡力 , 位该节点圆盘具

8、有的不平衡质量矩。则圆盘两端的状态2iiFUi矢量传递关系为: 22210001()R Ldipiisji iiXXAAIMMmKQQU 支承刚度增加,临界转速必然增高。在支承刚度的某些范围内,临界转速的变化很明显,有些范围内则不明显,变化比较平缓。对于多跨轴系,总体的临界转速和振型与单跨轴系存在一定的关系,即每阶振型都有一个转子主导。各向异性支承转子系统的计算分析需考虑转子铅垂和水平两个平面的运动耦合及支承阻尼作用。振动量有相位差,用复数表示。振动量 复振幅 复频率cos()tyYecsYiSi分别为余弦分量、正弦分量, 衰减指数、阻尼圆频率,cs ,振幅 相位角2csarctnsY常将 写

9、成 Y,或者 y。 2SRiccati 传递矩阵: TTi xyi ifeMQXAYB Z总支承刚度 212 S bbbSSSKCKCmKCK 油膜刚度系数矩阵,C 阻尼系数矩阵,下标 b 表示基座。质量矩阵 bxym若不考虑交叉刚度系数、交叉阻尼系数和系统阻尼时,有 , ,Si0sxysxK0b传递矩阵分块形式为:转子动力学读书笔记611i ilul2212 22()00()sxdsxypisyxpsydilmSKISlKISullm 23221 232616i ilEJlulEJl 3 23222 2223 3 222 2(1) (1)6 61(1) (1)661sxdsxypsx sxy

10、i syxp sydsyx syl lllmSKISKISJJEJEJEEul ll lISmSISJJJJKKEE 2i (系统频率方程式) ,为复数方程式,求解方法有121342N+341243N+1s=0sNewton-Raphson 法(切线法)和 Muller 法(抛物线法)转子系统复模态的特点:各结点轨迹为形状、方位不同的椭圆;转子轴线是一空间曲线;椭圆短半轴 b 为正,则正进动,正向涡动,小于零则反进动,反向涡动。存在混合进动(有正有负) ; 为负,运动是衰减的,为正是失稳的,为零是临界的。不为零时,结点涡动轨迹是不闭合的椭圆形螺旋线转子系统动不平衡响应计算即增加两个离心力分量

11、和 ,在薄圆盘上,两端截面状态矢量关系:2U2i转子动力学读书笔记72222100001001Rx xdpsxsxyy ypdsyx syiiM MIiIQ QmKKimX XA AY YB B 200LiiUiRiccati 变换 iiifSeP递推公式 11122i iiuuS1i fii+eiPFSF112122i i ieieeP其中 2fiiilUlF22()31()ei iillUlEJF传递矩阵法求转子系统的瞬态响应Riccati 传递矩阵法和 Newmark 法:根据 Newmark-,取 为 0.5,有22111t ttt t qqqtttt为广义坐标,若某瞬时,已知某结点的

12、位移 、速度 、加速度,并知 时的位移。q qt时刻构件瞬态传递矩阵:t12fei iiuffeF11ilul12131223221()00()lCKlDKull 转子动力学读书笔记8231231u312312321212 312()()CKlDKu l 1122()fAlEBlF2131212()()eAEBF其中 23231 322 2232()()11()t tt xytKKA Lt t 122 ()1xytt tt FtBxmmy 、 、 、 、 和 来自线性油膜力公式1C1D21E212Kt2dIt3pIKt1lEJ2lJ3316lEJ传递矩阵法阻抗耦合法基本思路:将复杂的转子、支承

13、系统分成若干子系统,对简单的子系统采用传递矩阵法建立方程,各子系统之间利用分割点处机械阻抗耦合关系形成一组齐次方程,从而求解整个系统的固有频率和临界转速。根据机械阻抗原理,系统在分割点处的机械阻抗总和为零。因此要分别计算系统在分割点处各子系统的机械阻抗。转子系统的机械阻抗特性图中 B 为特殊支承,将系统在 B 处分割成两个子系统。在传递矩阵方法中, 将起始端处理为自由端(若实际上并非自由端,则添加一段自由轴端,方便计算) ,即 和 。在分割0M0Q点处 和 , 为机械阻抗。0nMnnQZX转子子系统传递矩阵关系:转子动力学读书笔记90nXTMQ4123/nttZOABnOABC如果图中分割点在

14、中部,则将系统分割成两个子系统。在分割点处,还必须满足变形一致与弯矩现等的条件 ABCBACBXMQZ机械阻抗的一般形式为 2sKmiK 为静刚度,Zs 为机械阻抗,又称动刚度。传递矩阵分振型综合法将传递矩阵法与模态综合法相结合。基本思路同样是分割子系统,然后在分割点处约束子系统,用传递矩阵法计算约束振动模态与约束静态位移,然后进行模态综合。分振型综合法将整个系统的广义坐标分为内部坐标 和外部坐标 两部分(边界坐标代表全部分割点)iZb改写 为矩阵形式,其中各量按内部坐标和外部坐标分块。得约束00iMZGCK子系统一般方程 ,若假设约束子系统无阻尼,则0ii iii,其特征解为 ,其中 为约束

15、振动模态振型, 为0ii ii citeZc约束振动模态角频率。引入坐标变换 ,其中 为轮流释放约束边界坐标时,约束子系统内部坐标12ibZqI得静位移模态矩阵。引入 12stQeq得到复杂转子支承系统自由振动运动方程模态综合表达式转子动力学读书笔记101121220ABQ其系数矩阵便是频率方程。式中: ; ;210TciciAMG2bsbSACK; ;21 0TTiiSSMG2TiiSB;2iiB 200TiibbiG模态截阶处理复杂的转子支承系统,各转子系统的自由度很多,没有必要全部考虑。每个约束转子子系统仅保留若干必要的低阶约束振动模态。传递矩阵直接积分法对复杂的转子支承机匣系统,选定若

16、干各特征盘,对每条链的特征盘借传递矩阵建立运动方程,此即系统的运动方程式。对其次方程的求解,可以得到系统的固有频率和振型,对非其次方程求解,根据不平衡响应的峰值位置,可确定系统的临界转速;根据特征盘盘心轨迹,可确定系统是否失稳;如果在积分过程中采用瞬态传递矩阵,就可以得到系统的瞬态响应,所以这是一种通用的转子动力学计算方法。特征盘的选取,主要考虑转子主要模态的能观性和最易碰撞的部位。设,由链左端到第 i 个特征盘的传递矩阵为 ,从 i 个特征盘后到右端的传递矩阵为LiTRiT则不计耦合力时运动方程 11220RLdp xxITm1122RLdp yy联结结构的动柔度矩阵 12g第一列是由单位激

17、振力引起的相对位移和相对转角,第二列是单位激振力矩引起的相对位移和相对转角。该联结的动刚度为 1Kg例如处在不同链上的特征盘 j 和 k,联结的结构特征 ,则 j 盘受到来自 k 盘的耦合力K, kjjkxR特征盘的运动方程总特征盘 N 个。系统有 2N 个特征盘动力学方程,每个有 2 个代数方程式,故方程自由度数为转子动力学读书笔记114N 个。对于特征盘 i,到盘左侧传递矩阵 ,耦合力 ,它到该链右端的传递矩阵 ,耦合力 ,iALRiBR并记, 或 ,dpi iSxImiixedpi iSyImiiye系统运动微分方程: 111 1222, 2, 2,0RLRRLLi iiij ikiik

18、j kT TeTBAR受不平衡力作用的系统运动方程: 111 1222, 2, 2, ,12, 0RLRRLLi iiij ikiikj kRLLiiminiinimTC SeTDFF在 xoz 平面, ,20cosmUtF20cosnUt2cosiiUt在 yoz 平面, ,2sinmt2sinntF20siniitF临界转速的物理意义: 转子同步正进动的固有频率;转子不平衡产生的同步响应达到峰值时的转速xMxQnxMnx 1nxnln1n1nxM1nxQyxT i m o s h e n k o 梁轴有限元与轴单元关于盘、轴单元的传递矩阵及数值积分差分格式的讨论不计外力和惯性力,由微段平衡

19、条件可得 0xzxxQz转子动力学读书笔记12由 Timoshenko 梁假设2bxMEJz1xQzkAG是弯曲变形产生的挠度, 是剪力产生的剪切形变。b 1x轴单元传递矩阵: 2321/6/11tx xi ilEJllkAJMMlQQ 梁轴有限元形函数方程 01231 23060i teii txaEJkAGll q c两者都不计惯性力并采用 Timoshenko 梁假设,有限元法中考虑均匀分布的质量及转动惯量,而传递矩阵法中,只能简单的把轴元单元的质量和转动惯量分成两半,分别加于两端的圆盘上。改进的轴元单元传递矩阵不仅采用 Timoshenko 梁假设,还考虑惯性力的影响,并在集总参数时,

20、利用实际存在的折合关系。计惯性力,由微段平衡条件可得又有:xQAzxxMzsxtQzkAG得到统一方程式: ( )42()()()0tdXdTTXkEJ420dXabz计算稳态不平衡响应与临界转速时,转子作同步涡动 0cos)tt计算稳定性分析及瞬态响应时,略去 内进动速度变化,t()T则轴单元传递矩阵,即: 1 10010012212000tt tx xi itaEJJkAGkAGkTM MaEJQ Qak 转子动力学读书笔记13其中 中元素:T112cosh()cos()talal123142sinh()si()tall356 478 21122sin()sin()tll231142cos

21、()cos()tlal356haa478ha2211cos()cos()tll221314sin()sin()tll3562 4782 a3411sinh()sin()talal332114cosh()cos()tall3562 478221ac12c231ac214ac56712821ca1ca2b轴单元的功能为 220001()()2()l l liTAxydzJdzJdz 盘的动能为 2211()()()ddmII盘的瞬态传递矩阵基本关系:RLgRLixMQF由盘的动能方程导出陀螺力矩 ,gxdpyII又转子动力学读书笔记14cossinxryr则盘的瞬态传递矩阵 36324541727

22、6818500100R Lx xy yMMTTQQTT 其中, 232pdI36 2pdII241rTm45rTm7236763281458541T需要参量可以从轴单元传递矩阵方程式得到 1 0010011222000tt tx xi itaEJJkAGkAGkTM MaEJQ Qak 数值积分的差分格式关于数值积分方法,工程界常用的有 Houbolt 法、Newmark 法、尤拉后差法,对于转子系统,希望步长 t 大于一定数值,以保持数值稳定性。尤拉后差法稳定条件为 ,Houbolt 法稳定条件为20.1iipt0.53ipt:采用数值积分求瞬态响应时,尤拉后差法允许较小的步长,由于瞬态响应

23、分析总希望步长较小以充分显示细节,因此尤拉后差法更适宜与瞬态响应分析。Newmark 法的稳定性次于 Houbolt法,因此 Newmark 法求瞬态响应是不合适的。如果求系统的加速响应,可以用 Houbolt 法采用大步长以节省计算时间,精度优于尤拉后差法。对于一般结构,Newmark 法与 Houbolt 法具有同样的精度,但前者算法简单。若作转子稳定性分析,根据 Lyapunov 稳定性定义,最关心的是周期解受扰后瞬间内响应细节,用尤拉后差法为宜。转子动力学读书笔记15 第 2 章 轴承的动力特性通常是给定工况、轴承结构参数、润滑油特性等条件下,通过求解雷诺方程,确定系统静平衡位置时的油

24、膜压力场,并由此求得油膜力。将油膜力视为平衡点附近位移和速度的函数,油膜力增量的线性表达式为: xxyxyyFkc式中 8 个油膜力动力特性系数。固定瓦径向滑动轴承的油膜刚度和阻尼系数 WFxeyro径向滑动轴承简图径向滑动轴承非定常运动雷诺方程为: 3321612coshphphyxinrz瓦面为圆弧形,用 和 表示速度扰动,上式变为e3321612coshphpheinrz其中 为油膜厚度, 为油膜压力, 为润滑油动力粘度, 为轴向坐标。取 ,zhcH, , , , , , 。其中符号002pLzxcyec“ ”表求导, 为轴承半径间隙, 为轴承有效长度, 为轴径转动角速度, 为对应于进c

25、 0油温度和压力的动力粘度, 为偏心率。则圆弧瓦面雷诺运动方程写成无量纲形式为:转子动力学读书笔记16233126cosHpDHpsinL式中 为轴承直径。对于短轴承,非定常雷诺运动方程为 23126cossinL油膜刚度和阻尼系数的表达式 /2sinbaLxxFkprdz /2sinbaLxxyFkprdz/2cobaLyyx /2cobaLyy/2sinbaLxxFcprdz /2sinbaLxxyFcprdz/2cobaLyyx /2cobaLyy无限长轴承、有限长轴承和短轴承各自的油膜交叉阻尼系数相同;对于无限长轴承(Sommerfeld 边界条件)只有油膜交叉刚度和正阻尼系数,用一对

26、圆轴支承的刚性转子是不稳定的。短轴承油膜刚度和阻尼系数的表达式(用于计算) 222416xLkD22245/161xyk22245/3yxLkD 222434161yk 2223/ 261xLcD转子动力学读书笔记1722241616yLcD 2228161xy一般径向滑动轴承的油膜刚度系数和油膜阻尼系数的计算常用方法有差分法、偏导数法、小参数法及有限元法差分法基于 等,在平衡位置取 和 无量纲位移扰动,取 和xxFkxyx无量纲速度扰动。由几何关系计算新位置下偏心率和偏位角、无量纲油膜厚度 ,带入y 1H雷诺方程求解无量纲油膜压力 ,可以积分得到 和 (e.g. ) ,则响应的油膜刚度系数可

27、1p(1)xF(2)以按如下计算 。需要计算 8 次雷诺方程,若采用 和 (最后需要坐标变换)(1)(2)xFk,则只需计算 6 次雷诺方程。偏导数法是根据定义 (无量纲表达式) , 油膜破裂边界1sinbaxxkpd b随扰动而变化,经推导,得 ,只要计算出 、 、 和1ibaxxpxyx就可以求出 8 个所求系数,只需要计算 4 次雷诺方程且可避免差分代替导数产生得误差。py有限元法计算速度快,只需一次计算雷诺方程。m o e5 ( 1 , 1 )4 ( 1 , 1 )3 ( 1 , 1 )2 ( 1 , 1 )1 ( 1 , 1 )标准正方形单元设 为常数,静态雷诺方程为转子动力学读书笔

28、记18233pDpHHL由变分原理知,上式与下式泛函的极值等价 223() 6Jpdpd 采用自然坐标系下正方形标准单元,有:, , ,1,miik1,miikH1,miik1,iik其中, 为第 i 个节点的压力, 反映第 i 个节点影响的形变函数,m 为单元节点数, 为单ipi ik位布尔矩阵,k1,2,n,为节点总数。与 的关系:(,),12J式中, , , ,11miiJ121miiJ211miiJ21miiJ为 Jacobi 行列式,记做, , ,12kmk12emp1kmk1kmk若取 、 为常数,则有: ,3Hsin3sineeHCpWsisineJd其中 是单元刚度矩阵。3C总刚度矩阵和总线性方程可写为: sin3Kp式中 , , 。将边界条件带入,即可求得压31()neeKHsinsin1eW12,Tnp转子动力学读书笔记19力 p油膜力分量 1sinsin3TxFWK1cossin3TyFWK按非定常雷诺方程,设 为常数,系统的总刚度矩阵和总线性方程为sinsincos6pxy 第 3 章 转子系统的动力稳定性转子运动失稳的根本机理:扰动运动时受到一个切向力。在线性系统理论中,与转子扰动位移成正比与正交的切向力称为循环力,与扰动速度成正比与正交的力称为陀螺力稳定性理论和稳定性判据 0qMCK转子动力学读书笔记20

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 生活休闲 > 社会民生

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报