Mein MATLAB Forum

function DGL3 %DGL numerisch lösen

close all % schließt alle offenen Plots
clear all % löscht den Workspace
clc % löscht alle Einträge im command-window

%Variablen:
global x p2 y
global lk dk1 dk2 eps alpha roh CH CF mk AK AV FR FD p1 p3 AD y0

%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*x*sin(eps)*(dk1-x*sin(eps)*cos(eps)); %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-p2)/sqrt(p2-p3);



%Anfangsbed:
y0=[0;0;10];

%Simulationszeit:
t=0:0.1:10;
x=0:0.0001:0.0047;

[t,y]=ode23('VS1',t,y0); % aufruf der DGL über den solver ode45 --> siehe Hilfe

figure(1)
% Zeit und Verschiebung
plot(t,y(1))

figure(2)
% Zeit und Geschwindigkeit
plot(t,y(2))

figure(3)
% Zeit und Druck
plot(t,y(3))


end

function dy=VS1(y)
y(1)=x; %Zustandsvektor
y(2)=Dx;
y(3)=p2;

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