Нужна помощь в написании работы?

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; % Макс сила сопротивления, кН

Внимание!
Если вам нужна помощь в написании работы, то рекомендуем обратиться к профессионалам. Более 70 000 авторов готовы помочь вам прямо сейчас. Бесплатные корректировки и доработки. Узнайте стоимость своей работы.

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)); %температура рабочей жидкости

Получить выполненную работу или консультацию специалиста по вашему учебному проекту
Узнать стоимость
Поделись с друзьями