锭子临界转速及其振型的计算.docxVIP

锭子临界转速及其振型的计算 轴的工作旋转超过第一个第一骑力轴的速度，并落在第二个骑力轴的速度范围内。为了保证轴的稳定运行，应避免旋转噪声。 锭子的静动特性计算,主要是对它们本身的临界转速及其振型的计算.临界转速的计算,实际上是计算锭子的自由振动频率.锭子在纺纱工作过程中,在套筒管上纱的质量是随时间缓慢增加的.为此我们分别计算锭子在空锭、空筒管、满纱管三种状态下的第一、二阶临界转速. 1 简化无质量弹簧系统 转子———轴承系统由若干等截面的无重轴段,集中质量,园盘和支承系统组成.集中质量为与该点相邻的两轴段的质量和的一半,因而轴段成为无质量而有刚度的轴段,把实际的支承系统简化为单质量弹簧系统,上段弹簧为油膜刚度,下段为支承刚度,两弹簧的质量为参振质量. 2 模型的传递矩阵 我们把包含有无重轴段,集中质量,园盘以及支承的段做为轴系的典型单元,从而给出右端与左端的力学关系. 由于不考虑支承系统的阻尼,所以可以只考虑YZ面的关系,(XZ面类似),第i点的状态向量以Zi表示: 其中:y为挠度,θy为转角,Mx为弯矩,Vy为切力,1为强迫项. 第i点与第i-1点的状态向量之间的关系可写成:Zi=TiZi-1;称Ti为传递矩阵,下面给出其具体表达式,我们可把典型单元分为三部分:(1)支承,(2)集中质量和园盘,(3)无重轴段. 支承右端ZR(1)和左端ZL(1)的关系为: 其中:K=Kb(Kp-Mω2),称作动刚度;Kb为油膜刚度;Kp为支承静刚度;M为参振质量.集中质量和园盘的右端Zi(R2)和左端Zi L(2)的关系为: 其中:ε=Ip-Id,M0为集中质量m与园盘质量Md的和;Ip为园盘极惯性矩;Id为园盘直径惯性矩;Ω为涡动频率;ω为自转角速度;ΔVy为不平衡量ΔVy=me·e·Ω2,me为不平衡质量,e为不平衡半径,正向涡动时Ω=ω,反向涡动时Ω=-ω,静激振时Ω=0(程序中取Ω=ω). 无重轴段右端Zi(R3)和左端Zi L(3)的关系为: 其中:α=1/EI,β=12/2EI,γ=13/6EI;E为材料的杨氏模量;I为截面惯性矩;不难看出Zi L(1)=Zi R(2),Zi L(2)=Zi R(3);另外还有Zi=Zi(R1),Zi-1=Zi L(3);所以有: 利用(5)式,就可以把轴系两端的Z0与Zn联系起来: 称T为总体传递矩阵,有了(6)式就可以进行临界转速和强迫响应的计算了. 3 纺织机械体系的临界旋转法和振动型的计算 我们对空筒管式、满筒管及空锭状态下的锭子临界转速分别进行了计算,以下是几种不同情况下的计算结果. 3.1 各机构内的位移导纳 空筒管或满筒管状态在计算时没有什么本质不同,但与空锭状态下比较,在计算处理时却有较大的不同,关键在于筒管和键杆组装在一起,相互之间产生约束,在处理时带来一些麻烦.如图1所示: 在处理这种组合系统时我们采用机械阻抗的方法,将系统分成两个子系统即空锭杆上下轴承为第一子系统,筒管加纱卷为第二子系统,并在A、B两点组合.计算时假定AB两点为锭支计算步骤如下: 1)先算第一子系统:在A点加一单位激发力用TRP程序的不平衡响应计算部分,算出各个不同转速(ω)下A点位移导纳H1AA(ω),B点的位移导纳H1AB(ω). 2)用同样的方法在B点加一单位激发力算出各个不同转速(ω)下A点位移导纳H1BA(ω),B点位移导纳 3)用类似办法算出第二子系统的位移导纳HΠAA,HΠAB,HΠBA,HΠBB,它们也都是ω的函数.若在第一子系统A、B两点上有fA、fB激发力作用.由线性振动理论,可求出A、B两点的位移dAⅠ,dBⅠ,并由连接点A、B的位移协调和力的平衡条件,可得: dAⅠ=dAΠ;dBⅠ=dBΠ,fAⅠ+fAΠ=0;fBⅠ+fBΠ=0, (7)(8) 因此可得到: AHⅠAB+HΠABHⅠBB+HΠBB HⅠAA+HΠAAHⅠBA+HΠBA AfⅠBfⅠAA=0. (9) 要得到fAⅠ、fBⅠ的非零介,必须满足其导纳矩阵的行列式为0,即得到频率方程: (HⅠAA+HΠAA)(HⅠBB+HΠBB)-(HⅠAB+HΠAB)(HⅠBA+HΠBA)=0. (10) 4)解算频率方程(10)求出临界转速. 5)计算振型,先求导纳矩阵的伴随矩阵: dj(H)= -(HⅠBA+HΠBA)A.A-(HⅠAB+HΠAB)HⅠAA+HΠAA HⅠBB+HΠBB (11) 导纳以临界转速ωi值代入算出.伴随矩阵任一非零列为对应于临界转速ωi的特征向量. 故第一第二两个子系统在分别解除约束后只需在相连点A、B加上由(8)外力fA、fB,计算在临界转速ωi=ωi(i=1,2,…)下的强迫振动. 3.2 第二阶弯曲振型 用我们研制的计算机程序,分别对纺纱锭子组合系统各种状态下的临界转速及振型进行了大量计算,其计算结