
книги из ГПНТБ / Штагер, Е. А. Рассеяние волн на телах сложной формы
.pdfиых элементов и представлены в виде таблицы чисел:
аі(Т/); 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 - г р а н и ц ы р а с с ч и т а н н о г о и н т е р в а л а .