View difference between Paste ID: <a href="/gE5GjLYv">gE5GjLYv</a> and <a href="/ajFR7gKP">ajFR7gKP</a>

nodes = [ 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 ]
1	-	nodes = [ 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 ]
1	+	bottennoder för referens är 1 6 12 19 27 35 43 51 59 67 75 83 91 98 104
2
3	-	bottennoder för referens är 1 6 12 19 27 35 43 51 59 67 75 83 91 98 104
3	+
4		clear all
5
6		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 ];
7
8		plot(Coord(:,1), Coord(:,2),'x')
9
10		%konstanter
11		kasf = 30;
12		casf = 1200;
13
14		ep = 30; %mm
15		D = kasf*eye(2);
16
17		KE = zeros(3,3,178);
18		Edof = zeros(4,178);
19		Dof = [1:108]';
20
21
22		l = 0;
23		for j=1:49,
24		i = Coord(j,1);
25
26
27		if (i == 0)
28		nodejump = 5;
29		elseif (i == 20/3)
30		nodejump = 6;
31		elseif (i == 40/3)
32		nodejump = 7;
33		else
34		nodejump = 8;
35		end
36
37
38
39
40
41
42
43		if (Coord(j,2) < Coord(j+nodejump+1,2))
44		l = l + 1
45
46		x1 = Coord(j,1);
47		y1 = Coord(j,2);
48
49		x2 = Coord(j+nodejump,1);
50		y2 = Coord(j+nodejump,2);
51
52		x3 = Coord(j+nodejump+1,1);
53		y3 = Coord(j+nodejump+1,2);
54
55
56		Edof(:,l) = [l j j+nodejump j+nodejump+1 ];
57
58
59		ex = [x1 x2 x3];
60		ey = [y1 y2 y3];
61
62
63		KE(:,:,l) = flw2te(ex,ey,ep,D);
64
65		end
66
67		if (Coord(j,2) < Coord(j+1,2))
68		l = l + 1
69
70		x1 = Coord(j,1);
71		y1 = Coord(j,2);
72
73		x2 = Coord(j+1,1);
74		y2 = Coord(j+1,2);
75
76		x3 = Coord(j+nodejump+1,1);
77		y3 = Coord(j+nodejump+1,2);
78
79
80		Edof(:,l) = [l j j+1 j+nodejump+1 ];
81
82
83		ex = [x1 x2 x3];
84		ey = [y1 y2 y3];
85
86
87		KE(:,:,l) = flw2te(ex,ey,ep,D);
88
89		end
90
91
92		end
93
94
95		l = 89;
96		for j=51:102,
97		i = Coord(j,1);
98
99		if (i == 200/3)
100		nodejump = 7;
101		elseif (i == 220/3)
102		nodejump = 6;
103		else
104		nodejump = 8;
105		end
106
107
108		if (Coord(j,2) < Coord(j+1,2))
109		l = l + 1
110
111		x1 = Coord(j,1);
112		y1 = Coord(j,2);
113
114		x2 = Coord(j+nodejump,1);
115		y2 = Coord(j+nodejump,2);
116
117		x3 = Coord(j+1,1);
118		y3 = Coord(j+1,2);
119
120
121		Edof(:,l) = [l j j+nodejump j+1 ];
122
123
124		ex = [x1 x2 x3];
125		ey = [y1 y2 y3];
126
127
128		KE(:,:,l) = flw2te(ex,ey,ep,D);
129
130		end
131
132		if ((j < 102) && Coord(j,2) < Coord(j+nodejump+1,2))
133		l = l + 1
134
135		x1 = Coord(j+nodejump+1,1);
136		y1 = Coord(j+nodejump+1,2);
137
138		x2 = Coord(j+1,1);
139		y2 = Coord(j+1,2);
140
141		x3 = Coord(j+nodejump,1);
142		y3 = Coord(j+nodejump,2);
143
144
145		Edof(:,l) = [l j+nodejump+1 j+1 j+nodejump ];
146
147
148		ex = [x1 x2 x3];
149		ey = [y1 y2 y3];
150
151
152		KE(:,:,l) = flw2te(ex,ey,ep,D);
153
154		end
155
156
157		end
158		Edof = Edof'
159		[Ex,Ey] = coordxtr(Edof,Coord,Dof,3)