资源描述：
《传递矩阵-matlab程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、%main_critical.m%该程序使用Riccati传递距阵法计算转子系统的临界转速及振型%本函数中均采用国际单位制%第一步：设置初始条件（调用函数shaft_parameters）%初始值设置包括：轴段数N，搜索次数M%输入轴段参数：内径d,外径D，轴段长度l，支撑刚度K，单元质量mm，极转动惯量Jpp[N,M,d,D,l,K,mm,Jpp]=shaft_parameters;%第二步：计算单元的5个特征值（调用函数shaft_pra_cal）%单元的5个特征值：%m_k：:质量%Jp_k：极转动惯量%Jd_k：直径转动惯量%EI：弹
2、性模量与截面对中性轴的惯性矩的乘积%rr：剪切影响系数[m_k,Jp_k,EI,rr]=shaft_pra_cal(N,D,d,l,Jpp,mm);%第三步：计算剩余量（调用函数surplus_calculate）,并绘制剩余量图%剩余量：D1fori=1:1:Mptx(i)=0;pty(i)=0;endforii=1:1:Mwi=ii/1*2+50;[D1,SS,Sn]=surplus_calculate(N,wi,K,m_k,Jp_k,JD_k,l,EI,rr);D1;pty(ii)=D1;ptx(ii)=w1endylabel(‘剩余量
3、’)；plot(ptx,pty)xlabel(‘角速度red/s’);gridon%第四步：用二分法求固有频率及振型图%固有频率：Critical_speedwi=50;fori=1:1:4order=i[D1,SS,Sn]=surplus_calculate(N,wi,k,m_k,Jp_k,Jd_k,l,EI,rr);Step=1;D2=D1;kkk=1;whilekkk<5000ifD2*D1>0wi=wi+step;D2=D1;[D1,SS,Sn]=surplus_calculate(N,wi,K,m_k,Jp_k,Jd_k,l,EI,
4、rr);endifD1*D2<0wi=wi-step;step=step/2;wi=wi+step;[D1,SS,Sn]=surplus_calculate(N,wi,K,m_k,Jp_k,Jd_k,l,EI,rr);EndD1;Wi;Ifatep<1/2000Kkk=5000;endendCritical_speed=wi/2/pi*60figure;plot_mode(N,l,SS,Sn)wi=wi+2;end%surplus_calculate,.m%计算剩余量%（1）计算传递矩阵%（2）计算剩余量function[D1,SS,Sn]=
5、surplus_calculate(N,wi,K,m_k,Jp_k,Jd_k,l,EI,rr);%(1)计算传递矩阵%＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝%(a)初值设为0%＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝fori=1:1:N+1forj=1:1:2fork=1:1:2ud11(j,k.i)=0;ud12(j,k.i)=0;ud21(j,k.i)=0;ud22(j,k.i)=0;endendendfori=1:1:Nforj=1:1:2fork=1:1:2us11(j,k.i)=0;us12(j,k.i)=0;us21(j,k.i)=0;us22(j
6、,k.i)=0;endendendfori=1:1:Nforj=1:1:2fork=1:1:2u11(j,k.i)=0;u12(j,k.i)=0;u21(j,k.i)=0;u22(j,k.i)=0;endendend%＝＝＝＝＝＝＝＝＝＝＝＝%(b)计算质点上传递矩阵―――点矩阵的一部分！%＝＝＝＝＝＝＝＝＝＝＝＝fori=1:1:N+1ud11(1,1,i)=1;ud11(1,2,i)=0;ud11(2,1,i)=0;ud11(2,2,i)=1;ud21(1,1,i)=0;ud21(1,2,i)=0;ud21(2,1,i)=0;ud21(2
7、,2,i)=0;ud22(1,1,i)=1;ud22(1,2,i)=0;ud22(2,1,i)=0;ud22(2,2,i)=1;end%＝＝＝＝＝＝＝＝＝＝＝＝%(c)计算质点上传递矩阵―――点矩阵的一部分！%＝＝＝＝＝＝＝＝＝＝＝＝fori=1:1:N+1ud12(1,1,i)=0;ud12(1,2,i)=(Jp_k(i)-Jd_k(i))*wi^2;%%%考虑陀螺力矩ud12(2,1,i)=m_k(i)*wi^2-k(i);ud12(2,2,i)=0;end%＝＝＝＝＝＝＝＝＝＝＝＝%(d)以下计算的是无质量梁上的传递矩阵―――场矩阵%计
8、算的锥轴的us是不对的，是随便令的，在后面计算剩余量时，zhui中会把错误的覆盖掉%＝＝＝＝＝＝＝＝＝＝＝＝fori=1:1:Nus11(1,1,i)=1;us11
