故障轉(zhuǎn)子系統(tǒng)的非線性振動(dòng)分析與診斷方法-免費(fèi)閱讀
【正文】 tt=t((l+1):(lt1))。[I X(i)]]。nIMF%%該函數(shù)用來計(jì)算上下包絡(luò)的最大值和最小值function[Xa,Xb]=ProfileData(X)Ns=length(X)。h1=X。if SD h1=h1k。 [Ma,Na]=size(Xa)。 %原始信號(hào)左右延拓t=0(1/fs):((3*Ns1)/fs)。 xlabel(‘Frenquency/Hz’, ‘Fontsize’, 8)。%Hilbert時(shí)頻譜c=c’ 。for i=1:nIMF Hc1=imag(Hilbert(c(:,i)))。nIMFsubplot(nIMF+1,1,i)plot(t,c(:,i))。x=X(:,2)。ps=sum(p)。s(11,:)=wprocef(t,[4 12])。s(13,:)=wprocef(t,[4 8])。s(7,:)=wprocef(t,[4 4])。%基于小波包分解的能量比例計(jì)算s(1,:)=wprocef(t,[4 0])。dt1=1/fs1。T=dt*[0:N1]。 F (local*44+4,1)=T*sin(alfa)*cos(bita)。 Ati=fun(t,yi,w)。d(6)=(c11/m1)*y(6)(c12/m1)*y(8)(k11/m1)*y(5)(k12/m1)*y(7)+e1*omega*omega*cos(omega*t)+fy1/m1。fx2=kr2*(y(7)delta2)*Dirac2。c21=。 %碰摩摩擦系數(shù)1fr2=。e1=*0。endfor i=1:1:N yout(i)=b(i)+h*d(i)/6。 a(4)=h。end%function yout=ekutta21(t,h,y,A)N=length(y)。d=fun(t,y,w)。for i=1:1:N a(i)=0。Drds=eye(N)*1S。 for i=1:1:Fen*1 t=0。% 第二步：通過打靶法計(jì)算振幅的分岔特性shooting。dxn1=dxn+a6*ddxn+a7*ddxn1。endKj=KT+a0*M+a1*C。 K1(loc_rub*44+2,loc_rub*44+1)=k_rub*(rrdert)*miu_rub。k_rub=5e4。K1=K*0。 D(i*4+24,i*4+24)=D(i*4+24,i*4+24)+Dzy(i)。beita=2*(yita2*wi2yita1*wi1)/(1/wi1^2wi1^2)。endformat long g[Ax WW]=eig(inv(M)*K)。for i=1:1:N+1 K(i*4+14,i*4+14)= K(i*4+14,i*4+14)+Kzx(i)。 end endendfor i=2:1:Nfor j=1:1。 Ks(j,k,i)=Ks(k,j,i)。 Ks(8,8,i)=4*L(i)^2。Ks(6,5,i)=0。 Ks(7,2,i)=6*L(i)。 Ks(5,2,i)=0。Ks(3,1,i)=0。 Ge(7,7,i)=0。 Ge(8,2,i)=3*L(i)。Ge(6,2,i)=0。 Ge(4,2,i)=3*L(i)。 Msr(8,8,i)=4*L(i)^2。Msr(6,5,i)=0。 Msr(7,2,i)=3*L(i)。 Msr(5,2,i)=0。 Msr(2,2,i)=36。 Mst(6,6,i)=156。 Mst(7,3,i)=3*L(i)^2。 Mst(5,2,i)=0Mst(5,3,i)=0。Mst(3,1,i)=0。L=[ 30 20 ]/1000。loca=6。 xt1=[xt1。endt=0。a5=ht/*(gama/)。ht=2*pi/wi/Fen。for i=1:1:N+1 D1=D1*sign(det(Su(:,:,i))。 ud12= u12(:,:,i)。 u21(:,:,i)=us21(:,:,i)*ud11(:,:,i)+us22(:,:,i)*ud21(:,:,i)。cone_modify（16，wi）。 Us12(i,j)=0。 SS (2,1,i)=0。 Su (2,1,i)=0。 us21(2,1,i)=1(i)/EI(i)。 us11 (2,1,i)=0。 ud22 (2,1,i)=0。 ud11 (2,1,i)=0。 us21(j,)=0。plot_mode(N,l,SS,Sn) wi=wi+2。 end if D1*D20 wi=wistep。grid on% 第四步：用二分法求固有頻率及振型圖%固有頻率：Critical_speedwi=50。% 第二步：計(jì)算單元的5個(gè)特征值（調(diào)用函數(shù)shaft_pra_cal）%單元的5個(gè)特征值：%m_k：:質(zhì)量%Jp_k：極轉(zhuǎn)動(dòng)慣量%Jd_k：直徑轉(zhuǎn)動(dòng)慣量%EI：彈性模量與截面對(duì)中性軸的慣性矩的乘積%rr：剪切影響系數(shù)[m_k,Jp_k,EI,rr]=shaft_pra_cal(N,D,d,l,Jpp,mm)。 D1。while kkk5000 if D2*D10 wi=wi+step。Wi。 ud22(j,)=0。 u22(j,)=0。 ud21(2,2,i)=0。ud12(2,2,i)=0。 us12(2,2,i)=0。 us22(2,2,i)=1。 Sn(2,2,i)=0。 Ud21(i,j)=0。cone_modify（6，wi）。cone_modify（24，wi）。 u22(:,:,N+1)=ud22(:,:,N+1)。Sn(:,:,i)= inv(ud21* SS1+ ud22)。yt1=[]。a2=(beita*ht)。N3 yn(i:1)=0。 sitax(I,1)=yn(i+40,n)/pi*180 end u=F(loca*44+1,1)。save ’’ yt1 –ascii。fai=[30 30]/180*pi %%該程序?qū)effcott轉(zhuǎn)子系統(tǒng)進(jìn)行參數(shù)設(shè)置 function [N density Ef L R R0 miu]=rotor_parameters%N: 劃分的軸段數(shù)%density: 軸的密度 單位為：[kg/m^3]%Ef: 軸的彈性模量 單位為:[Pa]%L 每個(gè)軸段的長度 單位為:[m]%R 每個(gè)軸段的外半徑 單位為:[m]%R0: 每個(gè)軸段的內(nèi)半徑 單位為:[m]%miu: 每個(gè)軸段的單元質(zhì)量 [kg/m]N=11。NN1=NN0+1NN2=NN1+1for i=1:1:NN0 Mst(1,1,i)=156。 Mst(4,3,i)=0。 Mst(6,4,i)=0。 Mst(8,4,i)=3*L(i)^2。 Mst(8,6,i)=0 Mst(8,7,i)=0Mst(8,8,i)=4*L(i)^2。 Msr(4,3,i)=0。 Msr(6,3,i)= 3*L(i)。Msr(8,3,i)=0。 Msr(8,5,i)=3*L(i)。 Ge(3,2,i)=0。Ge(5,3,i)=3*L(i)。 Ge(7,3,i)=0。 Ge(6,6,i)=0。endfor i=1:1:NN0Ks(1,1,i)=12。 Ks(4,3,i)=0。 Ks(6,3,i)=6*L(i)。Ks(8,3,i)=0。 Ks(8,5,i)=6*L(i)。 Ks(j,k,i)=Ks(j,k,i)*EI/L(i)^3。 end endengfor i=1:1:N for j=1:1:8 for K=1:1:8 M(i*4+j4,i*4+k4)=Ms(j,k,i)。Kzy(1)=7e7。Kzy(12)=7e7。yita1=。Dzx(12)=7e3。if i==1xn=yn(:,n)。Pn1=F。if nFen*60amp。endKT=K0+K1。end[,RR]=qr(Kj)。 穩(wěn)定性分析程序% 第一步：設(shè)置初始條件[N Fen y d y2]=initial_conditionsfm=[]。 mod(i,Fen)==1 fprintf(f1,w/2/pi,y(1),y(2),y(3),y(4), y(5),y(6),y(7),y(8))。endends=y。w/2/pi floquet_mul]。 a(2)=h/2。 end tt=t+a(k)。enda(1)=h/2。 b(i)=b(i)+a(k+1)*d(i)/3。l1=。 %碰摩間隙2kr1=1e5。k21=k12。endif y(5)deltalDirac1=0。d(3)=y(4)。for i=1:
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1