balok baja, sama berat namun beda kekuatan lebih

.



.



.

LTB

.



.



.



.



.

Balok baja dua tumpuan jepit, beban merata, Mtump = 2*Mlap, optimasi dgn haunched, sliced section. kriteria kekuatan: tegangan von mises maksimum, stabilitas dan serviceability:lendutan. Questionable, LTB for haunced? Cb values? elastic or plastic, residual stress effect? etc. Studying advanced problem using Code_Aster as FE solver.

.

Validation - Calculation pad.

.

Lb = 1200
w = 1.00
q = w*20 = 20.00
M1=(1/12)*q*Lb^2 = 2400000.00
Z = 2480.62
Sigma1=M1/Z = 967.50
' Bot = 1351.9 kg/cm2
' Top = 1408.7 kg/cm2
M2=(1/24)*q*Lb^2 =1200000.00
Sigma2=M2/Z = 483.75
' Bot = 708.9 kg/cm2
' Top = 672.3 kg/cm2
Es=2000000
nu = 0.30
G = Es/(2*(1+nu)) = 769230.77
h = 60.00
b = 20.00
tf = 1.70
tw = 1.10
Iy = ((1/6)*tf*b^3)+((1/12)*(h-tf)*tw^3)  = 2273.13
d = h - 2*tf  = 56.60
Ix = ((1/12)*b*(d+2*tf)^3) - ((1/12)*(b-tw)*d^3) = 74418.64
J = (1/3)*((2*b*tf^3)+(h*tw^3)) = 92.13
Cw = (h^2*Iy)/4 = 2045819.80
w = 1.00
q = w*b = 20.00
As = h*tw = 66.00
delta = ((q * Lb^4) / (384*Es*Ix)) + ((q*Lb^2)/(8*G*As)) = 0.80
K = 1.00
Ms = Min(M1,M2)*(-1) = -1200000.00
Ml = Max(M1,M2) = 2400000.00
Cb = 1.75+1.05*(Ms/Ml)+0.3*(Ms/Ml)^2 = 1.30  < 2.3
Mcra = Cb*(Pi/(K*Lb))*Sqrt((Es*Iy*G*J)+((Pi*Es)/Lb)^2*Iy*Cw)  = 2282234.25
Sigmacra = Mcra / Z  = 920.03
Lamda = 1.21034
wcr = w * Lamda  = 1.21
qcr = wcr * b=24.21
Mcrfe = 1/24 * qcr * Lb^2=1452408.00
Rltb = Mcra/Mcrfe=1.57
'
Mmax =  2400000.00
rx = 0.25
x = rx * Lb= 300.00
Ma = (q/12)*(6*Lb*x - Lb^2 - 6*x^2) = 300000.00
rx = 0.50
x = rx * Lb= 600.00
Mb = (q/12)*(6*Lb*x - Lb^2 - 6*x^2) = 1200000.00
rx = 0.75
x = rx * Lb = 900.00
Mc = (q/12)*(6*Lb*x - Lb^2 - 6*x^2) = 300000.00
'
Cb = (12.5*Mmax) / (2.5*Mmax + 3*Ma + 4*Mb + 3*Mc) = 2.38
' Conservatively takes Cb = 1.00

.

all units are in (Kgf-cm)

.

Keterangan tambahan, dilain kesempatan :)

Read More

little bit of try GUI Salome

tahun lalu sya sdh coba Salome-Code_Aster, namun karena dokumentasi perintah analisa FE nya Code_Aster dlm bahasa france jadi mandeg, hanya sebatas prinsip analisa linear aja yg baru ngerti, lainnya seperti contact dan buckling dan juga nonlinearitas geometry&material belum paham perintahnya. Sedangkan GUI modeling dgn Salome sudah bagus dokumentasinya 'coz dlm bahasa english yg bisa sya phamin, dan lagi prinsip penggunaan seperti dulu sya pke AutoCAD 2000 buat bikin model 3D, ada extrude,revolve, path, boolean combine, substarct bedanya di CAD ngga terorganisir komponen pembentuknya. Dapat membuat solid model dgn kompleksitas tinggi karena komponen perangkai model sudah terorganisir. Kalo yg sya buat sperti dibwah ini sih masih sederhana, namun bisa dibilang kompleks kalo sya modeling di SAP2000 lsg, yg memungkinkan model sperti ini  import solid model DXF dari CAD.



Salome-Meca skrg sudah ada yg versi 2009, namun saat dinstall di pc sya Ubuntu 8 ngga bisa di-run Code_Aster nya, kemungkinan hubungannya dgn administrasi pd saat instalasi.



Ada khabar bagus juga :) dulu CAELinux menggunakan base pcLinuxOS (KDE) skrg yg versi 2009 beralih ke Ubuntu (GNOME) yg lebih stabil, namun khbar :( kurg baiknya (bgi sya) hanya bisa run di 64bit processor 'coz pc sya masih 32bit.



Diatas adalah tampilan dri meshing quality control, sehingga jika pengguna kurang puas dapat diubah lagi algorithma dan setting lainnya.

Khabar baik lainnya sudah mulai lebih banyak dokumentasi Code_Aster dlm bahasa english, namun itu tadi masih sebagian kecil sja yg sdah translate masih jauh belum mencakup bahsan kemampuan dri Code_Aster secara keseluruhan. Kalo dluar sih sudah banyak peminat serius dri praktisi industri, terlihat dri pelatihan yg ada sprti berikut.



kalo sya blum serius banget, lebih cenderung minat dan mungkin bisa dikatakan hobby banding2in textbook dgn program. Bagi sya yg penting modeling sma analysis ngga bener2 menyita waktu aja, hobby=ga serius=ga berarti?? ya ga gitu lah sempit bgt .. "i think i did better understand with FE"

Updates :

Code_Aster berhasil di-run dan menampilkan hasilnya, yg sya lakukan adalah me-remove yg versi 2009 diganti dgn versi 2008 krena kemungkinan besar dgn Ubuntu 8 OS sya terpasang, sebenarnya bisa aja upgrade ke Ubuntu 9 biar Salome-Meca versi 2009 bisa jalan namun belum sempet nanti lagian sya masih beginner belum yg advanced.



Model yg pertama kehapus juga, jadinya bikin lagi. tapi gpapa hanya butuh bbrapa menit aja buat modeling,meshing dan run-analysis nya.



Diatas adalah meshing tetrahedral, dimodelkan juga fillet radius pada balok WF comb.



Tampilan kontur lendutan arah -Z, untuk mengetahui nilai pada suatu titik koordinat dilakukan dgn cara pilih node kemudian display values. Boundary condition pada daerah ujung kanan adalah jepit all DOF restrained, sedangkan pada top flens diterapkan beban merata, ini belum sampe ke verifikasi :) namanya jga little bit of try GUI capabilities .



Kontur tegangan kriteria leleh von misses, ditampilkan dengan iso-surfaces.



Kontur tegangan arah sumbu-Y kontur ditampilkan dgn warna untuk membedakan tegangan tarik dan tekan, ini seperti pada umumnya post-processor program FE. Sedangkan gambar dibawah adalah tegangan ZY, sejajar bidang web karena pda penampang WF geser akan sebagian besar diterima oleh web.



Jenis element masih linear, kurang baik kalo linear bagusan quadratic atau element hexahedral. Sebenarnya sudah disediakan fasilitas pd Salome, pengguna tinggal pilih mesh dan convert aja ke quadratic. Namun hal ini meminta resource memory cukup besar, kalo gitu pc sya perlu tambah lagi 1GB RAM ini keliatannya, ato perlu sabar bntar saat nungguin run...

Tambahan, model struktur dibawah ini termasuk kategori rumit butuh software modeling+meshing yg advanced sekelas ANSYS atau ABAQUS. Hasil yang ditampilkan berikut adalah: mesh, deformed shape (deflection), stress distribution (von misses, axes).



.



.



fasilitas yg terlihat bermanfaat dgn cara visual adalah pilihan tampilan clipping planes, dengan cara tersebut pengguna dapat melihat distribusi tegangan suatu lapis ketebalan elemen pada potongan yg ditentukan sembarang.

Mungkin nanti sya akan bikin tutorial modeling, meshing & analysis serta perbandingannya dgn software FE yg sya udah familiar S2K, perlu pake software advanced meshing yg di-import ke S2K untuk model kompleks dgn shell atau solid element. Di awal jenis analysis linear elastis dulu mungkin, karena di Code_Aster baru ini yg saya paham, pengennya sih lebih ke jenis analisa masalah kontak, nonlinearitas material/geometry tapi tentunya akan sebanding dgn my learning curves using C_A, yah namanya juga keinginan.

from the spirit of free software, now its times to my pingus goes pro :)

Read More

mesin penyelesai persamaan

equation solver dlam suatu program FE merupakan hal yg paling menentukan hasil ketepatan dan kecepatan menyelesaikan persamaan matriks dgn jumlah yg sgt besar, pada prinsipnya semua program FE menylesaikan masalah persamaan yg sama yaitu {F} = {K} * {u} pda masalah mekanika statis linear, pada masalah dinamis selain matriks kekakuan akan adanya tambahan matriks massa {M} dan redaman {C}. Sedangkan pada masalah tekuk dan ragam getar akan mencari nilai eigen.



SAP2000 pada sejak versi 11 (setahu saya) sudah mempunyai solver advanced,  tahun lalu pertama kali pake kaget juga saat nge-run model dgn DOF's mencapai ~3500 dgn jenis analisa statis, dynamics & p-delta, diselesaikan hanya dgn waktu kurang dri 15 detik saja wahhh.



Ada rasa ingin tau aja beda bgt dgn yg saya pake jaman kuliah dulu yaitu versi 7,  dan jawabannya adalah seperti diatas yg terlampir juga dalam dokumentasi.



Ternyata itu bukan sesuatu yg baru dan bukan murni penelitian dari tiap2 developer, namun merupakan hasil dari penelitian suatu institusi/perorangan. Implementasinya aja baru belakangan seperti yg banyak dilakukan oleh pengembang software FE komersil. Setiap developer dri akademis boleh menggunakan math library dri hasil project tsb, kalo untuk developer komersil maka berbayar namun cheap (seharga motor 2nd jadul disni). Cheap sya pikir bagi pengembang besar sekelas Autodesk/Bentley atau lainnya untuk sebuah engine solver hebat seperti PARDISO, ... ini win-win solution ngga sih :(

Software FE yg mengunakan solver advanced tsb, diantaranya: MSC Nastran, CalculiX, Robot Structural Analysis, RAM Concept, NISA, TNO DIANA, COMSOL,  dll belum tau.



.



.



.



.

.



.



banyak program FE besar dan komersil mengunakan mesin penyelesai persamaan yg sama, kalo gitu hasil akan sama jika diberikan perintah yg sama. Jadi yg membedakan command & management data input/output saja? atau lain apa?

Note:

program dri negeri sendiri SANSPRO juga sedang dikembangkan menggunakan advanced solver, cuman dokumentasi ngga jelasin pake apa ato yang mana.

Read More

tambahan modul pd program glassboxDI #3

beberapa tambahan modul desain beton bertulang pada program glassboxDI v0.16 diantaranya adalah diagram interaksi (continued) dengan jumlah lapis tulangan saat ini sampai 6lapis.

.

.

.

.

.

.

Sebenarnya ingin bikin auto juga gambar2nya dimensi, plot Pn,Mn cuman waktu buat dedikasinya buat ngerjain karena ini bisanya di GUI bukan CLI yah nati ini lagi ngoprek bahasanya cukup ribet banyak aturan tpi bakalan lebih usefull nantinya :)



selain itu, ditambahkan pula modul desain pondasi tapak dgn posisi kolom center. keadaan posisi kolom tidak pada center pada akn dibuat lanjut.

[scribd id=25335503 key=key-zl4gtlo56zu0pbw7iha]

ada sedikit berbeda karena melibatkan proses iterasi mencari nilai lebar pondasi yg tergantung dari nilai P,Mx,My dan tegangan tanah ijin. mungkin nantinya dibuat lebih transparan dgn ditampilkan beberapa step saat mendekati nilai konvergensi seperti berikut,



.

.

snap shoot gambarnya  kekecilan, dibawah di-zoom biar lebih jelas maksud sya step diliatin saat proses.



.

work in progress ... juga untuk desain pilecap, next retaining wall mungkin untuk yg sub-structures.



.

Read More

freetime, mainan 3d lagi

bulan kemarin ada sodara di kampung yg minta di gambarin rumah tinggal. karena kelas rumah biasa/sederhana, banyak juga keterbatasan jadinya ngga bisa pake jasa arsitek beneran. yah itulah kenyataan.. keinginan juga baru diangan budget belum full, rencana bertahap, dll. Berbekal dari ngerjain tugas SBG1&2 dulu waktu kuliah dan my 2nd hobby gbr 3d realistik sya coba aja.. memang sudah beberapa kali sih sya bantuin sodara/temen yg mau bikin rumah, sya menyisihkan sedikit waktu luang buat bantuin.



sya hanya berpikir saat itu (skul) ngga bisa bantuin nyumbang matrial ya gppa dgn bantu desain aja sebisanya, termasuk penggunaan bahan yg dpat menekan biaya.



pake free software yg bagus efektif modelin 3d ngga menyita waktu, trus yg good juga di rendering. gppa emang khan bukan buat tujuan komersil. thanx Google SketchUp :)

Read More

tambahan modul pd program glassboxDI #2

Saat ini program glassboxDI mencapai versi 0.15 (alpha) dengan adanya tambahan modul pembuatan titik2 diagram interaksi kolom beton bertulang. Cukup banyak baris code yg dibuat, karena program ini bertujuan membuat sejenis  laporan desain yg transparan sebagai pengingat dan kontrol desain. Pembuatan diagram interaksi dipisah berdasarkan jumlah lapis tulangan, untuk yg sudah dibuat dibawah ini yaitu lapis-2 jumlah sampai ~1.000 list codes. Jadi jika akan dibuat sampai maks 10 lapis tul, maka akan mencapai > 10.000 list codes whohah jumlah yg tidak sedikit dan terlihat algorithma yg panjang.



Berikut adalah perbandingan dgn buku rujukan, ditambahkan titik tambahan diagram interaksi 1bh masing-masing pada daerah tekan menentukan dan tarik menentukan. Selain itu juga ditambahkan deteksi/perhitungan faktor aman terhadap beban P_u dan M_u yg bekerja.

[scribd id=24551922 key=key-1afzwg5tn23t64p6yps9]

Dibawah ini perbandingan dgn hasil program PCACOL,



Perbedaan disebabkan karena pembulatan atau lainnnya, terlihat pada keadaan tekan menentukan c = 350 mm > c_b. perlu dicheck ulang dgn berbagai konfigurasi jumlah tulangan.

Note:

ada kekeliruan (bugs) pada penentuan nilai c berikut snapshoot perbaikannya.



Perubahan notasi gaya aksial dan momen tersedia.

.

[scribd id=24579590 key=key-15xkkwm5zlvrjt97i044]

.

Berikut perbandingan lainnya,



Tambahan modul lainnya adalah koefisien momen pelat 2-arah kondisi tumpuan ideal, diambil dari tabel Marcus rujukan dari BS.



Selain itu juga sedang dilakukan proses digitizer data pada tabel koefisien momen dari M. Pigeaud rujukannya juga dari BS.



.

Read More