ccitonline.com

1. Pendahuluan

Komposit serat karbon memiliki sifat mekanik yang bergantung pada arah serat. Oleh karena itu, setiap ply perlu dimodelkan dengan sifat material yang sesuai, seperti elastisitas ortotropik. Arah serat dan urutan penyusunan ply memengaruhi kekakuan serta respons tarik laminat.

Tujuan proyek ini adalah membangun model elemen hingga uji tarik laminat T700/epoxy, mengamati responsnya terhadap displacement tarik, dan membandingkan respons numerik dengan teori laminat klasik. Hasil numerik baru dapat disajikan setelah job Abaqus dijalankan dan diperiksa.

2. Model dan Asumsi Analisis

Rancangan awal proyek menggunakan laminat delapan ply dengan susunan [0/90/+45/−45]_s. Parameter di bawah perlu dicocokkan dengan model Abaqus final dan sumber material yang digunakan.

Hal yang harus diperiksa pada model

• Konfirmasi bahwa sumbu displacement benar-benar sejajar dengan arah tarik spesimen.
• RP ujung tetap dan RP penarik harus terhubung ke permukaan grip yang sesuai.
• Periksa DOF coupling agar kondisi grip terepresentasi tanpa mengunci kontraksi lateral pada daerah ukur secara tidak sengaja.
• Orientasi material lokal dan sudut setiap ply harus sesuai dengan urutan layup.
• Konsistensi satuan geometri, material, dan displacement 5 mm perlu diperiksa.

3. Ilustrasi Uji Tarik dengan Reference Point

4. Metode Numerik: Finite Element Method

Metode numerik utama proyek ini adalah Finite Element Method (FEM) atau Metode Elemen Hingga. FEM membagi geometri menjadi elemen kecil yang terhubung pada node. Perpindahan dalam tiap elemen didekati dengan fungsi interpolasi, kemudian persamaan elemen dirakit menjadi persamaan keseimbangan global.

Untuk elastisitas linier dengan formulasi elemen yang sesuai, bentuk umum sistem diskret adalah:

K adalah matriks kekakuan global; u adalah vektor perpindahan nodal; F adalah vektor gaya nodal ekuivalen. Pada displacement control, perpindahan pada RP penarik ditentukan dan reaksi keseimbangan pada tumpuan/constraint dapat diekstrak.

Secara konseptual, matriks kekakuan elemen diperoleh dari integral berikut:

B menghubungkan perpindahan nodal dengan regangan; D adalah matriks konstitutif material; V_e adalah domain elemen. Untuk shell komposit, respons juga dipengaruhi orientasi ply dan integrasi section melalui ketebalan.

Penerapan pada laminat

Setiap ply memiliki sumbu material lokal. Kekakuan ortotropik didefinisikan pada sumbu material, lalu diarahkan sesuai sudut ply terhadap sumbu spesimen. Abaqus menggabungkan kontribusi tiap ply pada composite layup. Kesalahan dalam orientasi material atau arah stacking dapat mengubah hasil.

FEM adalah metode diskretisasi dan formulasi numerik. FEM bukan nama algoritma iterasi tertentu milik Abaqus.

5. Algoritma Abaqus: Newton–Raphson dan Iterative Solver

Model menggunakan Abaqus/Standard dengan prosedur Static, General, NLGEOM=ON, material elastik tipe lamina, dan automatic incrementation. Materialnya elastik, tetapi NLGEOM=ON mengaktifkan nonlinearitas geometri: perubahan konfigurasi struktur ikut diperhitungkan dalam keseimbangan. Karena itu, solusi dibangun secara inkremental dan iteratif.

Persamaan keseimbangan

Pada akhir setiap increment, Abaqus mencari perpindahan yang membuat gaya eksternal dan gaya internal seimbang. Residual keseimbangan secara umum dituliskan sebagai:

R adalah residual gaya; P_ext adalah gaya eksternal ekuivalen; P_int adalah gaya internal yang dirakit dari respons elemen; dan u adalah vektor perpindahan nodal. Perpindahan pada RP penarik diterapkan sebagai kondisi batas displacement.

Langkah Newton–Raphson dalam satu increment

1. Mulai increment. Abaqus menaikkan displacement RP-2 secara bertahap menuju +5 mm. RP-1 tetap ditahan. Perkiraan awal biasanya berasal dari keadaan increment sebelumnya.
2. Hitung respons elemen. Dari perkiraan perpindahan, Abaqus menghitung regangan dan tegangan tiap elemen/ply, lalu merakit gaya internal global.
3. Hitung residual. Abaqus membandingkan keseimbangan gaya eksternal dan internal. Residual yang belum cukup kecil berarti solusi belum seimbang.
4. Hitung koreksi Newton. Persamaan keseimbangan dilinierkan di sekitar perkiraan perpindahan saat ini:

K_T adalah matriks kekakuan tangent pada iterasi ke-k; Δu^k adalah koreksi perpindahan; dan u^k+1 adalah perkiraan perpindahan setelah dikoreksi.

5. Selesaikan persamaan koreksi. Abaqus menyelesaikan sistem linear di atas menggunakan equation solver yang dipilih. Pengaturan equation solver pada model ini adalah Iterative.
6. Perbarui konfigurasi. Abaqus menggunakan perpindahan yang telah dikoreksi, menghitung ulang respons elemen dan residual, serta memperhitungkan konfigurasi terdeformasi karena NLGEOM=ON.
7. Periksa konvergensi. Residual gaya dan besar koreksi perpindahan diperiksa terhadap toleransi konvergensi. Jika kriteria belum terpenuhi, iterasi Newton diulang pada increment yang sama.
8. Terima atau ulangi increment. Increment yang konvergen diterima. Jika solusi tidak konvergen, automatic incrementation dapat mengurangi ukuran increment dan mencoba kembali.
9. Lanjutkan sampai target tercapai. Proses increment dan iterasi berlanjut sampai displacement RP-2 mencapai +5 mm atau analisis berhenti karena masalah konvergensi.

Perbedaan Newton–Raphson dan Iterative Solver

Untuk model ini, Full Newton direkomendasikan sebagai solution technique, disertai automatic incrementation. Equation solver yang digunakan adalah Iterative. Untuk spesimen shell berukuran tugas, Direct dapat digunakan sebagai pembanding/baseline. Laporan perlu mencantumkan solver yang benar-benar aktif pada job analisis.

Automatic incrementation mengatur ukuran increment selama step. Tetap periksa ukuran increment awal, minimum, dan maksimum pada Step Editor. Periksa juga file .sta dan .msg untuk memastikan apakah increment diterima, dicoba ulang dengan ukuran lebih kecil, atau gagal konvergen.

Static, General mengabaikan efek inersia untuk analisis statik. Parameter “time” pada step ini berfungsi sebagai parameter kemajuan pembebanan, bukan otomatis waktu fisik eksperimen.

Rujukan: Dassault Systèmes, Solving Nonlinear Problems; Configuring a Static, General Procedure.

6. Verifikasi dengan Teori Klasik Laminat

Classical Laminate Theory (CLT) dapat digunakan sebagai pemeriksaan independen untuk respons membran elastik linier. CLT bukan solver Abaqus; dalam proyek ini CLT berfungsi sebagai verifikasi analitik/semianalitik pada daerah ukur yang jauh dari efek grip.

Langkah perhitungan

1. Ambil konstanta elastik tiap ply dari sumber data: E₁, E₂, G₁₂, dan ν₁₂ untuk asumsi tegangan bidang.
2. Bentuk matriks kekakuan ply Q pada sumbu material lokal.
3. Transformasikan Q menjadi Q̄ sesuai orientasi masing-masing ply.
4. Hitung matriks extensional A = Σ Q̄_k(z_k − z_k−1).
5. Gunakan resultan beban per satuan lebar dan kondisi batas yang sepadan dengan model.
6. Bandingkan regangan aksial atau gaya reaksi FEM pada daerah ukur dengan prediksi CLT.

N dan M adalah resultan gaya dan momen per satuan lebar; ε₀ adalah regangan bidang tengah; κ adalah kelengkungan; A, B, dan D adalah matriks kekakuan extensional, coupling, dan bending.

Untuk laminat simetris, B bernilai nol secara teori. Jika laminat juga balanced sehingga coupling tarik–geser dapat diabaikan, dan tepi spesimen bebas berkontraksi lateral (N_y=0), estimasi regangan aksialnya:

N_x adalah gaya tarik per satuan lebar; δ_x adalah pertambahan panjang daerah ukur; dan L adalah panjang ukur. Rumus ini mengasumsikan medan membran seragam dan sisi bebas berkontraksi lateral. Jika grip atau coupling mengubah kondisi tersebut, lokasi dan rumus perbandingan perlu disesuaikan.

7. Mesh Convergence

Mesh convergence diperiksa dengan menjalankan beberapa mesh yang semakin halus, sementara geometri, material, layup, kondisi batas, displacement, dan pengaturan solver tetap. Sumbu X grafik dapat menggunakan number of elements; sumbu Y berisi respons terpilih.

Untuk uji tarik displacement control, respons utama yang disarankan adalah gaya reaksi total pada RP penarik saat displacement mencapai 5 mm. Regangan aksial atau tegangan di area ukur dapat menjadi respons tambahan, selama lokasi dan jenis output-nya konsisten pada semua mesh.

R_i adalah respons pada mesh ke-i. Ambang konvergensi, misalnya 2%, merupakan kriteria yang dipilih untuk proyek ini, bukan aturan universal.

8. Hasil yang Perlu Dilaporkan Setelah Job Berjalan

9. Kesimpulan

Metode numerik proyek ini adalah FEM dengan Abaqus/Standard, Static, General, material elastik lamina, elemen S4R, coupling daerah grip ke reference point, dan displacement tarik +5 mm. NLGEOM=ON membuat analisis memperhitungkan nonlinearitas geometri. Abaqus menyelesaikan masalah secara inkremental; Full Newton direkomendasikan sebagai solution technique awal untuk mencapai keseimbangan pada tiap increment.

Equation solver method yang saat ini dipilih adalah Iterative. Untuk model shell spesimen ukuran tugas, Direct dapat dipakai sebagai baseline dan dibandingkan bila diperlukan. Keduanya adalah pilihan penyelesaian persamaan linear, bukan nama metode iterasi nonlinear. Keandalan hasil perlu diperiksa melalui keseimbangan reaksi, verifikasi CLT, pemeriksaan orientasi ply, status konvergensi, dan mesh convergence.

Referensi

1. Dassault Systèmes, Configuring a Static, General Procedure, dokumentasi Abaqus 2025. Buka dokumentasi.
2. Dassault Systèmes, Solving Nonlinear Problems, dokumentasi Abaqus 2025. Buka dokumentasi.
3. Dassault Systèmes, Including Nonlinearity in an Abaqus Analysis, dokumentasi Abaqus 2025. Buka dokumentasi.
4. Dassault Systèmes, Choosing a Shell Element, dokumentasi Abaqus 2025. Buka dokumentasi.
5. Dassault Systèmes, Linear Elastic Behavior, dokumentasi Abaqus 2025.