File Edit Ted Go Celi Tools Debug Desktop Window Help
|
I * % m *> o |
& ® - |
H |
Stack: |
k
|
♦tu r[3 E) Us |
ł-1 o |
+ |
,|,i |
X |
Ct |
|
35 |
% masa |
kg; |
m = 100;
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60 61 62
63
64
65
66
67
68
69
70
% wezly z masami [numer wezla z masa wm = [2 2*m;
4 m]
masa]
%% Black box — MES
% Zdefiniowanie przekątniowej macierzy mas M = zeros(max(max(Dof)));
E;'. for n = l:size(wm,1)
M ( Dof (wm (n, 1) ,1) M( Dof(wm(n,1) ,2)
L end
Dof (wm (n, 1) ,1) ) Dof(wm(n,1) ,2) )
wm (n, 2) ; wm (n, 2) ;
%----- Rozwiniecie macierzy Coord w macierze ex oraz ey
[ex,ey] = coordxtr(Edof, Coord, Dof, 2);
%-
A
I
E
— Charakterystyki geometryczne i materiałowe ----
7.64e-4; % pole przekroju IPE80 [m2]
80.1e-8; % moment bezwładności IPE80 [m4]
205e9; % moduł Young'a [Pa] dla stali
ep = [E A I];
% GL03ALNA MACIERZ SZTYWNOŚCI K = zeros(max(max(Dof)));
[nel,kol]=size(Edof) ;
0 for i=l:neł
[Ke] = beam2e(ex(i,:),ey(i,:),ep); K = assem(Edof(i,:),K,Ke);
L end
Ln 1
Col 1
□