Termal yönetim elektrik mühendisliğinde kritik bir endişedir. Her elektronik bileşen - basit bir dirençten karmaşık bir mikroişlemciye kadar - operasyon sırasında ısıyı optimize etmek, mühendisler risk erken bileşen başarısızlığı, performansı azaltmak veya felaket sistem arızası sağlar.Diferansiyel denklemler, zaman içinde sıcaklıkların nasıl geliştiğini tahmin etmek için matematiksel temel sağlar, tasarımcılara uygun malzemeleri seçmelerini sağlar ve uzun vadeli güvenilirlik sağlar.
Bu makale, diferansiyel denklemleri kullanarak termal modellemenin temel özelliklerini araştırıyor, temel ısı denklemlerinden mühendislik tasarımında pratik uygulamalar inşa edebilirsiniz. Elektronik bileşenler için termal modelleri nasıl oluşturacağınızı öğrenecek, yönetim denklemlerini yorumlayın ve bu teknikleri gerçek dünya problemlerine uygularsınız.Sonunda, matematiksel denklemlerin termal analizlerini nasıl pratik uygulamalarla ve neden modern elektronik elektronik elektronik sistemlerde vazgeçilmez olduklarını sağlam bir kavrayacaksınız.
Neden Termal Modelleme Maddeleri
Elektronik bileşenler elektrik direnci, geçiş kayıpları veya diğer inefficiler nedeniyle ısı olarak gücü dağıtır.Sıcak ısı kaldırıldığında, bileşen sıcaklığı yükselir. Elevated sıcaklıklar, ısı geçişi gibi başarısızlık mekanizmalarını hızlandırır ve uygun ısı parçaları, fanlar veya sıvı soğutma sistemleri, fiziksel prototipler olmadan 10°C artış gösterir.
Diferansiyel denklemler, sıcaklık zaman ile nasıl değiştiğini tanımlamak için doğal dildir. ısı nesli, depolama ve transfer arasındaki dengeyi yakalarlar, analitik veya sayısal olarak çözülebilir. Bu modelleme yaklaşımı endüstriler arasında kullanılır - otomotiv güç elektronik ve havacılık sürücülerine.
Heat Transfer Mechanisms
Heat üç birincil mekanizmayla hareket eder: iletim, konveksiyon ve radyasyon. Elektronik sistemlerde, sağlam malzemeler (örneğin silikon, bakır ve alüminyum) ve sıvı soğutucuları çevreleyen sıvıları (hava veya sıvı soğutucu) hakimdir. Radyasyon genellikle çok yüksek sıcaklıklar veya vakum ortamlardan başka bir şekilde ihmal edilebilir, ancak tamlık için dahil edilebilir.
- [FONT:0]Condüksiyon[[Dönetici:)[Dönetici: 0): 4B][/ ⁇ =) + [Düşük sıcaklık yüksek sıcaklık bölgelerinden sıcaklık gradient ve termal iletkenlik. üç boyutta, ısı iletkenliği: ısı iletkenliği = ⁇ T/ ⁇ t = ⁇ (k ⁇ T) + Q[D)
- [FONT=0)[Dönetici[Dönemli)[Dönemli))[Köpektif ısı transferi Newton’un serinleme yasası tarafından tanımlanır: [[QUYGÜN:2).
- [FONT=0)Radiation[[Döncükler: Katı olmayan radyasyon, Stefan-Boltzmann kanununu takip eder: 03:2).q = ⁇ A (T^4 - T surroundings^4)[DDDDDDönersizce)[Dönersiz)[Dönersiz)[Dönersiz)[Dönersiz)[Dönergesel radyasyonu genellikle küçük sıcaklık farklılıkları için lineerleştirilmiştir.
Birçok elektronik bileşen için, topak parametre modeli - bu varsayımın tek bir düğüm olarak ele alındığı gibi - yeterli doğruluk. Bu model, L c / k) ile normal bir diferansiyel denklemin sayısal bir denklemi, iç sıcaklık gradients tarafından kabul edilebilir olduğu varsayılır.
Temel Diferansiyel Denklemi Ayarlayın
Yüzeyin en güçlü dirençleri gibi küçük bir elektronik bileşeni düşünün, baskılı bir devre kuruluna (PCB) Kombine monte edilen ısı miktarındaki ısıyı C veya K olarak üretir.
Enerjinin korunması: depolamak için ısı enerjisinin değişimi ısı jenerasyonu eksi ısı kaybına eşittir:
[0]) c dT/dt = P - h A (T - T ambient)).
Bu ilk sipariş lineer sıradan diferansiyel denklemi yeniden yazılabilir:
[0])))
İlk koşul ile T **** = T 0 (genellikle ortam sıcaklığı), çözüm:
[0]T(t) = T ambient + (P / (h A) (1 - exp (t / ⁇ ) + (T 0 - T ambient)[Dış / ⁇ )[Dönem)[Dönemli)[Dönt / ⁇ )[Dönemli)
[FONT=0) ⁇ = m c / (h A) ) termal zaman sabitdir. Bu üstel yaklaşım, sabit devlet değeri T 0'dan 6,0'a doğru yükselir; T ambient + P/(hA) Sürekli sabit bir şekilde ⁇ .
Mühendisler bu modeli zirve ısıtımı tahmin etmek için kullanırlar, termal arayüz malzemeleri seçin ve pasif soğutmanın yeterli veya aktif soğutma olup olmadığını değerlendirin (fans, ısı boruları) gereklidir.Toplantılı modelin basitliği erken tasarım hesaplamaları için idealdir, ancak daha ayrıntılı analiz, iç sıcaklık gradyanları veya birden fazla ısı kaynakları ile bileşenler için gereklidir.
Basit Lumped Modelinin Ötesinde: Termal Direniş ve Kapacitance Networks
Gerçek elektronik bileşenler - özellikle entegre devreler (ICs) ve güç modülleri - birden fazla malzeme katmanı olan karmaşık iç yapılara sahiptir. Tek bir topakt düğümü, silikonun ölmesi, paketleme ve termo kapasiteleri arasındaki sıcaklık farklarını yakalayamaz.
Basit iki düğüm ağı için diferansiyel denklem (di ve vaka):
[0]C j dT j/dt + (T j - T c)/R jc = P).
[0]C c dT c/dt + (T c - T ambient)/R ca - (T j - T c)/R jc = 0).
Burada, C j ve C c, ilk kez ODE'lerin ısıtılması ve sayısal olarak analiz edilebilir. Elde edilen sıcaklık profilleri, başlangıçta sıcaklıktan daha hızlı yükseleceğini gösteriyor ve P, iki katın sürekli olarak belirlendiği güç değerleridir.
ICs ve güç transistörleri için standart veri tabloları genellikle R jc ve R ca için değerler sağlar (veya R j-a ek termal dirençler ve kapasitanslar için eklenir). ancak bu değerler belirli koşullarda ölçülmektedir (örneğin, tanımlanmış bir PCB ayak izi, hava akışı veya ısıtılır).
Kompleks Geometri için Numerical Method
Metal veya çok fazla-node yaklaşımı yetersiz olduğunda - örneğin, çok fazla çip modüllerini yakından uzaylı ısı kaynakları veya bir PCB ile anisotropic termal iletkenliği ile analiz ederken -mühendisler sonlu fark yöntemi (FDM) veya sonlu elemanlar yöntemi (FEM) gibi sayısal tekniklere dönüşürler.
Uzaysal alan ağına ve ısı denkleminde kısmi türevlere fark edilen bir fark yöntemi, x-kücut boyunca bir boyutlu problem için, açık bir ileri-zaman merkezi-uzay planı:
[0]T i^{n + 0 = T i^n + αt / (345x) ^2 (T {i + 0 {i i^n + T {i-1} ^n) + (345t / ( ⁇ c) Q i^n).
α = k/( ⁇ c) termal diffusivity, 345x ve 345t, uzaysal ve zaman adımları ve süperscript n zaman seviyesini gösterir. Stability, Courant-Friedrichs-Lewy (CFL) koşulunun karşılaştırılması gerekir: α {{t / (345x) ≤ 0,5x)
FEM, ANSYS gibi yazılımlarda uygulanmış, COMSOL ve OpenFOAM, bu simülasyonların arkasındaki temel hesaplama motorudur. Bu araçlar otomatik olarak ağlar oluşturur, malzeme özellikleri üretir ve ısı denklemlerini güçlüken, önemli hesaplama kaynakları ve uzmanlık gerektirir.
Birçok elektrik mühendisi ayrıca SPICE benzeri devre simülatörüni termal ağ modeline de kullanır. Sıcaklık direnci ve kapasitanslar eşdeğer elektrik bileşenleri olarak uygulanır ve simülatörün ortaya çıkan diferansiyel denklemleri çözer. Bu yaklaşım aynı anda ısı ve elektrik davranışını analiz eder - kısa devre dışı olayları sırasında.
Pratik Örnek: Bir Güç MOSFET'yi Bir Pulse sırasında modelleme
Bir geçiş dönüştürücüsinde bir güç MOSFET düşünün. Transistör kısa aralıklar için mevcut davranır, patlamalarda önemli ısıyı üretir.Kampiyon sıcaklığı üreticinin maksimum puanının altında kalmalıdır (tipik olarak 150°C 175°C'ye kadar) hasardan kaçınmak için.
MOSFET için basitleştirilmiş bir termal model üç düğümle inşa edilebilir: kavim, dava (veya sürmüş), ve ek bir ısısink.Diferansiyel denklemler oluşur:
[0]C j dT j/dt = P (t) - (T j - T c)/R jc).
[0]C c dT c/dt = (T j - T c) / R jc - (T c - T h) / R cs).
[0]C h dT h/dt = (T c - T h) / R cs - (T h - T ambient)/R ha[Dönem: 1)
R cs'in durum ısıtıcı ısıtıcı direnci (parat arayüz malzemesi dahil), R ha ısıtıcı direnç ve P (t) güç nabız dalgaformu (örneğin, bir dizi dikdörtgen puls).
Bu denklemleri çözün - Laplace dönüşümlerini analitik olarak bile kullanarak - hat sıcaklığının pul sırasında hızla yükseldiği ve puls arasındaki çürümeler sırasında hızla arttığını gösteriyor.Toplu sıcaklık sadece ortalama güçte değil aynı zamanda nabız genişliği ve görev döngüsü üzerinde de önemlidir.Bu tür bir analiz güvenli çalışma alanları (SOA) tasarlamak ve uygun ısısink termal kütle seçmek için çok önemlidir.
Bu eğrileri, süperpozisyon yöntemleri ile (örneğin, Cauer ağ gösterimleri) arasındaki hesaplı ısı sıcaklık artışını tahmin etmek için bu eğrileri kullanır.
Nonlinear Etkileri ve Değişken Özellikleri
Birçok gerçek sistemde, malzeme özellikleri sürekli değildir. Silikon ısı iletkenliği sıcaklıkla azalırken, belirli ısı geçişi katsayısı doğal konveksiyon için ısı transfer katsayısı sıcaklık farkı ve geometriye bağlıdır (Nsselt numara korelasyonu ile). Radyasyon, T4 terimini nadir hale getirir ve sayısal yöntemler gerekli hale gelir.
Örneğin, Newton'un soğutma yasasında katançlar şöyle ifade edilebilir:
[FONT=0)= C (345T / L) ^n ).
{{T = T surface - T ambient, L karakteristik bir uzunluktır ve C ve n, dikey plaka üzerinde laminar doğal konveksiyon için ⁇ 0.25. bunu diferansiyel denklemlere ekleyin:
[DÜDÜ:0) c dT/dt = P - A C (T - T ambient) / L ) ^n (T-T ambient)).
Bu, ilk sıra dışı bir ODE. Runge-Kutta yöntemlerini sayısal olarak kullanarak çözülebilir. Lineer olmayanlığı görmezden gelmek, özellikle de sıcaklık artışının büyük olduğu durumlarda, yüksek güç veya aşırı-environment uygulamaları için mühendisler hesabın gerektirdiğini varsayabilir.
Isif Modelleme ve Geçerlilik
Matematiksel model doğrulanmış ve doğrulanmış değildir. Doğrulama kontrolleri denklemlerin doğru şekilde çözüldüğünü kontrol eder (örneğin, programlama böcekleri), modelin deneysel verileri oynadığını kontrol eder.
- Bilinen bir güç döngüsü sırasında bir termo çift, kızılötesi kamera veya entegre sıcaklık sensörü ile durum veya kavim sıcaklığı ölçmek.
- Diferansiyel denklem veya ağ modelinin çözümüne yönelik olarak ölçülen geçici cevabı karşılaştırmak.
- Verilere uygun olarak fiziksel olarak makul olmayan aralıklarda (parametre tahminleri) belirsiz parametreleri (örneğin, h veya R jc) ayarlamayı sağlayın.
- Hangi parametrelerin sıcaklık tahmini üzerindeki etkisini belirlemek için hassasiyet analizi yapmak.
Mühendisler genellikle veri tablolarından ilk geçerlilik için termal dürtü eğrilerini kullanırlar, ancak bu eğriler gerçek uygulamadan farklı olabilecek standart koşullar altında ölçülmektedir. Özel bir doğrulama modeli inşa - sadece genel sayılara güvenmekten çok daha güvenilir.
Diferansiyel Equation-Based Termal Simülasyon için Yazılım Araçları
Birkaç yazılım paketi, termal diferansiyel denklemlerin elektronik bileşenler için oluşturulması ve çözümü basitleştirir:
- [FONT=0)SPICE tabanlı simülatörler[[Dönetici: LTspice, PSpice ve diğer SPICE varyantları, direnç ve kapasite kullanan termal R-C ağlarının doğrudan uygulanmasına izin verir.
- [FONTNT S:0]MATLAB/Simulink[DÜT:1): Mühendisler ODE'leri açıkça yazabilir ve bunları inşaatlı çözücülerle (Dörtüncü45, ode15s, vb.) Simulink, Simscape ürünü ailesindeki bloklar oluşturur.
- [FONT:0]Python ile SciPy[DÜT:1): "Ninve ivp' işlevi, sert ODE'lerin büyük sistemlerini ele alabilir. “CoolProp’un malzeme özelliklerini sağladığı gibi kütüphaneler ve "pytermal" modülleri özel elektronik soğutma analizi için mevcuttur.
- [FONT:0]CFD araçları : ayrıntılı 3D termal modelleme, ANSYS Buzpak, Flotherm ve OpenFOAM tam yönetim denklemlerini çözüyor (Navier-Stokes ve enerji). Bu araçlar hava akış modelleri ve yerel sıcak noktaların çözülmesi gereken temeldir, ancak daha fazla zaman ve hesaplama kaynakları gerektirir.
Her araç güçlü ve zayıf yönleri vardır. Seçim gerekli doğruluk, mevcut bütçeye bağlıdır ve mühendisin platformu ile aşinalığı. Ön tasarım için, analitik bir topak model veya SPICE'deki basit bir R-C ağı genellikle anahtar kararlar almak için yeterli bir fikir sunar.
Mühendislik Tasarımlarında Termal Modelleme için En İyi Uygulamalar
Diferansiyel denklem tabanlı termal modellerden güvenilir sonuçlar elde etmek için, bu yönergeleri takip edin:
- [FONT:0) Basite Başlayın[Dönemli Bir parametre modeliyle baskın zaman sabitlerini yakalamak için başlayın. Bu, sezgiyi inşa eder ve kritik parametreleri erken tanımlamaya yardımcı olur.
- [FONT:0)Ek karmaşıklığı yavaş yavaş): Basit model deneysel veri veya tasarım özelliklerini eşleştirmezken daha fazla düğüm veya doğrusal olmayanlığı tanıtın.
- [FONT:0) Güvenilir parametre değerleri [[Dönetici: Üretici veri tablolarından termal direnç ve kapasite elde etmek, ölçümlenen deneyler veya hakemli literatürden kaçınmak.
- [FONT:0] Ölçümlere karşı tarih[[Dönem: Her zaman simülasyon sonuçları en az birkaç sıcaklık ölçümleriyle karşılaştırır.Kesinlikle diskrepanzilerin% 5-10'u aştığının modelini %50'yi aşır.
- [FONT:0]En kötü durum koşullarını [Dönetici: Maksimum Orta Ortam Sıcaklıkta Simülasyonlar, en yüksek güç dağılımı ve yüksek soğutma (blok ve ilaçlar, toz ısısı) güvenlik marjlarını sağlamak için.
- [FONT:0)Belge varsayımları): Sınır koşulları, malzeme özellikleri ve basitleştirmeler hakkında tüm varsayımlar kaydetmek.Bu, tasarım tekrar gözden geçirildiğinde veya termal başarısızlıklar sorunken yardımcı olur.
Bu uygulamalara yemin ederek, mühendisler, gerçek dünya operasyonunun termal taleplerine dayanacak ürünleri temin edebilir, bağımsız tasarım kararlarını yapmak için diferansiyel denklem modellerine güvenebilirler.
Termal Modeling Future Trends in Thermal Modeling
Elektronikin artan güç yoğunluğu - miniaturizasyon ve daha yüksek performansa dayalı - daha sofistike termal modeller. Araştırmacılar, makine öğrenimini doğrudan ağlatmadan hızlandırırken, karmaşık geometrileri etkin bir şekilde ele almak için simülasyonu hızlandırıyor.
Uygulama mühendisi için, temel yetenek ısı akışını tanımlayan diferansiyel denklemleri nasıl ayarlayıp yorumlayacağımızı anlamaktadır. Matematik, araçlar geliştikçe bile.
Sonuç Sonuç Sonuç Sonuç Sonuç Sonuç Sonuç Sonuç
Elektriksel denklemlerle elektronik bileşenlerin termal davranışını modellemek, temel ısı denklemlerini geliştirmek ve uygun sayısal teknikleri uygulamak için temel ısı denklemlerini anlamak, enerji çözümlerinin belirlenmesi ve elektronik bileşenlerin operasyonel yaşamını genişletmek.