gambar diatas adalah portal sederhana dua dimensi dengan beban lateral & gravitasi, material beton dimensi balok kiri (30x50cm) dan balok kanan (30x80cm) kolom kiri & tengah (40x40cm) dan kolom kanan (40x100cm).
ditinjau model pertama (1) semua tanpa modifikasi
Nodal Displacements:
Dx = -7.114e-001 mm
Dy = -8.452e-002 mm
Rz = 6.627e-004 rad
Bending Moment Diagram (kN*m)
Shear Froce Diagram (kN)
Axial Force Diagram (kN)
ditinjau model kedua (2) modifikasi kekakuan dengan meningkatkan modulus elastisitas daerah end rigid joint menjadi 100x
Nodal Displacements:
Dx = -8.457e-001 mm
Dy = -7.713e-002 mm
Rz = 6.101e-004 rad
Bending Moment Diagram (kN*m)
Shear Froce Diagram (kN)
Axial Force Diagram (kN)
ditinjau model ketiga (3) semua tanpa modifikasi kecuali balok kanan, kedua ujungnya dilakukan modifikasi kekakuan dengan menurunkan modulus elastisitas daerah end release joint menjadi 0.01x atau 1/100
Nodal Displacements:
Dx = 6.877e-001 mm
Dy = -7.939e-002 mm
Rz = -1.680e-004 rad
Bending Moment Diagram (kN*m)
Shear Froce Diagram (kN)
Axial Force Diagram (kN)
terlihat agak rumit karena kondisi rigid end joint juga perlu dimodelkan pada kolom lebar/shearwall. dicoba contoh yg lebih sederhana modifikasi kekakuan untuk mewakilkan kondisi rotational release.
Max. Transv. Displ.: 3.255e+000 mm
Bending Moment Diagram (kN*m)
Shear Froce Diagram (kN)
modifikasi kekakuan ujung dengan merubah nilai modulus elastisitas sebesar 0.01x atau 1/100, panjang segmen balok ujung 0.01Lt atau 1/100
Max. Transv. Displ.: 1.149e+001 mm
Bending Moment Diagram (kN*m)
model lain dengan nilai modulus elastisitas 0.001x atau 1/1000 kondisi normal
Max. Transv. Displ.: 1.510e+001 mm
terlihat daerah pertemuan masih menerima momen sebesar ~5% namun ini sudah dapat dikategorikan jenis sambungan sendi atau simple pinned.
modifikasi kekakuan seperti diatas, hanya dapat digunakan dengan perhatian (using with caution) karena masalah numerik solver dalam komputasi. namun program FE saat ini sudah banyak yg mengggunakan advanced solver, lalu apakah hal tersebut masih menjadikan masalah (?)
**updates
analisa beam continum expansi 3D dengan elemen hingga, kondisi normal tanpa modifikasi end release
pada metode elemen hingga program bantu CalculiX elemen balok dilakukan expansi menjadi 3d continuum element, sehingga modifikasi properties penampang seperti inersia tidak dapat dilakukan. karena balok sudah merupakan element solid 3d maka yang dilakukan adalah properties materialnya yaitu modulus elastisitas, berbeda keadaannya pada analisa balok dengan program analisa struktur standard yg menggunakan one dimensional elemen perubahan material yg isotropis tidak terlihat berdampak pada stabilitas solver. pastinya tidak dengan FEM yg menggunakan element solid, perlu menggunakan material anisotropis (sumbu 1,2 &3) pada daerah end release untuk mensimulasikan tumpuan dapat menahan geser namun tidak dapat menahan rotasi.
analisa beam continum expansi 3D dengan elemen hingga, kondisi normal dengan modisikasi material daerah end release, dibuat 1/1000 konndisi material normal.
terlihat hasilnya sudah sesuai pada defleksi dan tegangan, dimana tegangan one dimensional adalah sebesar 15.62N/mm^2 (f=M/S) hanya ada sedikit yg terlihat berbeda pada koordinat tepat di tumpuan, mungkin ini dikarenakan penentuan multi point constraint dan penggunaan element quadratic.
***updates
ada yg terlewat pada penentuan hubungannya modulus geser, sehingga menjadikan node titik tumpuan tidak stay plane terhubung rigid links seperti tampilan diatas. berikut setelah diperbaiki, hasilnya sesuai yg diperkirakan.
deformasi sejajar serat balok sumbu axis/normal antara analisa dengan tumpuan jepit dan simulasi end release menunjukan nilainya sekitar 10 kali lipatnya pada serat bawah dan atas bagian dekat tumpuan.
perkiraan saya akan ada masalah dipemodelan portal diatas karena adanya gaya aksial pada balok. bagaimana melakukan modifikasi porperties material untuk simulasi hubungan dapat menerima tekan/tarik namun tidak dapat menerima momen (?)
**updates
model (1) kondisi material normal tanpa modifikasi
hasil beam 1D (Dx = -0.7114 mm)
hasil beam 1D (Dy = -0.0845 mm)
hasil beam 1D (f=6.86 N/mm^2)
hasil beam 1D (f=0.80 N/mm^2)
... to be added,
beam 3D continuum for model 2 & 3 of frames simulation in example above.
0 komentar:
Posting Komentar