Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

курсач / курсач 1 / курсовая_1_усвчилп_брв2201

.pdf
Скачиваний:
0
Добавлен:
13.05.2026
Размер:
3.28 Mб
Скачать

Министерство цифрового развития, связи и массовых коммуникаций Российской Федерации

Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования

Московский технический университет связи и информатики

___________________________________________________________________

Факультет

«Радио и телевидение»

Кафедра

«Техническая электродинамика и антенны (ТЭДиА)»

Курсовая работа по дисциплине «Устройства СВЧ и линии передачи» по теме: «Основные уравнения электродинамики»

Бригада 10, вариант 1

Выполнил Студент группы БРВ2201 _________________________

Проверил

 

К.т.н., доцент

_________________________ Гайнутдинов Т.А.

Москва 2024

1. СОДЕРЖАНИЕ

1. СОДЕРЖАНИЕ

2. ИСХОДНЫЕ ДАННЫЕ

3. РАСЧЁТНАЯ ЧАСТЬ

3.1. РАСЧЁТ СОСТАВЛЯЮЩИХ ВЕКТОРОВ ЭЛЕКТРИЧЕСКОГО И

МАГНИТНОГО ПОЛЕЙ 3.2. РАСЧЁТ ДИАПАЗОНА ЧАСТОТ, ПРИ КОТОРОМ ВОЛНА БЕЖИТ ВДОЛЬ ОСИ Z

3.3. РАСЧЁТ ВЫРАЖЕНИЙ ДЛЯ МГНОВЕННЫХ ЗНАЧЕНИЙ СОСТАВЛЯЮЩИХ ВЕКТОРОВ

3.3.1. ДЛЯ ЧАСТОТЫ f2

3.3.2. ДЛЯ ЧАСТОТЫ f1

3.4. ПРОВЕРКА ВЫПОЛНЕНИЯ ГРАНИЧНЫХ УСЛОВИЙ

3.5 РАСЧЁТ МАКСИМАЛЬНЫХ ЗНАЧЕНИЙ ПЛОТНОСТЕЙ ТОКОВ 3.6. РАСЧЁТ СРЕДНЕГО ЗА ПЕРИОД ПОТОКА ЭНЕРГИИ ЧЕРЕЗ

ПОПЕРЕЧНОЕ СЕЧЕНИЕ ВОЛНОВОДА

3.7. РАСЧЁТ РАЗЛИЧНЫХ СКОРОСТЕЙ 3.8 ЧЕРТЕЖИ СТРУКТУРЫ ВЕКТОРНЫХ ЛИНИЙ ПОЛЕЙ И ЭПЮРЫ

ТОКОВ НА СТЕНКАХ ВОЛНОВОДА

2. ИСХОДНЫЕ ДАННЫЕ

Исходные данные представлены ниже: рисунок 2.1 со всеми исходными данными.

Рисунок 2.1 – Исходные данные

Таким образом, исходные данные следующие: выражения для составляющих векторов поля (формулы 1-3); характеристики волновода, частоты волн, а также амплитуда вектора напряжённости магнитного поля (таблица 2.1).

Hxm

i β a

H0

 

2 π x

e

-i β z

 

 

2 π

sin

a

 

 

 

 

,

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

Hym 0,

 

 

 

 

 

 

 

(2)

Hzm H0

 

2 π x

e

-i β

z

.

 

(3)

cos

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.1 – характеристики волновода и поля

H0

εr μr σ

a

b

f1

f2

λкр

(A÷m) (mm) (mm) (GHz) (GHz) (mm)

10

4

1

0

30

40

1.2

8

30

В таблице 2.1 находятся следующие величины:

H0=10 mA – амплитуда вектора напряжённости магнитного поля;

εr=4 – относительная диэлектрическая проницаемость среды, заполняющей

волновод;

μr=1 – относительная магнитная проницаемость среды, заполняющей

волновод; σ=0 – удельная электропроводность среды, заполняющей волновод;

a=30 mm – ширина волновода; b=40 mm – высота волновода; f1=1.2 GHz – частота;

f2=8 GHz – частота;

λкр=30 mm – критическая длина волны в волноводе.

3.РАСЧЁТНАЯ ЧАСТЬ

3.1.РАСЧЁТ СОСТАВЛЯЮЩИХ ВЕКТОРОВ ЭЛЕКТРИЧЕСКОГО И МАГНИТНОГО ПОЛЕЙ

Необходимо рассчитать комплексные амплитуды составляющих векторов магнитного и электрического полей, не указанных в исходных данных, то есть, в рамках задачи, это комплексные амплитуды составляющих вектора

электрического поля E .

Так как комплексные амплитуды составляющих вектора магнитного поля H даны по условиям задачи, рассчитать комплексные амплитуды составляющих

вектора электрического поля E можно следующим образом:

Em

rot Hm

(4)

,

 

i ω εa

 

где Em – комплексная амплитуда вектора электрического поля; rot Hm – ротор от комплексной амплитуды магнитного поля;

i – мнимая единица;

ω – угловая частота волны;

εa 8.85 10-12 mF – абсолютная диэлектрическая проницаемость среды.

Ротор от комплексной амплитуды можно вычислить следующим образом:

Azm-

 

 

rot Am x

∂y

∂z

Aym

 

 

 

 

 

 

 

 

+ y

∂z

Aym-

∂x

Azm

 

(5)

 

 

 

 

 

Aym-

 

 

 

+ z

∂x

∂y

Axm

 

 

 

 

 

 

 

где x , y , z – единичные вектора соответствующих осей системы координат (вектора-направления);

Axm, Aym, Azm – соответствующие составляющие вектора поля.

Таким образом для вычисления ротора сперва необходимо рассчитать частные производные данных составляющих комплексной амплитуды вектора напряжённости магнитного поля. Однако стоит заметить, что игрековая составляющая равна нулю, а значит и её частные производные.

А иксовая и зетовая составляющая не зависят от игрека, так что частные производные по игрек этих составляющих равны нулю. В итоге необходимо рассчитать только игрековое слагаемое ротора, а значит только те частные производные, которые в нём находятся.

 

 

H0 β

2

a e

-(1i

β

z)

 

2 x π

 

 

 

 

 

 

sin

a

 

Hxm β,a,H0,x,z →

 

 

 

 

 

2 π

 

 

,

∂z

 

 

 

 

 

 

 

 

 

 

 

 

H0

π e

-(1i

β

z)

2 x π

- 2

 

 

 

 

sin

a

 

Hzm β,a,H0,x,z →

 

 

 

 

 

 

a

 

 

 

.

∂x

 

 

 

 

 

 

 

 

 

 

 

 

Тогда ротор будет вычисляться по следующей формуле:

 

H0 β

2

a e

-(1i

β

z)

 

 

2 x π

 

 

 

 

 

 

 

 

 

sin

 

a

 

 

 

 

rot Hm y

 

 

 

 

 

 

 

 

2 π

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x π

 

 

 

 

 

 

π

e

-(1i

β z)

 

 

 

 

2 H0

 

 

 

 

sin

a

 

 

 

 

+y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и будет равен:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H0 y π

2

+H0

β

2

a

2

 

 

 

-(1i

β z)

 

2 x π

4

 

 

 

 

y e

 

 

 

sin

a

 

rot Hm

 

 

 

 

 

 

 

 

 

2 a π

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6)

(7)

(8)

(9)

Так как в роторе присутствует только игрековая составляющая, иксовая и зетовая составляющие комплексной амплитуды напряжённости электрического поля также будут равны нулю. Следовательно, необходимо рассчитать только игрековую составляющую по формуле (4):

 

 

 

H0 π

2

+H0

β

2

a

2

e

-(1i

β z)

 

 

2 x π

 

 

1

4

 

 

 

 

 

sin

 

a

 

 

Eym

 

 

 

 

 

 

 

2 a π

 

 

 

 

 

,

(10)

 

i ω εa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eym -

H0 i 4 π2 2

 

a2

 

2 π x

e

-i

β z

,

 

(11)

 

2 a π ω εa

 

 

sin

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Eym – игрековая составляющая комплексной амплитуды вектора напряжённости электрического поля.

Итоговая матрица для всех составляющих комплексной амплитуды каждого поле имеет вид:

 

 

 

 

 

 

Exm 0

 

 

 

 

 

 

 

 

 

 

 

 

 

-H0 i 4 π

2

2

a

2

sin

2

π x

 

 

 

E

ym

 

 

 

e-i β z

 

 

 

2 a π ω εa

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ezm 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12)

 

 

 

i β a

 

 

 

2 π x

 

 

 

 

 

 

 

Hxm

 

H0 sin

e

-i β z

 

 

 

2 π

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hym 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 π x

 

 

-i

β z

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

Hzm H0 cos

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.2. РАСЧЁТ ДИАПАЗОНА ЧАСТОТ, ПРИ КОТОРОМ ВОЛНА БЕЖИТ ВДОЛЬ ОСИ Z

Для расчёта диапазона частот, в котором рассматриваемое поле представляет собой волну, бегущую вдоль оси z, коэффициент фазы β должен быть действительным числом. Сам коэффициент фазы можно рассчитать по следующей формуле:

β

2 π

 

 

λ

2

(13)

λ

1-

 

,

 

 

 

λкр

 

 

где λ – длина волны в волноводе на произвольной частоте.

Для расчёта величины λ необходимо рассчитать скорость распространения волны внутри волновода cср, что можно сделать по следующей формуле:

c

ср

c

= 1.499 108

 

m ,

(14)

 

εr μr

 

s

 

 

 

 

 

 

 

 

 

где c – это скорость света в вакууме.

 

 

 

Зная cср, можно выразить величину λ:

 

 

 

 

 

 

λ cср÷f.

 

 

(15)

Для того, чтобы коэффициент фазы β был равен действительному числу, отношение длины волны на произвольной частоте к критической длине волны должно быть меньше или равно 1:

λ ≤1.

λкр

Тогда, после подстановки в формулу формулу для f:

f≥ cср .

λкр

(16)

(16) выражения (15), можно получить

(17)

То есть, после подстановки в (17) критической частоты, диапазон частот можно вычислить по следующим формулам (λкр a указано в дано задания):

 

fкр

cср

,

(18)

 

λкр

 

 

 

 

 

fкр

cср

=4.997 GHz.

(19)

λкр

 

 

 

 

 

Тогда диапазон частот, при котором рассматриваемое поле есть волна, бегущая вдоль оси z:

fкр≤f<∞,

4.997 GHz<f<. (20)

3.3. РАСЧЁТ ВЫРАЖЕНИЙ ДЛЯ МГНОВЕННЫХ ЗНАЧЕНИЙ СОСТАВЛЯЮЩИХ ВЕКТОРОВ

Необходимо получить выражения для мгновенных значений составляющих векторов напряжённости электрического и магнитного полей. Для этого необходимо взять действительную часть от произведения составляющей поля по координате и множителя ei ω t. Получение зависимостей необходимо осуществить с учётом следующих условий:

x a÷3=10 mm,

y b÷3=13.333 mm, (21) 0≤z≤2 Λ,

где Λ

λ2

=23.992 mm – длина волны в волноводе при f2.

1- λ2÷λкр 2

 

 

Так же необходимо построить графики зависимости получившихся выражений для мгновенных значений напряжённостей полей от z в пределах из (20) при следующих значениях времени: t1 0 s, t2 T÷4.

3.3.1. ДЛЯ ЧАСТОТЫ f2

Так как частота f2=8 GHz больше, чем fкр=4.997 GHz, то величина

коэффициента фазы β будет действительной. Тогда выражения для составляющих напряжённостей полей будут иметь вид (с учётом того, что нулевые составляющие можно не рассчитывать).

Для Eym:

Eymf(t,z) Re Eym ei ω t ,

(22)

E

 

(t,z) -

H0 4 π2

2

a2

sin

 

2 π x

sin(β z-ω t),

ymf

2

εa ω a π

 

a

 

 

 

 

 

 

 

 

где ω 2 π f2=? rads – угловая частота для f2,

β

2 π

 

λ2

 

2

rad

– коэффициент фазы для частоты f2.

 

1-

 

 

=261.887

 

 

λ2

 

 

λкр

 

m

 

Ниже, на рисунке 3.3.1.1 приставлены зависимости напряжённости зетовой составляющей напряжённости электрического поля при двух значениях времени t1 и t2.

Диапазон значений z выбран следующим:

z 0 mm,0.1 mm 2 Λ. (23)

На графике сплошной линией обозначена зависимость игрековой составляющей напряжённости электрического поля при значении времени t1=0 s, а пунктиром – при значении времени t2= 3.125 10-11 s.

1.2 10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 10³

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 10³

 

 

 

 

 

 

 

 

 

E

 

t ,z V

2 10³

 

 

 

 

 

 

 

 

 

 

4 10³

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

ymf

1

 

 

-2 10³0

5

10

15

20

25

30

35

40

45

50

 

 

m

-4 10³

 

 

 

 

 

 

 

 

 

E

 

t ,z

V

-6 10³

 

 

 

 

 

 

 

 

 

 

-8 10³

 

 

 

 

 

 

 

 

 

 

ymf

2

 

 

-1 10

 

 

 

 

 

 

 

 

 

 

 

 

m

-1.2 10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z (mm)

Рисунок 3.2.1.1 – Графики зависимости игрековой составляющей напряжённости электрического поля от z при разных t.

Для Hxm:

 

 

(

)

Re

 

 

i ω t

,

 

Hxmf

t,z

 

Hxm e

 

 

 

 

H0 β a

 

2 π x

 

(24)

Hxmf

(t,z)

 

sin(β z-ω t).

2 π

sin

a

 

 

 

 

 

 

 

 

Для Hzm:

Hzmf(t,z) Re Hzm ei ω t ,

Hzmf(t,z) H0 cos 2 π x cos(β z-ω t).

(25)

a

Ниже, на рисунках 3.3.1.2 и 3.3.1.3 приставлены зависимости напряжённости иксовой и зетовой составляющих напряжённости магнитного поля при двух значениях времени t1 и t2.

Соседние файлы в папке курсач 1