bottennoder för referens är 1 6 12 19 27 35 43 51 59 67 75 83 91 98 104


clear all

Coord = [ 0 0; 0 5; 0 10; 0 15; 0 20; 20/3 0; 20/3 5; 20/3 10; 20/3 15; 20/3 20; 20/3 25; 40/3 0; 40/3 5; 40/3 10; 40/3 15; 40/3 20; 40/3 25; 40/3 30; 20 0; 20 5; 20 10; 20 15; 20 20; 20 25; 20 30; 20 35; 25 0; 25 5; 25 10; 25 15; 25 20; 25 25; 25 30; 25 35; 30 0; 30 5; 30 10; 30 15; 30 20; 30 25; 30 30; 30 35; 35 0; 35 5; 35 10; 35 15; 35 20; 35 25; 35 30; 35 35; 40 0; 40 5; 40 10; 40 15; 40 20; 40 25; 40 30; 40 35; 45 0; 45 5; 45 10; 45 15; 45 20; 45 25; 45 30; 45 35; 50 0; 50 5; 50 10; 50 15; 50 20; 50 25; 50 30; 50 35; 55 0; 55 5; 55 10; 55 15; 55 20; 55 25; 55 30; 55 35; 60 0; 60 5; 60 10; 60 15; 60 20; 60 25; 60 30; 60 35; 200/3 0; 200/3 5; 200/3 10; 200/3 15; 200/3 20; 200/3 25; 200/3 30; 220/3 0; 220/3 5; 220/3 10; 220/3 15; 220/3 20; 220/3 25; 240/3 0; 240/3 5; 240/3 10; 240/3 15; 240/3 20 ];

plot(Coord(:,1), Coord(:,2),'x')

%konstanter
kasf = 30;
casf = 1200;

ep = 30; %mm
D = kasf*eye(2);

KE = zeros(3,3,178);
Edof = zeros(4,178);
Dof = [1:108]';


l = 0;
for j=1:49,
    i = Coord(j,1);


    if (i == 0)
        nodejump = 5;
    elseif (i == 20/3)
        nodejump = 6;
    elseif (i == 40/3)
        nodejump = 7;
    else
        nodejump = 8;
    end







            if (Coord(j,2) < Coord(j+nodejump+1,2))
                 l = l + 1

                x1 = Coord(j,1);
                y1 = Coord(j,2);

                x2 = Coord(j+nodejump,1);
                y2 = Coord(j+nodejump,2);

                x3 = Coord(j+nodejump+1,1);
                y3 = Coord(j+nodejump+1,2);


                Edof(:,l) = [l j j+nodejump j+nodejump+1 ];


                ex = [x1 x2 x3];
                ey = [y1 y2 y3];


                KE(:,:,l) = flw2te(ex,ey,ep,D);

            end

            if (Coord(j,2) < Coord(j+1,2))
                l = l + 1

                x1 = Coord(j,1);
                y1 = Coord(j,2);

                x2 = Coord(j+1,1);
                y2 = Coord(j+1,2);

                x3 = Coord(j+nodejump+1,1);
                y3 = Coord(j+nodejump+1,2);


                Edof(:,l) = [l j j+1 j+nodejump+1 ];


                ex = [x1 x2 x3];
                ey = [y1 y2 y3];


                KE(:,:,l) = flw2te(ex,ey,ep,D);

            end


end


l = 89;
for j=51:102,
    i = Coord(j,1);

    if (i == 200/3)
        nodejump = 7;
    elseif (i == 220/3)
        nodejump = 6;
    else
        nodejump = 8;
    end


            if (Coord(j,2) < Coord(j+1,2))
                 l = l + 1

                x1 = Coord(j,1);
                y1 = Coord(j,2);

                x2 = Coord(j+nodejump,1);
                y2 = Coord(j+nodejump,2);

                x3 = Coord(j+1,1);
                y3 = Coord(j+1,2);


                Edof(:,l) = [l j j+nodejump j+1 ];


                ex = [x1 x2 x3];
                ey = [y1 y2 y3];


                KE(:,:,l) = flw2te(ex,ey,ep,D);

            end

            if ((j < 102) && Coord(j,2) < Coord(j+nodejump+1,2))
                l = l + 1

                x1 = Coord(j+nodejump+1,1);
                y1 = Coord(j+nodejump+1,2);

                x2 = Coord(j+1,1);
                y2 = Coord(j+1,2);

                x3 = Coord(j+nodejump,1);
                y3 = Coord(j+nodejump,2);


                Edof(:,l) = [l j+nodejump+1 j+1 j+nodejump ];


                ex = [x1 x2 x3];
                ey = [y1 y2 y3];


                KE(:,:,l) = flw2te(ex,ey,ep,D);

            end


end
Edof = Edof'
[Ex,Ey] = coordxtr(Edof,Coord,Dof,3)