Minggu, 16 Agustus 2015

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.

Kamis, 06 Agustus 2015

CAE - mesher solid brick linear hexahedral

postingan sebelumnya ingin mencoba konversi file *.inp ABAQUS hasil program mesher GMSH menjadikan format yang dapat dibaca SAP. program mesher tersebut sudah lamya sya tidak menggunakannya jadi agak lupa (sekalian posting disini buat pengingat sendiri juga), ini saya coba lagi pada object sederhana balok memanjang dan balok lengkung.

2015-08-07 07_53_22-Gmsh - C__Users_synt_Desktop_meshext.geo

2015-08-07 18_13_55-Statistics


2015-08-07 07_52_26-Gmsh - C__Users_synt_Desktop_meshext.geo 2015-08-07 07_53_44-Source 2015-08-07 08_06_39-Gmsh - C__Users_synt_Desktop_unsquadext.geo

2015-08-07 18_13_09-Statistics


2015-08-07 08_07_16-Gmsh - C__Users_synt_Desktop_unsquadext.geo 2015-08-07 08_07_42-Source

coba meshing yang agak rumit ketegori full hexa,

2015-08-07 20_02_14-Gmsh - C__Users_synt_Desktop_platehole.geo

2015-08-07 20_04_16-Statistics

2015-08-07 20_01_38-Gmsh - C__Users_synt_Desktop_platehole.geo

2015-08-07 21_09_21-Source

terilahat masih adanya element prismatiks/wedge (50 elmt.), ini karena boundary edges terdapat lima sisi. padahal jika meshing bidang 2D algorithma GMSH dapat berhasil full quad walau dengan jumlah boundari edges lebih banyak dari 4 sisi, saya coba lagi dgn partisi ulang:

2015-08-07 20_22_00-Gmsh - C__Users_synt_Desktop_platehole2.geo

2015-08-07 20_23_08-Statistics

2015-08-07 20_22_52-Gmsh - C__Users_synt_Desktop_platehole2.geo

2015-08-07 21_11_01-Source

penghalusan mesh dan penambahan tebal object,

2015-08-08 03_20_44-Gmsh - untitled.geo

2015-08-08 03_20_34-Statistics

terlihat hasil meshing masih full hexa. ya, cara partisi menjadikan boundary edges dibuat empat sisi dapat berhasil membuat meshing yg full hexahedral. perlu diperhatikan setting jumlah pembagi mesh pada edge boundary, jika tidak konsisten antar sisinya tetap akan menghasilkan element prism/wedge.

mengenai menjembatani import di SAP nanti akan saya coba pelajari ulang dulu untuk element solid tersebut, dikarenakan jenis 3D solid walau linear saja pengurutan posisi node berbeda tiap program apalagi yg quadratic. hal ini berbeda dengan element tria/quad yang biasanya untuk linear sama berlawanan jarum jam, kecuali yang quadratic.

tambahan catatan:

formulasi element solid brick untuk yang linear saja ada beberapa macam (full integration, reduced integration, incompatible mode dan hybrid) membaca theory manualnya ini hubungannya dengan masalah shear-locking dan lainnya masalah hour-glassing.

2015-08-08 02_16_05-graphics.cs.uh.edu_wp-content_papers_2015_2015-SIGGRAPH-HexmeshAlignment.pdf - C

(source: Deng Z., et al)


membaca dan mengikuti banyak paper, notes dan diskusi dari forum, pendapat para analyst yang sudah banyak berpengalaman (kebanyakan komentar dari tekhnik mesin mekanika) performance element ditentukan dari geometri dan kasus yang ditinjau (contact/plasticity) banyak yang merekomendasikan untuk menggunakan elemen quadratic brick/hexa (as possible) atau quadratic tetra (auto meshing), namun jika yang dihadapi adalah masalah contact, element yg direkomendasikan adalah linear brick/hexa incompatible (fine mesh) karena masalah konvergensi pada element quadratic dengan adanya reaksi pegas kontak yang terbalik di extra nodes, belum lagi tinjauan elemnt quadratic terhadap masalah dinamis?

berikut cuplikan dari pembuat SAP,

"Because of the serious problem associated with shear-locking, the classical compatible four-node quadrilateral and eight-node hexahedral elements should not be used to simulate the behavior of real structures. It has been demonstrated "that the addition of incompatible displacement modes, corrected to pass the patch test, significantly enhances the performance of quadrilateral and hexahedral isoparametric elements.

2015-08-07 15_19_59-edwilson.org_BOOK-Wilson_05-iso.pdf - Citrio

The nine-node quadrilateral and the 27-node hexahedral elements are accurate and can be improved by adding corrected incompatible modes. For example, cubic modes can be added to the nine-node plane element in which the exact results can be calculated, for tip shear loading, using only one element to model a cantilever beam."

element order and conectivity di SAP dan rekomendasi dari program theory manual,

2015-08-07 15_32_03-CSiRefer.pdf - Adobe Reader

"Solid element is an eight-node element for modeling three- dimensional structures and sol ids. It is based upon an isoparametric formulation that includes nine optional incompatble bending modes.

The locations of the joints should be chosen to meet the following geometric conditions:

  • The in side angle at each corner of the faces must be less than 180°. Best resultswill be ob tained when these an gles are near 90°, or at least in the range of 45° to 135°.

  • The aspect ratio of an element should not be too large. This is the ratio of the longest dimension of the element to its shortest dimension. Best results are obtained for aspect ratios near unity, or at least less than four. The aspect ratio
    should not exceed ten.


These conditions can usually be met with adequate mesh refinement."

dan berikut element order linear brick di ABAQUS,

2015-08-07 15_55_43-e3dsolid-node-face-num.png (521×575) - Citrio

"C3D8I 8-node linear brick, incompatible modes"


jadi di SAP menggunakan formula element SOLID yang sebanding lenear brick/hexa C38I nya ABAQUS. mengenai konversi input files kelihatan sudah, ini ada masalah kalo input ABAQUS mau dimodifikasi menjadikan terbaca di SAP, node element order posisinya berbeda. di ABAQUS masih menggunakan prinsip urutan berlawanan jarum jam sedangkan di SAP tidak seperti itu.