% Systemmodellierung 
% HS 2009
% Milestone 3
% 
% Anastasiades Bryan
% Schneider Melanie
% Weber Marco
% -------------------------------
% -------------------------------
% 
% Parameter
% 
% -------------------------------
% -------------------------------

close all
clear all
clc

%Physik
rho=1000;                % Dichte Wasser [kg/m^3]
g=9.81;                  % Erdbeschleunigung
p_0=1e5;                 % Vergleichsdruck / Umgebungsdruck


%________________________________________________
%Messwerte: **noch genauer zu bestimmen**
%________________________________________________

%Hydraulisches System

lambda_T=1e-1;              % Reibungskoeffizient des Tunnels
lambda_W=1e1;               % Verlustkoeffizient des Wassserschlosses
lambda_F=1e-1;              % Reibungskoeffizient des Fallrohrs
sigma_0=1e-15;              % Kompressibilitäts des Systems

%Ventil

c_d=1;                      % Verlustkoeffizient des Ventils

%Pelton Turbine

kappa=1;                    % Trägheitsmoment der Turbine

%Generator

Theta=1;                    % Generatorkonstante



%________________________________________________
%gegebene Werte:
%________________________________________________

%Hydraulische Systeme

h_R=413;        %Wasserstand Resevoir
l_T=2000;       %Länge des Tunnels
d_T=3;          %Durchmesser des Tunnels
A_W=5^2*pi;     %Querschnittsfläche des Wasserschlosses
h_W0=390;       %Höhe des Wasserschlossbodens
l_F=741;        %Länge des Fallrohrs
d_F=2;          %Durchmesser des Fallrohrs

%Ventil

tau=1;          % Zeit für einen Ventilhub
A_V0=0.518;     % Fläche des voll geöffneten Ventils


%Pelton Turbine

d_Turb=5.35;    % Durchmesser Turbine


%Generator

R_G=10;         % Widerstand des Generators
L_G=0.1;        % Impedanz des Generators

%Lokales Stromnetz

R_N=165;        % Widerstand des Stromnetzes
L_N=0.1;        % Impedanz des Stromnetzes

R_tot=R_G+R_N;
L_tot=L_G+L_N;

%________________________________________________
%Initiall Conditions
%________________________________________________

x_IC=0;         % Startposition Ventil
h_W_IC=413;     % Höhe des Wasserschlosses
v_T_IC=0;       % Start Tunnelgeschwindigkeit
v_F_IC=0;       %Geschwindigkeit im Fallrohr
V_F_IC=l_F*pi*d_F^2/4;        %Startvolumen des Fallrohrs
omega_IC=0;         %Drehgeschwindigkeit der Turbine (Beginn)
I_IC=0;             %Generierter Strom