Minggu, 18 Desember 2016

balok lengkung - perbandingan hasil FE berbagai jenis elemen

sebelumnya sya menulis untuk permasalahan pada keadaa balok lurus dan pendek,sehingga dominasi beban adalah lentur dan geser. pada kesempatan ini sya membuat perbandingan pada masalah balok lengkung.

2016-12-19-00_40_35

hasil dengan menggunakan element HEX20 - quadratic (structured hexahedral) sebagai acauan akurasi dari hasil elemen lainnya.

2016-12-19-00_39_14

2016-12-19-00_39_27

2016-12-19-00_39_37

2016-12-19-00_39_49

hasil dengan menggunakan elemen HEX8 - linear (structured hexahedral)

2016-12-19-00_34_55

2016-12-19-00_35_09

2016-12-19-00_35_23

2016-12-19-00_35_33

hasil dengan menggunakan elemen TET4 - linear,

2016-12-19-00_29_30

2016-12-19-00_30_07

2016-12-19-00_30_24

2016-12-19-00_30_39

hasil menggunakan elemen wedge linear,

model (1)

2016-12-19-01_04_00

2016-12-19-01_04_14

2016-12-19-01_04_28

2016-12-19-01_04_39

hasil menggunakan elemen wedge linear,

model (2)

2016-12-19-02_01_45

2016-12-19-02_02_01

2016-12-19-02_02_12

2016-12-19-02_02_27

mesh refinenement,

2016-12-19-02_11_06

2016-12-19-02_11_22

2016-12-19-02_11_39

2016-12-19-02_11_51

hasil menggunakan element unstructured HEX8 dan beberapa Wedge, didapat dengan membuat meshing 2D quad-dominate lalau dilakukan extrude

2016-12-19-01_07_09

 

2016-12-19-01_13_09

2016-12-19-01_13_20

2016-12-19-01_13_32

hasil dengan menggunakan HEX8 linear unstructured,

model (1)

2016-12-19-01_42_42

2016-12-19-01_43_20

2016-12-19-01_43_31

2016-12-19-01_43_43

model (2)

2016-12-19-12_12_57

2016-12-19-12_13_12

2016-12-19-12_13_25

2016-12-19-12_13_38

tabulasi hasil analisa berbagai jenis elemen tersebut,

.. to be add

  • perbandingan lainnnya pada keadaan permasalahan tegangan bidang dan pelat lentur, mungkin juga model yg agak rumit dgn konsentrasi tegangan.

  • diulang model dan permasalahan yg sama namun menggunakan jenis elemen quadratic.

  • grafik studi konvergensi dgn meshing penambahan jumlah element.

  • resume kesimpulan, number of nodes/equations dan resource dari computational times


 

 

 

Rabu, 14 Desember 2016

balok lurus - perbandingan hasil FE berbagai jenis elemen

pada penggunaan alat bantu elemene hingga untuk simulasi mekanika, pengguna diberikan pilihan pada jenis element diantaranya adalah sbb berikut:

  • tetrahedral

  • hexahedral

  • wedge/prism

  • pyramid


2016-12-15-02_40_56-tetrahedron

2016-12-15-02_40_56-hexahedron

2016-12-15-02_40_56-wedge

2016-12-15-02_40_56-pyramid

LISA-FET element type (sumber: Waluyo, 2010)


pemilihan jenis element tersebut disebabkan karena masalah atau kendala hardware resource, waktu penyelesaian, dan otomatisasi meshing. banyak software komersil atau free yg beredar, untuk kasus linear static hal yg membedakannya adalah kecepatan solver dan otomatisasi meshing. idealnya adalah menggunakan jenis element hexahedral linear atau quadratic, namun hal ini mempunyai kekurangan pada saat tahap meshing yg biasanya masih manual dgn partisi. meshing jenis element tetrahedral sudah banyak fully otomatis, sedangkan yg meshing fully automatis hexanedral masih jarang kalaupn ada jenisnya masih unstructured yg pelu peninjauan terhadapmesh quality. sisi lain dari keurangan menggunakan jenis element tetrahedral adalah masalah computational cost (hardware & times). belakangan banyak dikembangkan automatisasi meshing ntuk element kombinasi yg mempunyai dominasi besar jenis hexahedral, hal ini mengimbangi antara masalah computational cost dan meshing.

pada kesempatan ini sya membandingkan hasilberbagai element tersebut pada berbagai permasalahan sederhana (dominasi:lentur, geser, normal & torsi)

case (1) balok kantilever akibat beban merata dan berat sendiri

2016-12-15-03_25_07

 

sebagai acuan dibandingkan dgn jenis element hexahedral quadratic, yg mana jenis element tsb adalah paling akurat.

2016-12-15-03_23_42

2016-12-15-03_23_59

2016-12-15-03_24_14

2016-12-15-03_24_31

 

 

2016-12-15-02_26_03-testhexadom1-geo-scite

Modulus elastisitas material E=25000 210000MPa,  Berat jenis 2400kg/m3, Beban permukaan 0.5N/mm2

2016-12-15-03_01_09-gmsh-d__ccxstudy_testhexadom1a-geo

hasil element jenis tetrahedral,

2016-12-15-03_04_16

2016-12-15-03_05_54

2016-12-15-03_06_56

2016-12-15-03_07_10

dengan mesh refinement,

2016-12-15-07_06_48

 

2016-12-15-02_26_30-gmsh-d__ccxstudy_testhexadom1-geo

hasil analisa element wedge/prism linear,

2016-12-15-03_08_45

 

2016-12-15-02_34_21

2016-12-15-02_35_36

2016-12-15-02_35_06

dengan mesh refinement,

2016-12-15-06_20_57

2016-12-15-03_11_30-gmsh-d__ccxstudy_testhexadom1a-geo

hasil analisa element hexahedral unstructured linear,

2016-12-15-07_02_24

2016-12-15-07_03_07

2016-12-15-07_03_332016-12-15 07_03_20-.png

 

 

 

 

dengan mesh refinement,

2016-12-15-07_00_09

 

 

 

sedangkan hasil analisa element jenis hexahedral structured linear,

2016-12-15-03_30_42

2016-12-15-03_31_11

2016-12-15-03_31_26

2016-12-15-03_31_40

dengan mesh refinement,

2016-12-15-06_27_32

tabulasi perbandingan hasil diatas adalah sebagai berikut, terlihat element linear jenis hexahedral lebih baik dan dapat cepat mencapai konvergensi. element linear tetrahedral tidak akurat, sedangkan element wedge sudah cukup baik. jika menggnakan element kombinasi hexa-dominate sebaiknya mesh kategori coarse sampai medium namun dengan jenis quadratic, hal ini untuk mereduksi computational effort sehingga computer PC home biasa dapat digunakan tanpa komputer high end dgn cloud computing.

2016-12-15-07_51_59-hexdom-ods-openoffice-calc

hasil meshing jenis element hexa-dominate,

2016-12-15-06_56_58-gmsh-d__ccxstudy_testhexadom2-geo

object terlalu sederhana menjadikan meshing di GMSH dgn metode hexa-dominate tidak berhasil, atau sya ada yg terlewat (?)

gambar cukup banyak, nanti untuk kasus lain kondisi torsi & normal dibuat postingan terpisah. ditambahkan juga computational effort in milisecond