J. Math. and Its Appl. ISSN: 1829-605X Vol. 4, No. 1, May 2007, 27–37

## Analisis dan Perancangan Sistem Pengendali pada Sistem Pembangkit Listrik Tenaga Uap dengan Pengendali Robust Melalui Optimasi H ∞

## Mardlijah dan Aris Dian Permata
Jurusan Matematika FMIPA ITS Surabaya tidak ada email

## Abstrak
Naik turunnya kebutuhan konsumen terhadap energi listrik berpengaruh pada frekuensi sistem pembangkit listrik. Hal ini menyebabkan sistem pe- ngendali pada sistem pembangkit tersebut harus mengatur kembali parame- ter pengendalinya guna mempertahankan kestabilan sistem. Sehingga untuk lebih memperbaiki performansi sistem pembangkit listrik digunakan sistem pengendali optimal.
Pada penelitian ini dirancang sebuah sistem pengendali yang robust melalui optimasi yang diimplementasikan pada sistem pembangkit listrik tenaga uap dengan reheat steam turbine . Dari hasil simulasi yang dilakukan, overshoot sistem pengendali yang dirancang lebih kecil hingga 29 , 1667% dibandingkan sistem dengan gangguan dan waktu yang diperlukan untuk mencapai titik kestabilan lebih cepat hingga 21 , 3043%. Sistem pengendali yang dirancang memiliki sensitivitas yang lebih kecil hingga 30 , 0678%. Se- hingga sistem pengendali yang dirancang dapat dikatakan memiliki kestabil- an dan performansi yang robust melalui optimasi H ∞ .
Kata Kunci: Sistem PLTU dengan Reheat Steam Turbine, Pengendali Ro- bust, Optimasi H ∞

## 1. Pendahuluan
Energi listrik merupakan salah satu kebutuhan dasar manusia. Energi tersebut merupakan hasil dari sistem pembangkit listrik yang didistribusikan pusat pe- ngatur beban melalui jaringan transmisi. Perubahan jumlah energi listrik yang dibutuhkan konsumen berpengaruh terhadap frekuensi sistem pembangkit listrik. Oleh karena itu, sistem pengendali pada sistem pembangkit tersebut harus menga- tur kembali parameter pengendalinya guna mempertahankan kestabilan sistem.
Analisis terhadap kestabilan sistem pembangkit listrik dapat digunakan untuk mengetahui tingkat kestabilan sistem tersebut dan menjaga kelangsungan pro- ses operasi. Kestabilan sistem pembangkit listrik tenaga uap yang dikaji dengan memperhatikan parameter-parameter pada governor, turbin dan generator.
Sistem pengendali harus dirancang sedemikian hingga robust , baik terhadap perubahan gangguan yang timbul, maupun ketidakpastian dari sistem yang di- kendalikan. Ketidakpastian ( uncertainties ) dalam sistem nyata yang diubah ke dalam model matematika disebabkan oleh asumsi yang digunakan. Suatu sistem dapat dikatakan robust yaitu jika sistem tersebut mampu bekerja dengan baik dan mempertahankan kestabilan walaupun terdapat kesalahan pemodelan maupun gangguan yang tidak diketahui karakteristiknya.
Dalam penelitian ini akan dirancang suatu sistem pengendali robust melalui optimasi H ∞ sebagai metode alternatif memperbaiki performansi sistem pem- bangkit listrik tenaga uap dengan reheat steam turbine . Pengendali umpan-balik ini dirancang untuk meminimumkan overshoot dan waktu yang dibutuhkan untuk mencapai kestabilan. Manfaat dari penelitian ini adalah untuk memberikan gam- baran bagaimana merancang suatu sistem pengendali pada sistem pembangkit listrik tenaga uap agar keluaran sesuai dengan yang diinginkan.
Permasalahan yang dibahas disini dibatasi pada:
a. Masalah pengendali daya reaktif (generator) dari sistem pembangkit listrik tenaga uap dengan reheat steam turbine.
b. Variabel terkendali berupa frekuensi generator dan variabel termanipulasi berupa masukan ke generator.
c. Sinyal masukan, kondisi keadaan dan sinyal keluaran kontinu terhadap waktu.
d. Pemodelan sistem yang digunakan berdasarkan referensi.
e. Sistem bersifat linier, berdimensi hingga, tidak berubah terhadap waktu (steady-state), non-autonomous dan real.
f. Sistem berupa Single Input Single Output (SISO).
g. Galat pada sistem nominal dan sistem pengendali terbatas dan tak terstruktur.
h. Rancangan pengendali yang dibuat hanya bisa dijalankan pada program Mat- lab versi 7.0.
Dengan asumsi bahwa:
a. Gangguan yang dapat diatasi merupakan gangguan yang tidak bisa diprediksi (unpredictable) besarnya dan tidak diketahui karakteristiknya serta terbatas.
b. Keluaran tenaga turbin sama dengan keluaran generator.
Sedangkan metode yang dilakukan dalam penelitian ini antara lain:
• Mempelajari model matematika dari sistem pembangkit listrik tenaga uap.
• Melakukan simulasi dengan memasukkan nilai parameter pada simulink pro- gram Matlab 7.0 sehingga memperoleh bentuk ruang keadaannya.
• Menguji keminimalan sistem.
• Menentukan fungsi bobot (weighting function) yang sesuai.
• Menghitung pengendali robust melalui optimasi H ∞ .
• Menganalisis hasil dari simulasi yang dilakukan.

## 2. Pemodelan Sistem Pembangkit Listrik Tenaga Uap dengan Reheat Steam Turbine

## Gambar 1: EHC
Governor yang digunakan adalah Electro-Hydraulic Governor Control System ( EHC ). Bagian electric merupakan pemicu gerakan hidrolik ( hydraulic amplifier ).
Gerakan titik D dinyatakan sebagai:
∆ x D = k 6 Z x B 1 dt (1)
Sedangkan gerakan di titik A dinyatakan:
∆ x A = − k 1 ∆ ω (2)
Unit pembangkit dari governor diparalel dengan sistem yang besar, maka x B 2 juga dapat menyebabkan perubahan daya ∆ p dan tidak menimbulkan perubahan frekuensi, sehingga diperoleh:
∆ x B 2 = k 2 ∆ p (3)
Gerakan di titik B 1 dipengaruhi oleh titik A , maka:
∆ x B 1 = − k 3 k 1 ∆ ω + k 4 k 2 ∆ p (4)
Karena ada umpan balik dari titik D ke titik B 1 , maka:
∆ x B 1 = − k 3 k 1 ∆ ω + k 4 k 2 ∆ p + k
5 ∆ x D (5)
Dengan menggunakan transformasi Laplace dari bidang t ke bidang s , sehingga: ∆ ω menjadi ∆ F ( s ), ∆ p menjadi ∆ P ( s ) ∆ x C menjadi ∆ x C ( s ), ∆ x D menjadi ∆ x D ( s ) dan
∆ x B 1 ( s ) = − k 3 k 1 ∆ F ( s ) + k 4 k 2 ∆ P ( s ) + k 5 ∆ x D ( s ) (6)
∆ x D ( s ) = − k 6 s ∆ x b 1 ( s ) (7)
Maka dari persamaan (6) dan (7) diperoleh:
∆ x D ( s ) = − k 4 k 2 k 5 1 + s k 5 k 6 · ∆ P ( s ) − k 3 k 1 k 4 k 2 ∆ F ( s ) ¸ (8)
Jika: R = − k 4 k 2 k 3 k 1 , K G = k 4 k 2 k 5 , T G = 1 k 5 k 6 , G G ( s ) = K G 1+ sT G maka Persamaan (8) menjadi:
∆ x D ( s ) = G G ( s ) · ∆ P ( s ) − 1 R ∆ F ( s ) ¸ (9)
Uap keluaran dari governor (∆ x D ( s )) akan masuk ke dalam turbin, maka di- peroleh hubungan antara ∆ x D ( s ) dengan keluaran dari generator sebagai berikut:
G T ( s ) = ∆ P G ( s ) ∆ x D ( s ) = K T (1 + sF HP T RH (1 + sT CH )(1 + sT RH ) (10)
Saat terjadi penambahan beban sebesar ∆ P D , governor akan membuka lebih lebar. Penambahan beban ini akan menyebabkan penambahan daya sebesar ∆ P G , sehingga tambahan daya bersih yang diterima sistem adalah ∆ P G − ∆ P D , dengan:
∆ P G − ∆ P D = d dt W k + D ∆ ω (11)
Energi kinetik sebanding dengan kuadrat kecepatan sudut yang juga berarti meru- pakan kuadrat dari frekuensi, sehingga:
W k = ³ ω ω ∗ ´ 2 W ∗ k = µ ω ∗ + ∆ ω ω ∗ ¶ 2 W ∗ k (12)
Persamaan (12) menjadi:
W k = µ 1 + ∆ ω ω ∗ ¶ 2 W ∗ k ≈ µ 1 + 2 ∆ ω ω ∗ ¶ W ∗ k
karena kecil, sehingga diperoleh perubahan energi kinetik yaitu:
d dt W k = 2 W ∗ k ω ∗ d dt (∆ ω ) (dalam MW) (13)
Di lapangan, besaran energi kinetik sering dinyatakan dalam:
H = W k P S = ⇒ W k = HP S
Dengan substitusi persamaan diatas dan Persamaan (13) ke dalam Persamaan (11), maka diperoleh:
∆ P G − ∆ P D = 2 HP S ω ∗ d dt (∆ f ) + D ∆ ω = 2 H ω ∗ d dt (∆ ω ) + D ∆ ω (14)
Persamaan (14) di transformasi Laplace menjadi:
∆ P G ( s ) − ∆ P D ( s ) = µ 2 H ω ∗ s + D ¶ ∆ F ( s )
dengan
∆ F ( s ) = µ ω ∗ 2 Hs + Dω ∗ ¶ (∆ P G ( s ) − ∆ P D ( s )) (15)
Sehingga Persamaan (15) menjadi:
∆ F ( s ) ∆ P G ( s ) − ∆ P D ( s ) = 1 D + M s (16)
Dari uraian diatas diperoleh blok diagram untuk sistem PLTU dengan reheat steam turbine sebagai berikut:
dengan:
+ _ + _ 1 1+sT G 1 Da+Ms 1+sF T (1+sT )(1+sT ) HP RH CH RH 1 R Turbin dengan reheat Rotor inertia & Load (generator) Governor D P(s) D F(s) D x (s) D D P (s) G D P (s) D Dw
Gambar 2: Simulink Sistem PLTU dengan Reheat Steam Turbine
∆ P : perubahan daya ∆ x P : setting steam valve
∆ P G : perubahan keluaran generator ∆ P D : penambahan beban ∆ ω : perubahan frekuensi generator F HP : koefisien pemanas ulang (reheat)
R : konstanta pengaruh kerja governor terhadap frekuensi (speed droop) T G : time-constant dari mekanisme speed governing (detik) T CH : time-constant turbin bertekanan rendah (low pressure turbine) (detik) T RH : time-constant turbin bertekanan tinggi (high pressure turbine) (detik) H : konstanta inersia generator Da : konstanta load-damping

## 3. Analisis dan Pembahsan
Gambar 3: Simulink Sistem PLTU dengan Reheat Steam Turbine
Gambar 3 menunjukan model sistem PLTU dengan reheat turbine dalam ben- tuk simulink Matlab dengan:
T G = 0 , 09 det, R = 0 , 034 , M = 10 det, Da = 1 , 389 ,
F HP = 0 , 3 , T CH = 0 , 3 det, T RH = 7 det
Dari simulink diatas, diperoleh matriks ruang keadaan:
x = Ax + Bu y = Cx + Du
dengan:
A =     − 0 , 1389 1 0 , 1 0 0 − 3 , 3333 0 11 , 1111 0 3 , 3333 − 0 , 1429 0 − 2 , 9412 0 0 − 11 , 1111

   , B =     − 1 0 0 0     , C = ¡ 0 , 1 0 0 0 ¢ , D = ¡ 0 ¢
Dan matriks transfer G yaitu:
G =         
− 0 , 1389 1 0 , 1 0 − 1 4 0 − 3 , 3333 0 11 , 1111 0 0 0 3 , 3333 − 0 , 1429 0 0 0 − 2 , 9412 0 0 − 11 , 1111 0 0 0 , 1 0 0 0 0 0 0 0 0 0 0 ∞         
Sistem pembangkit yang diteliti bersifat minimal karena rank dari matriks keterkendalian dan keteramatannya sama dengan ordo matriks A yaitu n = 4. Kemudian dipilih fungsi bobot sedemikian hingga magnitudonya diatas sistem nominal dan menyebabkan matriks D 12 memiliki rank yang penuh, dengan:
w 1 = 9 0 , 1 s + 1 , w 2 = 0 , 17 1 = 0 , 17 w 3 = 0
Gangguan yang dipilih merupakan salah satu parameter dari sub blok pengen- dali daya reaktif (generator) sistem PLTU dengan reheat steam turbine, yaitu kon- stanta load-damping ( Da ) dengan asumsi bersifat tidak diketahui karateristiknya dan tidak dapat diprediksi besarnya. Sistem dengan gangguan akan memasukkan beberapa nilai Da , diantaranya Da = 1.
Dari perhitungan yang dilakukan dengan Matlab diperoleh gamma optimal sebesar 0 , 1543 pada iterasi ke delapan dengan toleransi kesalahan 0 , 01. Over- shoot dari sistem pengendali robust lebih kecil dibandingkan dengan sistem dengan gangguan. Perbandingan tingkat overshoot ini berkisar antara 27 , 1186% sampai
Gambar 4: Grafik Fungsi Bobot untuk Da plant with disturbance = 1
Gambar 5: Step response sistem nominal
dengan 29 , 1667%. Begitu pula waktu yang diperlukan sistem pengendali robust untuk mencapai titik kestabilan lebih cepat hingga 21 , 3043%.
Berdasarkan Gambar 5, overshoot sistem nominal adalah sebesar 0 , 0531 dan waktu untuk mencapai titik kestabilan sistem nominal tersebut adalah pada saat t = 9 , 41 det .
Overshoot sistem dengan gangguan lebih besar 0 , 0021 dibandingkan dengan sistem nominal dan waktu yang diperlukan untuk mencapai kestabilan lebih lama 2 , 9 detik.
Dari Gambar 7 dapat dilihat bahwa overshoot dari sistem pengendali cen- derung lebih kecil atau waktu yang diperlukan sistem pengendali robust untuk mencapai titik kestabilan lebih cepat. Overshoot sistem pengendali untuk nilai gangguan Da = 1 adalah sebesar 0 , 0391 dan waktu yang diperlukan sistem pen- gendali untuk mencapai kestabilan adalah gt = 9 , 05 sec .
Safonov, Stein, Doyle dan Sandell menyatakan bahwa dasar utama pengen- dalian robust adalah jika nilai singular sistem pengendali yang dirancang memiliki nilai singular jauh dari titik kritis atau jika nilai singular sistem pengendali yang
Gambar 6: Step Response Sistem dengan Gangguan Da = 1
Gambar 7: Step Response untuk Da plant with disturbance = 1
dirancang jauh dari titik kritis, maka sensitivitas sistem tersebut kecil.
Berdasarkan Gambar 8, nilai singular sistem pengendali untuk nilai gang- guan Da = 10 , 02784 pada frekuensi 0 , 1 rad/sec dan 0 , 00852 pada frekuensi 10 rad/sec . Nilai singular maksimumnya adalah sebesar 0 , 02462 pada frekuen- si 0 , 3784 rad/sec . Sedangkan nilai singular sistem dengan gangguannya adalah 0 , 03981 pada frekuensi 0 , 1 rad/sec dan 0 , 01019 pada frekuensi 10 rad/sec . Nilai singular maksimumnya adalah sebesar 0 , 1123 pada frekuensi 0 , 8214 rad/sec . Hal ini menyatakan bahwa sistem pengendali yang dirancang memiliki nilai singular dibawah titik kritis. Sehingga, sistem pengendali yang dirancang memiliki sen- stivitas yang kecil. Selain itu, ukuran kesalahan pemodelan sistem pengendali yang dirancang juga kurang dari satu, sehingga sistem pengendali yang dirancang memenuhi kestabilan robust . Dengan kata lain, sistem pengendali yang dirancang bersifat robust dengan optimasi H ∞ .
Gambar 8: Nilai singular sistem untuk Da plant with disturbance = 1

## 4. Kesimpulan dan Saran

## 4.1. Kesimpulan
a. Dari hasil simulasi dan analisa yang dilakukan telah diperoleh sistem pengen- dali yang memiliki kestabilan dan performansi yang lebih baik. Overshoot dari sistem pengendali robust lebih kecil hingga 29 , 1667% dibandingkan sis- tem dengan gangguan. Waktu yang diperlukan sistem pengendali robust untuk mencapai titik kestabilan lebih cepat hingga 21 , 3043%.
b. Sensitivitas sistem pengendali yang dirancang terhadap gangguan lebih kecil hingga 30 , 0678% dibandingkan sistem dengan gangguan.
c. Karena sistem pengendali yang dirancang memiliki performansi dan kestabilan yang robust , maka sistem tersebut dikatakan sistem pengedali robust dengan optimasi H ∞ .

## 4.2. Saran
a. Untuk mencegah terjadinya error pada pemrograman dengan H inf ty optimal, dapat dipilih fungsi bobot (weighting function) dengan (I adalah matriks iden- titas ordo satu).
b. Penelitian ini dapat dilanjutkan dengan menggunakan gangguan yang lain sesuai parameter dari sistem pembangkit listrik dengan reheat steam turbine.

## Pustaka
[1] Green, Michael, and David J. N. Limebeer, 1996, Linear robust Control, Pren- tice Hall International Inc., USA.
[2] Kreyzig, Erwin, 1978, Introductory Functional Analysis with Application, John Wiley and Sons, Inc., Canada.
[3] Kundur, Prabha S., M. Zwyno, M. Klieng, GJ. Rongers, 1989, Application of Power System Stabilizers for Enhancement of Overall System Stability, Mc Graw-Hill Inc., USA.
[4] Ogata, Katsuhito, 1984, Modern Control Engineering, Prentice Hall, USA.
[5] Olsder, G. J., prof. dr. and dr. J. W. van der Woude, 1994, Mathematical System Theory, delftse Uitgevers Maatschappij b.v., Netherlands.
[6] Permata, Aris Dian, 2007, Pemodelan Matematika dan Analisis Sifat-Sifat Sistem Pembangkit Listrik di PT Pembangkitan Jawa-Bali (PT PJB) Unit Pembangkitan Gresik, Kerja Praktek, Matematika, ITS, Surabaya.
[7] Rijanto, Estiko, 2000, robust Control Theory for Application, ITB Press, Ban- dung.
[8] Setiabudi, Andri, 1998, Perancangan Sistem Pengendali Moisture Pada Mesin Kertas Dengan Metode Kontrol robust , Tugas Akhir, Teknik Fisika, ITS, Surabaya.
[9] Tamaela, Gysber Jan, 1993, Analisis Dan Perancangan Sistem Pengendalian Tekanan Pada Deaerator Secara robust Melalui Optimasi , Tugas Akhir, Matematika, ITS, Surabaya.
[10] Zhou, Kemin, and John C. Doyle, 1997, Essentials of robust Control, Prentice Hall, USA.
[11] NN, Robust Control Toolbox, Matlab 7.0, The Mathworks Inc.
Analisis dan Perancangan Sistem Pengendali ...
.