مقدمه

در تحلیل سیستم‌های قدرت و به‌ویژه در مطالعه رفتار گذراهای الکترومغناطیسی، مدل‌سازی دقیق خطوط انتقال از اهمیت بسیار بالایی برخوردار است. در شرایط واقعی، پارامترهای الکتریکی خط انتقال شامل مقاومت سری، راکتانس سری و سوسپتانس شنت، همگی تابعی از فرکانس هستند. این وابستگی فرکانسی باعث می‌شود که مدل‌های ساده مانند مدل π-بخشی کلاسیک نتوانند در بازه‌های وسیع فرکانسی پاسخ دقیقی ارائه دهند.

در این شبیه‌سازی، یک مدل پیشرفته برای خط انتقال وابسته به فرکانس معرفی می‌شود که بر پایه استخراج پارامترهای فرکانسی، تقریب‌سازی تابعی با استفاده از ابزارهای RF Toolbox و پیاده‌سازی مدل Universal Line Model (ULM) در محیط Simscape توسعه یافته است. هدف اصلی این مدل، افزایش دقت در تحلیل گذراهای سریع و سیگنال‌های پهن‌باند است.

مدلسازی پارامترهای وابسته به فرکانس

در این رویکرد، ابتدا پارامترهای خط انتقال به‌صورت تابعی از فرکانس استخراج می‌شوند. این پارامترها شامل موارد زیر هستند:

– مقاومت سری واحد طول R(ω)
– راکتانس سری واحد طول X(ω)
– سوسپتانس شنت واحد طول B(ω)
– فرکانس متناظر freq
– طول کل خط len

این پارامترها معمولاً برای یک خط هوایی در ارتفاع مشخص و با فرض شرایط زمین ایده‌آل یا با مقاومت ویژه مشخص محاسبه می‌شوند. در این شبیه‌سازی فرض شده است که اثر پوست (Skin Effect) و مقاومت زمین ناچیز است تا تمرکز بر رفتار فرکانسی ذاتی خط حفظ شود.

با داشتن این داده‌ها، امپدانس و ادمیتانس سری و شنت به صورت زیر تعریف می‌شوند:

Z(ω) = R(ω) + jX(ω)
Y(ω) = jB(ω)

این روابط پایه‌ای‌ترین ساختار برای تحلیل انتشار موج در خط انتقال هستند.

پارامترهای انتشار و ادمیتانس مشخصه

یکی از مهم‌ترین مفاهیم در تحلیل خطوط انتقال وابسته به فرکانس، ادمیتانس مشخصه Yc و تابع انتشار H است. این پارامترها رفتار موج الکترومغناطیسی در طول خط را تعیین می‌کنند.

ادمیتانس مشخصه به صورت زیر تعریف می‌شود:

Yc = √(Y / Z)

که در آن Z و Y به ترتیب امپدانس و ادمیتانس واحد طول هستند.

همچنین ثابت انتشار Γ به صورت زیر تعریف می‌شود:

Γ = √(Y Z)

این مقدار شامل بخش حقیقی (تضعیف) و موهومی (فاز) است. سرعت انتشار موج نیز از رابطه زیر استخراج می‌شود:

V = ω / Im(Γ)

تابع انتقال یا تابع انتشار خط نیز به صورت زیر بیان می‌شود:

H = e^(-Γ · len)

این تابع نشان‌دهنده میزان تضعیف و تأخیر سیگنال در طول خط انتقال است و نقش کلیدی در مدل‌سازی دقیق رفتار گذرا دارد.

تقریب تابعی با استفاده از Rational Fitting

یکی از چالش‌های اصلی در مدل‌سازی خطوط انتقال وابسته به فرکانس، تبدیل داده‌های فرکانسی به یک مدل قابل پیاده‌سازی در حوزه زمان است. برای این منظور از روش Rational Fitting استفاده می‌شود که در RF Toolbox در دسترس است.

در این روش، تابع ادمیتانس مشخصه به صورت مجموعی از کسرهای جزئی نمایش داده می‌شود:

Yc(s) = Σ (Residue_i / (s – Pole_i))

در این رابطه:
– Pole_i قطب‌های سیستم هستند
– Residue_i ضرایب باقیمانده هستند
– n مرتبه تقریب است

در این شبیه‌سازی معمولاً از تقریب مرتبه هشتم استفاده می‌شود که تعادل مناسبی بین دقت و پیچیدگی محاسباتی ایجاد می‌کند.

نتیجه این تقریب آن است که رفتار فرکانسی پیچیده خط به مجموعه‌ای از سیستم‌های مرتبه اول قابل پیاده‌سازی در حوزه زمان تبدیل می‌شود.

تقریب تابع انتشار و در نظر گرفتن تأخیر زمانی

تابع انتشار H نیز مشابه ادمیتانس مشخصه باید به فرم قابل محاسبه تبدیل شود. برای افزایش دقت و کاهش مرتبه تقریب، ابتدا تأخیر زمانی τ از تابع استخراج می‌شود:

H(s) = Hk(s) · e^(-sτ)

در این رابطه Hk بخش بدون تأخیر است که برای تقریب استفاده می‌شود. سپس داریم:

Hk(s) = Σ (HkResidue_i / (s – HkPole_i))

این روش باعث می‌شود که بخش تأخیر به‌صورت یک بلوک مستقل در مدل در نظر گرفته شود و پیچیدگی تقریب کاهش یابد.

مدل Universal Line Model در Simscape

مدل Universal Line Model یا ULM یکی از دقیق‌ترین روش‌های تحلیل خطوط انتقال در حوزه گذراهای الکترومغناطیسی است. در این مدل، خط انتقال به صورت شبکه‌ای از جریان‌ها و ولتاژها در دو سر خط تعریف می‌شود.

در این ساختار، متغیرهای اصلی شامل موارد زیر هستند:

– V1, V2: ولتاژهای دو سر خط
– I1, I2: جریان‌های ورودی
– Ish: جریان شنت
– Iaux: جریان کمکی ناشی از انتشار موج
– Irfl: جریان بازتابی
– H: تابع انتشار

روابط بنیادی مدل به صورت زیر تعریف می‌شوند:

I1 = Ish,1 – Iaux,2
I2 = Ish,2 – Iaux,1

همچنین:

Ish = Yc · V
Iaux = 2H · Irfl
Irfl = (I + Ish) / 2

این روابط نشان می‌دهند که جریان در هر سر خط ترکیبی از مؤلفه‌های محلی و اثرات انتشار از سر دیگر خط است.

تبدیل به حوزه زمان و پیاده‌سازی عددی

برای پیاده‌سازی در Simscape، لازم است روابط حوزه لاپلاس به حوزه زمان تبدیل شوند. برای ادمیتانس مشخصه داریم:

Ish = Σ Wi
dWi/dt = Pi Wi + Ri V(t)

که در آن Pi قطب‌ها و Ri ضرایب تقریب هستند.

به طور مشابه برای تابع انتشار:

Iaux = 2 Σ Xi
dXi/dt = Pi Xi + Ri Irfl(t – τ)

این معادلات نشان می‌دهند که مدل نهایی شامل مجموعه‌ای از معادلات دیفرانسیل مرتبه اول با تأخیر زمانی است که به‌راحتی در Simscape قابل پیاده‌سازی هستند.

نتایج شبیه‌سازی

در اولین سناریوی شبیه‌سازی، منبع ولتاژ یک سیگنال سینوسی 60 هرتز تولید می‌کند. در این حالت، مدل کلاسیک π-بخشی که بر اساس فرکانس نامی 60 هرتز تنظیم شده است، رفتار مناسبی در حالت ماندگار نشان می‌دهد. با این حال، مدل وابسته به فرکانس دقت بالاتری در گذراهای اولیه و تغییرات سریع دارد.

در دومین سناریو، سیگنال ورودی شامل یک موج 60 هرتز به همراه مدولاسیون 10 کیلوهرتز است. در این حالت اختلاف عملکرد دو مدل به‌وضوح مشاهده می‌شود. مدل π-بخشی قادر به دنبال کردن تغییرات سریع نیست و دچار خطای فاز و دامنه می‌شود، در حالی که مدل وابسته به فرکانس رفتار بسیار دقیق‌تری در بازه پهن‌باند ارائه می‌دهد.

جمع‌بندی

در این شبیه‌سازی، نشان داده شد که استفاده از مدل‌های وابسته به فرکانس برای تحلیل خطوط انتقال، به‌ویژه در کاربردهای گذرا و سیگنال‌های پهن‌باند، ضروری است. با استفاده از روش Rational Fitting و پیاده‌سازی Universal Line Model در Simscape، می‌توان یک مدل دقیق و در عین حال قابل محاسبه برای خطوط انتقال ایجاد کرد.

این رویکرد نسبت به مدل‌های کلاسیک مانند π-Section، دقت بسیار بالاتری در تحلیل رفتار دینامیکی خطوط ارائه می‌دهد و برای کاربردهای پیشرفته در مهندسی قدرت و تحلیل EMC بسیار مناسب است.