PŁYTA IZOTROPOWA SZEŚCIOKĄTNA
Dla danej płyty stropowej, obciążonej na całej powierzchni siłami prostopadłymi
do powierzchni płyty o wartości 1322N/m2, obliczyć siły przekrojowe i przemieszczenia.
Należy wykorzystać symetrię konstrukcji.
Rys.1. a) Model płyty; b) ćwiartka płyty z zaznaczonymi węzłami
bazowymi do generatora gwe.
ROZPOCZĘCIE PRACY
$ feas - wejście do systemu feas,
FEAS> - system zgłasza gotowość do pracy,
OPIS KONSTRUKCJI
FEAS> ok - wejście do podsystemu OK - Opis Konstrukcji,
Podaj rodzaj konstrukcji pl-i - pl-i - symbol płyty izotropowej,
OK> - podsystem zgłasza gotowość do pracy,
wczytanie współrzędnych wezłów
OK> w a-1 10.2 0 - węzeł o nazwie a-1 ma współrzędne (10.2,0),
OK> w a-2 21.6 0 - węzeł o nazwie a-2 ma współrzędne (21.6,0),
OK> w a-3 10.8 18.72 - węzeł o nazwie a-3 ma współrzędne (10.8,18.72),
OK> w a-4 5.1 8.84 - węzeł o nazwie a-4 ma współrzędne (5.1,8.84),
OK> w b-1 0 8.84 - węzeł o nazwie b-1 ma współrzędne (0,8.84),
OK> w b-2 5.1 8.84 - węzeł o nazwie b-2 ma współrzędne (5.1,8.84),
OK> w b-3 10.8 18.72 - węzeł o nazwie b-3 ma współrzędne (10.8,18.72),
OK> w b-4 0 18.72 - węzeł o nazwie b-4 ma współrzędne (0,18.72),
generowanie węzłów i elementów
OK> gwe w-1 e-1 a-1 do a-4 10 18 - czworokąt o węzłach od
a-1 do a-4 zostanie podzielony na elementy trójkątne, odcinek między węzłami a-1 i a-2
będzie podzielony na 10 części, odcinek między a-2 a a-3 na 18,
liczba wygenerowanych : węzłów =209, elementów=360
- system informuje, że wygenerowanych zostało 209 węzłów o nazwach od w-1
do w-209, oraz 360 elementów o nazwach od e-1 do e-360,
OK> u a-1 do a-4 - usuwanie pomocniczych węzłów a-1, a-2, a-3, a-4,
USUWANIE WEZLA: A-1 [T/N/W/K]w
OK> gwe w-210 e-361 b-1 do b-4 9 10 - czworokąt o węzłach od
b-1 do b-4 zostanie podzielony na elementy trójkątne, odcinek między węzłami b-1 i b-2
będzie podzielony na 9 części, odcinek między b-2 a b-3 na 10,
liczba wygenerowanych : węzłów =110, elementów=180 - system informuje, że wygenerowanych zostało 110 węzłów o nazwach od w-210
do w-319, oraz 180 elementów o nazwach od e-361 do e-540,
OK> u b-1 do b-4 - usuwanie pomocniczych węzłów b-1, b-2, b-3, b-4,
USUWANIE WEZLA: B-1 [T/N/W/K] w OK> zw w-1 do w-319 - zszycie dwóch wygenerowanych siatek,
węzły o tych samych współrzędnych będą eliminowane,
Rys.2. Siatka poziału na węzły i elementy (ćwiartka konstrukcji).
wczytanie danych materiałowych
OK> m m1 - wczytanie danych materiałowych o nazwie m1,
MAT: Ro=0 ?> 0
MAT: E=0 ?> 36e6 - wpisujemy wartości odpowiednich stałych,
MAT: Ni=0 ?> 0.16667
MAT: AlfaT=0 ?> 0
geometria przekroju
OK> g g1 - wczytanie danych geometrycznych o nazwie g1,
GEOM:H=0 ?> 0.3 - wpisujemy grubość płyty,
elementy
OK> pe e-1 do e-540 m1 g1 - przyporządkowanie wszystkim elementom
danych materiałowych -- m1 i danych geometrycznych -- g1,
warunki brzegowe (w układzie globalnym)
OK> wb w-1 do w-199 co 11 uz fix fiy - w co 11 węźle od w-1
do w-199 odbieramy trzy stopnie swobody uz, fix, fiy,
OK> wb w-210 do w-218 uz fix fiy - węzly
od w-210 do w-218 nie mają możliwości przesuwu uz i obrotu fix i fiy}
OK> wb w-2 do w-10 fix - węzły od w-2 do w-10 mają odebraną
możliwość obrotu fix,
OK> wb w-220 do w-300 co 10 fiy - co 10--ty węzeł
pomiędzy w-220 a w-300 nie ma
możliwość obrotu fiy,
OK> wb w-310 uz fiy
OK> wb w-11 uz fix
OK> wb w-313 w-316 w-44 w-77 w-110 w-143 w-176 w-209 uz
obciążenie elementowe ciągłe na powierzchni elementu
OK> oe oe-1 /ln - wczytanie obiążenia elementowego o nazwie
oe-1,
OE:pz1=0 ?> -13.22
OE:mx1=0 ?> 0
OE:my1=0 ?> 0
OE:pz2=0 ?> -13.22
OE:mx2=0 ?> 0
OE:my2=0 ?> 0 - wczytanie wartości danego obciążenia,
OE:pz3=0 ?> -13.22
OE:mx3=0 ?> 0
OE:my3=0 ?> 0
OK> pe e-1 do e-540 oe-1 - przyporządkowanie elementom od e-1 do e-540 obciążenia
ciągłego o nazwie oe1,
WYJŚCIE Z PODSYSTEMU OK - ROZWIĄZANIE ZADANIA
OK> zap zad4 - zapisanie danych do zadania do pliku zad4,
TYTUL> Płyta izotropowa
OK> .. - wyjście z podsystemu OK,
FEAS> ro s - rozwiązanie zadania (statyka),
WYŚWIETLENIE WYNIKÓW
FEAS> WS - wejście do podsystemu ws,
WS> pm - wyświetlenie przemieszczeń,
WS> ne - wyświetlenie sił przekrojowych w elementach,
WS> r - wyświetlenie reakcji podporowych,
WS> dp - definiowanie pliku wynikowego
podaj nazwe pliku: a - plik wynikowy będzie miał nazwę A.PRN,
WS> pm /p - zapisanie przemieszczeń do pliku a, analogicznie
możemy zapisać przemieszczenia, siły przekrojowe (ne), reakcje (r)...,
WS> .. - wyjście z podsystemu,
GRAFICZNA PREZENTACJA WYNIKÓW
FEAS> GRAF - wejście do podsystemu GRAF,
GRAF> si - rysowanie siatki elementów,
GRAF> si /kn - rysowanie konturu siatki elementów,
GRAF> war pm uz - rysowanie warstwic przemieszczenia uz ,
GRAF> war nw mxx - rysowanie warstwic momentu mx,
GRAF> i - rysowanie w izometrii,
GRAF> pm - rysowanie siatki odkształconej w izometrii.
ZAKOŃCZENIE PRACY
FEAS> ..
Czy chcesz zakonczyc prace (t/n) t
$
Rys.3. Otrzymane wyniki, przemieszczenia, warstwice momentu Mx.
Dla danej płyty stropowej, obciążonej na całej powierzchni siłami prostopadłymi do powierzchni płyty o wartości 1322N/m2, obliczyć siły przekrojowe i przemieszczenia. Należy wykorzystać symetrię konstrukcji.
Rys.1. a) Model płyty; b) ćwiartka płyty z zaznaczonymi węzłami
bazowymi do generatora gwe.
Rys.2. Siatka poziału na węzły i elementy (ćwiartka konstrukcji).
Rys.3. Otrzymane wyniki, przemieszczenia, warstwice momentu Mx.
$ feas - wejście do systemu feas,
FEAS> - system zgłasza gotowość do pracy,
FEAS> ok - wejście do podsystemu OK - Opis Konstrukcji,
Podaj rodzaj konstrukcji pl-i - pl-i - symbol płyty izotropowej,
OK> - podsystem zgłasza gotowość do pracy,
OK> w a-1 10.2 0 - węzeł o nazwie a-1 ma współrzędne (10.2,0),
OK> w a-2 21.6 0 - węzeł o nazwie a-2 ma współrzędne (21.6,0),
OK> w a-3 10.8 18.72 - węzeł o nazwie a-3 ma współrzędne (10.8,18.72),
OK> w a-4 5.1 8.84 - węzeł o nazwie a-4 ma współrzędne (5.1,8.84),
OK> w b-1 0 8.84 - węzeł o nazwie b-1 ma współrzędne (0,8.84),
OK> w b-2 5.1 8.84 - węzeł o nazwie b-2 ma współrzędne (5.1,8.84),
OK> w b-3 10.8 18.72 - węzeł o nazwie b-3 ma współrzędne (10.8,18.72),
OK> w b-4 0 18.72 - węzeł o nazwie b-4 ma współrzędne (0,18.72),
OK> gwe w-1 e-1 a-1 do a-4 10 18 - czworokąt o węzłach od
a-1 do a-4 zostanie podzielony na elementy trójkątne, odcinek między węzłami a-1 i a-2
będzie podzielony na 10 części, odcinek między a-2 a a-3 na 18,
liczba wygenerowanych : węzłów =209, elementów=360
- system informuje, że wygenerowanych zostało 209 węzłów o nazwach od w-1
do w-209, oraz 360 elementów o nazwach od e-1 do e-360,
OK> u a-1 do a-4 - usuwanie pomocniczych węzłów a-1, a-2, a-3, a-4,
USUWANIE WEZLA: A-1 [T/N/W/K]w
OK> gwe w-210 e-361 b-1 do b-4 9 10 - czworokąt o węzłach od
b-1 do b-4 zostanie podzielony na elementy trójkątne, odcinek między węzłami b-1 i b-2
będzie podzielony na 9 części, odcinek między b-2 a b-3 na 10,
liczba wygenerowanych : węzłów =110, elementów=180 - system informuje, że wygenerowanych zostało 110 węzłów o nazwach od w-210
do w-319, oraz 180 elementów o nazwach od e-361 do e-540,
OK> u b-1 do b-4 - usuwanie pomocniczych węzłów b-1, b-2, b-3, b-4,
USUWANIE WEZLA: B-1 [T/N/W/K] w OK> zw w-1 do w-319 - zszycie dwóch wygenerowanych siatek,
węzły o tych samych współrzędnych będą eliminowane,
OK> m m1 - wczytanie danych materiałowych o nazwie m1,
MAT: Ro=0 ?> 0
MAT: E=0 ?> 36e6 - wpisujemy wartości odpowiednich stałych,
MAT: Ni=0 ?> 0.16667
MAT: AlfaT=0 ?> 0
OK> g g1 - wczytanie danych geometrycznych o nazwie g1,
GEOM:H=0 ?> 0.3 - wpisujemy grubość płyty,
OK> pe e-1 do e-540 m1 g1 - przyporządkowanie wszystkim elementom
danych materiałowych -- m1 i danych geometrycznych -- g1,
OK> wb w-1 do w-199 co 11 uz fix fiy - w co 11 węźle od w-1
do w-199 odbieramy trzy stopnie swobody uz, fix, fiy,
OK> wb w-210 do w-218 uz fix fiy - węzly
od w-210 do w-218 nie mają możliwości przesuwu uz i obrotu fix i fiy}
OK> wb w-2 do w-10 fix - węzły od w-2 do w-10 mają odebraną
możliwość obrotu fix,
OK> wb w-220 do w-300 co 10 fiy - co 10--ty węzeł
pomiędzy w-220 a w-300 nie ma
możliwość obrotu fiy,
OK> wb w-310 uz fiy
OK> wb w-11 uz fix
OK> wb w-313 w-316 w-44 w-77 w-110 w-143 w-176 w-209 uz
OK> oe oe-1 /ln - wczytanie obiążenia elementowego o nazwie
oe-1,
OE:pz1=0 ?> -13.22
OE:mx1=0 ?> 0
OE:my1=0 ?> 0
OE:pz2=0 ?> -13.22
OE:mx2=0 ?> 0
OE:my2=0 ?> 0 - wczytanie wartości danego obciążenia,
OE:pz3=0 ?> -13.22
OE:mx3=0 ?> 0
OE:my3=0 ?> 0
OK> pe e-1 do e-540 oe-1 - przyporządkowanie elementom od e-1 do e-540 obciążenia
ciągłego o nazwie oe1,
OK> zap zad4 - zapisanie danych do zadania do pliku zad4,
TYTUL> Płyta izotropowa
OK> .. - wyjście z podsystemu OK,
FEAS> ro s - rozwiązanie zadania (statyka),
FEAS> WS - wejście do podsystemu ws,
WS> pm - wyświetlenie przemieszczeń,
WS> ne - wyświetlenie sił przekrojowych w elementach,
WS> r - wyświetlenie reakcji podporowych,
WS> dp - definiowanie pliku wynikowego
podaj nazwe pliku: a - plik wynikowy będzie miał nazwę A.PRN,
WS> pm /p - zapisanie przemieszczeń do pliku a, analogicznie
możemy zapisać przemieszczenia, siły przekrojowe (ne), reakcje (r)...,
WS> .. - wyjście z podsystemu,
FEAS> GRAF - wejście do podsystemu GRAF,
GRAF> si - rysowanie siatki elementów,
GRAF> si /kn - rysowanie konturu siatki elementów,
GRAF> war pm uz - rysowanie warstwic przemieszczenia uz ,
GRAF> war nw mxx - rysowanie warstwic momentu mx,
GRAF> i - rysowanie w izometrii,
GRAF> pm - rysowanie siatki odkształconej w izometrii.
FEAS> ..
Czy chcesz zakonczyc prace (t/n) t
$