Penerapan Metode Numerik Finite Element Method (FEM) dalam Analisis Struktur Menara Turbin Angin – Final Project Komputasi Teknik – Ganang Adhitama (2606829613)

Assalamualaikum warahmatullahi wabarakatuh.

Bismillahirrahmanirrahim.

Pada kesempatan kali ini, saya ingin membahas mengenai metode numerik yang digunakan dalam proyek analisis struktur menara turbin angin, khususnya Finite Element Method (FEM). Pembahasan akan dimulai dari prinsip dasar metode numerik, persamaan diferensial yang mendasarinya, hingga algoritma matematis yang digunakan oleh software untuk menghitung perpindahan dan tegangan.

Fokus pembahasan: Memahami bagaimana software simulasi seperti Autodesk Inventor menghitung deformasi dan tegangan suatu struktur dengan geometri kompleks, persamaan matematika yang mendasarinya, serta algoritma FEM yang digunakan untuk menyelesaikan permasalahan tersebut.

1. Pendahuluan

Dalam bidang teknik mesin, banyak permasalahan yang dinyatakan menggunakan persamaan matematika, khususnya persamaan diferensial. Misalnya, aliran fluida, perpindahan panas, getaran, maupun deformasi struktur akibat pembebanan.

Untuk geometri yang sederhana, solusi dapat dihitung secara analitik. Sebagai contoh, perpindahan ujung suatu balok kantilever dengan penampang seragam dan beban sederhana dapat dihitung menggunakan teori balok. Namun, ketika geometri menjadi lebih kompleks, penyelesaian analitik sering kali sulit dilakukan.

Pada proyek ini, objek yang dianalisis adalah menara turbin angin yang dimodelkan sebagai geometri tiga dimensi menggunakan Autodesk Inventor Professional 2024. Tujuan utamanya adalah menentukan distribusi perpindahan (displacement) dan tegangan (stress) pada menara akibat pembebanan.

Oleh karena itu, digunakan metode numerik untuk mengubah persoalan matematika kontinu menjadi sejumlah persamaan yang dapat diselesaikan menggunakan komputer.

Gambar 1. Ilustrasi Analisis FEM Struktur Menara Turbin Angin

Pilih mode visualisasi untuk melihat representasi mesh, displacement, dan von Mises stress.

Mesh Displacement von Mises Stress

Interpretasi: Mode Mesh memperlihatkan pembagian menara menjadi elemen dan node sebagai dasar perhitungan FEM.

Ilustrasi konseptual yang menyerupai tampilan software Finite Element Analysis. Warna, bentuk deformasi, serta pola tegangan tidak berasal dari hasil perhitungan Autodesk Inventor dan tidak memiliki skala numerik.

Pada analisis FEM, model struktur terlebih dahulu dibagi menjadi elemen-elemen kecil melalui proses meshing. Software kemudian menyelesaikan perpindahan pada setiap node dan menggunakannya untuk menghitung regangan serta tegangan. Hasil tersebut dapat ditampilkan dalam bentuk kontur warna untuk membantu mengidentifikasi daerah yang mengalami deformasi atau tegangan relatif lebih tinggi.

2. Mengenal FDM, FVM, dan FEM

Terdapat beberapa pendekatan dalam menyelesaikan persamaan diferensial secara numerik. Tiga pendekatan yang umum digunakan dalam perangkat lunak teknik adalah Finite Difference Method (FDM), Finite Volume Method (FVM), dan Finite Element Method (FEM).

FDM

Finite Difference Method

Mengganti turunan diferensial dengan aproksimasi beda hingga pada titik grid.

Contoh: Persoalan temperatur pada domain sederhana.

FVM

Finite Volume Method

Menerapkan neraca konservasi pada setiap volume kontrol melalui fluks yang melewati batasnya.

Contoh: Aliran fluida dan perpindahan panas.

FEM

Finite Element Method

Membagi domain menjadi elemen dan memperkirakan variabel melalui fungsi interpolasi.

Contoh: Deformasi dan tegangan struktur.

2.1. Alasan Pemilihan FEM pada Proyek Menara Turbin Angin

Dalam proyek analisis struktur menara turbin angin ini, metode yang digunakan adalah Finite Element Method (FEM). Pemilihan tersebut didasarkan pada karakteristik permasalahan yang dianalisis, yaitu mekanika struktur dengan geometri tiga dimensi.

Menara turbin angin menerima berbagai pembebanan yang menghasilkan perpindahan dan distribusi tegangan pada seluruh bagian struktur. Melalui FEM, geometri menara dapat dibagi menjadi elemen-elemen kecil sehingga respons mekanis pada berbagai lokasi dapat dihitung secara numerik.

Dibandingkan dengan FDM, FEM memiliki fleksibilitas yang lebih baik dalam mengikuti geometri yang tidak beraturan dan menerapkan kondisi batas mekanika struktur. Sementara itu, FVM sangat umum digunakan pada permasalahan konservasi dan aliran fluida, meskipun metode tersebut juga dapat diterapkan pada mekanika padatan.

FEM juga memungkinkan evaluasi perpindahan dan tegangan pada lokasi tertentu, misalnya perpindahan puncak menara dan tegangan pada bagian dasar yang mengalami pembebanan lentur. Oleh karena itu, FEM merupakan pendekatan yang sesuai untuk proyek ini dan tersedia dalam fitur Stress Analysis pada Autodesk Inventor Professional 2024.

3. Pemodelan Menara Turbin Angin

Struktur menara dimodelkan sebagai 3D solid dengan bagian dasar ditahan secara tetap (fixed constraint). Dengan demikian, perilaku menara secara keseluruhan menyerupai struktur kantilever.

Pembebanan yang direncanakan mempertimbangkan gaya dorong rotor (thrust force), gaya hambat angin pada menara (drag force), dan berat turbin. Acuan pembebanan direncanakan menggunakan Simplified Load Model (SLM), khususnya Load Case H pada IEC 61400-2:2013 yang berkaitan dengan kondisi angin ekstrem.

Gambar 2. Ilustrasi Pembebanan Menara Turbin Angin

Thrust Force Drag Force Turbine Weight

Centang atau hilangkan pilihan gaya untuk melihat komponen pembebanan secara terpisah. Panjang panah drag force dibuat semakin besar menuju puncak menara untuk mengilustrasikan profil pembebanan angin yang meningkat terhadap ketinggian. Skema bersifat konseptual, bukan hasil simulasi maupun distribusi beban terhitung. Momen pada shaft tidak ditampilkan.

Pada ilustrasi tersebut, thrust force merepresentasikan gaya dorong akibat interaksi angin dengan rotor turbin. Turbine weight merupakan gaya vertikal ke bawah akibat berat turbin. Sementara itu, drag force merupakan pembebanan lateral akibat interaksi aliran angin dengan permukaan menara.

Panjang panah drag force pada ilustrasi dibuat meningkat ke arah atas untuk menggambarkan kecenderungan kecepatan angin bertambah terhadap ketinggian. Akan tetapi, distribusi gaya aktual tetap bergantung pada profil kecepatan angin, koefisien drag, dan dimensi menara pada setiap ketinggian.

3.1. Apa yang Dimaksud Linear Static Analysis?

Linear static analysis memiliki dua asumsi utama. Pertama, perilaku material diasumsikan elastis linier dan deformasinya cukup kecil. Kedua, beban dianalisis dalam kondisi kesetimbangan statik sehingga efek inersia dan redaman diabaikan.

Linear

Tegangan berhubungan linier dengan regangan melalui Hukum Hooke. Perubahan geometri akibat deformasi dianggap kecil.

Static

Gaya dan momen berada dalam kesetimbangan. Efek percepatan, getaran, dan respons dinamik tidak diperhitungkan.

Dengan pendekatan ini, hasil utama simulasi adalah perpindahan dan distribusi tegangan. Analisis ini belum mencakup respons dinamik, tekuk kritis, atau ketidaklinieran material dan geometri.

4. Persamaan Diferensial yang Diselesaikan

Setelah mengetahui model fisiknya, langkah berikutnya adalah memahami persamaan matematika yang mendeskripsikan deformasi struktur.

4.1. Persamaan Kesetimbangan

∇ · σ + b = 0

σ = tensor tegangan (Pa); b = gaya badan per satuan volume (N/m³); ∇ · = operator divergensi. Persamaan menyatakan keseimbangan gaya pada elemen material.

Secara sederhana, gaya internal akibat tegangan harus mampu menyeimbangkan gaya eksternal. Untuk kasus statik, percepatan tidak muncul dalam persamaan.

4.2. Hubungan Regangan dan Perpindahan

ε = ½ [ ∇u + (∇u)^T ]

ε = tensor regangan (tanpa satuan); u = vektor perpindahan (m); ∇u = gradien perpindahan (tanpa satuan); T = transpose matriks.

Regangan bukan hanya ditentukan oleh besarnya perpindahan, melainkan oleh bagaimana perpindahan berubah dari satu lokasi ke lokasi lainnya.

4.3. Hukum Hooke

σ = Dε

D = matriks elastisitas material (Pa); σ = tegangan (Pa); ε = regangan (tanpa satuan). Bentuk ini menggunakan representasi matriks dari hubungan tegangan-regangan.

Untuk material isotropik, sifat elastis dapat dinyatakan melalui modulus Young E dan rasio Poisson ν. Artinya, sifat material menentukan seberapa besar tegangan yang terbentuk akibat suatu regangan.

4.4. Persamaan Navier–Cauchy

Dengan menggabungkan persamaan keseimbangan, hubungan regangan-perpindahan, dan Hukum Hooke, untuk material homogen dan isotropik diperoleh persamaan Navier–Cauchy:

μ∇²u + (λ + μ)∇(∇ · u) + b = 0

u = vektor perpindahan (m); λ dan μ = konstanta Lamé (Pa); ∇²u = Laplasian perpindahan (1/m); b = gaya badan (N/m³). Persamaan ini berlaku untuk elastisitas linier homogen dan isotropik.

μ = E / [2(1 + ν)] λ = Eν / [(1 + ν)(1 – 2ν)]

E = modulus Young (Pa); ν = rasio Poisson (tanpa satuan); λ dan μ = konstanta Lamé (Pa).

Persamaan Navier–Cauchy merupakan sistem persamaan diferensial parsial orde dua dengan perpindahan sebagai variabel utama.

4.5. Laplasian Perpindahan 3D

Pada model 3D solid, perpindahan memiliki tiga komponen, yaitu u_x, u_y, dan u_z. Operator Laplasian tiga dimensi adalah:

∇² = ∂²/∂x² + ∂²/∂y² + ∂²/∂z²

Untuk perpindahan arah x:

∇²u_x = ∂²u_x/∂x² + ∂²u_x/∂y² + ∂²u_x/∂z²

u_x = perpindahan arah x (m); x, y, z = koordinat posisi (m); ∇²u_x = Laplasian perpindahan arah x (1/m). Komponen arah y dan z memiliki bentuk serupa.

Ketiga komponen tersebut saling berhubungan melalui persamaan Navier–Cauchy. Jadi, software tidak menyelesaikan ketiganya sebagai masalah yang terpisah, melainkan sebagai sistem persamaan yang saling terkait.

Konsep utama: Navier–Cauchy menjelaskan persoalan fisiknya, sedangkan FEM menyediakan cara numerik untuk mencari solusi perpindahannya.

5. Bagaimana Algoritma FEM Bekerja?

Bagian ini merupakan inti pembahasan. FEM tidak sekadar membuat mesh dan menjalankan perintah Solve. FEM mengubah persamaan diferensial menjadi sistem persamaan aljabar yang dapat diselesaikan menggunakan komputer.

Tahapan algoritma berikut ditampilkan secara berurutan. Setiap tahap menjelaskan masukan, proses perhitungan, serta keluaran yang diperlukan oleh tahap selanjutnya.

Tahap 1 — Strong Form menjadi Weak Form

Persamaan Navier–Cauchy termasuk strong form, karena menyatakan persamaan diferensial yang harus dipenuhi pada setiap titik domain.

Dalam FEM, persamaan kesetimbangan dikalikan dengan fungsi uji berupa perpindahan virtual, kemudian diintegralkan. Melalui integrasi parsial, diperoleh bentuk integral yang melibatkan turunan pertama perpindahan.

∫_Ω δε : σ dΩ = ∫_Ω δu · b dΩ + ∫_Γt δu · t dΓ

Ω = domain volume (m³); δε = regangan virtual; δu = perpindahan virtual (m); t = traction permukaan (Pa); Γt = bagian permukaan tempat traction diberikan. Simbol : menyatakan kontraksi tensor.

Input: Persamaan kesetimbangan, hubungan material, dan kondisi batas.

Output: Persamaan weak form yang cocok untuk diskretisasi FEM.

Tahap 2 — Diskretisasi Mesh

Domain kontinu menara dibagi menjadi elemen-elemen kecil yang saling terhubung melalui node. Untuk analisis solid 3D, Autodesk Inventor Stress Analysis menggunakan elemen tetrahedral kuadratik 10-node.

Setiap node memiliki tiga derajat kebebasan translasi: u_x, u_y, dan u_z. Dengan 10 node, suatu elemen mempunyai 30 komponen perpindahan nodal lokal.

Input: Geometri struktur dan parameter mesh.

Output: Bentuk elemen, koordinat node, dan konektivitas antar-elemen.

Tahap 3 — Interpolasi Shape Function

Perpindahan pada setiap titik elemen belum diketahui. FEM mengaproksimasinya dari perpindahan pada node-node elemen menggunakan shape function.

u^(e) = Nd^(e)

N = matriks fungsi bentuk (tanpa satuan); d^(e) = perpindahan nodal elemen (m); u^(e) = medan perpindahan elemen (m). Superskrip (e) berarti elemen tertentu.

Turunan shape function selanjutnya membentuk matriks B, yaitu matriks yang menghubungkan perpindahan nodal dengan regangan.

Input: Node dan geometri elemen.

Output: Matriks N dan B untuk formulasi kekakuan elemen.

Tahap 4 — Pembentukan Matriks Kekakuan Elemen

Matriks B dan sifat material D dimasukkan ke dalam weak form untuk memperoleh matriks kekakuan setiap elemen.

K^(e) = ∫_Ωe B^TDB dΩ

K^(e) = matriks kekakuan elemen (N/m); B = matriks regangan-perpindahan (1/m); D = matriks elastisitas (Pa); Ωe = volume elemen (m³).

Integral umumnya dievaluasi dengan metode numerik seperti Gaussian quadrature. Dari sini diperoleh hubungan antara gaya dan perpindahan pada elemen.

Input: Matriks B, sifat material D, dan volume elemen.

Output: Matriks kekakuan dan kontribusi gaya setiap elemen.

Tahap 5 — Assembly dan Boundary Conditions

Semua matriks kekakuan elemen dijumlahkan ke dalam matriks kekakuan global berdasarkan konektivitas node.

Selanjutnya diterapkan pembebanan dan kondisi batas. Pada dasar menara yang fixed:

u_x = u_y = u_z = 0

u_x, u_y, dan u_z adalah perpindahan translasi node dasar (m).

Input: Seluruh matriks elemen, beban, dan pengekangan.

Output: Sistem persamaan global yang siap diselesaikan.

Tahap 6 — Linear Equation Solver

Setelah assembly dan kondisi batas, FEM menghasilkan sistem persamaan aljabar:

KU = F

K = matriks kekakuan global (N/m); U = vektor perpindahan nodal (m); F = vektor gaya nodal (N).

Komputer menggunakan algoritma aljabar linear, misalnya faktorisasi atau metode iteratif, untuk memperoleh nilai perpindahan pada seluruh node yang tidak dikekang.

Input: Matriks K dan vektor F.

Output: Nilai perpindahan nodal U.

Tahap 7 — Post-processing: Strain dan Stress

Setelah nilai perpindahan nodal diketahui, regangan dan tegangan pada setiap elemen dapat dihitung:

ε^(e) = Bd^(e)
σ^(e) = Dε^(e)

ε^(e) = regangan elemen; σ^(e) = tegangan elemen (Pa); d^(e) = perpindahan nodal hasil solver (m).

Hasil tersebut kemudian digunakan untuk menampilkan kontur deformasi, regangan, tegangan, dan tegangan ekuivalen von Mises.

Input: Perpindahan nodal yang sudah diketahui.

Output: Distribusi displacement, strain, dan stress.

5.1. Ringkasan Algoritma FEM

Governing PDE
Persamaan Navier–Cauchy

↓

Weak Formulation
Persamaan diferensial menjadi integral

↓

Mesh Discretization
Elemen dan node

↓

Shape Function
u = Nd

↓

Element Stiffness
K^(e) = ∫ B^TDB dΩ

↓

Assembly + Boundary Conditions

↓

Linear Equation Solver
KU = F

↓

Displacement → Strain → Stress

Gambar 3. Hubungan antartahap algoritma FEM.

5.2. Contoh Sederhana: Penyelesaian Persamaan FEM

Sebagai ilustrasi, bayangkan model elastis yang telah disederhanakan menjadi dua derajat kebebasan. Matriks kekakuan dan gaya dapat berupa:

K = [ [3, −1], [−1, 2] ] kN/mm
F = [5, 0]^T kN

Penyelesaian persamaan KU = F menghasilkan:

U = [2, 1]^T mm

K = matriks kekakuan contoh; F = vektor gaya; U = vektor perpindahan. Nilai ini hanya contoh aljabar dua derajat kebebasan, bukan hasil simulasi menara.

Artinya, perpindahan pada dua derajat kebebasan yang dicontohkan adalah 2 mm dan 1 mm. Konsep yang sama diterapkan pada FEM 3D, tetapi jumlah persamaannya dapat mencapai ribuan hingga jutaan derajat kebebasan.

6. Verifikasi Numerik melalui Mesh Convergence

Karena FEM merupakan metode pendekatan, hasil perhitungan dapat dipengaruhi oleh ukuran dan kualitas mesh. Salah satu metode verifikasi numerik yang penting adalah mesh convergence study.

Prinsipnya adalah menjalankan simulasi menggunakan mesh yang semakin halus dan membandingkan hasil yang diperoleh. Jika perubahan hasil semakin kecil, maka solusi menunjukkan kecenderungan konvergen.

Perubahan (%) = |R_fine − R_coarse| / |R_fine| × 100%

R_fine = respons pada mesh lebih halus; R_coarse = respons pada mesh sebelumnya. Keduanya harus mewakili besaran dan lokasi yang sama.

6.1. Eksperimen Interaktif Mesh Refinement

Geser slider berikut untuk melihat ilustrasi bagaimana penghalusan mesh dapat menyebabkan hasil numerik semakin stabil.

Tingkat penghalusan mesh: Sedang

1.024 Jumlah elemen ilustratif

10.50 mm Perpindahan ilustratif

3.81% Perubahan dari mesh sebelumnya

Seluruh angka di atas adalah data contoh untuk menjelaskan konsep konvergensi, bukan hasil perhitungan Autodesk Inventor atau model menara. Grid yang ditampilkan adalah visualisasi 2D untuk mempermudah pemahaman.

6.2. Penerapan Mesh Convergence pada Menara

Untuk proyek menara turbin angin, mesh convergence dapat dilakukan dengan menggunakan beberapa tingkat mesh, misalnya kasar, sedang, halus, dan sangat halus. Parameter yang dibandingkan dapat berupa perpindahan puncak menara dan tegangan pada lokasi evaluasi tertentu.

Perhatian terhadap tegangan maksimum: Tegangan pada sudut tajam, titik pembebanan ideal, atau kondisi pengekangan tertentu dapat menunjukkan stress singularity. Nilainya mungkin terus meningkat meskipun mesh semakin halus. Karena itu, konvergensi tegangan harus dievaluasi dengan hati-hati, tidak hanya melihat satu angka tegangan maksimum global.

Mesh convergence merupakan bagian dari verification, yaitu memastikan bahwa solusi numerik tidak terlalu dipengaruhi diskretisasi. Sementara itu, validation menilai apakah model mampu menggambarkan kondisi fisik nyata, misalnya melalui perbandingan dengan data eksperimen.

7. Kesimpulan

Berdasarkan pembahasan tersebut, metode numerik yang dipilih dalam proyek analisis struktur menara turbin angin adalah Finite Element Method (FEM). Metode ini sesuai untuk analisis struktur 3D dengan geometri dan kondisi pembebanan kompleks.

Persoalan mekanikanya berawal dari persamaan kesetimbangan elastisitas, yang untuk material homogen isotropik dapat dinyatakan melalui persamaan Navier–Cauchy. Melalui weak formulation dan diskretisasi elemen, FEM mengubah PDE menjadi sistem persamaan aljabar linear KU = F.

Sistem tersebut diselesaikan untuk memperoleh perpindahan nodal, yang kemudian digunakan untuk menghitung regangan dan tegangan. Mesh convergence selanjutnya diperlukan untuk memeriksa kestabilan hasil numerik.

Inti yang saya pelajari: Software simulasi bukan sekadar alat untuk menghasilkan gambar kontur tegangan dan deformasi. Di baliknya terdapat formulasi mekanika kontinu, diskretisasi numerik, pembentukan matriks kekakuan, dan algoritma penyelesaian persamaan aljabar.

Referensi

International Electrotechnical Commission. (2013). IEC 61400-2:2013 — Wind turbines — Part 2: Small wind turbines.
IEC Webstore
Wood, D. (2011). Small wind turbines: Analysis, design, and application. Springer.
Hibbeler, R. C. (2011). Mechanics of materials (8th ed.). Pearson Prentice Hall.
Budynas, R. G., & Nisbett, J. K. (2011). Shigley’s mechanical engineering design (9th ed.). McGraw-Hill.
Karmouche, S. (2016). Tower design and analysis for a small wind turbine (Capstone design final report). Al Akhawayn University.
Autodesk Support. Which theoretical method (FEA) is used by Autodesk Inventor for stress analysis?
Autodesk Technical Article
Autodesk Support. Accuracy of stress analysis results and convergence in Inventor Stress Analysis.
Autodesk Technical Article

Demikian pembahasan mengenai metode numerik FEM dalam analisis struktur menara turbin angin. Semoga artikel ini dapat membantu memahami hubungan antara persamaan diferensial, algoritma numerik, dan penerapannya dalam software simulasi teknik.

Wassalamualaikum warahmatullahi wabarakatuh.

ccitonline.com

Quick Links