%-------------------------------------------------------------------------
%Initialisierungsdatei und
%Ausgabedatei zu gasverdichter_v5_s.mdl
%-------------------------------------------------------------------------

%Simulation des Eingangssignals

%Simulation (des Eingangssignals, das den unterschiedlich ankommenden
%Gasdruck, der zwischen 1.0 bar und 20.0 bar variiert,
%zum Gasverdichter aus dem Niedrig-Gasnetz der Stadt Münster
%symbolisiert)
sim('eingangssignal_s.mdl',700);

%Deklaration der Variablen/benoetigten
%Groessen

%Erdgas (Druck in bar)
inp1=11.5;
oup1=48.5;%gewuenschter Druck am Abfluss (zur Gasturbine)
oup2=48.5;%gewuenschter Druck am Abfluss (zur Gasturbine)

%Kenndaten des Hauptmotors
Un=10000;
In=92;
Pn=1400;
Mn=4493;
fn=50;
nnm=49.6;
J=0.875;
cos_phi=0.91;
Ia=552;
Ma=2920.45;
Mk=10109.25;
Ms=2156.64;

%Kenndaten des Getriebes
nng=60.5;

%Abmessungen/Konstante des Ventils/Schiebers
d=0.0365;
r=0.01825;
Kq1=oup1/inp1;
for i=1:60
    Kq2=oup2/esignal_array(i,1);
end

%Abmessungen und Kenndaten des Behälters/Tanks
Lt=3.5;
dt=1.0;
rt=dt/2;
La=0.811;
ha=0.811;
da=0.470;
ra=da/2;
hmax=0.8539;
ht1=0.8539;%100%
ht2=0.8539;%100%
Vt1=rt^2*Lt*(acos((rt-ht1)/rt)-(rt-ht1)*(sqrt(2*rt*ht1-ht1^2)/rt^2));
Vt2=rt^2*Lt*(acos((rt-ht2)/rt)-(rt-ht2)*(sqrt(2*rt*ht2-ht2^2)/rt^2));
A1=Vt1/ht1
A2=Vt2/ht2
Va=pi*ra^2*ha;
a=2*(Va/ha);

%Erdbeschleunigung
g=9.81;

%Simulation und Ausgabe des Modells

%Simulation 
sim('gasverdichter_v5_s.mdl',700);

%Ausgabe/Plot 
hold off
p21=plot(tsim,abfluss_konst,'b+');
set(p21,'LineStyle','-');
set(p21,'LineWidth',1.75);
hold on
p22=plot(tsim,abfluss_var,'g*');
set(p22,'LineStyle','-');
set(p22,'LineWidth',1.75);
grid on
grid MINOR
title('GASVERDICHTER - Druckerhoehung auf 45.5bar .. 51.5bar','FontSize',14,'FontWeight','Bold')
xlabel('Zeit in [sec] \rightarrow')
ylabel('Druck in [bar] \rightarrow')
legend('Ablussdruck (zur Gasturbine) (konst. Eingangsgroesse)','Ablussdruck (zur Gasturbine) (var. Eingangsgroesse)','location','SouthEast')