pz7 releika

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

# In[343]:


# 7
from math import *
def p(arr):
    for item in arr:
        print(item)
def t1():
    global Usnom,Usnom,Z1wud,cosf,uk,k1ots,k2ots,k3ots,kv,fmch,dt,S3kz
    global ST1nom,ST2nom,Pn1,Pn2,Pn3,Lw1,Lw2,Lw3,Lw4,Ureg,kz,Z3sz5,t3sz3
    global t3sz5,t3sz4
    Usnom = 110   ; Usnom = Usnom*10**3
    Z1wud = 0.4
    cosf = 0.85
    uk = 10.5
    k1ots = 0.85
    k2ots = k1ots
    k3ots = 1.1
    kv = 1.1
    fmch = 65
    dt = 0.5
    ###########
    S3kz = 1330;   S3kz = S3kz * 10**6
    ST1nom = 16;  ST1nom = ST1nom * 10**6
    ST2nom = 32;   ST2nom = ST2nom * 10**6
    Pn1 = 12.3; Pn1 = Pn1 * 10**6
    Pn2 = 21.9; Pn2 = Pn2 * 10**6
    Pn3 = 14.3; Pn3 = Pn3 * 10**6
    Lw1 = 25
    Lw2 = 11
    Lw3 = 10
    Lw4 = 14
    Ureg = 0.09
    kz = 1.14
    Z3sz5 = 16
    t3sz3 = 1.1
    t3sz4 = 0.8
    t3sz5 = 1.2
#################################
def t2():
    global Usnom,Usnom,Z1wud,cosf,uk,k1ots,k2ots,k3ots,kv,fmch,dt,S3kz
    global ST1nom,ST2nom,Pn1,Pn2,Pn3,Lw1,Lw2,Lw3,Lw4,Ureg,kz,Z3sz5,t3sz3
    global t3sz5,t3sz4
    Usnom = 110   ; Usnom = Usnom*10**3
    Z1wud = 0.4
    cosf = 0.85
    uk = 10.5
    k1ots = 0.85
    k2ots = k1ots
    k3ots = 1.1
    kv = 1.1
    fmch = 65
    dt = 0.5
    ###########
    S3kz = 1078;   S3kz = S3kz * 10**6
    ST1nom = 6.3;  ST1nom = ST1nom * 10**6
    ST2nom = 16;   ST2nom = ST2nom * 10**6
    Pn1 = 4.4; Pn1 = Pn1 * 10**6
    Pn2 = 13.2; Pn2 = Pn2 * 10**6
    Pn3 = 16.7; Pn3 = Pn3 * 10**6
    Lw1 = 10
    Lw2 = 20
    Lw3 = 10
    Lw4 = 30
    Ureg = 0.09
    kz = 1.18
    Z3sz5 = 16
    t3sz3 = 1.4
    t3sz4 = 1.9
    t3sz5 = 2
def t3():
    global Usnom,Usnom,Z1wud,cosf,uk,k1ots,k2ots,k3ots,kv,fmch,dt,S3kz
    global ST1nom,ST2nom,Pn1,Pn2,Pn3,Lw1,Lw2,Lw3,Lw4,Ureg,kz,Z3sz5,t3sz3
    global t3sz5,t3sz4
    Usnom = 110   ; Usnom = Usnom*10**3
    Z1wud = 0.4
    cosf = 0.85
    uk = 10.5
    k1ots = 0.85
    k2ots = k1ots
    k3ots = 1.1
    kv = 1.1
    fmch = 65
    dt = 0.5
    ###########
    S3kz = 6000;   S3kz = S3kz * 10**6
    ST1nom = 16;  ST1nom = ST1nom * 10**6
    ST2nom = 5.6;   ST2nom = ST2nom * 10**6
    Pn1 = 12.8; Pn1 = Pn1 * 10**6
    Pn2 = 4.2; Pn2 = Pn2 * 10**6
    Pn3 = 15; Pn3 = Pn3 * 10**6
    Lw1 = 30
    Lw2 = 13
    Lw3 = 18
    Lw4 = 25
    Ureg = 0.09
    kz = 1.25
    Z3sz5 = 16
    t3sz3 = 1.1
    t3sz4 = 0.9
    t3sz5 = 1.5
def t4():
    global Usnom,Usnom,Z1wud,cosf,uk,k1ots,k2ots,k3ots,kv,fmch,dt,S3kz
    global ST1nom,ST2nom,Pn1,Pn2,Pn3,Lw1,Lw2,Lw3,Lw4,Ureg,kz,Z3sz5,t3sz3
    global t3sz5,t3sz4
    Usnom = 110   ; Usnom = Usnom*10**3
    Z1wud = 0.4
    cosf = 0.85
    uk = 10.5
    k1ots = 0.85
    k2ots = k1ots
    k3ots = 1.1
    kv = 1.1
    fmch = 65
    dt = 0.5
    ###########
    S3kz = 2588;   S3kz = S3kz * 10**6
    ST1nom = 20;  ST1nom = ST1nom * 10**6
    ST2nom = 10;   ST2nom = ST2nom * 10**6
    Pn1 = 13.5; Pn1 = Pn1 * 10**6
    Pn2 = 5.9; Pn2 = Pn2 * 10**6
    Pn3 = 6.3; Pn3 = Pn3 * 10**6
    Lw1 = 9
    Lw2 = 23
    Lw3 = 14
    Lw4 = 32
    Ureg = 0.12
    kz = 0.99
    Z3sz5 = 16
    t3sz3 = 0.8
    t3sz4 = 1.4
    t3sz5 = 1
t4()


# In[344]:


# Расчет сопротивлений элементов
Zc = Usnom**2/S3kz
Zc


# In[345]:


Xw1 = Z1wud * Lw1
Xw2 = Z1wud * Lw2
Xw3 = Z1wud * Lw3
Xw4 = Z1wud * Lw4
p([Xw1,Xw2,Xw3,Xw4])


# In[346]:


Xtsr = uk/100 * Usnom**2 / ST1nom
Xtmin = Xtsr * (1 - Ureg)**2
Xtmax = Xtsr * (1 + Ureg)**2
Xt2 =uk/100 * Usnom**2 / ST2nom
p([Xtsr,Xtmin,Xtmax,Xt2])


# In[347]:


#Расчет минимального сопротивления нагрузки
PP = Pn3+Pn2+Pn1
Irabmax3 = Pn2/ 3**0.5 / Usnom / cosf
Irabmax2 = (Pn3+Pn2) / 3**0.5 / Usnom / cosf
Irabmax1 = (Pn3+Pn2+Pn1) / 3**0.5 / Usnom / cosf
p([Irabmax3,Irabmax2,Irabmax1])


# In[348]:


# Мин сопр нагрузки
ZrabMin3 = 0.9 * Usnom / 3**0.5 / Irabmax3
ZrabMin2 = 0.9 * Usnom / 3**0.5 / Irabmax2
ZrabMin1 = 0.9 * Usnom / 3**0.5 / Irabmax1
p([ZrabMin3,ZrabMin2,ZrabMin1])


# In[349]:


# Расчет ступени 1
Z1sz4 = k1ots * Xw4
Z1sz3 = k1ots * (Xw3 + Xt2)
Z1sz2 = k1ots * Xw2
Z1sz1 = k1ots * Xw1
p([Z1sz4,Z1sz3,Z1sz2,Z1sz1])


# In[350]:


# Коэффициент чувствительности I ступени защиты 3 в конце линии:
K1ch = Z1sz3 / Xw3; K1ch


# In[351]:


# Расчет II ступеней
# Расчет ступени 1
Z1sz4 = k1ots * Xw4
Z1sz3 = k1ots * (Xw3 + Xt2)
Z1sz2 = k1ots * Xw2
Z1sz1 = k1ots * Xw1

Z2sz2 = k2ots * (Xw2 + 0.9*Z1sz4)   #min([Z1sz4,Z1sz3,Z1sz2,Z1sz4]))
Z2sz1 = k2ots * (Xw1 + 0.9*Z1sz2)
Z2sz1_ = k2ots * (Xw1 + Xtmin)
Z2sz1__ = min([Z2sz1,Z2sz1_])
#Z2sz1
p([Z2sz2,Z2sz1,Z2sz1_,Z2sz1__])


# In[352]:


# Коэффициент чувствительности II ступени защиты 2:
K2ch2 = Z2sz2 / Xw2; K2ch2
#Коэффициент чувствительности II ступени защиты 1:
K2ch1 = Z2sz1__ / Xw1; K2ch1
p([K2ch2,K2ch1])


# In[353]:


#Расчет III ступеней
t3sz2 = max([t3sz3, t3sz4])+dt;t3sz2
t3sz1 = max([t3sz2,t3sz5]) + dt;t3sz1
p([t3sz2,t3sz1])


# In[354]:


def z(ZrabMin3):
    return (ZrabMin3 / k3ots / kz / kv)
# Сопротивление срабатывания III ступени защиты 3:
Z3sz3 = z(ZrabMin3)
#Сопротивление срабатывания III ступени защиты 2:
Z3sz2 = z(ZrabMin2)
#Сопротивление срабатывания III ступени защиты 1:
Z3sz1 = z(ZrabMin1)
p([Z3sz3,Z3sz2,Z3sz1])


# In[355]:


# Основной и резервный коэффициенты чувствительности III ступени защиты 3:

k3chOsn3 = Z3sz3 / Xw3
k3chRez3 = Z3sz3 / (Xw3 + Xt2)
k3chOsn2 = Z3sz2 / Xw2
k3chRez2 = Z3sz2 / (Xw2 + Xw3)
k3chOsn1 = Z3sz1 / Xw1
k3chRez1 = Z3sz1 / (Xw1 + Xw2)
k3chRezTr1 = Z3sz1 / (Xw1 + Xtmax)

p([k3chOsn3,k3chRez3, k3chOsn2, k3chRez2, k3chOsn1, k3chRez1, k3chRezTr1])


# In[356]:


#Пересчет III ступеней для обеспечения требуемой чувствительности
k1 = 0
k2=0
k3=0
if (k3chOsn3 >=1.5 and k3chRez3 >= 1.25):
    print('ZIIIс.з.3=0 оставим')
else:
    print('пересчет')
    Z3sz3 = ZrabMin1/k3ots/kz/cos(radians(fmch-acos(Pn1/(ST2nom+ST1nom)/0.85/((100+uk)/100))/pi*180))
    print('пересчет',Z3sz3)
    k3=1
if (k3chOsn2 >=1.5 and k3chRez2 >= 1.25):
    print('ZIIIс.з.2=0 оставим')
else:
    Z3sz2 = ZrabMin2/k3ots/kz/cos(radians(fmch-acos(Pn2/ST2nom/0.85/((100+uk)/100))/pi*180))
    print('пересчет',Z3sz2)
    k2=1
if (k3chOsn1 >=1.5 and k3chRez1 >= 1.25 and k3chRezTr1>=1.25 ):
    print('ZIIIс.з.1=0 оставим')
else:
    Z3sz1 = ZrabMin1/k3ots/kz/cos(radians(fmch-acos(Pn1/ST1nom/0.85/((100+uk)/100))/pi*180))
    print('пересчет',Z3sz1)
    k1=1

#    print(cos(radians(41.98)),sin(radians(41.98)))

#acos(Pn1/ST1nom/0.85/((100+uk)/100))/pi*180=


# In[357]:


# перерасчет kch
k3chOsn3 = Z3sz3 / Xw3
k3chRez3 = Z3sz3 / (Xw3 + Xt2)

k3chOsn2 = Z3sz2 / Xw2
k3chRez2 = Z3sz2 / (Xw2 + Xw3)

k3chOsn1 = Z3sz1 / Xw1
k3chRez1 = Z3sz1 / (Xw1 + Xw2)
k3chRezTr1 = Z3sz1 / (Xw1 + Xtmax)
if k3:
    print("#3осн/рез", k3chOsn3,k3chRez3)
else:
    print("#3 = 0")
if k2:
    print("#2осн/рез", k3chOsn2,k3chRez2)
else:
    print("#2 = 0")
if k1:
    print("#1осн/рез", k3chOsn1,k3chRez1,k3chRezTr1)
else:
    print("#1 = 0")


# In[358]:


#Пересчет III ступеней для обеспечения требуемой чувствительности
# Расчет зон срабатывания защит
L1_1 = Z1sz1 / Z1wud
L2_1 = Z2sz1 / Z1wud
L1_2 = Z1sz2 / Z1wud
L2_2 = Z2sz2 / Z1wud

p([L1_1, L2_1, L1_2,L2_2])


# In[359]:


# Охват II ступенью защиты 1 смежной линии:
lw2 = L2_1 - Lw1
lw4 = L2_2 - Lw2
p([lw2,lw4])


# In[360]:


# Оценка применения защиты 1 в качестве основной на линии W1
U3ost = Z1sz1 / ( Zc + Z1sz1 )
p([U3ost])
if U3ost > 0.6:
    print(1)
else:
    print(0)