def simulate(self, xpivot=0.0, vpivot=0.0, t1=1.5, dt=0.01):
# reset
self.system.q [:] = 0.0
self.system.qd[:] = 0.0
# initial conditions
self.system.q [self.pivot.istrain[0]] = xpivot # initial elevation
self.system.qd[self.pivot.istrain[0]] = vpivot # initial elevation spd
self.system.qd[self.axis .istrain[0]] = self.spin # initial rotation speed
# simulate
self.t = np.arange(0, t1, dt)
self.y = solve_system(self.system, self.t, outputs='vels')
def simulate(self, xpivot=0.0, vpivot=0.0, t1=1.5, dt=0.01):
# reset
self.system.q[:] = 0.0
self.system.qd[:] = 0.0
# initial conditions
self.system.q[self.pivot.istrain[0]] = xpivot # initial elevation
self.system.qd[self.pivot.istrain[0]] = vpivot # initial elevation spd
self.system.qd[
self.axis.istrain[0]] = self.spin # initial rotation speed
# simulate
self.t = np.arange(0, t1, dt)
self.y = solve_system(self.system, self.t, outputs='vels')
def sim_both(q0=None, qd0=None):
if q0 is None: q0 = np.zeros_like(w)
if qd0 is None: qd0 = np.zeros_like(w)
t = np.arange(0, 2, 0.005)
ylin = linsys.integrate(t, q0, qd0)
system.q[:] = 0.0
system.qd[:] = 0.0
system.q[system.iFreeDOF] = q0
system.qd[system.iFreeDOF] = qd0
y = solve_system(system, t)
ylin = ylin[:, :len(q0)]
ax = plt.figure().add_subplot(111)
ax.set_color_cycle(['r', 'g', 'b', 'c'])
ax.plot(t, ylin, ':', t, y, '--')
plt.legend(('Mode 1', 'Mode 2', 'Mode 3', 'Mode 4'))
return t, y, ylin
def test1():
# Initial conditions
system.qd[bearing.istrain[0]] = -10 * pi/180
# Prescribed DOFs
#system.prescribe(foundation.istrain, 0.0)
system.prescribe(tower.istrain, 0.0)
system.prescribe(bearing.istrain, 0.0)
system.prescribe(b1.istrain, 0.0)
system.prescribe(b2.istrain, 0.0)
system.prescribe(b3.istrain, 0.0)
# Solve
dt = 0.05
t = np.arange(0, 40, dt)
y = solve_system(system, t)
return t,y
def sim_both(q0=None, qd0=None):
if q0 is None: q0 = np.zeros_like(w)
if qd0 is None: qd0 = np.zeros_like(w)
t = np.arange(0, 2, 0.005)
ylin = linsys.integrate(t, q0, qd0)
system.q [:] = 0.0
system.qd[:] = 0.0
system.q [system.iFreeDOF] = q0
system.qd[system.iFreeDOF] = qd0
y = solve_system(system,t)
ylin = ylin[:,:len(q0)]
ax = plt.figure().add_subplot(111)
ax.set_color_cycle(['r','g','b','c'])
ax.plot(t,ylin,':',t,y,'--')
plt.legend(('Mode 1','Mode 2','Mode 3','Mode 4'))
return t,y,ylin
def simulate(system, t1, dt=0.05):
t = np.arange(0, t1, dt)
y = solve_system(system, t)
return t, y
def simulate(system, t1, dt=0.05):
t = np.arange(0, t1, dt)
y = solve_system(system, t)
return t,y
