% Erstellen der Verzeichnisstruktur
for n=1:length(m)
    node=num2str(m(n));
    mkdir('D:\sim\',node);
    cd(['D:\sim\',node]);
    copyfile('D:\ordner1\ordner2\kraft.key');
    copyfile('D:\ordner1\ordner2\modellcharakteristik.key');

    % Angabe des zu belasteten Knoten
    full_directory='kraft.key';
    fid = fopen((full_directory),'a');

      fprintf(fid,'     %s        3          2      -1.0         0',node);
      fprintf(fid,'\n*END');

    fclose(fid);

    % Vervollständigen der Datei D:\ordner1\ordner2\modellcharakteristik.key.key
    full_directory='D:\ordner1\ordner2\modellcharakteristik.key.key';
    fid = fopen((full_directory),'a');

      % NODOUT erzeugen
      fprintf(fid,'\n*DATABASE_HISTORY_NODE');
      fprintf(fid,'\n$HNAME OUTPUTBLOCKS       1nodout');
      fprintf(fid,'\n     %s',node);
      % Erzeuge *INCLUDE Befehl
      fprintf(fid,'\n*INCLUDE');
      fprintf(fid,'\nd:\\ordner1\\%s\\kraft.key',node);
      fprintf(fid,'\n*END');

    fclose(fid);


end


%---------------------------------------
%Auslesen der nodout Dateien
for n=1:length(m)
    node=num2str(m(n));
    full_directory=['D:\ordner1\',node,'\nodout'];

    % Öffen der Datei nodout
    fid = fopen(strcat(full_directory),'r');
    count1=0;

    % Messwerte einlesen
    while ~feof(fid)
        zeile = fgetl(fid);

        % nach Zeitschritt suchen
        if strfind(zeile,'time')
            strtime=zeile(106:116);
            time2=str2double(strtime);
        end

        werte = str2double(zeile);

        if numel(werte)~=0

            node2=str2double(node);

            if werte(1)==node2

                % translatorische Werte einlesen
                if numel(werte)==13
                    count1 = count1 + 1;
                    time(count1)=time2;
                    trans(count1,:)=[time2 werte];
                    save(['D:\sim\',node,'\read_nodout'],'trans','-ascii');

                    % Einlesen der verwendeten Kraft
                    full_directory=['d:\ordner1\',node,'\kraft.key'];
                    fid2 = fopen((full_directory),'r');

                    while ~feof(fid2)
                        zeile = fgetl(fid2);

                        % nach Kraft suchen
                        if strfind(zeile,'$HMCURVE')
                            textscan(fid2, '%*[^\n]', 3);
                            force=fscanf(fid2,'%*f %f');

                        end

                    end
                    fclose(fid2);

                    % Berechnung der Deformationsenergie
                    defenerg_1=[trans(:,1),trans(:,2),((trans(:,3).^2+trans(:,4).^2+trans(:,5).^2).^0.5).*force]; %Energie für Einzelknoten pro Zeitschritt
                    save(['D:\sim\',node,'\defenerg_1'],'defenerg_1','-ascii');
                end

            end

        end

    end
    fclose(fid);
    clear time trans defenerg;
end
%%