Import matplotlib pyplot as plt import numpy as np import

import matplotlib.pyplot as plt
import numpy as np
import math

print("Uch fazali Asinxron dvigatellarning elektr parametrlarini hisoblash")

A = "Maqsad : Uch fazali asinxron dvigatellarning parametrlarini hisoblashni o'rganish."
U = int(input("Kuchlanishni qiymatini kiriting = "))
d = int(input("Qisqa tutashtirilgan rotorli asinxron dvigatelning nominal quvvatini kiriting (Kvt) = "))
t = int(input("Dvigatelni ishlash vaqtini kiriting (soat) = "))
r = float(input("Rotorni aylanish tezligini kiriting (ayl/min) = "))
r1 = float(input("Rotorni nominal aylanish tezligini kiriting (ayl/min) = "))
n = float(input("Dvigatelni nominal foydali ish koeffitsientini kiriting (F.I.K) = "))
c = float(input("Quvvat koeffitsientini kiriting = "))
i = float(input("Ishga tushuruvchi tok karraligi = "))
l = float(input("Zo'riqish koeffitsientini kiriting = "))

# Aniqlash kerak:
# 1) Istemol qiladigan quvvat ;
# 2) Naminal va kritik aylantiruvchi momentlar;
# 3) Ishga tushiruvchi tok;
# 4)Nominal va kritik sirpanishi. Mexanik tavsifini qurish
# 5) Elektr energiya sarfi

# Istemol qiladigan quvvatni topamiz:

Pн = d * math.pow(10, 3) / n
print(f" Dvigatelni ist'emol qiladigan quvvati {Pн / math.pow(10, 3)} (Kvt) ga teng bo'ldi.")

# Nominal va kritik aylantiruvchi momentlarni quyidagicha topiladi:
Mн = (9550 * Pн) / r
print(f" Nominal momentning qiymati {Mн / math.pow(10, 3)} (N*m)ga teng bo'ldi. ")

# Nominal va ishga tushiruvchi toklari:
Iн = Pн / math.sqrt(3) * U * c
print(f" Nominal tokning qiymati {Iн}(A) ga teng bo'ldi. ")
# Nominal sirpanish
Sн = (r - r1) / r
print(f"Nominal sirpanishning qiymati {Sн} ga teng. ")

# Kiritik sirpanish
Sк = Sн * (l + (math.sqrt(l * l - 1)))
print(f"Kritik sirpanishning qiymati {Sк} ga teng ")
print(f"1-chi momentning qiymati {(2 * Mн / ((1 / (10)) / Sк + Sк / (1 / (10))))/1000} ga teng bo'ldi. ")
M1=((2 * Mн / ((1 / (10)) / Sк + Sк / (1 / (10))))/1000)
print(M1)
print(f"2-chi momentning qiymati {(2 * Mн / ((2 / (10)) / Sк + Sк / (2 / (10))))/1000} ga teng bo'ldi. ")
M2=((2 * Mн / ((2 / (10)) / Sк + Sк / (2 / (10))))/1000)
print(M2)
print(f"3-chi momentning qiymati {(2 * Mн / ((3 / (10)) / Sк + Sк / (3 / (10))))/1000} ga teng bo'ldi. ")
M3=(2 * Mн / ((3 / (10)) / Sк + Sк / (3 / (10))))/1000
print(M3)
print(f"4-chi momentning qiymati {(2 * Mн / ((4 / (10)) / Sк + Sк / (4 / (10))))/1000} ga teng bo'ldi. ")
M4=(2 * Mн / ((4 / (10)) / Sк + Sк / (4 / (10))))/1000
print(M4)
print(f"5-chi momentning qiymati {(2 * Mн / ((5 / (10)) / Sк + Sк / (5 / (10))))/1000} ga teng bo'ldi. ")
M5=(2 * Mн / ((5 / (10)) / Sк + Sк / (5 / (10))))/1000
print(M5)
print(f"6-chi momentning qiymati {(2 * Mн / ((6 / (10)) / Sк + Sк / (6 / (10))))/1000} ga teng bo'ldi. ")
M6=(2 * Mн / ((6 / (10)) / Sк + Sк / (6 / (10))))/1000
print(M6)
print(f"7-chi momentning qiymati {(2 * Mн / ((7 / (10)) / Sк + Sк / (7 / (10))))/1000} ga teng bo'ldi. ")
M7=(2 * Mн / ((7 / (10)) / Sк + Sк / (7 / (10))))/1000
print(M7)
print(f"8-chi momentning qiymati {(2 * Mн / ((8 / (10)) / Sк + Sк / (8 / (10))))/1000} ga teng bo'ldi. ")
M8=(2 * Mн / ((8 / (10)) / Sк + Sк / (8 / (10))))/1000
print(M8)
print(f"9-chi momentning qiymati {(2 * Mн / ((9 / (10)) / Sк + Sк / (9 / (10))))/1000} ga teng bo'ldi. ")
M9=(2 * Mн / ((9 / (10)) / Sк + Sк / (9 / (10))))/1000
print(M9)
print(f"10-chi momentning qiymati {(2 * Mн / ((10 / (10)) / Sк + Sк / (10 / (10))))/1000} ga teng bo'ldi. ")
M10=((2 * Mн / ((10 / (10)) / Sк + Sк / (10 / (10))))/1000)
print(M10)

plt.plot([M1, M2, M3, M4, M5, M6, M7, M8, M9, M10], [1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1], 'ro')

plt.xlabel('Moment')
plt.ylabel('Sirpanish')
plt.axis([min(M1, M2, M3, M4, M5, M6, M7, M8, M9, M10), max(M1, M2, M3, M4, M5, M6, M7, M8, M9, M10), 1, 0])
plt.show()

Natija