Symulacje komputerowe zjawisk fizycznych z zakresu mechaniki

krzywa2.append(pos=(2,1,0), color=color.green)
krzywa2.append(pos=(3,1,0), color=color.yellow)
# krzywa rysowana parametrycznie
n = 16 # ilosc pelnych cykli (zwojow) sprezyny na jednostke dlugosci
r = 2 # promien sprezyny :)
rd = 0.05 # promien drutu
krzywa3=curve(pos=[], radius=0.05)
for i in arange(0.,1.,0.001):
k=2*(1/(i+1)*i)
#print k
c=(k,1,1)
krzywa3.append(pos=(-4+10*i,-2+i*r*cos(n*2*math.pi*i),
i*r*sin(n*2*math.pi*i)), color=color.hsv_to_rgb(c))
#----- uklad trzech strzalek
strzalkax = arrow(pos=(-2,-2,-2), axis=(1,0,0), shaftwidth=1, color=color.red)
strzalkay = arrow(pos=(-2,-2,-2), axis=(0,1,0), shaftwidth=1, color=color.green)
strzalkaz = arrow(pos=(-2,-2,-2), axis=(0,0,1), shaftwidth=1, color=color.blue)
#-----
Jednym z wielu sposobów zastosowania
krzywych jest wizualizacja trajektorii
punktu materialnego lub specyficznych
aspektów ruchu bryły sztywnej. Kolejny
skrypt przedstawia animację kuli
poruszającej się w przestrzeni. Za każdym
razem gdy kula podąża w kierunku
rosnących wartości y do krzywej
(trajectory) jest przypisywany punkt o
współrzędnych równych położeniu kuli.
Zatem poruszająca się kula pozostawi za
sobą ślad w postaci krzywej. Dodatkowo na
płaszczyznach xy, yz, oraz zx są rysowane
rzuty tej trajektorii.
from visual import *
#scene.background=(1.0,1.0,1.0)
scene.range = 3.0
Rysunek 4: Wizualizacja krzywych.
rc=0 # promien krzywych
col=(0.1,0.7,0.7)
trajectory=curve(pos=[],color=col, radius=rc) #trajektoria kuli xyz
c1=curve(pos=[],color=(0.8,0.5,0.1), radius=rc) #rzut na plaszczyzne xy
c2=curve(pos=[],color=(0.2,0.6,0.2), radius=rc) #rzut na plaszczyzne yz
c3=curve(pos=[],color=(0.9,0.1,0.9), radius=rc) #rzut na plaszczyzne zx
a=1.0
b=1.0
c=0.0
kula=sphere(pos=vector(a,c,b), radius=0.05, color=col)
#p1,p2,p3 - prowadnice
szary=(0.7,0.7,0.7)
p1=curve(pos=[(kula.pos),(kula.pos[0],kula.pos[1],0.0)],color=szary, radius=rc)
p2=curve(pos=[(kula.pos),(kula.pos[0],0.0,kula.pos[2])],color=szary, radius=rc)
p3=curve(pos=[(kula.pos),(0.0,kula.pos[1],kula.pos[2])],color=szary, radius=rc)
r = 0.6 # promien okregu
#----- uklad trzech strzalek
p=(0.0,0.0,0.0)
k=0.01
strzalkax = arrow(pos=p, axis=(1,0,0), shaftwidth=k, color=color.red)
9