reynoldswending.m

来自「解析法建立挤压油膜阻尼器稳态数学模型并仿真特性关系」· M 代码 · 共 46 行

%
%
style_plot=menu('输入方式','手动','已经储存的数据','程序内固定的数据')
if style_plot==1
    %输入U L 等系数
    U=input('U');mu=input('mu');R=input('R');L=input('L');C=input('C');m=input('m');omega=input('omega');h=input('h');
    %输入delta epslion
    delta=input('delta');epsilon=input('epsilon');
elseif style_plot==1
    load jietie_ini_data
else
%     U=
%     mu=
%     R=
%     L=
%     C=
%     m=
%     omega=
%     h=
%     delta=
%     epsilon=  
end
%计算B
B=mu*R*L^3/(m*omega*C^3);

I_1=2*epsilon/(1-epsilon^2)^2;I_2=pi/(2*(1-epsilon^2))^(3/2);
Delta=epsilon-delta^2*U/sqrt((B*I_1*delta+1-delta)^2+(B*I_2*delta)^2);
%判断Delta
if Delta~=0 
    Delta_1=Delta;epsilon_1=epsilon;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
    epsilon=epsilon+h;
    I_1=2*epsilon/(1-epsilon^2)^2;I_2=pi/(2*(1-epsilon^2))^(3/2);
    Delta=epsilon-delta^2*U/sqrt((B*I_1*delta+1-delta)^2+(B*I_2*delta)^2);
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    while Delta>=0
        epsilon=epsilon+h;
        I_1=2*epsilon/(1-epsilon^2)^2;I_2=pi/(2*(1-epsilon^2))^(3/2);
        Delta=epsilon-delta^2*U/sqrt((B*I_1*delta+1-delta)^2+(B*I_2*delta)^2);        
    end
    
    Delta_2=Delta;epsilon_2=epsilon;epsilon=(epsilon_1+epsilon_2)/2;
end

delta,epsilon