Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS IJCCS.
Vol.
No.
July x, pp.
GEOMETRIC MORPHOMETRIC BANANA LEAF
MIDRIB AS WIND TURBINE AIRFOIL PROFILE
Firlya Rosa *1,2.
Priyono Soetikno1,3,4.
I Wayan Suweca1.
Mochammad Agoes Moelyadi 1 Mechanical Engineering Study Program.
Faculty of Mechanical and Aerospace Engineering.
Bandung Institute of Technology.
Bandung, 40132.
Indonesia Department of Mechanical Engineering.
Faculty of Engineering.
Universitas Bangka Belitung, 33172.
Indonesia Fluid Machinery Laboratory.
Faculty of Mechanical and Aerospace Engineering New and Renewable Research.
Bandung Institute of Technology.
Bandung, 40132.
Indonesia e-mail: *1 f105a@yahoo.
Abstrak Profil airfoil merupakan salah satu faktor yang mempengaruhi kinerja turbin angin sumbu Sejumlah penelitian telah mengadopsi bentuk dan gerak binatang baik binatang air maupun binatang yang terbang.
Anatomi alami pelepah daun pisang menjadi salah satu inspirasi untuk pemodelan profil airfoil.
Profil pelepah dari tujuh sampel daun pisang dipetakan dengan menggunakan teknik morfometri geometrik berbasis landmark.
Analisis sampel menggunakan metode statistik dan analisis aerodinamis airfoil menggunakan perangkat lunak QBLADE dengan perubahan bilangan Reynolds mulai dari 3E5 hingga 7E5 pada interval 1E5.
Ada tujuh sampel dengan variasi profil dan variasi aerodinamis yang konstan.
Dari hasil penelitian didapatkan bahwa semua model memiliki koefisien angkat tinggi, koefisien hambat rendah, dan glide ratio Sampel 1 pada bilangan Reynolds 7E5 memiliki glide ratio tertinggi dengan nilai 123,07.
Dari hasil analisa aerodinamis didapatkan bahwa dalam kondisi kecepatan angin rendah, profil airfoil berdasarkan morfologi pelepah pisang dapat diterima untuk turbin angin.
Kata kunciAiairfoil, pelepah daun pisang, glide ratio, morfologi Abstract The airfoil profile is one of the factors influencing the performance of horizontal-axis wind Numerous studies have adopted the form and motion of aquatic and aerial animals.
The natural anatomy of the midrib of banana leaves is one of the inspirations for modeling the airfoil The midrib profiles of seven samples of banana leaves were mapped using a landmarkbased geometric morphometric technique.
Sample analysis utilizing statistical methods and aerodynamic analysis of airfoils operating QBLADE software with Reynolds number changes ranging from 3E5 to 7E5 at 1E5 intervals.
There are seven samples with constant aerodynamic profiles and variants.
All models have high lift coefficients, low drag coefficients, and high glide ratios, according to the results.
Sample 1 at Reynolds number 7E5 had the highest glide ratio with a value of 123.
In conditions of low wind speed, an airfoil profile based on the morphology of the banana midrib may be accepTabel for wind turbines.
KeywordsAi airfoil, banana leaf midrib, glide ratio, morphology Received June 1st,2012.
Revised June 25th, 2012.
Accepted July 10th, 2012 Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS
ISSN: 1978-1520214
PENDAHULUAN
enomena energi yang dihasilkan dari sumber energi terbarukan semakin meningkat mengingat sumber daya fosil semakin berkurang.
Beberapa energi terbarukan mulai dikembangkan dalam rangka mendukung energi yang ramah lingkungan.
Pembangkit listrik yang berasal dari energi terbarukan mulai banyak dikembangkan, diantaranya bersumber pada angin, surya, biomassa dan lain sebagainya.
Performansi pembangkit listrik tenaga air/bayu (PLTB) dipengaruhi berbagai aspek yang salah satunya adalah profil airfoil.
Aerodinamis mempertimbangkan geometri airfoil, geometri sudu dan distribusi sudut puntir pada sudu serta massa sudu yang direncanakan.
Pertimbangan-pertimbangan tersebut diharapkan mampu menghasilkan gaya angkat tinggi dan gaya hambat rendah.
Berbagai geometri airfoil didapatkan dari berbagai profil yang berasal dari alam dan terus dikembangkan, terutama morfology airfoil dari pergerakan binatang.
Beberapa pergerakan hewan yang diadaptasi diantaranya gerakan hewan di air dan di udara.
Airfoil yang berasal dari morfologi dari hewan air diantaranya gerakan kepala ikan lumba-lumba yang menyerupai bentuk airfoil pada bagian ujungnya .
Begitu juga bentuk dari ikan sturgeon menyerupai profil airfoil .
Airfoil yang diadaptasi dari gerakan hewan di udara telah banyak diteliti, diantaranya morphology airfoil yang diadaptasi dari burung dan capung.
Diantaranya adalah morfologi burung hantu.
, kepakan sayap burung.
dan morfologi sayap capung Aeshna Cyanea.
Untuk penelitian profil airfoil yang berasal dari morphology tumbuhan masih sedikit yang dilakukan, padahal Martin Ryle telah meramalkan bahwa jika blade dapat bekerja seperti pohon palem dalam kecepatan angin tinggi sehingga efisiensi konversi energi dapat ditingkatkan hingga 50%, dan material dapat dihemat sebesar 2/3.
Bagian dari tumbuhan memiliki karakteristik akan bergerak mengikuti gerakan arah angin dibandingkan dengan gerakan hewan yang mengikuti pergerakan naluri hewan itu sendiri.
Pelepah daun pisang merupakan salah satu morfologi tumbuhan yang dapat dijadikan profil airfoil.
Beberapa penampang melintang daun yang diteliti, salah satunya adalah penampang berbentuk U yang sebanding dengan pelepah daun pisang dengan jenis Musa-x paradisiaca.
Dengan menggunakan metode numerik terhadap pengaruh bentuk penampang struktur penahan beban utama pada fungsionalitas mekanis, didapatkan susunan penampang pelepah daun pisang mempunyai kekakuan lentur yang lebih baik dibandingkan tangkai daun lainnya.
Selain itu, kombinasi bentuk dan struktur bagian dalam pelepah daun pisang juga mengurangi kekakuan torsi yang menyebabkan rasio puntir kelengkungan yang sangat tinggi dibandingkan dengan tangkai daun dari berbagai spesies pohon sehingga kemampuan pelepah daun pisang dalam menerima torsi berpotensi untuk dapat dijadikan bentuk profil airfoil.
Gambar 1 Arah penekungan pelepah pisang dikarenakan angin (A) pandangan atas tangkai daun (B) peningkatan kecepatan angin daun membelok menjauh dari angin (C) penekukan daun akibat penambahan kecepatan angin.
(D) rekonstruksi aliran angin pada penampang pelepah daun pisang IJCCS Vol.
No.
July 201x : first_page Ae end_page Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS IJCCS
ISSN: 1978-1520
Selain itu, pemilihan penampang pelepah daun pisang dikarenakan kecenderungan pelepah mengikuti pergerakan angin.
Kondisi pelepah angin yang mengalami penekukan mengikuti gerakan angin .
menyebabkan pelepah tidak mudah patah seperti terlihat pada Gambar 1.
Umumnya, daun lebar akan mengalami penekukan dan berkibaryang disebabkan oleh angin berkecepatan 5 m/s sebesar 46.
54Cdan 52.
56C pada kecepatan angin 1 m/s.
Oleh karena itu, penampang pelepah pisang yang digunakan dalam penelitian ini difokuskan menggunakan pelepah yang dipotong pada sudut 45C dengan asumsi gaya angin yang terjadi akan tegak lurus terhadap permukaan penampang pelepah.
Adapun penampang pelepah yang terpotong terlihat seperti pada Gambar 2.
Tidak semua morfologi pelepah daun pisang yang dijadikan profil airfoil.
Hal ini diasumsikan bahwa gaya angin hanya mengenai salah satu sisi pelepah daun pisang mulai dari tangkai daun pisang .
tem en.
sampai dengan bagian yang sejajar dengan tangkai daun seperti yang diilustrasikan pada Gambar 3.
Oleh karena itu, pemotongan profil yang dijadikan profil airfoil dimulai dari stem end sampai dengan daerah yang sejajar dengan sumbu pelepah .
itik A sampai dengan titik B) seperti yang ditunjukkan pada Gambar 3.
Gambar 2 .
Metode pemotongan 45C longitudinal profil airfoil berdasarkan morfologi pelepah pisang .
Profil penampang pelepah daun pisang Wind Direction Stem end Gambar 3 .
Profil airfoil pelepah daun pisang .
penentuan daerah pemotongan .
permukaan yang dijadikan airfoil
METODE PENELITIAN
Pemetaan profil dan ukuran tumbuhan umumnya menggunakan landmark-based geometric morphometric methods (GMM.
Ae.
dengan menggunakan 7 sampel dari pohon pisang yang sama dengan profil pelepah seperti pada Gambar 4.
Metode landmark akan memetakan titik-titik terluar pada bagian pelepah daun pisang yang berbentuk 3 dimensi untuk dijadikan profil 2 dimensi.
Penentuan tinggi profil penampang Title of manuscript is short and clear, implies research results (First Autho.
Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS
ISSN: 1978-1520216
pelepah pisang untuk semua sampel sama yang mulai dari stem end di titik A sampai dengan titik Dengan menempatkan koordinat x,y dijadikan sebagai acuan dengan titik arah koordinat y yang berjarak tetap.
Jarak masing-masing koordinat pada arah sumbu yang sama, memastikan tinggi dari profil akan sama mulai dari titik awal sampai dengan stem end.
Namun, jarak koordinat masing-masing profil arah sumbu x akan berbeda dikarenakan lebar profil pelepah daun pisang yang berbeda-beda.
Arah sumbu y akan menandakan panjang chord profil airfoil dan arah sumbu x akan mengindikasikan tinggi camber dari profil airfoil.
Ada 19 titik landmark yang terpilih pada pelepah daun pisang.
Landmark 1 dan 19 merupakan titik stem end pelepah daun pisang.
Landmark 2 Ae 18 merupakan titik pada jarak yang tetap pada koordinat y dengan tahapan proses pemetaan seperti pada Gambar 5.
Hasil pemetaan, didapatkan bentuk awal profil airfoil dari 7 sampel pelepah daun pisang seperti terlihat pada Gambar 6.
Koordinat-koordinat yang didapatkan terdata pada Tabel 1 yang kemudian akan dilanjutkan analisa hubungan antara satu sampel dengan sampel yang lain dengan menggunakan metode statistik.
Gambar 4 Ilustrasi sampel pelepah daun pisang .
Penamaan pelepah daun pisang .
Proses pemilihan dan penamaan sampel pelepah daun pisang Keterangan:
S1=Sample 1 S2=Sample 2 S3=Sample 3 S4=Sample 4 S5=Sample 5 S6=Sample 6 S7=Sample 7 .
Gambar 5 Tahapan proses pemetaan masing-masing morfologi pelepah daun pisang .
penyamaan tinggi profil .
hasil pemetaan profil .
penetapan koordinat y .
penampang yang terpotong Gambar 6 Profile 7 sampel airfoil IJCCS Vol.
No.
July 201x : first_page Ae end_page Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS IJCCS
ISSN: 1978-1520
Posisi Titik Tabel 1.
Koordinat .
hasil pemetaan Posisi Titik Koordinat Tetap HASIL DAN DISKUSI y-coordinate Analisa Data Penyebaran data untuk masing-masing koordinat titik tetap mulai posisi 1 sampai dengan 19 berdasarkan Tabel 1 yang terlihat pada grafik di Gambar 7.
Dari Gambar 7 terlihat perbedaan titik koordinat pada arah sumbu y yang menandakan adanya variasi bentuk profil masing-masing Umumnya puncak kurva berada pada posisi titik 11 dan 12 dengan camber tertinggi pada sampel 3 dan sampel 6 mempunyai nilai camber terendah.
0,15
0,05
10 11 12 13 14 15 16 17 18 19
x-coordinate Gambar 7 Sebaran titik koordinat arah sumbu y untuk masing-masing sampel Title of manuscript is short and clear, implies research results (First Autho.
Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS
ISSN: 1978-1520218
Langkah awal dengan melakukan uji korelasi menggunakan metode Pearson Correlation.
Hasil uji korelasi dengan membandingkan hasil hitung statistik dengan nilai koefisien korelasi yang dinotasikan dengan huruf r.
Uji korelasi akan diterima jika hasil uji korelasi, nilai rtabel 0.
untuk tingkat kepercayaan 99% lebih kecil dibandingkan dengan r hitung.
Selain itu dengan melihat nilai signikansi yang dihasilkan.
Jika nilai siginifikasi yang dihasilkan lebih kecil dibandingkan nilai tingkat kepercayaan sebesar 0.
01 maka uji korelasi tersebut dapat diterima.
Uji korelasi dalam penelitian ini menghasilkan nilai rtabel lebih kecil dibandingkan dengan nilai rhitung dan nilai siginifikansi lebih kecil dibandingkan dengan nilai tingkat kepercayaan.
Hasil uji korelasi dapat dilihat pada Tabel 2.
Sample1 Sample2 Sample3 Sample4 Sample5 Sample6 Sample7 Tabel 2 Test of correlation Sampl Sampl Sampl Pearson Correlation 983** .
Sig.
-taile.
Pearson Correlation .
Sig.
-taile.
Pearson Correlation Sig.
-taile.
Pearson Correlation Sig.
-taile.
Pearson Correlation Sig.
-taile.
Pearson Correlation Sig.
-taile.
Pearson Correlation Sampl Sampl Sampl 997** .
991** .
996** .
999** .
973** .
993** .
000* .
977** .
996** .
Sig.
-taile.
**.
Correlation is significant at the 0.
01 level .
-taile.
Sampl Uji normalitas dengan menggunakan uji Kolmogorov Smirnov didapatkan bahwa nilai signifikan berada di atas tingkat kepercayaan sehingga koordinat titik y di posisi tersebut terdistribusi normal untuk masing-masing sampel seperti ditunjukkan pada Tabel 3.
Dari uji kecukupan data dengan menggunakan uji nilai Kaiser-Mayer-Olkin of sampling Adequacy (KMO MSA) dan hubungan korelasi yang kuat antar variabel menunjukkan nilai KMO MSA lebih besar dibandingkan dengan nilai tingkat kepercayaan dan nilai BartlettAos test Sphericity lebih kecil dari 05 serta didapatkan nilai anti-image correlation antar variabel lebih besar dari nilai tingkat IJCCS Vol.
No.
July 201x : first_page Ae end_page Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS IJCCS
ISSN: 1978-1520
Dari hasil tersebut dapat dinyatakan bahwa data sampel memiliki nilai kecukupan dan berkorelasi antar variabel sehingga dapat dianalisa lebih lanjut.
Tabel 3 Normality distribution test Position Y Kolmogorov-Smirnova Shapiro-Wilk Statistic Df Sig.
Statistic Sig.
Position2 Position3 Position4 Position5 Position6 Position7 Position8 Position9 Position10 Position11 Position12 Position13 Position14 Position15 Position16 Position17 Position18 This is a lower bound of the true significance.
Lilliefors Significance Correction Selanjutnya, uji MANOVA digunakan untuk menganalisis garis sampel dengan variabel terikat y.
Uji statistik MANOVA (Multivariate Analysis of Varianc.
umumnya digunakan untuk menentukan keteraturan geometrik morfologi tumbuhan .
, .
Ae.
Berdasarkan hasil uji MANOVA pada Tabel 4, nilai signifikansi yang dicapai lebih besar dari nilai tingkat kepercayaan 0,05.
Hal ini menunjukkan bahwa koordinat sumbu y antar sampel memiliki hubungan bentuk Source Corrected Model Intercept Sample Error Dependent Variable Tabel 4 Tes MANOVA Tests of Between-Subjects Effects Type i Sum Mean Square of Squares Sig.
Partial Eta Squared Title of manuscript is short and clear, implies research results (First Autho.
Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS
Total
Corrected Total
ISSN: 1978-1520220
R Squared = .
000 (Adjusted R Squared = -.
R Squared = .
031 (Adjusted R Squared = -.
Analisa Aerodinamis Dari analisis data statistik didapatkan bahwa semua sampel dapat digunakan sebagai profil airfoil, dan gaya aerodinamis yang terukur.
Gaya aerodinamis diwakili oleh koefisien angkat, koefisien hambat, glide ratio, dan variabel distribusi tekanan pada permukaan airfoil.
Perangkat lunak QBLADE open source digunakan untuk melakukan analisis aerodinamis dengan bentuk profil yang akan dievaluasi, seperti yang ditunjukkan pada Gambar 8.
Parameter airfoil yang diperoleh tercantum pada Tabel 5.
Tabel 5 Specification of Samples Airfoil Thickness Camber Point
(%)
(%)
Analisa aerodinamis menggunakan XFOIL Direct Analysis dengan numerical set sebagai berikut bilangan mach sebesar 0, free transition menggunakan method Ncrit: 9 dan force transition sebesar 1 pada top dan bottom.
Untuk penentuan bilangan Reynolds didasarkan pada beberapa variabel yang diasumsikan, yaitu panjang chord airfoil sebesar 1500 mm, freestream density (A) 225 kg/m3, dynamic viscosity (A) sebesar 0.
17894x10-4 kg/ms.
Bilangan Reynolds yang didapatkan pada kecepatan angin mulai dari 3 m/s sampai dengan 7 m/s berdasarkan rumus:
ycIyce = yuUycOya yuN Gambar 8 The output of each sample QBLADE software profile IJCCS Vol.
No.
July 201x : first_page Ae end_page .
Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS IJCCS
ISSN: 1978-1520
Tabel 6 menyajikan bilangan Reynolds yang diperoleh pada kecepatan angin berkisar antara 2,92 m/s hingga 6,82 m/s, sedangkan Tabel 7 menyajikan hasil uji dinamis.
Tabel 6 Convert wind speed to the Reynolds number.
Freestream Velocity (V) Reynolds Numbers (R.
Sample Tabel 7 Aerodynamic Result of Samples Criteria Reynolds Number 00E 05 4.
00E 05 5.
00E 05 6.
00E 05 7.
00E 05
Glide Ratio
AOA
Glide Ratio
AOA
Glide Ratio
AOA
Glide Ratio
AOA
Glide Ratio
AOA
Glide Ratio
AOA
Glide Ratio
AOA
Title of manuscript is short and clear, implies research results (First Autho.
Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS
ISSN: 1978-1520222
Koefisien hambat Koefisien Angkat Simulasi aerodinamis menunjukkan bahwa koefisien angkat rata-rata yang terjadi lebih besar dari 1.
Kecuali untuk sampel nomor 3, nilai koefisien angkat rata-rata untuk bilangan Reynolds 4E5 hingga 7E5 berkisar antara 0,96 dan 0,98.
Seperti digambarkan pada Gambar 9.
, nilai rata-rata koefisien angkat pada sampel 7 dengan bilangan Reynolds 3E5 adalah 0,9736.
Gambar menunjukkan bagaimana koefisien angkat maksimum bervariasi dengan bilangan Reynolds yang berbeda.
Profil airfoil sampel 4 menghasilkan koefisien angkat tertinggi, sedangkan profil airfoil sampel 6 menghasilkan koefisien angkat terendah.
Namun, nilai koefisien angkat sampel 1 lebih konsisten dibandingkan engan koefisien angkat sampel yang lain.
Berdasarkan Gambar 9.
, koefisien hambat kurang dari 0,25.
Sampel 1, 4, dan 5 memiliki distribusi nilai koefisien hambat yang lebih seragam.
Sampel nomor 6 memiliki koefisien hambat terendah pada bilangan Reynolds 6E7.
Hasil koefisien angkat dan koefisien hambat menunjukkan bahwa glide ratio akan lebih kecil jika selisih antara nilai koefisien angkat dan nilai koefisien hambat kecil.
Glide ratio akan semakin baik jika semakin besar signifikan rentang perbedaan antara koefisien gaya angkat dan gaya hambat.
3E 05
4E 05
5E 05
6E 05
7E 05
0,25
0,15
0,05
3E 05 4E 05 5E 05 6E 05 7E 05
Reynolds Numbers Reynolds Numbers .
Glide Ratio .
140,00 120,00 100,00 80,00 60,00 40,00 20,00 0,00 3E 05 4E 05 5E 05 6E 05 7E 05 Reynolds Numbers .
Gambar 9 Hasil analisa aerodinamis .
koefisien angkat tiap sampel .
koefisien hambat tiap sampel .
glide ratio tiap sampel Nilai glide ratio meningkat seiring dengan kenaikan bilangan Reynolds, dengan sudut serang biasanya berkisar antara 2 dan 4 derajat.
Glide ratio tertinggi didapatkan pada sampel 1 sebesar 123,07 dengan bilangan Reynolds 7E5.
Sampel dengan glide ratio terendah sebesar 70,02 pada bilangan Reynolds 3E5.
Distribusi tekanan akan menjadi lebih ekspansif dari glide ratio IJCCS Vol.
No.
July 201x : first_page Ae end_page Jurnal Mekanova : Mekanikal.
Inovasi dan Teknologi Vol 8 No.
Oktober 2022
P-ISSN : 2477-5029
E-ISSN : 2502-0498IJCCS IJCCS
ISSN: 1978-1520
yang tinggi.
Hasil glide ratio menunjukkan bahwa semua sampel cocok untuk digunakan sebagai profil airfoil turbin angin sumbu horizontal untuk mencapai nilai torsi yang tinggi.
KESIMPULAN
Beberapa studi biologi mendapatkan bahwa profil airfoil dari binatang dan tumbuhan dapat diadaptasi menjadi profil airfoil turbin angin.
Geometri morfologi tumbuhan memberikan solusi alternatif untuk profil sudu turbin angin sumbu horizontal.
Uji statistik menunjukkan bahwa profil pelepah daun pisang memiliki bentuk yang seragam.
Dengan panjang chord yang sama yaitu 1500 mm dan rentang bilangan Reynolds yang mewakili kecepatan angin dari 2,92 hingga 6,82 m/s, diperoleh nilai koefisien angkat yang tinggi, nilai koefisien hambat yang rendah, dan nilai glide ratio yang tinggi untuk semua sampel pelepah daun pisang.
Menurut analisis aerodinamis yang dilakukan dengan perangkat lunak QBLADE, glide ratio pada kecepatan angin 6,82 m/s dengan nilai di atas 100 lebih baik dibandingkan kecepatan angin yang lebih rendah untuk enam sampel.
Berdasarkan pemodelan airfoil dari sampel morfologi pelepah daun pisang, ditemukan bahwa morfologi tanaman memiliki profil yang dapat digunakan sebagai profil airfoil turbin angin sumbu horizontal.
UCAPAN TERIMA KASIH
Kami berterima kasih kepada LPDP Kementerian Keuangan RI yang telah memberikan beasiswa penuh kepada penulis pertama sebagai mahasiswa pasca sarjana di Institut Teknologi Bandung.
Juga ucapan terima kasih kepada Institut Teknologi Bandung yang telah mendanai sepenuhnya penelitian ini.
DAFTAR PUSTAKA