clear all
clc



R_G=15.74;      %Nocken-Grundkreisradius, mm
r=42.5;         %Krümmungsradius der Druckfläche, welche mit dem Nocken zusammenarbeitet, mm
i=1.2;          %Schlepphebelverhältnis
Z=5/8;

%Wählbare Größen
H_0=0.26;       %Vornockenhöhe, mm
H=5.5;          %Hauptnockenhöhe, mm
f_0=26;          %Abschnittlänge vom Vornocken, Grad
f_1=18;          % Grad
f_2=4;           % Grad
f_3=38;          % Grad





F_0=pi*f_0/180   ;  %Abschnittlänge vom Vornocken im Winkelmass, rad bei 26°
F_1=pi*f_1/180    ;%Abschnittlänge im Winkelmass, rad bei 18°
F_2=pi*f_2/180    ;%;Abschnittlänge im Winkelmass, rad bei 4°
F_3=pi*f_3/180    ;%Abschnittlänge im Winkelmass, rad bei 38°


phi_0=0:0.0001:F_0;
phi_1=F_0:0.0001:F_0+F_1;
phi_2=F_0+F_1:0.0001:F_0+F_1+F_2;
phi_3=F_0+F_1+F_2:0.0001:F_0+F_1+F_2+F_3;

phi=[phi_0 phi_1 phi_2 phi_3];

%Berechnung der Konstanten:
k_1=0.506606*F_2^2
k_2=0.9375*F_3^2
k_3=1.75*F_3
K_1=k_1+k_2+k_3*F_2
K_2=k_3+0.795775*F_2
ds_OE=H_0*pi/2/F_0                          %, mm
c_11=(K_1*ds_OE+K_2*H)/(K_1*2+K_2*F_1)      %, mm
c_12=(c_11-ds_OE)*F_1/pi                   %, mm
c_32=(2*c_11-ds_OE)/K_2                     %, mm
c_21=c_32*k_3                               %, mm
c_22=c_32*k_1                               %, mm
c_31=c_32/(16*F_3^2)                        %, mm
c_33=c_32*k_2                               %, mm
ds_max=2*c_11-ds_OE                         %, mm

%Abschnittfunktionen (f_s=F_3-(phi_3-F_0-F_1-F_2)):

f_s=F_3-(phi_3-F_0-F_1-F_2);
s_1=c_11*F_1-c_12*sin(pi/F_1*F_1)+H_0;       %, mm
s_2=c_21*F_2+c_22*sin(pi/2/F_2*F_2)+s_1;    %, mm


%Erhebung
s0=H_0*(1-cos(pi/2/F_0*phi_0)) ;              %, mm
s1=c_11*(phi_1-F_0)-c_12*sin(pi/F_1*(phi_1-F_0))+H_0;      %, mm
s2=c_21*(phi_2-F_0-F_1)+c_22*sin(pi/2/F_2*(phi_2-F_0-F_1))+s_1;     %, mm
s3=c_31*f_s.^4-c_32*f_s.^2+c_33+s_2;          %, mm

y=[s0 s1 s2 s3];
y1=R_G+y;
plot(phi*180/pi,y1);