2

Haris Hunafa Hanifan / 1610631160064

PENYELESAIAN DENGAN MENGGUNAKAN METODE

GAUSS-SHEIDL

1 2
0.02 + 𝑗0.04
256.6
MW

110.2
0.01 + 𝑗0.03 0.0125 + 𝑗0.025 Mvar
Slack Bus 3
𝑉1 = 1.05∠0°
138.6 45.2
MW Mvar

Konversikan saluran impedansi ke dalam admittance

1
𝑦12 = = 10 − 𝑗20
0.02 + 𝑗0.04

1
𝑦13 = = 10 − 𝑗30
0.01 + 𝑗0.03

1
𝑦23 = = 16 − 𝑗32
0.0125 + 𝑗0.025

Ubahlah nilai P-Q bus kedalam satuan 𝒑𝒖

(256.6 + 𝑗110.2)
𝑆2𝑠𝑐ℎ = − = −2.566 − 𝑗1.102 𝑝𝑢
100

(138.6 + 𝑗45.2)
𝑆3𝑠𝑐ℎ = − = −1.386 − 𝑗0.452 𝑝𝑢
100

Nilai tegangan pada V2 dan V3 tidak diketahui, maka estimasikan

(0) (0)
𝑉2 = 1.0 + 𝑗0.0 𝑑𝑎𝑛 𝑉3 = 1.0 + 𝑗0.0
Sehingga dapat diperoleh:

𝑃2𝑠𝑐ℎ − 𝑗𝑄2𝑠𝑐ℎ (0)

∗(0) + 𝑦12 𝑉1 + 𝑦23 𝑉3
(1) 𝑉2
𝑉2 =
𝑦12 + 𝑦23
 Haris Hunafa Hanifan / 1610631160064

1 2
𝑦12 = 10 − 𝑗20
256.6
MW

110.2
𝑦13 = 10 − 𝑗30 𝑦23 = 16 − 𝑗32
Mvar
Slack Bus 3
𝑉1 = 1.05∠0°
138.6 45.2
MW Mvar
−2.566+𝑗1.102
+(10−𝑗20)(1.05+𝑗0)+(16−𝑗32)(1.0+𝑗0)
(1) 1.0−𝑗0
𝑉2 = (26−𝑗52)

= 𝟎. 𝟗𝟖𝟐𝟓 − 𝒋𝟎. 𝟎𝟑𝟏𝟎

𝑃𝑠𝑐ℎ 𝑠𝑐ℎ
3 −𝑗𝑄3 +𝑦 𝑉 +𝑦 𝑉 (1)
∗(0) 13 1 23 2
(1) 𝑉3
𝑉3 =
𝑦13 +𝑦23

−1.386 + 𝑗0.452
+ (10 − 𝑗30)(1.05 + 𝑗0) + (16 − 𝑗32)(0.9825 − 𝑗0.0310)
1 − 𝑗0
=
(26 − 𝑗62)

= 𝟏. 𝟎𝟎𝟏𝟏 − 𝒋𝟎. 𝟎𝟑𝟓𝟑

𝑃𝑠𝑐ℎ 𝑠𝑐ℎ
2 −𝑗𝑄2 +𝑦 𝑉 +𝑦 𝑉 (1)
∗(0) 12 1 23 3
(2) 𝑉2
𝑉2 = 𝑦12 +𝑦23

−2.566 + 𝑗1.102
1.0 − 𝑗0 + (10 − 𝑗20)(1.05 + 𝑗0) + (1.0011 − 𝑗0.0353)
(1)
𝑉2 =
(26 − 𝑗52)
= 𝟎. 𝟗𝟖𝟏𝟔 − 𝒋𝟎. 𝟎𝟓𝟐𝟎
 Haris Hunafa Hanifan / 1610631160064

𝑃𝑠𝑐ℎ 𝑠𝑐ℎ
3 −𝑗𝑄3 +𝑦 𝑉 +𝑦 𝑉 (2)
∗(0) 13 1 23 2
(2) 𝑉3
𝑉3 = 𝑦13 +𝑦23

−1.386 + 𝑗0.452
+ (10 − 𝑗30)(1.05 + 𝑗0) + (16 − 𝑗32)(0.9816 − 𝑗0.0520)
1 − 𝑗0
=
(26 − 𝑗62)

= 𝟏. 𝟎𝟎𝟎𝟖 − 𝒋𝟎. 𝟎𝟒𝟓𝟗

𝑃𝑠𝑐ℎ 𝑠𝑐ℎ
2 −𝑗𝑄2 +𝑦 𝑉 +𝑦 𝑉 (2)
∗(0) 12 1 23 3
(3) 𝑉2
𝑉2 = 𝑦12 +𝑦23

−2.566 + 𝑗1.102
+ (10 − 𝑗20)(1.05 + 𝑗0) + (1.0008 − 𝑗0.0459)
(3) 1.0 − 𝑗0
𝑉2 =
(26 − 𝑗52)

= 𝟎. 𝟗𝟖𝟎𝟖 − 𝒋𝟎. 𝟎𝟓𝟕𝟖

𝑃𝑠𝑐ℎ 𝑠𝑐ℎ
3 −𝑗𝑄3 +𝑦 𝑉 +𝑦 𝑉 (3)
∗(0) 13 1 23 2
(3) 𝑉3
𝑉3 = 𝑦13 +𝑦23

−1.386 + 𝑗0.452
+ (10 − 𝑗30)(1.05 + 𝑗0) + (16 − 𝑗32)(0.9808 − 𝑗0.0578)
1 − 𝑗0
=
(26 − 𝑗62)

= 𝟏. 𝟎𝟎𝟎𝟒 − 𝒋𝟎. 𝟎𝟒𝟖𝟖

(4)
𝑉2 = 0.9803 − 𝑗0.0594
(4)
𝑉3 = 1.0002 − 𝑗0.0497
(5)
𝑉2 = 0.9801 − 𝑗0.0598
(5)
𝑉3 = 1.0001 − 𝑗0.0499
(6)
𝑉2 = 0.9801 − 𝑗0.0599
(6)
𝑉3 = 1.0000 − 𝑗0.0500
(7)
𝑉2 = 0.9800 − 𝑗0.0600
(7)
𝑉3 = 1.0000 − 𝑗0.0500
 Haris Hunafa Hanifan / 1610631160064

Sehingga iterasi yang diambil adalah pada iterasi ke-7

𝑉2 = 0.9800 − 𝑗0.0600 = 0.98183∠ − 3.5035˚ 𝑝𝑢

𝑉3 = 1.0000 − 𝑗0.0500 = 1.00125∠ − 2.8624˚ 𝑝𝑢

Sehingga dapat diperoleh nilai P-Q:

𝑃1 − 𝑗𝑄1 = 𝑉1∗ [𝑉1 (𝑦12 + 𝑦13) − (𝑦12 𝑉2 + 𝑦13 𝑉3 )]

= 1.05[1.05(20 − 𝑗50) − (10 − 𝑗20)(0.98 − 𝑗. 06) − (10 − 𝑗30)(1.0 − 𝑗0.05)]
= 4.095 − 𝑗1.890 𝑝𝑢

Sehingga nilai

P1 = 4.095 Q1 = 1.890.
𝑝𝑢 = 409.5 MW 𝑝𝑢 = 189 Mvar

Mencari nilai arus pada setiap jalur

𝐼12 = 𝑦12 (𝑉1 − 𝑉2 )
= (10 − 𝑗20)[(1.05 + 𝑗0) − (0.98 − 𝑗0.06)]
= 1.9 − 𝑗0.8

𝐼21 = −𝐼12
= −1.9 + 𝑗0.8

𝐼13 = 𝑦13 (𝑉1 − 𝑉3 )

= (10 − 𝑗30)[(1.05 + 𝑗0) − (1.0 − 𝑗0.05)]
= 2.0 − 𝑗1.0

𝐼31 = −𝐼13
= −2.0 + 𝑗1.0

𝐼23 = 𝑦23(𝑉2 − 𝑉3 )
= (16 − 𝑗32)[(0.98 + 𝑗0.06) − (1.0 − 𝑗0.05)]
= −0.64 + 𝑗0.48

𝐼32 = −𝐼23
= 0.64 − 𝑗0.48
 Haris Hunafa Hanifan / 1610631160064

Aliran daya:

∗
𝑆12 = 𝑉1𝐼12
= (1.05 + 𝑗0.0)(1.9 + 𝑗0.8)
= 1.995 + 𝑗0.84 𝑝𝑢
= 199.5 𝑀𝑊 + 𝑗84.0 𝑀𝑣𝑎𝑟

∗
𝑆21 = 𝑉2 𝐼21
= (0.98 − 𝑗0.06)(−1.9 − 𝑗0.8)
= −1.91 − 𝑗0.67 𝑝𝑢
= −191.0 𝑀𝑊 + 𝑗67.0 𝑀𝑣𝑎𝑟

∗
𝑆13 = 𝑉1𝐼13
= (1.05 + 𝑗0.0)(2.0 + 𝑗1.0)
= 2.1 + 𝑗1.05 𝑝𝑢
= 210.0 𝑀𝑊 + 𝑗105.0 𝑀𝑣𝑎𝑟

∗
𝑆31 = 𝑉3 𝐼31
= (1.0 + 𝑗0.05)(−2.0 − 𝑗1.0)
= −2.05 − 𝑗0.90 𝑝𝑢
= −205.0 𝑀𝑊 − 𝑗90.0 𝑀𝑣𝑎𝑟

∗
𝑆23 = 𝑉2 𝐼23
= (0.98 − 𝑗0.06)(−0.656 + 𝑗0.48)
= −0.656 − 𝑗0.432 𝑝𝑢
= −65.6 𝑀𝑊 + 𝑗43.2 𝑀𝑣𝑎𝑟

∗
𝑆32 = 𝑉3 𝐼32
= (1.0 − 𝑗0.05)(0.64 + 𝑗0.48)
= 0.664 + 𝑗0.448 𝑝𝑢
= 66.4 𝑀𝑊 + 𝑗44.8 𝑀𝑣𝑎𝑟

Sehingga losses daya yang dihasilkan adalah:

𝑺𝑳 𝟏𝟐 = 𝑺𝟏𝟐 + 𝑺𝟐𝟏 = 𝟖. 𝟓 𝑴𝑾 + 𝒋𝟏𝟕. 𝟎 𝑴𝒗𝒂𝒓
𝑺𝑳 𝟏𝟑 = 𝑺𝟏𝟑 + 𝑺𝟑𝟏 = 𝟓. 𝟎 𝑴𝑾 + 𝒋𝟏𝟓. 𝟎 𝑴𝒗𝒂𝒓
𝑺𝑳 𝟐𝟑 = 𝑺𝟐𝟑 + 𝑺𝟑𝟐 = 𝟎. 𝟖 𝑴𝑾 + 𝒋𝟏. 𝟔𝟎 𝑴𝒗𝒂𝒓
2
Haris Hunafa Hanifan / 1610631160064

PENYELESAIAN DENGAN MENGGUNAKAN METODE

NEWTON-RHAPSON

1 2
0.02 + 𝑗0.04
400
MW

250
0.01 + 𝑗0.03 0.0125 + 𝑗0.025 Mvar
Slack Bus 3
𝑉1 = 1.05∠0°
200 | 𝑉3 | = 1.04
MW

Konversikan saluran impedansi ke dalam admittance

1
𝑦12 = = 10 − 𝑗20
0.02 + 𝑗0.04

1
𝑦13 = = 10 − 𝑗30
0.01 + 𝑗0.03

1
𝑦23 = = 16 − 𝑗32
0.0125 + 𝑗0.025

Kemudian ubahlah kedalam bentuk matrix:

𝑌11 𝑌12 𝑌13
𝑌𝑏𝑢𝑠 = [𝑌21 𝑌22 𝑌23 ]
𝑌31 𝑌32 𝑌33

20 − 𝑗50 −10 + 𝑗20 −10 + 𝑗30

𝑌𝑏𝑢𝑠 = [−10 + 𝑗20 26 − 𝑗52 −16 + 𝑗32]
−10 + 𝑗30 −16 + 𝑗32 26 − 𝑗62
Konversikan matrix diatas ke dalam bentuk radian:
53.85165∠ − 1.9029 22.36068∠2.0344 31.62278∠1.8925
𝑌𝑏𝑢𝑠 = [ 22.36068∠2.0344 58.13777∠ − 1.1071 35.77709∠2.0344 ]
31.62278∠1.8925 35.77709∠2.0344 67.23095∠ − 1.1737
 Haris Hunafa Hanifan / 1610631160064

Mencari nilai daya aktif pada bus 2 dan 3 serta daya reaktif pada bus 2 menggunakan metode
power balance equation terhadap 𝛿2 , 𝛿3 , 𝑑𝑎𝑛 |𝑉2 |: asumsikan bahwa |V1|=1 dan 𝛿1 = 0, 𝛿3 =
0

𝑷𝟐 = |𝑽𝟐 ||𝑽𝟏 ||𝒀𝟐𝟏 | 𝐜𝐨𝐬(𝜽𝟐𝟏 − 𝜹𝟐 + 𝜹𝟏 ) + |𝑽𝟐𝟐 ||𝒀𝟐𝟐 |𝒄𝒐𝒔𝜽𝟐𝟐 + |𝑽𝟐 ||𝑽𝟑 ||𝒀𝟐𝟑 |𝐜𝐨𝐬(𝜽𝟐𝟑 − 𝜹𝟐 + 𝜹𝟑 )
𝜕𝑃2
= |𝑉2 ||𝑉1 ||𝑌21 | sin(𝜃21 − 𝛿2 + 𝛿1 ) + |𝑉2 ||𝑉3 ||𝑌23 |sin(𝜃23 − 𝛿2 + 𝛿3 )
𝜕𝛿2
= | − 10 + 𝑗20||𝑉2 | sin(0.9994 − 𝛿2 ) + | − 16 + 𝑗32||𝑉2 ||1.04|sin(0.9994 − 𝛿2 )
= | − 26.64 + 𝑗53.28|sin(0.9994 − 𝛿2 )

𝜕𝑃2
= −|𝑉2 ||𝑉3 ||𝑌23 | sin(𝜃23 − 𝛿2 + 𝛿3 )
𝜕𝛿3

= −| − 16 + 𝑗32||1.04||𝑉2 |sin(0.9994 − 𝛿2 )
= −| − 16.64 + j33.28||𝑉2 |sin(0.9994 − 𝛿2 )

𝜕𝑃2
= |𝑉1 ||𝑌21 | cos(𝜃21 − 𝛿2 + 𝛿1 ) + 2|𝑉2 ||𝑌22 | cos 𝜃22 + |𝑉3 ||𝑌23 | cos( 𝜃23 − 𝛿2 + 𝛿3 )
𝜕|𝑉2 |
= |−10 + 𝑗20| cos(2.0344 − 𝛿2 ) + 2|26 − 𝑗52||𝑉2 | cos(−1.1071) + |1.04|| − 16 + 𝑗32|cos(2.0344
− 𝛿2 )
= |−26.64 + 𝑗53.28| cos(2.0344 − 𝛿2 ) + |52 − 𝑗104||𝑉2 | cos(−1.1071)

𝑷𝟑 = |𝑽𝟑 ||𝑽𝟏 ||𝒀𝟑𝟏 | 𝐜𝐨𝐬(𝜽𝟑𝟏 − 𝜹𝟑 + 𝜹𝟏 ) + |𝑽𝟑 ||𝑽𝟐 ||𝒀𝟑𝟐 | 𝐜𝐨𝐬(𝜽𝟑𝟐 − 𝜹𝟑 + 𝜹𝟐 ) + |𝑽𝟐𝟑 ||𝒀𝟑𝟑 |𝒄𝒐𝒔𝜽𝟑𝟑
𝜕𝑃3
= −|𝑉3 ||𝑉2 ||𝑌32 | sin(𝜃32 − 𝛿3 + 𝛿2 )
𝜕𝛿2

= −|1.04|| − 16 + j32||V2|sin(0.9994 − 𝛿3 + 𝛿2 )
= −| − 16.64 + 𝑗33.28||V2|sin(0.9994 − 𝛿3 + 𝛿2 )

𝜕𝑃3
= |𝑉3 ||𝑉1 ||𝑌31 | sin(𝜃31 − 𝛿3 + 𝛿1 ) + |𝑉3 ||𝑉2 ||𝑌32 |sin(𝜃32 − 𝛿3 + 𝛿2 )
𝜕𝛿3
= | − 10.4 + j31.2| sin(0.99945 − 𝛿3 ) + | − 16.64 + 𝑗33.28||V2|sin(0.9994 − 𝛿3 + 𝛿2 )

𝜕𝑃3
= |𝑉3 ||𝑌32 | cos(𝜃32 − 𝛿3 + 𝛿2 )
𝜕|𝑉2 |
= | − 16.64 + 𝑗33.28|cos(2.0344 − 𝛿3 + 𝛿2 )

𝑸𝟐 = −|𝑽𝟐 ||𝑽𝟏 ||𝒀𝟐𝟏 | 𝐬𝐢𝐧(𝜽𝟐𝟏 − 𝜹𝟐 + 𝜹𝟏 ) − |𝑽𝟐𝟐 ||𝒀𝟐𝟐 |𝒔𝒊𝒏𝜽𝟐𝟐 − |𝑽𝟐 ||𝑽𝟑 ||𝒀𝟐𝟑 |𝐬𝐢𝐧(𝜽𝟐𝟑 − 𝜹𝟐 + 𝜹𝟑 )
𝜕𝑄2
= |𝑉2 ||𝑉1 ||𝑌21 | cos(𝜃21 − 𝛿2 + 𝛿1 )+|𝑉2 ||𝑉3 ||𝑌23 |cos(𝜃23 − 𝛿2 + 𝛿3 )
𝜕𝛿2

= |−10 + j20||𝑉2 | cos(0.9994 − 𝛿2 ) + | − 16.64 + 𝑗33.28||𝑉2 |cos(0.9994 − 𝛿2 + 𝛿3 )

 Haris Hunafa Hanifan / 1610631160064

𝜕𝑄2
= −|𝑉2 ||𝑉3 ||𝑌23 | cos(𝜃23 − 𝛿2 + 𝛿3 )
𝜕𝛿3

= −| − 16.64 + 𝑗33.28||𝑉2 |cos(0.9994 − 𝛿2 + 𝛿3 )

𝜕𝑄2
= −|𝑉1 ||𝑌21 | sin(𝜃21 − 𝛿2 + 𝛿1 ) − 2|𝑉2 ||𝑌22 | sin(𝜃22 ) − |𝑉3 ||𝑌23 |𝑠𝑖𝑛(𝜃23 − 𝛿2 + 𝛿3 )
𝜕|𝑉2 |
= −|−10 + j20| sin(2.0344 − 𝛿2 ) − |52 − j104||𝑉2 | sin(−1.1071) − | − 16.64 + j33.28|sin(2.0344
− 𝛿2 + 𝛿3 )

Ubah ke bentuk Jacobian matrix:

𝝏𝑷𝟐 𝝏𝑷𝟐 𝝏𝑷𝟐
𝝏𝜹𝟐 𝝏𝜹𝟑 𝝏|𝑽𝟐 |
𝝏𝑷𝟑 𝝏𝑷𝟑 𝝏𝑷𝟑 ∆𝜹𝟐 ∆𝑷𝟐
[ ∆𝜹𝟑 ] = [∆𝑷𝟑 ]
𝝏𝜹𝟐 𝝏𝜹𝟑 𝝏|𝑽𝟐 |
∆|𝑽𝟐 | ∆𝑸𝟐
𝝏𝑸𝟐 𝝏𝑸𝟐 𝝏𝑸𝟐
[ 𝝏𝜹𝟐 𝝏𝜹𝟑 𝝏|𝑽𝟐 | ]

−𝟐𝟔. 𝟔𝟒 + 𝒋𝟓𝟑. 𝟐𝟖|𝑽𝟐 |𝒔𝒊𝒏(𝟎. 𝟗𝟗𝟗𝟒 − 𝜹𝟐 ) 𝟏𝟔. 𝟔𝟒 − 𝐣𝟑𝟑. 𝟐𝟖||𝑽𝟐 |𝐬𝐢𝐧(𝟎. 𝟗𝟗𝟗𝟒 − 𝜹𝟐 ) |−𝟐𝟔. 𝟔𝟒 + 𝒋𝟓𝟑. 𝟐𝟖| 𝐜𝐨𝐬(𝟐. 𝟎𝟑𝟒𝟒 − 𝜹𝟐 ) + |𝟓𝟐 − 𝒋𝟏𝟎𝟒||𝑽𝟐 |𝐜𝐨𝐬(−𝟏. 𝟏𝟎𝟕𝟏) ∆𝜹𝟐 ∆𝑷𝟐
[ 𝟏𝟔. 𝟔𝟒 − 𝒋𝟑𝟑. 𝟐𝟖||𝑽𝟐 |𝐬𝐢𝐧(𝟎. 𝟗𝟗𝟗𝟒 − 𝜹𝟑 + 𝜹𝟐 ) −𝟏𝟎. 𝟒 + 𝐣𝟑𝟏. 𝟐| 𝐬𝐢𝐧(𝟎. 𝟗𝟗𝟗𝟒𝟓 − 𝜹𝟑 ) + −𝟏𝟔. 𝟔𝟒 + 𝒋𝟑𝟑. 𝟐𝟖||𝐕𝟐|𝐬𝐢𝐧(𝟎. 𝟗𝟗𝟗𝟒 − 𝜹𝟑 + 𝜹𝟐 ) −𝟏𝟔. 𝟔𝟒 + 𝒋𝟑𝟑. 𝟐𝟖|𝐜𝐨𝐬(𝟐. 𝟎𝟑𝟒𝟒 − 𝜹𝟑 + 𝜹𝟐 ] [ ∆𝜹𝟑 ] = [∆𝑷𝟑 ]
|−𝟏𝟎 + 𝐣𝟐𝟎||𝑽𝟐 |𝐜𝐨𝐬(𝟎. 𝟗𝟗𝟗𝟒 − 𝜹𝟐 ) + | − 𝟏𝟔. 𝟔𝟒 + 𝒋𝟑𝟑. 𝟐𝟖||𝑽𝟐|𝐜𝐨𝐬(𝟎. 𝟗𝟗𝟗𝟒 − 𝜹𝟐 + 𝜹𝟑 ) −| − 𝟏𝟔. 𝟔𝟒 + 𝒋𝟑𝟑. 𝟐𝟖||𝑽𝟐 |𝐜𝐨𝐬(𝟎. 𝟗𝟗𝟗𝟒 − 𝜹𝟐 + 𝜹𝟑 ) −|−𝟏𝟎 + 𝐣𝟐𝟎| 𝐬𝐢𝐧(𝟐. 𝟎𝟑𝟒𝟒 − 𝜹𝟐 ) − |𝟓𝟐 − 𝐣𝟏𝟎𝟒||𝑽𝟐 |𝐬𝐢𝐧(−𝟏. 𝟏𝟎𝟕𝟏) − | − 𝟏𝟔. 𝟔𝟒 + 𝐣𝟑𝟑. 𝟐𝟖|𝐬𝐢𝐧(𝟐. 𝟎𝟑𝟒𝟒 − 𝜹𝟐 + 𝜹𝟑 ) ∆|𝑽𝟐 | ∆𝑸𝟐

(400+𝑗250)
𝑆2𝑠𝑐ℎ = − = −4.0 − 𝑗2.5 𝑝𝑢
100
200
𝑃3𝑠𝑐ℎ = 100 = 2.0 𝑝𝑢
250
𝑄2𝑠𝑐ℎ = 100 = 2.5 𝑝𝑢

Sehingga didapatkan:
(0) (0)
∆𝑃2 = 𝑃2𝑠𝑐ℎ − 𝑃2 = −4.0 − (−1.14) = −2.8600
(0) (0)
∆𝑃3 = 𝑃3𝑠𝑐ℎ − 𝑃3 = 2.0 − (0.5616) = 1.4384
(0) (0)
∆𝑄2 = 𝑄2𝑠𝑐ℎ − 𝑄2 = −2.5 − (−2.28) = −0.2200
Matrik jacobian dengan menggunakan estimasi awal (iterasi pertama):
(0)
−2.8600 54.28000 −33.28000 24.86000 ∆𝛿2
[ 1.4384 ] = [−33.28000 66.04000 −16.64000] [ ∆𝛿3(0) ]
−0.2200 −27.14000 16.64000 49.72000 ∆|𝑉 (0)|
2
 Haris Hunafa Hanifan / 1610631160064

Sehingga didapatkan :
(0)
∆𝛿2 = −0.045263
(0)
∆𝛿3 = −0.00718
(0)
∆|𝑉2 | = −0.026548

Dengan mengasumsikan nilai |V1|=1 dan 𝜹𝟏 = 𝟎, 𝜹𝟑 = 𝟎

Maka
(1)
𝛿2 = 0 − 0.045263 = −0.045263
(1)
𝛿3 = 0 − 0.00718 = −0.00718
(1)
|𝑉2 | = 1 − 0.026548 = 0.97345

Sudut Phase tegangan dalam radian pada iterasi kedua:

(1)
−0.099218 51.724675 −31.765618 21.302567 ∆𝛿2
[ 0.021715 ] = [−32.981642 65.656383 −15.379086] [ ∆𝛿3(1) ]
−0.050914 −28.538577 17.402838 48.103589 ∆|𝑉 (1) |
2

Sehingga didapatkan:
(1)
∆𝛿2 = −0.001795
(1)
∆𝛿3 = −0.000985
(1)
∆|𝑉2 | = −0.001767

Maka
(2)
𝛿2 = −0.045263 − 0.001795 = −0.04706
(2)
𝛿3 = −0.00718 − 0.000985 = −0.008703
(2)
|𝑉2 | = 0.97345 − 0.001767 = 0.971684

Sehingga pada iterasi ketiga didapatkan:

(2)
−0.000216 51.596701 −31.693866 21.147447 ∆𝛿2
[ 0.000038 ] = [−32.933865 65.597585 −15.351628] [ ∆𝛿3(2) ]
−0.000143 −28.548205 17.396932 47.954870 ∆|𝑉 (2)|
2
 Haris Hunafa Hanifan / 1610631160064

Sehingga diperoleh:
(2)
∆𝛿2 = −0.0000038
(2)
∆𝛿3 = −0.0000024
(2)
∆|𝑉2 | = −0.0000044

Maka
(3)
𝛿2 = −0.04706 − 0.0000038 = −0.04706
(3)
𝛿3 = −0.008703 − 0.0000024 = −0.008705
(3)
|𝑉2 | = 0.971684 − 0.0000044 = 0.971684

Dalm solusi dari 3 iterasi terdapat ketidakcocokan daya sebesar 2.5x10 -4 dengan
V2=0.97168∠-2.696˚ dan V3=1.04∠ − 𝟎. 𝟒𝟗𝟖𝟖˚.
Mencari nilai daya aktif pada bus 3 serta daya reaktif pada bus 3 menggunakan metode
power balance equation:
𝑸𝟑 = −|𝑽𝟑 ||𝑽𝟏 ||𝒀𝟑𝟏 | 𝐬𝐢𝐧(𝜽𝟑𝟏 − 𝜹𝟑 + 𝜹𝟏 ) − |𝑽𝟑 ||𝑽𝟐 ||𝒀𝟑𝟐 | 𝐬𝐢𝐧(𝜽𝟑𝟐 − 𝜹𝟑 + 𝜹𝟐 ) − |𝑽𝟑 |𝟐 |𝒀𝟑𝟑 | 𝐬𝐢𝐧 𝜽𝟑𝟑

= −|𝑉3 ||−10 + 𝑗30| sin(1.8925 − 𝛿3 + 𝛿1 ) − |𝑉3 ||−16 + 𝑗32||𝑉2 | sin(2.0344 − 𝛿3 + 𝛿2 ) − |𝑉3 |2 |26
− 𝐽62|sin(−1.1737)

= −|𝑉3 ||−10 + 𝑗30|(0.95) cos(𝛿3 ) − |𝑉3 ||−16 + 𝑗32|(0.894) cos(𝛿3 ) − |𝑉3 |2 | − 23.98 + 𝐽57.18|

= −|𝑉3 ||−9.5 + 𝑗28.5| cos(𝛿3 ) − |𝑉3 ||−14.304 + 𝑗28.608| cos(𝛿3 ) − |𝑉3 |2 | − 23.98 + 𝐽57.18|

𝑷𝟏 = |𝑽𝟏 |𝟐 |𝒀𝟏𝟏 | 𝐜𝐨𝐬(𝜽𝟏𝟏 ) − |𝑽𝟏 ||𝑽𝟐 ||𝒀𝟏𝟐 | 𝐜𝐨𝐬(𝜽𝟏𝟐 − 𝜹𝟏 + 𝜹𝟐 ) − |𝑽𝟏 ||𝑽𝟑 ||𝒀𝟏𝟑 | 𝐜𝐨𝐬(𝜽𝟏𝟑 − 𝜹𝟏 + 𝜹𝟑 )

= |𝑉1 |2 |20 − 𝑗50| cos(−1.9029) − |𝑉1 ||−10 + 𝑗20| cos(2.0344 − 𝛿1 + 𝛿2 )

− |𝑉1 ||𝑉3 ||−10 + 𝑗30| cos(1.8925 − 𝛿1 + 𝛿3 )

= |𝑉1 |2 |−6.5 + 𝑗16.3| + |𝑉1 ||4.4 − 𝑗8.8| sin(𝛿1 ) + |𝑉1 ||𝑉3 ||3.16 − 𝑗9.46| sin(𝛿1 )

𝑸𝟏 = |𝑽𝟏 |𝟐 |𝒀𝟏𝟏 | 𝐬𝐢𝐧(𝜽𝟏𝟏 ) − |𝑽𝟏 ||𝑽𝟐 ||𝒀𝟏𝟐 | 𝐬𝐢𝐧(𝜽𝟏𝟐 − 𝜹𝟏 + 𝜹𝟐 ) − |𝑽𝟏 ||𝑽𝟑 ||𝒀𝟏𝟑 | 𝐬𝐢𝐧(𝜽𝟏𝟑 − 𝜹𝟏 + 𝜹𝟑 )

= |𝑉1 |2 |20 − 𝑗50| sin(−1.9029) − |𝑉1 ||−10 + 𝑗20| sin(2.0344 − 𝛿1 + 𝛿2 )

− |𝑉1 ||𝑉3 ||−10 + 𝑗30| sin(1.8925 − 𝛿1 + 𝛿3 )

= |𝑉1 |2 |−18.91 + 𝑗47.27| − |𝑉1 ||−8.944 + 𝑗17.89| cos(𝛿1 ) − |𝑉1 ||𝑉3 ||9.49 + 𝑗28.46| cos(𝛿1 )
 Haris Hunafa Hanifan / 1610631160064

Ubah ke bentuk jacobian matrix:

𝝏𝑸𝟑 𝝏𝑸𝟑 𝝏𝑸𝟑
𝝏𝜹𝟐 𝝏𝜹𝟑 𝝏|𝑽𝟐 |
∆𝜹𝟐 ∆𝑸𝟑
𝝏𝑷𝟏 𝝏𝑷𝟏 𝝏𝑷𝟏
[ ∆𝜹𝟑 ] = [∆𝑷𝟏 ]
𝝏𝜹𝟐 𝝏𝜹𝟑 𝝏|𝑽𝟐 | ∆𝑸𝟏
∆|𝑽𝟐 |
𝝏𝑸𝟏 𝝏𝑸𝟏 𝝏𝑸𝟏
[ 𝝏𝜹𝟐 𝝏𝜹𝟑 𝝏|𝑽𝟐 |]

−|𝑽𝟑 ||−𝟏𝟒. 𝟑𝟎𝟒 + 𝒋𝟐𝟖. 𝟔𝟎𝟖| 𝐜𝐨𝐬(𝜹𝟐 ) −|𝑽𝟑 ||−𝟗. 𝟓 + 𝒋𝟐𝟖. 𝟓| 𝐜𝐨𝐬(𝜹𝟑 ) − |𝑽𝟑 ||−𝟏𝟒. 𝟑𝟎𝟒 + 𝒋𝟐𝟖. 𝟔𝟎𝟖| 𝐜𝐨𝐬(𝜹𝟑 ) −|𝑽𝟑 ||−𝟏𝟔 + 𝒋𝟑𝟐||𝑽𝟐 | 𝐬𝐢𝐧(𝟐. 𝟎𝟑𝟒𝟒 − 𝜹𝟑 + 𝜹𝟐 ) ∆𝜹𝟐 ∆𝑸𝟑
[ |𝑽𝟏 ||𝟒. 𝟒 − 𝒋𝟖. 𝟖| 𝐬𝐢𝐧(𝜹𝟏 ) |𝑽𝟏 ||𝑽𝟑 ||𝟑. 𝟏𝟔 − 𝒋𝟗. 𝟒𝟔| 𝐬𝐢𝐧(𝜹𝟏 ) |𝑽𝟏 ||−𝟏𝟎 + 𝒋𝟐𝟎| 𝐜𝐨𝐬(𝟐. 𝟎𝟑𝟒𝟒 − 𝜹𝟏 ) ] [ ∆𝜹𝟑 ] = [∆𝑷𝟏 ]
−|𝑽𝟏 ||−𝟖. 𝟗𝟒𝟒 + 𝒋𝟏𝟕. 𝟖𝟗| 𝐜𝐨𝐬(𝜹𝟏 ) −|𝑽𝟏 ||𝑽𝟑 ||𝟗. 𝟒𝟗 + 𝒋𝟐𝟖. 𝟒𝟔| 𝐜𝐨𝐬(𝜹𝟏 ) −|𝑽𝟏 ||−𝟏𝟎 + 𝒋𝟐𝟎| 𝐬𝐢𝐧(𝟐. 𝟎𝟑𝟒𝟒 − 𝜹𝟏 ) ∆|𝑽𝟐 | ∆𝑸𝟏

Jadi, nilai P-Q bus 3 yang didapatkan adalah:

Q = 1.4617 pu
P1 = 2.1842 pu
Q1 = 1.4085 pu