🔬 applied physics

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

تقدم هذه الورقة أداة "طريقة العزوم" (Method of Moments) مستقرة عددياً لمحاكاة معاملات S (S-parameters) للوحات الدوائر المطبوعة الصندوقية متعددة الطبقات، وذلك عبر الجمع بين صيغة مصفوفة S لاستخلاص دالة "جرين" ثنائية الالتواء (dyadic Green's function) كاملة، وبين دوال أساسية متنوعة لنمذجة التيارات المستعرضة والطولية على حد سواء.

المؤلفون الأصليون: A. O. Makarenko, P. Zheglova, R. Gaponenko, R. V. Salimov, R. I. Tikhonov, A. A. Shcherbakov

نُشر 2026-05-19

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: A. O. Makarenko, P. Zheglova, R. Gaponenko, R. V. Salimov, R. I. Tikhonov, A. A. Shcherbakov

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

تخيل أنك مهندس معماري تصمم ناطحة سحاب. ولكن بدلاً من الخرسانة والفولاذ، يتكون مبناك من طبقات من البلاستيك وشرائح رقيقة من النحاس، مكدسة فوق بعضها البعض مثل الشطيرة. هذه هي لوحة الدوائر المطبوعة (PCB)، وهي بمثابة العقل المدبر لكل جهاز إلكتروني تقريباً.

قبل أن تبدأ في بناء هذه الناطحة بالفعل، تريد أن تعرف: هل ستتدفق الكهرباء بسلاسة من الطابق العلوي إلى الطابق السفلي؟ هل ستتعثر أو ترتد بطرق غريبة؟

في العالم الحقيقي، سيتعين عليك بناء نموذج أولي، واختباره، وإذا فشل، عليك هدمه والبدء من جديد. هذا أمر مكلف وبطيء. لذلك، يستخدم المهندسون عمليات المحاكاة الحاسوبية لـ "اختبار" التصميم افتراضياً. تقدم هذه الورقة البحثية طريقة جديدة وأكثر ذكاءً لتشغيل عمليات المحاكاة تلك.

إليك تفصيل طريقتهم باستخدام تشبيهات بسيطة:

1. المشكلة: غرفة صاخبة ومزدحمة

تخيل أنك تحاول سماع همسة في غرفة مليئة بالصدى.

اللوحة (PCB) هي الغرفة. تتكون من طبقات عديدة (مواد عازلة) وصفائح معدنية (موصلات).
الإشارة هي الهمسة (الكهرباء).
التحدي: عندما تتحرك الكهرباء عبر هذه الطبقات، فإنها ترتد عن الجدران والصفائح المعدنية. للتنبؤ بكيفية سلوك الإشارة بدقة، يجب عليك حساب كيفية تفاعل كل موجة مع كل موجة أخرى.

تقليدياً، حساب هذه "الأصداء" (والتي تسمى رياضياً دوال غرين - Green's functions) يشبه محاولة عد كل حبة رمل على الشاطئ. إنها تتطلب قدرة حوسبة هائلة ووقتاً طويلاً، خاصة عندما يكون مصدر الإشارة والمستمع قريبين من بعضهما البعض. تصبح الرياضيات معقدة، وغير مستقرة، وبطيئة.

2. الحل: "مصفوفة التشتت" (المرآة السحرية)

يقترح المؤلفون طريقة جديدة للتعامل مع هذه الأصداء باستخدام ما يسمى بطريقة مصفوفة S (S-matrix/Scattering Matrix).

تخيل لوحة الدوائر المطبوعة كسلسلة من المرايا والنوافذ المكدسة فوق بعضها البعض.

الطريقة القديمة: تقوم بتتبع مسار شعاع الضوء من خلال تتبع كل ارتداد عن كل سطح بشكل فردي. إنه أمر مرهق.
الطريقة الجديدة (S-matrix): بدلاً من تتبع كل ارتداد، تعامل كل طبقة كـ "صندوق أسود" له كتاب قواعد محدد.
- إذا اصطدمت موجة بأعلى الطبقة (أ)، يخبرك كتاب القواعد بالضبط كم سيبلغ مقدار الارتداد وكم سينفذ إلى الطبقة (ب).
- تقوم بدمج كتب القواعد الخاصة بالطبقة (أ)، والطبقة (ب)، والطبقة (ج) للحصول على كتاب القواعد للمبنى بأكمله.

هذا يشبه لعب لعبة "الهاتف المكسور" حيث لا تحتاج لمعرفة القصة كاملة؛ بل تحتاج فقط لمعرفة كيف يغير كل شخص في السلسلة الرسالة. ومن خلال استخدام هذه "كتب القواعد" (S-matrices)، تصبح الرياضيات أكثر استقراراً وأسهل في الحساب، حتى بالنسبة للهياكل المعقدة متعددة الطبقات.

3. "الأسطح" و"النبضات" (لبنات البناء)

لمحاكاة الكهرباء، يحتاج الكمبيوتر إلى تقسيم الصفائح والأسلاك المعدنية إلى قطع صغيرة.

الصفائح المعدنية المسطحة: يستخدم المؤلفون أشكالاً تشبه الأسطح (سطح مستوٍ بجوانب مائلة) لتمثيل التيار المتدفق عبر الطبقات المعدنية المسطحة.
الأسلاك الرأسية (Vias): غالباً ما تحتوي لوحات الدوائر المطبوعة على أسلاك صغيرة تخترق الطبقات لربط الأعلى بالأسفل. يستخدم المؤلفون أشكال النبض (pulse) والخطية (linear) (مثل كتلة مسطحة أو منحدر) لتمثيل التيار المتدفق صعوداً وهبوطاً في هذه الأسلاك.

لقد توصلوا إلى الصيغ الرياضية الدقيقة لحساب كيفية تفاعل هذه "الأسطح" و"النبضات" مع "الأصداء" (قواعد S-matrix). وهذا يسمح للكمبيوتر ببناء معادلة ضخمة تتنبأ بسلوك اللوحة بأكملها.

4. دفعة السرعة (تحويل فوريه السريع)

حتى مع طريقة "كتاب القواعد" الجديدة، لا يزال على الكمبيوتر إجراء ملايين الحسابات.

التشبيه: تخيل أن لديك جدول بيانات ضخم حيث يجب ملء كل خلية فيه. القيام بذلك واحدة تلو الأخرى يستغرق وقتاً طويلاً جداً.
الحل: يستخدم المؤلفون تقنية تسمى FFT (تحويل فوريه السريع). فكر في هذا كآلة فرز فائقة السرعة. بدلاً من فحص كل خلية على حدٍ، تقوم الآلة بتجميعها بطريقة ذكية لإيجاد الإجابة بشكل فوري تقرياً. وهذا يجعل المحاكاة سريعة بما يكفي لتكون عملية للتصاميم الواقعية.

5. الإثبات: هل نجح الأمر؟

اختبر المؤلفون طريقتهم الجديدة على مثالين:

المرشح (Filter): مكون إلكتروني قياسي يحتوي على ثلاث طبقات من البلاستيك وست شرائط من المعدن. قارنوا نتائج الكمبيوتر الخاصة بهم مع بيانات معروفة من دراسات أخرى، وكانت الأرقام متطابقة تماماً.
المرشح مع الأسلاك: أضافوا سلكين رأسيين (vias) يربطان الطبقات. هذه مشكلة أصعب لأنها تتضمن حركة التيار صعوداً وهبوطاً، وليس فقط من جانب إلى آخر. تعاملت طريقتهم مع هذا بنجاح، مما أظهر كيف غيرت الأسلاك الإشارة.

الخلاصة

هذه الورقة البحثية لا تخترع نوعاً جديداً من لوحات الدوائر المطبوعة. بدلاً من ذلك، هي تخترع آلة حاسبة أفضل للمهندسين.

من خلال استخدام نهج "كتاب القواعد" (S-matrix) للتعامل مع الأصداء المعقدة داخل اللوحة، وباستخدام أشكال "الأسطح" لرسم خريطة الكهرباء، أنشأوا أداة محاكاة تتميز بـ:

أكثر استقراراً: لا تتعطل أو تعطي أرقاماً غريبة عندما تصبح الأمور معقدة.
أكثر وضوحاً: من الأسهل فهمها وبرمجتها مقارنة بالطرق السابقة.
أسرع: تستخدم حيل تسريع لحل المشكلات بسرعة.

هذا يساعد المهندسين على تصميم أجهزة إلكترونية أفضل وأكثر موثوقية دون الحاجة إلى بناء الكثير من النماذج الأولية المادية.

ملخص تقني: محاكاة معاملات S للوحات الدوائر المطبوعة (PCBs) متعددة الطبقات العامة ذات الصندوق المغلق باستخدام طريقة العزوم وطريقة مصفوفة التشتت

بيان المشكلة
تتناول الورقة البحثية التحدي الحسابي لنمذجة لوحات الدوائر المطبوعة (PCBs) متعددة الطبقات والمحمية لمحاكاة معاملات S قبل التصنيع الفعلي. وبينما تتميز لوحات الـ PCB بعناصر رقيقة وعالية الموصلية (رقع التعدين، والشرائط، والفتحات/vias) مغمورة داخل وسط عازل طبقي، فإن طرق التجزئة الحجمية التقليدية تكون غير فعالة. وتعد طريقة العزوم (MoM) المطبقة على المعادلات التكاملية السطحية والحجمية هي النهج المفضل، لأنها تحصر التجزئة في الموصلات فقط. ومع ذلك، تواجه صياغات طريقة العزوم ثنائية الأبعاد (2.5D) الحالية للهياكل المحمية صعوبتين رئيسيتين:

تعقيد دالة غرين (Green's Function): إن اشتقاق دالة غرين ثنائية المصفوفة (dyadic Green's function) الكاملة لموجّه موجي متعدد الطبقات، والذي يأخذ في الاعتبار التيارات المستعرضة (الأفقية) والطولية (الرأسية)، هو أمر معقد تحليلياً. وتركز العديد من الأعمال الحالية فقط على التيارات المستعرضة أو تتطلب اشتقاقات معقدة لعدد محدد من الطبقات.
الاستقرار العددي والتقارب: تتضمن دالة غرين في الموجّه الموجي المحمي عادةً سلاسل مزدوجة لانهائية تتقارب ببطء، خاصة عندما تكون نقاط المصدر والمراقبة قريبة من بعضها البعض. علاوة على ذلك، فإن بناء مصفوفة طريقة العزوم يتطلب حساب تكاملات التداخل بين مكونات دالة غرين وهذه الدوال الأساسية، وهي عملية تفتقر غالباً إلى إطار خوارزمي موحد ومستقر وبديهي.

المنهجية
يقترح المؤلفون إطاراً موحداً يجمع بين طريقة العزوم (MoM) وصيغة مصفوفة التشتت (S-matrix) لحل معادلة التكامل للمجال الكهربائي (EFIE) للوحات الـ PCB المحمية.

اشتقاق دالة غرين عبر مصفوفة التشتت (S-Matrix): بدلاً من اشتقاق سلاسل لانهائية صريحة لدالة غرين، يعبر المؤلفون عن مكونات المجال من خلال مصفوفات تشتت مجمعة. يتم التعامل مع الموجّه الموجي الطبقي كسلسلة من الطبقات العازلة المتجانسة والواجهات.
- المصفوفات الأولية: تحدد الطريقة مصفوفات S أولية للواجهات بين الطبقات العازلة، وللانتقال عبر الطبقات المتجانسة، وللشروط الحدية (الموصل المثالي للكهرباء PEC أو الموجّه المفتوح).
- المصفوفات المجمعة: باستخدام ضرب "ريدهيفر" (Redheffer star product)، يبني المؤلفون ثلاث مجموعات من مصفوفات S المجمعة: المصفوفات "العلوية" ( $S^A$ )، و"السفلية" ( $S^B$ )، و"الزوجية" ( $S^P$ ). وتمثل هذه المصفوفات مقاطع الموجّه الموجي فوق وتحت وبين نقاط المصدر والمراقبة، على التوالي.
- الاستقرار: تضمن هذه الصيغة الاستقرار العددي من خلال التعبير عن الحل باستخدام عوامل أسية ذات ثوابت طول موجبة، مما يتجنب عدم الاستقرار العددي المرتبط غالباً بالأنماط المتلاشية (evanescent modes) عند استخدام الجمع المباشر للسلاسل.
الدوال الأساسية وتكاملات التداخل:
- التعدينات المستوية: يتم نمذجة التيارات على الشرائط والرقع الأفقية باستخدام دوال "روفتوب" السطحية (الكاملة، والمنحدر الصاعد، والمنحدر الهابط) على شبكة منتظمة.
- الوصلات الرأسية (Vias): يتم نمذجة التيارات على الفتحات الرأسية باستخدام دوال نبضية وخطية حجمية.
- التكامل التحليلي: توفر الورقة تعبيرات تحليلية صريحة لتكاملات التداخل (تكاملات التفاعل) بين مكونات دالة غرين المشتقة وهذه الدوال الأساسية. ويشمل ذلك تكاملات لنقاط المصدر والمراقبة في نفس الطبقة، وفي طبقات مختلفة، وكذلك للفتحات (vias) التي تعبر طبقات متعددة.
بناء المصفوفة والتحجيم:
- يتم بناء نظام طريقة العزوم الخطي عن طريق إسقاط المعادلة التكاملية على الدوال الأساسية المختارة.
- للتعامل مع الأبعاد الفيزيائية المختلطة لتيارات السطح (A/m) والتيارات الحجمية (A/m²)، يقدم المؤلفون إعادة تحجيم محددة لعناصر المصفوفة ومتجهات المصدر. يضمن ذلك أن مصفوفة النظام النهائية بوحدات الأوم (Ohms)، ومتجه الجهد بوحدات الفولت (Volts)، ومتجه التيار بوحدات الأمبير (Amperes).
التسريع:
- تقر الورقة بأن السلاسل اللانهائية في دالة غرين تتقارب ببطء. ولتسريع ملء المصفوفة وحلها، يستخدم المؤلفون تقنيات تحويل فوريه السريع (FFT). ومن خلال استغلال هيكل الشبكة المنتظمة، يتم التعبير عن عناصر المصفوفة كمجموع من مصفوفات "تويتز" (Toeplitz) و"هانكل" (Hankel) الكتلية، مما يسمح بإجراء ضرب المصفوفات في عمليات قدرها $O(n \log n)$ .

المساهمات الرئيسية

صيغة مصفوفة التشتت الموحدة: توفر الورقة خوارزمية شفافة وعامة لبناء دالة غرين ثنائية المصفوفة الكاملة لأي نقاط مصدر ومراقبة في موجّه موجي متعدد الطبقات ومحمي. يتم التعبير عن الدالة حصرياً من خلال ثلاث مجموعات محسوبة مسبقاً من مصفوفات S ( $S^A, S^B, S^P$ )، مما يسهل التنفيذ مقارنة بنهج خطوط النقل أو التوسع النمطي التقليدي.
دعم شامل للدوال الأساسية: تشتق هذه الدراسة صراحةً تكاملات التداخل المطلوبة لنمذجة كل من التعدينات الأفقية (باستخدام دوال الروفتوب) والوصلات الرأسية (باستخدام الدوال النبضية والخطية الحجمية) ضمن نفس إطار عمل مصفوفة التشتت.
الاستقرار العددي: أظهرت التعبيرات المشتقة أنها مستقرة عددياً لجميع تكوينات الطبقات، بما في ذلك الأنماط المتلاشية، وذلك بسبب الترتيب المحدد للمصطلحات الأسية.
التنفيذ والتحقق: تم تنفيذ الطريقة بلغة ++C والتحقق من صحتها مقابل الأدبيات الموجودة وحالة اختبار معدلة.

النتائج
تحقق المؤلفون من صحة الطريقة باستخدام مثالين عدديين:

مرشح (Filter) بثلاث طبقات عازلة: تمت محاكاة هيكل يحتوي على ستة شرائط معدنية داخل موجّه موجي. أظهرت معاملات S المحسوبة ( $S_{11}$ و $S_{21}$ ) عبر النطاق 1-4 جيجاهرتز توافقاً ممتازاً مع نتائج دراسة منشورة سابقاً (المشار إليها بـ [82])، مما يؤكد دقة اشتقاقات دالة غرين وتكاملات التداخل.
مرشح (Filter) مع وصلات رأسية: تم تعديل المثال الأول ليشمل فتحتين (vias) رأسيتين تربطان طبقات معدنية مختلفة. نجحت المحاكاة في إظهار قدرة الطريقة على نمذجة تأثير هذه الوصلات الرأسية على معاملات S، مما أثبت صحة دمج الدوال الأساسية الحجمية مع دالة غرين الخاصة بمصفوفة التشتت.

الأهمية والادعاءات
تدعي الورقة أنه بينما ليست مشكلة نمذجة الـ PCB عامة أمراً جديداً، فإن المساهمة المحددة تكمكمن في طريقة الحل البديلة التي تتميز بـ:

الشفافية والأناقة: يوفر نهج مصفوفة التشتت طريقة أكثر بديهية وإيجازاً لبناء دالة غرين مقارنة بالاشتقاقات المعقدة الموجودة.
العمومية: الطريقة قابلة للتطبيق على أي عدد من الطبقات ويمكنها نمذجة التيارات الأفقية والرأسية دون الحاجة إلى أطر نظرية منفصلة ومتميزة لكل منها.
سهولة التنفيذ: إن الطبيعة الخوارزمية لبناء مصفوفة التشتت (الدمج المتكرر للمصفوفات الأولية) تجعل النهج سهلاً في التنفيذ برمجياً.
القابلية للتوسع: على الرغم من أن التنفيذ الحالي يفترض توزيع تيار منتظم على طول الفتحات وطبقات التعدين الرقيقة، إلا أن المؤلفين يشيرون إلى أن تعبيرات دالة غرين المشتقة يمكن دمجها مع دوال أساسية أخرى (مثل RWG أو الدوال الجيبية) لنمذجة أجسام ذات أشكال متنوعة.

يتبنى المؤلفون موقفاً متواضعاً، ملاحظين أن التنفيذ الحالي يفترض توزيع تيار منتظم على طول الفتحات وطبقات التعدين الرقيقة، لكن الطريقة قابلة للتوسع بسهولة لتشمل التيارات غير المنتظمة والطبقات ذات السماكة المحددة. لا تدعي الورقة حل مشكلة دوال غرين للوسط المفتوح (غير المحمي)، بل تركز حصرياً على بيئة الموجّه الموجي المحمي (الصندوق).