TRAIETTORIA

cd99f305-32ae-482d-9ace-9d70b7b29e83 5 4 5 false false 0 decimal

&[DATE] &[TIME] - &[FILENAME]

V 10 :

g 9.8 :

y.0 0.8 :

p 14 :

TRAIETTORIA

t 0.001 (5 (0.01 range :

n t length :

n 556

k 1 n range x k el V t k el * : y k el y.0 π V * t k el * p / sin * : 21 line for

M x y augment :

t x V t * :

t y y.0 π V * t * p / sin * :

Grandezze geometriche caratteristiche della troiettoria

curvatura traiettoria

t c t x t diff t y t 2 diff * (t y t diff t x t 2 diff * (- (t x t diff (2^t y t diff (2^+ (32 /^/ :

x 0 > ( x 5 < ( & x c not defined if

raggio di curvatura

t R.c 1 t c / :

x 0 > ( x 5 < ( & x R.c not defined if Parametri caratteristici motocicletta

Aprilia RSV 1000

m 189 :

p 1.418 :

I.r 0.6 :

r.r 0.45 :

I.xG 6.73 :

h.G 0.43 :

ω.r V r.r / :

I.zG 36.4 :

b.G 0.703 :

EDO velocità di rollio RK23

Runge - Kutta 2° e 3° ordine

2 D\032E\ to\032E\ xo\032E\ N\032E\ E\032E\ RK23 x\032E\ xu\032E\ UIn x\032E\ xu\032E\ UnitsOf * : x\032E\ UOut x\032E\ x\032E\ UnitsOf / vectorize : 12 mat k\0324\ 1 xo\032E\ length range : x\0324\ xo\032E\ k\0324\ el : X\0324\ x\0324\ transpose : 13 mat n\0324\ 1 : t\0324\ to\032E\ 1 el : T\0324\ t\0324\ 11 mat : h\0324\ to\032E\ 2 el t\0324\ - N\032E\ / eval : ho\0324\ h\0324\ : 15 mat t\0324\ x\0324\ k\0324\ h\0324\ K\0324\ h\0324\ t\0324\ x\0324\ D\032E\ * eval h\0324\ E\032E\ ho\0324\ UIn x\0324\ E\032E\ x\0324\ UIn + D\032E\ * eval try 11 line : t\0324\ to\032E\ 2 el ≤ k\0324\.1 t\0324\ x\0324\ x\0324\ h\0324\ K\0324\ : \0324\k.2 t\0324\ h\0324\ + x\0324\ k\0324\.1 + k\0324\.1 h\0324\ K\0324\ : \0324\k.3 t\0324\ 0.5 h\0324\ * + x\0324\ 0.25 k\0324\.1 * + 0.25 \0324\k.2 * + \0324\k.2 h\0324\ K\0324\ : t\0324\ trace d\0324\ k\0324\.1 2 \0324\k.3 * - \0324\k.2 + 3 / UOut normi : \0324\q E\032E\ x\0324\ UOut normi ho\0324\ UOut 12 mat max * : d\0324\ \0324\q ≤ t\0324\ t\0324\ h\0324\ + : T\0324\ n\0324\ n\0324\ 1 + : el t\0324\ : 12 mat x\0324\ x\0324\ k\0324\.1 4 \0324\k.3 * + \0324\k.2 + 6 / + eval : X\0324\ n\0324\ k\0324\ el x\0324\ k\0324\ el : 12 mat 21 line if d\0324\ 0 ≠ h\0324\ ho\0324\ 0.9 h\0324\ * \0324\q d\0324\ / 3 nthroot * 12 mat min eval : if 81 line while T\0324\ X\0324\ augment eval 61 line :

t Clear 1

m V 2^* h.G R.c / * φ cos * m g * h.G * φ sin * - I.xG m h.G 2^* + (φ'' * - 0 ≡

ODE

soluzione equazione differenziale φ

t φ D φ 2 el m V 2^* h.G t R.c / * φ 1 el cos * m g * h.G * φ 1 el sin * - I.xG m h.G 2^* + / 21 mat :

t.o t.end 12 mat 0512 mat :

φ.o φ'.o 12 mat 0 0.5 12 mat :

φsol t φ D t.o t.end 12 mat φ.o φ'.o 12 mat 50103 -^RK23 :

φsol rows 146

t φ φsol 1 col φsol 2 col t cinterp :

t φ' φsol 1 col φsol 3 col t cinterp :

IMBARDATA

t α b.G t R.c / atan :

t Δ p t R.c / atan :

t F.sf t α sin t Δ sin / m * V 2^t R.c / * :

M.gr 2 I.r ω.r * (* t φ' * :

M.gr t φ cos *

t Clear 1

t ψ'' t F.sf p * m V 2^* b.G * t R.c / t α cos * - I.zG m b.G 2^* + / :

x 0 > ( x 5 < ( & x ψ'' not defined if

x Y D x ψ'' :

init 0 :

end 4 :

t VelImb init 0 end 100 x Y D dn_GearsBDF :

t ψ' t VelImb 1 col t VelImb 2 col t cinterp :

angolo di rollio φ(t)

x 0 > ( x 5 < ( & x φ 360 2 π * / * not defined if

velocità di rollio φ'(t)

x 0 > ( x 5 < ( & x φ' not defined if

velocità di imbardata ψ'(t)

x 0 > ( x 4 < ( & x ψ' not defined if

t x' t x t diff :

t y' t y t diff :

t m t y' t x' / :

t ψ π 2 / t m atan - (:

t ϕ t ψ' t φ' / atan :

angolo asse di mozzi

coordinate traccia di mozzi nel sistema solidalle alla motocicletta

t x.m 0 :

t y.m t ψ' V * t ψ' (2^t φ' (2^+ / :

coordinate centro di istantanea rotazione nel sistema solidale alla motocicletta

t x.t 0 :

t y.t V t ψ' / :

traccia mozzi e centro istatanea rotazione x 0 > ( x 4 < ( & x y.m not defined if x 0 > ( x 4 < ( & x y.t not defined if

t 0.001 4 0.01 range : j 1 t rows range x.M j el t j el y t j el ψ sin t j el y.m * (- eval : y.M j el t j el x t j el ψ cos t j el y.m * (+ eval : 21 line for 21 line

t 0.001 4 0.01 range : j 1 t rows range x.T j el t j el y t j el ψ sin t j el y.t * (- eval : y.T j el t j el x t j el ψ cos t j el y.t * (+ eval : 21 line for 21 line

Z y.M x.M augment :

T y.T x.T augment :

(coordinate punti allo stesso istante)

indice 144589 range :

Z indice 1212 mat el

M indice 1212 mat el

l t indice el :

l l 1 round vectorize :

funzione plot

data char size clr plot k 1 data rows range r3 k el . : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

data leg size clr plot2 k 1 data rows range l2 k el l k el var2str : for data l2 size clr augment 21 line :

R Z indice 1212 mat el : L R . 15 black plot : LL R l 14 black plot2 : 31 sys

E M indice 1212 mat el : L1 E . 15 black plot : LL1 E l 14 black plot2 : 31 sys

sulla traccia di Mozzi

sulla traiettoria

Z M T 3 1 sys

angolo asse di Mozzi

x 0 > ( x 4 < ( & x ϕ 360 2 π * / * not defined if Z M T L LL L1 LL1 7 1 sys