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

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

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

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

Добавлен: 19.10.2024

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

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

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

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

. ,

 

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


Для иллюстрации рассмотрим пример построения номограммы (рис. 6-5) для определения оптимальной периодичности замены эле­ ментов электронных систем по фактическому состоянию. По оси абсцисс откладывают значения экономического показателя ср(гр—

—сп) -1, характеризующего степень опасности отказа изделия вовре­

мя оперативного использования — чем

больше ср (ср—сп) _1, тем

меньше ущерб,

обусловленный отказом. По оси ординат отклады­

вают значения

нормированного оптимального времени г]Топт, где

г\= Ш іА х ~ 1, mi — скорость ухудшения

параметра, А х — интервал

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

из

Правило

пользования

этой

номограммой

следующее. Исходя

схемного назначения элемента, определяют

его

номинальное

х0

и

критическое

х2 значения,

рассчитывают интервал

прогноза Ах=

= (xo+x2)2~t.

 

 

 

 

 

 

 

 

По

результатам k измерений

(k'— 10 -н- 20)

'оценивают среднюю

 

 

 

 

 

 

k

 

 

 

скорость ухудшения параметра т*х

(xt—xt+1) ( Ш ) - 1 , где At—

 

 

 

 

 

 

I- 1

 

 

 

интервал

времени измерения,

й

хі+і — значения параметра

X,

измеренные в

моменты времени

f,-

и ti+i = ti+At.

Ориентировочно

определяют ущерб ср, обусловленный ремонтом и простоем изделия из-за отказа рассматриваемого элемента, и стоимость с„ ПЗ элемен­ та. Вычисляют ср(ср—Сп) _1 и определяют вероятность р*0 (если р*о неизвестно, выбирают р*0= 1). Используя ср(ср—сп) _1, находят точ­ ку пересечения прямой с абсциссой ср(ср—сп) _1, и кривой, для кото­ рой ро= р*о. Ордината точки пересечения и определяет нормирован­ ную оптимальную периодичность ПЗ. Для определения периодичности в реальном масштабе времени полученную величину необходимо разделить на г|* = т*іАх-1.

Пример 6 -8 . Пользуясь номограммой рис. 6-5, найдем оптималь­ ную периодичность замены электровакуумных приборов в самолетной радиосвязной аппаратуре с учетом того, что ср= (50н-100) руб., а Сп«1 руб. Выходным параметром ламп обычно служит крутизна характеристики или анодный ток. Как правило, допускается ухудше­ ние крутизны на 20—30%, следовательно, Дх=15%. Средняя ско­ рость ухудшения крутизны m * ,« 10 - 4 % ■ч~1, поэтому

^ ~

10 —/о-«-1

= 0,667-

15%

Дисперсия крутизны ламп невелика, если учесть отбраковку, ' проводимую в начальный момент времени, поэтому р*о~0,985.

Рассчитаем величину отношения ср(ср—сп) _і для верхнего и нижнего значения диапазона изменения ср, тогда

Ср .м а к с ( С р . м а к с — С п ) І=Ч00-99 1 ä 5 1 , 0 1 ;

Др.мин(Др.мин—Дд)—~ 50 *49~* ~ 1,02.

22Ѳ