metode partisi hexahedral mesh di SALOME

prinsip dari metode partisi untuk menghasilkan jenis mesh yg full hexahedral terstruktur didasarkan dari pembagian suatu object menjadikan kubus pesegi yg terbentuk dari 6 (enam) bidang permukaan. bidang tersebut dapat berupa kelengkungan atau curvature sembarang, sedangkan object pemisah/partisi terbuat dari suatu bidang segi-empat yg dapat bidang/plane atau lengkung/curved.

pembagi sisi tebal didasarkan dari penentuan sub-mesh pada garis tepian (edge) yg didapat dari pemisahan secara explode dari main object solid. hal tersebut ditempuh ntuk menghasilkan rasio mesh yg cenderung merata uniform antara sisi tebal dan sisi lebarnya, yg mana banyak merekomendasikan nilai rasio tersebt tidak lebih dari 10.

berikut adalah beberapa contoh sederhana suatu object yg dilakukan meshing dgn metode partisi pada program SALOME-CAD/CAE

contoh (1) pelat berlubang









adapun langkah dari pembuatannya diperlihatkan pada Object Browser berikut,



langkah yg sya tempuh, secara urut adalah:

• object kubus (Box) dengan penentuan dimensi lebar dan ketebalan


• object tabung solid (Cylinder) dengan penentuan radius dan ketinggian


• operasi Boolean object solid dengan cara Cut object kubus terhadap tabung solid


• penentuan point/vertex sebagai titik bantu pembuatan garis


• pembuatan bidang permukaan (plane) dengan cara extrude  object garis


• lakukan partisi object solid pelat berlubang dengan bidang permukaan bantu (plane) yg telah dibuat


• meshing algorithma dgn metode Hexahedron (i,j,k) penentuan edge pembagi uniform untuk global mesh


• penentuan local mesh pada garis tepi (edge) group dari garis ketebalan

 

contoh (2) silinder









contoh (3) cylinder with cap/head









contoh (4) segitiga prismatic (wedge)









contoh (5) pelat berlubang dgn chamfer









contoh (6) lug









contoh (7) embossed & fillet of circular plate - quarter models











contoh (8) curved plate with holes











 

contoh (9) bracket







karena sifat dari metode partisi adalah manual, maka perlu beberapa cara/strategi yg diterapkan oleh pengguna (user) seperti local mesh dan letak bidang pembaginya untuk mencapai hasil kwalitas mesh yg baik. sedangkan kriteria dari kwalitas mesh hexahedral element jenis Aspect Ratio 3D dijelaskan dalam dokumentasi sebagai berikut,



contoh (10) rod connected to ring with chamfer





saat dilakukan meshing hexahedral ada yg tertingga bada sub-bagian pertemuan  padahal model adalah symetris, ini mungkin dikarenakan akurasi akibat penentuan sya (user) yg kurang teliti (?) atau perlu dicari lain penyebabnya.



 



kesalahan diatas sya cari lagi ini diakibatkan oleh alat bantu object partisi yg berupa banyak bidang, setelah diperbaiki dan disatukan berhasil dalam meshing hexahedral.





untuk project yg cukup serius proses meshing perlu yg fully atomatic agar dapat mempercepat waktu. interface dengan program automatic hexahedral meshing disediakan dengan integrasi program meshing commercial MeshGems dari DISTENE, namun pengguna perlu license untuk memiliki dan menggunakannya.



sedangkan meshing automatic yg sudah tersedia di SALOME hanya tetrahedral algorithma Netgen, mengapa sya masih melakukan manual meshing dgn partisi padahal sudah ada automatic tetrahedral meshing? bagi sya adalah masalah computational cost (hardware & times) due to large numbers of nodes/DoF's.

(sumber: SALOME Platform, 2015)

(sumber: DISTENE, 2013)

hal yg memungkinkan bagi sya yg hanya perorangan dalam mereduksi computational cost adalah dgn automatic Hexa-Dominate meshing dari GMSH seperti yg sdah sya posting sebelumnya, namun ini perlu study perbandingan sebelumnya terhadap hasil solver FE. nanti dilain waktu dan kesempatan sya akan menulisnya.

ada kemampuan advanced dari program SALOME-CAD/CAE yaitu parametric modeling & meshing dengan tanpa menggunakan API atau Scripting Python melainkan cukup hanya dengan NoteBook sebagai nilai masukan dan definisi parameter. keseluruhan object masih dalam kendali parameter tersebut, sehingga akan dilakukan updates secara otomatis terhadap model dan hasil meshing. Perhatian khusus diperlukan pada penentuan parameter tersebut terhadap ketergantungan (dependency) object lain pembentuknya.



kemampuan tersebut akan sangat membantu pengguna, dikarenakan tidak perlunya re-modeling dan re-meshing suatu object yg sedang di analisa/review untuk parametric study. seperti contoh pada study perilaku sambungan baja, untuk parameter perubahan: tebal pelat, lubang baut, diameter baut, kepala baut/mur, dimensi profil WF/HB/L/UNP beserta nilai takikan (fillet radius) dan lainnya.



NoteBook pada SALOME-CAD/CAE adalah fetures yg sangat berguna karena nilai suatu parameter dapat bergantung kepada nilai parameter lain yg sudah ditentukan sebelumnya sehingga dapat membuat object sembarang dgn cara fully parametric modeling & meshing.

 

Read More

perbedaan hasil FE analisa berbagai jenis element pada balok pre-twisted

[draft]

contoh kasus analisa menggunakan element shell biasa/standard yg umumnya banyak digunakan program FE general dibandingkan terhadap hasilnya dengan element shell continuum 3D yg hanya tersedia pada beberapa software FE spesifik seperti Abaqus/CalculiX,



diketahui dimensi balok:

• panjang, L=12.00 m


• lebar, b=1.10 m


• tebal, t=0.32 m

material properties:

• Modulus elastisitas, E=29000 kPa


• Angka banding poison, v=0.22

geometri balok dilakukan pre-twisted sebesar 90deg dengan perubahan linear tiap segment elemennya. beban yg diterapkan berupa point load sebesar 1.0 N pada arah sumbu-Y global dan akan dibandingkan hasil analisa defleksi arah tersebut, study konvergensi  menggunakan berbagai element serta kehalusan mesh:

• conventional 1D beam element, both Euler and Timoshenko beams


• expanded beam 3D continuum


• linear shell element (quad 4), both for Mindlin-Reissner thick and Kirchof-Love thin


• quadratic shell element (quad 8), both for Mindlin-Reissner thick and Kirchof-Love thin


• continuum 3D shell element


• brick/solid element, linear & quadratic

seperti yg sudah banyak diteliti, penggunaan conventional shell element gagal mencapai konvergensi pada kasus diatas walau digunakan refined dengan mesh yg sangat halus.

(source: LISA-FET/Kemp V.,2013)

tampilan hasil analisa menggunakan jenis element quadratic shell (quad 8) Mindlin-Reissner thick plate theory memperhitungkan deformasi geser.

(16 element, 69 nodes)

(64 element, 223 nodes)

berikut perbandingannya dengan hasil analisa menggunakan element shell continuum (S8),

(16 element, 165 nodes)

(64 element, 551 nodes)

terlihat penggunaan element shell continuum dapat mencapai konvergensi dan bahkan dengan cukup cepat diterapkannya mesh refinement, dan ini merupakan indikasi yg menunjukan penggunaaan element tersebut lebih akurat. dibawah adalah tampilan kontur distribusi tegangan von-Mises



 

to be added,..

yg belum beam 1D element conventional/standard dengan beam continuum 3D, Kirchof thin shell & Timoshenko thick beam serta element solid true. nanti ada waktu dibuat menyusul kemudian perbandingannya,

lainnya analisa nonlinear plastisitas mencari nilai tegangan residu dan deformasi permanen akibat pemutaran penampang balok sebesar 90deg pada jung yg sebelumnya tegak seperti posisi pada tumpuan.

 

*updates

hasil analisa dengan menggunakan conventional 1D beam element, Euler mengabaikan deformasi geser.

(8 element, 9 nodes)

penghalusan jumlah pias pembagi (mesh refinement),

(16 element, 17 nodes)

(32 element, 33 nodes)

kontur tegangan arah sumbu utama pada titik sudut tepi penampang balok,

hasil analisa menggunakan elemen beam continuum 3d quadratic (B32),

(16 element, 200 nodes)

tampilan kontur distribusi tegangan von-Mises,



penghalusan mesh transisi kelengkungan,





 

**Note: Finite Element Solver are LISA-FET/MW and CalculiX

Read More