rez 13 32 1 5 1

  • Published on
    07-Mar-2016

  • View
    215

  • Download
    0

DESCRIPTION

this is an useless file.

Transcript

  • % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

  • A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

  • % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

  • dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B vorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%TROUGLOVI%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%%

    trougloviSelektor=[1 2 4; 2 3 4; 1 3 4; 1 2 3];% [ABD,BCD,ACD,ABC]

    trougloviIDs = Elementi(:,[1 2 4]);

    for i=2:4

    trougloviIDs = [trougloviIDs; Elementi(:,trougloviSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, trougloviIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%%%%%%CENTRALNE

    LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, Elementi);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%%%%%%%%%%KRAJ

    MERDzOVANjA%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

  • AB_pola = (A+B)/2;

    % AB_pola = AB_pola*dt;

    % AB_pola = fix(AB_pola);

    % AB_pola = AB_pola/ dt;

    AC_pola = (A+C)/2;

    % AC_pola = AC_pola*dt;

    % AC_pola = fix(AC_pola);

    % AC_pola = AC_pola/ dt;

    AD_pola = (A+D)/2;

    % AD_pola = AD_pola*dt;

    % AD_pola = fix(AD_pola);

    % AD_pola = AD_pola/ dt;

    DB_pola = (D+B)/2;

    % DB_pola = DB_pola*dt;

    % DB_pola = fix(DB_pola);

    % DB_pola = DB_pola/ dt;

    = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

  • % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

  • C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

  • % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

  • arsMerdzovaniCvorovi=[arsMerdzovaniC

    % BC_pola = BC_pola*dt;

    % BC_pola = fix(BC_pola);

    % BC_pola = BC_pola/ dt;

    DC_pola = (D+C)/2;

    % DC_pola = DC_pola*dt;

    % DC_pola = fix(DC_pola);

    % DC_pola = DC_pola/ dt;

    %%%%%%%%%%%%%%%%

    tezisteTrouglaABC = (A+B+C)/3;

    % tezisteTrouglaABC = tezisteTrouglaABC*dt;

    % tezisteTrouglaABC = fix(tezisteTrouglaABC);

    % tezisteTrouglaABC = tezisteTrouglaABC/ dt;

    tezisteTrouglaBCD = (D+B+C)/3;

    % tezisteTrouglaBCD = tezisteTrouglaBCD*dt;

    % tezisteTrouglaBCD = fix(tezisteTrouglaBCD);

    % tezisteTrouglaBCD = tezisteTrouglaBCD/ dt;

    tezisteTrouglaADC = (A+D+C)/3;

    % tezisteTrouglaADC = tezisteTrouglaADC*dt;

    % tezisteTrouglaADC = fix(tezisteTrouglaADC);

    % tezisteTrouglaADC = tezisteTrouglaADC/ dt;

  • tezisteTrouglaABD = (A+B+D)/3;

    % tezisteTrouglaABD = tezisteTrouglaABD*dt;

    % tezisteTrouglaABD = fix(tezisteTrouglaABD);

    % tezisteTrouglaABD = tezisteTrouglaABD/ dt;

    C_pola(end)+brTetri] ;

    idDB_pola = [idAD_pola(end)+1:idAD_pola(end)+brTetri] ;

    idDC_pola = [idDB_pola(end)+1:idDB_pola(end)+brTetri] ;

    idtezisteTrouglaABD = [idDC_pola(end)+1:idDC_pola(end)+brTetri] ;

    idtezisteTrouglaBCD = [idtezisteTrouglaABD(end)+1:idtezisteTrouglaABD(end)+brTetri];

    clear('idtezisteTrouglaABD');

    clear('idcentarTetre');

    clear('A');

    clear('B');

    clear('C');

    clear('D');

    clear('AB_pola');

    clear('BC_pola');

    c

    materijaliPoHeksama=materijaliPoHeksama(indexHekse);

    rez.Hekse =

    Hekse;clear('Hekse');clear('BHekse');clear('indexHekse');clear('bckpHekse');clear('idOrigMerdzCvorovi');

    rez.materijaliPoHeksama=materijaliPoHeksama;clear('materijaliPoHeksama');

    %%%%%%%%%%%%%%%%%%%

    rez.Cvorovi = merzdovaniCvorovi; clear('merzdovaniCvorovi');

    %%% NADjI Obodne cvorove

  • selektor = [1 2 3 4;

    4 3 7 8;

    8 7 6 5;...3

    5 6 2 1;

    6 7 2 3;

    8 5 1 4];

    fejsovi=[0 0 0 0];

    %%%Hekse1 =Hekse(1:100,:);

    idElementFejs = [0];

    brElemenata = numel(rez.Hekse(:,1));

    for i = 1:6

    element = rez.Hekse(:,selektor(i,:)) ;

    fejsovi = [fejsovi;element(:,[1 2 3 4])];

    idElementFejs = [idElementFejs,1:brElemenata] ;

    f=unique(idElementFejs((i-1)*brElemenata+1:i*brElemenata));

    end%

    f

    indexElemenataNaZidu = idElementFejs(indexFejsovi);

    indexElemenataNaZidu = unique(indexElemenataNaZidu);

    clear('ia');clear('ia2');clear('indexBB');clear('indexIA');clear('index');clear('BB');clear('idElementFejs');

    upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(indexElemenataNaZidu,:), ['ElementiNaZidu3Dmreza.vtk']);

    rez.indexElemenataNaZidu= indexElemenataNaZidu ;

    rez.idFejsoviNaZidu = indexFejsovi ;

    fejsoviNaZidu = fejsovi(indexFejsovi,:) ;

    rez.fejsoviNaZidu = fejsoviNaZidu ;

  • %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %NACRTAJ FEJSOVE

    upisiSTL(rez.Cvorovi, Trouglovi,['obodSTLVene.stl']);

    % % %%%%%%

    % % upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(1:end,:), [filename '3Dmreza.vtk']);

    % % % arsVTK.saveVTKfromCvoroviElementi(rez.Cvorovi,rez.Hekse, '3Dmreza.vtk');

    end%func% rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

  • % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

  • % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

  • % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B vorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%TROUGLOVI%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%%

    trougloviSelektor=[1 2 4; 2 3 4; 1 3 4; 1 2 3];% [ABD,BCD,ACD,ABC]

    trougloviIDs = Elementi(:,[1 2 4]);

    for i=2:4

    trougloviIDs = [trougloviIDs; Elementi(:,trougloviSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, trougloviIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

  • %%%%%%%%%%%%%%%%%%%%CENTRALNE

    LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, Elementi);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%%%%%%%%%%KRAJ

    MERDzOVANjA%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    AB_pola = (A+B)/2;

    % AB_pola = AB_pola*dt;

    % AB_pola = fix(AB_pola);

    % AB_pola = AB_pola/ dt;

    AC_pola = (A+C)/2;

    % AC_pola = AC_pola*dt;

    % AC_pola = fix(AC_pola);

    % AC_pola = AC_pola/ dt;

    AD_pola = (A+D)/2;

    % AD_pola = AD_pola*dt;

    % AD_pola = fix(AD_pola);

    % AD_pola = AD_pola/ dt;

    DB_pola = (D+B)/2;

    % DB_pola = DB_pola*dt;

    % DB_pola = fix(DB_pola);

  • % DB_pola = DB_pola/ dt;

    = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

  • end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

  • for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

  • linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC

    % BC_pola = BC_pola*dt;

    % BC_pola = fix(BC_pola);

    % BC_pola = BC_pola/ dt;

    DC_pola = (D+C)/2;

    % DC_pola = DC_pola*dt;

    % DC_pola = fix(DC_pola);

    % DC_pola = DC_pola/ dt;

    %%%%%%%%%%%%%%%%

    tezisteTrouglaABC = (A+B+C)/3;

    % tezisteTrouglaABC = tezisteTrouglaABC*dt;

    % tezisteTrouglaABC = fix(tezisteTrouglaABC);

    % tezisteTrouglaABC = tezisteTrouglaABC/ dt;

    tezisteTrouglaBCD = (D+B+C)/3;

    % tezisteTrouglaBCD = tezisteTrouglaBCD*dt;

    % tezisteTrouglaBCD = fix(tezisteTrouglaBCD);

  • % tezisteTrouglaBCD = tezisteTrouglaBCD/ dt;

    tezisteTrouglaADC = (A+D+C)/3;

    % tezisteTrouglaADC = tezisteTrouglaADC*dt;

    % tezisteTrouglaADC = fix(tezisteTrouglaADC);

    % tezisteTrouglaADC = tezisteTrouglaADC/ dt;

    tezisteTrouglaABD = (A+B+D)/3;

    % tezisteTrouglaABD = tezisteTrouglaABD*dt;

    % tezisteTrouglaABD = fix(tezisteTrouglaABD);

    % tezisteTrouglaABD = tezisteTrouglaABD/ dt;

    C_pola(end)+brTetri] ;

    idDB_pola = [idAD_pola(end)+1:idAD_pola(end)+brTetri] ;

    idDC_pola = [idDB_pola(end)+1:idDB_pola(end)+brTetri] ;

    idtezisteTrouglaABD = [idDC_pola(end)+1:idDC_pola(end)+brTetri] ;

    idtezisteTrouglaBCD = [idtezisteTrouglaABD(end)+1:idtezisteTrouglaABD(end)+brTetri];

    clear('idtezisteTrouglaABD');

    clear('idcentarTetre');

    clear('A');

    clear('B');

    clear('C');

    clear('D');

    clear('AB_pola');

    clear('BC_pola');

  • c

    materijaliPoHeksama=materijaliPoHeksama(indexHekse);

    rez.Hekse =

    Hekse;clear('Hekse');clear('BHekse');clear('indexHekse');clear('bckpHekse');clear('idOrigMerdzCvorovi');

    rez.materijaliPoHeksama=materijaliPoHeksama;clear('materijaliPoHeksama');

    %%%%%%%%%%%%%%%%%%%

    rez.Cvorovi = merzdovaniCvorovi; clear('merzdovaniCvorovi');

    %%% NADjI Obodne cvorove

    selektor = [1 2 3 4;

    4 3 7 8;

    8 7 6 5;...3

    5 6 2 1;

    6 7 2 3;

    8 5 1 4];

    fejsovi=[0 0 0 0];

    %%%Hekse1 =Hekse(1:100,:);

    idElementFejs = [0];

    brElemenata = numel(rez.Hekse(:,1));

    for i = 1:6

    element = rez.Hekse(:,selektor(i,:)) ;

    fejsovi = [fejsovi;element(:,[1 2 3 4])];

    idElementFejs = [idElementFejs,1:brElemenata] ;

    f=unique(idElementFejs((i-1)*brElemenata+1:i*brElemenata));

    end%

    f

    indexElemenataNaZidu = idElementFejs(indexFejsovi);

  • indexElemenataNaZidu = unique(indexElemenataNaZidu);

    clear('ia');clear('ia2');clear('indexBB');clear('indexIA');clear('index');clear('BB');clear('idElementFejs');

    upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(indexElemenataNaZidu,:), ['ElementiNaZidu3Dmreza.vtk']);

    rez.indexElemenataNaZidu= indexElemenataNaZidu ;

    rez.idFejsoviNaZidu = indexFejsovi ;

    fejsoviNaZidu = fejsovi(indexFejsovi,:) ;

    rez.fejsoviNaZidu = fejsoviNaZidu ;

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %NACRTAJ FEJSOVE

    upisiSTL(rez.Cvorovi, Trouglovi,['obodSTLVene.stl']);

    % % %%%%%%

    % % upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(1:end,:), [filename '3Dmreza.vtk']);

    % % % arsVTK.saveVTKfromCvoroviElementi(rez.Cvorovi,rez.Hekse, '3Dmreza.vtk');

    end%func% rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

  • % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

  • % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

  • % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B vorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%TROUGLOVI%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%%

    trougloviSelektor=[1 2 4; 2 3 4; 1 3 4; 1 2 3];% [ABD,BCD,ACD,ABC]

    trougloviIDs = Elementi(:,[1 2 4]);

  • for i=2:4

    trougloviIDs = [trougloviIDs; Elementi(:,trougloviSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, trougloviIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%%%%%%CENTRALNE

    LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, Elementi);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%%%%%%%%%%KRAJ

    MERDzOVANjA%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    AB_pola = (A+B)/2;

    % AB_pola = AB_pola*dt;

    % AB_pola = fix(AB_pola);

    % AB_pola = AB_pola/ dt;

    AC_pola = (A+C)/2;

    % AC_pola = AC_pola*dt;

    % AC_pola = fix(AC_pola);

    % AC_pola = AC_pola/ dt;

    AD_pola = (A+D)/2;

  • % AD_pola = AD_pola*dt;

    % AD_pola = fix(AD_pola);

    % AD_pola = AD_pola/ dt;

    DB_pola = (D+B)/2;

    % DB_pola = DB_pola*dt;

    % DB_pola = fix(DB_pola);

    % DB_pola = DB_pola/ dt;

    = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

  • arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

  • idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

  • % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC

    % BC_pola = BC_pola*dt;

    % BC_pola = fix(BC_pola);

    % BC_pola = BC_pola/ dt;

    DC_pola = (D+C)/2;

    % DC_pola = DC_pola*dt;

    % DC_pola = fix(DC_pola);

    % DC_pola = DC_pola/ dt;

    %%%%%%%%%%%%%%%%

    tezisteTrouglaABC = (A+B+C)/3;

  • % tezisteTrouglaABC = tezisteTrouglaABC*dt;

    % tezisteTrouglaABC = fix(tezisteTrouglaABC);

    % tezisteTrouglaABC = tezisteTrouglaABC/ dt;

    tezisteTrouglaBCD = (D+B+C)/3;

    % tezisteTrouglaBCD = tezisteTrouglaBCD*dt;

    % tezisteTrouglaBCD = fix(tezisteTrouglaBCD);

    % tezisteTrouglaBCD = tezisteTrouglaBCD/ dt;

    tezisteTrouglaADC = (A+D+C)/3;

    % tezisteTrouglaADC = tezisteTrouglaADC*dt;

    % tezisteTrouglaADC = fix(tezisteTrouglaADC);

    % tezisteTrouglaADC = tezisteTrouglaADC/ dt;

    tezisteTrouglaABD = (A+B+D)/3;

    % tezisteTrouglaABD = tezisteTrouglaABD*dt;

    % tezisteTrouglaABD = fix(tezisteTrouglaABD);

    % tezisteTrouglaABD = tezisteTrouglaABD/ dt;

    C_pola(end)+brTetri] ;

    idDB_pola = [idAD_pola(end)+1:idAD_pola(end)+brTetri] ;

    idDC_pola = [idDB_pola(end)+1:idDB_pola(end)+brTetri] ;

    idtezisteTrouglaABD = [idDC_pola(end)+1:idDC_pola(end)+brTetri] ;

    idtezisteTrouglaBCD = [idtezisteTrouglaABD(end)+1:idtezisteTrouglaABD(end)+brTetri];

    clear('idtezisteTrouglaABD');

  • clear('idcentarTetre');

    clear('A');

    clear('B');

    clear('C');

    clear('D');

    clear('AB_pola');

    clear('BC_pola');

    c

    materijaliPoHeksama=materijaliPoHeksama(indexHekse);

    rez.Hekse =

    Hekse;clear('Hekse');clear('BHekse');clear('indexHekse');clear('bckpHekse');clear('idOrigMerdzCvorovi');

    rez.materijaliPoHeksama=materijaliPoHeksama;clear('materijaliPoHeksama');

    %%%%%%%%%%%%%%%%%%%

    rez.Cvorovi = merzdovaniCvorovi; clear('merzdovaniCvorovi');

    %%% NADjI Obodne cvorove

    selektor = [1 2 3 4;

    4 3 7 8;

    8 7 6 5;...3

    5 6 2 1;

    6 7 2 3;

    8 5 1 4];

    fejsovi=[0 0 0 0];

    %%%Hekse1 =Hekse(1:100,:);

    idElementFejs = [0];

    brElemenata = numel(rez.Hekse(:,1));

    for i = 1:6

  • element = rez.Hekse(:,selektor(i,:)) ;

    fejsovi = [fejsovi;element(:,[1 2 3 4])];

    idElementFejs = [idElementFejs,1:brElemenata] ;

    f=unique(idElementFejs((i-1)*brElemenata+1:i*brElemenata));

    end%

    f

    indexElemenataNaZidu = idElementFejs(indexFejsovi);

    indexElemenataNaZidu = unique(indexElemenataNaZidu);

    clear('ia');clear('ia2');clear('indexBB');clear('indexIA');clear('index');clear('BB');clear('idElementFejs');

    upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(indexElemenataNaZidu,:), ['ElementiNaZidu3Dmreza.vtk']);

    rez.indexElemenataNaZidu= indexElemenataNaZidu ;

    rez.idFejsoviNaZidu = indexFejsovi ;

    fejsoviNaZidu = fejsovi(indexFejsovi,:) ;

    rez.fejsoviNaZidu = fejsoviNaZidu ;

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %NACRTAJ FEJSOVE

    upisiSTL(rez.Cvorovi, Trouglovi,['obodSTLVene.stl']);

    % % %%%%%%

    % % upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(1:end,:), [filename '3Dmreza.vtk']);

    % % % arsVTK.saveVTKfromCvoroviElementi(rez.Cvorovi,rez.Hekse, '3Dmreza.vtk');

    end%func% rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

  • brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

  • B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

  • function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

    % A = A/ dt;

    B = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    % rez.idFejsoviNaZidu

    % rez.Cvorovi

    % rez.Hekse

    function [rez] = Backup_of_prebaciIzTetriUHekse(Cvorovi, Elementi,materijaliIDs,filename)

    brTetri = numel(Elementi(:,1)) ;

    materijaliPoTetrama = materijaliIDs;

    %%%%%%%%%%%%%%

    dt=100;

    A = Cvorovi(Elementi(:,1),:);

    % A = A*dt;

    % A = fix(A);

  • % A = A/ dt;

    B vorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%TROUGLOVI%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%%

    trougloviSelektor=[1 2 4; 2 3 4; 1 3 4; 1 2 3];% [ABD,BCD,ACD,ABC]

    trougloviIDs = Elementi(:,[1 2 4]);

    for i=2:4

    trougloviIDs = [trougloviIDs; Elementi(:,trougloviSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, trougloviIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%%%%%%CENTRALNE

    LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, Elementi);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniCvorovi;merdzovanniCvorovi];

    arsIdMerdzCvorovi =[arsIdMerdzCvorovi;idMerdzovanihCvorovaIzID+brCvorova];

    %%%%%%%%%%%%%%%%%%%%%%%%KRAJ

    MERDzOVANjA%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    AB_pola = (A+B)/2;

    % AB_pola = AB_pola*dt;

    % AB_pola = fix(AB_pola);

    % AB_pola = AB_pola/ dt;

  • AC_pola = (A+C)/2;

    % AC_pola = AC_pola*dt;

    % AC_pola = fix(AC_pola);

    % AC_pola = AC_pola/ dt;

    AD_pola = (A+D)/2;

    % AD_pola = AD_pola*dt;

    % AD_pola = fix(AD_pola);

    % AD_pola = AD_pola/ dt;

    DB_pola = (D+B)/2;

    % DB_pola = DB_pola*dt;

    % DB_pola = fix(DB_pola);

    % DB_pola = DB_pola/ dt;

    = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

  • D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

    % C = C*dt;

    % C = fix(C);

  • % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC = Cvorovi(Elementi(:,2),:);

    % B = B*dt;

    % B = fix(B);

    % B = B/ dt;

    C = Cvorovi(Elementi(:,3),:);

  • % C = C*dt;

    % C = fix(C);

    % C = C/ dt;

    D = Cvorovi(Elementi(:,4),:);

    % D = D*dt;

    % D = fix(D);

    % D = D/ dt;

    %merdzuj A B C D

    idABCD = [Elementi(:,1);Elementi(:,2);Elementi(:,3);Elementi(:,4)];

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, idABCD);

    arsMerdzovaniCvorovi=merdzovanniCvorovi;

    arsIdMerdzCvorovi =idMerdzovanihCvorovaIzID;

    %%%%%%%%%%%%%%LINIJE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %%%%%%%%%%%%%

    linijeSelektor=[1 2; 2 3; 1 3; 1 4; 2 4; 3 4];% [AB,BC,AC,AD,BD,DC]

    linijeIDs = Elementi(:,[1 2]);

    for i=2:6

    linijeIDs = [linijeIDs; Elementi(:,linijeSelektor(i,:))];

    end

    [merdzovanniCvorovi, idMerdzovanihCvorovaIzID] = merdzujCvorove(Cvorovi, linijeIDs);

    brCvorova = numel(arsMerdzovaniCvorovi(:,1));

    arsMerdzovaniCvorovi=[arsMerdzovaniC

    % BC_pola = BC_pola*dt;

    % BC_pola = fix(BC_pola);

    % BC_pola = BC_pola/ dt;

  • DC_pola = (D+C)/2;

    % DC_pola = DC_pola*dt;

    % DC_pola = fix(DC_pola);

    % DC_pola = DC_pola/ dt;

    %%%%%%%%%%%%%%%%

    tezisteTrouglaABC = (A+B+C)/3;

    % tezisteTrouglaABC = tezisteTrouglaABC*dt;

    % tezisteTrouglaABC = fix(tezisteTrouglaABC);

    % tezisteTrouglaABC = tezisteTrouglaABC/ dt;

    tezisteTrouglaBCD = (D+B+C)/3;

    % tezisteTrouglaBCD = tezisteTrouglaBCD*dt;

    % tezisteTrouglaBCD = fix(tezisteTrouglaBCD);

    % tezisteTrouglaBCD = tezisteTrouglaBCD/ dt;

    tezisteTrouglaADC = (A+D+C)/3;

    % tezisteTrouglaADC = tezisteTrouglaADC*dt;

    % tezisteTrouglaADC = fix(tezisteTrouglaADC);

    % tezisteTrouglaADC = tezisteTrouglaADC/ dt;

    tezisteTrouglaABD = (A+B+D)/3;

    % tezisteTrouglaABD = tezisteTrouglaABD*dt;

    % tezisteTrouglaABD = fix(tezisteTrouglaABD);

    % tezisteTrouglaABD = tezisteTrouglaABD/ dt;

  • C_pola(end)+brTetri] ;

    idDB_pola = [idAD_pola(end)+1:idAD_pola(end)+brTetri] ;

    idDC_pola = [idDB_pola(end)+1:idDB_pola(end)+brTetri] ;

    idtezisteTrouglaABD = [idDC_pola(end)+1:idDC_pola(end)+brTetri] ;

    idtezisteTrouglaBCD = [idtezisteTrouglaABD(end)+1:idtezisteTrouglaABD(end)+brTetri];

    clear('idtezisteTrouglaABD');

    clear('idcentarTetre');

    clear('A');

    clear('B');

    clear('C');

    clear('D');

    clear('AB_pola');

    clear('BC_pola');

    c

    materijaliPoHeksama=materijaliPoHeksama(indexHekse);

    rez.Hekse =

    Hekse;clear('Hekse');clear('BHekse');clear('indexHekse');clear('bckpHekse');clear('idOrigMerdzCvorovi');

    rez.materijaliPoHeksama=materijaliPoHeksama;clear('materijaliPoHeksama');

    %%%%%%%%%%%%%%%%%%%

    rez.Cvorovi = merzdovaniCvorovi; clear('merzdovaniCvorovi');

    %%% NADjI Obodne cvorove

    selektor = [1 2 3 4;

    4 3 7 8;

    8 7 6 5;...3

    5 6 2 1;

  • 6 7 2 3;

    8 5 1 4];

    fejsovi=[0 0 0 0];

    %%%Hekse1 =Hekse(1:100,:);

    idElementFejs = [0];

    brElemenata = numel(rez.Hekse(:,1));

    for i = 1:6

    element = rez.Hekse(:,selektor(i,:)) ;

    fejsovi = [fejsovi;element(:,[1 2 3 4])];

    idElementFejs = [idElementFejs,1:brElemenata] ;

    f=unique(idElementFejs((i-1)*brElemenata+1:i*brElemenata));

    end%

    f

    indexElemenataNaZidu = idElementFejs(indexFejsovi);

    indexElemenataNaZidu = unique(indexElemenataNaZidu);

    clear('ia');clear('ia2');clear('indexBB');clear('indexIA');clear('index');clear('BB');clear('idElementFejs');

    upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(indexElemenataNaZidu,:), ['ElementiNaZidu3Dmreza.vtk']);

    rez.indexElemenataNaZidu= indexElemenataNaZidu ;

    rez.idFejsoviNaZidu = indexFejsovi ;

    fejsoviNaZidu = fejsovi(indexFejsovi,:) ;

    rez.fejsoviNaZidu = fejsoviNaZidu ;

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    %NACRTAJ FEJSOVE

    upisiSTL(rez.Cvorovi, Trouglovi,['obodSTLVene.stl']);

  • % % %%%%%%

    % % upisiVTK_HEXE(rez.Cvorovi,rez.Hekse(1:end,:), [filename '3Dmreza.vtk']);

    % % % arsVTK.saveVTKfromCvoroviElementi(rez.Cvorovi,rez.Hekse, '3Dmreza.vtk');

    end%func

Recommended

View more >