Gazi Üniv. Müh. Mim. Fak. Der. J. Fac. Eng. Arch. Gazi Univ. Cilt 20, No 4, 433-442, 2005 Vol 20, No 4, 433-442, 2005 SABİT MIKNATISLI DEMİR NÜVELİ KÜRESEL EYLEYİCİNİN BİLGİSAYAR DESTEKLİ 3 BOYUTLU STATİK MANYETİK ANALİZİ VE UYGULAMASI Osman GÜRDAL ve Yusuf ÖNER* Elektrik Eğitimi Bölümü, Teknik Eğitim Fakültesi, Gazi Üniversitesi, Teknikokullar, 06500, Ankara *Elektrik ve Elektronik Mühendisliği Bölümü, Mühendislik Fakültesi, Pamukkale Üniversitesi, Denizli ogurdal@gazi.edu.tr, yoner@pamukkale.edu.tr (Geliş/Received: 23.07.2004; Kabul/Accepted: 20.09.2005) ÖZET Bu çalışmada, NdFeB sabit mıknatıs dilimlerinden yapılmış sabit mıknatıslı küresel rotor ve yumuşak demir ve basit sargı düzenlemesine sahip bir statordan oluşan küresel eyleyici tanımlanmaktadır. Küresel rotor düşük sürtünme katsayılı teflon yüzey kaplamalı bir stator içine yerleştirilmiştir. Stator üzerindeki üç ortogonal sargı (8 adet bağımsız enerjilendirilebilen sargı) halka akısının artırılması için bir dış küresel demir kabukla çevrelenmiştir. Stator demir nüvesi yokken, rotorun analitik olarak formülasyonu yapılmış manyetik alan dağılımı değerleri sonlu elemanlar analizi ve deneysel ölçümler ile doğrulanmış ve stator demir nüvesi varken rotorun ürettiği tork değerleri sonlu elemanlar metodu (Ansoft Maxwell 3D ile) ile tahmin edilerek deneysel olarak doğrulanmıştır. Tasarımı yapılan ve uygulanan sabit mıknatıslı demir nüveli küresel eyleyici açık çevrimli olarak 360° dönme hareketinin yanında ![]() Anahtar Kelimeler: Küresel eyleyici, sabit mıknatıs, 3 boyutlu manyetik analiz. COMPUTER AIDED THREE DIMENSIONAL MAGNETOSTATIC ANALYSIS AND APPLICATION OF IRON CORED PERMANENT MAGNET SPHERICAL ACTUATOR ABSTRACT In this study, a spherical actuator which has spherical permanent magnet rotor made of NdFeB slices and stator comprising mild steel and simple winding arrangement is described. The spherical rotor is housed within the spherical stator on a low friction constant teflon surface coating. Accommodated on the stator are three orthogonal windings (8 coils to be energized independently) enclosed by an outer spherical iron shell to increase the flux linkage. Without stator iron core, the values of magnetic field distribution formulated analytically are validated by finite element analysis and experimentally, and with stator iron core the torque values generated by rotor are validated by finite element method (by Ansoft Maxwell 3D) and experimentally. The permanent magnet rotor spherical actuator designed and applied is capable of ![]() Keywords: Spherical actuator, permanent magnet, 3D magnetic analysis.
Robotik, büro otomasyonu ve akıllı esnek üretim ve montaj sistemleri çok serbestlik dereceli küçük hassas eyleyicilerin geliştirilmesini zorunlu hale getirmektedir. Genelde çok serbestlik dereceli hareket her eksen için hemen hemen ayrı bir motor/eyleyicinin kullanılmasıyla gerçekleştirilir ki bu karmaşık geçiş sistemleri ve oldukça ağır yapılarla sonuçlanmaktadır. Bundan dolayı dişlilerin ataleti, geri kayması, doğrusal olmayan sürtünme ve elastik deformasyonundan dolayı yüksek bir dinamik performansın elde edilmesi zor olmaktadır. Üç serbestlik derecede kontrollü hareket yeteneğine sahip eyleyiciler bu problemlerin üstesinden gelebilmektedir. Böyle eyleyiciler mevcut durumda araştırma konusu olma özelliklerini sürdürmekle beraber, karmaşık yapıları ve elektromanyetik davranışlarının modellenmesi ve tasarımlarının optimize edilmesinde karşılaşılan zorluklardan dolayı ticari boyuta ulaşmakta zorlanmaktadırlar [1]. Yüksek performanslı hareket kontrol sistemleri için güç elektroniği bilimi ve dijital kontrol teknolojisi süregelen pazar talepleri sonucunda küresel eyleyici uygulamalarına daha fazla eğilimin olmasını sağlamıştır. Böylece çoğu küresel eyleyici tipleri ve topolojileri çok değişken işlevsel özellikli, bir eksen etrafında dönebilme veya doğrusal hareket yapabilme özelliklerine sahip olmuştur. Küresel eyleyiciler bir eksen etrafında dönebilen, doğrusal hareket yapabilen, hareketli sargısı, hareketli demir nüvesi veya hareketli daimi mıknatısı ile 3 boyutta sınırlı hareket yapabilen elektromekanik elemanlardır. Günümüze kadar farklı yapılarda küresel eyleyiciler tasarlanmış ve performans analizleri yapılmıştır. Williams ve Laithwaite’in 1950’lerde değişken hızlı indüksiyon motorlar üzerine çalışması, indüksiyon motorun küresel geometri ile geliştirilmesine yol açmıştır [2]. Williams ve Laithwaite’in bu çalışması üzerine Vachtsevanos ve arkadaşları robot bileği için üç serbestlik dereceli küresel indüksiyon motor kavramını ortaya atmışlar [3] ve detaylı bir analizle indüksiyon motor torku ve faydalı akısı arasındaki ilişkiyi incelemişlerdir [4]. İncelenen indüksiyon tipi küresel motor üç boyutlu uzayda basit bir tork üretimine izin vermektedir. Foggia ve arkadaşları küresel indüksiyon tip motor üzerinde hareket karakteristik oranının 60º konik olması için çalışmışlardır [5]. Oldukça geniş hareket oranına rağmen deneysel sonuçlar bu motorun prototipinin önemli bir gürültüsü olduğunu ve 5 saniyeden daha fazla tepkime zamanına sahip olduğunu göstermiştir. Küresel indüksiyon motorun farklı tasarımları karmaşık mekanik ve sargı tasarımına sahip olmasına rağmen Lee ve arkadaşları küresel adım motor kavramını keşfederek küresel yüzey üzerine 20 adet kutup yerleştirmişlerdir. Lee’nin küresel step motor tasarım kavramı ile yüksek kararlılıkta hareketler, gerçekte daha az sayıda rotor ve stator kutuplarıyla başarılabilmiştir [6]. Wang ve arkadaşları çok basit manyetik rotor montajı ve sargı düzenlemesiyle iki ve üç serbestlik derecesinde hareket edebilecek basit bir küresel motor geliştirmişler [7] ve motoru PD kontrolörle kontrol etmişlerdir [8]. John Hopkins Üniversitesi’nden Chirikjan ve arkadaşları rotorun hareketleri için kavramsal bir potansiyel enerji fonksiyonu geliştirerek rotorun denge durumunu sağlayacak manyetik alana bir yaklaşımda bulunmuşlardır [9]. Alman araştırmacı Kahlen ve arkadaşları son zamanlarda Lee tarafından önerilen benzer bir küresel adım kavramını rapor etmişlerdir [10]. Kullanılan rotora yerleştirilen sabit mıknatıslar ve statora monte edilmiş, tek tek kontrol edilebilen 96 adet elektromıknatıs ile 40 Nm’lik bir tork üretimi başarılmıştır. Kore’de Yang ve arkadaşları Lee’nin önerdiği benzer bir küresel eyleyici için PD kontrolör geliştirmişlerdir [11]. Hollis ve arkadaşları temel doğru akım prensiplerinde çeşitli serbestlik derecelerinde dönebilen çok hassas hareket tahrikli eyleyici geliştirmişlerdir [12]. Benzer bir DC küresel motor Kaneko ve arkadaşları tarafından yapılmıştır [13]. Bu motor sürekli dönebilmekte ve maksimum eğilmesi de ![]() ![]()
Küresel eyleyici teorisi, stator sargıları ve rotorun oluşturduğu manyetik kuvvetlerin bileşkesine göre küresel yapılı rotorun hareketi olarak tanımlanabilir. Bundan dolayı rotor ve stator sargılarının oluşturacakları kuvvetler küresel eyleyicinin çalışmasını etkilemektedir. Küresel eyleyicinin hareket serbestliğinin sağlanabilmesi ve enerji sarfiyatının az olması nedeniyle küresel rotor sabit mıknatıstan seçilmiştir. Böylece rotor hiçbir elektriksel bağlantı olmadan istenilen manyetik kuvveti sağlamış olacaktır. Stator sargılarına uygulanacak uygun bir akım ile oluşacak kuvvet rotorun oluşturacağı kuvvete karşı koyacak ve bu kuvvetlerin bileşkesi yönünde rotor hareket edecektir. Şekil 1’de sabit mıknatıslı küresel eyleyicinin yapısı şematik olarak gösterilmiştir. Tasarlanan küresel eyleyici stator, rotor, yataklar, sargılar ve stator kapağı olmak üzere beş parçadan oluşmaktadır [27]. Stator olarak küresel yapıda ve eyleyicinin dış kabuğunu oluşturan kalıcı mıknatıslığı az yumuşak demir malzeme kullanılmıştır. Yumuşak demir malzemenin B-H eğrisi Şekil 2’de gösterilmiştir. Statorun dış çapı 160 mm, iç çapı 140 mm, boyunduruk kalınlığı 10 mm ve stator kutupları küresel rotoru daha iyi kavrayabilecek şekilde 32 mm iç yarıçapta CNC tezgahında hassas bir şekilde işlenmiştir. ![]() Küresel eyleyicinin en önemli parçalarından biri olan rotor, eyleyicinin istenilen 3 serbestlik dereceli hareetini yapabilecek şekilde manyetik alan yönleri dikkate alınarak dört kutuplu olarak seçilmiştir. İdeal küresel rotor dört çeyrek küre diliminden oluşmaktadır. İdeal küresel rotorun yapısı ve manyetik alan yönleri Şekil 3’de gösterilmiştir. Bu tasarımın üç boyutlu manyetik analizleri yapılmış ve istenilen hareket yönlerinde tork oluşturduğu görülmüştür. Ancak tasarlanan yapıda bir küresel rotor mıknatıs imalatçıları tarafından temin edilemediğinden, rotor yüzeyinde aynı manyetik etkiyi oluşturabilecek farklı rotor tasarımları geliştirilmiştir. Bir, üç, beş parça silindirik sabit mıknatıs dört kutup oluşturacak şekilde küresel bir teflon yapının içerisine gömülmüş ve bu şekilde manyetik analizleri yapılmıştır. Yapılan manyetik analizler sonucunda beş parça silindirik sabit mıknatıstan oluşan rotor yüzeyindeki manyetik akı yoğunluğu ideal dilimli rotorun yüzey manyetik akı yoğunluğuna en yakın değerleri vermiştir. Bu nedenle uygulamada, beş parça silindirik sabit mıknatıstan oluşan küresel rotor kullanılmıştır. ![]() Kutupların mıknatıslanma vektörü yarıçap doğrultusundadır. Silindirik sabit mıknatıslar küresel teflon içine gömülerek küresel yapı kazandırılmıştır. Mıknatısların üzerindeki teflonun kalınlığı 1 mm’dir. Şekil 4’de görüldüğü gibi statorun demir nüvesinin açıklık kısmına karşılık gelen rotorun üst kısmında mekanik enerji alınabilecek bir mil çıkartılmıştır. Kullanılan sabit mıknatısın kalıcılık değeri Br=1,234 T ve manyetik geçirgenliği μr=1,15 dir [27]. Şekil 4’de beş parça silindirik sabit mıknatıstan tasarlanan küresel rotor yapısı ve Şekil 5’de ise rotor sabit mıknatısının B-H eğrisi görülmektedir. ![]() ![]() Küresel eyleyicinin yataklaması üzerinde önemle durulması gereken bir konudur. Eyleyici üç serbestlik derecede hareket yapacağı için henüz bu hareketi sağlayabilecek bir yatak tasarımı mevcut olmadığından çeşitli problemlerle karşılaşılabilir. Bu çalışmada yataklama sistemi olarak sürtünme katsayısı çok az olan teflon malzeme seçilmiştir. Rotorun boyutlarından 0,2 mm daha büyük bir teflon malzemeden iç küre işlenmiş ve rotor bu iç kürenin içine yerleştirilmiştir. Teflondan yapılan iç küre içerisine kayganlaştırıcı yağ konularak sürtünme azaltılmaya çalışılmıştır. Teflondan oluşan yataklama malzemesine sargıların ve stator kutuplarının gireceği yuvalar açılarak yatağın stator merkezinde sabitlenmesi sağlanmıştır. Tasarlanan sabit mıknatıslı küresel eyleyicide birbirinden bağımsız sekiz adet sargı kullanılmıştır. Sargı yapıları küresel statora daha iyi yerleştirilebilmesi amacıyla konik şekilde yapılarak stator yüzeyine tam bir dağılım sağlanmıştır. Sargılar 0,50 mm çapında bakır iletkenden 1000 sipir olarak sarılmıştır. ![]() Yumuşak demirden yapılmış stator kapağında rotor milinin çıkartıldığı bir boşluk bulunmaktadır. Stator kapağındaki bu açıklık, rotorun ±45º lik hareket yapabilmesi için tasarlanmıştır. Kapaktaki açıklığın sınırlarının değiştirilmesi rotorda oluşan döndürme torkunda da değişim yapmaktadır [25]. Şekil 6’da sabit mıknatıslı küresel eyleyiciyi oluşturan parçalar gösterilmiştir. 3. KÜRESEL EYLEYİCİDE MANYETİK ALAN DAĞILIMI ![]() Küresel eyleyicinin manyetik alan dağılımının hesaplanması eyleyicinin tasarım optimizasyonu ve küresel eyleyicinin doğru modelini kurmanın temelini oluşturur. Hava nüveli küresel eyleyicinin manyetik alan dağılımı küresel harmonikler içerdiğinden küresel koordinat sistemi kullanılarak analitik yolla hesaplanabilir. Şekil 7’de küresel koordinat sistemi görülmektedir. Manyetik alan iki alt bölgeye ayrılabilir; μ0 manyetik geçirgenliğine sahip hava aralığı veya sargı bölgesi diğeri ise μ0μr manyetik geçirgenliğine sahip sabit mıknatıs bölgesidir. Buna göre: ![]() Burada μr mıknatısın bağıl geçirgenliği ve M artık mıknatıslanmadan kaynaklanan mıknatıslanma vektörüdür. Doğrusal demanyetizasyon karakteristiğine sahip bir kalıcı mıknatısta μr sabittir ve mıknatıslanma vektörü, M ve kalıcı mıknatısiyet, Br arasında ![]() Mıknatıslanma vektörünün radyal bileşeni, Mr küresel harmonikler şeklinde açılarak aşağıdaki gibi ifade edilebilir [26]. ![]() Burada ![]() ![]() Burada ![]() ![]() eşitliği ile hesaplanabilir. Burada; ![]() ![]() ![]() Şekil 8’de stator nüvesi yokken rotorun Rm=0,030 m, r= 0,031 m ve θ=0˚ pozisyonu için α=-π ile α=+π aralığında (4) eşitliğine göre hesaplanmış manyetik akı yoğunluğunun radyal bileşenlerinin analitik, sonlu elemanlar ve uygulama sonuçlarının bir karşılaştırılması verilmiştir. Analitik metotla hesaplanan (ideal) akı yoğunluğu (r=0,031 mm yarıçapındaki bir çember üzerindeki normal bileşenleri, Bn) ile sonlu elemanlar metoduyla tahmin edilen akı yoğunluğu arasındaki fark Ansoft Maxwell 3D yazılımının istenilen bir çizgi üzerindeki akı yoğunluklarının normal (Bn) ve teğet (Bt) bileşenlerinin bileşkesini ![]() 4. MANYETİK SİMÜLASYONLAR VE UYGULAMA SONUÇLARI ![]() ![]() Sabit mıknatıslı küresel eyleyicinin kendisi oldukça basit bir geometriye sahip olmasına rağmen rotorun dönüşü süresince manyetik alanın belirlenmesi için üç boyutlu manyetik alan analizinin kullanılması gerekmiştir. Ansoft Maxwell 3D ile üç boyutlu manyetik simülasyon yapılmış ve buradan alınan sonuçlar ışığında eyleyicinin prototipi gerçekleştirilmiştir. Hava nüveli bir küresel eyleyici için rotorda oluşan torkun analitik yolla hesaplanması mümkündür (Ek.A). Fakat uygulaması yapılan demir nüveli küresel eyleyicide stator sargıları enerjili iken rotorun açısal yer değişimleri için relüktans torku ve kuvvetlerin (problem simetrilikten uzak 3 boyutlu bir özelliğe sahip olduğundan) analitik yolla hesaplanması zorlaşmakta ve bu Ansoft Maxwell 3D tarafından otomatik olarak Zahiri İş Metodu ile hesaplanmaktadır. Küresel motorun alan çözümünde sınır şartı (Statorun iki katı yarıçapa sahip en dış küresel bölge yüzeyinde Hn=0 dır. Hava nüveli küresel eyleyicinin sonlu elemanlar analizinde, doğru sonuçların elde edilmesi için bu yarıçapın 10-100 kat arası bir değerde alınması gerekir) ve bir dönüşe ait sargı akım değerleri Şekil 10’daki çizelgede görülmektedir. Simülasyon sonuçlarına göre rotorun farklı pozisyonlarına göre farklı torklar oluşmaktadır. Rotor pozisyonu α ve θ açısal değişkenlerine göre küresel koordinat sisteminde α = 0-2π ve θ = 0-π/2 aralığında değişmektedir. Rotorun hareket ettirilmesi için örneğin A-A' sargılarına 1000 amper-sarımlık akım verildiğinde rotorun S kutup ekseni ile A sargısının manyetik eksenleri aynı doğrultuya gelmek ister ve A' sargısı da rotorun N kutbunu iterek bu harekete yardım eder. Böylece akımların miktarına bağlı olarak rotor yönünde 0˚- 45° dönme gerçekleştirebilir. Bu dönmenin akabinde A-A' sargılarının akımı kesilerek B-B' sargılarına aynı miktarda akımlar verilerek daha ileri bir 45°’lik dönme gerçekleşir. Bu işlem böylece devam ettirilerek istenildiği kadar dönme hareketi yaptırılabilir. Rotorun eksenindeki (aşağı-yukarı) hareketleri için de uygun sargı kombinasyonunun enerjilendirilmesi gerekir. Tüm bu aralıklar içinde kalan rotor pozisyonları için 5° ’lik aralıklarla toplam 648 adet üç boyutlu statik manyetik simülasyon yapılmıştır [27]. Küresel eyleyici üç boyutlu geometriye sahip olduğundan yapılan 3 boyutlu simülasyonlar oldukça büyük miktarda bilgisayar hesaplama zamanı gerektirmektedir. Örneğin sadece bir pozisyona ait manyetik akı, indüktans ve tork değerlerinin bulunması için gereken zaman P4 1,7 Ghz bilgisayar için 360 dakikadır. ![]() Sabit mıknatıslı küresel eyleyicinin θ ve α değerlerine göre 1 amper sargı akımında ürettiği tork değerleri Şekil 12’de üç boyutlu grafikte verilmiştir. Grafikte θ ve α konum değişkenlerine göre torkun farklı değerler alması rotordaki sabit mıknatısların yönlerinden kaynaklanmakta ve bundan dolayı kontrol sisteminde bu değişimin dikkate alınması önemlidir. ![]() Sabit mıknatıslı küresel eyleyicinin y ve x eksenlerinde α=0˚ iken θ =-π/4 ile θ =+π/4 aralığında ve 1A akım uygulandığı durumda simülasyon sonucu elde edilen tork eğrileri ve bu konumlara ait tork algılayıcısı (Şekil 13) ile ölçülen tork değerlerinin karşılaştırılması Şekil 14’de görülmektedir. Simülasyon ve uygulama sonucu elde edilen tork grafikleri kabul edilebilir bir yakınlıkta örtüşmektedir. 5. KÜRESEL EYLEYİCİNİN AÇIK ÇEVRİM KONTROLÜ ![]() Sabit mıknatıslı küresel eyleyicinin sekiz adet sator sargısı bulunmaktadır. Stator sargılarının her biri ayrı ayrı enerjilendirilebilmektedir. Bu yönüyle küresel eyleyici adım motoru prensibine yakın bir prensipte çalıştırılabilir. Sargılar istenilen sırada yeteri kadar bir akımla uyarıldığında küresel eyleyicinin hareketi sağlanmış olur [27]. Şekil 15a’da sabit mıknatıslı küresel eyleyici ve kontrol düzeneğinin fotoğrafı görülmektedir. Kontrol elemanı olarak PIC 16F877 mikro denetleyicisi kullanılmıştır. PIC 16F877’e yazılan programla, seçilecek değerlere göre eyleyicinin hareketi sağlanmış olur. Dört adet buton kullanılmıştır. Butonlardan birincisi tab butonu, ikincisi değerleri artırma, üçüncüsü değerleri azaltma ve dördüncü butonda start butonudur. Program çalıştırıldığında LCD display’de kutup seçimi, hareket yönü seçimi, adım seçimi ve hız seçimi gibi seçenekler vardır. Bu seçenekleri uygun konumlara getirerek start butonuyla PIC 16F877’nin çıkış vermesi sağlanır. Bu sinyaller 4N25 ile yapılmış yalıtım devresine gelir. Yalıtım devresine gelen sinyal sargı akım yönlerinin değiştirildiği H tipi köprü sürücü devresine gelir. Şekil 15b’de H tipi köprü sürücü devre gösterilmiştir. H tipi köprü sürücü devre çıkışı eyleyici sargı uçlarına bağlanmıştır. Böylece PIC 16F877’de üretilen bir sinyal küresel eyleyiciye ulaştırılmış olacaktır [27]. 6. SONUÇLAR Bu çalışmada üç serbestlik dereceli sabit mıknatıs rotorlu demir nüveli bir küresel eyleyicinin statik manyetik analizi ve uygulaması yapılmıştır. Uygulaması yapılan küresel eyleyici açık çevrimli olarak 360° dönme hareketinin yanında ![]()
KAYNAKLAR ![]()
Ek.A: Hava nüveli bir küresel eyleyicide rotorda oluşan torkun Zahiri İş Metodu ile analitik yolla hesaplanması [24]: Stator sargılarından geçen akımın manyetik alanı ile rotor mıknatıslarının manyetik alanının etkileşimi sonucu rotor üzerinde aşağıdaki eşitlik ile verilen bir tork oluşur. ![]() Burada J, Şekil E.1’de görüldüğü gibi v hacmindeki sargının akım yoğunluk vektörünü temsil eder. Şekilde görüldüğü gibi r R0 ile Rs arasında 0 ile 1 arasında değişim göstermektedir. Yüzeyi ![]() ![]() Blr ’nin temel bileşenleri (l=2, m=1) baskın olduğundan yüksek dereceli harmonikler ihmal edilerek tork aşağıdaki gibi verilir. ![]() Burada ![]() ![]() ![]() ![]() ![]() (E.4) ile bulunur. Eşitlikten görüleceği gibi Tm torku Brem ve J ’ye bağımlı olduğu gibi rotor ve sargıların geometrik parametreleri ile de ilgilidir. Rotor üzerindeki toplam tork bütün sargı akımlarının bileşke etkisi ile elde edilir. |
![]() | ![]() |