clear all;
close all;
clc;

%% Eingabeparameter
% Systemparameter
argTimeIncrease = 1e-3;     %[s] Zeitschritt
argFrequencyStart = 1;      %[Hz] Anregungsfrequenz zu Beginn
argFrequencyEnd = 20;       %[Hz] Anregungsfrequenz zum Ende
argFrequencyIncrease = 1;       %[Hz/s] Steigerung der Frequenz
argDuration = (argFrequencyEnd-argFrequencyStart)/argFrequencyIncrease; %[s] Simulationsdauer

argForceAmplitude = 40;    %[N] Anregungskraftamplitude an der Balkenspitze
argPreStress = 500;        %[MPa] Vorspannung

argCe = 0;                  % Dämpfung des Systems [N-s/m]                
argK = 1000;                   % Steifigkeit des Systems [N/m]

argl_0 = 0.012;             %[m] Länge der Stabfedern im entspannten Zustand
argW = 0.008;               %[m] Breite der Stabfeder
argT = 0.008;               %[m] Balkentiefe
argRho = 7850;              %[kg/m^3] Dichte (Stahl) des Balkens
argM = argW*argW*argT*argRho;   % [kg] Masse des eingespannten Balkenstücks
argArea = (argW*argW)/4*pi();   % [m^2] kreisförmiger Probenquerschnitt

% Anfangsbedingungen
argF0 = 0;              %[N] Anregungskraft zum Zeitpunkt t-2
argF1 = 0;              %[N] Anregungskraft zum Zeitpunkt t-1 
argU0 = 0;              %[m] Auslenkung zum Zeitpunkt t-2
argU1 = 0;              %[m] Auslenkung zum Zeitpunkt t-1 


%% Werkstoffdaten      
argE_A = 63300;             % [N/mm^2]=[MPa]
argE_M = 22700;             % [N/mm^2]=[MPa]
arg_eps_Tr_gr_r = 0.0521;   % maximum transformation strain [-]
arg_eps_Tr_gr_l = 0.0521;   % maximum transformation strain [-]
arg_sig_fh = 475;           % [N/mm^2]=[MPa]	
arg_sig_fr = 474.9;         % [N/mm^2]=[MPa]      

argE_T = 1500;              % durch Ablesen; attention when calculating eps_n_1 (maybe: divison by 0) -------------------------NEU-------------------------
argE_n_K = (argE_A + argE_M)/2;  % E_n zur Bestimmung von K -------------------------NEU-------------------------
argKprob = argE_n_K*argE_T/(argE_n_K-argE_T); % -------------------------NEU-------------------------

%argH = 7;                  % war = 7; Spannungsanstieg durch Temperaturerhöhung (Clausius-Clapeyron-Faktor)	
argH = 0;        

%argC = 700;               % war = 700; "frei waehlbarer Parameter" Achtung: Kelvin oder Grad Celsius?
argC = 0;   

%arg_delta = 10;	
arg_delta = 0;	

argT_0 = 20;	
argT_R = 20;

argA_s = 34.5; 		
argA_f = 49; 		
argM_s = 18.4; 		
argM_f = 9;    

%% Objekterstellung, Methodenaufruf & Übergabe von Eigenschaften

s = Solver(argDuration,argCe,argK,argl_0,argArea,argM,argTimeIncrease,argForceAmplitude,argPreStress,argE_A,argE_M,argE_T,argE_n_K,argKprob,arg_eps_Tr_gr_r,arg_eps_Tr_gr_l,arg_sig_fh,arg_sig_fr,argH,argC,arg_delta,argT_0,argT_R,argA_s,argA_f,argM_s,argM_f,argFrequencyStart,argFrequencyIncrease,argFrequencyEnd);