clc;
Pnom=16e+6; %Номинальное давление, Па
a1=62;
b1=32;
c1=23 ;%коэф нагрузок на штоках ГЦ, кН, кН/м, кН/м2
x=0.75 ;% ход штоков, м
hmax=x;
Vmax=0.07; % максимальная скорость штоков, м/с
m=55 ;% масса грузов, кг
M=240 ;% момент на выходном валу ГМ, Нм
n=120 ;%частота врашение вала ГМ, об/мин
J=2.8 ;%момент инерции нагрузки на валах ГМ, кгм2
%Выбор ГЦ
Fmax=a1+b1*x+c1*x^2; % Макс сила сопротивления, кН
P1=0.95*Pnom; %давление в напорной полости
P2=0.05*Pnom; %давление в сливной полости
fi=0.49; %коэф асимметрии рабочих полостей цилиндра
etamc=0.95; % мех КПД
%dp=sqrt(4*Fmax*1000/(pi*(P1-P2)*(1+fi)*etamc)) %диаметр поршня
dp=0.080;
Apn=pi*dp^2/4;
%dsh=dp*sqrt(1-fi) %диаметр штока
dsh=0.063;
Aps=pi*(dp-dsh)^2/4;
%Выбор ГМ
etap=0.95 ;%КПД редуктора
Ptm=M*2*pi*n*60/etap; %полезная мощность ГМ
%Принимаем ГМ 210.16
qm=28.1e-6; %рабочий обьём гидромотора, м3
etamm=0.95 ;%мех КПД ГМ
Mm=qm*(P1-P2)*etamm/(2*pi); %момент на валу выбранного ГМ
up=M/(Mm*etap); %передаточное число редуктора
ntm=n*up; % требуемая частота вращения вала ГМ, об/мин
%Определение параметров и выбор насоса
Aef=pi*dp^2/4; %Эфективная площадь порщня
etaoc=0.98 ;%обьёмный КПД ГЦ
etaoa=0.94 ;%обьёмный КПД гидроаппаратов
Qc=Aef*Vmax/(etaoc*etaoa); %Расход ГЦ, м3/с
etaom=0.965; %обьёмный КПД ГМ
Qm=qm*ntm/60/(etaom*etaoa); %Расход ГМ, м3/с
Qmax=Qc+Qm; %макс расход, м3/с
etaon=0.95; %обьёмный КПД насоса
nn=1920/60; % частота вращения вала насоса, об/с
qnt=Qmax/(nn*etaon); %рабочий обьём насоса, м3
%Принимаем насос 210.16
etam=0.95; %мех КПД насоса
Qnt=52.1/60/1000;%теоретическая подача насоса,м3/с
qn=28.1e-6; %рабочий обьём выбранного насоса, м3
Qn=qn*nn*etaon; %подача насоса,м3/с
%Выбор параметров трубопроводов
vg1=1.2; %скорость жидкости для всасывающих трубопроводов
vg2=2 ;%скорость жидкости для сливных трубопроводов
vg3=5 ;%скорость жидкости для напорных трубопроводов
% d1=sqrt(4*Qmax/(pi*vg1)) %диаметр трубопровода,м
% d2=sqrt(4*Qc/(pi*vg3))
% d3=d2
% d4=sqrt(4*Qc/(pi*vg2))
% d5=sqrt(4*Qm/(pi*vg3))
% d6=sqrt(4*Qm/(pi*vg2))
% d7=d4
% d8=d5
% d9=d6
d1=0.025 ;
d2=0.01 ;
d3=d2 ;
d4=0.016;
d5=0.008;
d6=0.012 ;
d7=d4 ;
d8=d5 ;
d9=d6;
%Определение обьёма гидробака
Vnom=Qn*60/3; %номинальная ёмкость бака
%Определение параметров гидрораспределителя
Frmax=5000; %Н
alfa=0.1;
mu1=100;
xi=0.6;
%dz1=sqrt(4*Qc/(pi*(1-xi^2)*vg3)) %диаметр золотника распределителя, м
%dz2=sqrt(4*Qm/(pi*(1-xi^2)*vg3))
dz1=0.013;dz2=0.011;
% dsh1=xi*dz1 %диаметр шейки распределителя, м
% dsh2=xi*dz2
dsh1=0.008; dsh2=0.007;
vd=2*vg3; %скорость потока жидкости в дроссельном окне, м/с
Amax1=Qc/vd; %макс площадб проходного сечения
Amax2=Qm/vd;
y1=0.003; y2=0.002; %величина перекрытия, мм
xp1=Amax1/(pi*dz1)+y1; %рабочее смещение золотника
xp2=Amax2/(pi*dz2)+y2;
%Определение параметров переливного клапана
pkl=Pnom; %давление настройки клапана
Vkmax=20 ;% макс скорость течения через клапан, м/с
%dkl=sqrt(4*Qn/(pi*Vkmax)) %диаметр золотника клапана, м
dkl=0.007;
Fpr0=pkl*pi*dkl^2/4 ;%усилие пружины клапана,Н
l0=1e-3; %предворительная деформация,м
l=3e-3; % перемещение соот началу открытия,м
x0=l0+l ;%деформация пружины в момент открытияпроходного сечения клапана,м
cpr=Fpr0/x0; %коэф жёсткости пружины,Н/м
Rez=;
end
figure(1)
plot(Rez(:,1),Rez(:,2))
xlabel('Q')
ylabel('P')
grid on
%Вычисление коэффициентов уравнений модели
nu = 0.15e-6; %кинематическая вязкость
E=1.5e+9;
l2=1;
l31=0.75;l32=0.75;
l71=0.75;l72=0.75;
A2=pi*d2^2/4;
A3=pi*d3^2/4;
A7=pi*d7^2/4;
mugl_1=25.2*po*nu*(l2/A2^2); %коэф учит линейные потери
mugl_2=25.2*po*nu*((l31+l32)/A3^2);
mugl_3=25.2*po*nu*((l71+l72)/A7^2);
lamdat=0.025;
zetta2=1;zetta3=0.75+0.75;zetta7=0.75+0.75;
mugn_1=0.443*lamdat*po*(l2/(A2^2.5))+po/2*(zetta2/A2^2);%коэф учит нелинейные потери
mugn_2=0.443*lamdat*po*((l31+l32)/(A3^2.5))+po/2*(zetta3/A3^2);
mugn_3=0.443*lamdat*po*((l72+l72)/(A7^2.5))+po/2*(zetta7/A7^2);
V2=l2* A2; V3=(l31+l32)*A3; V7=(l71+l72)*A7;
mg_1=po*(V2/A2^2); %коэф масс
mg_2=po*(V3/A3^2);
mg_3=po*(V7/A7^2);
m1=2;m2=10;mu2=100;
cst=1e+6;must=2500;
%Анализ статических состояний гидропривода
Pnymax=pkl+(kkl*Qnt*(Pnom-(1-etaon)*pkl))/(Pnom+kkl*(1-etaon)*Qnt); %макс давление созд насосной установкой
Qnk=(1-pkl/Pnom*(1-etaon))*Qnt; %подача насоса в момент срабатывания клапана
Rez1=;Rez3=;
p=F/Q;
P_sum=del_Pt1+del_Pt2+del_Pm1+del_Pm2+p;
Q_graf=;
P_garf=;
Rez1=;
Rez2=;
end
v=Q_graf*pi*dp^2/4*1e+4; %скороть движения штока ГЦ
F1=P_garf*(pi*dp^2/4);
Pvih=Fmax*v; %мощность на штоке ГЦ
eta_gp=Pvih/Pn; %КПД гидропривода
Rez3=;
figure(2)
hold on
plot(Rez1(:,1),Rez1(:,2),Rez2(:,1),Rez2(:,2),Rez2(:,1),Rez2(:,3),Rez2(:,1),Rez2(:,4),Rez2(:,1),Rez2(:,5))
xlabel('Q')
ylabel('Pny')
grid on
figure(3)
plot(F1, v)
xlabel('F')
ylabel('v')
grid on
figure(4)
plot(v, eta_gp)
xlabel('v')
ylabel('кпд')
grid on
%Анализ тепловой нагруженности гидропривода
F=(1-eta_gp)*Pn; %тепловой поток
alfa1=15;
Ths=35+max(F)/(alfa1*0.065*0.8*Vnom^(2/3)); %температура рабочей жидкости
Поможем написать любую работу на аналогичную тему