机械振动大作业

机械振动项目2017-5-21沈阳航空航天大学学院：创新创业学院班级：1434200101班组员1：邵增德学号：143406010621组员2：哈哲远学号：1434070103301目录一、significanceoftheproject............................................................2二、MatrixIterationmethod........................................................22.1方法的一般性说明或介绍....................................................22.2方法的应用说明....................................................................32.3计算结果...............................................................................5三、Transfermatrixmethod..........................................................63.1方法的一般性说明或介绍......................................................63.2方法的应用说明....................................................................73.3计算结果..................................................................................8四、近似计算结果的验证..................................................................9五、两种近似方法及其计算结果的比较分析.................................10附录.......................................................................................................11附录ATheMatlabprogramofmatrixiterationmethod...............11附录BTheMatlabprogramoftransfermatrixmethod................122一、significanceoftheprojectAeroenginemodernpursuitofhighperformanceandhighthrustweightratio,thestructureismorecomplex,moreandmoreharshworkingconditions,whichleadstotheexcessivevibrationfactorsgraduallyincreased.Vibrationisthedecisivefactorthataffectsthelifeofaero-engineandflightsafety.Thefollowingfigureshowsthemainshaftofanaero-enginecompressorandturbine,whichissimplifiedasamultidegreeoffreedomvibrationsystem.Inordertoensureitssafeoperation,itisnecessarytodetermineitsmainvibrationmodes,naturalfrequenciesandmainvibrationlaws.二、MatrixIterationmethod2.1方法的一般性说明或介绍系统的主振型方程可分别表示为}]{[][}{12AKMA（1）}]{[}{12AMA。（2）矩阵迭代法是从假设主振型出发，对以上两式进行矩阵迭代运算。对（1）式进行迭代时，是从系统的最高阶固有频率与主振型开始，逐步求得系统的全部或一部分固有频率和主振型。而对于（2）进行迭代时，则以系统的最低阶固有频率与主振型开始，逐步求得系统的全部或一部分固有频率和主振型。因为工程上通常对系统的最低阶或较低的几阶固有频率和主振型比较重视，所以下面只考虑的（2）式的迭代运算。引进动力矩阵]][[][MD,(3)则（2）变为aaaaaa6aa3}]{[}{12ADA(4)当（4）式中的{A}取为系统的各阶主振型}{)(iA时，该式成为}]{[}{12ADA(5)当（5）式中的{A}取为系统各阶主振型}{)(iA时，该式成为}]{[}{1)()(2iiADA（i=1,2...，r）（6）随着迭代次数的增加，{𝐴}𝐾越来越接近于第一阶主阵型，即,{𝐴}𝐾={𝐴}𝐾−1，因此我们将得到关于固有频率的关系式，从而计算出n阶固有频率。2.2方法的应用说明1、本方法为矩阵迭代法，所需要的公式如下所示：转动惯量𝐽=12𝑚𝑅2刚度𝐾=G∗IpL惯性矩Ip=π𝐷4322、根据所给条件，进行编程计算刚度k和转动惯量𝐽，为了计算的方便，程序如下：//D1=0.4;h=0.04;G=7.69e10;rou=7800;D2=0.04;a=0.12;V=pi*(D1^2)*0.25*h;m=rou*V;j=0.5*m*(D1^2)*0.25;Ip=pi*D2^4*(1/32);k=(G*Ip)/a;//3、对系统各段的转动惯量以及刚度进行赋值并分别列出对应的矩阵，具体程序如下所示：//J1=j;J2=j;J3=j;J4=j;J5=j;J6=j;J7=j;J8=j;k1=k;k2=k;k3=k;k4=k;k5=k;k6=k;k7=k/6;k8=k;M=[J10000000;0J2000000;00J300000;000J40000;0000J5000;400000J600;000000J70;0000000J8];K=[k2-k2000000;-k2k2+k3-k300000;0-k3k3+k4-k40000;00-k4k4+k5-k5000;000-k5k6+k5-k600;0000-k6k6+k7-k70;00000-k7k7+k8-k8;000000-k8k8];//4、然后进行编写矩阵迭代方法的程序，首先根据我们的编程思路画出流程框图如下图2-1所示图2-1具体程序如下所示：/*/K=K+alf*M;D=inv(K)*M;forj=1:nA=[10000000]';fori=1:10000B=D*A;a=B(n);B=B/a;ifmax(abs(B-A))0.00000001om(j)=1/sqrt(a);om(j)=sqrt(abs(om(j)^2-alf));break5endA=B;end/*/5、最后计算出结果，并画出主阵型D=D*(eye(n)-B*B'*M/(B'*M*B));MM(:,j)=B;endMMMMF=[00000000;MM;00000000];plot(MMF)2.3计算结果1、根据所编写的程序计算出八阶固有频率如下图2-2所示图2-22、作出的八阶主阵型图如下图2-3所示：图2-36三、Transfermatrixmethod3.1方法的一般性说明或介绍对于盘轴扭振系统，可以将该整体系统分割成结构形式相同的单元，即一段轴连接一个圆盘，将盘视为刚性质量元件、轴段视为柔性无质量元件，每个单元的扭转运动取决于轴段梁福安的转角和扭矩，只有这两个主要因素，因此取各元件的转角和扭矩作为状态变量，表示为[],TXT，对整个盘轴扭振系统可以理解为状态变量是通过各个单元从系统的左侧传递到系统的右侧，我们通过对一个单元进行分析，得出各个单元的传递矩阵，进而可以求出整个系统的传递矩阵，然后再根据边界条件就可以求得关于2的方程，由于方程过于复杂，我们可以将该表达式的图形作出来，找出使该表达式为0的点，该点即是我们要求的值，从而求出该系统的各阶固有频率。取单元中的第i个单元来来进行分析：（1）：对刚性质量元件盘进行分析可得：2RLiiRLiiiiTTJ，写成矩阵的形式为：111/01LRPiiiikSTT,该矩阵称为点传递矩阵。（2）：对柔性无质量元件轴进行分析可得：111()LRRiiiLRiikTTT，写成矩阵的形式为：2101RLFiiiiSTJT，该矩阵称为场传递矩阵。于是得到第i单元两端状态的传递矩阵（从第i-1个盘右侧到第i个盘右侧）22111/(1)RRPFiiiiiikSSSJTTJk，该矩阵为单元传递矩阵。对于有n个盘的轴系，将各个单元的单元传递矩阵反顺序连乘就可得到整个系统两端状态的矩阵形式：1100RRRiiiXSSSXSX,即：22011122202122()()()()RRnRRnSSTTSS再根据轴系两边的边界条件，在该系统中00LRnTT，可以得到关于2的方程，7对该关系式作图求出与x轴的交点就可以得到系统的固有频率。3.2方法的应用说明1、本方法所使用的的数据和公式同2.2.1步骤，赋值程序也是相同的，故这部分不再陈列出来。2、对系统各段的转动惯量以及刚度进行赋值并分别列出对应的矩阵，具体程序如下所示：J1=j;J2=j;J3=j;J4=j;J5=j;J6=j;J7=j;J8=j;k1=k;k2=k;k3=k;k4=k;k5=k;k6=k;k7=k/6;k8=k;J=[J1J2J3J4J5J6J7J8];k=[k1k2k3k4k5k6k7k8];3、然后进行编写传递矩阵方法的程序，首先根据我们的编程思路画出流程框图如下图2-4所示图2-4具体程序如下所示：step=0.1;mn=1;T21p=0;fori=1:9000om=(i-1)*step;T=eye(2);8omx(i)=om;z(i)=0;forj=1:nT=[11/k(j);-J(j)*om^21-J(j)*om^2/k(j)]*T;endT21(i)=T(2,1);ifT(2,1)==0omr(mn)=om;mn=mn+1;endifT(2,1)*T21p0omr(mn)=om-step/2;mn=mn+1;endifmnnbreakendT21p=T(2,1);endomr3、最后计算出结果，并画出主阵型fori=1:nom=omr(i);SV(:,1)=[10]';forj=1:nSV(:,j+1)=[11/k(j);-J(j)*om^21-J(j)*om^2/k(j)]*SV(:,j)endMM(:,i)=SV(1,:);endplot(MM)3.3计算结果1、根据所编写的程序计算出八阶固有频率如下图2-5所示图2-52、作出的八阶主阵型图如下图2-6所示：9图2-6四、近似计算结果的验证1、采用特征根特征向量法进行精确计算并验证，具体实现程序如下所示alf=1;[Clam]=eig(inv(M)*K);fori=1:nomeig(i)=sqrt(abs(lam(i,i)));endforw=1:nC(:,w)=C(:,w).*(1/C(1,w));endomeigCx=[1234561213];plot(x,C)gridon2、运算结果结果以及图像如下图4.1和图4.2所示10图4.1图4.2五、两种近似方法及其计算结果的比