机械振动--盘轴扭振系统固有频率和主振型的计算.docxVIP

机械振动--盘轴扭振系统固有频率和主振型的计算 机械振动大作业 （盘轴扭振系统固有频率和主振型的计算） 学院：航空航天工程学部 班级 姓名：李根 学号：2021040402093 2021年 5 月 12号 盘轴扭振系统 固有频率和主振型的计算 一：简化分析 该系统为非约束性盘轴扭振系统，并简化分析： 1．忽略轴的质量; 2．轴的刚度对盘的影响不做考虑; 3．将圆盘的质量集中于圆盘中心，不考虑圆盘厚度对系统的影响; 4．系统为线弹性系统，盘为刚体。 对于非约束系统，其只存在刚度矩阵，不存在柔度矩阵，即不能对刚度矩阵求逆。 二：条件 圆盘： 1. 几何尺寸：直径d1?0.4m，厚度h?0.02m； 2. 材料：杨氏模量E?2?1011(N/m2)，剪切模量G?7.69?1010(N/m2) ) 密度??7800(kg/m3轴： 1. 几何尺寸：直径d2?0.04m，a?0.1m 2. 材料：杨氏模量E?2?1011(N/m2)，剪切模量G?7.69?1010(N/m2) ) 密度??7800(kg/m3 三（1）：矩阵迭代法 1.概述 ?1?1???（1）：系统主振型方程为?2?A???，引入动力矩阵D?MMKA?????????K?。任??- 2 - 取一个经过归一化的假设阵型?A?0,用动力矩阵[D]前乘它，并对通过乘法运算新得到的阵型矢量进行归一化，则得：[D]?A?0?a1{A}1， 式中a1为新振型矢量归一化后的系数。 （2）若{A}1??A?0，从{A}1开始，重复上述步骤得：[D]?A?1?a1{A}2， 式中a2为新振型矢量归一化后的系数。 （3）：若{A}2??A?1，继续重复上述步骤，进过K次矩阵乘法运算后，得到 [D]?A?k?1?ak{A}k，在规定的有效位数内，{A}k??A?k?1时停止运算，此时的?A?k?1即 为系统第一阶主振型{A(1)}的近似值，即：{A(1)}??A?k?1，而这时的系数ak即是系统第一阶固有频率平方倒数的近似值，即：?12?1/ak。 该方法的精确度不依赖于假设阵型，假设阵型的好坏只影响迭代的次数。即使假设的固有频率域一阶主振型相差很远，经过充分的迭代运算，仍可求得足够精确的基频值。 求得第一阶主振型以后，利用主振型的正交性来清除掉假设阵型中的分量，然后再进行迭代求解可以是结果收敛于第二阶主振型。同理，如果我们在假设阵型中清除掉所有前s阶这阵型分量，那么迭代的结果将得到第s+1阶固有频率及主振型。 (j)（j）Ts引进清型矩阵：[Q]?[I]??{A}{A}?M?。由于实际计算中舍入误差的存在，每次迭 j?1Mj代后，所得的主振型中还包含前面几阶的主振型分量，因此每次计算前都要进行清型才能保证最后收敛的主振型。 2.计算程序 clc clear n=8; d1=0.4; %圆盘直径 d1=0.4 d2=0.04; %轴直径 d2=0.04 a=0.1; %轴几何尺寸 den=7800; %密度 （轴和圆盘） G=7.69e+10; %剪切模量 h=0.02; %圆盘厚度 J1=0.5*pi*den*h*(d1/2)^4; %转动惯量 - 3 - Ip=pi*d2^4/32; %质心 k=G*Ip/a; %刚度 %假设阵型中刚度 k1=k; k2=k; k3=k; k4=k; k5=k; k6=k; k7=k/5; k8=k; J2=J1; J3=J1; J4=J1; J5=J1; J6=J1; J7=J1; J8=J1; J=diag([J1,J2,J3,J4,J5,J6,J7,J8]) %刚度矩阵 K=[k2,-k2,0,0,0,0,0,0; -k2,k2+k3,-k3,0,0,0,0,0; 0,-k3,k3+k4,-k4,0,0,0,0; 0,0,-k4,k4+k5,-k5,0,0,0; 0,0,0,-k5,k5+k6,-k6,0,0; 0,0,0,0,-k6,k6+k7,-k7,0; 0,0,0,0,0,-k7,k7+k8,-k8;