IndexError: индекс назначения списка вне диапазона.Как мне это решить?

Question

Я пробую метод постоянного среднего ускорения Ньюмарка.Я получаю эту ошибку.Как мне восстановиться после этой ошибки?

IndexError Traceback (последний последний вызов)

41 для i в диапазоне (len (t)):

42pn [i + 1] = p [i + 1] + a1 * u [i] + a2 * v [i] + a3 * a [i]

43 u [i + 1] = pn [i + 1] / kn

44 v [i + 1] = y * (u [i + 1] -u [i]) / (b * dt) + (1-y / b) v [i] + dt (1-й / (2 * b)) * a [i]

IndexError: индекс назначения списка вне диапазона

y = 1/2
b = 1/4

u = []
v = []

t = []
p = [0,25,43.3013,50,43.3013,25,0,0,0,0,0,0]

a = []

pn = []
pn.append(0)



x = 0.0
for i in range(11):
    z = 0.0 + x
    t.append(z)
    x = x + 0.1

m = 0.45594
k = 18
c = 0.2865

u.append(0)
v.append(0)


a.append((p[0]-c*v[0]-k*u[0])/m)

dt = 0.1

a1 =(m/(b*dt*dt)+y*c/(b*dt))
a2 = (m/(b*dt)+(y/b-1)*c)
a3 = (((1/(2*b))-1)*m + dt*((y/(2*b))-1)*c)
kn = k + a1


for i in range(len(t)-1):
    pn[i+1] = p[i+1]+ a1*u[i] + a2*v[i] + a3*a[i]
    u[i+1] = pn[i+1]/kn
    v[i+1] = y*(u[i+1]-u[i])/(b*dt) + (1-y/b)*v[i] + dt* (1-y/(2*b))*a[i]
    a[i+1] = (u[i+1]-u[i])/(b*dt*dt) - v[i]/(b*dt)-(1/(2*b)-1)*a[i]

Answer 1

Ваш pn, a, u, v определен как список длиной 1, поэтому нет индекса, такого как pn[1].Вы можете использовать append или определить список необходимой длины.

for i in range(len(t)):
    pn.append(p[i+1] + a1*u[i] + a2*v[i] + a3*a[i])
    u.append(pn[i+1]/kn)
    v.append(y*(u[i+1]-u[i])/(b*dt) + (1-y/b)*v[i] + dt* (1-y/(2*b))*a[i])
    a.append((u[i+1]-u[i])/(b*dt*dt) - v[i]/(b*dt)-(1/(2*b)-1)*a[i])

или

pn, a, u, v = [0]*11, [0]*11, [0]*11 [0]*11
pn[0], u[0], v[0] = 0, 0, 0
a[0] = (p[0]-c*v[0]-k*u[0])/m

...

for i in range(len(t)-1):
    pn[i+1] = p[i+1] + a1*u[i] + a2*v[i] + a3*a[i]
    u[i+1] = pn[i+1]/kn
    v[i+1] = y*(u[i+1]-u[i])/(b*dt) + (1-y/b)*v[i] + dt* (1-y/(2*b))*a[i]
    a[i+1] = (u[i+1]-u[i])/(b*dt*dt) - v[i]/(b*dt)-(1/(2*b)-1)*a[i]