CC BY
88
2008
Доклады БГУИР
№ 6 (36)
УДК 621.382.019.3
ПРОГНОЗИРОВАНИЕ НАДЕЖНОСТИ ИЗДЕЛИЙ ЭЛЕКТРОННОЙ ТЕХНИКИ НА ОСНОВЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ДЕГРАДАЦИИ ФУНКЦИОНАЛЬНОГО ПАРАМЕТРА
СМ. БОРОВИКОВ, А.В. ШАЛАК, А.И. БЕРЕСНЕВИЧ, А.В. ЕМЕЛЬЯНОВ, Е.Н. ШНЕЙДЕРОВ
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 27 декабря 2007
Решалась задача прогнозирования параметрической надежности изделий электронной техники с помощью математической модели деградации функционального параметра. Предлагается получать модель в виде условной плотности распределения рассматриваемого параметра. Это позволяет применить метод группового прогнозирования.
Ключевые слова: изделия электронной техники, прогнозирование, надежность, функциональный параметр, физико-статистическая модель деградации.
Введение
Постепенные отказы изделий электронной техники (ИЭТ) в технической литературе нередко называют деградационными отказами. Критериями их возникновения являются уровни изменения функциональных параметров, при которых, согласно технической документации, изделие считается вышедшим из строя. Скорость протекания деградационных процессов в нормальных (лабораторных) условиях широкого класса ИЭТ незначительна. Однако в современных электронных устройствах общее количество ИЭТ велико и аддитивное влияние изменений функциональных параметров отдельных элементов вследствие медленной деградации может оказать заметное влияние на поведение выходного параметра электронного устройства в целом. При повышенных электрических, температурных, климатических и прочих нагрузках физико-химические процессы деградации ускоряются, и при этом возрастает количество постепенных отказов.
По мере совершенствования технологии изготовления ИЭТ причины возникновения внезапных отказов могут быть в значительной мере устранены. Постепенные отказы, отражающие свойства, внутренне присущие материалам ИЭТ, в частности старение, исключить невозможно. Однако эти отказы и, следовательно, параметрическую надежность ИЭТ можно прогнозировать [1, 2], что позволит выполнять отбор партий ИЭТ с требуемым уровнем надежности. Этим и вызван растущий интерес к постепенным (деградационным) отказам ИЭТ.
Актуальность разработки
При работе ИЭТ его функциональный параметр (обозначим через у) может рассматриваться как функция времени Замечено, что в большинстве случаев параметр у изменяется в нежелательном направлении и в определенный момент времени может достичь предельного уровня, указанного в технической документации или установленного потребителем. Параметрическая надежность характеризует способность ИЭТ сохранять уровень функционального па-32
раметра у(0 в пределах норм (от а до Ь), указанных в технической документации или установленных потребителем, в течение заданного времени ¿З при выбранных режимах и условиях работы. Математическим критерием постепенного отказа в данном случае является нарушение неравенства
а<у t <Ь, (1)
а в качестве количественной меры параметрической надежности может использоваться вероятность выполнения этого условия в течение времени ¿З, т.е.
Рп =Р а< у I 0...*з =Р а < у I О...Г3 . (2)
Для получения достоверного прогноза о постепенном отказе и, следовательно, параметрической надежности ИЭТ надо располагать количественной моделью надежности в виде зависимости деградации функционального параметра ИЭТ от времени, температуры, электрической нагрузки и других эксплуатационных факторов [3]. Такая модель базируется на изучении поведения ИЭТ не только в момент отказа, но и в ходе изменения функционального параметра ИЭТ, т.е. на исследовании кинетики отказов и может быть получена с помощью вероятностно-статистических методов. Построенную таким способом модель деградации функционального параметра будем называть физико-статистической.
/0 Г] и 4 Время г
Рис. 1. Изменение плотности распределения функционального параметра у при работе ИЭТ: ..., ^ — моменты времени (временные сечения); т0, т1, ..., тк — средние значения у во временных сечениях (штриховая линия)
Интерес для практики представляют физико-статистические модели в виде условной плотности распределения параметра у(0. Смысл этих моделей понятен из рисунка 1, на котором величины м(у^) представляют собой плотности распределения параметра у(0 для времени ti (/=0, 1, 2, ..., к), названные временными сечениями. Имея физико-статистическую модель в виде функции м>(у/{), можно решать задачу группового прогнозирования параметрической надежности ИЭТ. Решение задачи сводится к нахождению вероятности того, что для заданного будущего момента времени ^ функциональный параметр выборки ИЭТ будет отвечать условию (1). Согласно работе [2], вероятность Р{утт<у(0<утах, t=t3) может быть найдена по выражению
Утях
Р утт < у(0 < >>тах, Г = Г3 = \м>(у/Г = Г3)с1у. (3)
Записанный интеграл численно представляет собой заштрихованную площадь, ограниченную графиком функции w(y), соответствующей времени t=tЗ, и ординатами у=утт, у=утах, задаваемыми с учетом технической документации или с позиций потребителя ИЭТ (рис. 2).
Рис. 2. К определению вероятности нахождения параметра в заданных пределах
Теоретический анализ
Получение физико-статистической модели деградации функционального параметра ИЭТ облегчается с помощью физического эксперимента, заключающегося в моделировании наиболее типичных условий возникновения механизмов отказов и процессов физико-химической деградации функционального параметра. Поясним некоторые математические аспекты физического моделирования деградации функциональных параметров и прогнозирования параметрической надежности ИЭТ на основе этих моделей.
Во многих случаях количественная характеристика параметрической надежности -Рп(^з), определяемая выражением (2), может быть получена на основе знания закона распределения функционального параметра у(^ в начальный момент времени, например условной плотности распределения w(y/t=0), а также функции изменения у(^ во времени у(0=ф(у0, 0, где ф — символ функциональной связи; у0 — значение у(0 в момент времени t=0. Значение РП(^) является результатом изменения статистического распределения параметра у(0 за время работы
В случае монотонности функции у(0 выборка ИЭТ, характеризуемых параметрами уг(0, в первом приближении сохраняет вид начального распределения в любом временном сечении [4]. Скорость деградационных процессов в ИЭТ при нормальных эксплуатационных условиях невысока, dyi/dt^0 и эффект "переплетания" с течением времени функций у() малозаметен. В этом случае можно говорить о сохранении в течение длительного времени не только вида закона распределения — условной плотности распределения функционального параметра ^(у/0, но и тесной корреляции параметра у(0 для различных временных сечений (рис. 1). Тогда приближенное аналитическое выражение условной плотности распределения для любого
интересующего времени ^ может быть получено путем математических преобразований начального распределения w(y/t=0):
W(y/t = t3) = ц, w(y/t = t3Xt3 , (4)
где у — символ функциональной зависимости.
Физико-химические характеристики деградации функционального параметра у(0, получаемые усреднением по исследуемой выборке ИЭТ, войдут в виде коэффициентов в правую часть выражения (4). Точный вид условной плотности распределения w(y/t=tз) в общем случае будет зависеть не только от вида функций ф(у0, 0 относительно аргумента у0, являющегося источником начального статистического распределения, но и от параметров самой функции ф(у0, 0, которые становятся причинами дополнительных разбросов, со временем накладывающих отпечаток на окончательный вид плотности распределения w(y/t=tз). Нахождение точных аналитических выражений для w(y/t=tз) сопряжено со значительными математическими трудностями. Поэтому идеализация процесса у(0 и упрощения, допускаемые при определении w(y/t=tз), оправдывают себя, поскольку дают возможность хотя бы приближенно определить
количественную характеристику параметрической надежности РП(^) по принятым в теории вероятностей правилам нахождения вероятности вида РП(^)=Р{а<у(0<Ь, t=tЗ}, используя закон распределения случайных величин [5]:
ь
РпЮ = = Р(ЪН3)-Р{аН3), (5)
а
где F(Ь/tЗ) и F(a/tЗ),— значения условной (для времени ^ функции распределения функционального параметра у, подсчитанные для значений у=а и у=Ь.
Разработка метода получения модели деградации
При налаженном технологическом процессе изготовления ИЭТ наиболее часто наблюдается нормальное распределение параметров изделий. Причинами возникновения разброса параметров являются различные флюктуации технологического процесса, неоднородность применяемых материалов, различие в контроле из-за погрешности измерительных приборов, способностей операторов и другие факторы [3].
Возьмем за основу получения модели деградации нормальный закон распределения функционального параметра у. Условная плотность распределения у для рассматриваемого временного сечения t в этом случае
( 1 I У~т(у/ О
НУ'Ч = ,— ^ ехр^ —
,--^Лр >
Л/27га(у/() 2 стО/О
(6)
где m(y/t), o(y/t) — параметры (характеристики) закона распределения.
Величины т(у/0 и o(y/t) представляют собой, соответственно, среднее значение и среднее квадратическое отклонение функционального параметра у во временном сечении t и в неявном виде включают физико-химические характеристики деградации функционального параметра ИЭТ для интересующего времени t. В соответствии с выражением (4) значения т(у/() и а(у/0 определяются как функции времени t и величин m(y/t=0) и о(у^0), являющихся параметрами нормального закона в начальный момент времени ^0):
т у!х = 0 ,ст = 0 ]; (7)
а у!г =фу!1 = 0 ,а у!1 = 0 ], (8)
где ф1, ф2 — операторы функциональных зависимостей, подлежащие определению.
Условная плотность распределения (6), полученная с учетом выражений (7) и (8), может рассматриваться в качестве физико-статистической модели деградации функционального параметра у. Для получения этой модели нужны предварительные исследования выборки интересующего типа ИЭТ. Выборку будем называть обучающей. Ее объем п должен быть не менее 80-100 экземпляров. Процедура построения модели включает следующие основные этапы:
1) измерение значения функционального параметра каждого экземпляра обучающей выборки в начальный момент времени ^0);
2) получение плотности распределения функционального параметра для времени t=0;
3) физическое моделирование деградации функционального параметра каждого экземпляра обучающей выборки в течение интересующего времени
4) нахождение операторов ф1 и ф2 выражений (7) и (8);
5) получение модели деградации в виде условной плотности распределения (6).
Для физического моделирования деградации функционального параметра пригодны ускоренные форсированные испытания ИЭТ, проводимые обычно по типовым методикам. Эти испытания позволяют для каждого экземпляра обучающей выборки получить экспериментальное изменение функционального параметра у в течение интересующего времени
Нахождение операторов ф1 и ф2 выражений (7) и (8) является наиболее трудоемким этапом. Авторами предлагается метод, описанный ниже. Используя результаты физического моде-
лирования деградации функционального параметра экземпляров обучающей выборки, формируют массив данных в виде табл. 1, в которой приняты следующие обозначения: п1, п2, • ••, Щ — номер экземпляра обучающей выборки, отвечающий последнему элементу соответствующей группы; т0®, о0® — среднее значение и среднее квадратическое отклонение функционального параметра для экземпляров/-й группы, /=1, 2, • .., /; ^ — момент времени (временное сечение), для которого контролировались значения функциональных параметров при физическом моделировании их деградации, 7=1, 2, •.., к.
В обозначениях т0®, о0® нижний индекс "нуль" подчеркивает то, что эти характеристики относятся к начальному моменту времени (¿=0).
При решении практических задач число групп / рекомендуется выбирать из условия />(4-5). Количество экземпляров в каждой группе, согласно [2, 5], следует брать не менее 20-30. Число точек контроля (временных сечений) к выбирается с учетом процедуры физического моделирования деградации функциональных параметров экземпляров обучающей выборки и должно отвечать условию к>(4-5).
Табл. 1 можно рассматривать в качестве результатов пассивного факторного эксперимента, в котором факторами являются значения т0&), о0(/) и откликами — интересующие величины т(уИ) и о(у//). При соблюдении указанных рекомендаций число опытов эксперимента N составит 16-25 при объеме обучающей выборки «>(80-120).
Таблица 1. Массив данных для получения операторов ф! и ф2 выражений (7) и (8):
Номер экземпляра обучающей выборки Номер группы Характеристика га(у/г=0) Характеристика а(у/Т=0) Временное сечение
1...П1 1 т0(1) п О ¿1
2 т0(2) п (2) ¿1
¿1
(И/_1+1)-И/ / т0(/) п (/) ¿1
1...п1 1 т О т0 п О ¿2
(И1+1>..И2 2 та2 п (2) а0 ¿2
¿2
(И/-1+1)-И/ / т0(/) п (/) ¿2
1...п1 1 т О т0 п (1) ¿к
(И1+1>..И2 2 та2 п (2) а0 ¿к
¿к
(И/-1+1)-И/ / т0(/) п (/) ¿к
Групповое прогнозирование параметрической надежности
Полученную физико-статистическую модель деградации функционального параметра в виде условной плотности распределения (6) можно использовать на практике для группового прогнозирования параметрической надежности новых выборок ИЭТ исследуемого типа. При таком прогнозировании прогноз получают в виде вероятности того, что функциональный параметр у выборки ИЭТ в течение заданного времени работы не выйдет за пределы норм (а, Ь), указанных в технической документации или установленных потребителем ИЭТ [2]. Эта вероятность, согласно формуле (5) и гипотезе о нормальном законе распределения функционального параметра у во временных сечениях может быть определена по выражению
Ъ-т(уи 3)
-Ф
а-т{у! t3)
Чу/О
(9)
где Ф[-.] — табличная нормальная функция распределения [2, 5], найденная для аргумента, указанного в скобках; т(у/^З) и о(у/^З) — параметры нормального закона распределения, подсчитанные по выражениям (7) и (8) для времени ¿=З.
Оценка достоверности прогнозирования
При решении задач прогнозирования возникает вопрос, в какой степени можно доверять прогнозу, оправдается ли прогноз. Для ответа на этот вопрос пользуются понятием достоверности прогнозирования. О достоверности прогнозирования параметрической надежности на интервале наработок (0. tз) предлагается судить с помощью средней ошибки прогнозирования Лср, определяемой по выражению
Д =
ср
1 к -Т.
к Л2
100 %, (10)
где к — число моментов времени (временных сечений), для которых контролировались значения функциональных параметров при физическом моделировании их деградации; [РП(^)]пр, — прогнозное значение параметрической надежности ИЭТ контрольной выборки, соответствующее временному сечению Ь; [РП(^)]Э — экспериментальная оценка вероятности нахождения параметра в заданных пределах (параметрической надежности ИЭТ контрольной выборки), соответствующая временному сечению ^.
Значения [РП(^)]Э, соответствующие временным сечениям определяют, используя результаты физического моделирования экземпляров контрольной выборки. Пользуются классическим выражением оценки вероятности по частоте события [5], которое применительно к определению [РП(^)]Э принимает вид
кс,)]'='■>. (и,
т
где т(а<у<Ь, t=t) — количество экземпляров контрольной выборки, для которых функциональный параметр у в момент времени ^ находится в пределах указанных норм от а до Ь; т — общее количество экземпляров в контрольной выборке (объем контрольной выборки).
Физико-статистическую модель деградации функционального параметра можно считать пригодной для использования на практике, если средняя ошибка прогнозирования Аср не превышает 5-7%.
Экспериментальные исследования
Предложенный метод получения физико-статистической модели деградации функционального параметра и выполнения группового прогнозирования параметрической надежности ИЭТ был апробирован на биполярных транзисторах типа КТ872А. В качестве функциональных параметров рассматривались И2\Э — статический коэффициент передачи тока базы в схеме с общим эмиттером и £/КЭнас — напряжение насыщения коллектор-эмиттер. Электрические режимы измерения параметров соответствовали требованиям технической документации транзисторов. Для проведения физического моделирования деградации указанных параметров было сформировано две выборки: обучающая объемом «=200 экземпляров, и контрольная объемом 300 экземпляров. Обучающая выборка использовалась для получения моделей деградации функциональных параметров. Контрольная выборка предназначалась для оценки достоверности группового прогнозирования. Применительно к ней для четырех временных сечений t (3840, 8320, 12800 и 17280 ч) решалась задача группового прогнозирования, а на этапе физического моделирования деградации функциональных параметров контролировались их значения в этих временных сечениях. Физическое моделирование деградации состояло в проведении для ИЭТ по типовым методикам [6-8] ускоренных форсированных испытаний, эквивалентных времени 17280 ч с точки зрения функционирования транзисторов в нормальных условиях работы.
При формировании массива данных, аналогичного табл. 1, обучающая выборка разбивалась на пять групп по 40 экземпляров в каждой из них. Выражения вида (7) и (8) строились по данным этой таблицы в виде уравнений множественной линейной регрессии. Для получения уравнений использовалось приложение MS Excel. Вид выражений (7), (8) для параметра Л21Э:
т И21ЭИ = 0,9684да0 + 0,0121ст0 - 0,0092^;
а к2ХЭИ = 0,0134да0 + 0,9083а0 -0,0024^, для параметра иКэнас:
т икэнас /1 = 1,1175да0 + 0,0774ст0 - 0,4118^; а икЭнасИ = -0,0187да0 +1,4812а0 + 0,0704Л
(12)
(13)
(14)
(15)
Величины т(Н21Э/1) и а(Л21Э/0, определяемые по выражениям (12) и (13), являются характеристиками физико-статистической модели деградации параметра Л21Э, а величины т(иКЭнас/0 и о(иКЭнас/0, определяемые по выражениям (14), (15) — характеристиками физико-статистической модели деградации параметра иКЭнас.
В табл. 2 приводятся результаты группового прогнозирования параметрической надежности в виде вероятности выполнения условия (1) для транзисторов контрольной выборки. Приняты следующие обозначения: [Р^)]^, [Рп(^)]э — соответствующее временному сечению прогнозное значение и экспериментальная оценка вероятности нахождения параметра в заданных пределах.
При определении [^П(^г)]пр условие (1) для параметра Н21Э было выбрано в виде Л21Э>18, а для параметра иКЭнас — в виде иКЭнас<660 мВ. Оценки вероятностей вида [РП(^З)]Э подсчиты-вались по выражению (11).
Таблица 2. Результаты прогнозирования параметрической надежности по параметрам Н21Э и иКЭнас транзисторов контрольной выборки
Параметр Вероятность РП(Г,) для времени 11
3840 ч 8320 ч 12800ч 17280ч
[РпЛ)]ио [РП&)]Э [РП(Г,)]ИО [РП&)]Э [РП(Г,)]ИО [РП&)]Э [Рп(Г,)]ИО [Рп&)Ъ
^21Э 0,716 0,719 0,689 0,682 0,656 0,656 0,623 0,631
иКЭнас 0,682 0,678 0,593 0,554 0,506 0,508 0,424 0,423
Средняя ошибка прогнозирования Лср определялась по формуле (10). В случае оценки параметрической надежности ИЭТ по функциональному параметру Н21Э эта ошибка составила 0,55%, в случае использования функционального параметра иКЭнас — 1,96%.
Заключение
Предложенный в работе метод позволяет экспериментально с использованием физического моделирования деградации функционального параметра ИЭТ и статистического анализа результатов моделирования получать физико-статистическую модель деградации этого параметра. Ее получают один раз для рассматриваемого типа ИЭТ, исследуя обучающую выборку. С помощью построенной физико-статистической модели деградации можно в начальный момент времени решать задачу группового прогнозирования применительно к другим выборкам. Прогноз дают в виде определения вероятности того, что функциональный параметр рассматриваемой выборки ИЭТ в течение интересующего времени работы ¿З будет находиться в пределах заданных норм.
Получение физико-статистических моделей деградации и их использование для выполнения группового прогнозирования на примере функциональных параметров биполярных транзисторов типа КТ872А показало оправданность принятия гипотезы о нормальном распределении функциональных параметров ИЭТ во временных сечениях, эффективность метода построения физико-статистических моделей деградации функциональных параметров и высокую достоверность прогнозирования (табл. 2).
THE ELECTRONIC TECHNOLOGY DEVICES RELIABILITY FORECASTING BY THE INSTRUMENTALITY OF THE MATHEMATIC MODEL OF FUNCTIONAL PARAMETERS DEGRADATION
S.M. BARAVIKOU, A.V. SHALAK, A.I. BERASNEVICH, A.V. EMELYANOV, E.N. SHNEIDEROV
Abstract
The problem of forecasting of parametrical reliability of products of electronic techniques by means of mathematical model of degradation of functional parameter was solved. It is offered to receive model in the form of conditional density of distribution of the considered parameter. It has allowed to apply a method of group forecasting.
Литература
1. European Organization of the Quality Control Glossary. Bern: EOQC. 1988.
2. Боровиков С.М. Теоретические основы конструирования, технологии и надежности. Минск, 1998.
3. СыноровВ.Ф., ПивоваровР.П., ПетровБ.К. и др. Под ред. Ю.Г. Миллера. Физические основы надежности интегральных схем. М., 1976.
4. МихайловА.В. Эксплуатационные допуски и надежность в радиоэлектронной аппаратуре. М., 1970.
5. Вентцель Е.С. Теория вероятностей. М., 1969.
6. Peattie D.S., Adams J.D., Carreli S.L. et al. // Proc. IEEE. 1974. Vol. 62, № 2. P. 149-168.
7. Quick Logic Reliability Report / рASIC, Vialink and Quick Logic Corp. Orleans, 1998.
8. Reliability Audit Report 1999. Life Test Data. ON Semiconductor L.L.C., Formerly a Division of Motorola, 1999.
