Simulation of S-parameters of general multilayer boxed PCBs with the method of moments and the scattering matrix algorithm

Dieser Beitrag stellt ein numerisch stabiles Werkzeug der Momentenmethode zur Simulation von S-Parametern von mehrschichtigen, abgeschirmten Leiterplatten vor, das eine S-Matrix-Formulierung zur Herleitung der vollständigen dyadischen Greenschen Funktion mit verschiedenen Basisfunktionen zur Modellierung sowohl transversaler als auch longitudinaler Ströme kombiniert.

Das Wesentliche

Stellen Sie sich vor, Sie sind ein Architekt, der einen Wolkenkratzer entwirft. Doch statt aus Beton und Stahl besteht Ihr Gebäude aus Schichten aus Kunststoff und dünnen Kupferblechen, die wie ein Sandwich gestapelt sind. Dies ist eine Leiterplatte (PCB), das Gehirn fast jedes elektronischen Geräts.

Bevor Sie diesen Wolkenkratzer tatsächlich bauen, möchten Sie wissen: Wird der Stromfluss vom obersten zum untersten Stockwerk reibungslos verlaufen? Wird er stecken bleiben oder auf seltsame Weise zurückprallen?

In der realen Welt müssten Sie einen Prototyp bauen, ihn testen und falls er versagt, ihn abreißen und von vorne beginnen. Das ist teuer und langsam. Daher nutzen Ingenieure Computersimulationen, um das Design virtuell zu „testen". Diese Arbeit stellt eine neue, intelligentere Methode vor, um diese Simulationen durchzuführen.

Hier ist die Aufschlüsselung ihrer Methode unter Verwendung einfacher Analogien:

1. Das Problem: Ein lauter, überfüllter Raum

Stellen Sie sich vor, Sie versuchen, ein Flüstern in einem Raum voller Echos zu hören.

• Die Leiterplatte ist der Raum. Sie besteht aus vielen Schichten (Dielektrika) und Metallblechen (Leitern).
• Das Signal ist das Flüstern (Elektrizität).
• Die Herausforderung: Wenn Elektrizität durch diese Schichten wandert, prallt sie von den Wänden und den Metallblechen ab. Um genau vorherzusagen, wie sich das Signal verhält, müssen Sie berechnen, wie jede einzelne Welle mit jeder anderen Welle interagiert.

Traditionell ist das Berechnen dieser „Echos" (mathematisch Green-Funktionen genannt) wie der Versuch, jedes einzelne Sandkorn an einem Strand zu zählen. Es erfordert enorme Rechenleistung und Zeit, besonders wenn sich die Signalquelle und der Empfänger nahe beieinander befinden. Die Mathematik wird unübersichtlich, instabil und langsam.

2. Die Lösung: Die „Streuungs-Matrix" (Der magische Spiegel)

Die Autoren schlagen eine neue Methode vor, um diese Echos mithilfe der sogenannten S-Matrix (Streuungs-Matrix) zu behandeln.

Stellen Sie sich die Leiterplatte als eine Reihe von Spiegeln und Fenstern vor, die übereinander gestapelt sind.

• Der alte Weg: Sie berechnen den Weg eines Lichtstrahls, indem Sie jeden einzelnen Abprall von jeder einzelnen Oberfläche einzeln nachverfolgen. Das ist mühsam.
• Der neue Weg (S-Matrix): Anstatt jeden Abprall nachzuverfolgen, behandeln Sie jede Schicht als eine „Blackbox" mit einem spezifischen Regelwerk.
  • Wenn eine Welle auf die Oberseite von Schicht A trifft, sagt Ihnen das Regelwerk genau, wie viel zurückprallt und wie viel zu Schicht B weitergeht.
  • Sie kombinieren die Regelwerke von Schicht A, Schicht B und Schicht C, um das Regelwerk für das gesamte Gebäude zu erhalten.

Das ist wie ein Spiel „Stille Post", bei dem Sie nicht die ganze Geschichte kennen müssen; Sie müssen nur wissen, wie jede Person in der Kette die Nachricht verändert. Durch die Verwendung dieser „Regelwerke" (S-Matrizen) wird die Mathematik viel stabiler und leichter zu berechnen, selbst für komplexe, mehrschichtige Strukturen.

3. Die „Dächer" und „Impulse" (Die Bausteine)

Um die Elektrizität zu simulieren, muss der Computer die Metallbleche und Drähte in winzige Stücke zerlegen.

• Flache Metallbleche: Die Autoren verwenden Formen, die wie Dächer aussehen (eine flache Oberseite mit geneigten Seiten), um den Strom darzustellen, der über die flachen Metallschichten fließt.
• Vertikale Drähte (Vias): Leiterplatten haben oft winzige Drähte, die durch die Schichten hindurchstoßen, um die Oberseite mit der Unterseite zu verbinden. Die Autoren verwenden Impuls- und lineare Formen (wie einen flachen Block oder eine Rampe), um den Strom darzustellen, der durch diese Drähte auf und ab fließt.

Sie haben die exakten mathematischen Formeln entwickelt, um zu berechnen, wie diese „Dächer" und „Impulse" mit den „Echos" (den S-Matrix-Regeln) interagieren. Dies ermöglicht es dem Computer, eine riesige Gleichung zu erstellen, die das Verhalten der gesamten Platine vorhersagt.

4. Der Geschwindigkeitsschub (Die schnelle Fourier-Transformation)

Selbst mit der neuen „Regelwerk"-Methode muss der Computer immer noch Millionen von Berechnungen durchführen.

• Die Analogie: Stellen Sie sich vor, Sie haben eine riesige Tabellenkalkulation, bei der jede Zelle ausgefüllt werden muss. Das eine nach dem anderen zu tun, dauert ewig.
• Die Lösung: Die Autoren verwenden eine Technik namens FFT (Fast Fourier Transform). Stellen Sie sich dies als eine superschnelle Sortiermaschine vor. Anstatt jede einzelne Zelle einzeln zu prüfen, gruppiert die Maschine sie auf clevere Weise, um die Antwort fast augenblicklich zu finden. Dies macht die Simulation schnell genug, um für reale Designs praktisch nutzbar zu sein.

5. Der Beweis: Hat es funktioniert?

Die Autoren testeten ihre neue Methode an zwei Beispielen:

1. Ein Filter: Eine Standard-Komponente mit drei Kunststoffschichten und sechs Metallstreifen. Sie verglichen ihre Computerergebnisse mit bekannten Daten aus anderen Studien, und die Zahlen stimmten perfekt überein.
2. Der Filter mit Drähten: Sie fügten zwei vertikale Drähte (Vias) hinzu, die die Schichten verbinden. Dies ist ein schwierigeres Problem, da es Stromfluss auf und ab beinhaltet, nicht nur seitlich. Ihre Methode bewältigte dies erfolgreich und zeigte, wie die Drähte das Signal veränderten.

Das Fazit

Diese Arbeit erfindet keinen neuen Typ von Leiterplatte. Stattdessen erfindet sie einen besseren Rechner für Ingenieure.

Durch die Verwendung eines „Regelwerk"-Ansatzes (S-Matrix) zur Behandlung der komplexen Echos innerhalb der Platine und durch die Verwendung von „Dach"-Formen zur Abbildung der Elektrizität haben sie ein Simulationswerkzeug geschaffen, das:

• Stabiler ist: Es stürzt nicht ab oder liefert seltsame Zahlen, wenn die Dinge kompliziert werden.
• Intuitiver ist: Es ist leichter zu verstehen und zu programmieren als frühere Methoden.
• Schneller ist: Es nutzt Geschwindigkeits-Boost-Tricks, um Probleme schnell zu lösen.

Dies hilft Ingenieuren, bessere, zuverlässigere elektronische Geräte zu entwerfen, ohne so viele physische Prototypen bauen zu müssen.

Technische Zusammenfassung

Technische Zusammenfassung: Simulation von S-Parametern allgemeiner mehrschichtiger, abgeschirmter PCBs mit der Momentenmethode und dem Streumatrix-Algorithmus

Problemstellung
Der Beitrag adressiert die rechnerische Herausforderung der Modellierung abgeschirmter mehrschichtiger Leiterplatten (PCBs) zur Simulation von S-Parametern vor der physischen Fertigung. Während PCBs durch dünne, hochleitfähige Elemente (Metallisierungsflächen, Streifen und Durchkontaktierungen) gekennzeichnet sind, die in einen geschichteten dielektrischen Wirt eingebettet sind, sind traditionelle Volumen-Diskretisierungsmethoden ineffizient. Die auf Oberflächen- und Volumenintegralgleichungen angewandte Momentenmethode (MoM) ist der bevorzugte Ansatz, da sie die Diskretisierung auf die Leiter beschränkt. Allerdings stehen bestehende 2,5D-MoM-Formulierungen für abgeschirmte Strukturen vor zwei Hauptproblemen:

1. Komplexität der Greenschen Funktion: Die analytische Herleitung der vollen dyadischen Greenschen Funktion für einen geschichteten Hohlleiter, der sowohl transversale (horizontale) als auch longitudinale (vertikale) Ströme berücksichtigt, ist umständlich. Viele bestehende Arbeiten konzentrieren sich nur auf transversale Ströme oder erfordern komplexe Herleitungen für spezifische Schichtanzahlen.
2. Numerische Stabilität und Konvergenz: Die Greensche Funktion in einem abgeschirmten Hohlleiter umfasst typischerweise unendliche Doppelreihen, die langsam konvergieren, insbesondere wenn Quellen- und Beobachtungspunkte nahe beieinander liegen. Darüber hinaus erfordert die Konstruktion der MoM-Matrix die Berechnung von Überlappungsintegralen zwischen diesen Komponenten der Greenschen Funktion und Basisfunktionen, ein Prozess, dem oft ein einheitlicher, stabiler und intuitiver algorithmischer Rahmen fehlt.

Methodik
Die Autoren schlagen einen einheitlichen Rahmen vor, der die Momentenmethode mit dem Streumatrix-Formalismus (S-Matrix) kombiniert, um die Integralgleichung für das elektrische Feld (EFIE) für abgeschirmte PCBs zu lösen.

• Herleitung der Greenschen Funktion via S-Matrix: Anstatt explizite unendliche Reihen für die Greensche Funktion abzuleiten, drücken die Autoren die Feldkomponenten in Form von aggregierten Streumatrizen aus. Der geschichtete Hohlleiter wird als Folge homogener dielektrischer Schichten und Grenzflächen behandelt.

  • Elementarmatrizen: Die Methode definiert elementare S-Matrizen für Grenzflächen zwischen dielektrischen Schichten, Ausbreitung durch homogene Schichten und Randbedingungen (Perfekter elektrischer Leiter oder offener Hohlleiter).
  • Aggregierte Matrizen: Unter Verwendung des Redheffer-Sternprodukts konstruieren die Autoren drei Sätze aggregierter S-Matrizen: „Oben" ($S^A$), „Unten" ($S^B$) und „Paar" ($S^P$). Diese repräsentieren die Hohlleiterabschnitte oberhalb, unterhalb und zwischen den Quellen- und Beobachtungspunkten.
  • Stabilität: Die Formulierung gewährleistet numerische Stabilität, indem sie die Lösung in Form von Exponentialfaktoren mit positiven Längenkonstanten ausdrückt und so die numerische Instabilität vermeidet, die oft mit evaneszenten Moden bei der direkten Reihensummation verbunden ist.

• Basisfunktionen und Überlappungsintegrale:

  • Planare Metallisierungen: Ströme auf horizontalen Streifen und Flächen werden mittels Oberflächen-Dachziegel-Basisfunktionen (voll, Auf-Rampe und Ab-Rampe halbe Dachziegel) auf einem regulären Gitter modelliert.
  • Vertikale Verbindungen (Vias): Ströme auf vertikalen Durchkontaktierungen werden mittels Volumen-Impuls- und Linear-Basisfunktionen modelliert.
  • Analytische Integration: Der Beitrag liefert explizite analytische Ausdrücke für die Überlappungsintegrale (Reaktionsintegrale) zwischen den abgeleiteten Komponenten der Greenschen Funktion und diesen Basisfunktionen. Dies umfasst Integrale für Quellen- und Beobachtungspunkte in derselben Schicht, in verschiedenen Schichten sowie für Vias, die mehrere Schichten durchqueren.

• Matrixkonstruktion und Skalierung:

  • Das MoM-Linearsystem wird durch Projektion der Integralgleichung auf die gewählten Basisfunktionen konstruiert.
  • Um die gemischten physikalischen Dimensionen von Oberflächenströmen (A/m) und Volumenströmen (A/m²) zu handhaben, führen die Autoren eine spezifische Reskalierung der Matrixelemente und Quellvektoren ein. Dies stellt sicher, dass die finale Systemmatrix Einheiten in Ohm, der Spannungsvektor Einheiten in Volt und der Stromvektor Einheiten in Ampere hat.

• Beschleunigung:

  • Der Beitrag erkennt an, dass die unendlichen Reihen in der Greenschen Funktion langsam konvergieren. Um die Matrixfüllung und -lösung zu beschleunigen, setzen die Autoren Fast-Fourier-Transformations-(FFT-)Techniken ein. Durch Ausnutzung der regulären Gitterstruktur werden die Matrixelemente als Summen aus Block-Toeplitz- und Hankel-Matrizen ausgedrückt, was eine Matrix-Vektor-Multiplikation in $O(n \log n)$-Operationen ermöglicht.

Hauptbeiträge

1. Einheitlicher S-Matrix-Formalismus: Der Beitrag liefert einen transparenten und allgemeinen Algorithmus zur Konstruktion der vollständigen dyadischen Greenschen Funktion für beliebige Quellen- und Beobachtungspunkte in einem mehrschichtigen, abgeschirmten Hohlleiter. Die Funktion wird ausschließlich in Form von drei vorab berechneten Sätzen von S-Matrizen ($S^A, S^B, S^P$) ausgedrückt, was die Implementierung im Vergleich zu traditionellen Ansätzen mit Leitungs- oder Modalerweiterung vereinfacht.
2. Umfassende Unterstützung von Basisfunktionen: Die Arbeit leitet explizit die Überlappungsintegrale her, die erforderlich sind, um sowohl horizontale Metallisierungen (mittels Dachziegel-Funktionen) als auch vertikale Verbindungen (mittels Impuls- und Linear-Volumenfunktionen) innerhalb desselben S-Matrix-Rahmens zu modellieren.
3. Numerische Stabilität: Die abgeleiteten Ausdrücke erweisen sich für alle Schichtkonfigurationen, einschließlich solcher mit evaneszenten Moden, als numerisch stabil, bedingt durch die spezifische Anordnung der Exponentialterme.
4. Implementierung und Validierung: Die Methode wurde in C++ implementiert und gegen bestehende Literatur sowie einen modifizierten Testfall validiert.

Ergebnisse
Die Autoren validierten die Methode anhand zweier numerischer Beispiele:

1. Filter mit drei dielektrischen Schichten: Eine Struktur mit sechs Metallstreifen in einem Hohlleiter wurde simuliert. Die berechneten S-Parameter ($S_{11}$ und $S_{21}$) im Bereich von 1–4 GHz zeigten eine hervorragende Übereinstimmung mit Ergebnissen aus einer zuvor veröffentlichten Studie (referenziert als [82]), was die Genauigkeit der Herleitungen der Greenschen Funktion und der Überlappungsintegrale bestätigt.
2. Filter mit vertikalen Verbindungen: Das erste Beispiel wurde so modifiziert, dass es zwei vertikale Vias enthält, die verschiedene Metallschichten verbinden. Die Simulation demonstrierte erfolgreich die Fähigkeit der Methode, den Einfluss dieser vertikalen Verbindungen auf die S-Parameter zu modellieren, und validierte die Integration von Volumen-Basisfunktionen mit der S-Matrix-Greenschen Funktion.

Bedeutung und Behauptungen
Der Beitrag behauptet, dass, obwohl das allgemeine Problem der PCB-Modellierung nicht neu ist, der spezifische Beitrag in der alternativen Lösungsmethode liegt, die sich durch Folgendes auszeichnet:

• Transparenz und Eleganz: Der S-Matrix-Ansatz bietet eine intuitivere und prägnantere Möglichkeit zur Konstruktion der Greenschen Funktion im Vergleich zu bestehenden komplexen Herleitungen.
• Allgemeingültigkeit: Die Methode ist auf beliebige Schichtanzahlen anwendbar und kann sowohl horizontale als auch vertikale Ströme modellieren, ohne separate, unterschiedliche theoretische Rahmenwerke für jeweils zu erfordern.
• Einfache Implementierbarkeit: Der algorithmische Charakter des S-Matrix-Aufbaus (rekursive Kombination von Elementarmatrizen) macht den Ansatz straightforward in Software zu implementieren.
• Erweiterbarkeit: Während die aktuelle Implementierung Dachziegel- sowie Impuls-/Linear-Basen verwendet, stellen die Autoren fest, dass die abgeleiteten Ausdrücke der Greenschen Funktion mit anderen Basisfunktionen (z. B. RWG, sinusförmig) kombiniert werden können, um Objekte verschiedener Formen zu modellieren.

Die Autoren nehmen eine bescheidene Haltung ein und stellen fest, dass die aktuelle Implementierung eine gleichmäßige Stromverteilung entlang der Länge der Vias und dünne Metallisierungsschichten annimmt, die Methode jedoch leicht auf nicht gleichmäßige Ströme und Schichten mit endlicher Dicke erweiterbar ist. Die Arbeit beansprucht nicht, das Problem der Greenschen Funktionen für offene (nicht abgeschirmte) Medien zu lösen, sondern konzentriert sich strikt auf die abgeschirmte (eingehüllte) Hohlleiterumgebung.