Re: projectielbeweging met kwadratische wrijvingskracht
Geplaatst: di 24 sep 2019, 12:05
Vermoedelijk dat door k zo te definiëren de massa m van het projectiel niet ter zake doet.
Ik heb twijfels over jouw Python code.
In regels 21 en 22 staat bijvoorbeeld:
vx = v0*np.sin(th0)
vy = v0*np.cos(th0)
Ik weet weliswaar niet wat np betekent maar het lijkt me dat hier sin en cos verwisseld zijn.
Ik neem hierbij aan dat de x-as de horizontale as is en de y-as de verticale, en dat θ de hoek is tussen de bewegingsrichting en de horizontale x-as.
Zolang theta 45 graden is maakt deze fout toevallig niet uit maar bij andere hoeken wel.
Ook schreef je eerder: Fx= - k.vt2.sin(θ) waarbij sin en cos ook verwisseld zijn.
Een ander probleem is de tekenconventie voor Fd, v en a waar ukster en jij je, voor zover ik kan zien, niet om lijken te bekommeren.
Nadat het projectiel het hoogste punt bereikt dienen de tekens (+/-) van diverse variabelen te veranderen, echter door het kwadrateren van een negatieve v of door twee negatieve v's met elkaar te vermenigvuldigen ontstaat een positieve Fd terwijl die na de top negatief moet zijn want dan is Fd omhoog gericht.
De bewegingsvergelijking, in het allereerste bericht, voor de y-richting heeft dat probleem al want vt en vy hebben altijd hetzelfde teken. Na de top, als beiden negatief zijn, levert dat toch een positieve Fd op, terwijl die nu negatief zou moeten zijn. En dat heeft vanzelfsprekend invloed op versnelling a.
Bij gebruik van Pythagoras (in de Python code) loop je vanwege het kwadrateren ook tegen dat probleem aan.
ik heb zo als je kunt zien in het opzetten van de vergelijkingen in x- en y richting wel degelijk rekening gehouden met de gangbare tekenafspraken.