'float' object is not subscriptable

@AidaR_D · Регистрация: 31.01.2021

Студворк — интернет-сервис помощи студентам

Пытаюсь смоделировать гидроудар в трубопроводе
на строчке J_B = p_1[j-1][i+1] - density*c/F*q[j-1][i+1] + dx*fi_B выходит ошибка
не понимаю в чем дело
p и q одномерные массивы, поэтому приходится сначала для p_1=[] q_1=[] дважды вставить одномерный, чтобы получить двумерный, в конце удаляю ненужные списки

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
p_1=[]
q_1=[]
 
p_2=[]
q_2=[]
for j in range(m+1):
    if j==0:
        p_1.append (p)
        q_1.append (q)
        p_1.append (p)
        q_1.append (q)
    else:
        for i in range(n+1):
            if i == 0:
                speed_B = 4 * q_1[j-1][i+1] /(3.14159*(diameter/1000)**2)
                Re_B = speed_B*diameter/(1000*viscosity)
                # Вычисляем лямбда 
                if Re_B < 2300:
                    lymda_B = 64/Re_B
                elif Re_B < 10*diameter/0.2:
                    lymda_B = 0.3164/Re_B**0.25
                elif Re_B < 500*diameter/0.2:
                    lymda_B = 0.11*(0.2/diameter+68/Re_B)**0.25
                else:
                    lymda_B = 0.11*(0.2/diameter)**0.25
 
                fi_B = - lymda_B * density * speed_B * abs(speed_B)/(diameter*2)
 
                J_B = p_1[j-1][i+1] - density*c/F*q[j-1][i+1] + dx*fi_B
                
                q_3 = (p0 - J_B)*F/(density*c) 
                p_2.append (p0)
                q_2.append (q_3)
 
 
 
        if j == 1:
            p_1.pop()
            q_1.pop()
        p_1.append (p_2)
        q_1.append (q_2)

@Catstail · 31.01.2021, 16:31

Такое впечатление, что p и q - не массивы

@Viktorrus · 31.01.2021, 16:36

Catstail, Код выложен не полный, поэтому остается только гадать. По этой же причине для поиска ошибки здесь не получится использовать отладчик. ТС хочет, что бы подробно изучили его фрагмент кода и нашли ошибку не используя средства отладки.
Я в такие игры не играю.

@AidaR_D · 31.01.2021, 18:18 **[ТС]**

ОЙ, я первый раз здесь, сори
вот весь код

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
import matplotlib.pyplot as plt
 
density = float(input('Введите плотность нефти: '))
viscosity = float(input('Введите вязкость нефти в мм2/с: '))
pipe_flow = float(input('Введите расход нефти в трубопроводе: '))
 
T = 300         # секунд - время моделирования
diameter = 1000 # мм
delta = 10      # мм
length = 4000   # м
x_valve = 3000  # м, координата задвижки 
density_water = 1000
 
pipe_flow /= 3600
diameter -= 2*delta
viscosity /= 1000000
K = 1.5*10**9           #  расчет констант 
E = 2*10**11
G = density/density_water
N = 0.865/(3600*(10**5)**0.5)
F = (3.14159*(diameter/1000)**2)/4  # м2
c = 1/(density/K + density*diameter/(E*delta))**0.5     # м/с
 
dx = 1000   # шаг сетки по x  
dt = dx/c    # шаг по времени 
n = int(length/dx)  
m = int(T//dt)      # n*m - сетка на плоскости t0x
 
# начачальные условия 
# 1) определяем начальное давление в трубопроводе 
speed = 4*pipe_flow/(3.14159*(diameter/1000)**2)
Re = speed*diameter/(1000*viscosity)
lymda = 0.11*(0.2/diameter)**0.25
length_loss = lymda*length*speed**2/(2*diameter/1000*9.81)
p0 = 10**5 + length_loss*density*9.81
 
# 2) распределение давления по длине трубоповода 
p=[]
x=[]
for i in range(n+1):
    x.append (dx * i)
    p.append (p0 - lymda*x[i]*speed**2/(2*diameter/1000*9.81)*density*9.81)
 
# 3) Распределение расхода по длине трубопровода 
q=[]
for i in range(n+1):
    q.append (pipe_flow)
 
# вычисляем давление по длине трубопровода по слоям времени по методу характеристик
t=[]
for j in range(m+1):
    t.append (dt*j)   #сетка по времени 
 
p_1=[]
q_1=[]
valve_opening = 0.1
Cv = 0.50894308*valve_opening**3 - 0.45105658*valve_opening**2 + 0.145726857*valve_opening - 0.002619972
 
p_2=[]
q_2=[]
for j in range(m+1):
    if j==0:
        p_1.append (p)
        q_1.append (q)
        p_1.append (p)
        q_1.append (q)
    else:
        for i in range(n+1):
            if i == 0:
                speed_B = 4 * q_1[j-1][i+1] /(3.14159*(diameter/1000)**2)
                Re_B = speed_B*diameter/(1000*viscosity)
                # Вычисляем лямбда 
                if Re_B < 2300:
                    lymda_B = 64/Re_B
                elif Re_B < 10*diameter/0.2:
                    lymda_B = 0.3164/Re_B**0.25
                elif Re_B < 500*diameter/0.2:
                    lymda_B = 0.11*(0.2/diameter+68/Re_B)**0.25
                else:
                    lymda_B = 0.11*(0.2/diameter)**0.25
 
                fi_B = - lymda_B * density * speed_B * abs(speed_B)/(diameter*2)
 
                J_B = p_1[j-1][i+1] - density*c/F*q[j-1][i+1] + dx*fi_B
                
                q_3 = (p0 - J_B)*F/(density*c) 
                p_2.append (p0)
                q_2.append (q_3)
 
 
 
        if j == 1:
            p_1.pop()
            q_1.pop()
        p_1.append (p_2)
        q_1.append (q_2)
 
 
print (p_1)

@Welemir1 · 31.01.2021, 18:29

Сообщение от AidaR_D

q[j-1][i+1]

q может и список, но содержит то он не списки (не путать с q_1)

@AidaR_D · 31.01.2021, 19:23 **[ТС]**

Такая дурацкая ошибка.. Спасибо !!

@Viktorrus · 31.01.2021, 19:40

AidaR_D, Я не знаю какие у Вас входные данные, поэтому взял ориентировочно следующие

Python
1
2
3
Введите плотность нефти: 500
Введите вязкость нефти в мм2/с: 2000000
Введите расход нефти в трубопроводе: 7200

Строка 84

Python
1
J_B = p_1[j-1][i+1] - density*c/F*q[j-1][i+1] + dx*fi_B

выдает ошибку, так как список q одномерный ( в моем примере он [2.0, 2.0, 2.0, 2.0, 2.0] ) и поэтому к нему нельзя применять двойную индексацию [j-1][i+1].
Поэтому или разбирайтесь, почему элементами списка q являются числа float, а не списки или нужно оставить только один индекс q[...].
Если не использовать отладчик, то убедится в моих словах можно вставив перед указанной строкой
print(q), и увидеть, что находится в списке q .

Кликните здесь для просмотра всего текста

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
#import matplotlib.pyplot as plt
 
density = float(input('Введите плотность нефти: '))
viscosity = float(input('Введите вязкость нефти в мм2/с: '))
pipe_flow = float(input('Введите расход нефти в трубопроводе: '))
 
T = 300         # секунд - время моделирования
diameter = 1000 # мм
delta = 10      # мм
length = 4000   # м
x_valve = 3000  # м, координата задвижки 
density_water = 1000
 
pipe_flow /= 3600
diameter -= 2*delta
viscosity /= 1000000
K = 1.5*10**9           #  расчет констант 
E = 2*10**11
G = density/density_water
N = 0.865/(3600*(10**5)**0.5)
F = (3.14159*(diameter/1000)**2)/4  # м2
c = 1/(density/K + density*diameter/(E*delta))**0.5     # м/с
 
dx = 1000   # шаг сетки по x  
dt = dx/c    # шаг по времени 
n = int(length/dx)  
m = int(T//dt)      # n*m - сетка на плоскости t0x
 
# начачальные условия 
# 1) определяем начальное давление в трубопроводе 
speed = 4*pipe_flow/(3.14159*(diameter/1000)**2)
Re = speed*diameter/(1000*viscosity)
lymda = 0.11*(0.2/diameter)**0.25
length_loss = lymda*length*speed**2/(2*diameter/1000*9.81)
p0 = 10**5 + length_loss*density*9.81
 
# 2) распределение давления по длине трубоповода 
p=[]
x=[]
for i in range(n+1):
    x.append (dx * i)
    p.append (p0 - lymda*x[i]*speed**2/(2*diameter/1000*9.81)*density*9.81)
 
# 3) Распределение расхода по длине трубопровода 
q=[]
for i in range(n+1):
    q.append (pipe_flow)
 
# вычисляем давление по длине трубопровода по слоям времени по методу характеристик
t=[]
for j in range(m+1):
    t.append (dt*j)   #сетка по времени 
 
p_1=[]
q_1=[]
valve_opening = 0.1
Cv = 0.50894308*valve_opening**3 - 0.45105658*valve_opening**2 + 0.145726857*valve_opening - 0.002619972
 
p_2=[]
q_2=[]
for j in range(m+1):
    if j==0:
        p_1.append (p)
        q_1.append (q)
        p_1.append (p)
        q_1.append (q)
    else:
        for i in range(n+1):
            if i == 0:
                speed_B = 4 * q_1[j-1][i+1] /(3.14159*(diameter/1000)**2)
                Re_B = speed_B*diameter/(1000*viscosity)
                # Вычисляем лямбда 
                if Re_B < 2300:
                    lymda_B = 64/Re_B
                elif Re_B < 10*diameter/0.2:
                    lymda_B = 0.3164/Re_B**0.25
                elif Re_B < 500*diameter/0.2:
                    lymda_B = 0.11*(0.2/diameter+68/Re_B)**0.25
                else:
                    lymda_B = 0.11*(0.2/diameter)**0.25
 
                fi_B = - lymda_B * density * speed_B * abs(speed_B)/(diameter*2)
 
 
                print(q)
 
 
                J_B = p_1[j-1][i+1] - density*c/F*q[j-1][i+1] + dx*fi_B
                
                q_3 = (p0 - J_B)*F/(density*c) 
                p_2.append (p0)
                q_2.append (q_3)
 
 
 
        if j == 1:
            p_1.pop()
            q_1.pop()
        p_1.append (p_2)
        q_1.append (q_2)
 
 
print (p_1)

В моем примере выводит:

Python
1
2
3
4
Введите плотность нефти: 500
Введите вязкость нефти в мм2/с: 2000000
Введите расход нефти в трубопроводе: 7200
[2.0, 2.0, 2.0, 2.0, 2.0]

Добавлено через 7 минут
Вот он и ругается, потому, что к числу float не применим индекс. float[i] выдает ошибку.

Добавлено через 1 минуту
Пока я разбирался, Вам уже все объяснили.

Новые блоги и статьи Все статьи Все блоги /
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Access VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .

@Catstail Супер-модератор 38161 / 21096 / 4306 Регистрация: 12.02.2012 Сообщений: 34,680 Записей в блоге: 14
	31.01.2021, 16:31
	Такое впечатление, что p и q - не массивы 0

@Viktorrus 1732 / 970 / 199 Регистрация: 22.02.2018 Сообщений: 2,693 Записей в блоге: 6
	31.01.2021, 16:36
	Catstail, Код выложен не полный, поэтому остается только гадать. По этой же причине для поиска ошибки здесь не получится использовать отладчик. ТС хочет, что бы подробно изучили его фрагмент кода и нашли ошибку не используя средства отладки. Я в такие игры не играю. 3

@AidaR_D 0 / 0 / 0 Регистрация: 31.01.2021 Сообщений: 4
	31.01.2021, 19:23 [ТС]
	Такая дурацкая ошибка.. Спасибо !! 0

'float' object is not subscriptable

Решение