Symulacje komputerowe zjawisk fizycznych z zakresu mechaniki

More documents

Recommendations

Info

do poziomu wystarczy sporządzić jej wykres (patrz plik sd_wyniki_wykres.py). Rysunek 24: Zależność drogi hamowania od kąta nachylenia jezdni do poziomu. Przecięcie zielonych linii obrazuje rozwiązanie dla wprowadzonych przez użytkownika parametrów początkowych (patrz rys. poniżej) Rysunek 25: Wyznaczenie drogi hamowania dla bryły sztywnej. 33
Ruch postępowo obrotowy na równi Jeżeli na równi pochyłej ustawimy walec [3],[8],[9] i puścimy swobodnie, zacznie się obracać i zjeżdżać. Będzie się tak działo tylko wtedy gdy współczynnik tarcia będzie miał odpowiednią wartość. Jeżeli nie byłoby tarcia lub byłoby ono zbyt małe, walec nie staczał by się, lecz zsuwał po równi a jego prędkość kątowa u podnóża równi mogłaby być nawet równa zeru. Przyjmijmy lokalny układ współrzędnych w którym oś x-ów jest równoległa po płaszczyzny równi a oś y-ów jest do niej prostopadła. Wówczas równania ruchu dla walca znajdującego się na równi pochyłej przyjmują następującą postać. ∑ F y local =mg cos−mg cos=m a y ∑ F xlocal =mg sin −T =m a x gdzie T - wartość siły tarcia a x - przyspieszenie liniowe środka ciężkości walca Gdy siła sprężystości równoważy siłę nacisku m a y =0 . Zakładamy także , że walec stacza się bez poślizgu. Ruch obrotowy dookoła środka masy opisany jest równaniem = I śr.m. Obrotu walca nie powoduje jednak ani siła sprężystości F s ani F c , gdyż obie przechodzą przez środek masy i ich moment wypadkowy jest równy zeru. Moment siły różny od zera daje tylko siła tarcia T , której ramię jest równe promieniowi walca. Możemy więc zatem napisać: F R=I śr.m. moment bezwładności walca jest równy I śr.m. = 1 2 m R2 , Rysunek 26: Rozkład sił działających na walec staczający się po równi pochyłej. a przyspieszenie kątowe = a R więc F =I śr.m. R = m a 2 Przyspieszenie walca na równi obliczymy podstawiając powyższe równanie do drugiego równania ruchu (dla osi x-ów). W rezultacie otrzymamy a= 2 3 g sin . Prędkość środka masy przy założeniu zerowej prędkości początkowej jest określona równaniem v= 4 3 g h gdzie h jest wysokością równi. 34
Page 2 and 3: Temat: Symulacje komputerowe zjawis
Page 4 and 5: zeczywiście dotyczącą interesuj
Page 6 and 7: Obiekty, nazwy i atrybuty. Obiekty
Page 8 and 9: p=(0.0,0.0,0.0) k=0.01 strzalkax =
Page 10 and 11: krzywa2.append(pos=(2,1,0), color=c
Page 12 and 13: ziemia=sphere(pos=vector(1.5,0.0,0.
Page 14 and 15: i maksimum dla każdej z osi. bbs=0
Page 16 and 17: zaokrągleń niewiele wpływają na
Page 18 and 19: efekt Coriolisa oraz przyjąć sta
Page 20 and 21: Rysunek 12: Symulacja ruchu punktu
Page 22 and 23: Tabela. Wyniki symulacji dla ruchu
Page 24 and 25: dla składowych y v y2 = 1 C d { s
Page 26 and 27: 4. Modelowanie obiektów (część
Page 28 and 29: Jeżeli nie określimy atrybutu dla
Page 30 and 31: Rysunek 20: Scena renderowana z uż
Page 32 and 33: Wykresy - znakowanie danych Przy ry
Page 36 and 37: Rysunek 27: Symulacja ruchu walca n
Page 38 and 39: 6. Modelowanie obiektów (część
Page 40 and 41: W połączeniu między punktami dzi
Page 42 and 43: Rysunek 33: Histogram ocen uzyskany
Page 44 and 45: 7. Wiele ciał w przestrzeni Ruchy
Page 46 and 47: Wyniki symulacji można porównać
Page 48 and 49: Prawo pól Keplera Ruch w czasie kt
Page 50 and 51: Rysunek 39: Symulacja trajektorii r
Page 52 and 53: vel = V_posr_inalgorithm vel_posred
Page 54 and 55: a prędkość pośrednia dla chwili
Page 56 and 57: Kolejna symulacja ilustruje trudnie
Page 58 and 59: Wahadło torsyjne. Wahadłem takim
Page 60 and 61: pędu jak i zmianę pędu równą i
Page 62 and 63: umieszczamy aktualne położenia i
Page 64 and 65: V rn =v 1init ⋅n #predkosc wzgled
Page 66 and 67: Inne sytuacje jakie mogą być symu
Page 68 and 69: • Przesuwanie obiektów Aby móc
Page 70 and 71: all_list.append(ball) • Konstrukc
Page 72 and 73: exchange=vj-vi #exchange momentum b
Page 74 and 75: kroków czasowych przy zachowaniu d
Page 76 and 77: Możemy zatem się spodziewać, że
Page 78 and 79: # Plot the velocity histogram v_lis
Page 80 and 81: Kinetyczny model gazu doskonałego.
Page 82 and 83: pos = array(poslist) p = array(plis
Page 84 and 85:
hitlist = sort(nonzero(hit.flat)).t
Page 86 and 87:
v CM i M v śr M M =v i M m i m i v
Page 88 and 89:
Rysunek 54: Rozkłady prędkości M
Page 90 and 91:
9. Zakończenie Eksperyment kompute
show all

Symulacje komputerowe zjawisk fizycznych z zakresu mechaniki

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?