%% Standard
clear all       %alles löschen
close all       %alles schließen
clc             %Commander Fenster säubern

%% Anzeigen der Aufgabe

disp('Aufgabe sowieso')

%% Festlegen von Konstanten

Hbem = input('Hbem :');     %Aufforderung der Eingabe nach Hbem
rw = 1.025;                 %t/m³ %Meerwasser
g = 9.81;

    %cot(a)|Anz. Lagen| P(%|  kd| We.b.flanke|We.b.kopf| rs
                                    %b |n.b     %b |n.b
%Spalte 1      2         3    4     5    6       7  8     9
                                       
N = [1.5       2        38   1.02  1.2  2.4     1.1 1.9  2.6]; %Naturstein Glatt
B = [1.5       2        37   1.15   2    4      1.9 3.2  2.6]; %Bruchsteine scharkantig
T = [1.5       2        50   1.04   7    8      5   6    2.6]; %Tetrapoden
%% Berechnung Naturstein Glatt
rs = N(9);
kdn = N(4);
c = N(1);
n = N(2);
rs = B(9);
    
Wn = (rs*g*Hbem.^3)./(kdn.*(rs/rw-1)^3*c);    %erforderliches Steingewicht
    
fprintf('W für Natursteinglatt (kd=1.02): %.2f kN\n', Wn)

plot(kdn,Wn,'+-','LineWidth',2.5,'MarkerEdgeColor','r','MarkerFaceColor','r','MarkerSize',10)
hold on
grid on

%% Berechnung Bruchstein scharfkantig
rs = B(9);
kdb = B(4);
c = B(1);
n = B(2);
rs = B(9);
    
Wb = (rs*g*Hbem.^3)./(kdb.*(rs/rw-1)^3*c);  %erforderliches Steingewicht

fprintf('W für Bruchstein scharfkantig (kd=1.15): %.2f kN\n', Wb)
    
plot(kdb,Wb,'+-','LineWidth',2.5,'MarkerEdgeColor','g','MarkerFaceColor','g','MarkerSize',10)
hold on
grid on
%% Berechnung Tetrapoden
rs = T(9);
kdT = T(4);
c = T(1);
n = T(2);
rs = T(9);
    
WT = (rs*g*Hbem.^3)./(kdT.*(rs/rw-1)^3*c);    %erforderliches Steingewicht

fprintf('W für Tetrapoden (kd=1.04): %.2f kN\n', WT)
    
plot(kdT,WT,'+-','LineWidth',2.5,'MarkerEdgeColor','m','MarkerFaceColor','m','MarkerSize',10)
hold on
grid on

%% Variieren von kd
kd1 = [0.8:0.1:2];
rs = 2.6;
c = 1.5;
g = 9.81;

W = (rs*g*Hbem^3)./(kd1.*(rs/rw-1)^3*c);

disp(' ')

fprintf('W (kd=0.8): %.2f kN\n', W(1))
fprintf('W (kd=0.9): %.2f kN\n', W(2))
fprintf('W (kd=1.0): %.2f kN\n', W(3))
fprintf('W (kd=1.1): %.2f kN\n', W(4))
fprintf('W (kd=1.2): %.2f kN\n', W(5))
fprintf('W (kd=1.3): %.2f kN\n', W(6))
fprintf('W (kd=1.4): %.2f kN\n', W(7))
fprintf('W (kd=1.5): %.2f kN\n', W(8))
fprintf('W (kd=1.6): %.2f kN\n', W(9))
fprintf('W (kd=1.7): %.2f kN\n', W(10))
fprintf('W (kd=1.8): %.2f kN\n', W(11))
fprintf('W (kd=1.9): %.2f kN\n', W(12))
fprintf('W (kd=2.0): %.2f kN\n', W(13))


plot(kd1,W,'bo-','LineWidth',2.5,'MarkerEdgeColor','k','MarkerFaceColor','k','MarkerSize',10)
hold on
grid on

%% Bezeichnung Ausdruck
title('Variieren von kd') %Titel anzeigen
xlabel('kd')
ylabel('W')                                          %Achsenbeschriftung
legend('Naturstein Glatt','Bruchstein scharfkantig','Tetrapoden','Variieren von kd',1) %Legende