Kamis, 21 Juli 2016

tinjauan singkat FE dgn material kemampuan terbatas

masalah nonlinear yg paling sederhana adalah adanya batasan elemen (truss) dari kemampuan yang hanya menerima tarik atau tekan, dapat juga ditentukan elemen hanya tarik namun mempunyai sedikit batasan tekan atau kondisi sebaliknya. pada penggunaan software sya sebelumnya SAP yaitu dikenal dengan frame with tension only atau compression only. elemen frame jenis tension only biasa digunakan pada komponen gedung untuk bracing, namun pengguna dapat memberikan suatu nilai kecil batasan jika diinginkan mempunyai kemampuan tekan. saya tertarik mencari rujukannya, misal untuk jenis brace tsb atau frame dgn tension only maka analisa adalah nonlinear dilakukan secara bertahap atau iterasi, jika pada tahap pertama terdeteksi adanya suatu brace yg mengalami nilai tekan diatas yg ditentukan maka akan dianggap tidak ada pada tahap analisa berikutnya dst.

2016-07-21 21_38_04-Clausen2006

criteria for Modified MC materials (source: Clausen etal, 2006)


berbeda dgn program FE nonlinear, umumnya yg diutamakan adalah perilaku nonlinear materialnya. untuk material yg dominasi kemampuannya tekan dan getas maka akan banyak pilihan pada program advanced FE seperti adanya jenis nonlinear material Mohr Coulomb atau Drucker Prager yg biasa diaplikasikan untuk pemodelan pondasi atau dinding penahan tanah dan pemodelan beton bertulang. jenis material elasto-plastis yg lebih sederhana dari kedua material elasto-plastis diatas adalah Compression_Only material yg sudah tersedia pada program elemen hingga CalculiX versi 2.8 keatas.

2016-07-23 23_53_59-Kupfer1973

perbandingan kurva diatas dengan material beton hasil pengujian tegangan biaxial dari Kupfer (1973) yg sya ambil dari paper Como etal (2007) ada sedikit perbedaan dimana pada tegangan tekan biaxial maka kemampuannya akan meningkat sampai sekitar 25% dari kondisi tegangan uniaxial.

2016-07-21 21_59_24-Source

disini saya mencoba meninjau secara singkat dan belum sampai ke perbandingan pendekatan secara teoritis, contoh adalah balok beton bertulang (bxh=20x40cm) kantilever (l=2.0m) dengan beban merata pada permukaan atas. mutu tekan beton (fc=30MPa) kapasitas tegangan tarik material ditentukan sebesar 10% dari kapasitas tegangan tekan umumnya. mutu tekan beton tsb dalam hal ini hanya sebagai acuan tidak berpengaruh pada analisa dgn elemen hingga menggunakan material Compression_Only.

 

2016-07-21 21_54_02-Calculix Graphix

tulangan dimodelkan dengan elemen truss (T3D2) mengunakan material yg masih elastis, hubungan dengan node elemen solid beton dianggap monolit sempurna.

2016-07-21 21_50_45-Calculix Graphix

deformasi,

2016-07-21 21_11_06-Calculix Graphix

2016-07-21 21_11_41-Calculix Graphix

tegangan lentur,

2016-07-21 21_10_20-Calculix Graphix

2016-07-21 21_09_23-Calculix Graphix

tegangan geser,

2016-07-21 21_13_39-Calculix Graphix

walau belum dilakukan validasi atau verifikasi namun sekilas hasil dan perilaku menunjukkan seperti yg diperkirakan. tegangan tarik sebagian besar diterima oleh tulangan  atas, sedangkan tegangan geser hanya diterima pada balok yg mengalami tekan (tidak retak). masih ada beberapa yg belum sesuai, diantaranya adalah sisi tumpuan bagia atas yg mengalami tarik seharusnya sudah retak dan kemudian tegangan tarik diterima penuh oleh tulangan.

dibawah dicoba dengan mengabaikannya kapasitas tarik material pada elemen solid balok beton, ditentukan nilai hanya sebesar 1% kuat tekannya. dan hasilnya menunjukan sebagai berikut,

2016-07-21 22_40_38-Calculix Graphix

terlihat saat mengabaikan kemampuan tarik material beton menjadikan tegangan tarik baja tulangan atas dan tegangan tekan balok beton sisi bawah meningkat cukup signifikan.

*updates

model balok diatas menggunakan material beton elasto-plastis kriteria Modified Mohr Coulomb (Abaqus UMAT by J. Clausen, 2006)

defleksi,

2016-07-23 02_44_56-Calculix Graphix

tegangan lentur,

2016-07-23 02_46_24-Calculix Graphix

2016-07-23 02_50_14-Calculix Graphix

tegangan geser,

2016-07-23 02_46_46-Calculix Graphix

2016-07-23 02_49_01-Calculix Graphix

perlu model yg lebih realistis yaitu representasi tulangan dengan elemen solid, ditambahkannya pengaruh tahanan geser dan kekangan dari tulangan sengkang. juga pelu diterapkan nonlinearitas materil baja tulangan tersebut. lainnya yg tidak boleh terlewat adalah  perbandingan dengan analisa secara pendekatan teoritis perhitungan tangan sebagai acuan awal.

 

**updates

model diatas agak rumit terutama saat pemodelan tulangan sengkang, sya coba model yg lebih sederhana untuk mencari perkiraan awal perbedaan pemodelan tulangan pokok dengen elemen truss (t3D2) dan full solid. ditinjau balok 50x150mm panjang 500mm.  material hanya tekan (kapasitas 10%), beban merata.

lendutan,

2016-07-23 01_12_44-Calculix Graphix

tegangan lentur,

2016-07-23 01_11_26-Calculix Graphix

2016-07-23 01_16_37-Calculix Graphix

tegangan geser,

2016-07-23 01_13_13-Calculix Graphix

2016-07-23 01_17_39-Calculix Graphix

terlihat hasilnya mendekati, dapat diartikan model tulangan pokok dgn elemen truss  (T3D2) dapat digunakan untuk memprediksikan awal. perbaikan model tetap diperlukan seperti model tulangan dgn elemen solid, diikut sertakannya model tulangan sengkang dan  kriteria elasto-plastis materialnya.

***updates

model balok beton bertulang sederhana diatas menggunakan material elasto-plastis kriteria Modified Mohr Coulomb (Abaqus UMAT by J. Clausen, 2006)

2016-07-23 02_22_16-Source

defleksi,

2016-07-23 02_10_46-Calculix Graphix

tegangan lentur,

2016-07-23 02_11_35-Calculix Graphix

2016-07-23 02_17_10-Calculix Graphix

tegangan geser,

2016-07-23 02_12_01-Calculix Graphix

2016-07-23 02_18_11-Calculix Graphix

tegangan tarik pada balok beton yg terjadi pada sisi atas masih cukup besar, ini kemungkinan besar disebabkan dari penentuan nilai parameter dari material modified MC yaitu angle of internal friction, dilatancy dan cohesion. perlu dicari rujukan dalam penentuan nilai tersebut hubungannya dengan mutu tekan beton (?)

 

2016-07-23 03_09_19-MCforConcrete

source: Ardiaca, 2009


 

menggunanakan parameter nilai kohesi (c=3.0 N/mm^2) sedangkan nilai lainnya tetap.

lendutan,

2016-07-23 21_50_22-Calculix Graphix

tegangan lentur,

2016-07-23 21_49_35-Calculix Graphix

penekanan tampilan warna batas/kontur sisi bawah - tegangan tekan,

2016-07-23 21_59_16-Calculix Graphix

penekanan tampilan warna batas/kontur sisi atas - tegangan tarik,

2016-07-23 22_34_25-Calculix Graphix

tegangan geser,

2016-07-23 21_49_52-Calculix Graphix

2016-07-23 21_59_52-Calculix Graphix

terlihat nilai kohesi erat hubungannya dengan batasan tarik Rankine atau "tension cut-off" yg mana nilai tersebut adalah batasan maksimum tegangan tarik yg diperbolehkan terjadi. berdasarkan rujukan Eurocode nilai cohesi beton cukup kecil yaitu hanya 0.55 Nmm^2 untuk beton mutu f'c 30MPa dan pada Abaqus CAE menjadikan masalah stabilitas numerik seperti yg disampaikan oleh Chiruta G. (2014) dgn supervisor prof. J. Clausen (2015) bahkan pada analisa sejenis hanya digunakan nilai kohesi sebesar 7.5N/mm^2 dan ini jauh lebih besar dari yg sya gunakan hanya 3.0N/mm^2 pada CalculiX.

sya sudah mencobanya di CalculiX dengan parameter nilai kohesi tersebut serta penggunaan sudut geser friksi sebesar 9deg rekomendasi dari Eurocode, konvergensi solver cukup lama dan belum sampai selesai.. nanti saya coba lagi.

***updates

menggunakan nilai kohesi (c=1.5Nmm^2) mendekati nilai yg direkomendasikan oleh Montoya (1971) yg diambil dari paper Ardiaca (2009), nilai sudut geser friksi dan dilatansi dibuat tetap.

2016-07-24 16_29_52-Iss25 Art2 - Montoya1971

running cukup lama dan solver tidak mencapai konvergensi, hanya sekitar 81% tahap beban diterapkan. artinya hasil tidak dapat digunakan, namun sebagai gambaran saja sebagai perbandingan dgn sebelumnya.

2016-07-24 16_11_31-rcbeam2.inp - SciTE

defleksi pada step ~81%

2016-07-24 16_17_27-Calculix Graphix

tegangan lentur,

2016-07-24 16_19_59-Calculix Graphix

 

 

tegangan geser,

2016-07-24 16_20_31-Calculix Graphix

 

kegagalan solver tersebut yg tidak mencapai konvergensi perlu dicari penyebabnya (?) dan ini ada kemungkinan dari penerapan beban sya yg terlalu tinggi, diperkirakan sudah gagal karena berkisar 3x lipat kapasitas geser balok.

sebagai informasi tambahan mengenai rekomendasi nilai kohesi untuk beton dari :

  • CEB-FIB Model Code (2010) menerangkan nilai tersebut diantara 1.5 - 3.5 N/mm^2

  • Mattock (1974) menerangkan nilai tersebut berkisar 2.8 N/mm^2

  • AASHTO LRFD (2007) dapat diambil dari hasil uji geser langsung, tau_v = c * A_cv

  • Vecchio etal (2001) dapat diambil dgn rumusan hubungan dgn sudut geser friksi,


2016-07-24 18_18_23-c

jadi penggunaan nilai yg digunakan terakhir sebesar 3.0 N/mm^2 pada data masukan analisa cukup reliable untuk kepemilikan beton pada umumnya.

sedangkan untuk nilai sudut geser friksi atau sudut friksi dalam (internal angle friction) hasil penelitian Suenaga (1998) untuk beton normal (Fc~26.5Nmm^2) menunjukan nilainya berkisar antara 37 - 55 deg dan nilai maksimum tsb mendekati dgn nilai sesuai Wegde Theory.

2016-07-23 22_50_13-5.1.5 Modified Mohr-Coulomb - Opera

source: DIANA FEA (2010)


biasanya analisa model dgn material beton menggunakan kriteria Drucker-Prager namun pada CalculiX belum diimplementasikan UMAT tersebut. terlihat dari rujukan grafik diatas, penggunaan nilai sudut geser friksi yang kecil akan menjadikan bentuk kurva modifikasi Mohr-Coulomb mendekati kurva Drucker-Prager.

 

 

0 komentar:

Posting Komentar