Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лабораторный практикум АиРРВ.pdf
Скачиваний:
283
Добавлен:
09.02.2015
Размер:
1.33 Mб
Скачать

Вклад каждого вибратора в полное поле излучения ВК определяется амплитудой тока, наведенного на элементе. На рис. 6.3 приводится распределение амплитуды тока по излучателям для 15-элементной антенны. Видно, что наиболее сильно возбуждаются два первых директора и активный вибратор.

В результате синфазного сложения полей в осевом направлении антенна имеет довольно узкий лепесток ДН. На рис. 6.4 приводится диаграмма направленности 15-элементной ВК, на рис. 6.5 ее сечения в Е- и Н- плоскостях.

6.2. Математическое описание вибраторных антенных решеток

Для проектирования антенных решеток и оптимизации их характеристик необходимо наличие точной и высокоэффективной в вычислительном отношении математической модели, поскольку многоэлементные ВК имеют большое число степеней свободы (длины вибраторов, расстояние между ними). Математическое описание можно свести к системе связанных интегральных уравнений относительно неизвестных функцийраспределения тока в каждом элементе.

6.2.1. Система связанных интегральных уравнений для многоэлементной антенной решетки ВК

Математическая модель ВК сводится к системе связанных интегральных уравнений, порядок которой совпадает с количеством элемента. Например, для пятиэлементной антенны, представленной на рис. 6.6. искомыми являются пять функций распределения токов в рефлекторе, активном элементе и трех директорах.

Основываясь на тех же допущениях, что и при выводе уравнения Поклингтона одиночного вибратора (т. е. наличие сильного скин-эффекта, приближение тонкого провода) [1], можно сформулировать систему уравнений относительно неизвестных токов на каждом элементе ВК:

65

z

Рис. 6.6

 

 

1(z) =

L1

 

LN

IN (z)K1N (z, z)dz;

Eстz

I1(z)K11(z, z)dz′+... +

 

 

L1

 

LN

 

 

 

 

 

 

 

2(z) =

L1

 

I1(z)K12(z, z)dz′+... +

LN

IN (z)K2N (z, z)dz;

Eстz

 

 

 

 

L1

LN

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

L1

 

LN

 

z

 

 

N (z) =

 

I1(z)K1N (z, z)dz′+... +

IN (z)KNN (z, z)dz

Eст

 

 

 

 

L1

 

LN

 

 

 

 

 

 

 

(6.1)

где Eстz n (z) – стороннее электрическое поле (продольная компонента), создаваемое внешним источ-

xниками поля, Im (z) – ток, протекающий по оси m- го элемента. Ядра системы (6.1) определяются на основе соотношений:

1

d 2...

 

2

eikRij

 

Kij (z, z ) = −

 

 

dz2

+k

 

...

.

(6.2)

 

 

 

iωε

 

 

 

R ji

 

 

 

 

 

 

 

 

 

 

 

 

Здесь – Rii

= (z

 

2

+a

2

 

– расстояние между точкой интегрирования на

 

z )

 

 

оси

и точкой

 

наблюдения

на поверхности i-го вибратора, Rij = R ji =

=

2

+dij

2

 

– расстояние между точкой интегрирования на оси i-го

(z z )

 

 

вибратора и точкой наблюдения на поверхности j-го вибратора, из (6.2) вид-

 

 

но, что: Kij (z, z ) = K ji (z, z ).

 

 

В случае одного активного элемента левые части системы уравнений

(6.1) имеют вид:

 

 

 

 

 

0

 

 

 

 

 

 

Eстz k(z)

 

 

 

 

 

 

...

 

 

 

0

 

 

 

 

66

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

Физический смысл системы интегральных уравнений (6.1) заключается в том, что каждая строчка – это запись граничного условия для касательной компоненты электрического поля на поверхности вибратора. Интегралы в правой части любой строки системы (6.1) – это вклад в поле на поверхности одного из вибраторов от каждого из N элементов антенной решетки, таким образом, общее число интегралов в системе равно квадрату числа элементов решетки.

Полученная система (6.1) описывает в самосогласованной постановке систему из N вибраторов с учетом взаимного влияния их друг на друга. В частном случае, при отсутствии взаимной связи (например, в случае значительного междуэлементного расстояния, т. е. при d → ∞), ядра системы (6.1) стремятся к нулю. Легко убедиться, что в этом случае система распадается на N независимых систем относительно тока в каждом уединенном вибраторе (т.е. на N независимых уравнений Поклингтона).

6.2.2. Решение системы связанных ИУ

Система (6.1) является обобщением уравнения Поклингтона для одиночного элемента и решается в пять этапов:

1-й этап. Выбираем систему базисных функций f1n (z), f2n (z), ..., fNn (z) по которым раскладывается предполагаемое решение на первом, втором и N-м вибраторах.

При этом возможна ситуация, в которой число базисных функций на каждом вибраторе неодинаково:N1 NN , а их вид на разных вибраторах разный. Получается система уравнений

 

N1

 

I1(z) = I1n f1n (z)

 

 

n=1

 

 

N2

 

 

 

(z) = I2n f2n (z)

 

I2

(6.3)

 

n=1

 

 

...

 

 

 

 

NN

 

 

 

 

IN (z) = INn fNn (z)

 

 

n=1

 

67

2-й этап. Подставляем разложение (6.3) в исходную систему интегральных уравнений (6.1) и, меняя порядок суммирования и интегрирования, получаем:

 

 

N

L1

Eстz

1(z) = 1

I1n

 

 

n=1

L

 

 

 

1

 

 

N

L1

Eстz

2

(z) = 1

I1n

 

 

n=1

L

 

 

 

1

...

 

 

 

 

 

N1

L1

 

 

Eстz

N(z) = I1n

 

 

n=1 L

 

 

 

1

f1n(z ')K11(z, z ')dz '+...+NN IN n

L

 

 

N fNn(z ')K1N (z, z ')dz '

n=1

LN

 

f1n(z ')K21(z, z ')dz '+...+NN IN n

L

 

N

fNn(z ')K2N (z, z ')dz '

n=1

LN

 

f1n(z ')KN1(z, z ')dz '+...+NN IN n

L

 

N

fNn(z ')KNN (z, z ')dz '

n=1

 

LN

 

(6.4)

3-й этап. Формируем систему линейных алгебраических уравнений на основе метода Галеркина. На этом этапе конкретизируем и упростим ситуацию, будем считать, что ток на каждом вибраторе описывается тремя базисными функциями. Далее, последовательно умножив правую и левую части уравнений системы (6.4) на проекционные функции fm (z) , получим систему линейных алгебраических уравнений (СЛАУ) M = 3N порядка относительно неизвестных базисных коэффициентов на каждом вибраторе:

68

L1

Eстz

1(z) fm(z)dz =

N =3

I1n

 

 

L1

 

L1

 

 

 

fn

(z ') fm(z)K11(z, z ')dz 'dz +...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

N =3

 

 

 

 

 

 

 

L1

 

 

LN

 

 

n=1

 

 

 

L1L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

... +

I

 

 

 

 

 

 

 

 

f

 

 

(z ') f

 

 

 

(z)K

 

 

 

 

 

 

(z, z ')dz 'dz ;

 

1m 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2

n=1

 

 

NnL1LN

 

 

 

 

 

 

 

 

 

 

 

 

1N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E z

 

 

 

 

 

(z) f

 

 

(z)dz

=

 

N =3

I

 

 

 

 

L2

 

 

 

L1

 

f

 

 

(z ') f

 

 

 

(z)K

 

 

 

 

(z, z ')dz 'dz +...

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

m

 

 

 

 

L2

ст 2

 

 

 

 

 

 

 

LN

 

 

 

 

 

n=1

1nL2

 

L1

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

 

 

 

 

 

 

 

 

... +

N =3

INn

 

L2

 

 

 

 

fn(z ') fm(z)K2N (z, z ')dz 'dz;

 

1m 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

n=1

 

 

 

 

 

 

L2

LN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

LN

E z

 

 

 

 

 

(z)

f

 

 

 

(z)dz

=

N =3

I

 

 

 

LN

 

 

 

 

L1

 

f

 

 

(z ') f

 

 

(z)K

 

 

 

(z, z ')dz 'dz +...

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

m

N1

 

 

 

стN

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

1n

LN L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

LN

 

 

 

 

 

 

 

 

 

LN

 

LN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

... +

N =3

INn

 

 

 

 

 

fn(z ') fm(z)KNN (z, z ')dz 'dz, 1m 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

LN LN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.5)

 

 

Систему (6.5) можно переписать в более наглядном виде, введя сквоз-

ную нумерацию базисных функций и токов:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u1 = i1z1,1 +i2z1,2 +i3z1,3 +i4z1,4 +i5z1,5 +i6z1,6 +... +iM z1,M ;

 

 

 

 

 

 

 

 

u

 

= i z

2,1

+i

 

z

2,2

+i

z

2,3

 

+i

z

2,4

 

 

+i

 

 

z

2,5

+i

 

z

2,6

+... +i

M

z

2,M

;

 

 

 

 

 

2

 

1

 

 

 

2

 

3

 

 

 

 

4

 

 

 

 

 

5

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

= i z

3,1

 

+i

 

 

z

 

 

 

 

+i

z

3,3

+i z

 

 

 

+i

z

3,5

+i

 

 

z

3,6

 

+... +i

M

z

 

 

;

 

 

 

 

 

 

3

 

1

 

 

 

2

3,2

 

3

 

 

 

 

4

3,4

 

 

 

5

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

3,M

 

 

 

 

 

 

 

 

u4 = i1z4,1 +i2z4,2 +i3z4,3 +i4z4,4 +i5z4,5 +i6z4,6 +... +iM z4,M ;

 

 

 

 

(6.6)

u

 

= i z

 

 

 

 

+i

 

 

z

 

 

 

 

+i

z

 

 

 

 

 

+i z

 

 

 

+i

z

 

 

 

 

+i

 

 

z

 

 

 

+... +i

M

z

 

 

;

 

 

 

 

 

5

 

1 5,1

 

 

 

2

5,2

 

3

5,3

 

 

 

4

5,4

 

 

 

5

 

 

5,5

 

 

 

6

 

5,6

 

 

 

 

 

 

 

 

5,M

 

 

 

 

 

 

 

 

u

 

= i z

 

 

 

 

+i

 

z

 

 

 

 

+i

z

 

 

 

 

 

+i

z

 

 

 

+i

 

z

 

 

 

 

+i

 

z

 

 

 

+... +i

M

z

 

 

 

;

 

 

 

 

 

 

6

 

1 6,1

 

 

 

2

 

6,2

3

 

6,3

 

 

 

4

 

6,4

 

 

 

5

 

 

6,5

 

 

 

6

6,6

 

 

 

 

 

 

 

 

6,M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

M

= i z

M

,1

+i z

M ,2

+i z

M ,3

+i

 

 

z

M ,4

+i z

M

,5

+i z

M ,6

+... +i

M

z

M

,M

.

 

1

 

 

 

2

 

 

 

 

3

 

 

 

 

4

 

 

 

5

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или в матричном виде – [U]=[Z][I]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь ui

 

(1 i M ) – амплитуды напряжения,

 

выражаемые через интегралы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ui

 

 

 

L

Eстz

i(z) fi(z)dz ,

от напряженности стороннего электрического поля:

 

=

 

L

причем вектор-столбец левых частей для ВК имеет вид

69

0

 

 

0

 

 

 

 

 

0

 

 

 

 

 

[U ]= 0

 

,

2U01

 

 

 

 

 

0

 

 

...

 

 

 

 

 

0

 

 

где 2U01 – напряжение дельта-источника, включенного в центре активного элемента (центральный сегмент активного элемента имеет 5-й порядковый номер), остальные элементы вектора [U] равны нулю; элементы матрицы [Z] –

L

L

обобщенных взаимных импедансов Zij =

fi (z ') f j (z) Kij (z, z ')dz 'dz ; [I] –

L L

вектор неизвестных базисных коэффициентов.

4-й этап. Решаем систему уравнений (6.6) любым известным способом, например методом Гаусса. В результате решения будут найдены базисные коэффициенты токов.

 

 

y

 

 

 

0

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

60

 

 

 

 

–60

 

 

 

–10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

–20

 

120

 

–120

180

Рис. 6.7

5-й этап. По найденным коэффициентам на основе (6.3) восстанавливается ток в каждом вибраторе, затем находится входное сопротивление активного вибратора и ДН всей системы.

Результаты решения системы и вычисления ДН 5-элементной ВК (рис.6.6) представлены на рис. 6.7.

70