# University of Florida/Eml4507/s13.team4ever.Wulterkens.R3.6

Problem 3.6 Determine Truss System Deformation

## Part 1
### Case A

EDU>> Edof=[1 1 2 3 4;2 3 4 5 6; 3 5 6 7 8;4 7 8 9 10;5 9 10 11 12; 6 11 12 13
14; 7 15 16 17 18;8 17 18 19 20;9 19 20 21 22; 10 21 22 23 24; 11 23 24 25 26;
12 25 26 27 28; 13 1 2 15 16;14 3 4 17 18; 15 5 6 19 20;16 7 8 21 22;17 9 10 23
24; 18 11 12 25 26; 19 13 14 27 28; 20 1 2 17 18; 21 3 4 19 20; 22 5 6 21 22; 23
7 8 23 24;24 9 10 25 26;25 11 12 27 28];
EDU>> L = 0.3;
EDU>> E = 100000000000;
EDU>> A = 0.0001;
EDU>> L2 = L*sqrt(2);
EDU>> K = zeros(28);
EDU>> f = zeros(28,1);
EDU>> f(13) = 10000; f(27) = 10000;
EDU>> ep = [E A];
EDU>> Ex = [0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;0 0; L L; 2*L 2*L;3*L 3*L;4*L 4*L;5*L 5*L;6*L 6*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L];
EDU>> Ey = [0 0;0 0;0 0;0 0;0 0;0 0;L L;L L;L L;L L;L L;L L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L];
EDU>> for i = 1:25
Ke = bar2e(Ex(i,:),Ey(i,:),ep);
K = assem(Edof(i,:),K,Ke);
end;
EDU>> bc = [1 0;2 0;15 0;16 0];
EDU>> [a,r]=solveq(K,f,bc)
a =
r =

### Case B

EDU>> Edof=[1 1 2 3 4;2 3 4 5 6; 3 5 6 7 8;4 7 8 9 10;5 9 10 11 12; 6 11 12 13 14; 7 15 16 17 18;8 17 18 19 20;9 19 20 21 22; 10 21 22 23 24; 11 23 24 25 26; 12 25 26 27 28; 13 1 2 15 16;14 3 4 17 18; 15 5 6 19 20;16 7 8 21 22;17 9 10 23 24; 18 11 12 25 26; 19 13 14 27 28; 20 1 2 17 18; 21 3 4 19 20; 22 5 6 21 22; 23 7 8 23 24;24 9 10 25 26;25 11 12 27 28];
EDU>> L = 0.3;
EDU>> E = 100000000000;
EDU>> A = 0.0001;
EDU>> L2 = L*sqrt(2);
EDU>> K = zeros(28);
EDU>> f = zeros(28,1);
EDU>> f(14) = 10000; f(28) = 10000;
EDU>> ep = [E A];
EDU>> Ex = [0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;0 0; L L; 2*L 2*L;3*L 3*L;4*L 4*L;5*L 5*L;6*L 6*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L];
EDU>> Ey = [0 0;0 0;0 0;0 0;0 0;0 0;L L;L L;L L;L L;L L;L L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L];
EDU>> for i = 1:25
Ke = bar2e(Ex(i,:),Ey(i,:),ep);
K = assem(Edof(i,:),K,Ke);
end;
EDU>> bc = [1 0;2 0;15 0;16 0];
EDU>> [a,r]=solveq(K,f,bc)
a =
r =

### Case c

EDU>> Edof=[1 1 2 3 4;2 3 4 5 6; 3 5 6 7 8;4 7 8 9 10;5 9 10 11 12; 6 11 12 13 14; 7 15 16 17 18;8 17 18 19 20;9 19 20 21 22; 10 21 22 23 24; 11 23 24 25 26; 12 25 26 27 28; 13 1 2 15 16;14 3 4 17 18; 15 5 6 19 20;16 7 8 21 22;17 9 10 23 24; 18 11 12 25 26; 19 13 14 27 28; 20 1 2 17 18; 21 3 4 19 20; 22 5 6 21 22; 23 7 8 23 24;24 9 10 25 26;25 11 12 27 28];
EDU>> L = 0.3;
EDU>> E = 100000000000;
EDU>> A = 0.0001;
EDU>> L2 = L*sqrt(2);
EDU>> K = zeros(28);
EDU>> f = zeros(28,1);
EDU>> f(13) = 10000; f(27) = -10000;
EDU>> ep = [E A];
EDU>> Ex = [0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;0 0; L L; 2*L 2*L;3*L 3*L;4*L 4*L;5*L 5*L;6*L 6*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L];
EDU>> Ey = [0 0;0 0;0 0;0 0;0 0;0 0;L L;L L;L L;L L;L L;L L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L];
EDU>> for i = 1:25
Ke = bar2e(Ex(i,:),Ey(i,:),ep);
K = assem(Edof(i,:),K,Ke);
end;
EDU>> bc = [1 0;2 0;15 0;16 0];
EDU>> [a,r]=solveq(K,f,bc)
a =
r =

## Part 3
### Case A

EDU>> Edof=[1 1 2 3 4;2 3 4 5 6; 3 5 6 7 8;4 7 8 9 10;5 9 10 11 12; 6 11 12 13 14; 7 13 14 15 16; 8 15 16 17 18; 9 19 20 21 22; 10 21 22 23 24; 11 23 24 25 26; 12 25 26 27 28; 13 27 28 29 30; 14 29 30 31 32; 15 31 32 33 34; 16 33 34 35 36; 17 1 2 19 20;18 3 4 21 22; 19 5 6 23 24; 20 7 8 25 26; 21 9 10 27 28; 22 11 12 29 30; 23 13 14 31 32; 24 15 16 33 34; 25 17 18 35 36; 26 1 2 21 22; 27 3 4 23 24; 28 5 6 25 26; 29 7 8 27 28; 30 9 10 29 30; 31 11 12 31 32; 32 13 14 33 34; 33 15 16 35 36]
Edof =
EDU>> L = 0.3;
EDU>> E = 100000000000;
EDU>> A = 0.0001;
EDU>> L2 = L*sqrt(2);
EDU>> K = zeros(36);
EDU>> f = zeros(36,1);
EDU>> f(17) = 10000; f(35) = 10000;
EDU>> ep = [E A];
EDU>> Ex = [0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L;0 0; L L; 2*L 2*L;3*L 3*L;4*L 4*L;5*L 5*L;6*L 6*L;7*L 7*L;8*L 8*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L];
EDU>> Ey = [0 0;0 0;0 0;0 0;0 0;0 0;0 0;0 0;L L;L L;L L;L L;L L;L L;L L;L L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L];
EDU>> for i = 1:33
Ke = bar2e(Ex(i,:),Ey(i,:),ep);
K = assem(Edof(i,:),K,Ke);
end;
EDU>> bc = [1 0;2 0;19 0;20 0];
EDU>> [a,r]=solveq(K,f,bc)
a =
r =

### Case B

EDU>> Edof=[1 1 2 3 4;2 3 4 5 6; 3 5 6 7 8;4 7 8 9 10;5 9 10 11 12; 6 11 12 13 14; 7 13 14 15 16; 8 15 16 17 18; 9 19 20 21 22; 10 21 22 23 24; 11 23 24 25 26; 12 25 26 27 28; 13 27 28 29 30; 14 29 30 31 32; 15 31 32 33 34; 16 33 34 35 36; 17 1 2 19 20;18 3 4 21 22; 19 5 6 23 24; 20 7 8 25 26; 21 9 10 27 28; 22 11 12 29 30; 23 13 14 31 32; 24 15 16 33 34; 25 17 18 35 36; 26 1 2 21 22; 27 3 4 23 24; 28 5 6 25 26; 29 7 8 27 28; 30 9 10 29 30; 31 11 12 31 32; 32 13 14 33 34; 33 15 16 35 36]
Edof =
EDU>> L = 0.3;
EDU>> E = 100000000000;
EDU>> A = 0.0001;
EDU>> L2 = L*sqrt(2);
EDU>> K = zeros(36);
EDU>> f = zeros(36,1);
EDU>> f(18) = 10000; f(36) = 10000;
EDU>> ep = [E A];
EDU>> Ex = [0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L;0 0; L L; 2*L 2*L;3*L 3*L;4*L 4*L;5*L 5*L;6*L 6*L;7*L 7*L;8*L 8*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L];
EDU>> Ey = [0 0;0 0;0 0;0 0;0 0;0 0;0 0;0 0;L L;L L;L L;L L;L L;L L;L L;L L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L];
EDU>> for i = 1:33
Ke = bar2e(Ex(i,:),Ey(i,:),ep);
K = assem(Edof(i,:),K,Ke);
end;
EDU>> bc = [1 0;2 0;19 0;20 0];
EDU>> [a,r]=solveq(K,f,bc)
a =
r =

### Case C

EDU>> Edof=[1 1 2 3 4;2 3 4 5 6; 3 5 6 7 8;4 7 8 9 10;5 9 10 11 12; 6 11 12 13 14; 7 13 14 15 16; 8 15 16 17 18; 9 19 20 21 22; 10 21 22 23 24; 11 23 24 25 26; 12 25 26 27 28; 13 27 28 29 30; 14 29 30 31 32; 15 31 32 33 34; 16 33 34 35 36; 17 1 2 19 20;18 3 4 21 22; 19 5 6 23 24; 20 7 8 25 26; 21 9 10 27 28; 22 11 12 29 30; 23 13 14 31 32; 24 15 16 33 34; 25 17 18 35 36; 26 1 2 21 22; 27 3 4 23 24; 28 5 6 25 26; 29 7 8 27 28; 30 9 10 29 30; 31 11 12 31 32; 32 13 14 33 34; 33 15 16 35 36]
Edof =
EDU>> L = 0.3;
EDU>> E = 100000000000;
EDU>> A = 0.0001;
EDU>> L2 = L*sqrt(2);
EDU>> K = zeros(36);
EDU>> f = zeros(36,1);
EDU>> f(17) = 10000; f(35) = -10000;
EDU>> ep = [E A];
EDU>> Ex = [0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L;0 0; L L; 2*L 2*L;3*L 3*L;4*L 4*L;5*L 5*L;6*L 6*L;7*L 7*L;8*L 8*L;0 L; L 2*L; 2*L 3*L;3*L 4*L; 4*L 5*L;5*L 6*L;6*L 7*L;7*L 8*L];
EDU>> Ey = [0 0;0 0;0 0;0 0;0 0;0 0;0 0;0 0;L L;L L;L L;L L;L L;L L;L L;L L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L;0 L];
EDU>> for i = 1:33
Ke = bar2e(Ex(i,:),Ey(i,:),ep);
K = assem(Edof(i,:),K,Ke);
end;
EDU>> bc = [1 0;2 0;19 0;20 0];
EDU>> [a,r]=solveq(K,f,bc)
a =
r =
EDU>>