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

Dit artikel presenteert een numeriek stabiel Method of Moments-hulpmiddel voor het simuleren van S-parameters van meerschichtige ingekaste PCB's door een S-matrixformalisme te combineren om de volledige dyadische Green-functie af te leiden met diverse basisfuncties om zowel transversale als longitudinale stromen te modelleren.

De kern

Stel je voor dat je een architect bent die een wolkenkrabber ontwerpt. Maar in plaats van beton en staal, is je gebouw gemaakt van lagen plastic en dunne koperplaten, opgestapeld als een sandwich. Dit is een Printed Circuit Board (PCB), het brein van bijna elk elektronisch apparaat.

Voordat je deze wolkenkrabber daadwerkelijk bouwt, wil je weten: Zal de elektriciteit soepel stromen van de bovenste verdieping naar de onderste? Zal het vastlopen of op vreemde manieren terugkaatsen?

In de echte wereld zou je een prototype moeten bouwen, het testen, en als het faalt, het moeten afbreken en opnieuw beginnen. Dat is duur en traag. Daarom gebruiken ingenieurs computersimulaties om het ontwerp virtueel te "testen". Dit artikel presenteert een nieuwe, slimmere manier om die simulaties uit te voeren.

Hier is de uitleg van hun methode met eenvoudige analogieën:

1. Het Probleem: Een Ruim met Lawaai en Drukte

Stel je voor dat je probeert een fluistering te horen in een kamer vol echo's.

• De PCB is de kamer. Het heeft vele lagen (diëlektrica) en metalen platen (geleiders).
• Het Signaal is de fluistering (elektriciteit).
• De Uitdaging: Wanneer elektriciteit door deze lagen beweegt, kaatst het af tegen de wanden en de metalen platen. Om precies te voorspellen hoe het signaal zich gedraagt, moet je berekenen hoe elke enkele golf met elke andere golf interacteert.

Traditioneel is het berekenen van deze "echo's" (wiskundig Green-functies genoemd) als proberen elk zandkorreltje op een strand te tellen. Het vereist enorme rekenkracht en tijd, vooral wanneer de signaalbron en de luisteraar dicht bij elkaar zitten. De wiskunde wordt rommelig, instabiel en traag.

2. De Oplossing: De "Strooiingsmatrix" (De Magische Spiegel)

De auteurs stellen een nieuwe manier voor om met deze echo's om te gaan met behulp van de S-matrix (Strooiingsmatrix) methode.

Stel je de PCB voor als een reeks spiegels en ramen die op elkaar zijn gestapeld.

• Oude Manier: Je berekent het pad van een lichtstraal door elke enkele botsing tegen elk oppervlak afzonderlijk te traceren. Dat is vermoeiend.
• De Nieuwe Manier (S-matrix): In plaats van elke botsing te traceren, behandel je elke laag als een "zwarte doos" met een specifiek regelboek.
  • Als een golf de bovenkant van Laag A raakt, vertelt het regelboek je precies hoeveel er terugkaatst en hoeveel er naar Laag B doorgaat.
  • Je combineert de regelboeken van Laag A, Laag B en Laag C om het regelboek voor het hele gebouw te krijgen.

Dit is als het spel "telefoon" waarbij je niet het hele verhaal hoeft te kennen; je hoeft alleen te weten hoe elke persoon in de keten het bericht verandert. Door deze "regelboeken" (S-matrices) te gebruiken, wordt de wiskunde veel stabieler en makkelijker te berekenen, zelfs voor complexe, meerlagige structuren.

3. De "Daken" en "Pulsen" (De Bouwblokken)

Om de elektriciteit te simuleren, moet de computer de metalen platen en draden in kleine stukjes breken.

• Vlakke Metalen Platen: De auteurs gebruiken vormen die eruitzien als daken (een plat dak met hellende zijden) om de stroom voor te stellen die over de vlakke metalen lagen stroomt.
• Verticale Draden (Via's): PCB's hebben vaak kleine draden die door de lagen heen prikken om de bovenkant met de onderkant te verbinden. De auteurs gebruiken puls en lineaire vormen (zoals een vlak blok of een helling) om de stroom voor te stellen die omhoog en omlaag door deze draden stroomt.

Ze hebben de exacte wiskundige formules bedacht om te berekenen hoe deze "daken" en "pulsen" interageren met de "echo's" (de S-matrix regels). Hierdoor kan de computer een enorme vergelijking opbouwen die het gedrag van het hele bord voorspelt.

4. De Snelheidssprint (De Snelle Fourier-transformatie)

Zelfs met de nieuwe "regelboek"-methode moet de computer nog steeds miljoenen berekeningen uitvoeren.

• De Analogie: Stel je voor dat je een enorme spreadsheet hebt waarbij elke cel moet worden ingevuld. Het één voor één doen duurt eeuwen.
• De Oplossing: De auteurs gebruiken een techniek genaamd FFT (Snelle Fourier-transformatie). Denk hierbij aan een supersnelle sorteermachine. In plaats van elke enkele cel afzonderlijk te controleren, groepeert de machine ze op een slimme manier om het antwoord bijna direct te vinden. Dit maakt de simulatie snel genoeg om praktisch toepasbaar te zijn voor real-world ontwerpen.

5. Het Bewijs: Heeft het Gewerkt?

De auteurs hebben hun nieuwe methode getest op twee voorbeelden:

1. Een Filter: Een standaard elektronisch component met drie lagen plastic en zes metalen strips. Ze vergeleken hun computerresultaten met bekende data uit andere studies, en de cijfers kwamen perfect overeen.
2. Het Filter met Draden: Ze voegden twee verticale draden (via's) toe die de lagen met elkaar verbonden. Dit is een moeilijker probleem omdat het stroom omvat die omhoog en omlaag beweegt, niet alleen zijwaarts. Hun methode slaagde hierin, en liet zien hoe de draden het signaal veranderden.

De Conclusie

Dit artikel bedenkt geen nieuw type printplaat. In plaats daarvan bedenkt het een betere rekenmachine voor ingenieurs.

Door een "regelboek"-benadering (S-matrix) te gebruiken om de complexe echo's binnen het bord te behandelen, en door "dak"-vormen te gebruiken om de elektriciteit in kaart te brengen, hebben ze een simulatiehulpmiddel gecreëerd dat:

• Stabieler is: Het crasht niet of geeft geen rare getallen als dingen complex worden.
• Intuïtiever is: Het is makkelijker te begrijpen en te programmeren dan eerdere methoden.
• Sneller is: Het gebruikt snelheidssprint-trucs om problemen snel op te lossen.

Dit helpt ingenieurs betere, betrouwbaardere elektronische apparaten te ontwerpen zonder dat ze evenveel fysieke prototypes hoeven te bouwen.

Technische samenvatting

Technische Samenvatting: Simulatie van S-parameters van Algemene Multilayer Geboxte PCB's met de Methode der Momenten en het Scattering-Matrixalgoritme

Probleemstelling
Het artikel behandelt de computationele uitdaging van het modelleren van afgeschermde multilayer Printed Circuit Boards (PCB's) om S-parameters te simuleren voordat er fysiek wordt gefabriceerd. Hoewel PCB's worden gekenmerkt door dunne, hooggeleidende elementen (metallisatievlekken, strips en vias) die zijn ingebed in een gestratificeerde diëlektrische host, zijn traditionele volumediskretisatiemethoden inefficiënt. De Methode der Momenten (MoM) toegepast op oppervlakte- en volume-integraalvergelijkingen is de voorkeursbenadering, omdat deze de diskretisatie beperkt tot de geleiders. Echter, bestaande 2.5D MoM-formuleringen voor afgeschermde structuren staan voor twee primaire moeilijkheden:

1. Complexiteit van de Groenfunctie: Het afleiden van de volledige dyadische Groenfunctie voor een gelaagde golfgeleider die rekening houdt met zowel transversale (horizontale) als longitudinale (verticale) stromen is analytisch omslachtig. Veel bestaande werken focussen alleen op transversale stromen of vereisen complexe afleidingen voor specifieke aantallen lagen.
2. Numerieke Stabiliteit en Convergentie: De Groenfunctie in een afgeschermde golfgeleider omvat doorgaans oneindige dubbele reeksen die langzaam convergeren, met name wanneer bron- en observatiepunten dicht bij elkaar liggen. Bovendien vereist het construeren van de MoM-matrix het berekenen van overlapintegralen tussen deze componenten van de Groenfunctie en basisfuncties, een proces dat vaak een gebrek vertoont aan een unificerend, stabiel en intuïtief algoritmisch raamwerk.

Methodologie
De auteurs stellen een unificerend raamwerk voor dat de Methode der Momenten combineert met de Scattering-Matrix (S-matrix) formaliteit om de Elektrisch-Veld-Integraalvergelijking (EFIE) op te lossen voor afgeschermde PCB's.

• Afleiding van de Groenfunctie via S-Matrix: In plaats van expliciete oneindige reeksen voor de Groenfunctie af te leiden, drukken de auteurs de veldcomponenten uit in termen van geaggregeerde scattering-matrices. De gelaagde golfgeleider wordt behandeld als een reeks homogene diëlektrische lagen en interfaces.

  • Elementaire Matrices: De methode definieert elementaire S-matrices voor interfaces tussen diëlektrische lagen, voortplanting door homogene lagen en randvoorwaarden (Perfecte Elektrische Geleider of open golfgeleider).
  • Geaggregeerde Matrices: Met behulp van het Redheffer-sterproduct construeren de auteurs drie sets geaggregeerde S-matrices: "Boven" ($S^A$), "Onder" ($S^B$) en "Paar" ($S^P$) matrices. Deze vertegenwoordigen respectievelijk de golfgeleidersecties boven, onder en tussen de bron- en observatiepunten.
  • Stabiliteit: De formulering zorgt voor numerieke stabiliteit door de oplossing uit te drukken in termen van exponentiële factoren met positieve lengteconstanten, waardoor de numerieke instabiliteit wordt vermeden die vaak gepaard gaat met evanescente modi bij directe reekssommatie.

• Basisfuncties en Overlapintegralen:

  • Planaire Metallisaties: Stromen op horizontale strips en vlekken worden gemodelleerd met oppervlakte-rooftop-basisfuncties (volledige, opwaartse en neerwaartse halve rooftops) op een regelmatig rooster.
  • Verticale Interconnecties (Vias): Stromen op verticale vias worden gemodelleerd met volume-puls- en lineaire basisfuncties.
  • Analytische Integratie: Het artikel levert expliciete analytische uitdrukkingen voor de overlapintegralen (reactieintegralen) tussen de afgeleide componenten van de Groenfunctie en deze basisfuncties. Dit omvat integralen voor bron- en observatiepunten in dezelfde laag, verschillende lagen, en voor vias die meerdere lagen kruisen.

• Matrixconstructie en Schaling:

  • Het MoM-lineaire systeem wordt geconstrueerd door de integraalvergelijking te projecteren op de gekozen basisfuncties.
  • Om de gemengde fysieke dimensies van oppervlaktstromen (A/m) en volumestromen (A/m²) te hanteren, voeren de auteurs een specifieke herschaling in van de matrixelementen en bronvectoren. Dit zorgt ervoor dat de uiteindelijke systeemmatrix eenheid Ohm heeft, de spanningsvector eenheid Volt en de stroomvector eenheid Ampère.

• Versnelling:

  • Het artikel erkent dat de oneindige reeksen in de Groenfunctie langzaam convergeren. Om het vullen van de matrix en de oplossing te versnellen, passen de auteurs Fast Fourier Transform (FFT)-technieken toe. Door gebruik te maken van de structuur van het regelmatige rooster, worden de matrixelementen uitgedrukt als sommen van blok-Toeplitz- en Hankel-matrices, waardoor matrix-vectorvermenigvuldiging kan worden uitgevoerd in $O(n \log n)$ operaties.

Belangrijkste Bijdragen

1. Unificerende S-Matrix Formaliteit: Het artikel biedt een transparant en algemeen algoritme om de volledige dyadische Groenfunctie te construeren voor willekeurige bron- en observatiepunten in een multilayer afgeschermde golfgeleider. De functie wordt uitsluitend uitgedrukt in termen van drie vooraf berekende sets S-matrices ($S^A, S^B, S^P$), wat de implementatie vereenvoudigt in vergelijking met traditionele transmissielijn- of modale expansiebenaderingen.
2. Uitgebreide Ondersteuning voor Basisfuncties: Het werk leidt expliciet de overlapintegralen af die nodig zijn om zowel horizontale metallisaties (met behulp van rooftop-functies) als verticale interconnecties (met behulp van puls- en lineaire volumefuncties) binnen hetzelfde S-matrixraamwerk te modelleren.
3. Numerieke Stabiliteit: De afgeleide uitdrukkingen blijken numeriek stabiel te zijn voor alle laagconfiguraties, inclusief die met evanescente modi, vanwege de specifieke rangschikking van exponentiële termen.
4. Implementatie en Validatie: De methode is geïmplementeerd in C++ en gevalideerd tegen bestaande literatuur en een gewijzigde testcase.

Resultaten
De auteurs hebben de methode gevalideerd met twee numerieke voorbeelden:

1. Filter met Drie Diëlektrische Lagen: Een structuur met zes metalen strips in een golfgeleider werd gesimuleerd. De berekende S-parameters ($S_{11}$ en $S_{21}$) over het bereik van 1–4 GHz toonden uitstekende overeenkomst met resultaten uit een eerder gepubliceerd onderzoek (verwezen naar als [82]), wat de nauwkeurigheid bevestigt van de afleidingen van de Groenfunctie en overlapintegralen.
2. Filter met Verticale Interconnecties: Het eerste voorbeeld werd gewijzigd om twee verticale vias op te nemen die verschillende metalen lagen verbinden. De simulatie toonde succesvol het vermogen van de methode aan om de impact van deze verticale interconnecties op de S-parameters te modelleren, waardoor de integratie van volumebasisfuncties met de S-matrix-Groenfunctie wordt gevalideerd.

Betekenis en Beweringen
Het artikel beweert dat hoewel het algemene probleem van PCB-modellering niet nieuw is, de specifieke bijdrage ligt in de alternatieve oplossingsmethode die wordt gekenmerkt door:

• Transparantie en Elegantie: De S-matrixbenadering biedt een meer intuïtieve en beknopte manier om de Groenfunctie te construeren in vergelijking met bestaande complexe afleidingen.
• Generaliteit: De methode is toepasbaar op willekeurige aantallen lagen en kan zowel horizontale als verticale stromen modelleren zonder dat er aparte, distincte theoretische raamwerken voor elk nodig zijn.
• Gemak van Implementatie: Het algoritmische karakter van de S-matrixconstructie (recursieve combinatie van elementaire matrices) maakt de benadering rechttoe rechtaan te implementeren in software.
• Uitbreidbaarheid: Hoewel de huidige implementatie rooftop- en puls/lineaire bases gebruikt, merken de auteurs op dat de afgeleide uitdrukkingen van de Groenfunctie kunnen worden gecombineerd met andere basisfuncties (bijvoorbeeld RWG, sinusvormig) om objecten van verschillende vormen te modelleren.

De auteurs houden een bescheiden standpunt aan, waarbij zij opmerken dat de huidige implementatie een uniforme stroomverdeling langs de lengte van vias en dunne metallisatielagen veronderstelt, maar dat de methode gemakkelijk uitbreidbaar is naar niet-uniforme stromen en lagen met eindige dikte. Het werk claimt niet het probleem van open (ongeachermde) media-Groenfuncties op te lossen, maar focust strikt op de afgeschermde (geboxte) golfgeleideromgeving.