混沌.zip_动力学_混沌 动力学_混沌动力学_齿轮_齿轮 动力学

function dy=Chaos_3(t,x) % m_g1=1; % m_g2=1; %齿轮重量 % r_b1=0.04; % r_b2=0.04; %齿轮半径 % I_1=0.0008; % I_1=0.5*m_g1*(r_b1^2); % I_2=0.0008; % I_2=0.5*m_g2*(r_b2^2); % % g_1=0; % g_1=F_1/[m_g1*b_3*(k_n^2)]; % g_2=0; % g_2=F_2/[m_g2*b_3*(k_n^2)]; % g_m=0.1; % g_m=F_m/[m_e*b_3*(k_n^2)]; % g_a=0.05; % g_a=e_a/b_3; % % a_13=0.0125; % a_23=0.0125; % a_33=0.05; % a_11=0.1; % a_22=0.1; % % m_e=(I_1*I_2)/[I_1*(r_b2^2)+I_2*(r_b1^2)]; % % % k_m=2*10^8; % X=0.2; % k_a=k_m*X; % % K_n=k_m/m_e; % K=1.53; % 无量纲啮合频率 % K_h=K*K_n; % %% % % t=0:0.01:100;%时间范围 % h=1; % k_h=k_m+k_a*sin(K_h*t+h);%% 需要编自编函数 % k_b1=5*k_h; % k_b2=k_b1; % % %% % k_11=(k_b1*m_e)/(m_g1*k_m); % k_22=(k_b2*m_e)/(m_g2*k_m); % k_33_q=1-X*cos(K_h*t); % k_13_q=k_33_q*m_e/m_g1; % k_23_q=k_33_q*m_e/m_g2; % % dy=zeros(6,1); %% 微分方程 % dd_x_1=-g_1-2*a_11*d_x_1-2*a_13*d_x_3-k_11*f(x_1)-k_13_q*f(x_3); % dd_x_2=g_2-2*a_22*d_x_2+2*a_23*d_x_3-k_22*f(x_2)+k_23_q*f(x_3); % dd_x_3=g_m+g_a*sin(K_h*t)*K^2+dd_x_1-dd_x_2-2*a_33*d_x_3-k_33_q*f(x_3); %% % dy(1)=y(2); % dy(2)=-g_1-2*a_11*y(2)-2*a_13*y(4)-k_11*f(y(1))-k_13_q*f(y(3)); % dy(3)=y(4); % dy(4)=g_2-2*a_22*y(4)+2*a_23*y(4)-k_22*f(y(2))+k_23_q*f(y(3)); % dy(5)=y(6); % dy(6)=g_m+g_a*sin(K_h*t)*K^2+dy(2)-dy(4)-2*a_33*y(5)-k_33_q*f(y(3)); a_11=0; a_13=0.0125; % global K; K=1.53; h=0; k_h=4*K*2*10^8;%----------------------------K无量纲啮合频率 k_11=1.25+0.25*sin(k_h*t+h); k_13=0.25-0.05*cos(k_h*t); %------------------------- a_22=0; a_23=0.0125; k_22=k_11; k_23=k_13; %---------------------------- a_33=0.05; k_33=1-0.2*cos(k_h*t); g_1=0; g_2=0; g_m=0.1; g_a=0.05; %% 另外一种微分方程方法 f1=x(2); f2=-g_1-2*a_11*x(2)-2*a_13*x(6)-k_11*f(x(1))-k_13*f(x(5)); f3=x(4); f4=g_2-2*a_22*x(4)+2*a_23*x(6)-k_22*f(x(3))+k_23*f(x(5)); f5=x(6); f6=g_m+g_a*sin(k_h*t)*K^2+f2-f4-2*a_33*x(6)-k_33*f(x(5)); dy=[f1;f2;f3;f4;f5;f6]; % function out=f(in) b_1=0.0001; b_2=0.0001; b_3=0.0001; % 轴承径向间隙 suitiao % b_3= 齿轮副侧隙 B=[b_1/b_3,b_2/b_3,b_3/b_3]; % var_nam=inputname(1); % switch var_nam % case 'y(1)' % i=1; % case 'y(2)' % i=2; % case 'y(3)' % i=3; % end i=1; if in>B(i) out=in-1; elseif in<-B(i) out=in+1; else out=0; end