Difference between revisions of "Evi Elisa Ambarita"
(→UTS) |
(→Artikel Kelompok, OSCILLATING 1-D DYNAMIC SYSTEM) |
||
Line 322: | Line 322: | ||
Dari kedua grafik diatas, dapat dilihat bahwa metode numerik ini (grafik kuning) memiliki tingkat keakurasian yang sangat baik dengan pemverifikasian terhadap metode analitik (grafik hijau). Dengan demikian, disimpulkan bahwa a magic fix for a numerical method telah terverifikasi terhadap metode analitik yang mana dapat mewakili metode analitik. untuk menentukan nilai u dan v. | Dari kedua grafik diatas, dapat dilihat bahwa metode numerik ini (grafik kuning) memiliki tingkat keakurasian yang sangat baik dengan pemverifikasian terhadap metode analitik (grafik hijau). Dengan demikian, disimpulkan bahwa a magic fix for a numerical method telah terverifikasi terhadap metode analitik yang mana dapat mewakili metode analitik. untuk menentukan nilai u dan v. | ||
− | ==Artikel | + | ==Artikel 1, OSCILLATING 1-D DYNAMIC SYSTEM== |
===Pendahuluan=== | ===Pendahuluan=== | ||
Line 550: | Line 550: | ||
b) Hasil perhitungan numerik dan simulasi berbeda dengan metode analitik. Seperti sudah disebutkan pada materi perkuliahan sebelum ini, bahwa pada metode Forward Euler terdapat hal yang harus diperhatikan pada perkembangan amplitudonya. | b) Hasil perhitungan numerik dan simulasi berbeda dengan metode analitik. Seperti sudah disebutkan pada materi perkuliahan sebelum ini, bahwa pada metode Forward Euler terdapat hal yang harus diperhatikan pada perkembangan amplitudonya. | ||
c) Hasil metode forward euler telah terverifikasi dengan metode Artificial Neural Networks (ANN) yang dapat dilihat dari nilai R sebesar 0,99996 (mendekati 1) dari persamaan output metode ANN menggunakan software matlab. | c) Hasil metode forward euler telah terverifikasi dengan metode Artificial Neural Networks (ANN) yang dapat dilihat dari nilai R sebesar 0,99996 (mendekati 1) dari persamaan output metode ANN menggunakan software matlab. | ||
+ | |||
+ | ==Artikel 2, Finite Element Simplified == |
Revision as of 22:43, 10 May 2020
Contents
- 1 BIODATA
- 2 Pengenalan Diri
- 3 Pendahuluan Komputasi Teknik
- 4 Sinopsis Skripsi
- 5 Presentasi Sinopsis Skripsi
- 6 Quiz
- 7 Model Komputasi Teknik
- 8 Extended Abstract
- 9 Optimasi Kebutuhan Energi Manusia
- 10 9 Maret 2020
- 11 UTS
- 12 Quiz Oscillating dynamic system
- 13 Artikel 1, OSCILLATING 1-D DYNAMIC SYSTEM
- 14 Artikel 2, Finite Element Simplified
BIODATA
Nama : Evi Elisa Ambarita
NPM : 1906324044
Jurusan : Konversi energi, Teknik mesin, Universitas Indonesia
Email : evi.elisa91@ui.ac.id
Pengenalan Diri
Baseline of my knowledge and skill related to computing
I am a master student of mechanical engineering in Universitas Indonesia. I got my undergraduate degree from Marine Engineering of Universitas Indonesia in Januari 2016. For my bachelor thesis, i designed a mini container with 1 per 12 of the 12 feet container size using autodesk inventor software which applies finite element approach for stress analysis.
After got my bachelor degree, I worked for a year in a Shipyard company in Jakarta. I was in PPC department working on new ship building design using autodesk autocad 2D and 3D software. From May 2018 until August 2020, I moved to Nagasaki, Japan to work with engineering consultant company. In 2 years, I finished 3 projects about ship building, bridge construction, and main foundation of tidal turbine. The company utilized a software known as LS-DYNA whose concept is Incompressible Computational Fluid Dynamics (ICFD) . With that software, i tested the construction strength of tidal turbine foundation. Moreover, I was registered as a research student in Nagasaki University. I learned about Ansys software with Linux program to design and simulate the tidal turbine.
After worked for 3 years in total, I decided to continue my study in Universitas Indonesia. I have a big desire to construct the real tidal turbine in Indonesian's ocean, since Indonesia has great potential of ocean current to generate the electricity around 75 GW. Now for my papers and master thesis, I have been working on design development on diffuser angle of horizontal axis tidal turbine for efficiency optimization by controlling water flow passively. This utilizes CFD solver of Ansys software.
In the conclusion, I have been working on engineering designing software since 2015 until now. Hopefully I'll be able to apply all these knowledge about engineering computation on the technology development for Indonesia.
Pendahuluan Komputasi Teknik
I. Pengertian
Komputasi teknik adalah suatu kegiatan menambah, kurang, kali, bagi angka-angka yang disebut numerik dari eksak dengan menggunakan komputer. Inti berfikir dalam pemograman adalah input-proses-output. Dimana algoritma manusia dalam dipindahkan dalam bentuk algoritma computer, yang disebut dengan bahasa pemograman.
II. Tujuan Pembelajaran
1. Mampu memahami konsep-konsep dan prinsip-prinsip didalam pelajaran komputasi teknik, seperti konsep iterasi, konsep error (floating point error), dsb.
2. Mampu menerapkan pemahaman tersebut dalam menyelesaikan masalah didalam masyarakat
3. Memiliki pengenalan diri yang lebih tentang muhasabah, skill, value, dsb.
Sinopsis Skripsi
Judul: Perancangan Mini Kontainer dengan Sistem Bongkar Pasang untuk Angkutan Kapal Serbaguna
Ringkasan:
Indonesia masih memiliki banyak pelabuhan kecil dengan fasilitas yang tidak memadai layaknya pelabuhan besar. Fasilitas tersebut diantaranya crane yang mampu mengangkut kontainer 20 ft dan mobil truk kontainer. Dengan keadaan tersebut, maka pelabuhan-pelabuhan tersebut tidak dapat digunakan untuk transportasi barang dari pulau satu ke pulau lainnya. Oleh karena itu, diperlukan suatu benda dengan dimensi dan massa yang lebih kecil dibandingkan kontainer biasa. Pada skripsi ini, untuk tujuan pelayaran rakyat, penulis mendesain mini kontainer dengan dimensi dan massa 12 kali lebih kecil dibandingkan kontainer biasa, sesuai spesifikasi bak mobil pick up mitsubishi colt L300. Disamping itu, mini kontainer ini dirancang khusus dengan sistem bongkar pasang untuk mengurangi ruang dan titik gravitasi pada saat membawa mini kontainer dalam keadaan kosong. Penelitian ini dilakukan dengan metode numerik yakni mendesain dan mensimulasi kekuatan dari mini kontainer tersebut menggunakan pendekatan finite element yang terdapat pada software autodesk inventor. Berdasarkan simulasi stress analysis menggunakan software tersebut, maka mini container ini dikategorikan kuat dan mampu menjadi standar nasional untuk diaplikasikan di pelabuhan-pelabuhan kecil.
Korelasi terkait sinopsis skripsi dengan komputasi teknik
Desain mini kontainer pada skripsi ini dirancang menggunakan software autodesk inventor dengan pendekatan finite element pada stress analysis.
Presentasi Sinopsis Skripsi
Quiz
Model Komputasi Teknik
Pengerjaan skripsi saya menggunakan perangkat lunak inventor yang didasarkan pada Finite Element Method (FEM) untuk stress analysis. FEM adalah model perhitungan numerik yang menemukan perkiraan solusi untuk masalah persamaan differensial parsial dan integrasi lainnya yang dihasilkan dari hasil diskritisasi benda kontinum. Apa yang dilakukan FEM adalah membagi sebuah masalah besar menjadi bagian-bagian yang lebih kecil dan sederhana. Dalam Autodesk Inventor Stress Analysis, itu membutuhkan struktur yang complex dan mengubahnya menjadi bagian-bagian lebih kecil (meshing process) dan kemudian menyelesaikan persamaan dengan sistem persamaan yang memiliki input yang berbeda-beda, seperti constraints, materials, and loads.
In stress analysis, the weak form is called the principle of virtual work.
dimana,
w = deflection
u = displacement
A = area
E = young’s modulus
b = axial loading
Untuk persamaan lainnya terkait Finite Element Analysis dapat dilihat pada tabel dibawah ini,
source: https: //www. simscale.com/docs/content/simwiki/fea/whatisfea.html
Extended Abstract
Nama: Evi Elisa Ambarita Judul: Mini Container Design with Knock-Down System for Multipurpose Vessel Transportation.
In fact, Indonesia still has a lot of small ports whose facilities are not as compatible as big ports’ such as a container truck and crane with small power which is not able to lift a 20 ft container. Accordingly, those harbours cannot be categorized as a cargo port. Therefore, this study proposed a mini container design with dimension and mass of 12 times smaller than the regular container so that it can be carried only by a pick-up car and forklift. Moreover, this mini container was especially designed with known-down system to minimize the space and decrease the centre or gravity while toting it in empty condition. The inventor software was utilized to design and simulate mini container, also to test its material strength well known as stress analysis. The software applies finite element method (FEM) to subdivide a large system into smaller and simpler parts. To be able to run simulations, a mesh, consisting of up to millions of small elements that together form the shape of the structure needs to be created. Calculations are made for every single element. Combining the individual results gives the result of the structure. Afterwards, the simulation was conducted by inputting the stress distribution of force value and the fixed constraint. The result of this numerical method was represented the von mises stress, 1st and 3rd principal stress to obtain the occurred displacement when there was bending on the object. Based on the simulation and calculation, this mini container was concluded to meet the ISO standard about size and material of container to carry cargo for multipurpose vessel transportation.
Optimasi Kebutuhan Energi Manusia
Sebagai unit dari energi, kalori mengacu pada energi yang diperoleh dari makanan dan minuman yang kita konsumsi setiap hari. Kalori sangat penting bagi tubuh manusia. Tanpa kalori, sel-sel dalam tubuh akan mati, dan organ-organ vital seperti jantung dan paru-paru tidak akan mampu melaksanakan fungsi dasarnya.
Di Indonesia, berdasarkan rekomendasi Angka Kecukupan Gizi (AKG) dari Kementerian Kesehatan RI, rata-rata kebutuhan untuk pria usia 30 – 49 tahun adalah 2625 kkal per hari. Sedangkan perempuan usia 30 – 49 tahun adalah 2150 kkal per hari, demikian seperti dikutip dari laman depkes.go,id.
Menghitung Kalori
Terdapat dua jenis kalori, yaitu kalori kecil yang ditulis dalam satuan “kal” dan kalori besar atau “kilokalori” (kkal), dengan perhitungan 1 kkal sama dengan 1.000 kal. Cara menghitung kalori yang kita butuhkan menurut P2PTM Kemenkes RI sangat sederhana, karena hanya berdasarkan jenis kelamin dan tinggi badan. Sebelum menghitung, kita harus mengetahui terlebih dahulu tinggi badan (TB) dalam sentimenter dan Berat Badan Ideal (BBI) dengan rumus sebagai berikut:
BBI = (TB-100) – (10% x (TB – 100))
Selanjutnya, kita dapat hitung Kebutuhan Kalori Basal (KKB) atau Basal Metabolic Rate (BMR). KKB merupakan kebutuhan kalori yang dibutuhkan oleh tubuh untuk metabolisme basal, yakni metabolisme yang wajib dilakukan mahluk hidup walaupun tidak membutuhkan energi. Ini berarti metabolisme tubuh akan tetap berjalan meskipun dalam keadaan tidur atau tidak melakukan apa-apa. Berikut adalah rumus menghitung angka KKB:
KKB Laki-laki = 30 kkal x BBI
KKB Perempuan = 25 kkal x BBI
Menghitung Kebutuhan Kalori Total (KKT)
Kebutuhan kalori total adalah jumlah kebutuhan kalori tubuh ditambah dengan jumlah kalori saat melakukan aktivitas fisik. Kita mengenal tiga jenis aktivitas, yaitu
1. Aktivitas ringan seperti membaca (10%), menyetir mobil (10%), kerja kantoran (10%), mengajar (20%), berjalan (20%).
2. Aktivitas sedang: kerja rumah tangga (20%), jalan cepat (30%), bersepeda (30%).
3. Aktivitas berat: aerobik (40%), mendaki (40%), dan jogging (40%)
Rumus KKT = KKB + Aktivitas Fisik - Faktor Koreksi
Faktor koreksi adalah sebagai berikut:
Usia 40 - 59 tahun, nilai koreksinya minus 5%
Usia 60 - 69 tahun, nilai koreksinya minus 10%
Usia >70 tahun, nilai koreksinya minus 20%
Resource: https:// www.anlene.com/id/ms/berapa-banyak-kalori-yang-anda-butuhkan-setiap-hari.html
9 Maret 2020
Berikut tahapan dalam menganalisa suatu masalah:
1. Initial thinking (analisis) untuk mengidentifikasi masalah (objektif) yang akan diselesaikan, dapat berupa masalah terstruktur yakni terkait SOP, dan unstructured issue terkait modifikasi.
2. Mengembangkan rumusan (permodelan matematis) yang mengandung asumsi (tergantung pengetahuan, ilmu, dan pengalaman kita), dimana asumsi sama dengan batasan masalah.
3. Pengumpulan data yang dibutuhkan untuk menyelesaikan permodelan matematis diatas.
4. Simulasi = menjalankan atau mengeksekusi model (constraint) yang kita buat dengan suatu variabel bebas
5. Verifikasi = resolve the equation right, menguji apakah model yang dilakukan tidak ada kesalahan numerik, contohnya kesalahan ngitung, pada proses meshing
6. Validasi = resolve the right equation, menguji kebenaran / keaktualan dari hasil simulasi tadi
7. Result dan discussion
8. Recommendation
UTS
Video presentasi
Video presentasi hasil belajar terkait pengetahuan (konsep/teori) dan keterampilan (menggunakan komputasi teknik)
Laporan tugas optimasi
Laporan hasil tugas optimasi kebutuhan energi manusia.
Berdasarkan hasil diskusi kelas pada ( http://air.eng.ui.ac.id/index.php?title=Optimasi_Sistem_Energi_pada_Manusia ) ini, berikut optimasi kebutuhan energi saya pribadi selama seminggu (March 16th - 22nd),
Kebutuhan kalori yang saya butuhkan dalam 1 minggu:
Kebutuhan listrik yang saya butuhkan dalam 1 minggu:
Sehingga diperoleh total energi dan total biaya yang saya butuhkan dalam seminggu:
Dan menghasilkan grafik sebagai berikut:
Draft paper
Design Development of Strength Testing on Mini Container Construction with Knock-Down System for Multipurpose Vessel
Background
Indonesia has thousands of ports including small ports. A lot of small ports cannot be categorized as a cargo port since there is no standard facilities to be a cargo port, such as crane and container car. Whereas the citizens dwelling there still need food storage even just a little. Therefore, in order to support their needs, this study proposed a mini container design with dimension and mass of 12 times smaller than the regular container so that it can be carried only by a pick-up car and forklift. Moreover, this mini container was especially designed with known-down system to minimize the space and decrease the center or gravity while toting it in empty condition.
The design has been presented on author's bachelor thesis in 2016. Autodesk inventor software was utilized to simulate and test the construction of mini container for stress analysis. It was concluded that the container is safe to be used in small cargo port since the safety factor was above 1, the deflection was under 3 mm, and the von mises stress of mini container was less than the material yield stress of carbon steel.
Due to time constraints, the strength testing can only be done using inventor software for stress analysis. Even though, there are still many aspects that need to be investigated further, such as temperature, pressure, corrosion, etc. As known, factors as wind flow, water density, food types, can affect the material utilized for the mini container. Therefore, the current study submits a validation of material strength and feasibility tests of mini container construction using finite element method (FEM) on ansys software, so that it can be made as an international standard and commercialized.
Objectives
1. Further examine the strength material of mini container construction
2. Validate the simulation results of stress and strength analysis on mini container construction using FEM on ansys software
3. Define the best material for mini container
Problem Formulation
1. What other aspects influence the material strength of mini container construction can be shown in FEM on ansys software?
2. What kind of material is suitable for this mini container?
3. What kind of storage is suitable for this mini container?
Methodology
This study will utilize numerical method. Since the design has been performed on previous study shown in figure 1, so this study will straight away use ansys software to simulate and test the material applying FEM. The trial and error will be done to receive the best results of material by adjusting the carbon content of the material and a number of load the mini container can be carried.
Hypothesis
The validation will support the statement of previous study which is to utilize carbon steel AISI 1060 with 0.6% of carbon with these properties.
Quiz Oscillating dynamic system
Memverifikasi sebuah sistem dinamik berosilasi tanpa gesekan (fig. 4.15) dengan menggunakan metode numerik dan metode analitik. Selain itu, terdapat 2 metode numerik yang akan digunakan, yakni standard numerical solution (persamaan 4.47 dan 4.48) dan a magic fix of a numerical method (persamaan 4.49 dan 4.50).
Gambar diatas menjelaskan suatu benda dengan massa m melekat/dikaitkan pada pegas dan bergerak sepanjang garis tanpa gesekan. Ketika pegas diregangkan (atau dikompresi), gaya pegas menarik (atau mendorong) body (penampang m) kembali dan bekerja "melawan" gerakan, misalkan x (t) adalah posisi body pada sumbu x, dimana body bergerak. Pegas tidak direntangkan ketika x= 0, sehingga gaya adalah nol, dan x= 0 karenanya posisi keseimbangan bodi. Gaya pegas adalah -kx, dimana k adalah konstanta yang diukur. Kami berasumsi bahwa tidak ada gaya lain (mis., Tidak ada gesekan). Hukum kedua Newton menyatakan F=ma kemudian memiliki F=-kx dan a=x",
-kx = mx" (4.41)
sehingga dapat dituliskan
x" + ω^2 x = 0 (4.42)
Persamaan (4.42) adalah persamaan diferensial orde kedua, dan oleh karena itu kita memerlukan dua kondisi awal, satu pada posisi x(0) dan satu pada kecepatan x’(0).
x (0) = 1,5 dan x’(0) = 0
Untuk mengendalikan turunan orde dua menggunakan metode numerik, cara yang dapat kita tuliskan pada persamaan 4.42 sebagai sebuah sistem orde pertama dari 2 persamaan turunan. Kita memperkenalkan u=x dan v=x'=u' sebagai 2 fungsi baru yang tidak diketahui. Dua persamaan yang sesuai muncul dari definisi v=u' dan persamaan asal (4.42):
u' = v (4.43)
v' = -ω²u (4.44)
maka u"=v' untuk menghilangkan turunan orde kedua dari hukum kedua newton.
Dengan menggunakan metode numerik pertama, yakni standard numerical solution dengan persamaan-persamaan berikut:
u^(n+1) = u^n + ∆t v^n (4.47)
v^(n+1) = v^n - ∆t ω^2 u^n (4.48)
Selain itu, terdapat juga metode numerik kedua yang anggap lebih akurat, disebut a magic fix of the numerical method dengan mengganti u^n untuk persamaan kedua menjadi u^(n+1) dari persamaan pertama, sehingga:
u^(n+1) = u^n + ∆t v^n (4.49)
v^(n+1) = v^n - ∆t ω^2 u^(n+1) (4.50)
dimana diketahui ω = 2,5 rad/s dan ∆t = 0,05 s, maka dilakukan perhitungan sebagai berikut.
Selanjutnya dilakukan verifikasi antara setiap metode numerik diatas terhadap metode analitik, dimana persamaan-persamaan untuk metode analitik adalah sebagai berikut:
u = x cos ωt
v = -x ω cos ωt
yang perhitungan detail nya dapat juga dilihat pada lampiran (ms.excel).
Sehingga diperoleh hasil yang merupakan grafik sebagai berikut:
Dari kedua grafik diatas dapat dilihat jelas bahwa metode numerik yang standard dengan menggunakan persamaan 4.47 dan 4.48 jauh dari kata akurat, dikarenakan hasil yang didapat dari metode ini (grafik biru) sangat berbeda dengan hasil dari metode analitik (grafik abu-abu). Oleh karena itu, metode numerik yang disebut a magic fix for a numerical method perlu dilakukan, dengan hasil sebagai berikut:
Dari kedua grafik diatas, dapat dilihat bahwa metode numerik ini (grafik kuning) memiliki tingkat keakurasian yang sangat baik dengan pemverifikasian terhadap metode analitik (grafik hijau). Dengan demikian, disimpulkan bahwa a magic fix for a numerical method telah terverifikasi terhadap metode analitik yang mana dapat mewakili metode analitik. untuk menentukan nilai u dan v.
Artikel 1, OSCILLATING 1-D DYNAMIC SYSTEM
Pendahuluan
Artikel ini akan membahas penyelesaian permasalahan oscillating 1-D dynamic system pada system damping sesuai gambar di bawah. Pada system di bawah terdapat sebuah massa yang dihubungkan dengan 3 spring, dan 1 dashpot. System tersebut fix pada 2 boundary sisi kiri dan kanan.
Penyelesaian dari permasalahan tersebut akan dilakukan dengan 3 metode, yang kemudian akan dibandingkan hasil dari perhitungan pada tiap metode tersebut. Metode yang digunakan yaitu: analitikal, numerik, dan simulasi. Pada metode analitikal dilakukan tinjauan gaya aksi – reaksi pada free body diagram (FBD) yang muncul dari eksitasi yang diberikan. Kemudian untuk penyelesaian metode numerik akan dilakukan iterasi dengan bantuan software Microsoft Excel menggunakan Forward Euler. Sedangkan untuk proses simulasi akan dilakukan dengan bantuan software MATLAB.
Metode
Metode Analitik
Free Body Diagram (FBD):
• Tinjauan massa:
• Tinjauan titik: (diantara c1, k2, dan k3)
Diketahui:
Maka persamaan 1 dan 2 dapat dituliskan kembali menjadi;
Persamaan 1
Persamaan 2
Penyelesaian
Asumsi
Karena pada persamaan 5 masih terdapat ̇ y maka subtitusi persamaan 6 ke persamaan 5, sehingga didapatkan nilai ̇ z sebagai berikut;
Penyelesaian menggunakan matrix;
Maka didapatkan nilai λ sebesar
Sehingga didapatkan persamaan
Jika diasumsikan nilai C1,C2,dan C3 adalah 1, 2 dan 3 maka;
Nilai x untuk t = 0 s, t = 1 s, dan t = 10 s adalah
Mencari nilai yh dengan mensubtitusi nilai xh ke dalam persamaan ̇ z=-x+3y, di mana nilai ̇ z = ̈ x.
Mencari nilai ̇x;
Nilai ̇ x untuk t = 0 s, t = 1 s, dan t = 10 s adalah
Mencari nilai ̈ x;
Persamaan ̈ x di atas dapat disederhanakan menjadi;
Nilai ̈ x untuk t = 0 s, t = 1 s, dan t = 10 s adalah
Mencari nilai y;
Persamaan yh dapat disederhanakan menjadi;
Nilai y untuk t = 0 s, t = 1 s, dan t = 10 s adalah
Mencari nilai ̇ y;
Nilai ̇ y ketika t = 0 s, t = 1 s, dan t = 10 s adalah
Metode Numerik
Metode ini dilakukan dengan menggunakan software Microsoft Excel. Parameter input yang harus diketahui yaitu x pada t = 0 s, y pada t = 0 s, dan selisih penambahan waktu (∆t). Parameter input x dan y ini dapat juga disebut dengan initial condition. Lalu kemudian dilakukan iterasi sesuai dengan ∆t sampai waktu yang ditentukan. Nilai parameter input tersebut didapat dengan cara analitis seperti sudah disebut di atas, yaitu x((t=0))=4 m; y((t=0))=5.9 m; dan selisih waktu ∆t=0.1 second.
Seperti sudah disebutkan sebelumnya bahwa penyelesaiannya menggunakan Forward Euler, yaitu sebagai berikut:
Simulasi
Simulasi dilakukan dengan bantuan software MATLAB. Adapun agar simulasi dapat berjalan, maka perlu diberikan input coding pada software MATLAB yang isinya sebagai berikut.
% Calculates the position, velocity, and acceleration as a function of time % of a system of carts connected by springs and dashpots. Euler's Method is % used to solve the equations of motion numerically. clear all; close all; clc; tic
% Problem parameters k1=1; % cart 1 spring constant (N/m) k2=2; % cart 1 spring constant (N/m) b1=1; % cart 1 viscous damping coefficient (kg/s) m1=1; % cart 1 mass (kg) k3=3; % cart 1 spring constant (N/m) x10=4; % cart 1 initial position (m) x20=5.9; % y movement initial position (m) v10=-2.135; % cart 1 initial velocity (m/s) v20=-35.36; % cart 2 initial velocity (m/s)
% Set time step stuff simTime=60; % simulation time (s) tStep=0.1; % simulation time step iterations=simTime/tStep; t=0:iterations;
% Pre-allocate variables for speed and add initial conditions x1=zeros(iterations,1); x1(1,:)=x10; x2=zeros(iterations,1); x2(1,:)=x20; v1=zeros(iterations,1); v1(1,:)=v10; v2=zeros(iterations,1); v2(1,:)=v20; a1=zeros(iterations,1); a1(1,:)=((0-(k1*x10)-b1*(v10+v20)-k2*(x20+x10))/m1);
% Solve the ODE's with Euler's Method for n=2:(iterations+1) x1(n,:)=x1(n-1,:)+v1(n-1,:)*tStep; % cart 1 position x2(n,:)=x2(n-1,:)+v2(n-1,:)*tStep; % cart 2 position v1(n,:)=v1(n-1,:)+a1(n-1,:)*tStep; % cart 1 velocity v2(n,:)=-(2*x1(n,:)+5*x2(n,:)+v1(n,:)) % Find cart accelerations a1(n,:)=((0-k1*x1(n,:))-b1*(v1(n,:)+v2(n,:))-k2*(x2(n,:)+x1(n,:))/m1); end
% Plot results subplot(3,1,1) hold on; plot(t',x1,'r') plot(t',x2,'m') ylabel('Position (m)') title('Position, Velocity, & Acceleration as a Function of Time') legend('Pergerakan ke arah X','Pergerakan ke arah Y') subplot(3,1,2) hold on; plot(t',v1,'b') plot(t',v2,'c') ylabel('Velocity (m/s)') legend('Pergerakan ke arah X','Pergerakan ke arah Y') subplot(3,1,3) hold on; plot(t',a1,'g') ylabel('Acceleration (m/s^2)') xlabel('time (1 = 0.01 detik)') legend('Pergerakan ke arah X','Pergerakan ke arah Y')
toc
Hasil
Gambar 2. Grafik velocity (dx/dt) perbandingan analitik, numerik, dan simulasi
Gambar 3. Grafik acceleration (d^2x/dt^2) perbandingan analitik, numerik, dan simulasi
Gambar 4. Grafik velocity (dy/dt) perbandingan analitik, numerik, dan simulasi
Selanjutnya, untuk memverifikasi hasil numerik menggunakan metode euler yang dijelaskan diatas, digunakan metode Artificial Neural Networks (ANN) pada perangkat lunak matlab. Berikut sedikit pengenalan tentang ANN yang dapat kita temui dari beberapa sumber di internet. ANN adalah beberapa lapisan jaringan saraf (neural network) yang terhubung secara keseluruhan yang dapat dilihat dari gambar dibawah ini. Setiap node pada setiap lapisan terhubung ke setiap node dari lapisan selanjutnya. Jaringan dibuat lebih dalam dengan menambah jumlah lapisan tersembunyi (hidden layer).
Jika salah satu hidden layer atau node output diperbesar, maka akan ditemukan gambar dibawah ini.
Node tertentu mengambil jumlah bobot dari inputnya, dan meneruskannya melalui fungsi aktivasi non-linear. Ini merupakan output dari node yang kemudian menjadi input dari node yang lain pada layer selanjutnya. Sinyal mengalir dari kiri ke kanan, dan output akhir dihitung dengan cara melakukan prosedur ini untuk semua node. Melatih jaringan saraf berarti mempelajari bobot yang terkait dengan semua edge.
Metode ANN ini digunakan untuk melatih data-data yang didapat dari numerik, dimana nilai posisi pada sumbu x dan y digunakan sebagai nilai input, dan nilai kecepatan pada sumbu x (u) dan sumbu y (v) yang merupakan hasil diferensial dari posisi digunakan sebagai nilai target. Data yang di training sebanyak 400 data.
Dengan demikian, diperoleh hasil yang dapat dilihat pada gambar diatas, dengan nilai R sebesar 0,99996 yang berarti data yang didapat dari metode Euler terverifikasi secara baik dengan persamaan yang dihasilkan dari metode ANN pada matlab.
Kesimpulan
Dari hasil yang sudah didapat untuk metode analitik, numerik, dan simulasi berikut beberapa hal yang dapat disimpulkan: a) Hasil perhitungan numerik terhadap simulasi memiliki nilai yang berimpit. Hal ini dikarenakan, pada dasarnya simulasi yang dilakukan dengan MATLAB, menggunakan cara yang sama dengan numerik, hanya saja pada numerik persamaan-persamaannya harus dirubah menjadi 3 persamaan dengan orde 1. Sedangan pada simulasi, persamaan dari FBD langsung menjadi input pada MATLAB. b) Hasil perhitungan numerik dan simulasi berbeda dengan metode analitik. Seperti sudah disebutkan pada materi perkuliahan sebelum ini, bahwa pada metode Forward Euler terdapat hal yang harus diperhatikan pada perkembangan amplitudonya. c) Hasil metode forward euler telah terverifikasi dengan metode Artificial Neural Networks (ANN) yang dapat dilihat dari nilai R sebesar 0,99996 (mendekati 1) dari persamaan output metode ANN menggunakan software matlab.