Materyel Mekaniği - 2
Kirişin Yatay Kesme Stresi
Yatay kesme stresinin mevcut olduğunun belki de en iyi ispatı alttaki şekli göstermek. İki tane tıpatıp aynı kirişi üst üste koysak ve üstten bir P yükü uygulasak, eğer kirişler arası sürtünme çok az ise kirişler birbirinden ayrı şekilde büküleceklerdi, (b) şeklinde görüldüğü gibi.
Şimdi hayal edelim ki kirişler birbirine tutkal ile yapıştırıldı, bu şekilde iki parça tek bir parça haline geldi. Bu birleşik kiriş yüklendiğinde o yapıştırılan yatay yüzeyde stresler oluşmalıdır ki yapıştırılmış yüzeyin üst (b) şekildeki gibi kayması engellensin. Bu kesme stresleri sayesinde tek birleşik kiriş ayrı ayrı iki kirişten daha katıdır / serttir (stiff). Aynı kavram genel olarak moleküler yapışma ile tek parça nesnelere de uygulayanabilirdi. Bir vidaya dikey yönde uygulanan yük kesme yönünde stres yaratır, ki bu stresler moleküler bağlantılar üzerinden ortaya çıkar.
Türetmeye gelelim [2, sf. 388]. Birörnek bükülmeye maruz olan bir kirişi düşünelim, yine üstteki (a) resmini referans alıyoruz, bu kirişte dx genişliğindeki ufak bir parçaya odaklanalım, bu ufak parçanın sol kısmına etki eden kesme kuvveti ve bükülme momenti V ve M. Eksende sağa gittikçe bu değerler değişik olabileceği için parçanın sağında V+dV ve M+dM olacaktır.
σ1 ve σ2 formüllerini yazalım, bükülme normal stres (flexure) formülünden σ=My/I olduğunu biliyoruz,
σ1=−MyI,σ2=(M+dM)yI
Üstteki (b) resminde daha da ufak bir bölgeye odaklanalım, alttaki (c) resminde daha net gösteriliyor,
Bu parçanın en üst kısmı kirişin en üstü, orada yatay kesme stresi yok. Parçanın alt kısmındaki strese τ diyelim. (c) figürünün sol taraftan bakılan hali (d) resminde. Kuvvetleri düşünürsek, F1,F2,F3 diyelim alttaki gibi olur.
Denge açısından
F3=F2−F1
olmalıdır.
Bu kuvvetleri yerine koyalım o zaman, mesela σ1'den hareketle,
σ1dA=MyIdA
Tüm F1 kuvveti için tanımladığımız alan üzerinden entegral alalım,
F1=∫σ1dA=∫MyIdA
F2,F3 için benzer şekilde,
F2=∫σ2dA=∫(M+dM)yIdA
F3=∫(M+dM)yIdA−∫MyIdA=∫(dM)yIdA
F3=dMI∫ydA
F3'e farklı bir açıdan yaklaşalım, eğer b boyunca kesme stresi τ değişmiyor ise, o zaman F3'ü alttaki gibi de belirtebilirdik,
F3=τbdx
Kuvvet eşittir stres çarpı alandır, bdx ile belirtilen alan üstteki resimde görülüyor, kirişin altındaki eni b boyu dx olan bölgeden bahsediyoruz. Burada τ sabit ise üstteki çarpım yapılabilir, tabii ki τ büyüklüğü y'ye bağlı olduğu için aşağı, yukarı değişimde τ değişirdi.
Devam edelim, son iki formülü birleştirince
dMI∫ydA=τbdx
Tekrar düzenleyince
τ=dMdx(1Ib)∫ydA
dM/dx büyüklüğü kesme kuvveti V'ye eşittir.
τ=VIb∫ydA
∫ydA entegrali Q ile gösterilir, alansal bir momenttir, yine kiriş yan yüzey şekli ile alakalı, standart şekiller için bilinen formüller vardır, o zaman nihai formül
τ=VQIb
Kirişin Yatay Kesme Stresi - Alternatif Anlatım
Kesme stresi τ'yu bulmak için yine kirişin ufak bir kısmına odaklanalım,
Tüm etki eden kuvvetlerin toplamı sıfır olmak zorundadır [3],
−P+(P+dP)+τbdx=0
ki b kesme stresinin uygulandığı noktadaki kiriş derinliğidir.
−dP/dx=τb(1)
P'yi bulmak için A bölgesindeki stresleri entegre ediyoruz,
∫AdP=∫AσbdA
Fakat daha önce bulduk ki σb=−My/I, yerine koyunca,
P=∫A−MyIdA
M ve I sabittir, entegral dışına çıkartılabilir,
P=−MI∫AydA=−MQI
ki görülen entegral bir alanın kütle merkezini bulmak için kullanılan standart bir entegraldir, Q=∫AydA. Devam edelim üstte bulunan P'yi (1)'e sokunca,
−ddx(−MQI)=τb
QIdMdx=τb
Şimdi hatırlarsak dM/dx türevi yatay kesme / teğetsel yükü V'ye eşittir. O zaman
QIV=τb
Nihai yatay kesme stres denklemi,
τ=VQIb
Problemler
Altta kiriş odaklı bazı örnek problemleri çözeceğiz. Bir kirişe yük uygulandığında dengenin muhafaza edilmesi için kiriş içinde kuvvetler oluşur. Bu iç kuvvetler kirişin destek yapısına göre farklı şekillerde ortaya çıkabilir [1].
Üstteki soldaki resimde mesela iki boyutta pimli destek dönüşe izin verir, tekerlekli yatay sağ, sol hareketi ve dönüşü serbest bırakır. Sabit destekte hiç harekete izin yoktur. Hangi harekete izin verilmediğine göre yük uygulanması ardından üst sağdaki iç kuvvetler ortaya çıkacaktır, bunlar pimli durumda dikey ve yatay kuvvetler, tekerlekli durumda dikey kuvvet, sabit durumda ise her üç mümkün tepkilerdir, yani moment, dikey ve yatay.
Yükler noktasal ya da dağıtık şekilde uygulanabilir, altta noktasal kuvvet, dağıtık kuvvet ve noktasal moment örneklerini görüyoruz.
Tipik olarak problemin beklediği kesim kuvveti ve bükülme momenti grafikleridir, bu grafiklerde x ekseni yatay olarak kirişin kendisi, y ekseni ise o noktada etki eden kesim ya da moment büyüklüğüdür.
Çözme yöntemi olarak iki yaklaşım mevcut, biri her kritik noktada kirişin hissettiği içsel kuvvetler ve momentleri hesaplamak için o noktalarda denge denklemlerini kullanmak, ki bu denklemlere (ve temel fiziğe göre) kirişe uyguladığımız hayali bir kesitte etki eden tüm kuvvetler ve momentler birbirini dengelemeli. Ardından bu kesit tüm kiriş boyunca kaydırılır ve gereken kuvvetler aynı denge üzerinden hesaplanır. Eğer tüm yükler noktasal ise bu yaklaşım iyi işler.
Bir diğer yaklaşım Calculus kullanmak. Bu yaklaşım temelde sürekli bazda çözüm verdiği için dağıtık yük durumunda daha kolay işler, kritik noktalara odaklanmak yerine pür formulsel düşünebiliriz . Daha önce görmüştük ki kesim kuvvet formülünün eğimi (türevi) o noktadaki uygulanan yükün negatifidir dV/dx=−w, ve bükülme moment grafiğinin eğimi ise o noktadaki kesim kuvvetine eşittir, dM/dx=V. Şimdi ters yönde gidersek, ilk formülü entegre edince mesela iki nokta arasındaki kesme kuvveti farkının yükleme eğrisinin altında kalan aynı noktalar arasındaki alanın olduğunu görebilirdik.
V2−V1=−∫x2x1wdx
İkinci formülü entegre edince iki nokta arasındaki bükülme moment farkının kesme kuvveti eğrisinin altında kalan alanın olduğunu bulurduk [1].
M2−M1=∫x2x1Vdxa
Problem 1
Üstüne birörnek yük q uygulanan bir dirsekli (cantilever) kiriş için kesme kuvveti ve bükülme momenti diyagramı çizin [2, sf. 334].
Çözüm
V−VA=V−0=V=−∫x0qdx=−qx
M−MA=M−0=M=∫x0−qxdx=−qx22
Problem 2
Alttaki kiriş için kesme ve bükülme moment diyagramlarını oluşturun [4, sf. 179].
Çözüm
Önce hareketsizliğin olduğu noktalardaki reaksiyon kuvvetlerini bulmak gerekir. Bu kuvvetler denge noktasında dikey ve yatay kuvvetlerin toplamının sıfır olma zorunluluğu üzerinden hesaplanabilirler.
Fakat bazı durumlarda zorluk çıkabilir, mesela eğer üstteki örnekte tek bir reaksiyon noktası olsa aşağı yöndeki kuvvetleri toplayıp tek reaksiyona eşitleyebilirdik. Fakat iki reaksiyon noktası var, etki iki tane reaksiyona nasıl yayılır? Bu noktada bükülme momentlerini dahil ederek ortaya bir formül çıkartabiliriz. Reaksiyonlar bilinmiyor olsa da birer kuvvettirler, ve bu kuvvetler, yük kuvvetleri ile beraber diğer sabit noktalar üzerinde bir moment oluştururlar. Genel bir figür üzerinde görelim [4, sf. 86],
Reaksiyonlar Ay ve By nasıl hesaplanır? Moment dengesi MA=0 ve MB=0 olmalı,
MA=0=ByS−Pa→By=PaS
MB=0=−AyS+Pb→Ay=PbS
Burada gördüğümüz a/S ve b/S basit birer orantıdan ibaret, yani tek P yükünü mesafe oranlarına göre iki reaksiyon arasında paylaştırabiliriz.
Problemimiz icin
Ay=60(6/9)+150(3/9)=90kN
Bir tane reaksiyonu bulunca aslında diğerini dikey kuvvetlerin toplamının sıfır olma şartından elde edebiliriz,
90−60−150+Dy=0
Dy=120kN
Kesme diyagramı için kesme kuvvetinin yükün entegrali olduğunu hatırlayalım. Entegral demek alan demektir, bu problemde elimizde bir analitik sürekli fonksiyon yok, bazı yerlerde noktasal yükler var. Eğer bir alan olarak yük uygulaması var ise, bu durumda dikdörtgensel alanları hesaplamak yeterli olurdu (bir problemde hem noktasal, hem alansal yük uygulaması aynı anda olabilir).
A üzerinde noktasal +90 kN yük var, kesme kuvvetine ne olur? Burada noktasal veri ile çalışmamız gerekiyor. Tam yükün uygulandığı noktada çok ufak bir Δx alanına odaklanalım, bu ufak parçanın serbest-kütle (free-body) diyagramı şuna benzer [5, sf. 208], alt solda,
Dikey yöndeki kuvvetleri toplayıp sıfıra eşitlersek, ve ΔV için çözersek,
P0+V−(V+ΔV)=0→ΔV=P0
Yük noktasında kesme kuvveti farkı yüke eşitmiş, o zaman kesme kuvveti 90 kN'e zıplar.
Daha ilerlemeden bir noktasal bükülme momenti uygulansa ne olur onu da görelim. Yine üstteki resimde sağdaki figüre bakalım, buradaki ufak parçanın serbest-gövde diyagramında tüm momentlerin toplamını sıfıra eşitleyelim, etrafında dönme olan referans noktası O olsun, ve ΔM için çözelim,
M0−M−VΔx+M+ΔM=0
ΔM=M0+VΔx
Δx sıfıra yaklaştıkça (noktasal bu demek) üstteki,
ΔM=M0
eşitliğine yaklaşır.
Problemimizle devam edelim, 90 kN'e zıplama oldu, kesme kuvveti grafiği o noktaya çıkar, B noktasına kadar düz devam eder, 60 kN ters yönde yük alınca 30 kN'ye düşüş olur, C'ye gelince yine aşağı yönde 150 kN yük alınca -120'ye düşüş olur. Yani grafik şuna benzer,
Moment grafiğine geldik, analitik versiyondaki alan hesabı burada ayrıksal bağlamda uygulanabilir, üstteki figürde alanların içine büyüklüklerini yazdık, 270, 90, -360 diye. Tabii her noktada o noktanın solundaki alan dikdörtgen toplamları hesaba alınmalı, bu da bize alttaki şekli verir,
Problem 3
Alttaki destekli kiriş için kesme kuvveti ve bükülme moment diyagramlarını çizin. Maksimum bükülme momentinin hangi noktada olduğunu bulun.
Reaksiyonları bulalım, hareketi engelleyen bir nokta var, kirişin duvara bağlandığı yer, orada dikey ve bükülme engelleniyor, o zaman iki tane reaksiyon ortaya çıkar, Ay ve My diyelim, denge denklemleri üzerinden,
∑Fy=Ay+180kN−50kN=0→Ay=−130kN
Bükülme moment dengesi üzerinden,
∑MA=(180kN)(1.5m)−(50kN)(5m)−140kN⋅m−MA=0→MA=−120kN
Kesme kuvveti diyagramını alttaki gibi çizeriz,
Daha önce türettiğimiz ΔV=Ay eşitliği üzerinden grafik hemen -130 kN'ye düşecektir. Dağıtık yükün altındaki alanın, entegralinin kesme kuvvetine eşit olması bilgisini kullanacağız şimdi, bu örnekte analitik formül yok, alan hesabını basit dikdörtgen hesabı ile yapabiliriz, başlangıç kesme kuvvetinden 3 metre sonra 60 kN/m çarpı 3 m = 180 kN eklenmiş olacak, bu bizi 50 kN seviyesinde getirir. Bu geçiş lineer olur, bir eğim üzerinden, hesaplanan alan sağa gittikçe gittikçe büyüyen dikdörtgenler vasıtasıyla. Devam ediyoruz, D noktasına gelene kadar değişim yok, o noktada ters yönde 50 kN var, böylece kesme kuvveti 0 seviyesine iniyor.
g noktasını bulmak için -130 kN'i nötralize etmek için kaç metre boyunca 60 kN/m dağıtık yük eklemek gerekir sorusunu cevaplamak lazım, basit bir bölme işlemi 130 / 60 = 2.1668 metre.
Bükülme moment diyagramı kesme kuvveti diyagramının alan hesabı üzerinden yapılacak,
Başlangıç bükülme momenti ilk anda −120kN⋅m seviyesine düşer, çünkü orada noktasal bir moment var. Daha sonraki eklenen moment için son iki diyagramdaki g−j eşlemesine bakalım, oraya gelirken bir üçgen alanı var (yeşil renkli) bu alan -130 kN çarpı 2.1667 m bolu 2 ile hesaplanır, ve -140.8355 sonucunu verir. Bu alanı -120 kN⋅m'ye ekleyince -260.8355 kN⋅m elde edilir. Dikkat oraya giderken ortaya çıkan eğri bir parabölün parçası, hatırlarsak y=x'in entegrali x2, yani bir parabol.
Geri kalan noktalar benzer şekilde, l−m geçişine dikkat, orada bir 140 kN⋅m büyüklüğünde noktasal moment uygulandığı için direk bir zıplama var, gerisi lineer şekilde devam ediyor.
Kaynaklar
[1] The Efficient Engineer, Understanding Shear Force and Bending Moment Diagrams, https://youtu.be/C-FEVzI8oe8
[2] Gere, Mechanics of Materials, 7th Edition
[3] Gramoll, Mechanics, http://www.ecourses.ou.edu/cgi-bin/ebook.cgi?topic=me
[4] Kassimali, Structural Analysis, 6th Ed
[5] Philpot, Mechanics of Materials
Yukarı