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

книги из ГПНТБ / Дракин, И. И. Основы проектирования беспилотных летательных аппаратов с учетом экономической эффективности

.pdf
Скачиваний:
39
Добавлен:
21.10.2023
Размер:
8.98 Mб
Скачать

параметров процессов и конструкций. Этот метод использует не только способы уменьшения рассеивания параметра, но и актив­ ную компенсацию отрицательного влияния рассеивания.

Идея этого метода заключается в следующем. Рассеивание любых параметров БЛА влияет в конечном счете на те или иные баллистические характеристики (параметры, связанные непо­ средственно с ошибками наведения или эффективностью боевой части, здесь не рассматриваются). Но баллистические характе­ ристики ЛА определяются при заданном законе баллистического коэффициента величиной относительного веса топлива и законом его расхода (законом регулирования тяги). Следовательно, уве­ личивая запас топлива, можно компенсировать отрицательное влияние рассеивания.

Если в предыдущей постановке задачи (см. рис. 5. 3) в каче­ стве критерия необходимо было принимать затраты на выполне­ ние целевой задачи с заданной вероятностью, то в последней постановке задачи можно ограничиться для критерия стоимостью Qo одного БЛА. Действительно, в этом случае параметрическая надежность по рассматриваемому параметру для всего БЛА бу­ дет доводиться до і?£=1, поэтому потребное количество БЛА для выполнения целевой задачи не меняется.

На этапе проектирования целесообразно определять не опти­ мальную параметрическую надежность, а оптимальные допуски на параметры. Следует заметить, что оптимальная параметриче­ ская надежность в конечном счете сводится к определению вели­ чины допусков.

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

В целях контроля допусков обычно интересуются предельны­ ми отклонениями параметров, в данном случае — весом конструк­ ции, который обозначим AnpGK. Схема определения оптимальных допусков следующая. Пусть максимальное рассеивание веса де­ тали или агрегата

Дтах°к~За(Ок).

Ущерб от этого рассеивания в пределах AnpGK можно компенси­ ровать увеличением запаса топлива. Пользуясь формулой (2. 17), находим увеличение полетного веса БЛА, вызванное увеличением веса конструкции на AnpGK,

Д ^ О ~ f f Д пр^к •

При этом баллистические характеристики БЛА останутся преж­ ними. Оставшееся некомпенсированным рассеивание (AmaxGK — AnpGK) можно устранить путем отбраковки или доработки конст­ рукции по весу.

200

Очевидно, что допуск AnpGK будет оптимальным в том случае, если сумма ущербов от увеличения веса А G0 и от браковки или доработки конструкции будет минимальна. На рис. 5. 4 показан характер зависимостей следующих величин: AQKOM — затраты на компенсацию предельного допуска AnpGK, вызванные увеличени­ ем полетного веса; AQ6p — экономический ущерб от браковки ана­ лизируемых изделий, для которых A6K>AirpGK; А(2дор— затра­

ты на весовую доработку изделия до величины AnpGK; (ДОк)™? — оптимальный допуск на вес, в случае доработки отбракованных

изделий; (Д(?к)опт — оптимальный допуск на вес в случае браков­ ки изделий с AGK>AnpGK.

Для определения AQK0M нахо­ дим увеличение веса БЛА при увеличении веса GK анализируе­ мого изделия (тот или иной кон­ структивный элемент) на вели­ чину AnpGK. Согласно формуле (2.17)

G K

Дпр^к

Д

Go

^рДпр^К-

 

Здесь коэффициент роста следует определять с учетом изменения веса как первой, так и второй сту­ пеней, если анализируемая кон­ струкция находится на 2-й сту­ пени. Аналогично следует опре­ делять /р для трехступенчатого БЛА.

Зная относительные весовые характеристики первой и второй ступеней, определяем изменения весов. Для случая нахождения анализируемой конструкции на второй ступени получим

Д02 — /ргДпр^к,

Дб?і—(/pi

1)Д02,

здесь AGi — изменение веса

энергоблока

(ускорителя) первой

ступени.

 

 

Для проведения экономических оценок разбиваем общие из­ менения весов на изменения веса конструкции и веса топлива.

В связи с тем, что величина AnpGK небольшая

(единицы процен*

тов), то (см. разд.\ П

1 гл. I)f*g2 +

а 2Мт2

f Д Q

 

(5. 37)

 

 

/м4"'

■ '

 

w

+ <! + « !) «

Д 0 т2= P g 2

+ (1 + а 2) Мт-2

/ р г Д іц А о

(5.38)

 

 

I", .

 

201

ДОкТ =

HÄT+ аіЦт1

■ (/,pl'

1) f ргДцр^к»

(5. 39)

д

Q.

Bgn +

(1 +

«О Н-ті

 

(5.40)

ТІ~

И.?1 +

(1 +

a l) Hrl

-(/pi ^ /ргДігр^к-

Зная удельные стоимости конструкции qK, емкостей qa и топ­ лива qT первой (энергоблока) и второй ступеней, находим по­ требные затраты A Q K O M на компенсацию ущерба от допущения предельного отклонения веса конструкции AnpGK от номинала

AQKOM

V

- g l Q i a

+

а

2 \ > - Л Ч а 2 т2<7т2

/р^ +

 

а іНті9аі,

+

+

(1 +

а2) Нт2

^

(5.41)

! !У§ч9к1 + .

Нті9т1

/ ^

^

 

ДцрОк •

+

(1 + а і) Цгі

(-/P I

 

 

Р2

I

 

ч

 

 

 

 

 

 

Как видно из формулы

(5. 41), величина AQK0M линейно зависит

ОТ ВеЛИЧИНЫ AnpGK.

Определяем ущерб от забраковки всех анализируемых дета­ лей, имеющих отклонение веса от номинального больше AnpGK. Пусть общее количество деталей будет «д, а количество забра­ кованных деталей Пбр. тогда относительное количество забрако­

ванных деталей

(5.42)

£ = —«л .

Если стоимость одной детали <3Д, то распределяя стоимость забракованных деталей на оставшиеся пл — «бр, получим увели­ чение стоимости одной детали

А<2бр

С?длбр

k

(5.43)

фд

фд (/2д ^бр)

^ ^

Рассмотрим распределение плотности вероятности <p(GK) по величине GK. Для случайных отклонений AGK от номинального значения G„ можно принять нормальный закон распределения плотности вероятности (рис. 5. 5). Интересующая нас величина k будет равна заштрихованной площади, которая при нормаль­ ном законе распределения плотности вероятности равна (см., на­ пример, [26]) :

k =

—^=r

Г

dv,

 

у/

 

 

где

2rt

J

 

 

 

Дпр^К

 

(5.44) ’

 

а (Д О к)

 

здесь o(AGK) — среднеквадратичное значение величины

АGK*

устанавливаемое по статистике аналогичных образцов.

 

202

Используя функцию Лапласа

' С

—1

 

у 1

 

Ф(г) = — —

о е

dv,

(5.45)

у2п,

J

 

 

можно написать

& = ф(оо) — Ф(г).

Так как

Ф (оо) =0,5,

(см. [26]), то

 = 0,5 — Ф(г).

(5.46)

Применяем в качестве критерия оптимальности

величины

AnpGK суммарный ущерб

 

(5.47)

QS= AQK0M4 -AQ6p.

Оптимальное значение AnpGK находим из уравнения

 

dQz

0.

(5.48)

^ОпрОк)

Уравнение (5. 48) является трансцендентным и в конечном виде аналитически не решается. Поэтому целесообразно его ре­ шать графически. Еще более целесообразно графически строить функцию (5. 47) и определять минимум по графику. В этом слу­ чае можно обоснованно сделать отступление от формального оптимума.

Как показывает опыт численных расчетов, величина k, вхо­ дящая в выражение (5.43), невелика: тысячные или сотые доли

203

единицы. Поэтому без

большой

погрешности можно

принять

в выражении (5. 43)

ДОбр

48)

можно разрешить

(5.49)

В этом случае уравнение

(5.

в конеч­

ном виде до конца. Для упрощения выкладок обозначим:

 

-X

ДпрОк,

 

(5. 50)

 

A

Q K O M

 

(5.51)

 

q " P

A„pGK

*

 

 

 

о==з(дGJ.

 

(5.52)

Выражение (5. 47) при этих обозначениях будет

Qs ~j-fcQд.

Дифференцируем по х в соответствии с уравнением (5.48),

 

^пр + Фд

dk

0 .

 

 

dx =

 

Учитывая значение

величины

k

по

формулам (5.46), (5.45) и

(5. 44), получаем

 

Qx

 

 

~2

°ä= o .

 

 

?„р----- — е

 

 

1 Д

 

Отсюда, логарифмируя,

- / 2зт

 

 

 

 

 

 

 

(5. 53)

 

-^011'

Г2 ln

 

 

Т

 

 

 

L

 

т/2я9пра J

 

Формулу (5. 53) целесообразно применять в тех случаях, ког­ да деталь или агрегат нельзя доработать по весу, или когда этонецелесообразно, например, при массовом производстве. В боль­ шинстве случаев, вместо браковки, целесообразно доработать агрегат до приемлемого веса.

Доработка по весу каркасных агрегатов является сравнитель­ но несложной операцией и сводится к фрезерованию отдельных мест. Стоимость доработки можно представить в виде

Дфдор—BQx-

При браковке е=1. Доработка же, обычно не будет превышать 10% стоимости агрегата или детали, т. е. е<70,1.

Величина е главным образом зависит от типа агрегата, ос­ новного конструкционной материала и технологии его произ­ водства. Величина е слабо зависит от величины снимаемого веса, так как большая часть затрат при доработке идет на установку

204

доводимого агрегата или детали, на установку и регулировку инструмента, на замеры, взвешивания и т. д.

Например, согласно приводимым в работе [62] графикам по затратам на снятие материала фрезерованием, стоимость сня­ тия дюралюминия за одну установку детали и инструмента со­ ставляет: для снятия 0,5 г материала 6 коп., 1 г — 7 коп., 2 г — 8 коп. Как видим, вес снимаемого материала увеличился в 4 ра­ за, а затраты на снятие увеличились только на 30%.

Как нетрудно видеть из вывода формулы (5.53), для случая

доработки оптимальная величина

Е<?дПри

 

 

 

 

 

 

 

 

 

-^опт

2 In-

 

 

 

 

 

 

(5.54)

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/2я<7

 

 

 

 

 

 

 

 

Эта формула включает и формулу (5.53), которая соответст­

вует е=1.

 

 

 

 

 

 

 

 

 

 

0 , 0 2

вес

 

qKp = qa.=

Пример.

Определить оптимальный предельный допуск на

крыла одно­

ступенчатого БЛА при следующих

условиях: сг(ДО„р) =

 

 

GKp,

 

= q„ = 50 qT,

pg=Ü,12, a = 0,12,

 

pT = 0,5, fp = 2,5,

e= 0,08.

 

 

 

 

 

По формулам (5.41) и (5.51)

 

 

 

 

0

 

 

 

 

 

 

 

 

0,12 + .0,12-0,5 + 0,5-0,02

ft =

0 ,7

-qk.

 

 

 

^np — 2,5

0,12 +

1,12-0,5

 

 

 

 

 

 

 

По формуле (5. 54)

 

 

______ 0,08Qfl______

 

 

 

 

 

 

 

Дпр^кр

 

 

 

 

 

 

 

 

 

−−−−−−−− − −−−−−−−−−−−= 0,02

21n

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что

GKp

 

К 2я ■0 ,

7q$

,02GKp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q n = G KpqK,

 

 

 

 

 

 

 

 

 

 

 

 

 

^пр^кр

0,0256.

GKp

Относительное количество дорабатываемых крыльев

6 = 0,5 — Ф(г),

где по формуле (5. 44)

0,0256

0,02 = 1,28.

Функция Лапласа Ф (г) по таблицам (см., например, [26])

Следовательно,

 

Ф( 1,28) =0,4.

 

 

6=0,1.

Дополнительное топливо по формуле (5. 38)

Увеличение полетного веса

Д

G r

= 0,0471.

кр

G

 

 

 

 

- ^ р° кр- = 0,064.

ДG0

G

кр =

кр

 

 

 

1

205

5.ОПТИМАЛЬНАЯ НАДЕЖНОСТЬ СИЛОВЫХ КОНСТРУКЦИЙ

5.1.Детерминированный и статистический методы обеспечения

прочности силовых конструкций

Потребную прочность и устойчивость силовых конструкций ЛА обычно обеспечивают, исходя из детерминированной несу­ щей способности конструкции, за которую принимается разру­ шающая нагрузка. Разрушающая нагрузка іѴр определяется из максимальной эксплуатационной нагрузки Дуэтах, при этом

■ ^p= f N aгааХ)

где / — коэффициент безопасности.

Для самолетов коэффициент безопасности берется в интер­ вале 1,5—2,0, для БЛА 1,2— 1,3 [19].

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

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

Для БЛА одноразового применения картина получается иная. Действительно, накопления остаточных деформаций практически нет, так как повторение максимальных эксплуатационных нагру­ зок, или не имеет место, или эти повторения исчисляются едини­ цами раз. Заметное же снижение разрушающих напряжений возникает при числе циклов нагрузки, исчисляемом сотнями. Поэтому для БЛА можно допускать в полете напряжения, не­ сколько превосходящие предел текучести. Практическое отсут­ ствие накопления деформаций у БЛА повело к снижению коэф­ фициента безопасности для БЛА, по сравнению с самолетами.

Еще более существенным фактором для прочности БЛА яв­ ляется то, что для БЛА допустимы разрушения в полете, что не допустимо для самолетов. Действительно, для БЛА допускается надежность меньше единицы. Например, при тренировочных запусках в ВВС СШ А многих тысяч управляемых ракет средняя их надежность составляла 0,9 [74]. Большинство отказов БЛА бывает вследствие отказов системы управления и двигательной установки. Очевидно, что если небольшое снижение надежности конструкции планера или двигательной установки, например,

206

в 10—1Ö0 раз меньше общего снижения надежности, поведет к снижению веса БЛА и соответственно его стоимости, то такое снижение надежности может быть целесообразным.

Для иллюстрации высказанной мысли на рис. 5. 6 приведен график зависимости надежности R от коэффициента безопасно­ сти /. График построен для случая среднеквадратичного суммар­ ного рассеивания эксплуатационной нагрузки и несущей способ­ ности конструкции, отнесенной к математическому ожиданию

эксплуатационной

нагрузки

ан. э=0,2. Предполагается

нормаль­

ный

закон

рассеивания.

 

 

 

 

 

Одна

кривая

соответствует

R

 

 

 

 

случаю

принятия

для пре­

0,99

 

 

 

 

дельно

допустимых

напря­

бр

/

Ор

 

 

жений

 

Од значений,

равных

0,98

 

-"

 

 

ал/ар =

0,97

 

 

разрушающему

напряже­

 

 

 

 

нию сгр. Другая кривая соот­

 

 

 

 

 

ветствует

значению

 

 

 

Он.э 0,2

 

 

= 0,8.

случае

<тд/сгр = 1, оче­

 

 

 

 

 

В

О М

 

 

 

 

видно,

 

что без

заметного

0,95

 

 

 

 

ущерба

 

для

надежности

 

 

 

 

 

коэффициент

безопасности

0,9b

 

 

 

 

может быть принят

равным

0,6 0,7

0,8

0,9 1,0

1,1

/

/«0,9, так как при этом на­

 

 

 

 

дежность

очень

.близка

к

 

Рис. 5. 6.

 

 

единице

(Д « 0,9999). В слу­

 

 

 

чае ад/ор = 0,8,

при

той

же

 

 

 

 

 

примерно

надежности мож­

 

 

 

 

 

но принять /«1,1. Эти коэффициенты безопасности значительно ниже обычно принимаемых для БЛА /= 1,2—1,3 [19].

Заметим, что принятое на рис. 5.6 значение ан.э = 0,2 не яв­ ляется большим, так как относительное среднеквадратичное зна­ чение рассеивания только перегрузки при выполнении целевой задачи на разных высотах может иметь такое значение (см., на­ пример, рассеивание перегрузок у самолетов [19]). Увеличение

од. э поведет к уменьшению коэффициента безопасности. Стохастический подход к определению коэффициента безо­

пасности, проведенный по рис. 5. 6, базируется на том, что мак­ симальные нагрузки встречаются очень редко, и поэтому приня­ тие их за основу формирования прочности конструкция БЛА яв­ ляется нецелесообразным.

Допуская отказ БЛА в полете вследствие разрушения или потери устойчивости конструкции, проектировочную задачу мож­ но поставить так: обеспечить такую прочность и устойчивость конструкции БЛА, при которых соответствующие вес и надеж­ ность конструкции приводили бы к минимальным затратам на выполнение целевой задачи. При этом учитывается, что при сни-

207

женин надежности конструкции вес и ее стоимость уменьшается, а количество БЛА для выполнения целевой задачи увеличи­ вается.

Критерием оптимальности надежности при указанной поста­ новке задачи должны являться затраты на выполнение целевой задачи Qq.3. Уравнение оптимизации в этом случае будет (5.1), следствием которого является уравнение (5.3).

5. 2. Методика оптимизации надежности силовых конструкций по экономической эффективности

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

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

Левую часть уравнения (5. 3') можно представить в виде

dQ B

__

dQa d G

д ___

dQ„

dG0 dG^

dR R

дС?д

d R

Ö

G

Q

dG^

dR Ä

 

 

;l

 

 

 

 

 

 

 

 

 

 

 

Здесь учтено, что надежность влияет на стоимость вылета толь­ ко через вес детали, то есть такие операции как отбраковка по надежности отсутствуют, то есть д ^ в/дЯя = 0.

Первый множитель правой части найденного выражения мож­ но представить в виде

 

 

dQB

 

dQg

I

dQa

ßQT

pQ3

 

 

 

dG0

 

ÖGQ

 

Ö GQ

 

G

 

 

 

 

 

 

 

Ö

 

Q

ÖG Q

 

где

Qg

— стоимость конструкций БЛА, вес которых пропорцио­

нален полетному весу;

Q

a — стоимость

конструкций,

вес кото­

рых

пропорционален

весу

топлива;

 

QT — стоимость

топлива;

Q0 — стоимость эксплуатации БЛА,

включая затраты на соору­

жения и оборудование стартовой позиции, приходящиеся на один

БЛА. Для авиационных БЛА

сюда входит стоимость носителя

и его эксплуатация.

 

 

 

Представляем стоимости конструкций и топлива в виде (2. 1)

Qg

--

Q*

-- a<*Gaa,

Q

T-- аТ^дТ, Qg-- QgO-pÖgG^,

 

 

 

здесь Qao — эксплуатационные затраты, не зависящие от G0.

208

Учитывая, что

 

Ga=aw TO0,

GT= p TG0,

ö„=!J-„O0,

Я*

Q g

я*-- G . ■ .

Яг

 

Qr

Gr

 

Gr

 

n Qa Q-эо ___

 

V—1

>

получаем

Чэ

“'3lJo3

 

 

 

vTqr)

 

3

ÖQB

Ѵ Л

+ (M .a +

 

 

(Іт + ѴЭ<? .

dGn

 

 

Здесь, в связи с небольшим влиянием на результат, не учтены весовой и аэродинамический масштабные эффекты.

Согласно формуле (2. 17)dG0 _

,

 

 

 

 

дал

 

 

 

где /р — коэффициент

 

Ур’

 

уравнение

(5. 3)

роста. Следовательно

можно теперь представить в виде

ѵ э ^ э]

=

.

(5. 55)

/ р М Л + М «

а + ѵ т<7т) I S +

 

К

В этом уравнении

величины vg и qg — средние значения для

всех конструкций БЛА, вес которых пропорционален полетному весу. Величина QB равна затратам на вылет с учетом эксплуата­ ционных затрат, приходящихся на один БЛА. Следует заметить, что, вследствие сравнительно слабого влияния Q B на оптималь­ ную надежность конструкции, эксплуатационные затраты могут учитываться приближенно или ориентировочно.

В уравнение (5.55) коэффициент безопасности непосредст­ венно не входит. Однако вес оптимизируемой детали или агре­ гата зависит от коэффициента безопасности. Эту зависимость

с приемлемой точностью можно представить в виде

(5.56)

Оя= А + В / .

 

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

Справедливость формулы (5.56) для несущих поверхностей следует из многочисленных весовых формул для крыльев. См., например, [69], [70], а также формулу (3. 57).

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

209

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