scenario A implementatie

This commit is contained in:
Stein_van_vliet 2024-02-28 15:53:46 +01:00
parent d59e813517
commit 85c8c1bd80
1 changed files with 50 additions and 8 deletions

View File

@ -9,23 +9,44 @@ m = 6 #[kg]
r = 0.075 #[m]
Fz = m*g/4
dt = 0.01
#SCENARIO 2a
F_aan = 0 #[N]
Fn = Fz
F_rol = Fn*mu_r
F_eff = 50
v_kar = 0 #[m/s]
F_eff = 0
v_kar = 2.1 #[m/s]
s_kar = 0 #[m]
time = 20
dt = 0.01
F_list = []
v_list = []
s_list = []
t_list = []
"""
#SCENARIO 1a
Fn = Fz*mp.cos(mp.radians(20))
Fzt = Fz*mp.sin(mp.radians(20))
F_rol = Fn*mu_r
for i in range(round(10/dt)):
a_kar = F_eff/(m/4)
v_kar = v_kar+a_kar*dt
s_kar = s_kar+v_kar*dt
F_eff = (F_aan-F_rol-Fzt)
F_list.append(F_eff)
v_list.append(v_kar)
s_list.append(s_kar)
t_list.append(dt*i)
#"""
# """
#SCENARIO 2a
Fn = Fz
F_rol = Fn*mu_r
for i in range(round(time/dt)):
if(v_kar < 0):
F_rol=0
F_eff=0
@ -39,6 +60,27 @@ for i in range(round(10/dt)):
v_list.append(v_kar)
s_list.append(s_kar)
t_list.append(dt*i)
#"""
"""
#SCENARIO 3a
Fn = Fz*mp.sin(mp.radians(20))
Fzt = Fz*mp.cos(mp.radians(20))
F_rol = Fn*mu_r
for i in range(round(10/dt)):
a_kar = F_eff/(m/4)
v_kar = v_kar+a_kar*dt
s_kar = s_kar+v_kar*dt
F_eff = (F_aan-F_rol+Fzt)
F_list.append(F_eff)
v_list.append(v_kar)
s_list.append(s_kar)
t_list.append(dt*i)
#"""
plt.plot(t_list,F_list)
plt.plot(t_list,v_list)