titik berat penampang bidang segi-empat sembarang

[ draft ]

metode sebelumnya yg sya gunakan adalah dengan bantuan luasan batas luar (bounding box) ternyata saya coba cukup panjang dan melibatkan banyak pencabangan, unfinished yet. sya coba cara yg lain yaitu dengan membagi bidang segi-empat sembarang tersebut menjadi dua buah bangun segi-tiga bidang. sebenarnya ini adalah cara grafis yg disimulasikan dengan metode analitis, metode grafis tersebut biasanya menggunakan penggaris manual untuk mecari titik berat masing-masing segi-tiga tersebut.

pencarian titik garis perpotongan dari pusat segi-tiga tidak ditampilkan lagi disini karena ini juga sudah pernah dibahas sebelumnya. kelebihan dari metode ini adalah dapat digunakan pada mesh yg triangular hanya dengan sedikit saja modifikasi algorithma.

example case (1)

2015-11-15 23_41_20-center_arbitary_quad.ods - OpenOffice Calc

2015-11-15 23_49_26-center_arbitary_quad.ods - OpenOffice Calc

segment pembentuk segi-empat sembarang (1)

2015-11-15 23_49_44-center_arbitary_quad.ods - OpenOffice Calc

segment pembentuk segi-empat sembarang (2)

perbandingan dgn metode perhitungan langsung dari Zienkiewicz (1977)

2015-11-16 04_58_20-SMath Studio Desktop - [segitiga_luascenter.sm]

2015-11-16 04_58_41-SMath Studio Desktop - [segitiga_luascenter.sm]

terlihat metode langsung tsb lebih cepat dibandingkan cara grafis dgn mencari titik perpotongannya. ada yg beda pada perhitungan luas mengapa selalu bernilai negatif (?) apa cukup dengan merubah menjadi nilai absolutnya, selain itu juga pada segment (2) hasil nilainya selisih -3.2%. perlu banyak test case sebagai pembanding.

example case (1)

2015-11-15 23_45_09-center_arbitary_quad.ods - OpenOffice Calc

example case (2)

2015-11-15 23_46_08-center_arbitary_quad.ods - OpenOffice Calc

example case (3)

2015-11-15 23_46_56-center_arbitary_quad.ods - OpenOffice Calc

example case (4)

2015-11-15 23_47_26-center_arbitary_quad.ods - OpenOffice Calc

dari langkah diatas didapat luasan dan posisi titik berat beratnya sehingga selanjutnya akan didapat properties dan stress dari penampang sembarang (arbitrary section),

  • sectional areas, normal stress
  • sectional inertias (strong/weak) axes – top/bottom & left/right, bending stress
  • combined stress (normal +/- bending) at any point

to be added,

  • more example case and algorithm validation/compared with available FEM programs
  • reading coordinate from meshing programs and make further calculation as describe previously.
  • shear center, effective shear areas (?)
  • torsional properties & stress (St. Venant)
  • warping properties & stress (VZ. Vlasov)
  • plot stress contour with color ranges (?)
  • study meshing algorithm of Delaunay Triangulation, known less complexity copared with quad algorithm – implemented in personal mesher programs(?)

2015-11-16 01_46_22-Gmsh - C__Users_synt_Downloads_gmsh-2.10.1-Windows32_gmsh-2.10.1-Windows_ysect.g

example of coarse mesh unsructured quad (GMSH)

2015-11-16 01_48_28-Gmsh - C__Users_synt_Downloads_gmsh-2.10.1-Windows32_gmsh-2.10.1-Windows_ysect.g

zoom detail in flanges and display node numbers

2015-11-16 01_45_48-Gmsh - C__Users_synt_Downloads_gmsh-2.10.1-Windows32_gmsh-2.10.1-Windows_ysect.g

triangular coarse mesh result for same boundary edges (GMSH)

shewchuk2012

triangular meshing steps (source:Shewchuck, 2012)

One thought on “titik berat penampang bidang segi-empat sembarang

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