![]() |
|
|
Sabit Kaynaklı RC Devresi Analizi / 20 Şubat 2008
Düzenlemeye çalışıyorum(yazının eski versiyonu için tıklayınız) Merhaba arkadaşlar bu yazıda RC devrelerinin analizini yapmaya çalışacağız, genellikle kitaplara formül olarak düşülen denklemlerin aslında basit bir birinci mertebeden diferansiyel denklemin çözümünden türediğini göreceğiz. RC devrelerinde adından da anlaşıldığı üzere ve yukarda da resmi görüldüğü üzere seri bağlı direnç ve kapasiteden oluşan bir devrenin davranışı incelenir. Biz bu basit devrenin analizinde devre analizinin en temel kanunlarından olan KVL yani kirchoff voltaj yasasını kullanacağız. I yönünde bir KVL dönersek; 1) Kapasitenin tanım denkleminin Q= C. V olduğunu biliyoruz. Her iki tarafın t ye göre türevini aldığımızda: dq/dt = C . dV(t)/dt olur. Yani t anındaki kapasite akımının voltajı cinsinden değeri; I(t) = C . dV(t)/dt olur. Buna göre KVL den yazdığımız ilk eşitlik 2) – ε + C. dv(t)/dt . R + V(t) = 0 haline gelir. kolaylık olması açısından dV(t)/dt ifadesi yerine V’ yazarsak eşitliğimiz: 3) C.R.V’ + V = ε haline gelir ki bu denklem görmeye alışık olduğumuz diferansiyel denklem modellerindendir; birinci mertebeden, doğrusal, homojen olmayan dif. denklem(first order linear non-homogen diff. equ.). Bu denklemi birinci mertebeden doğrusal diferansiyel denklem(first order linear dif. equ) olarak düşünüp çözebileceğimiz gibi(integral faktörü bularak) yüksek mertebeli diferansiyel denklem çözüm yöntemlerinden birini kullanarak da çözebiliriz. Tabi burada bir mesele de ε un yani devremizin kaynağının nasıl bir davranış gösterdiğidir. Kaynağımızın zamana bağlı değişim gösteriyor olup olmaması çözümümüzü değiştirecektir. Bu yazıda sabit kaynaklı devreleri değerlendirdiğimiz için kaynağımızın sabit olduğunu düşünerek 3 nolu denklemimizi çözelim. Diferansiyel denklemimizin çözümünün Y=Yh + Yp şeklinde olduğunu biliyoruz. Burda Yh , V(t) nin geçici(transient response) tepkisi, Yp V(t) nin kalıcı(steady-state response) tepkisini göstermektedir.Çözümü bulduğumuzda t sonsuza giderken geçici etkinin ortadan kalktığını, kalıcı tepkinin devam ettiğini göreceğiz. Önce Yh i bulalım: C.R.V’ + V = 0 Geçici hal çözümümüzün aşağıdaki formda olduğunu biliyoruz, bulmamız gerekenler c1 ve r değeri. V’ = c1 . r . e^(r*t) bu iki ifadeyi yukarda yerine yazdığımızda: c1 . e^(r *t) [ C. R .r + 1 )]= 0 olur, bu durumda eşitliğin sağlanması için C.R.r = -1 olmalıdır. r = – 1/ R.C bu durumda Yh (homojen kısmın çözümü): Yh = c1. e^(-t/RC) olur.(c1 sabit bir sayıyı temsil ediyor, kapasite sığası C ile bir ilgisi yok, c1 değeri kapasitenin ilk voltajıyla ilgilidir) Sıra Yp: ın çözümüne geldi, yöntem olarak bilinmeyen katsayılar yöntemi (method of undetermined coefficients) kullanabiliriz, zaten kaynağımız da sabit olduğu için işimiz oldukça kolay. Sağ tarafta 0. dereceden bir polinom gördüğüm için yöntem gereği önereceğim Yp= A olacaktır. Yerine yazıp A yı çektiğimizde(Yp‘ 0 olur ve Yp direkt ε a eşit olur) Yp: V(t) = c1. e^(-t/RC) + ε olur. Burda c1 sabitini bilmiyoruz, bulabilmemiz için bir başlangıç değerine sahip olmamız gerekir( yani herhangi bir t anındaki V(t) yi biliyor olmalıyız ). V(0) = V0 olarak bildiğimizi varsayalım, genelde 0 anındaki değer verilir veya bulunur. V(0) = c1 + ε = V0 olacağından c1 = V0 – ε olur. c1 i yerine yazdığımızda: V(t) = (V0 – ε)* e^(-t/RC) + ε olduğunu buluruz. Böylece bu devrede kapasitemizin geriliminin zamana bağlı davranışını ifade etmiş olduk. Şimdi de bu sonucu biraz yorumlamaya çalışalım: t=0 anında kapasite gerilimimiz V(0), eşitliğimizde t yerine 0 yazarsak: V(0) = (V0 – ε). e^(-0/RC) + ε = V0 – ε + ε = V0 bu sonuç bizi şaşırtmadı çünkü biz hesaplama yaparken kapasitenin 0 anında bir V0 gerilimine sahip olduğunu varsaymıştık. V0 gerilimimiz 0 da olabilir, yani kapasitemizin başlangıç gerilimi olmak zorunda değil. V(t) nin gidişatına baktığımızda eksponensiyel azalmanın olduğu (V0 – ε) ifadesinin giderek etkisinin azalacağını ve t = sonsuz olduğunda sadece ε değerinin kalacağını görürüz. Nitekim t yerine 0 yazdığımız gibi sonsuz yazarsak V(sonsuz)=ε olduğunu görürüz. Anlaşıldığı üzere ekponensiyel kısım sadece bir süre etkisini gösterdikten sonra etkisini kaybedecek ve kapasite gerilimi -kaynak varlığını sürdürdüğü müddetçe- ε olarak kalacaktır. Ne kadar bir süre etkisini gösteriyor? Vampire in Brooklyn movie download Yukarda da söylediğimiz gibi ekponensiyel kısım ancak sonsuzda 0 olur ve etkisinin tamamen yitmesi için sonsuza kadar beklememiz gerekir . “Zaman sabitinin birimi nasıl saniye olabilir? ” sorusu akıllara gelmiş olabilir, şöyle: kapasite tanım bağıntımız Q=C.V idi, bu durumda C nin birimi coulomb / V dir. RC nin birimi ise ohm*coulomb/V dir. Ohm/V nin 1/A (bir bölü amper) olduğunu ohm yasasından biliyoruz. Bu durumda RC nin birimi coulomb/A oldu. Akımın tanım bağıntısın dQ/dt olduğunu biliyoruz, yani akımın birimi aynı zamanda coulomb/t dir. Bu durumda RC nin birimi coulomb/(coulomb/t) olacağından RC miz saniye cinsinden bir değer olacaktır(akım bağıntısında dQ/dt deki Q coulomb t de saniye cinsindendir). Hatırlayacağınız üzere V(t) nin ekponensiyel azalmanın olduğu kısmını ilk yazdığımız diferansiyel denklemin(3 nolu denklem) homojen çözümden; ε u bulduğumuz kısmı da diferansiyel denklemimizin particular(kısmi) çözümünden bulmuştuk. Homojen ve particular ifadeleri matematiksel ifadelerdir, bu ifadelerin devremizdeki karşılıkları ise sırasıyla, geçici ve sürekli hal tepkisi dir. Burda geçici hal: (V0 – ε). e-t/RC ifadesinden gelmekte, kalıcı hal ise ε ifadesinden gelmektedir, umarım neden geçici ve kalıcı dediğimiz anlaşılmıştır. Şimdi isterseniz sonucumuzu görselleştirebilmek için kapasitemizin gerilimin zamana göre değişimini veren ifadenin yani V(t) = (V0 – ε). e-t/RC + ε matlab yardımıyla grafiğini çizelim. Grafiğimizi çizdirmeden önce nasıl bir grafikle karşılaşırız sorusuna cevap aramakta fayda var.İlk durumda kapasitemizdeki gerilimin 0 olduğunu yani olmadığını düşünelim, her şeyden önce ilk değerimizin 0 ve son değerimizin ε olacağını denklemimize baktığımızda görebiliyoruz ve 5RC zamana kadar bir geçiş(artış) olacak(5RC zamandan sonra çok az bir artış sonsuza kadar devam edecek ancak biz o artışı ihmal ediyoruz) ve sonrasında kapasite gerilim ε olarak kalacaktır. Fiziksel yapısı gereği kapasite gerilimi hiç bir zaman sıçrama yapamaz yani ani bir şekilde 3V gibi bir gerilimden 5V gibi bir gerilime atlayamaz. Kapasite geriliminin ani değişememesi sebebiyledir ki kapasite elemanının durum değişkeni(state variable) voltajıdır, durum değişkeni herhangi bir sistemin durumuyla ilgili bize öz bir bilgi verir. Mesela bir kapasite elemanı üzerinde 4V luk gerlim varsa deriz ki demek ki bir kaynak bu kapasiteye bir süre uygulanmış. Kapasite belleği-hafızası- olan bir elemandır ve belleğinde voltaj değerini tutar. Eğer kapasitenin durum değişkenine akımıdır demeye kalkacak olursak aşağıda göreceğimiz gibi bir anda 0 A den 5 A gibi bir değere çıkabilen veya 5A den bağlantıyı açtığımızda bir anda 0A a düşebilen bir durum değişkeni bize ne ifade edebilir ki! V0= 0, kaynağımız ε=10V, R=10k ohm C=0.1 uF (RC zaman sabitimiz= 10000ohm * 0.1 uF = 1 ms oldu) Beraber çizdirelim, Matlab’ı açalım: t=0:0.0001:10*0.001; %yazarak 0 dan 10 milisaniyeye kadar 100 değer oluşturalım. %(5. ms de kapasitemizin hemen hemen dolacağını biliyoruz ) %Formülümüzü yazalım: V=(0-10)*exp(-t/10^-3)+10; %ve çizdirelim: plot(t,V); %Eksenlerimize isimlerini verelim: xlabel('zaman(saniye)'); ylabel('Gerilim(V)'); grid Dediğimizde grafiğimiz: Neticeten kapasite olarak adlandırdığımız devre elemanımız nasıl davranıyor? Görüldüğü gibi kare dalgamız max değerine aniden çıkmasına rağmen kapasite gerilimimiz hemen yetişemiyor(5RC zaman gerekiyor) ve kare dalgamız 0 a düştüğünde de kapasite gerilimimiz hemen sıfıra düşemiyor. Peki kare dalgamızın periyodunu biraz arttırsaydık ne olurdu? Kapasite dolmak için yeterlli zamanı(5*R*C) bulacağından rahat rahat dolabilecekti. Seçtiğimiz örnek değerlere göre biliyoruz ki 5*R*C çarpımı 5ms, o zaman periyodu 20 ms seçelim bakalım ne oluyor: Kapasitemizin dolabilmesi(veya boşalabilmesi) için hiç zaman tanımayan 300us periyotlu kare dalga uygularsak: Yukarıdaki denemelerimizden çıkaracağımız önemli bir sonuç da şudur ki: Kapasite elemanımız üzerine uygulanan AC gerilimin frekansı düştükçe(periyodu arttıkça) kapasitemiz bu gerilimi takip edebildiği için daha yüksek tepkinlik gösterir, yani varlığını daha iyi hissettirir. Frekans arttıkça kapasite elemanı varlığını hissettirememeye başlar, yani düz bir tel gibi davranmaya başlar. Yukarıda ikinci yaptığımız 300us periyodu biraz daha düşürecek olsaydık(frekansı arttırsaydık) artık kapasite uçlarında neredeyse hiç gerilim okuyamayacaktık. “Klişe” bir cümleyle kapasite elemanı yüksek frekanslarda kısa devre gibi, düşük frekanslarda açık devre gibi davranır. Bu yazıda yaptığımız analizlerin bize kalan bir diğer yegane pratik bilgisi de şudur: Bir kapasite ve bir direnç seri bağlı ve biz bu devreye A V büyüklüğünde DC bir gerilim uyguladık, ne kadar süre sonra kapasite gerilimimiz A V olur, cevap basit: 5*R*C kadar süre sonra. Peki, B V a dolmuş bir kapasitenin kaynak bağlantısını kesip bir R direnci üzerinden boşaltmak istiyoruz, ne kadar sürede boşalır, 0 V olur? 5*R*C sürede. Kapasitemizin geriliminin zamana göre değişimini bir üstte çizdirdik, peki akımı? Yukarda ilk anda max olduğunu sonra 0 a kadar azaldığını söyledik. Şimdi de yine Matlab yardımıyla (yukarda verdiğimiz değerlerle) kapasitemizin akımının zaman göre değişimini çizdirelim. Değerlerimiz şöyleydi: Bunun için önce akımın zaman bağlı değişimini veren ifadeyi bulalım. Kapasitemizin akımının: I(t) = C . dV(t)/dt olduğunu biliyoruz, o zaman önceden yazmış olduğumuz V(t) ifadesinin t ye göre bir defa türevini alıp C ile çarptığımızda i(t) ifadesini bulabiliriz. V(t) miz i(t)=C*-10*exp(-t/10^-3)/ (-1/10^-3) olur , yani: i=0.1*10^-6*-10*exp(-t/10^-3)*(-1/10^-3) ifadesini matlaba verip t ye bağlı değişimini çizdirmek istediğimizde; plot(t,i) grid dediğimizde grafiğimiz: Akımımız 10 mA den başlayıp(0 anında 0 amperden 10 mA e ani bir çıkış var, demek ki kapasite akımı ani değişim gösterebiliyormuş.) 5RC saniye sonra -kapasite gerilimimiz 10V a ulaştığı için- 0 a düşüyor. Sorular: “Sabit Kaynaklı RC Devresi Analizi” için 5 YorumYorum Yapın |
|
| |
25 Ağustos 2008, 09:46
şimdi baktım.bravo helal olsun.yubil de güzel.tebrikler böyle çalışan insanlara her daim ihtiyacımız var Allah kolylık versin
25 Ağustos 2008, 09:46
tebrik ederim çok iyi bir şekilde hazırlamışsınız yazıyı. sayenizde bir şeyler öğrendim bu konu hakkında çünkü derste hocalarımız ya tam açıklamıyorlardı ya da hızlı geçiyorlardı ama siz de bizim gibi bir öğrenci olduğunuz için aynı sorular oluşmuş olsaki kafamızda tam istediğim gibi anlatmışsınız. TEŞEKKÜRLER…
21 Eylül 2008, 15:01
yazılarınız çok başarılı siteniz tam bir kaynak.teriyle pratiğin birleştiiği bilgiler bulunuyor. elinize sağlık gerçekten.başarılarınızın devamını dileriz..
18 Nisan 2009, 01:32
çok teşekkür ederim ,sayenizde derste anlamadığım bu konuyu matlabda deneyerek öğrendim, allah razı olsun.
09 Mayıs 2009, 18:45
Müthiş bir çözümleme. Bu tür kaynak ancak konunun çeşitli boyutlarına vakıf bilim gönüllüleri tarafından hazırlanabilir. Bu tür kaynaklar Türkiye’nin geleceğini oluşturacak gençlere çok faydalı olacaktır. Eline sağlık…