Файл: Игнатов, В. А. Статистическая оптимизация качества функционирования электронных систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 19.10.2024

Просмотров: 132

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Распределение нормированной частоты

при /= 0 близко к дель­

та-функции, что отражает условие точной

настройки генератора

в начальный момент времени на номинальную частоту—параметры ФЦ в точности соответствуют номинальным. В случае необходимо­ сти введением соответствующих начальных вероятностей нетрудно описать наличие разрегулировки генератора в начальный момент времени. Графический анализ показывает, что с течением времени

Рис. 6-2. Графики зависимо­

Рис. 6-3. Графики зави­

стей от и а частоты от вре­

симостей коэффициентов

мени.

асимметрии

и эксцесса

 

частоты от

времени.

Рис. 6-4. Вид плотностей вероятностей нормирован­ ной частоты при ^=0 (/); 326 (2); 653 г (3).

210

математическое

ожидание, коэффициенты

асимметрии и эксцесса

нормированной

частоты Fa= (RHCK)~l

генератора

убывают, а дис­

персия возрастает.

 

 

 

 

При регулировке йС-генератора, как правило, изменяют вели­

чину лишь одного из сопротивлений, поэтому полагая среднее

вре­

мя регулировки

/иРЕ=0,3 ч и учитывая,

что <?е = 1,

получим т р=

= т рЛ= 0,3 ч, следовательно,

 

 

 

 

 

'WO

 

 

 

 

 

kr = ggg g= 0,999541.

 

 

Таким образом, располагая лишь

средними значения /Пін,

т1С

и іП\R трех Э'СВ, мы определили все

основные вероятностные

ха­

рактеристики качества ÄC-генератора. В этом одно из главных до­ стоинств вероятностных расчетов с помощью ЭСВ.

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

6-4. Оптимизация синтеза качества ^восстанавливаемых элементов

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

ляемых

переменных — математическое

ожидание,

дисперсия

на­

чального

значения и

стабильность выходного параметра; в

роли

ограничений — основные изменяющиеся при оптимизации

слагаемые

приведенных часовых

расходов (ПЧР):

оптовая

цена

и ущерб,

обусловленный отказом устройства во время использования по наз­ начению. Экономические показатели, необходимые для оптимизации, рассчитывают методами экономического анализа производства и эксплуатации изделий, излагаемыми, например, в [Л. 12].

Начнем с рассмотрения простейшей задачи оптимизации ста­ бильности выходного параметра изделия при фиксированном мате­ матическом ожидании и дисперсии

 

2

 

 

 

 

2

 

 

 

Щ = тх (0) =

2

х *і Рі ,

4 =

(0) =

Ц

(x*t )2Pi — т0,2

(6-35)

 

;=о

 

 

 

 

і=о

 

 

где, как обычно,

х * і

— і-ое

квантованное

значение параметра;

р, —

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

в

і-м

интервале

квантования.

Выберем в роли целевой функции ПЧР. Для невосстанавлива­

емых изделий приведенные часовые расходы

 

 

 

 

 

П = (2+ ГдГа) Т

1o + ElZ.

 

 

 

Применим простейшую

степенную

аппроксимацию

и получим

г яалитическое выражение для оптовой

цены

в виде г = са~чі,

тогда

211


с учетом того, что Та=Ах(2ро+рР)а~1и для математического ожи­ дания ПЧР получим:

П = (cm\-q+ гАТАпи) [4х (2р0+ > ,)] ~1 +

е,спцч.

(6-36)

Дифференцирование уравнения (6-36) по mi и приравнивание

нулю производной дает уравнение

 

 

гАТА,пЧі + і + С(1 — q) Щ qeiChX (2р 0+

p t) = 0,

(6-37)

решение которого определяет /и1оп*. Характеристики экономических ЭСВ с и q, которые определяют зависимость оптовой цены невосстанавливаемого элемента от стабильности его выходного параметра, находят путем экономического прогнозирования [Л. 12].

Пример 6-5. Рассмотрим простейший невосстанавливаемый эле­

мент

типа микромодуля.

Пусть

Ах — 0,1;

ро ~ 0,9; рі =

0,1; Гд =

=200

руб -ч -1; ГА= 2 ч;

еі = 2,28-10-5

ч~1;

т с= 2 • 10-6;

mq = l.

При 9 = 1 уравнение

(6-37)

имеет

простое решение

 

УМ * (2р0 + р,) =* 1,62-10-’ ч -1;

 

 

V

ГА Т А

 

 

 

 

 

 

 

 

 

 

 

 

 

7'ооп*=1,17-10е

ч;

г 0пт = 12,3 руб.;

ПМин—6,82-10- 4

руб-(ч-изде­

л и е ) =5,97 руб (год • изделие)-1.

 

устройств

mq = 2 и

Для

более

сложных невосстанавливаемых

в рассматриваемом

случае

решение

(6-37) гпІ0ПТ«0,973 • ІО- 2

ч~1,

ГОп т»1 950 ч,

2 0nT« 2 1 1

руб., Пмин=Э,318 руб-(ч-изделие)-1=

=278 руб (год ■изделие).

 

 

 

 

 

 

Таким образом, уравнение (6-37)

позволяет

определять

опти­

мальную

стабильность параметров невосстанавливаемых

элементов

с учетом

амортизационных

отчислений, эффективности

капитальных

вложений в данной отрасли и ущерба, обусловленного отказами. Перейдем к рассмотрению более сложной задачи. Найдем опти­

мальные математическое ожидание и дисперсию а0, а также мате­ матическое ожидание' ßi. При выбранных квантованных значениях

вероятности

рі однозначно определяют

моменты

а0 [см. формулу

(6-35)], эти

вероятности

удобно рассматривать

как управляемые

переменные. Тогда целевая функция

 

 

 

 

 

7’0

Ах (2ра+ рі) шД1

 

 

(6-38)

включает три управляемых переменных, а ограничения

 

2

 

 

 

 

 

 

 

Т і

p t == I > 2 cam[

-f c0p*>+

схр ^

+

сгр ^ яК

(6-39)

і=0

 

 

 

 

 

 

 

Записывая вспомогательную

функцию Лагранжа

 

 

(Ра >Р \і P it

j

^-2)&X (2p0

P \ )

 

( 6 -4 0 )

212


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

2 Д

— Л2 — Âsc0<70jc>g0

1 =

0 ; Дх п ц 1— +

^гсі9 і/’Г (,1 + 1> =

— X t +

Л2с2<7 г ^ (?,+ 1) =

0;

 

 

■— (2р0+ р г) Дx m f2 + X2caq1Щ '<’а+Х) = 0;

 

2

 

 

 

 

2 Pt =-■1; c0« 7 ?a + c0/J?3 +

CjpJ-'b + с .р ^ » =

г.

i=0

 

 

 

 

(6-41)

Эту систему лучше всего решать с помощью стандартных под­ программ ЦВМ.

Пример 6-6. Рассмотрим относительно сложное невосстанавли-

ваемое устройство. Пусть іДх=20%,

математическое ожидание са=

= Co= Ci = ca= l ;

qa = q 0 = q l =

q 2=:\-t

z = 200

руб.; **о = 90%;

x*i =

= 70%; **2=50%.

 

 

 

 

 

 

Оптимальное решение ро=0,839;

p i=0,094;

р2=0,066; /Яі»5,37Х

Х І0- 3 % • ч~1;

Лі=74,1; Л2=0,355; m0^ 85,5%;

а20 « ( П ,7 % ) 2

найде­

но с максимальной погрешностью 5%

 

 

 

 

 

20 (2 -0 ,8 3 9 9 + 0,0937)

: 6,63-10* ч.

 

‘ омакс ' :

5,37 -10_3

 

 

 

 

 

 

Коэффициент Л2 =

— дг К

показывает,

что увеличение

макси­

мального среднего времени безотказной работы на один час влечет повышение оптовой цены на 2,82 руб. Множитель Лі показывает, что ошибка в оценке Г0Мако на 1 ч возникает тогда, когда начальные вероятности определены с погрешностью примерно 1,35%.

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

2

2 Рі = К Т 01--= Хх {2р0 + Рі ) т ~ 1.

(6-42)

і=о

Оптимальные параметры определяет решение системы

со9оРо0_1 = + 2ДхХ2т^и, c1qtp\~4l+i'>+ Хг — Дх \2т у 1 = 0;

Са<72Р р Чг+1) + Л, = 0; caqam l (?а+1) = (2р 0 + p t ) Л2Дх т ~ 2;

2

2 А = 1; (2л + Рі)Алда7"'= 7'01.

і=о

(6-43)

15—385

213


При параболической аппроксимации ухудшения определяющего параметра

2р0(X , — д 0)

,

-J/raf —4a2(*і-

•«„)

2 (/70 + /7,) Ах

(6-44)

 

| / Я] — 4аг (х2 — а0) + ] / я? — 4д2 (*, — а»)

Соотношение (6*44) позволяет оптимизировать не только то, а 2о» т ь но и т 2, о22, сг2ь а также и поле допуска [то, х2]. Роль ограни­ чений играют

2

S Pt ==■ 1 . 2 =

Са\ПНч' -+- Ca2n q qi + catsmq° +

 

( = 0

 

 

+

c oPo + С , Р Г Ѵі + c i P 2 4 -

( 6 ‘ 4 5 )

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

Как правило, из опыта проектирования и изготовления нетрудно ориентировочно оценить расходы, связанные с повышением ста­ бильности параметров элементов, с улучшением других технических показателей и т. и. (Л. 12]. В этой связи важное значение имеют развиваемые в настоящее время методы технико-экономического прогнозирования. Так как на ранних этапах проектирования исход­ ные данные оптимизации и сама оптимизация носят приближенный характер, во многих случаях целесообразно применять линейную аппроксимацию ухудшения определяющего параметра— это сущест­ венно упрощает вычислительные процедуры.

6-5. Оптимизация синтеза качества обслуживаемых систем

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

214


мальное качество изделий. Обратная задача ОС в принципе является многоцелевой. Целевыми функциями могут быть оптовая цена, капи­ тальные вложения, расходы на ТО, ущерб от простоя из-за отказов и т. п. Классический подход к ее решению заключается в использо­ вании принципа Парето [Л.- 59]. Более удобным является применение единой целевой функции (1-3) или ее отдельных слагаемых.

В этом параграфе мы рассмотрим задачу комплексной оптимиза­ ции стабильности выходного параметра, ремонтопригодности и экс­ плуатационной технологичности. В роли целевой функции выберем коэффициент технического использования Кт (6-18); в роли управ­ ляемых переменных — т р, т п, /щ, ѵ; в роли ограничения — ПГР.

При изменении т%, т р, ѵ и тп изменяются оптовая цена z изде­ лия, расходы на ПО (спо) и АР ( c a p ) и стоимость ЗИП (сзп). По­

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

Оптовая цена обслуживаемого изделия

 

 

 

г =

ср ѣ ^ 1-f- сгу~Чг+ с3т~Яі + сріц 4*,

руб (изделие)-1,

(6-46)

где Сі

и qt — экономические ЭСВ. Предполагая

для простоты все ПО

однотипными, выразим число N ПО в^год

через ѵ и /«„, тогда

7ѴП0=

= c5v?s -j- ty n j, профилактик-год-1. Учтем,

что средние расходы

на одно ПО складываются из стоимости' деталей гп1, основной и дополнительной заработной платы (гп2), приходящейся на одно ПО (если за изделием не закреплен постоянный обслуживающий персо­ нал), и стоимости ущерба из-за простоя гпз, тп, тогда

Спо =-- (Пи + Сцг + rnama) (csyqs + cem j) +

 

+

стт™, руб (год-изделие)-1,

(6-47)

где Съ, се, ст, <?5, <7в,

<7т — экономические ЭСВ.

с учетом

Так как среднее

число отказов в год 7ѴАр = 7’гГ_1о, то

того, что Тъ=Ах(2ро+рі)ті,

а ТА = Щ, получим:

 

 

 

 

 

Ттш1

 

С А Р = ('Т + '-а +

'ѴИр) А х { 2 р 0 + Р і ) +

 

— Яѵ

 

— Я

п,

руб (год-изделие)-1.

(6-48)

+ cv V

+ cama

Стоимость ЗИП

 

 

 

 

 

Сзп =

с ?ѵ?7 +

csmi8> руб (изделие)-1.

(6-49)

Таким образом, мы получили аналитические выражения для тех составляющих ПГР, которые существенно изменяются при измене­ нии тр mp, /Лп п V. В общем случае необходимо учитывать и изме­ нение расходов на контроль работоспособности, на вспомогательные материалы, изменение стоимости потребляемой энергии и т. д.

Если ограничения заданы раздельно на г, Спо, Сдр и Сзп, то оптимальное решение определяет система из восьми нелинейных

15*

2 15