stabilisation efficiency

gistfile1.txt

ssm = 100
S = ssm / (100**2)
sr = 0.03+math.sqrt(ssm)/2/100

r = 0.03
m = 0.3

reff=r/2*3

import math

V = 30
n = 10
rv = 3.14*2*reff*n
dt = 0.0001
angle = 0
t = 0.0
while angle < 3.14:
    cx = 1
    f = cx*1.2*(V**2)*S/2
    f *= math.sin(angle)
    f *= sr/reff
    drv = f/m*dt
    rv -= drv
    dangle = dt*rv / reff
    angle += dangle
    t += dt
    if rv<0.01:
      print('dangle')
      break