-
Notifications
You must be signed in to change notification settings - Fork 0
/
case30_cssr.py
85 lines (82 loc) · 2.29 KB
/
case30_cssr.py
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
import numpy as np
# tg = np.tan(np.arccos(0.9))
CASEDATA = {
'Vs': 11, # voltage at node 1 (kV)
'Vmin': 10.5, # minimum voltage (kV)
'alpha': 0.601, # load factor
'g': 0.04, # capital recovery rate 2% kamata na 30 godini
'cl': 0.06, # cost of energy loss ($/kWh)
'delta': 50, # precentage deviation in load Smax*(1'+/-' delta)
'branch': np.array([
# from to length(km) P(kW) Q(kVA)
[0, 1, 1.5, 35.0, 30],
[1, 2, 1.0, 25.0, 15],
[2, 3, 0.5, 10.0, 6.0],
[3, 4, 2.5, 25, 15],
[4, 5, 0.5, 38, 28],
[5, 6, 0.5, 10, 5],
[6, 7, 1, 15, 12],
[7, 8, 1.5, 30, 30],
[3, 9, 2, 40, 30],
[9, 10, 0.5, 54, 30],
[10, 11, 0.25, 30, 15],
[11, 12, 0.25, 10, 5],
[5, 13, 0.5, 25, 20],
[13, 14, 1, 40, 30],
[14, 15, 1, 40, 30],
[15, 16, 0.5, 100, 90],
[15, 17, 0.5, 60, 30],
[17, 18, 0.5, 30, 30],
[4, 19, 0.5, 10, 5],
[19, 20, 1, 40, 40],
[20, 21, 1.5, 40, 30],
[21, 22, 0.5, 27, 20],
[20, 23, 1, 80, 70],
[23, 24, 1.1, 45, 39],
[24, 25, 1, 25, 20],
[25, 26, 1, 20, 10],
[10, 27, 0.5, 40, 30],
[27, 28, 0.5, 30, 20],
[28, 29, 1.5, 22, 12],
]),
'cable': np.array([
# area(mm^2) r(ohm/km) x(ohm/km) Imax(A) instalation($/km)
[12.9, 1.376, 0.3896, 115, 2600], #Squirrel
[31.6, 0.9108, 0.3797, 150, 4000],#Weasel
[50, 0.5441, 0.3973, 208, 6400],#Rabbit
[80, 0.3675, 0.3579, 270, 9600],#Raccoon
]),
'xy': np.array([
# x(mm) y(mm)
[0, 4], #1
[1, 4],#2
[2, 4],#3
[3, 4],#4
[4, 4],#5
[5, 4],#6
[6, 4],#7
[7, 4],#8
[8, 4],#9
[3, 3],#10
[3, 2],#11
[3, 1],#12
[3, 0],#13
[5, 3],#14
[6, 3],#15
[7, 3],#16
[8, 3],#17
[7, 2],#18
[7, 1],#19
[5, 5],#20
[5, 6],#21
[5, 7],#22
[5, 8],#23
[5, 6],#24
[6, 6],#25
[7, 6],#26
[8, 6],#27
[4, 2],#28
[5, 2],#29
[6, 2],#30
])
}