grafKulaga

#!/usr/bin/env python
# coding: utf-8

# In[75]:


import math
Xd = 1.8
Xq = 1
Xid = 0.3
Ug = 1
s = 1
cosfi = 0.85
# naiti Id Iq Eq Eiq EQ
#Xd = 2.04
#Xq = 0.189
#Xid = 0.271
#Ug = 1
#s = 1
#cosfi = 0.8


# In[76]:


Pg = s * cosfi
Qg = s * (1-cosfi*cosfi)**0.5
EQ = ( (Ug + Qg*Xq/Ug)**2 + (Pg*Xq/Ug)**2 )**0.5


# In[77]:


z = Pg*Xq/(Ug**2+Qg*Xq)
Lamda = math.atan(z)/math.pi*180
print(Lamda)


# In[78]:


Id =( Ug*math.cos(Lamda/180*math.pi) - EQ) / Xq
Iq = U*math.sin(Lamda/180*math.pi) / Xq
Eq = EQ - (Xd - Xq)*Id
Eiq = EQ - (Xid - Xq)*(Id)


# In[79]:


print(Pg,Qg,EQ,Lamda,Id, Iq, Eq, Eiq)


# In[94]:


import numpy as pn
import matplotlib.pyplot as plt
import seaborn as sns


# In[95]:


sns.set()


# In[160]:


ss=-0.01

x=-Eq*math.cos(Lamda); y=-Eq*math.sin(Lamda);  plt.text(x, y, "Eq")
plt.quiver(0,0,x,y, scale_units='xy', angles='xy', scale=1, color='r')

x=-EQ*math.cos(Lamda); y=-EQ*math.sin(Lamda);  plt.text(x, y, "EQ")
plt.quiver(0,0,x,y, scale_units='xy', angles='xy', scale=1, color='g')

x = Ug; y = 0; plt.text(x*1, y*1.3, "U")
plt.quiver(0,0,x,y, scale_units='xy', angles='xy', scale=1, color='b')
x = Ug*math.cos(Lamda)**2; y = Ug*math.cos(Lamda)*math.sin(Lamda);  plt.text(x*0.8, y*1.3, "Uq")
plt.quiver(0,0,x,y, scale_units='xy', angles='xy', scale=1, color='b')
x = Ug*math.cos(math.pi/2 - Lamda)*math.sin(Lamda); y =- Ug*math.cos(math.pi/2-Lamda)*math.cos(Lamda);  plt.text(x*0.8, y*1.3, "Ud")
plt.quiver(ss,ss,x+ss,y+ss, scale_units='xy', angles='xy', scale=1, color='b')


I = (Iq**2+Id**2)**0.5
x=-Iq*math.cos(Lamda); y=-Iq*math.sin(Lamda);  plt.text(x*0.8, y*1.3, "Iq")
plt.quiver(ss,ss,x+ss,y+ss, scale_units='xy', angles='xy', scale=1, color='y')
x=Id*math.sin(Lamda); y=-Id*math.cos(Lamda);  plt.text(x*0.8, y*1.3, "Id")
plt.quiver(0,0,x,y, scale_units='xy', angles='xy', scale=1, color='y')
#y=I*(1-cosfi*cosfi)**0.5*math.cos(Lamda); x=I*cosfi;  plt.text(x, y, "I")
#plt.quiver(0,0,x,y, scale_units='xy', angles='xy', scale=1, color='y')
x = Id*math.sin(Lamda)-Iq*math.cos(Lamda); y=-Iq*math.sin(Lamda)-Id*math.cos(Lamda);  plt.text(x, y, "I")
plt.quiver(0,0,x,y, scale_units='xy', angles='xy', scale=1, color='y')


plt.title("Иллюстративная диаграмма синхронной машины")
plt.xlim(-0.1,2)
plt.ylim(-1,2)
plt.show