Mein MATLAB Forum

function [dy]= VentilSimulation1(t,y)
y=[0;0;0]; % Vorinitialisierung

% Kenngrössen des Öls
alpha=0.7; % Durchflusszahl
roh=0.89*10^3; % Dichte [kg/m^3]
CH=5.11*10^(-15); % Hydraulische Kapazität des Öls [m^5/N]

% Feder
CF=82731.7; % Federkonstante [N/m]

% Ventilschieberdaten
dk1=0.005; % Ventilschieberaussendurchmesser [m]
% dk2=0.0034; % Ventilschieberinnendurchmesser [m]
% lk=0.004; % Länge des Kegels [m]
mk=0.003; % Ventilschiebermasse [kg]
AK=pi/4*dk1^2; % Ventilschieberfläche [m^2]
% eps=atan(((dk1-dk2)/2)/lk); % Kegelwinkel [rad]
% AV=pi*y(1)*sin(eps)*(dk1-sin(eps)*cos(eps)*y(1)); % Eingangsquerschnitt des Stromregelventil

% Reibung
FR=0.0015; % Coulombsche Gleitreibung [N]
FD=0.003; % Dämpfungskraft [N]

% DrosselQuerschnitt
p1=250*10^5; %Eingangsdruck [N/m^2]
p3=228*10^5; %Lastdruck [N/m^2]
% AD=AV*sqrt(p1-y(3))/sqrt(y(3)-p3);

dy=[0;0;0]; %Vorinitialisierung

dy(1)=y(2);
dy(2)=-CF*y(1)/CH+AK*y(3)/mk+CF*(y(1)-0.0047)-AK*p3/mk-FR-FD;
dy(3)=0.25*y(2)/CH-alpha*0.25*sqrt(2/roh)/CH*sqrt(p1-y(3))+alpha*0.25*sqrt(2/roh)*sqrt(p1-y(3))/CH;

end

function DGLVentil %DGL numerisch lösen
close all % schließt alle offenen Plots


% Anfangsbed:
y0=[0;0;10*10^5];

% Simulationszeit:
t=0:0.1:10;

[t,Loesung]=ode45(@VentilSimulation,t,y0); % aufruf der DGL über den solver ode45 --> siehe Hilfe

figure(1)
% Zeit und Verschiebung
plot(t,Loesung(:,1),'r')

figure(2)
% Zeit und Geschwindigkeit
plot(t,Loesung(:,2),'g--')


figure(3)
% Zeit und Druck
plot(t,Loesung(:,3),'mo')

end