clear all
    close all
    mueR = [0.1 0.14 0.14 0.14 0.14 0.14 0.14 0.14];                          % Gleitreibwert ?Gleit Fügefläche innen [-]
    mueH = [0.12 0.17 0.17 0.17 0.17 0.17 0.17 0.17];                         % Haftreibwert ?Haft Fügefläche innen [-]
%     Da=[21 28.01 36.014 45.008 54.028 70.018];                          % Aussendurchmesser Da[mm]
    %     Da=[32.003 39.994 47.992 ];
    Pi=pi;                                                                    % Kreiszahl [-]
    %     Di=[20.898	31.845	39.93 47.992 ];
%     Di=[11.999	20.950	27.920	35.918	44.913	53.916 70.018];               % Innendurchmesser [mm]
    E=[530 208 208 208 208 208 208 208]*10^3;                                 % Elastizitätsmodul[N/mm^2]
    nu=[0.23 0.28 0.28 0.28 0.28 0.28 0.28 0.28];                             % Querdehnzahl [-]
    anz =7;                                                                  % Anzahl der kern + Ring !!!!
    n=[];
    for  i=2:anz   % i=anz:-1:2
        n(i-1)=i;
    end
    n=n(end:-1:1);
    gFa=[1 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2];                              % Glättungsfaktor gF außen [-]
    gFi=[0 0.5 0.2 0.2 0.8 0.2 0.2 0.2 0.2 0.2];                              % Glättungsfaktor gF innen [-]
    Rza=[3 3 3 3 3 3 3 3 3];                                                  % Oberflächenrauheit Rza außen [-]
    Rzi=[0 3 3 3 3 3 3 3 3];                                                  % Oberflächenrauheit Rzi innen [-]
    delta_r=0;                                                                % Bereücksichtigung von Glättung 'Ohne=0 & Mit=1'
    m= [0.000001563 -5.09e-10 0 0 ];                                          % Parameter der Polynomfunktion (Fließkurve)
    b= -6.13e-5;
    
%     lamda=[96 39.88 39.88 39.88 39.88 39.88 39.88 39.88 39.88 39.88 39.88];  %[0.01 0.035] % [96 40 40 40 40 40 40 40 40 40 40 40 40 40 40];                     % Wärmeleitfähigkeit [W/m.K]
    alpha=[1e16 50]; %[13*10^-6,13*10^-6] %[0.0000048 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114]; %[13*10^-6,13*10^-6];
    R=[5.9995	10.475	13.96	17.959	22.4565	26.958	35.009];%[5.9995	10.4750	13.9600	17.9590	22.4565	26.9580	35.0090];
    gamma = [0.0000048 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114 0.0000114];
    l=0.1;  %103;                                                                    % fiktive Matrizenlänge 1 [m], kürzt sich später raus [mm]
    TiI=200;   %200+273.15                                                             % Innere Temperatur [K]
    TaA=80;
%     D=[26.75 45 54 70 83 100.1 110 100.008 130.015 150 200];    % Di(1) + Da Innere Durchmesser der Kern + Alle Außere Durchmesser + Da9
    lamda_kern =[0 0 -0.035 96]; %[96 -0.035 0 0] % Wärmeleitfähigkeit [W/(m*K)] Kern
    lamda_Ring = [-3.019e-10 4.941e-7 -0.0002965  0.051 39.88 ]; % Wärmeleitfähigkeit [W/(m*K)] Kern
    lamda_W = 37;   % Wärmeleitfähigkeit [W/(m*K)] Werzeug-Aufnahme
     Da=[24 34.01 40.99 50.22 70 80 100.2 120];
      Di=[10.82 24 34.01 40.99 50.22 70 80 100];
       R=[12 17.0050 20.4950 25.1100 35 40 50.1000 60];
       D=[12 24 34.01 40.99 50.22 70 80 100.2 120];
       lamda=[96	39.8800	39.8800	39.8800	39.8800	39.8800	39.8800	37];
         
       lamda_w=37;
    % Mantelfläche
  
    for j=1:anz+1
        
        
        for i=1:anz
            
            Ai(i) = pi*D(i)/1000*l;
            
        end
        
        for i=1:anz+1 %anz+2
            
            T_Start(i) = 1;
            
            
        end
        
        for i=1:anz+2 %anz+2
            
            T_(i) = 0;
            
        end
        
        for i = 1:anz+1
            
            T(i) = T_(i+1) - (T_(i+1)-T_(i))/2 ;
            
        end
        
        while or(T_Start(j) > T(j) + 0.5 , T_Start(j) < T(j) - 0.5)
            
            T_Start(j) = T(j);
            
            % Wärmeleitfähigkeit [W/(m*K)] Kern bei vorliegender Temperatur im Kern
            
            lamda_kern = polyval (lamda_kern,T(1));
            
            for i=1:anz-1
                
                lamda_Ring(i) = polyval (lamda_Ring,T(i+1)) ;
                
            end
            
            lamda = [lamda_kern lamda_Ring lamda_w]
%             lamda=[96	39.8800	39.8800	39.8800	39.8800	39.8800	39.8800	37];
            
            for i = 1:anz
                
                Summand(i) = 1/(2*pi*l*lamda(i)*log(D(i+1)/D(i)));
                
            end
            
            
            Summand_i= 1 / (alpha(1) * Ai(1));
            %Summand_i= 0;
            Summand_a= 1 / (alpha(2) * Ai(end));
            % Summand_a= 0;
            
            Summand_w= 1/(2*pi*l*lamda_w*log(D(end)/D(anz)));
            
            Summand = [Summand_i Summand  Summand_a  Summand_w];
            
            Sum = sum(Summand);
            
            % Wärmestrom [W]
            
            Q = (TiI - TaA)/(Sum);
            
            % Temperaturen an Fugendurchmessern [°C]
            %          T_(1) = TiI ;
            T_(1) = TiI - Q / (alpha(1) * Ai(1));
            for i=1:anz+1
                
                T_ (i+1) = T_(i) - Q/(2*pi*l*lamda(i)*log(D(i+1)/D(i)));  % Außer Ta0 da lamda_kern = 0
                
            end
            
            TaA = T_(end) - Q / (alpha(2)*Ai(end)) ;  % mit Konvektion an Matrizen-Außendurchmesser oder TaA=T_(end)=Ta9_
            %          TaA = T_(end);
            for i = 1:anz+1
                
                T(i) = T_(i+1) - (T_(i+1)-T_(i))/2 ;
                
            end
            
        end
        %      Ta_mk(1) = TiI
        %     Ta_mk(j+1) = T_(j+1)
        Ta_mk(j) = T_(j)
        
    end