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

книги из ГПНТБ / Игнатов, В. А. Статистическая оптимизация качества функционирования электронных систем

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

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

при /= 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

у р а в н е н и й , п о р о ж д а е м а я в с п о м о г а т е л ь н о й ф у н к ц и е й Л а г р а н ж а

. ,

 

mp, тп, V ,

А,, А2, А3,

А4)

(2/Й! +

ѵДх)

L (m it

= д7-------------------

 

 

 

 

 

 

Л ' ( / И і , V , тѵ,тп)

К

(С11% 42 +

с гѵ9а + с 3дар 9з +

c4m j

ъ _

г*) _

■А2 І(га1 + гп2 +

гпз/кп) (c5v?s +

c6m j) +

сшт\т -

С*по]

 

 

— А» ( r 1 + r 2 + r 3m p )

 

ТѴлг,

 

 

 

 

, д х ( 2 л + л )

+

 

+ c v V

9 ѵ +

с п и и

D —

С * AP -А* (<ѵ9г +

С./Я?* -

с*зп), (6-50)

где Кт —

Ах (2/я, -f- ѵДх)

 

 

 

 

 

 

Ѵ; mjb mп)

 

 

 

 

 

Если ограничены ПГР в целом, необходимо выбрать суммарную норму а=Сі + а2 амортизационных отчислений на возобновление и капитальный ремонт а= (0,05ч-0,2), год-1, учесть, что норма эффек­ тивности капиталовложений в приборостроительной промышленности е=0,2 год-1, тогда ПГР

Л ~ (&і + а2+ е) (г + Сзп) + Спо + САР, руб (год -изделие)-1.

(6-51)

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

(

,,

 

 

 

 

 

dN

xN (піі,

Щ< v> т п) — Дх (2/Ил + ѵДх)

 

 

 

N2 (mt, mp, т п,

v)

-= А

dz

I dC3I!

1Г1-по _ дС,

öm! "Л d/и,

' дт.1

dm!

 

mf (2от, -j- ѵДх) Дх

+ А

 

dz

ÖCAP

IV2 (/я,, ?йр,

mn, v)

 

 

= o ;

 

 

 

 

Дх2 (2/и, -(- ѵДх) (2/И] 4- ѵ) V

 

Д12 (mlt

mv, ma, v)

 

+

(6-52)

 

 

 

dz

 

дСдо

ÖCAP

 

 

 

= 0 >

 

dm,, л

0/Ип

Дх LIV (m,,

mp, ma, v)

ov (2«j + ѵДх)

 

 

IV2 (Wj,

mp,

w.n, v)

 

 

/ C72Ö2

 

, дСЧ/ЬздТЛ1 .,

d C ~

4--dC

 

 

 

 

 

d\

 

öv

ei (z +

Сзп) +

Cn0 +

CAp —=П*,

 

2 1 6

п о р о ж д а е м а я в с п о м о га т е л ьн о й

ф ун кц ией

Л а г р а н ж а

 

 

 

L (ш ,,

дар ,

тп,

ч, X)

Дх (2да, -j- ѵДх)

 

 

У (да,, дар, да„, ѵ)

 

 

 

 

 

 

 

 

 

 

— А [г, (г + С3„) + Спо +

Сдр — П*],

(6-53)

где

г, = а , - |- а 2-f е;

jV (/«!,

дар,

Гдап,

ѵ) = Дх (2да, +

ѵДх) +

 

 

-f

да^дар +

да,ѵдадДх +

ѵДхдап (да, +

v);

 

 

 

dN

9

 

dN

 

 

 

 

 

 

-дщ=т';

â ^ =Ax(2mi + v)v;

 

dN

 

 

 

 

 

 

dN

 

mPm, + ѵдапДх);

- ^ = Ax (Ах + 2т1тш+

2дапѵ); g ^ - = 2 (Дх +

<?С„

-1

дСп'=0,

ÖCn' ~ СэЯь ('"пі +

ГП2 +

'■дз'Яп) v7s 1;

ÖOT,

СтЯт.т~

,

д/П]?

 

дѵ

 

 

 

 

ÖC

= /-пз (c6v9s + c6m j) +

dm-а

 

 

dCAP

 

dv

(r„, + rn2 + гпзда„) св^да^ 1 ;

(?» +1).

-W

dCAP _

(r, +

r2 4- räWp) T r

^

dCAP =

ггТТтг

 

âmj ~

Аx(2p0-j-pI)

ÖOTp

Дх (2p0+ jp,) '

 

 

dC,

 

 

 

 

 

 

 

 

"AP

-СпЯпта

 

 

 

 

Ода,,

 

 

 

 

 

 

n„—1.

 

 

 

fl,

1.

^^пз

 

 

öCV, <ЭС3

 

 

 

 

--

= 0 .

" ä r =<w *

- ш г ~ с*я^? ’

 

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

Вспомогательная функция Лагранжа

L (да,, дар, X) =

Ах (2да, + ѵДх)

 

\ (с3тр ?3

ciml Qi — z*)

(6-54)

N (да,, дар)

 

 

 

 

 

\ .

 

порождает систему

 

 

 

 

2 [AxN (да,, дар) — Дх (2да, +

ѵАх) (Дх +

2да,дар -f- 2ѵдапДх)]

 

 

Д/г (да,, дар)

 

 

 

 

 

 

С4?4

о;

 

 

 

+ Хда^+'

 

 

 

 

 

 

 

Дх (2да, +

VДх) т\ _

^

c3? 3

(6-55)

 

^ 2(т,,

даР)

*

да^5—1 ’

 

 

 

с3да“ ?3 +

e4m f ?4 = г*.

 

217

Пример

6-7. Пусть Дл;= 0,2;

ѵ=10“ 2 ч~'\

mn= 10 ч.

= 1, г*=2Х

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

ожидания

с3 = 1 ,

Сі =

2;

<73 =

 

ХІО3 руб (изделие)-'.

 

 

 

 

 

 

 

 

Решение (6-55), полученное с максимальной погрешностью 2%,

имеет

вид:

т І0ПТ «

1,00 • 10_3

ч~'\

отр.0пт «

0,49

ч;

X = 0,273 X

ХІО“ 3

руб-'. Максимальное значение коэффициента

технического

использования

Кт макс «0,953.

 

 

 

 

 

 

X,

 

 

X

д%т.макс

показывает,

 

 

 

у

Множитель

К=

— gj-----

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

цены,

например, на Лг=П00

руб. дает увеличение Кт на 0,0273, т. е.

/Ст = /'Ст.макс+ЛЛг=0,9803.

Выберем а = 0,1,

тогда

еі=0,3 год-' и

доля

ПГР,

обусловленная

г,

П1= еі2=600

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

Отсюда следует, что увеличение ПГР на 5%, вызванное ростом z, приводит к увеличению Кт всего лишь на 2 ,8 6 %.

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

При использовании раздельных ограничений на г, Спо, Сар и Сзп возникает практически очень важная и интересная возможность определения наиболее целесообразного дополнительного вложения средств для улучшения оптимальных ТЭХ качества изделий. Зная

 

дКг. ,

дКт

дК.

дКт

 

к'

=- rdz ’ Лз -

дСп0'

К

0СЗП’

 

дСАР

 

выполнив ранжировку этих

множителей и получив, например,

что

А.4>Лі>А,з>Л2,

тогда нетрудно заключить, что в данном случае

до­

полнительные

средства целесообразно

вкладывать

в комплектование

более полного ЗИП. Эта возможность — одно из самых главных преимуществ применения метода Лагранжа.

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

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

Использование вероятностных характеристик качества позволяет прогнозировать оптимальное ТО проектируемых изделий. Отличи­ тельной особенностью является то, что из-за неопределенности исход­ ных данных применяется параметрический анализ — рассчитывается

218

целый диапазон оптимальных режимов ТО. Например, при оптими­ зации периодичности ПО методами, изложенными в § 4-3 и 4-4, приходится рассчитывать топт для целого диапазона изменения интенсивности ухудшения параметра г\=тіАх~1 при вариации т\, Ах, ро, £р и Си- По таким данным можно построить ряд номограмм,

Р0~0,55

0,79 0,91 0,97

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

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

219

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