Analisis Radius Aman Akibat Kegagalan. (Sulistyo Atmadi et al. ANALISIS RADIUS AMAN AKIBAT KEGAGALAN STRUKTUR SUDU SKEA 50 KW PADA SAAT BEROPERASI Sulistyo Atmadi*). Ahmad Jamaludin Fitroh**) *) Peneliti Bidang Konversi Energi Dirgantara. LAPAN **) Peneliti Aerodinamika. LAPAN e-mail: sulistyoa@hotmail. ABSTRACT Analysis on the safe radius of a blade structural failure during operation, is studied on this paper. The particular blade under analysis is the SKEA 50 kW LAPAN. The rotor is installed at 18 m height. The analysis is based on the modified straight motion and rotational momentum formulation. The result showed that at the maximum design rotational speed of 80 rpm, the blade would be thrown out of it centre to 73 m horizontally, and 35 m vertically above the ground. On the other hand, the failure at the root of the blade would be thrown and spinned at its centre of mass at the speed of 145 rpm. Key words: Radius of safety. Rotating blade, 50 kW WECS ABSTRAK Telah dilakukan analisa radius aman akibat kegagalan struktur sudu atau sudu patah pada saat beroperasi. Analisis dilakukan pada sudu SKEA 50 kW LAPAN. Rotor dipasang pada menara setinggi 18 m. Dengan menggunakan Modifikasi dari Gerak Lurus Berubah Beraturan dan Momentum Putar, hasil analisa menunjukkan bahwa jika sudu mengalami kegagalan struktur atau patah pada putaran rancangan atau maksimum 80 rpm, maka sudu akan terlempar maksimum sejauh 73 m ke samping dan 35 m ke atas dihitung dari permukaan tanah. Akibat patah di pangkal, sudu akan terlempar dan berputar pada titik beratnya dengan putaran 145 rpm. Kata kunci: Radius aman. Sudu berputar. SKEA 50 kW PENDAHULUAN Turbin angin yang segera akan dipasang oleh LAPAN adalah SKEA 50 Sesuai dengan namanya, maka kapasitas daya maksimum 50 kW (Ahmad J. Fitroh, 2009. Sulistyo Atmadi & Ahmad J. Fitroh, 2. Faktor keamanan adalah unsur utama dalam pengoperasian Turbin Angin. Faktor keamanan yang dimaksud adalah apabila terjadi kegagalan operasi turbin angin, maka pada radius daerah tertentu dari turbin angin dinyatakan sebagai daerah yang harus dihindari pada saat pengoperasian turbin angin. Salah satu komponen utama SKEA 50 kW yang dapat menjadi penyebab kecelakaan adalah rotor dengan tiga buah sudu. Penelitian teoritis menunjukkan bahwa dengan sudu rotor yang modern akan terjadi kecepatan aerodinamika yang kritis kirakira tiga kali kecepatan ujung rancangan. Pada kecepatan angin yang konstan, ini berarti tiga kali kecepatan rotor (Milborrow,D. Apabila pada saat pengoperasian terjadi kegagalan struktur sudu, maka sudu akan terlempar. Dengan berat sekitar 230 kg per sudu, maka dampak yang akan diakibatkan apabila tertimpa sudu yang terlepas atau patah tersebut pasti sangat besar. Oleh karena itu perlu diperhitungan daerah atau radius aman selama pengoperasian SKEA 50 kW sehingga Jurnal Teknologi Dirgantara Vol. 8 No. 1 Juni 2010:57-63 apabila terjadi sudu lepas atau patah, maka dalam radius tersebut tidak ada korban jiwa maupun harta. 1 Tujuan Perhitungan dilakukan untuk mendapatkan radius aman selama pengoperasian SKEA 50 kW. Selain itu mendapatkan karakteristik sentrifugal failure sudu. Radius aman diperlukan agar tidak terjadi korban jiwa dan harta selama pengoperasian. Karakteristik mengetahui berapa kecepatan putar sudu terhadap titik beratnya pada saat 2 Batasan Masalah SKEA 50 kW yang sedang dikembangkan oleh LAPAN tersebut dirancang agar mampu menghasilkan daya sebesar 50 kW pada kecepatan angin 10 m/det dan putaran 80 rpm (Ahmad J. Fitroh, 2009. Sulistyo Atmadi & Ahmad J. Fitroh, 2. Dengan kata lain generator turbin angin yang maksimum 80 rpm. Oleh karena itu dalam perhitungan ini, putaran rotor yang dipilih adalah 80 rpm. Setiap sudu mempunyai berat 230 kg. Sudu tersebut diukur mempunyai titik berat pada posisi 2,407 m dihitung dari pangkal sudu. Apabila memperhitungkan jarak antara pangkal sudu dan sumbu putar rotor, maka titik berat sudu menjadi 2,907 m. Rotor dipasang di atas menara setinggi 18 m. Dengan demikian titik ketinggian 18 m di atas permukaan Pada saat sudu terlempar, terjadi gesekan antara permukaan sudu dengan Analisis perhitungan panjang lintasan dengan mempertimbangkan gaya angkat dan hambat aerodinamika menghasilkan bahwa suatu koefisien hambat udara rata-rata cD = 0,25 (Sorensen. Selain itu terjadi gaya hambat lainnya yang bekerja pada sudu sehingga dalam penelitian ini koefisien gaya hambat sudu, cD diasumsikan sebesar 0,6. Asumsi tersebut dipilih karena pada gerakan yang terdiri banyak manuver, harga cD menjadi cukup besar (Anderson. John D. Jr. Mc. Cormick. Barnes W. , 1. Penampang Untuk menghitung distribusi massa sudu, maka distribusi luas airfoil perlu dihitung terlebih dahulu. Untuk perhitungan luas airfoil diasumsikan sama dengan perhitungan luas belah ketupat, yaitu hanya memperhitungkan chord dan tebal airfoil. DASAR TEORI DAN METODE Dalam perhitungan radius aman akibat kegagalan rotor, titik berat sudu disajikan sebagai acuan. Sketsa turbin angin beserta sudunya disajikan dalam Gambar 2-1 dan 2-2 sebagai berikut: Gambar 2-1: Sketsa SKEA 50 kW tampak Analisis Radius Aman Akibat Kegagalan. (Sulistyo Atmadi et al. Pada saat terlempar, sudu akan mendapatkan gaya hambat udara. Besar gaya hambat tersebut secara langsung dipengaruhi oleh koefisien gaya hambat itu sendiri dan kecepatan, seperti pada persamaan . di bawah ini (Anderson. John D. Jr, 1985. Mc. Cormick. Barnes , 1. dA Gambar 2-2: Sketsa vektor pada sudu Gambar 2-2 menjelaskan vektor kecepatan pada sudu untuk tiga posisi, yaitu posisi 0, 120, dan 240 derajat. V tan merupakan kecepatan tangensial yang mana merupakan fungsi dari putaran dan posisi titik berat terhadap sumbu putar rotor, sebagaimana dituliskan sebagai berikut: Vtan A A R dengan A dan A masing Ae masing adalah kerapatan udara dan luas sudu. Kerapatan udara disesuaikan dengan kondisi pengujian. Luas sudu yang digunakan adalah luas permukaan Dalam pengerjaan ini luas AA Ac r A ct Ab cr = chord di pangkal . ct = chord di ujung . b = panjang sudu O = putaran . ad/de. R = posisi titik berat . Untuk menentukan jangkauan perpindahan sudu dalam arah horisontal dan vertikal akibat patah pada saat operasional, maka kecepatan tangensial tersebut juga harus diuraikan dalam arah horisontal dan vertikal sebagai Vx A Vtan * cosA V y A V tan * sin A A V 2 A cD . Vx = kecepatan tangensial dalam arah horisontal . /de. Vy = kecepatan tangensial dalam arah vertikal . /de. = posisi sudu . Notasi V dalam persamaan . merupakan kecepatan translasi sudu pada saat terlempar. Harga V diperoleh dari persamaan . Selain dalam bentuk satuan sudut, posisi sudu pada saat patah juga dapat disajikan dalam bentuk ketinggian. Rotor dipasang pada menara setinggi 18 m sehingga posisi sudu yang sebenarnya merupakan penjumlahan dari ketinggian menara dan posisi sudu pada saat Posisi dituliskan sebagai berikut: h A hm A R * cos A . h = ketinggian sudu dihitung permukaan tanah . hm = tinggi menara . Metode yang digunakan untuk menghitung ketinggian dan jarak terlempar Jurnal Teknologi Dirgantara Vol. 8 No. 1 Juni 2010:57-63 maksimum sudu setelah patah adalah Gerak Lurus Berubah Beraturan (GLBB). Persamaan tersebut diaplikasikan dalam arah vertikal dan mendatar sebagai berikut: 1 E mg A D y E 2 EE t y A h A V y t A g EE 1 Dx 2 x A Vx t A . = jarak antara elemen ke c. Harga x pada persamaan . di atas dapat berharga positif maupun negatif, tergantung dari arah referensinya. Momentum putar sudu diperoleh dari penjumlahan momentum putar per elemennya sebagai berikut: M A dM y = ketinggian maksimum . x = perpindahan ke samping maksimum . t = waktu yang diperlukan untuk mencapai ketinggian maksimum . g = percepatan gravitasi bumi . /det. Dy = gaya hambat dalam arah vertikal (N) Dx = gaya hambat dalam arah mendatar (N) Harga y dan x di atas berlaku untuk setiap harga pada saat sudu patah. Dengan memvariasikan harga , maka akan diperoleh radius maksimum untuk daerah aman pada saat pengujian sudu. Akibat adanya distribusi kecepatan tangensial di sepanjang sudu sesaat setelah sudu patah pada pangkalnya, akan membentuk momentum putar pada sudu. Adanya momentum putar tersebut menyebabkan sudu berputar terhadap titik beratnya sendiri pada saat Sketsa distribusi kecepatan tangensial tersebut disajikan dalam Gambar 2-3. Momentum putar per elemen sudu dipengaruhi oleh distribusi massa, kecepatan tangensial, dan posisi elemen terhadap titik berat sudu. Perumusan momentum putar per elemen sudu dapat dituliskan sebagai berikut: dM A dm * Vtan * x Gambar 2-3: Sketsa distribusi kecepatan HASIL ANALISA DAN PEMBAHASAN 1 Radius Aman Dalam analisa ini posisi sudut sudu divariasikan dari 0 hingga 330 derajat dengan kelipatan 15 derajat. Hasil mendapatkan radius aman disajikan dalam Gambar 3-1 berikut ini: Sudu SKEA 50 kW . 80 rpm . /de. Posisi . dM = momentum putar per elemen dm = massa per elemen Gambar 3-1a: Hasil perhitungan kecepatan awal terlempar Analisis Radius Aman Akibat Kegagalan. (Sulistyo Atmadi et al. Sudu SKEA 50 kW . 80 rpm (N) Posisi . Gambar 3-1b: Hasil perhitungan Sudu SKEA 50 kW . 80 rpm . Posisi . Gambar 3-1c: Jarak lontar Gambar 3-1d: Trayektori sudu terlempar Gambar 3-1a menyajikan hasil perhitungan kecepatan awal sudu Harga positif dan negatif untuk Vx hanya merupakan perbedaan gerak ke kanan dan kiri. Demikian juga pada harga Vy, harga positif dan negatif hanya menunjukkan perbedaan gerak ke atas dan ke bawah pada saat pangkal sudu patah. Dari hasil tersebut dapat dilihat bahwa misalnya sudu patah saat berada pada posisi 0 derajat yaitu tepat pada posisi vertikal ke atas, maka sudu akan mempunyai kecepatan awal terlempar ke arah kanan atau mendatar sebesar 24 m/det. Gambar 3-1b menyajikan hasil perhitungan gaya hambat sudu pada saat terlempar. Dari hasil tersebut dapat dilihat bahwa misalnya sudu patah saat berada pada posisi 90 derajat yaitu tepat pada posisi mendatar ke kanan, maka sudu hanya mendapatkan gaya hambat dalam arah vertikal saja sebesar 206 N. Gambar 3-1c menyajikan hasil perhitungan jarak lontar sudu pada saat Dari hasil tersebut dapat dilihat bahwa misalnya sudu patah saat berada pada posisi 180 derajat yaitu tepat pada posisi vertikal ke bawah, maka sudu akan terlontar ke kiri sejauh 41 m. Harga 15 m pada kolom y . rah vertika. bukan berarti sudu terlontar sejauh 15 m ke atas, namun angka tersebut merupakan posisi sudu terhadap tanah mengingat rotor dipasang pada menara setinggi 18 m. Posisi sudu pada saat terjadi kerusakan/patah, mempengaruhi jarak sudu terlempar. Dengan menggunakan diketahui bahwa posisi sudu patah pada posisi 3150 akan menghasilkan jarak lempar yang terjauh (Erich Hau, 2. Gambar 3-1d memperlihatkan bahwa jarak lempar dalam arah mendatar maksimal sebesar 73 m ke kanan pada posisi 315 derajat saat sudu Dengan disimpulkan bahwa radius aman SKEA 50 kW pada putaran desain rotor 80 rpm adalah sebesar 73 m dan hal tersebut terjadi jika sudu patah pada posisi 315 derajat sesuai dengan prediksi. Radius aman sejauh 73 m sebenarnya tergolong kecil. Hal tersebut disebabkan oleh putaran rotor yang Jika sebuah turbin angin tidak menggunakan sistem roda gigi . ear bo. sebagai penghubung antara rotor dan generator maka semakin besar kapasitas turbin angin tersebut, semakin besar diameter rotor yang diperlukan. Dengan adanya batasan harga Tip Speed Ratio maka semakin besar diameter rotor, semakin lambat putarannya (David M. Eggleston & Forrest S. Stoddard, 1987. Erich Hau, 2006. Tony Burton, et al. Dengan kata lain semakin besar Jurnal Teknologi Dirgantara Vol. 8 No. 1 Juni 2010:57-63 kapasitas turbin angin, maka semakin kecil radius aman yang diperlukan. Namun yang perlu diingat bahwa semakin besar kapasitas turbin angin, biasanya sudu yang digunakan semakin berat sehingga momentum yang dibawa oleh setiap sudu juga menjadi besar. Distribusi Massa . Nodal . 2 Putaran Sudu Saat Terlempar C Diameter rotor = 17,00 m C Chord di pangkal = 0,70 m C Chord di ujung = 0,35 m Hasil perhitungan luas airfoil persegmen dan distribusi volumenya disajikan dalam Gambar 3-2 berikut ini. Distribusi Penampang 0,05 0,04 0,03 Posisi . 0,00 Gambar 3-2a: Distribusi 0,035 0,020 0,015 0,010 0,005 Nodal . 0,000 Gambar 3-2b: Distribusi volume sudu Dengan menjumlahkan distribusi volume di sepanjang sudu, maka diperoleh volume sudu sebesar 0,287 m3. Dengan menggunakan data berat sudu 230 kg, maka diperoleh rapat jenis sudu sebesar 802 kg/m3. Dengan mengasumsikan rapat jenis sudu konstan di sepanjang sudu, maka distribusi berat sudu menjadi seperti pada Gambar 3-3 berikut ini: Nodal . Gambar 3-4: Distribusi . /de. 0,025 Distribusi Tangensial Distribusi Volume 0,030 Pada paragraf sebelumnya telah diasumsikan bahwa rapat jenis sudu konstan di sepanjang sudu. Oleh karena itu distribusi massa pada Gambar 3-3 akan sebanding dengan distribusi volume pada Gambar 3-2b. Pada kondisi putaran konstan, kecepatan tangensial di setiap segmen hanya bergantung pada posisi nodal segmen terhadap sumbu putar rotor. Distribusi kecepatan tangensial di setiap segmen disajikan dalam Gambar 3-4 di bawah ini. 0,02 0,01 Gambar 3-3: Distribusi massa Sudu SKEA 50 kW secara umum mempunyai data geometri: 0,06 Pada Sub Bab Batasan Masalah telah disebutkan bahwa titik berat sudu berada pada posisi 2,907 m dihitung dari sumbu putar rotor. Dengan demikian jarak antara nodal setiap segmen terhadap titik berat sudu adalah sebagai berikut: -1 0 Jarak Nodal ke c. Nodal . Gambar 3-5: Jarak nodal ke titik berat Analisis Radius Aman Akibat Kegagalan. (Sulistyo Atmadi et al. Jarak nodal ke titik berat sudu pada Gambar 3-5 di atas ada yang berharga positif dan negatif. Perbedaan berpengaruh pada arah momentum putar. Perumusan momentum putar yang digunakan dalam Persamaan . merupakan fungsi dari distribusi massa, kecepatan tangensial, dan jarak antara nodal ke titik berat sudu. Ketiga parameter tersebut telah disajikan masing-masing dalam Gambar 3-3. Gambar 3-4, dan Gambar 3-5. Dengan memasukkan ketiga parameter tersebut dalam Persamaan . , maka diperoleh distribusi momentum putar seperti pada Gambar 3-6 di bawah ini. dM . /de. Distribusi Momentum Putar Nodal . Gambar 3-6:Distribusi momentum putar Sesuai dengan kondisi operasional sistem turbin angin maka apabila dilihat dari arah depan, rotor akan berputar searah jarum jam (David M. Eggleston & Forrest S. Stoddard, 1987. Erich Hau. Tony Burton, et al. , 2. Distribusi momentum putar pada Gambar 3-6 menunjukkan bahwa pada daerah pangkal sebelum titik berat sudu, momentum putarnya berharga negatif. Hal tersebut tersebut setiap elemen sudu memaksa sudu untuk berputar berlawanan arah jarum jam, demikian juga sebaliknya. Dengan diperoleh momentum putar untuk setiap sudu pada saat patah di pangkal sebesar 633 kg. /de. Harga tersebut bersesuaian dengan putaran sudu terhadap titik beratnya sendiri searah jarum jam sebesar 145 rpm. KESIMPULAN Kesimpulan yang dapat ditarik dari analisa ini adalah bahwa radius aman untuk pengujian SKEA 50 kW cukup kecil, yaitu hanya sekitar 73 m. Selain itu kecepatan jatuhnya sudu juga tidak terlalu besar karena putaran rotor yang kecil. Namun tetap perlu diingat bahwa sudu rotor SKEA 50 kW sangat berat sehingga dampak yang diakibatkan jika menabrak sesuatu dapat berakibat Selain radius aman, juga dihitung putaran sudu pada saat terlempar akibat patah pada pangkalnya. Hasil perhitungan menunjukkan bahwa sudu berputar terhadap titik beratnya sendiri sebesar 145 rpm. DAFTAR RUJUKAN Ahmad J. Fitroh, 2009. Perancangan dan Analisis Aerodinamika Sudu SKEA 50 kW. LAPAN. Bogor. Anderson. John Jr. Fundamentals of Aerodynamics. Mc. Graw Hill company. Singapore. David M. Eggleston & Forrest S. Stoddard, 1987. Wind Turbine Engineering Design. Van Nostrand Reinhold. New York. Erich Hau. Wind Turbine Fundamentals. Technologies. Application. Economics. Springer. Berlin. Mc. Cormick. Barnes Aerodynamics. Aeronautics, and Flight Performance. John Wiley & Sons. Canada. Milborrow. Wind Turbine Runaway Speed, 5. BWEA Wind Energy Conference. Reading. Sorensen. , 1980. Calculation of Trajectories of Detached Wind Turbine Blades and Prediction of Side Risk Levels Associated with Failures of HAWTAoS. IEA-ExpertMeeting. Munich. Sulistyo Atmadi & Ahmad J. Fitroh. Efek Penambahan Sudut Pitch Terhadap Kurva Daya SKEA 50 kW. LAPAN. Bogor. Tony Burton, et al. , 2001. Wind Energy Hand Book. John Wiley & Sons.