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

книги из ГПНТБ / Штагер, Е. А. Рассеяние волн на телах сложной формы

.pdf
Скачиваний:
58
Добавлен:
19.10.2023
Размер:
7.41 Mб
Скачать

иых элементов и представлены в виде таблицы чисел:

аі(Т/); 52(т/); •••

«МТ/); ': = 1>2,

(4.1)

где оп (у.) — величина

ЭПР п-го элемента

в направле­

нии 7,..

 

 

Задача заключается в том, .как объединить все а„(у,-) для получения значения ЭПР тела сложной формы в це­ лом. Существует по крайней мере три подхода к реше­ нию этой задачи. Два из них мы рассмотрим ниже, тре­ тий подход был изложен в § 2.4 и § 3.1.

Метод суммирования с учетом фаз отраженных полей (метод относительной фазы)

Допустим, что по каким-либо соображениям выбра­ но р направлении из сектора углов наблюдения ±уо- В простейшем случае величину сектора углов 2у0 делят на р частей и определяют ЭПР тела сложной формы лишь в направлениях, отличающихся друг от друга на величину 2у0/р. Обозначая одно из указанных направле­

ний в секторе углов наблюдения через у м о ж н о

запи­

сать ЭПР тела сложной формы в виде

 

Л '

 

2

 

с (Т/) = £

°Я 2 (т/) ехР [ — УФя (т/.)] *• = 1 , 2 ,

( 4 . 2 )

л = 1

 

 

где on(7;)

и

(7,-) — величина ЭПР и фаза поля,

отра­

женного от

 

элемента в направлении 7;, УѴ— общее чис­

ло элементов.

 

Расчеты по формуле (4.2) достаточно просто выпол­ няются, если ЭПР тела слабо изменяется в выбранном секторе углов наблюдения. В этом случае число на­ правлений р может быть выбрано равным одному, двум, а требуемая точность определения фазы отраженного поля — сравнительно невысокой.

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

9 З а к а з Л5’ 166

129

рактер. Это прежде всего относится к определению разно­ сти фаз полей, отраженных от отдельных отражателей. Разность фаз пропорциональна отношению расстояния между элементами к длине волны поля. Если ограни­ читься точностью определения разности фаз в пределах 10°, то расстояния между элементами будет необходимо измерять с точностью до сотых долей длины волны. Как правило, это технически невозможно, если длина волны поля выбрана в сантиметровом или дециметровом диапа­ зонах волн. Неточность же определения расстояний между элементами, а следовательно, и разности фаз отраженных от них полей і|->„— фгР приводит к угловому смещению выбросов ЭПР тела сложной формы. В ре­ зультате ошибка в определении ЭПР тела с выбранного направления может достигнуть сотен процентов, что недопустимо.

Выше был указан лишь один источник погрешно­ стей вычисления і|:„ (у,-). На самом деле их несколько,, причем некоторые слабо зависят от длины волны поля. К числу таковых в первую очередь относятся погрешно­ сти, связанные с приближенным характером определения фазы отраженного поля при помощи методов коротко­ волновой ассимптотикп. В отдельных случаях ошибки в определении і):„ (у,) могут достигать десятков градусов, что приводит к значительным погрешностям расчета ЭПР тела в целом.

Метод случайной фазы

Если в формуле (4.2) полагать, что а) фазы отра­ женных полей от элементов тела сложной формы яв­ ляются взаимно независимыми случайными величинами; б) распределение вероятности фп(у,) для всех п яв­ ляется равномерным в интервале 0 -т- 2л, то среднее зна­ чение ЭПР такого тела

N

(4.3)

где і = 1, 2, . . ., р.

При выполнении указанных условий может бытьопределена и величина среднеквадратичного отклонения

130

величины

ЭПР

тела сложной

формы от ее среднего

значения

(о):

N

о

 

 

 

 

 

 

д т Е

с«(п)

 

Так как

трудно

1 = 1 ,2 ,

(4.4)

определяемая

величина

фазы фя (у,)

полей, отраженных от отдельных элементов, в расчетные формулы не входит, расчеты по методу случайной фазы оказываются несложными, и этот метод нашел доста­ точно широкое распространение в инженерной практике

[63, 88, 89]. І\ тому же, если равенства (4.3) и (4.4) рас­ сматривать как приближенные, то условия а) и б) мо­ гут быть заменены более слабыми. В первом условии вместо независимости фаз отраженных полей достаточно потребовать их некоррелированность. Во втором же условии можно отказаться от требования одинаковости распределений ф„(уі) и их равномерности в интервале О— 2я; достаточно лишь, чтобы среднеквадратичное от­ клонение случайной величины фп(уі) превышало 2л радиан [61].

Однако вместе с упрощением расчетной стороны ме­ тода случайной фазы по сравнению с методом относи­ тельной фазы возникает ряд трудностей. Предположим, что по методу случайной фазы определены среднее и среднеквадратичное значения ЭПР конкретного тела сложной формы, наблюдаемого с направления у*. Если затем мы измерим ЭПР этого тела в лабораторных ус­ ловиях, то получим какое-то значение этой величины, которое обозначим стоЕсли случайные факторы на трассе распространения волны отсутствуют, то во будет представлять не случайную величину, которая с какой-то

долей вероятности может

содержаться в интервале

в + Л 4 2 1 / 2 { ст} , рассчитанном

по методу случайной фазы.

Значение уровня этой вероятности в каждом конкрет­ ном случае неизвестно, что затрудняет интерпретацию результатов расчета по-методу случайной фазы.

Большое значение имеет знание разброса величины ЭПР тела сложной формы, определяемого методом слу­ чайной фазы. Некоторое представление о его величине можно получить, полагая, что ЭПР элементов одинако­

Р*

131

вы (ст„ (у,) = сто) • Тогда искомый интервал достаточно просто вычисляется по формуле (4.4) и оказывается равным

с ± M V2{а } = с 0 [7V + / Л Д Л Г ^ Т У ] .

( 4 . 5 )

Как следует из (4.5), интервал возрастает пропорцио­ нально количеству элементов N. При N = 2 этот интер­ вал составляет сг0(0,59 ч- 3,4), а для N = 3 соответствен­ но сто (0,55 -г- 5,45) и т. д.

Для элементов, различающихся величиной своей ЭПР и состоящих из двух и трех разных участков локально­ го отражения, величины интервалов указаны в табл. 3 и 4. Там же для сравнения приведены минимаксные значе­ ния ЭПР тела сложной формы, рассчитанные по методу относительной фазы.

Данные табл. 3 и 4 свидетельствуют о том, что ме­ тод случайной фазы может быть использован лишь для оценки порядка величины ЭПР тела сложной формы. К такому заключению приходят авторы книги [88], от­ куда заимствованы эти таблицы.

Как было показано выше, объединение величин ЭПР элементов с целью вычисления ЭПР тела в целом, про­

водимое на основе

метода

относительной

и случайной

 

 

 

 

 

Т а б л и ц а

3

 

Р а с ч е т ы э ф ф е к т и в н о й п л о щ а д и р а с с е я н и я

 

 

 

 

д в у х э л е м е н т н о г о т е л а

 

 

 

 

 

 

Р е з у л ь т а т ы р а с ч е т а п о м е т о д у

 

 

 

 

 

 

с л у ч а й н о й ф а з ы

 

 

 

Э П Р

Э П Р

 

 

 

 

М и н и м а к с н ы с

п е р в о г о

в т о р о г о

 

 

 

 

 

 

И н т е р в а л з н а ч е н и и Э П Р

з н а ч е н и я

э л е м е н т а

э л е м е н т а

С р е д н я я Э П Р

 

(*з)

 

 

Э П Р

 

 

 

и О,))

Оі) + Г и , ’ р ( Т , ) ]

 

 

 

 

 

 

 

1

1

 

2

0,59-^3,4

 

0 1

О

 

 

 

 

 

 

1

 

1

1,5

 

2,5

0,77-4-4,2

 

0,05^4,9

1

2

 

3

1,0-4-5,0

 

0,17-4-5,8

1

4

 

5

2,2н-7,8

 

1,0-4-9

 

1

9

 

10

5,8-і-14,2

 

4,0-4-16

132

Т а б л и ц а 4

Р а с ч е т ы э ф ф е к т и в н о й п л о щ а д и р а с с е я н и я т р е х э л е м е н т н о г о т е л а

о г)

о 3)

г о 3 )

о г ( а

о г ( о

ь е (сг

е р в т а

т о р т а

р е т т а

п н

в н

т н

Э П Р э л е м е

Э П Р э л е м е

Э П Р э л е м е

1

1

1

1

1

9

1

1

49

1

2

2

1

2

9

1

2

49

1

4

4

1

4

9

1

4

49

Р е з у л ь т а т ы р а с ч е т а

п о м е т о д у с л у ч а й н о й ф а з ы

С р е д н я я

И н т е р в а л з н а ч е н и й Э П Р

Э П Р

 

- ы

° ( і ;) + I7 ли [" 0;)]

30,55+5,45

п4,84+17,2

51

36,9-Т-65 1

5

О •I- СО о

12

4,38+19,6

52

34,7+69,3

9

2,07+15,9

14

4,10+23,9

54

31,5+76,5

=0} CJ

5 * X S

S о

0+9-

1+25

25+81

0+15

0,36+29

21+89

0+25

0+36

16+100

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

133

4.2.Пример расчета

По методике, изложенной в предыдущем параграфе, был проведен расчет эффективной площади рассеяния реактивного самолета «Конвер 990». Ниже воспроизве­ ден расчет ЭПР для углов наблюдения самолета с «но­ совых» ракурсов в плоскости крыльев. Источником ис­ ходных данных о размерах и форме отдельных частей

самолета служит его

модель, выполненная

в

масштабе

1 : 25

(рис. 36) [88].

порядком

расчета

на

первом его

В

соответствии с

этапе

следует произвести расчленение

тела

сложной

формы на отдельные элементы.

Для углов

наблюдения

в небольшом секторе носовых ракурсов можно выделить следующие элементы: 1) двигатели, 2) крылья, 3) фюзе­ ляж, 4) хвостовое оперение.

На каждом из перечисленных элементов следует ука­ зать участки локального отражения, т. е. произвести еще дополнительное расчленение поверхности на более мел­ кие части. Результаты расчленения и последующей аппроксимации как самих элементов, так и их частей телами простой формы представлены в табл. 5. Данные этой таблицы показывают, что большинство элементов удается аппроксимировать телами простой формы, при­ ближенные аналитические выражения ЭПР которых известны. В силу этого следует считать известными ве­ личины ст5 -f- о14 для каждого значения угла наблюдения из сектора носовых курсовых углов. Значения щ н- сц

не выражаются

в аналитической форме,

однако

для

их определения

имеются подпрограммы

вычислений

на

цифровой вычислительной машине, составленные на основе данных работы [63] н др.

Заключительный этап расчетов состоит в объедине­ нии величин ЭПР отдельных элементов для определения результирующей ЭПР самолета. Если объединение про­ изводить по методу случайной фазы, то среднее значе­ ние ЭПР модели самолета в соответствии с (4.3) опре­ деляется как

11

 

°(Тг) = - j j S °Я(Т<): / = 1 , 2 , . . . , / ; .

(4.6)

п-1

 

134

 

 

 

Т а б л и ц а 5

 

Э П Р

Н п п м е н о о л н и е о т д е л ь н ы х

 

Ы а н м е н о в а н н е

о т д е л ь н ы х

 

ч а с т е й э л е м е н т а и т е л ,

П р и м е ч а н и е

э л е м е н т а

ч а с т е й

н х з а м е н я ю щ и х

 

 

э л е м е н т а

 

 

 

 

Двигатели

Со

3-1

 

3 0

Крылья

а 7

~9

° і о

Первыіі двигатель Втором двигатель Третий двигатель Четвертый двигатель

Часть передней кромки левого крыла заменяется круговым цилиндром, длиной 250 м м и радиусом

6мм

Такая же часть перед­

ней

кромки

правого

крыла

 

кромки

Средняя часть

левого

крыла заменяется

круговым цилиндром, дли­ ной 270 м м и радиусом

5м м

Такая же часть кромки

правого крыла Оконечная часть перед­

ней кромки левого крыла заменяется кругойым ци­ линдром, длиной 270 м м

ирадиусом 3 мм Такая же часть перед­

ней

кромки

правого

крыла

 

 

Величины определиются при помощи подпрограмм для ЦВМ

° 5 =

° 7 = ° 8

° » — а 1»

Фюзеляж

J I1

Передняя часть корпуса

13ß

Продолжение табл. 5

 

Э П Р

Н а и м е н о в а н и е о т д е л ь н ы х

 

Н а и м е н о в а н и е

о т д е л ь н ы х

 

ч а с т е й э л е м е н т а и т е л ,

П р и м е ч а н и е

э л е м е н т а

ч а с т е й

и х з а м е н я ю щ и х

 

 

э л е м е н т а

 

 

 

 

Хвостовое

оперение

3 12

С 13

3 І-1

самолета заменяется вытя­

 

нутым

сфероидом с полу­

 

осями,

равным

77

м м

 

и 330 м м

 

 

 

 

 

Передняя

кромка левой

3 12 = 3 13

горизонтальной

плоскости

 

заменяется

круговым

ци­

 

линдром, длиной 260

м м

 

м радиусом

4

м м .

Ци­

 

линдр ориентирован

под

4

углом

133°

к осп

фюзе­

ляжа

же

кромка

пра­

 

Такая

 

вой горизонтальной

пло­

 

скости

 

 

 

 

 

 

Передняя кромка верти­

 

кальной

плоскости

заме­

 

няется

круговым

цилин­

 

дром,

длиной

3 0 0

м м

 

ирадиусом 4 м м

Вто же время интервал предполагаемых значений ЭПР модели самолета будет равен

а(7,-)±{тѴИС« (Т ,)-

М

 

Н І

С«(Т/)

Л-1

Л= 1

1, 2, р. (4.7)

 

/ =

В формулах (4.6) и (4.7) символом у,- обозначены углы наблюдения из рассматриваемого сектора «носовых» углов наблюдения.

Значения ЭПР, определяемые по формулам (4.6) и (4.7), нанесены в виде графиков на рис. 37—39. Эти графики, заимствованные из работы [88], включают не только «носовые» углы наблюдения самолета, но и дру-

136

Рис. 37. Зависимость ЭПР самолета „Конвер 990м от азимутального угла наблюдения при вертикальной поляризации падающего поля:

/ — э к с п е р и м е н т а л ь н а я к р и в а я ; і?

г р а н и ц ы р а с с ч и т а н н о г о и н т е р в а л а .

Рис. 38. Зависимость ЭПР самолета „Конвер 990“ от азимутального угла наблюдения при горизонтальной поляризации падающего поля:

1— э к с п е р и м е н т а л ь н а я к р и в а я ; 2 - г р а н и ц ы р а с с ч и т а н н о г о и н т е р в а л а .

Соседние файлы в папке книги из ГПНТБ