Skip to content

Commit c95e221

Browse files
PMeiraPMeira
committed
Add a simple example for the IEEE 13 bus system
1 parent 6a9e215 commit c95e221

File tree

4 files changed

+502
-0
lines changed

4 files changed

+502
-0
lines changed

‎examples/13Bus/IEEE13Node_BusXY.csv

Lines changed: 18 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,18 @@
1+
SourceBus, 200, 400
2+
650, 200, 350
3+
RG60, 200, 300
4+
646, 0, 250
5+
645, 100, 250
6+
632, 200, 250
7+
633, 350, 250
8+
634, 400, 250
9+
670, 200, 200
10+
611, 0, 100
11+
684, 100, 100
12+
671, 200, 100
13+
692, 250, 100
14+
675, 400, 100
15+
652, 100, 0
16+
680, 200, 0
17+
18+

‎examples/13Bus/IEEE13Nodeckt.dss

Lines changed: 177 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,177 @@
1+
Clear
2+
Set DefaultBaseFrequency=60
3+
4+
!
5+
! This script is based on a script developed by Tennessee Tech Univ students
6+
! Tyler Patton, Jon Wood, and David Woods, April 2009
7+
!
8+
9+
new circuit.IEEE13Nodeckt
10+
~ basekv=115 pu=1.0001 phases=3 bus1=SourceBus
11+
~ Angle=30 ! advance angle 30 deg so result agree with published angle
12+
~ MVAsc3=20000 MVASC1=21000 ! stiffen the source to approximate inf source
13+
14+
15+
16+
!SUB TRANSFORMER DEFINITION
17+
! Although this data was given, it does not appear to be used in the test case results
18+
! The published test case starts at 1.0 per unit at Bus 650. To make this happen, we will change the impedance
19+
! on the transformer to something tiny by dividing by 1000 using the DSS in-line RPN math
20+
New Transformer.Sub Phases=3 Windings=2 XHL=(8 1000 /)
21+
~ wdg=1 bus=SourceBus conn=delta kv=115 kva=5000 %r=(.5 1000 /)
22+
~ wdg=2 bus=650 conn=wye kv=4.16 kva=5000 %r=(.5 1000 /)
23+
24+
! FEEDER 1-PHASE VOLTAGE REGULATORS
25+
! Define low-impedance 2-wdg transformer
26+
27+
New Transformer.Reg1 phases=1 bank=reg1 XHL=0.01 kVAs=[1666 1666]
28+
~ Buses=[650.1 RG60.1] kVs=[2.4 2.4] %LoadLoss=0.01
29+
new regcontrol.Reg1 transformer=Reg1 winding=2 vreg=122 band=2 ptratio=20 ctprim=700 R=3 X=9
30+
31+
New Transformer.Reg2 phases=1 bank=reg1 XHL=0.01 kVAs=[1666 1666]
32+
~ Buses=[650.2 RG60.2] kVs=[2.4 2.4] %LoadLoss=0.01
33+
new regcontrol.Reg2 transformer=Reg2 winding=2 vreg=122 band=2 ptratio=20 ctprim=700 R=3 X=9
34+
35+
New Transformer.Reg3 phases=1 bank=reg1 XHL=0.01 kVAs=[1666 1666]
36+
~ Buses=[650.3 RG60.3] kVs=[2.4 2.4] %LoadLoss=0.01
37+
new regcontrol.Reg3 transformer=Reg3 winding=2 vreg=122 band=2 ptratio=20 ctprim=700 R=3 X=9
38+
39+
40+
!TRANSFORMER DEFINITION
41+
New Transformer.XFM1 Phases=3 Windings=2 XHL=2
42+
~ wdg=1 bus=633 conn=Wye kv=4.16 kva=500 %r=.55 XHT=1
43+
~ wdg=2 bus=634 conn=Wye kv=0.480 kva=500 %r=.55 XLT=1
44+
45+
46+
!LINE CODES
47+
redirect IEEELineCodes.dss
48+
49+
// these are local matrix line codes
50+
// corrected 9-14-2011
51+
New linecode.mtx601 nphases=3 BaseFreq=60
52+
~ rmatrix = (0.3465 | 0.1560 0.3375 | 0.1580 0.1535 0.3414 )
53+
~ xmatrix = (1.0179 | 0.5017 1.0478 | 0.4236 0.3849 1.0348 )
54+
~ units=mi
55+
New linecode.mtx602 nphases=3 BaseFreq=60
56+
~ rmatrix = (0.7526 | 0.1580 0.7475 | 0.1560 0.1535 0.7436 )
57+
~ xmatrix = (1.1814 | 0.4236 1.1983 | 0.5017 0.3849 1.2112 )
58+
~ units=mi
59+
New linecode.mtx603 nphases=2 BaseFreq=60
60+
~ rmatrix = (1.3238 | 0.2066 1.3294 )
61+
~ xmatrix = (1.3569 | 0.4591 1.3471 )
62+
~ units=mi
63+
New linecode.mtx604 nphases=2 BaseFreq=60
64+
~ rmatrix = (1.3238 | 0.2066 1.3294 )
65+
~ xmatrix = (1.3569 | 0.4591 1.3471 )
66+
~ units=mi
67+
New linecode.mtx605 nphases=1 BaseFreq=60
68+
~ rmatrix = (1.3292 )
69+
~ xmatrix = (1.3475 )
70+
~ units=mi
71+
72+
/*********** Original 606 Linecode *********************
73+
74+
You have to use this to match Kersting's results:
75+
76+
New linecode.mtx606 nphases=3 BaseFreq=60
77+
~ rmatrix = (0.7982 | 0.3192 0.7891 | 0.2849 0.3192 0.7982 )
78+
~ xmatrix = (0.4463 | 0.0328 0.4041 | -0.0143 0.0328 0.4463 )
79+
~ Cmatrix = [257 | 0 257 | 0 0 257] ! <--- This is too low by 1.5
80+
~ units=mi
81+
82+
Corrected mtx606 Feb 3 2016 by RDugan
83+
84+
The new LineCode.606 is computed using the following CN cable definition and
85+
LineGeometry definition:
86+
87+
New CNDATA.250_1/3 k=13 DiaStrand=0.064 Rstrand=2.816666667 epsR=2.3
88+
~ InsLayer=0.220 DiaIns=1.06 DiaCable=1.16 Rac=0.076705 GMRac=0.20568 diam=0.573
89+
~ Runits=kft Radunits=in GMRunits=in
90+
91+
New LineGeometry.606 nconds=3 nphases=3 units=ft
92+
~ cond=1 cncable=250_1/3 x=-0.5 h= -4
93+
~ cond=2 cncable=250_1/3 x=0 h= -4
94+
~ cond=3 cncable=250_1/3 x=0.5 h= -4
95+
96+
****End Comment******/
97+
98+
New Linecode.mtx606 nphases=3 Units=mi
99+
~ Rmatrix=[0.791721 |0.318476 0.781649 |0.28345 0.318476 0.791721 ]
100+
~ Xmatrix=[0.438352 |0.0276838 0.396697 |-0.0184204 0.0276838 0.438352 ]
101+
~ Cmatrix=[383.948 |0 383.948 |0 0 383.948 ]
102+
New linecode.mtx607 nphases=1 BaseFreq=60
103+
~ rmatrix = (1.3425 )
104+
~ xmatrix = (0.5124 )
105+
~ cmatrix = [236]
106+
~ units=mi
107+
108+
109+
!LOAD DEFINITIONS
110+
New Load.671 Bus1=671.1.2.3 Phases=3 Conn=Delta Model=1 kV=4.16 kW=1155 kvar=660
111+
New Load.634a Bus1=634.1 Phases=1 Conn=Wye Model=1 kV=0.277 kW=160 kvar=110
112+
New Load.634b Bus1=634.2 Phases=1 Conn=Wye Model=1 kV=0.277 kW=120 kvar=90
113+
New Load.634c Bus1=634.3 Phases=1 Conn=Wye Model=1 kV=0.277 kW=120 kvar=90
114+
New Load.645 Bus1=645.2 Phases=1 Conn=Wye Model=1 kV=2.4 kW=170 kvar=125
115+
New Load.646 Bus1=646.2.3 Phases=1 Conn=Delta Model=2 kV=4.16 kW=230 kvar=132
116+
New Load.692 Bus1=692.3.1 Phases=1 Conn=Delta Model=5 kV=4.16 kW=170 kvar=151
117+
New Load.675a Bus1=675.1 Phases=1 Conn=Wye Model=1 kV=2.4 kW=485 kvar=190
118+
New Load.675b Bus1=675.2 Phases=1 Conn=Wye Model=1 kV=2.4 kW=68 kvar=60
119+
New Load.675c Bus1=675.3 Phases=1 Conn=Wye Model=1 kV=2.4 kW=290 kvar=212
120+
New Load.611 Bus1=611.3 Phases=1 Conn=Wye Model=5 kV=2.4 kW=170 kvar=80
121+
New Load.652 Bus1=652.1 Phases=1 Conn=Wye Model=2 kV=2.4 kW=128 kvar=86
122+
New Load.670a Bus1=670.1 Phases=1 Conn=Wye Model=1 kV=2.4 kW=17 kvar=10
123+
New Load.670b Bus1=670.2 Phases=1 Conn=Wye Model=1 kV=2.4 kW=66 kvar=38
124+
New Load.670c Bus1=670.3 Phases=1 Conn=Wye Model=1 kV=2.4 kW=117 kvar=68
125+
126+
!CAPACITOR DEFINITIONS
127+
New Capacitor.Cap1 Bus1=675 phases=3 kVAR=600 kV=4.16
128+
New Capacitor.Cap2 Bus1=611.3 phases=1 kVAR=100 kV=2.4
129+
130+
!Bus 670 is the concentrated point load of the distributed load on line 632 to 671 located at 1/3 the distance from node 632
131+
132+
!LINE DEFINITIONS
133+
New Line.650632 Phases=3 Bus1=RG60.1.2.3 Bus2=632.1.2.3 LineCode=mtx601 Length=2000 units=ft
134+
New Line.632670 Phases=3 Bus1=632.1.2.3 Bus2=670.1.2.3 LineCode=mtx601 Length=667 units=ft
135+
New Line.670671 Phases=3 Bus1=670.1.2.3 Bus2=671.1.2.3 LineCode=mtx601 Length=1333 units=ft
136+
New Line.671680 Phases=3 Bus1=671.1.2.3 Bus2=680.1.2.3 LineCode=mtx601 Length=1000 units=ft
137+
New Line.632633 Phases=3 Bus1=632.1.2.3 Bus2=633.1.2.3 LineCode=mtx602 Length=500 units=ft
138+
New Line.632645 Phases=2 Bus1=632.3.2 Bus2=645.3.2 LineCode=mtx603 Length=500 units=ft
139+
New Line.645646 Phases=2 Bus1=645.3.2 Bus2=646.3.2 LineCode=mtx603 Length=300 units=ft
140+
New Line.692675 Phases=3 Bus1=692.1.2.3 Bus2=675.1.2.3 LineCode=mtx606 Length=500 units=ft
141+
New Line.671684 Phases=2 Bus1=671.1.3 Bus2=684.1.3 LineCode=mtx604 Length=300 units=ft
142+
New Line.684611 Phases=1 Bus1=684.3 Bus2=611.3 LineCode=mtx605 Length=300 units=ft
143+
New Line.684652 Phases=1 Bus1=684.1 Bus2=652.1 LineCode=mtx607 Length=800 units=ft
144+
145+
146+
!SWITCH DEFINITIONS
147+
New Line.671692 Phases=3 Bus1=671 Bus2=692 Switch=y r1=1e-4 r0=1e-4 x1=0.000 x0=0.000 c1=0.000 c0=0.000
148+
149+
Set Voltagebases=[115, 4.16, .48]
150+
calcv
151+
Solve
152+
BusCoords IEEE13Node_BusXY.csv
153+
154+
!---------------------------------------------------------------------------------------------------------------------------------------------------
155+
!----------------Show some Results -----------------------------------------------------------------------------------------------------------------
156+
!---------------------------------------------------------------------------------------------------------------------------------------------------
157+
158+
159+
// Show Voltages LN Nodes
160+
// Show Currents Elem
161+
// Show Powers kVA Elem
162+
// Show Losses
163+
// Show Taps
164+
165+
!---------------------------------------------------------------------------------------------------------------------------------------------------
166+
!---------------------------------------------------------------------------------------------------------------------------------------------------
167+
! Alternate Solution Script
168+
! To force the taps to be same as published results, set the transformer taps manually and disable the controls
169+
!---------------------------------------------------------------------------------------------------------------------------------------------------
170+
/*
171+
Transformer.Reg1.Taps=[1.0 1.0625]
172+
Transformer.Reg2.Taps=[1.0 1.0500]
173+
Transformer.Reg3.Taps=[1.0 1.06875]
174+
Set Controlmode=OFF
175+
176+
Solve
177+
*/

0 commit comments

Comments
 (0)