Symulacje komputerowe zjawisk fizycznych z zakresu mechaniki

More documents

Recommendations

Info

p=(0.0,0.0,0.0) k=0.01 strzalkax = arrow(pos=p, axis=(1,0,0), shaftwidth=k, color=color.red) strzalkay = arrow(pos=p, axis=(0,1,0), shaftwidth=k, color=color.green) strzalkaz = arrow(pos=p, axis=(0,0,1), shaftwidth=k, color=color.blue) #----------------- while b
przykład moja_krzywa.pos[0] jest pozycją pierwszego punktu w moja_krzywa. Możemy również otrzymać krzywą bezpośrednio z listy współrzędnych [1]. Listę taką tworzymy podobnie jak krotkę (sekwencję) zamykając listę współrzędnych w okrągłych nawiasach. I tak na przykład aby narysować kwadrat wystarczy napisać: kwadrat = curve(pos=[(0,0),(0,1),(1,1),(1,0),(0,0)]) Oczywiście (1,1) jest skrótem (1,1,0). Jednakże, nie możemy mieszać ze sobą punktów 2D oraz 3D w jednej liście. Krzywa posiada grubość, definiowaną poprzez parametr radius. Zatem średnica krzywej jest równa podwójnemu promieniowi. krzywa=curve(pos=[(0,0,0), (1,0,0)], radius=0.05) Domyślny promień krzywej wynosi 0, wtedy jest rysowana bardzo cienka krzywa. Jeżeli podamy bardzo małą niezerową wartość, krzywa stanie się prawie niewidoczna dlatego zaleca się używać wartości zero ponieważ wtedy średnica krzywej jest domyślnie równa jednemu pikselowi. Krzywą narysować możemy również z wygenerowanej listy punktów, w tym celu możemy posłużyć się funkcją arange( [start,] stop[, step]) , która zwraca sekwencję liczb. Poniżej przedstawiony jest fragment skryptu pozwalającego narysować sprężynę. # Rysowanie sprężyny c = curve( x = arange(0,20,0.1) ) c.y = sin( 2.0*c.x ) c.z = cos( 2.0*c.x ) Parametryczne rysowanie także nie powinno sprawiać problemu. Jak wspomniane było wczśniej niektóre atrybuty, na przykład color mogą być różne dla każdego punktu w krzywej. Kiedy będziemy chcieli nagle zmienić kolor krzywej, po prostu należy dodać kolejny punkt w miejscu zmiany koloru, lecz o innym kolorze. from visual import * c = curve( pos=[(0,0,0), (1,0,0)], color=color.red ) c.append( pos=(1,1,0) ) # dodanie punktu w kolorze czerwonym c.append( pos=(1,1,0), color=color.blue) # ten sam punkt, lecz kolor niebieski c.append( pos=(0,1,0) ) # dodanie punktu w kolorze niebieskim Poniżej przedstawiony jest pełny skrypt obrazujący niektóre z wymienionych operacji na krzywych. Pierwsza krzywa (krzywa1) jest rysowana z punktów które mogą pochodzić na przykład z obliczeń, kolejna krzywa (krzywa2) zawiera jedynie dwa punkty w kolorze szarym. Do tak stworzonej krzywej dodajemy punkty przy pomocy funkcji append(). Ostatnia krzywa (krzywa3) jest rysowana parametrycznie. Do barwienia krzywej użyto pośrednio schematu HSV (Hue, Saturation, Value – barwa, nasycenie, wartość), który następnie zamieniony został na schemat RGB. from visual import * scene.background=(1.0,1.0,1.0) a=(0,1,0) b=(1,1,0) c=(2,2,0) d=(3,2,0) krzywa1=curve(pos=[a,b,c,d],color=(0.5,0.5,0.5), radius=0.05) krzywa2=curve(pos=[(0,0,0), (1,0,0)],color=(0.5,0.5,0.5), radius=0.05) # dodanie punktow do krzywej 8
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 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 34 and 35: do poziomu wystarczy sporządzić j
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?