clc; clear all;



x=0:0.0005:0.264;

Volumen_Hub=1.1e-4       
t_Hub=0.264
d_Schlauch=0.223;
Q_Pumpe_pro_Hub=Volumen_Hub/t_Hub                          
c_Pumpe_pro_Hub=4*Q_Pumpe_pro_Hub/(d_Schlauch^2*pi)        

W_kin=0.5*1380*Volumen_Hub*c_Pumpe_pro_Hub^2

Beschleunigung_Hub=W_kin/(0.5*1380*Volumen_Hub*(c_Pumpe_pro_Hub/t_Hub)*t_Hub^2)

V=(d_Schlauch^2*pi/4)*Beschleunigung_Hub*x.^2

t_Pumpe_Impuls=[-0.05 0 x 0.264 0.3];
Volumen_Impuls=[0 0 V 0 0];

figure(1);
subplot(2,1,1), plot(t_Pumpe_Impuls, Volumen_Impuls);grid;
title('Pumpenhub');
xlabel('Zeit [Sek]');  
ylabel('Volumen [m³]');
axis([-0.02, 0.3, -0.2e-4, 1.5e-4])
set(gca,'FontName','calibri','FontSize',20);
ti=get(gca,'Title'); set(ti,'FontName','calibri','FontSize',22);
xl=get(gca,'xlabel'); set(xl,'FontName','calibri','FontSize',22);
yl=get(gca,'ylabel'); set(yl,'FontName','calibri','FontSize',22);


F=Beschleunigung_Hub*1380*d_Schlauch^2*pi/4*Beschleunigung_Hub*x.^2;

Kraft_Impuls=[0 0 F 0 0];
figure(1);
subplot(2,1,2), plot(t_Pumpe_Impuls, Kraft_Impuls);grid;
xlabel('Zeit [Sek]');  
ylabel('Kraft [N]');
axis([-0.02, 0.3, -0.3, 0.8])
set(gca,'FontName','calibri','FontSize',20);
ti=get(gca,'Title'); set(ti,'FontName','calibri','FontSize',22);
xl=get(gca,'xlabel'); set(xl,'FontName','calibri','FontSize',22);
yl=get(gca,'ylabel'); set(yl,'FontName','calibri','FontSize',22);