齿轮动力学.rar_MATLAB齿轮仿真_动力学仿真_齿轮 matlab_齿轮动力学_齿轮动力学‘

function dx=RK_fun(t,x,W) b(1)=0.2;b(2)=0.2;b(3)=0.2;b(4)=0.2;b(5)=1; %量纲一化后轴承径向游隙及齿侧间隙 rp=0.033829;rg=0.050744; %基圆半径 mp=0.676;mg=1.084; %齿轮质量 Jp=0.0407;Jg=0.1168; %齿轮惯量 Cxp=4860;Cyp=4860;Cxg=6160;Cyg=6160;Ch=6000; %支承阻尼及啮合阻尼 Kxp=3.5e8;Kyp=3.5e8;Kxg=3.5e8;Kyg=3.5e8;Km=3.04e8;Kr=0.4e8;%支承刚度及啮合刚度、波动值 wp=15;Tp=540;Tg=810; %齿轮转速、转矩 Kh=Km+Kr*cos(16*wp*t); %时变啮合刚度 aa=0.324;ac=0.413;ad=0.717;bb=0.392; %齿轮展角及啮合角 lp=rp*(aa+mod(wp*t,(ad-aa))); %小齿轮摩擦力臂 lg=(rp+rg)*tan(bb)-lp; %大齿轮摩擦力臂 nn=sign((ac-aa)-mod(wp*t,ad-aa));u=0.1; %摩擦力方向系数 me=Jp*Jg/(rp*rp*Jg+rg*rg*Jp); %齿轮等效质量 Wn=sqrt(Km/me); %齿轮固有频率 cxp=Cxp/(2*mp*Wn);cyp=Cyp/(2*mp*Wn);cxg=Cxg/(2*mg*Wn);cyg=Cyg/(2*mg*Wn); chp=Ch/(2*mp*Wn);chg=Ch/(2*mg*Wn);ch=Ch/(2*me*Wn); %量纲一化后的阻尼系数 g=1+me*u*nn*rp*lp/Jp+me*u*nn*rg*lg/Jg; %传动误差的部分系数 kxp=Kxp/(mp*Wn*Wn);kyp=Kyp/(mp*Wn*Wn);kxg=Kxg/(mg*Wn*Wn);kyg=Kyg/(mg*Wn*Wn); khp=Kh/(mp*Wn*Wn);khg=Kh/(mg*Wn*Wn);kh=1+Kr*cos(W*t)/(me*Wn*Wn); %量纲一化后的刚度系数 fp=rp*Tp/(Jp*0.00005*Wn*Wn);fg=rg*Tg/(Jg*0.00005*Wn*Wn);fe=0.2*W*W*cos(W*t); %量纲一化后的转矩及误差激励 M=[1,0,0,0,0;0,1,0,0,0;0,0,1,0,0;0,0,0,1,0;0,-1,0,1,1]; C=[(2*cxp),0,0,0,(-2*u*nn*chp);0,(2*cyp),0,0,(-2*chp);0,0,(2*cxg),0,(2*u*nn*chg);0,0,0,(2*cyg),(2*chg);0,0,0,0,(2*g*ch)]; K=[(kxp),0,0,0,(-u*nn*khp);0,(kyp),0,0,(-khp);0,0,(kxg),0,(u*nn*khg);0,0,0,(kyg),(khg);0,0,0,0,(g*kh)]; P=[0;0;0;0;(fp+fg-fe)]; F=zeros(5,1); if x(1)>b(1) F(1)=x(1)-b(1); elseif x(1)<(-b(1)) F(1)=x(1)+b(1); else F(1)=0; end if x(2)>b(2) F(2)=x(2)-b(2); elseif x(2)<(-b(2)) F(2)=x(2)+b(2); else F(2)=0; end if x(3)>b(3) F(3)=x(3)-b(3); elseif x(3)<(-b(3)) F(3)=x(3)+b(3); else F(3)=0; end if x(4)>b(4) F(4)=x(4)-b(4); elseif x(4)<(-b(4)) F(4)=x(4)+b(4); else F(4)=0; end if x(5)>b(5) F(5)=x(5)-b(5); elseif x(5)<(-b(5)) F(5)=x(5)+b(5); else F(5)=0; end dx=[x(6:10);inv(M)*(P-C*x(6:10)-K*F)];