Uwaga: np. IM XYi - jest komendą w systemie FEAS- KAM, dX dY - są parametrami
I) OPIS GEOMETRII
6. Inicjacja wektorów geometrii prętów - opisujących ich długość oraz orientację w stosunku do globalnego układu współrzędnych (dX, dY - przyrosty współrzędnych: koniec pręta- początek)
IM XYi 2 1
dX dY
7. Inicjacja wektorów własności materiałowych : E, Ro i alfaT - odpowiednio moduł Younga , gęstość właściwa, współczynnik rozszerzalności termicznej materiału
IM MATi 3 1
E Ro alfaT
8. Inicjacja wektorów własności geometrycznych : A, Jz - odpowiednio pole przekroju, moment bezwładności przekroju
IM GEOi 2 1
A Jz
9. Inicjacja wektorów alokacji Ai - opisujących rozmieszczenie niewiadomych w kolejnych prętach: n1, n2 n3 - nr niewiadomych na początku, n4n5n6- na końcu pręta.
IMI Ai 6
n1 n2 n3 n4 n5 n6
10. Budowa macierzy sztywności prętów Ki na podstawie biblioteki R2 - XYi, MATi, GEOi - poprzednio zdefiniowane macierze.
R2-ST Ki XYi MATi GEOi
11. Zdefiniowanie i wyzerowanie globalnej macierzy sztywności konstrukcji K - n - jest liczbą niewiadomych metody
DM K n n
12. Wypełnienie globalnej macierzy sztywności K - na podstawie macierzy sztywności prętów Ki i wektorów alokacji Ai
DAL K Ki Ai
13. Inicjacja wektorów pomocniczych DSi - określających obciążenie elementowe (międzywęzłowe) prętów
- dla sił skupionych - wektor DSi, w którym Sx, Sy - składowe lokalne siły skupionej , Mz - moment skupiony, ksi - bezwymiarowa odcięta położenia obciążenia
IM DSi 4 1
Sx Sy Mz ksi
- dla obc. rozłożonego - wektor DLi, w którym qx1, qy1, mz1 oraz qx2, qy2, mz2 rzędne składowych obciążenia rozłożonego w ukl. lokalnym na początku i końcu, ksi1, ksi2 - bezwymiarowe odcięte wyznaczające początek i koniec obciążenia
IM DLi 481
qx1 qy1 mz1 qx2 qy2 mz2 ksi1 ks2
14. Budowa wektorów sił przywęzłowych QLi, QSi dla poszczególnych prętów- sprowadzających obciążenie elementowe do końców prętów (na podstawie biblioteki R2)
R2-SS QSi XYi DSi - dla obciążeń skupionych elementowych
R2-LN QLi XYi DLi - dla obciążeń rozłożonych elementowych
15. Zdefiniowanie i wyzerowanie sumarycznego wektora sił przywęzłowych konstrukcji od obciążeń elementowych QEO - n - liczba niewiadomych
DM QEO n 1
16. Wypełnienie sumarycznego wektora sił przywęzłowych konstrukcji od obciążeń elementowych QEO - na podstawie wektorów obciążeń poszczególnych prętów QSi lub QLi i wektorów alokacji Ai
DWA QEO QSi Ai
DWA QEO QLi Ai
17. Inicjacja wektora obciążeń węzłowych konstrukcji QW0 - pochodzących od obciążeń konstrukcji działających na jej węzły; w1, w2, w3, ... wn - składowe obciążeń węzłowych w układzie globalnym w kolejności odpowiadanej kolejnym niewiadomym metody lub zero
IM QWO n 1
w1 w2 w3 ..wi ... wn
18. Zdefiniowanie sumarycznego wektora wyrazów wolnych QO
DM QO n 1
19. Obliczenie sumarycznego wektora wyrazów wolnych Q0 układu równań metody przemieszczeń jako sumy wektora sił przywęzłowych QEO i obciążeń węzłowych QWO.
D QEO QWO QO
20 Kopiowanie utworzonej macierzy sztywności konstrukcji K i wektora wyrazów wolnych QO:
KP QO Z
KP K KZ
21. Rozwiązanie układu równań kanonicznych metody przemieszczeń KZ* X = Z Rozwiązanie (niewiadome geometryczne) umieszczone jest przez program w wektorze wyrazów wolnych Z.
RO KZ Z
22. Wyznaczenie wektorów przemieszczeń przywęzłowych Qi dla poszczególnych prętów - na podstawie globalnego wektora przemieszczeń układu Z wyznaczonego powyżej oraz wektora alokacji.
TWA Z Ai Qi
(23. Pominięto)
24. Wyznaczenie wektorów zawierających rzędne wykresów sił przekrojowych na końcach oraz w `k' przekrojach pośrednich - od przemieszczeń węzłów pręta NQi , obciążeń międzywęzłowych.NSi, NLi, NMi
R2-NP NQi XYi MATi GEOi Qi k - od przemieszczeń końców pręta
R2-NS NSi XYi DSi k - od sił skupionych międzywęzłowych pręta
R2-NS NMi XYi DMi k - od momentów skupionych międzywęzłowych pręta
R2-NL NLi XYi DLi k - od obciążeń rozłożonych międzywęzłowych pręta
25. Wyznaczenie sumarycznych macierzy sił przekrojowych d dla poszczególnych prętów jako suma wektorów wyznaczonych w punkcie poprzednim, np
D NQi NSi NTMi - dla pręta obciążonego siła skupioną
D NQi NLi NTMi - dla pręta obciążonego obciążeniem rozłożonym
KP NQi NTMi - dla pręta nieobciążonego
26. Wyświetlenie wyników - macierzy NTMi , zapisanie do pliku
WS NTMi
ZAP NTMi nazwa_pliku
Wyniki umieszczone zostaną w pliku nazwa_pliku.mat w katalogu FEAS/BIN
27. Sporządzenie wykresów sił wewnętrznych na podstawie otrzymanych wyników
Fragment pliku z danymi (niezwiązany z pod