Я создал длинный код, состоящий из нескольких списков внутри циклов for.Там нет ничего плохого в расчете. Он получает результаты, как ожидается.Код и построение списков работают нормально.Проблема в том, что когда он выполняется, я определил разрыв, когда определенное условие соответствует.Но он не прерывается при первом запуске, продолжается и проходит через все значения в пределах функции диапазона первого цикла.Я хочу добиться возврата истинного значения, когда условие соответствует и останавливается и не будет продолжать расти в соответствии со значениями диапазона в первом цикле.
Я бы объяснил, как будет работать код!
Код: Первая часть соответствует и является входными данными
import math
import numpy as np
Ned = -500
fcd = 20
fyd = 435
E = 2e5
h = 200
cb = 35
ct = 35
ca = 35
b= 150
y = 12
d = h - cb
ds = ct
a = 25
yb = 8
ecu = 0.0035
rebarnumber = math.floor((b-(2*cb+2*yb+y))/a)
PI_Y2_4 = int(math.pi/4*(y)**2)
disc = []
dis = []
Asi = []
Asci = []
Esc = []
Esci = []
Sc = []
Sci =[]
#############################
# Calculation starts here
#############################
for n in range(0,10): # <------- First for-loop
cbb = cb + yb + y/2
ctt = ct + yb + y/2
if 0 < n <= rebarnumber:
Asi.append(PI_Y2_4)
dis.append( h - cbb)
Asci.append(PI_Y2_4)
disc.append( ctt )
if rebarnumber < n <= (2 * rebarnumber):
Asi.append(PI_Y2_4)
dis.append( h - cbb - ca)
Asci.append(PI_Y2_4)
disc.append(cbb + ca)
if (2*rebarnumber) < n <= (3 * rebarnumber):
Asi.append(PI_Y2_4)
dis.append( h - cbb - 2*ca)
Asci.append(PI_Y2_4)
disc.append(cbb + 2*ca)
if (3*rebarnumber) < n <= (4 * rebarnumber):
Asi.append(PI_Y2_4)
dis.append( h - cbb - 3*ca)
Asci.append(PI_Y2_4)
disc.append(cbb + 3*ca)
if (4*rebarnumber) < n <= (5 * rebarnumber):
Asi.append(PI_Y2_4)
dis.append( h - cbb - 4*ca)
Asci.append(PI_Y2_4)
disc.append(cbb + 4*ca)
if (5*rebarnumber) < n <= (6 * rebarnumber):
Asi.append(PI_Y2_4)
dis.append( h - cbb - 5*ca)
Asci.append(PI_Y2_4)
disc.append(cbb + 5*ca)
if (6*rebarnumber) < n <= (7 * rebarnumber):
Asi.append(PI_Y2_4)
dis.append( h - cbb - 6*ca)
Asci.append(PI_Y2_4)
disc.append(cbb + 6*ca)
if (7*rebarnumber) < n <= (8 * rebarnumber):
Asi.append(PI_Y2_4)
dis.append( h - cbb - 7*ca)
Asci.append(PI_Y2_4)
disc.append(cbb + 7*ca)
for i in range(0,len(dis)):
Esc.insert(i, dis[i])
Esci.insert(i, disc[i])
Sc.insert(i, dis[i])
Sci.insert(i, disc[i])
for x in np.linspace(1,h,10000): # <-------- Second for-loop
for k, _ in enumerate(Esc):
try:
if x < dis[k]:
Esc[k]=( ecu/x*(dis[k]-x) )
else:
Esc[k]=(- ecu/x*(x-dis[k] ) )
if x < disc[k]:
Esci[k]=( -ecu/x*(x-disc[k]) )
else:
Esci[k]=(- ecu/x*(x-disc[k]) )
except (ZeroDivisionError, RuntimeWarning):
Esc[k]=( 0 )
Esci[k]=( 0 )
for k, _ in enumerate(Sc): # <-------- Third for-loop
ss = Esc[k]*E
if ss <= -fyd:
Sc[k]= -fyd
elif ss >= -fyd and ss < 0:
Sc[k]=ss
else:
Sc[k]=min(ss,fyd)
for k, _ in enumerate(Sci):
sci = Esci[k]*E
if sci <= -fyd:
Sci[k]= -fyd
elif sci >= -fyd and sci < 0:
Sci[k]=sci
else:
Sci[k]=min(sci,fyd)
FS = 0
FSC = 0
for a, _ in enumerate(Sc):
FS += Sc[a]*Asi[a]
FSC+=Sci[a]*Asci[a]
MS = 0
MSC = 0
for m, _ in enumerate(Sc):
MS += Sc[a]*Asi[a]*(dis[m]-h/2)
MSC+= Sci[a]*Asci[a]*(h/2-disc[m])
Nrd = 0
Mrd = 0
Nrd = int((-0.8*x*b*fcd+FSC+FS)/1000)
Mrd = (0.8*x*b*fcd*(h/2-0.4*x)+MS-MSC)/1000000
if 0 <= (float(Nrd) - Ned) <= 1:
print(Nrd, x, Asi)
break
break
Как это работает?
Первый цикл for создает список с индексом 0, т. Е. [value1]
, второй цикл for создает приращение значения топора (диапазон), третий цикл forсоздает списки в зависимости от первого созданного списка [value1]
.затем значение Nrd
определяется на основе приращения значений x.если условие 0 <= (float(Nrd) - Ned) <= 1:
выполнено, то вычисление остановится и вернет значение Nrd
.если не совпадает, он возвращается и берет индекс 1 из первого цикла for, [value1, value2]
создается, снова достигает Nrd, если условие выполнено, в противном случае прерывание продолжается до тех пор, пока оно не будет найдено.
Моя проблема в том, что когда код запускается, я получаю эти выводы.
Nrd x Asi
---------------------------------------------------
-499 181.84938493849384 [113]
-499 162.36533653365336 [113, 113]
-499 147.3990399039904 [113, 113, 113]
-499 137.48784878487848 [113, 113, 113, 113]
-499 130.72117211721172 [113, 113, 113, 113, 113]
-499 126.10391039103911 [113, 113, 113, 113, 113, 113]
-499 122.7006700670067 [113, 113, 113, 113, 113, 113, 113]
-499 120.01390139013901 [113, 113, 113, 113, 113, 113, 113, 113]
-499 119.71537153715371 [113, 113, 113, 113, 113, 113, 113, 113, 113]
Выше приведены все истинные множественные решения.Но я хочу остановиться (разбить) первым матчем, вместо того, чтобы давать все решения.
Я использовал return True
, но он здесь не совсем работает, когда жалуется, что он не работает.
Второй вопрос здесь - второй цикл for for x in np.linspace(1,h,10000):
Я действительно хочу запустить его со многими десятичными знаками, чтобы получить лучший результат, но он замедляется и требует много времени для вычисления.Есть ли способ ускорить?
Возможно, определение вышеупомянутых строк кода в функции будет работать более эффективно.