3.8 Devre Kapalılık Denklemlerinin Sayısal Çözümü Elde etmiş olduğumuz denklemler ve benzer yöntemlerle başka mekanizmalar için elde edeceğimiz denklemler bir mekanizmanın her hangi bir konumunu veya mekanizmanın olası tüm konumlarını belirlememizi sağlar. Bunun için uzuv boyutları bilindiğinde bu denklemlerin sayısal olarak çözülmesi gereklidir. Mühendislikte analitik çalışmalar sonucunda elde edilenler sayısal değerlere dönüştürülmeden yapılanların uygulamaya indirilebilmesi mümkün değildir. Mekanizmaların sayısal çözümü uzun yıllar çizim masasında pergel ve cetvel kullanılarak geometrik yöntemlerle yapılmıştır. Bu yöntemin olumlu yanı kişinin yapmakta olduğu işi çok iyi canlandırmasına olanak tanımasıdır. Benzer geometrik çözüm bilgisayarlarda günümüzde AUTOCAD®, CADKEY® ve benzeri çizim programları kullanılarak yapılabilir. Geometrik yöntem açıklanırken bu tür çözümün örnekleri verilmiştir. Ancak geometrik yöntemler çok zaman aldığı gibi, hassasiyetleri sınırlı olup uygulayan kişinin becerilerine çok bağlı olduğundan mühendislikte sayısal çözüme daima ihtiyaç duyulmuştur. Bilgisayarların uygulamaya girdiği ilk yıllarda mekanizmaların bilgisayar destekli analizi ve sentezi söz konusu olmuştur. 1950'li yıllardan itibaren ilk olarak elde edilmiş olan analitik denklemler ve algoritma FORTRAN programlama dili ile bilgisayarlara uyarlanmıştır. Bilgisayarlar teknolojisinde gelişmelerle birlikte BASIC, PASCAL ve C lisanlarıda bu algoritmaların bilgisayarlarda uyarlanmasında kullanılmıştır. Bu uygulamaların sonucu ne yazıkki görsel açıdan çok zayıf olup genellikle çok sayıda rakamlardan oluşan uzun bilgisayar çıktıları olmuşdur. Bu nedenle işin uzmanları tarafından çok özel uygulamalar için kullanılabilmiştir. Bu bilgisayar lisanlarına günümüzde grafik komutlar kazandırılmış olmalarına rağmen, genellikle grafik uygulamalar için grafik özellikleri daha uygun ve kullanıcı ile etkileşimli program yazılmasına uygun olanaklar ve ortam oluşturabilen VisualBasic, VisualC++ ve Delphi gibi program yazma platformları geliştirilmiştir. Kişisel bilgisayarlarla birlikte bilgisayar kullanımının daha yaygınlaşabilmesi için kullanıcılara öğrenilmesi programlama dillerinden nisbeten daha kolay, genel amaçlı kullanıcı ile iletişimli paket programlar geliştirilmiştir. En yaygın olan paket programlar bilgisayarın bir yazıcı olarak çalışmasını sağlayan paket programlardır. Daktilodan çok daha farklı olarak, değişik yazı karakterleri ve boyutları kullanmamıza olanak tanıyan, satır sonlarını düzenleyen bu paket programlar günümüzde daktilolardan daha fazla kullanılmaktadır. İkinci yaygın paket program ise, çarşaf liste hazırlamak için kullanılan programlardır. Bu programlar ilk olarak muhasebe işlemlerini kolaylaştırmak için düşünülmüştür. Örneğin masrafları gösteren bir kolondaki sayılar kolayca toplanabilir veya maliyeti gösteren bir kolon belirli bir genel gider yüzdesi ile çarpıldıktan sonra elde edilmesi hedeflenen yüzde kar ile çarpılarak bir başka kolona satış fiyatı kolonu olarak elde edilmesi açısından bu tür liste programları çok kullanılışlı olmuştur. Ancak zamanla basit çarpma, bölme toplama gibi hesaplamalardan çok daha fazla fonksiyonları (her türlü istatistik işlemleri, trigonometrik fonksiyonları) içerir hale getirilen bu programlar günümüzde mühendisler için her türlü analizi yapabilecek, sonucu grafik olarak gösterebilecek vazgeçilemiyen bir paket olmuştur. Günümüzde Microsoft Excel®, Lotus 123® ve Corel Quattro Pro® değişik firmalar tarafından üretilen üç önemli çarşaf listeleme paket programlarıdır. Günümüzde bu paket programların çok yaygın olarak bilgisayarlarda olması, bu paket programların her mühendis tarafından öğrenilmesini gerekli kılmaktadır. Bilgisayar alanında gelişmelerle birlikte mühendislik, matematik ve fen bilimlerinde uğraşan kişiler için son yıllarda her türlü matematiksel işlemi yapabilecek matematik uygulama programlarının geliştirilmesidir. Maple®, Mathematica®, MATLAB®, Mathcad®, SCILAB®, TKsolver®, Eureka® bunlardan önemli olanlardır. Mekanizmaların kinematik ve dinamik analizi için ayrıca özel paket programlar geliştirilmiştir (ADAMS®, WorkingModel®, vb ).. Bu paket programlar kullanılarak düzlemsel veya uzaysal mekanizmaların kinematik ve dinamik analizi mümkündür. Ancak konunun temeli öğrenilirken bu paket programların kullanılması sakıncalı olacaktır. Konunun uzmanı olunduktan sonra bu programlar zaman kazandırdığından dolayı kullanılabilir. Analizi yapılabilecek mekanizmaların burada vereceğimiz örneklerle sınırlı olması düşünülemez. Ayrıca her geçen gün farklı çözüm imkanlarının ortaya çıktığıda bir gerçektir. Ancak temel matematiksel yaklaşımlar mekanizmaların temel geometrilerine bağlı kaldığından, değişmesi beklenmemelidir. Burada amaç bilgisayar programlamayı öğretmek veya bir paket programı kullanmayı öğretmek olmadığından anlatılacak olan daha çok söz konusu mekanizma analizinin bilgisayarlarda sayısal yöntemlerle en kolay bir şekilde nasıl yapılacağını açıklamaktır. Bunun için Excel ve Mathcad paket programları kullanılarak yapılan çözüm örnekleri verilecektir. Kullanılacak olan denklemler ve algoritma aynı olacağından açıklanmış olan yöntem kolaylıkla başka bir paket program kullanarak veya bir bilgisayar lisanı kullanılarak o problem için özel bir program yazılarak çözülebilir. Bilgisayar kullanmayı çok iyi bilenler, burada analtılanlardan çok daha fazlasını daha güzel bir görünüm verecek şekilde yapılabileceğini kolaylıkla söyleyebilirler. Örnek Gösterilen krank-biyel mekanizmasında uzuv boyutları: C noktasının koordinatlarını Krank Biyel mekanizması için elde ettiğimiz denklemler yeniden yazıldığında:
verilmiş olan sayısal değerler kullanıldığında:
Çözüm
C noktasının konum vektörü A0C = A0B + BC olarak yazılır, x ve y bileşenleri ayrı ayrı ele alınır ise:
Hesap makinası programlanabilir ise, bir defa yapılan bu işlemler makina
hafızasına alınarak değişik Örnek
Yukarıda gösterilen dört-çubuk mekanizmasında uzuv boyutları: a1=70 mm; a2= 35 mm; a3= 62.3 mm; a4= 56 mm; BC= c3 = 84.1 mm; AC=b3=126.6 mm dir. 0o<
Bu denklemde
Bu örneği EXCEL® paket programını kullanarak çözelim. Excel bir genel tablolama programıdır. Genel olarak tablolama programlarında çok sayıda hücreler bulunmaktadır. Hücreler rakamla belirtilen yatay satırlar ve harflerle gösterilen dikey kolonlar üzerinde bulunduğundan her hangi bir hücre kolon harfi ve satır rakamı ile gösterilir. Yani her bir hücrenin bir adı vardır. C25 (C kolonunda, 25inci satırda bulunan hücre ), A1 (birinci satır birinci kolon) AH55… gibi. Bir hücrenin adını isterseniz sizde etiketleme yaparak değiştirebilirsiniz. Bu hücreleri içine her hangi bir rakam, veya isim yerleştirebiliriz. Bu şekilde isterseniz güzel görünümlü tablolar veya listeler oluşturabiliriz. Bu tür paket programların mühendislik açısından en önemli avantajı, her türlü matematiksel kuralları kullanarak hücrelere yazılmış olan rakamlar üzerinde gerekli işlemlerin yapabilmesidir. Basit dört işlemin yanında bu programlarda her türlü trigonemetrik fonksiyon, istatistik fonksiyonları (maksimum, minimum, ortalama bulma,…) ve benzer fonksiyonlar bulunduğu gibi matris inversi alınması, matris çarpımı veya bir lineer olmayan denklemin köklerinin bulunması dahil olmak üzere çeşitli işlemler basit komutlarla yapılabilir. Gerektiğinde paket programda bulunmayan ancak kullanıcının devamlı kullanacağı fonksiyonlar ise MACRO olarak adlandırılan işlem ile veya VISUAL BASIC lisanı kullanarak yazılıp eklenebilir ve kullanılabilecek fonksi-yon sayısı artırabilir. Örneğin bir hücredeki sayının karesini alıp diğer bir hücredeki sayının kosinüsü ile çarpıp bir başka hücreye bu değerin kare kökünü yazdırabilirsiniz. Bu özellikler mekanizmaların analizi için yeterlidir.
Sabit uzuv boyutları a1, a2,
a3, a4
değerlerini A1 den A4 e kadar olan hücrelere yerleştirelim. A7 hücresine
Şimdi A8 hücresine farklı bir Sonuç değerler kolonlar şeklinde elde edildi ise de, bu değerlerin anlaşılması
oldukça zor olabilir. Hareketi gözlemliyebilmek için Excel paket programının
Chart menusunda bulunan grafik çizme özellikleri kullanılabilir. Aşağıda,
Excel paket programının bu özellikleri kullanılarak bağımlı konum değişkenleri
Krank Açısını göre çıkış kolu, biyel ve bağlama açıları
C biyel noktası eğrisi Sayfamızda yapmak istediğimiz işlemler girildikten sonra çok önemli ve farklı bir konuya girme imkanımızda olacaktır. Mekanizmanın analizi için baştan verildiğini kabul ettiğimiz uzuv boyutlarını sorgulamamız mümkün olacaktır. Örnek olarak 70 mm olarak verilmiş olan a1 boyutunu ele alırsak, acaba bu boyut 65 mm veya 60 mm olur ise sonuçları nasıl etkiler diye bir araştırmaya girilebilecektir. Benzer sorgulama her uzuv için yapılabilir ve sonuç bütün farklı değerler için anında Excel sayfamızda hesaplanarak görülecektir ve grafikler yeni duruma göre anında değişecektir. Yukarıda gösterilen işlemlerin Excel kütüğü için (four-bar.xls) tıklayınız. Bu tür paket programlarda eğer giriş parametresine göre diğer değerleri elde etmek için uygun ve doğru bir algoritma geliştirmiş isek, her hangi bir formülasyon yöntemini kullanabiliriz. Örnek Şekilde gösterilmiş olan mekanizma iki mekanizmanın seri bağlantısı olarak düşünülebilir. Birinci mekanizma (A0AB0, 1,2 ve 3 uzuvları), salınan kol-kızak mekanizması olup 2 krankının dönmesi ile 3 kolu salınım yapacaktır. İkinci mekanizma ise bir dört-çubuk mekanizması olup(B0BCC0 , 1,3,4 ve 5 uzuvları) 3 giriş uzvunun salınımı, 5 uzvunda farklı bir salınıma dönüşecektir. Önceden görüldüğü gibi, uygulamada mekanizmalar çoğunlukla belirli basit mekanizmaların seri bağlantısı olup bir basit mekanizmanın çıkış uzvu diğer mekanizmanın çıkış uzvu olmaktadır. Temel olarak karşılaşacağımız basit mekanizmalar dört-çubuk, kol-kızak, krank-biyel mekanizmalarıdır. Bu basit mekanizmaları tek bir fonksiyon kullanarak çözebilmek için Excel içinde program yazmamız ile mümkündür. Bunun için TOOLS menüsünden Macro seçilir ve VisualBasic Editörüne girilebilir ve VisualBasic programlama komutları kullanılarak Çeşitli fonksiyonlar ve alt programlar yazılabilir. Yazılan bu programlar "FILE- EXPORT"komutu ile diskte saklanabildiği gibi, başka bir Excel kütüğünde kullanmak istediğimizde "FILE - IMPORT" komutu ile çağırılabilir. Bu işlem basit mekanizmalar için bir defa hazırlanıp bir kütüphane oluşturulur ise, mekanizma analizi yazılmış olan bu fonksiyonlar ile bir iki komuta indirgenebilecektir. - EK 2 için tıklayınız - Ek 2 de Visual
Basic® ile hazırlanmış, Excel paket programı veya VisualBasic
programlama lisanı içinde kullanılabilecek temel bazı fonksiyonlar verilmektedir.
Bu örnekte verilmiş olan mekanizma için bu hazır programlar kullanılacaktır.
Bir önceki örnekten farklı olarak hücre adlarını kolon harfi ve sıra sayısı
olarak kullanmaktansa, bilhassa sabit uzuv boyutlarını yerleştirdiğimiz
hücrelere isimler verelim. Bu isimler aşağıda verilen excel tablosunda
gösterilmektedir. Sabit uzuv değerleri bu gösterilen hücrelere girildikten
sonra,
İsim verilen hücreler
Grafik çıktı Yukarıda gösterilen işlemlerin Excel kütüğü için (fivelink.xls) tıklayınız.
Örnek Bu örnekte MathCad® paket programı kullanılacaktır. Bu
çeşit programlar her türlü matematik denklemi (sayısal veya sembolik olarak)
çözmek için geliştirilmişlerdir. Excel gibi programlarda genel olarak
sonucun yer aldığı veya kendi yerleştirdiğimiz değerlerin olduğu hücrelerin
her biri ayrı ayrı görülmekte, bir hücreye koyduğumuz değerleri oluştururken
yaptığımız işlemler perde arkasında kalmaktadır. Mathcad türü programlarda
ise bu hücreleri kolon harfi ve satır sayısı ile adlandırmaktansa biz
hücre kapağına bir etiket yerleştirip üstüne bu hücreyi temsil edecek
istediğimiz herhangi bir isim koyabiliriz. Bu durumda o hücrenin içinde
ne olduğunu etiketteki bu isimden biliriz. Ancak hücrenin kapağını açmadan
içinde ne kadar olduğunu bilemeyiz. Yani Excelde kapaksız olan hücrelerin
her birine kapak ve kapakların üzerinede birer etiket konmuştur (Excelde
bir hücreye ad verseniz bile etiketin üzerinde bulunan adı ancak imleci
hücreye götürdüğünüz zaman görebilirsiniz). MathCadde örneğin eğer : a:=120
, b:=50, Denklemini görürsünüz. Bu işlemle, etiketinde xa yazılı yeni bir hücre
oluşturulmuş ve denklemin sağında yazılmış olan operasyonlar ile a,b ve
Yukarıda ayarlı bir pompa mekanizması görülmektedir. Eksantrik olarak tasarlanmış olan krank mili bir elektrik motorundan sonsuz vida - dişli sistemi ile devir sayısı gerektiği kadar düşürüldükten sonra döndürülmektedir. Mekanizmanın şematik görünümü Şekil yukarıda gösterildiği gibidir. Pompanın kurs boyu ayarlaması B0 mafsal noktasının bir ayar vidası ile yukarı aşağı hareket ettirilerek, s1 uzunluğu değiştirilerek, yapılmaktadır. A0ABB0 devresi (1,2,3 ve 4 uzuvları) bir dört çubuk mekanizmasını oluşturmaktadır. Önceki örneklerden farklı olarak A0B0 sabit mafsal noktaları yatay değildirler, B0 noktası A0 noktasının solunda ve yukarısındadır. Yine de vektör devre denklemi (A0A + AB = A0B0 + B0B) aynı şekilde yazılabilir: Vektör denklemi B0A
= B0A0
+ A0A eğer reel ve sanal kısımlar ayrı ayrı eşitlenir ise: s cos Önceden gösterildiği gibi s ve Pompanın ikinci kısmını oluşturan ve yukarıda gösterilen çift kızaklı mekanizma için (1,4 ve 5 uzuvları), vektör devre denklemi kompleks sayılarla: c4
MathCAD ile çözüm yapılırken bu denklemleri aynen yazmak yeterlidir.
Ancak her hangi bir işlem yapmadan önce sayısal parametrelerin (bağımsız
değişken veya sabit parametrelerin) değerleri mutlaka programda yazılmalıdır.
Program tanımlanmamış veya içi boş hücreler ile işlem yapamaz. MathCad Sayfası - 1 Verilen bir krank açısına göre çıkış uzvu konumu: A0A = a2:=55 AB = a3:=240 BB0 = a4:=165 b1:=185 s1:=90 c4:=70 xa
:= b1 + a2
.cos(
Eğer hareket analizini mekanizmanın tüm konumları için yapmak ister isek;
örneğin krank açısının her 50 aralığında
00 dan 3600
aralıkda değerlerin hesaplanması gerekli ise, bu durumda hücrelere yazacağımız
etiketlerde isim bulmakta zorluk çekebileceğimizden indis kullanmamız
gerekli olacaktır. Oluşacak 73 krank açısının her birini MathCad Sayfası - 2 Tüm devir analizi: Yukarıda s15 in krank açısına göre
değişimi görülmektedir. Ancak s1 mesafesi
bir ayar vidası ile değiştirilmektedir. Eğer bu mesafenin 5 uzvu hareketine
etkisi incelenmek istenir ise, s15 konum
parametresinin s1 e göre nasıl değiştiği
belirlenmelidir. Bunun için değişkenlerde iki indis birlikte kullanılmalıdır
(bir indis değişik MathCad Sayfası - 3 Tüm bir devir ve değişik ayar durumları için analiz Görüldüğü gibi, s1 mesafesi azaldıkça, piston stroku azalacak, arttıkca strok artacaktır.
Örnek Bir önceki MathCad ile çözülmüş olan denklemler Excel paket programı
kullanılarak aynen çözülebilir. Tabii ara basamaklar daha fazla olduğundan,
kullanılacak olan kolon sayısı fazla olur. Normal Excel Komutları yerine
bu işlemlerin tümünü yapan, verilen Global Const PI = 3.1415926 AyarlıPompa() fonksiyonundaki denklemler, analiz sonucunda elde ettiğimiz denklemlerin BASIC lisanı ile yazılımıdır. Ayrıca bu fonksiyon içinde, Boy(X,Y) (X ve Y boyutları verildiğinde hipotenüs bulma), Açı(X,Y) (X ve Y değerleri verildiğinde açı bulma), AçıCos(U1,U2,U3) (Kosinüs teoremine göre üç kenarı verilen bir üçgenin U3 kenarının karşısındaki iç açıyı bulma) gibi üç ayrı fonksiyonu kullanmaktadır. Kullanılan bu fonksiyonlar EK-2 de verilmiştir. Bu fonksiyon kullanılmadan önce sabit uzuv boyutlarını B1-B5 hücrelerine sırası ile yerleştirelim ve bu hücrelere sırası ile Krank, Biyel, Arakol, SabitU_Y, ve Eksantrik isimlerini verelim. B7 hücresine ise vidanın bir ayar boyutuna gelen s1 uzunluğunu yerleştirelim (s1=20mm). İkinci olarak A8 - A26 hücrelerine krank açısını her 200 için yerleştirelim (ilk iki hücre doldurulduktan sonra, fare ile sağ alt köşeden tutup, A26 hücresine kadar uzatılabilir. İmleci B8 hücresine getirelim. Basic lisanı ile yazmış olduğumuz fonksiyonu
kullanmak için bu hücreye = AyarlıPompa(Krank; Biyel; Arakol; SabitU_Y;
Eksantrik; B$7; $A8 ) komutunu tuşlayalım. Bu komut hazırlanmış olan
program devreye girer ve A8 hücresindeki AyarlıPompa fonksiyonu artık Excel de hazır herhangi bir fonksiyon programı gibi kullanıma hazırdır. s1 mesafesi bir ayar vidası ile değiştirilmektedir. Eğer bu ayarın 5 uzvunun hareketine etkisi incelenmek istenir ise, s15 konum parametresinin s1 e göre nasıl de-ğiştiği belirlemek gerekecektir. Bunun için B7, C7, D7, E7, F7 ve G7 hücrelerine değişik s1 değerleri yerleştirilmiştir. B8 hücresinde yazmış olduğumuz formülü şimdi B8-G27 arasında bulunan tüm hücrelere kopyalamamız gerekecektir. C kolonuna gelindiğinde s1 in bulunduğu hücreler B$7 şeklinde yazıldığından, B yerine C kolon harfi gelecek ancak hangi satırda olursak olalım, satır sayısının önüne $ işareti bulunduğundan (B$7) C kolonunun her hangi bir hücresine formül kopyalandığında, C7 hücresindeki s1 değeri kullanılacaktır. Benzer bir şekilde C kolonuna gelindiğinde krank açısı olarak $A8 yazıldığından krank açısı değeri mutlaka A kolonunda bulunan hücre olacak ve bulunduğumuz satır A kolonunda kullanılacak olan krank değerinin bulunduğu satır ile aynı olacaktır. Sonuç Şekil aşağıda de görülmektedir ve Mathcad ile yapılmış olan çözüm ile aynıdır. Excel kütüğü için (ayarlıpompa.xls) tıklayınız.
Hesaplama hassasiyeti ile ilgili not: Bilgisayarlarda sayılar belirli sayıda rakam ile depo edilir. Örneğin Excel de sayılar 15 hane olarak depolanır. Hane sayısı programa ve bilgisayara göre 6 ila 64 arasında olabilir. Mekanizmaların analizi veya sentezi sırasında bu hassasiyet gereğinden fazladır ve çok hassas sonuçlar, eğer doğru yöntem kullanılır ise, alınabilir. Bu önlemlerden birincisi yapılacak hesaplama miktarının minimuma indirilmesi ve bilhassa bölme işlemlerinde paydanın sıfıra yakın olmamasına dikkat edilmesidir. Ayrıca mekanizmanın özel konumlarında (ölü konumlar, kilitlenme konumları gibi uzuvlarda mafsalların birkaçının aynı doğru üzerinde olması durumu) hassasiyet azalabilir ve ayrıca bilgisayardan hata (sıfıra bölme veya negatif bir sayının kökü gibi) mesajları alınabilir. Yapılan işlemler doğru ise, Bu aynı zamanda o mekanizmanın pratikte sorunu olacağını da gösterir. Mühendislikte hiçbir sonuç başlangıç değerlerinden daha hassas olamaz. Bu nedenle elde edilen sonuç rakamlar sınırlı sayıda hane ile gösterilmelidir. Örneğin ±.1mm toleransla üretilmiş uzuvlardan oluşan bir mekanizmada sonuç s konum değişken değerini x=26.231456789mm olarak göstermek yanlıştır (mühendislik açısından mantıksızdır) x=26.2 ±.1mm den daha hassas elde edilmiş olamaz. Bu nedenle ekranda görülen rakamlar mutlaka kullanıcı tarafından incelenmelidir. |