jenis element linear brick pada masalah pelat berlubang

Jenis element brick hexahedral yang linear ada beberapa macam diantaranya full integration (C3D8), reduced integration (C3D8R) dan incopatible mode (C3D8I). pada tulisan ini ditinjau masalah pelat tebal berlubang dengan jenis mesh unstructured hexahedral, modelasisasi menerapkan kondisi batas simetris. pemodelan geometry dan meshing menggunakan program bantu GMSH, solver FEA yang digunakan adalah CalculiX.

2015-08-17 05_56_22-Gmsh - C__Users_synt_Desktop_platehole2.geo

adapun pertimbangannya pemilihan elment tsb adalah thd jenis analisa yg akan ditempuh,

menurut Dhondt G., 1998 pembuat program yg saya gunakan ini,

elemen delapan nodal brick dengan full integration (C38) cukup sederhana dan dapat banyak digunakan, namun perlu dihindari untuk permasalahan lentur seperti model balok dan pelat karena kekakuan lentur untuk elemen ini cukup tinggi. keadaan ini biasa disebut shear-locking.

selian itu juga berperilaku kurang sesuai pada kondisi isochoric (tegangan bertambahan namun volume tetap), ini terjadi pada kepemilikan material dengan nilai banding poison tinggi dan keadaaan lain masalah plastisitas.

elemen delapan nodal brick dengan reduced integration  (C38R) mengurangi kekakuan lentur dan masalah shear-locking, namun disisi lain menjadikan terlalu lentur. integrasi point nya hanya satu ditengah, maka perlu penghalusan mesh pada derah konsentrasi seperti tumpuan dan beban. masalah hour-glassing (?) namun sejak program keluaran tahun 2011 sudah ditanamkan otomatisasi untuk mengendalikan ini.

elemen delapan nodal brick dengan incompatible mode (C38I) sudah menghilangkan masalah shear-locking dan volumetric locking dengan menggunakan shape function standar beserta bubble function (?). untuk jenis elemen linear C3DI ini yg direkomendasikan untuk digunakan pada banyak permasalahan, namun disisi lain tidak dapat digunakan untuk masalah dinamis karena DOF’s tambahan pada element ini tidak ada keterkaiatan hubungan dengan massa yang menyebabkan singularitas pada matrix massa.

walaupun elemen C3DI sudah baik untuk jenis lement linear, tetap direkomendasikan menggunakan element quadratic karena elment ini biasanya menghasilkan analisa yang baik.

2015-08-18 16_21_55-https___www.michigan.gov_documents_mdot_MDOT_Research_Report_R1517_248445_7.pdf

(souce: Michigan DOT’s)

Berikut ditampilkan perbedaan hasil keluaran dari berbagai lement brick linear tersebut pada permasalahan pelat tebal berlubang.

Acuan awal dengan perhitungan tangan dari buku rujukan,

2015-08-18 00_16_33-PTC Mathcad Express Prime 3.0 - C__Users_synt_Documents_platehole.mcdx

Data masukan karakteristik material,

2015-08-18 01_01_47-platehole2.inp _ SciTE

Penentuan kondisi batas/restraint (hijau fixed ux=0, biru fixed uy=0)2015-08-18 08_53_29-Calculix Graphix

Beban merata sebesar diterapkan pada face element solid.

2015-08-18 09_05_46-Calculix Graphix

Deformasi pelat dan dan vektor resultannya,

2015-08-18 10_23_29-Calculix Graphix

2015-08-18 10_18_28-Calculix Graphix

Tinjauan (1) pertama adalah element brick hexahedral yang jenis full integration (C3D8)

2015-08-17 05_41_23-platehole2.inp _ SciTE

hasil keluaran program – besarnya defleksi arah x,y (mm)

2015-08-18 09_09_59-Calculix Graphix

2015-08-18 09_10_34-ParaView 4.3.1 32-bit

2015-08-18 09_10_59-Calculix Graphix

Tegangan arah x,y, xy dan von Mises (N/mm^2)2015-08-18 09_11_15-Calculix Graphix2015-08-18 09_11_29-ParaView 4.3.1 32-bit2015-08-18 09_11_51-ParaView 4.3.1 32-bit2015-08-18 09_12_06-ParaView 4.3.1 32-bit

Tinjauan (2) pertama adalah element brick hexahedral yang jenis reduced integration (C3D8R)

2015-08-17 05_41_39-platehole2.inp _ SciTE

hasil keluaran program – besarnya defleksi arah x dan y (mm)

2015-08-18 09_53_13-Calculix Graphix

2015-08-18 09_53_41-Calculix Graphix

2015-08-18 09_53_54-platehole2x.inp - SciTE

Tegangan arah x,y, xy dan von Mises (N/mm^2)

2015-08-18 09_54_10-Calculix Graphix

2015-08-18 09_54_26-platehole2x.inp - SciTE

2015-08-18 09_54_45-Calculix Graphix

2015-08-18 09_54_57-Calculix Graphix

Tinjauan (3) pertama adalah element brick hexahedral yang jenis incompatible mode (C3D8I)

2015-08-17 05_40_53-platehole2.inp - SciTE

hasil keluaran program – besarnya defleksi arah x dan y (mm)

2015-08-18 09_59_41-Calculix Graphix

2015-08-18 10_00_16-platehole2x.inp - SciTE

2015-08-18 10_00_29-platehole2x.inp - SciTE

Tegangan arah x,y dan kriteria von Mises (N/mm^2)

2015-08-18 10_00_59-Calculix Graphix

2015-08-18 10_01_16-Calculix Graphix

2015-08-18 10_01_32-platehole2x.inp - SciTE

2015-08-18 10_01_43-platehole2x.inp - SciTE

perbedaan tersebut ditabelkan berikut:

2015-08-18 11_05_02-Untitled 1.ods - OpenOffice Calc

terlihat perbedaan masih cukup besar, memang masalah pelat berlubang dengan konsentrasi tegangan sangat sensitif terhadap kehalusan mesh pada titik tersebut. perlu dibahas tersendiri untuk tinjauan beberapa kondisi pias elemen hingga dari berbagai model tingkatan.

hal yg perlu ditinjau ulang juga mengenai perbedaan pada panjang pelat karena persamaan dari Roark tidak menunjukannya, apakah ini berlaku juga untuk pelat yg pendek seperti diatas? hal yg penting juga perlu perbandingan dengan hasil experiment.

lainnya mungkin, tambahan analisa terhadap pengaruh plastisitas material karena terlihat ada bagian yg leleh (fy>240MPa) merujuk pustaka Neuber dan ini ada hubungannya dengan sifat kelelahan (fatique) material baja.

Leave a Reply

Please log in using one of these methods to post your comment:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s