590 likes | 813 Views
ZADATAK 1. Izvesti matricu elementa i vektor desne strane koristeći metodu funkcije zaokreta za prizmatični štap pravokutnog poprečnog presjeka prema crtežu. Koristi se osamčvorni konačni element. REALNE KOORDINATE ČVOROVA:. DETERMINANTA JAKOBIJANE. ČLANOVI MATRICE KRUTOSTI.
E N D
ZADATAK 1 Izvesti matricu elementa i vektor desne strane koristeći metodu funkcije zaokreta za prizmatični štap pravokutnog poprečnog presjeka prema crtežu. Koristi se osamčvorni konačni element. REALNE KOORDINATE ČVOROVA:
MATRIČNA JEDNADŽBA RAVNOTEŽE PRAVOKUTNOG KONAČNOG ELEMENTA matrica konačnog elementa vektor rješenja vektor opterećenja (vektor aktivnog opterećenja i reakcija u točkama sa zadanim kinematičkim rubnim uvjetom).
ZADATAK 2 Za štap kvadratnog poprečnog presjeka stranice duljine 2a opterećenog momentom torzije Mt (tako da je zasukan za kut ) potrebno je odrediti funkciju vitoperenja (x, y), komponente posmičnih naprezanja te torzijsku krutost presjeka. Materijal je izotropan s modulom posmika G. - nepoznati samo pomaci čvorova 4, 5 i 6; ostali čvorovi nemaju pomaka
FUNKCIJA VITOPERENJA NA K.E. PRIBLIŽNO RJEŠENJE
5 6 4 7 ξ 8 3 2 1 ZADATAK 3 Varijacijska metoda minimuma potencijalne energije koristeći osi simetrije i izoparametrički osmočvorni konačni element.
5 6 4 7 ξ 8 3 2 1 Rješenje tražimo u obliku linearne kombinacije baznih funkcija: Bazne funkcije na osmočvornom izoparametričkom konačnom elementu u punom obliku su: NE ZADOVOLJAVAJU RUBNE UVJETE SPREGNEMO U JEDNU BAZNU FUNKCIJU
Naše bazne funkcije: Bazna funkcija Nina rupi =1 =2Ak Član koji se dodaje bilo kojem čvoru (ali samo jednom) na rubu rupe. U našem slučaju imamo samo jednu baznu funkciju pa ćemo ovaj član dodati pripadajućem i-tomindeksu koji ima i ta bazna funkcija. Negativno zbog smjera integracije!
Moment torzije: Posmična naprezanja:
Štap poprečnog presjeka u obliku istostraničnog trokuta opterećen je momentom torzije . Koristeći metodu naprezanja, potrebno je odrediti polje naprezanja. y x ZADATAK 4 SAINT – VENANTOV PROBLEM Primjer: Stranice:
y 4 +1 -1 x 1 3 ξ +1 -1 2 Geometrija područja:
y x . . 7 6 8 2 1 5 1 2 2 4 ξ 3 Rješenje će se odrediti pomoću istih baznih funkcija kao i za kvadratno područje.
Rješenje tražimo u obliku linearne kombinacije baznih funkcija:
Pogreška iznosi Za Točno rješenje je:
Pogreška iznosi . Usporedba točnosti momenta torzije dobivenog analitičkim i numeričkim postupkom: Točno rješenje: Granice područja: po “x” po “y” Za: Numerički:
Pogreška iznosi NAPREZANJE: Za Točno: Za:
y x 1. ZADATAK: Zadan je prizmatični štap opterećen momentom torzije Mt (zasukan za kut ). Poprečni presjek štapa prikazan je na crtežu. Modul posmika iznosi . Koristeći jedan 8.-čvorni izoparametarski konačni element treba odrediti funkcije rješenja i , komponente naprezanja i torzijsku krutost, te međusobno usporediti dobivene rezultate.
1.1.METODA FUNKCIJE ZAOKRETA Funkcional potencijalne energije u ovom slučaju glasi: Funkcija se može izraziti kao linearna kombinacija: VEKTOR POMAKA ČVOROVA NA KRAJEVIMA KE VEKTOR BAZNIH FUNKCIJA DEFINIRANIH NA POMOĆNOM KOORDINATNOM SUSTAVU
0 0 0 0 0 0 KAKO NAM IZGLEDAJU NAŠE BAZNE FUNKCIJE? 0 0 0 0 0 0 0 0
Deriviranjem izraza za funkciju po koordinatnim smjerovima x,y, uvrštenjem u funkcional te njegovim minimiziranjem dobije se sustav jednadžbi: j=1,2...,8 odnosno matrično zapisan: =? =? KAKO BI VI ODREDILI ? =? =? SLIJEDE IZ J-1
MATRICA KRUTOSTI U REALNIM KOORDINATAMA PRELAZI U
y 8 4 x 3 1 2 7 6 5 EKVIVALENTNE ČVORNE SILE: ...
IDEMO SADA RIJEŠITI NAŠ ZADATAK! y 8 + 4 x 3 1 2 7 6 5 PRVO ĆEMO DEFINIRATI ČVOROVE KE!
y 8 + (-1,1) (0,1) (1,1) 4 (-1,0) (1,0) (0,0) x 3 1 2 7 6 5 (-1,-1) (1,-1) (0,-1)
y + x 1
REALNE KOORDINATE ČVOROVA: y 8 8 1 1 2 2 4 4 x 3 3 3 1 3 1 2 2 7 6 7 5 6 5 4 4 Koordinate točaka (x,y) izražene preko koordinata (,) su: 5 5 6 6 7 7 8 8
Jacobijan je : a njegova determinanta Sada se mogu izračunati potrebne derivacije:
(VEKTOR AKTIVNOG OPTEREĆENJA I REAKCIJA U TOČKAMA SA ZADANIM KINEMATIČKIM RUBNIM UVJETOM). POZNATI POMACI RUBNI UVJETI: y 8 8 Koji članovi matrice množe poznate pomake? Koji članovi matrice ‘ostaju’? Koji članovi matrice koji odgovaraju poznatim ‘reakcijama’? 4 4 x 3 1 3 1 2 2 7 6 7 5 6 5
Suma svih reakcija je jednaka nuli što znači da je u smjeru osi štapa uspostavljena ravnoteža.
Funkcija vitoperenja ima oblik: Torzijska krutost poprečnog presjeka računa se po izrazu: Za cijeli poprečni presjek se dobije: