Выпуск #2/2024
А. Строгонов
МЕТОДЫ ПРОГНОЗИРОВАНИЯ ДОЛГОВЕЧНОСТИ ИС ПО ПАРАМЕТРИЧЕСКИМ ОТКАЗАМ
МЕТОДЫ ПРОГНОЗИРОВАНИЯ ДОЛГОВЕЧНОСТИ ИС ПО ПАРАМЕТРИЧЕСКИМ ОТКАЗАМ
Просмотры: 976
DOI: 10.22184/1992-4178.2024.233.2.88.96
В статье рассмотрены основные методологические подходы, принятые в зарубежной и отечественной инженерных школах в области прогнозирования долговечности ИС по параметрическим отказам.
В статье рассмотрены основные методологические подходы, принятые в зарубежной и отечественной инженерных школах в области прогнозирования долговечности ИС по параметрическим отказам.
Теги: aging effects degradation durability parametric failures prediction reliability деградация долговечность надежность параметрические отказы прогнозирование эффекты старения
Методы прогнозирования долговечности ИС по параметрическим отказам
А. Строгонов, д.т.н.
Современные радиоэлектронные системы для ответственных применений рассчитаны на эксплуатацию в активном или дежурном режимах на протяжении 15–20 лет. Продление срока службы этих дорогостоящих систем требует научного обоснования, которое базируется на прогнозировании долговечности используемых в них высоконадежных ИС. В настоящее время в России разрабатывают новый подход к оценке надежности ИС и полупроводниковых приборов с помощью математических методов прогнозирования. Минпромторг России в 2023 году выделил на эти цели более 450 млн руб. Проводимые исследования должны базироваться на новейших российских и мировых достижениях в таких областях, как анализ характеристик надежности, теория вероятностей, математическая статистика и автоматизированные методы проектирования изделий микроэлектроники с учетом эффектов старения. В статье рассмотрены основные методологические подходы, принятые в зарубежной и отечественной инженерных школах в области прогнозирования долговечности ИС по параметрическим отказам.
Введение
Спектр методов прогнозирования долговечности ИС по параметрическим отказам очень широк: методы теории катастроф и нечетких множеств, прогнозирование методами распознавания образов, имитационная физико-статистическая теория надежности, методы, основанные на анализе и прогнозировании рядов деградации контролируемых параметров, методы параметрической теории надежности (статистические и вероятностно-статистические методы прогнозирования), термодинамические методы прогнозирования, энергетически-временная модель деградации параметров ИС, эволюционно-энтропийная модель процесса деградации, феноменологический метод физического анализа и прогнозирования надежности, методы идентификации систем (ARX и ARMAX-модели) и др.
Многие теоретические исследования не были доведены до практического их применения. Почти все методы требуют разработки дополнительного программного обеспечения.
Единственный подход, который широко используется на практике и имеет строгое математическое обоснование, это подход с использованием статистических и вероятностно-статистических методов прогнозирования. Однако, для современных ИС с интенсивностью отказов порядка 10–9 1/ч и менее, характеризующихся преобладанием параметрических отказов над катастрофическими, более эффективным оказался подход с использованием имитационной физико-статистической теории надежности.
Под «старением» в статистической теории надежности понимают явление постепенной утраты одного или нескольких свойств в расчетных условиях эксплуатации [1]. Существует и другое определение термина «старение». Например, в [2, 3] под старением подразумевают ухудшение показателей надежности с течением времени, то есть старение может означать возрастание с течением времени интенсивности отказов.
Старение характеризуется скоростью потери работоспособности, определяемой как внутренними индивидуальными особенностями конструкции ИС, так и внешними условиями эксплуатации. В отличие от старения, процесс постепенной утраты работоспособности в результате неблагоприятного сочетания внешних условий называется износом [2]. Если старение характеризуется какой-то усредненной скоростью потери работоспособности, то износ может быть определен как старение со скоростью, превышающей номинальную. Считается, что износ для ИС не характерен, так как при эксплуатации ИС в радиоэлектронной аппаратуре не наблюдается период старения, характерный для кривой интенсивности отказов ИС, который выявляется только по результатам ускоренных испытаний [4].
Установлено, что деградация является отражением процессов перестройки дефектной структуры ИС. Кинетика деградации определяется условиями и режимом эксплуатации, при этом дефектная структура перестраивается, приобретая свойства, отвечающие измененным внешним условиям [5]. Параметры ИС являются прямым отражением состояния дефектной структуры и ее эволюции. Вносимые в дефектную структуру изменения и соответственно наблюдающиеся изменения параметров при эксплуатации могут быть либо обратимыми, либо необратимыми [5]. Считается, что после изготовления полупроводникового прибора его дефектная структура находится в потенциальной яме на вершине холма, высота которого характеризует степень неравновесности дефектной структуры, а крутизну склонов определяют условия эксплуатации [5].
На рис. 1 условно показан потенциальный рельеф, описывающий жизненный цикл полупроводникового прибора [5]. Срок службы прибора определяется эффективной длиной радиуса-вектора, проведенного из центра потенциальной ямы на уровень основания холма. При этом ориентация радиуса-вектора характеризует особенности эксплуатации полупроводникового прибора в выбранном режиме. Срок службы прибора будет максимальным при хранении в комфортных условиях (вектор a) или эксплуатации в щадящем режиме (вектор b), когда для перехода в равновесное состояние дефектная структура должна преодолеть высокий потенциальный барьер. При эксплуатации в штатном режиме (вектор c) величина потенциального барьера ниже. Это, однако, не означает, что в данном случае срок службы прибора будет меньше.
Для необратимого изменения параметров необходимо накопление определенной концентрации дефектов, присущих выбранному режиму эксплуатации. Происходящая при этом трансформация дефектной структуры может перевести ее в другую потенциальную яму, которых на пути эволюции может быть несколько [5].
Считается, что естественное старение ИС происходит в процессе их нормальной эксплуатации и в условиях хранения при оптимальных параметрах среды. Вынужденная деградация происходит под воздействием дестабилизирующих факторов, таких как радиация, механические и тепловые нагрузки и др. [6, 7].
В работе [6] развита концепция естественного старения ИС, основанная на учете электрон-решеточных взаимодействий, вызванных электронными возбуждениями. Выявлены механизмы, приводящие к старению двух типов: макроскопическое размытие примесных профилей (p-n-переходов) и микроскопическая перестройка дефектно-примесных комплексов.
Получение оценок долговечности ИС по результатам изучения процесса деградации параметров ИС представляется достаточно сложной задачей. Эволюционные процессы, протекающие в дефектной структуре, могут сопровождаться и улучшением параметров ИС. Зачастую многие важные электрические параметры ИС не изменяют своего значения вплоть до отказа, другие же не несут информацию, необходимую для прогнозирования долговечности, а лишь указывают на то, что схема в какой-то момент времени является работоспособной. Поэтому для предсказания долговечности приходится выбирать параметры, лишь косвенно характеризующие естественное старение, или те параметры, контроль которых является обязательным по ТУ.
На механизмы и скорость деградационных процессов оказывают влияние различные случайные возмущения (дефекты), исходно имеющиеся или возникающие в элементах конструкции ИС на различных этапах их жизненного цикла. В зависимости от степени влияния (размера) этих дефектов характер проявления отказов может носить детерминированный, полуслучайный или случайный характер. Если предположить, что дефекты отсутствуют или не оказывают влияния на скорость или механизмы деградационных процессов, то проявления отказов будут носить детерминированный характер, а изучение надежности ИС при известных механизмах отказов может проводиться аналитическими методами [8].
Если размеры дефектов малы, то есть их влияние сказывается только на скорости деградации, а не на механизме деградации, то проявление отказов может быть описано методами теории функций случайных аргументов [8]. Если в изделии имеются «грубые» дефекты, то есть дефекты, изменяющие механизмы и скорость деграда-
ционных процессов, то проявление отказов будет носить чисто случайный характер, и поэтому изучение надежности должно проводиться на основе методов теории вероятности и математической статистики [8].
Согласно работам [8–11] в теории надежности ИС выделяют три методологических подхода – физический, физико-статистический и статистический (вероятностный), которые позволяют прогнозировать количественные показатели надежности, такие как интенсивность отказов, вероятность безотказной работы и др., а также время наступления параметрических отказов в выборке. Однако существуют и другие подходы. Это термодинамический подход [12], подход с использованием теории катастроф и нечетких множеств [13] и др. Из всего многообразия предложенных подходов прогнозирования долговечности ИС физико-статистический подход, основанный на решении физико-статистического уравнения, является наиболее законченным.
Разделяют краткосрочные и долгосрочные прогнозы технического состояния РЭА. Краткосрочное прогнозирование основано на изучении эволюции контролируемых технических характеристик исследуемого объекта. Техническую характеристику (параметр) рассматривают как случайный процесс. Для оценки вероятности выхода параметра за поле допуска используют различные математические методы, основанные на отыскании распределения экстремальных значений, анализе многомерных плотностей распределений значений контролируемого параметра, канонических разложениях случайного процесса, марковской аппроксимации (дискретной модели) поведения контролируемых параметров. При долгосрочном прогнозировании технического состояния сложных РЭА, как правило, рассчитывают только показатели надежности. Это связано с трудностями прогнозирования поведения контролируемых параметров РЭА. Для ИС такой классификации не проводят.
Физико-статистический подход
при прогнозировании долговечности ИС
по параметрическим отказам
В настоящее время в теории надежности уделяется большое внимание поиску выражений для основных количественных показателей (физико-статистический подход [8–11]), сущность которого заключается в предположении о доминирующем влиянии на надежность естественных деградационных процессов.
Пусть работоспособная ИС характеризуется N-мерным вектором технических характеристик X с компонентами X1, X2, …, XN. Компонентами вектора технических характеристик X ИС могут быть, например, входное напряжение высокого уровня (UIH), входное напряжение низкого уровня (UIL), выходное напряжение низкого уровня (UOL) и др. Параметры режимов и условия эксплуатации, деградационные процессы и статистические характеристики дефектов, на которые может оказываться прямое воздействие, рассматриваются как управляющие переменные. Обозначим эти переменные вектором U. Поскольку влияние возмущающих параметров проявляется фактически в виде случайных отклонений технических характеристик изделий от их номинальных значений, то эти параметры учитываются путем придания переменным состояния вероятностно-статистического смысла. Возмущения приводят к тому, что значения переменных состояния становятся случайными величинами. Для описания последних вводится зависящая от времени и управляющих переменных N-мерная функция плотности их совместного распределения f (X, t, U). В качестве целевой функции выберем вероятность безотказной работы P (t). Выражение для целевой функции имеет следующий вид [8]:
P(t,U) = ∫V f(X,t,U)dV,
где V – область допустимых значений технических характеристик изделий с элементом объема dV = dX1dX2…dXN.
Необходимо найти такие значения управляющих переменных U, при которых значения количественных показателей надежности были бы в пределах заранее установленных границ. При заданных времени наработки t и минимально допустимой вероятности безотказной работы α необходимо отыскать область значений управляющих переменных D(U), удовлетворяющую следующему критерию эффективности [9]:
P(t,U) = ∫V f(X,t,U)dV ≥ α для всех U ∈ D(U).
Подынтегральная функция f(X, t, U) является решением физико-статистического уравнения непрерывности для значений технических характеристик изделий [9]:
f(X,t,U) + div(f(X,t,U)C(X,t,U)) = 0
с начальным условием f(X, 0, U) = f0(X, U), где f0(X, U) – функция плотности распределения исходных (в момент времени t = 0) значений технических характеристик X; div – оператор дивергенции N-мерный; C(X, t, U) = dX / dt – вектор скоростей деградации.
В общем случае решение задачи возможно получить на основе априорной информации аналитическими или численными методами. Сформулированная задача относится к классу задач нелинейного программирования. Однако общего решения уравнения для прогнозирования количественных показателей надежности серийно выпускаемых современных БИС нет. Тем не менее, для некоторых частных случаев, когда исходное распределение и скорость деградации можно аппроксимировать достаточно простыми аналитическими функциями, можно получить искомое решение.
В работе [10] показано применение физико-статистического подхода на примере деградации тока потребления КМОП БИС в течение наработки 30 тыс. ч. Исходная плотность распределения значений параметра Iпит имеет вид:
f0(Y) = , 0 ≤ Y ≤ 1,
где Y – вектор технических характеристик, состоящий из параметра Iпит; [0, 1] – границы нормированных значений Y; L – значение параметра распределения, равное 0,27. Скорость деградации этой характеристики аппроксимируется зависимостью вида:
,
где E – случайная, асимптотически нормально распределенная величина с характеристиками M(E) = 0, σ(E) = 0,14. В общем случае для достижения поставленной цели решается уравнение непрерывности [10]:
df / dt + div(fC) = 0,
f(Y, 0) = f0(Y).
Решение этого уравнения методом Коши есть функция плотности [10]:
f(Y,t,E) = .
Вероятность безотказной работы и интенсивность отказов имеют вид [10]:
P(t,E) = 1 – exp(–Z / L) / [1 – exp(–1 / L)],
λ(t,E) = Z2 exp(–Z / L) / [L(1 + E)(t + e)(1 – exp(–1 / L))],
Z = 1 + E / [E + ln(t + e)].
Таким образом, решение уравнения непрерывности позволяет получить важный вывод: скорость деградации является неявной функцией времени, то есть деградация носит автономный характер. При t→∞ λ(t) не зависит от исходного распределения значений технических характеристик f0(Y). Следовательно, влияние исходных статистических распределений на интенсивность отказов весьма значительно при малых временах наработки и существенно уменьшается при больших. Напротив, разброс скоростей деградационных изменений в материалах ИС гораздо сильнее проявляется при больших временах наработки, чем при малых [10].
В работе [8] авторы приводят вероятность безотказной работы и выражение для интенсивности отказов λ с учетом внезапных и постепенных отказов:
P(t) = exp(–∫0tq(t)dt) ∫V0 f0 (Z)(dZ / dX)dX1 dX2 … dXN,
,
где q(t) – имеет смысл интенсивности внезапных отказов, обусловленных грубыми дефектами; V0 – элемент объема технических характеристик ΔV = ΔX1, ΔX2…, ΔXN; Z = (Z1, Z2, …, ZN) – N-мерный вектор частных интегралов, удовлетворяющих условию Z(X, 0) = X системы уравнений в обыкновенных производных: dX / dt = C(X, t); S0 – часть поверхности S, охватывающей объем V0. Интегрирование проводится по поверхности S0, охватывающей объем V0, для которого скалярное произведение CdS положительно, что отражает условие о необратимости отказов. Для высоконадежных изделий предполагается отсутствие «грубых» дефектов и q(t) = 0. Тогда первые члены выражений для P(t) и λ(t) равны 1 и 0 соответственно, а сами выражения дают физико-статистические определения количественных показателей надежности.
Вероятностный и вероятностно-физический методы прогнозирования надежностных характеристик ИС
Долговечность ИС в условиях эксплуатации может быть оценена с помощью функций параметрической надежности [14]. Параметрическая надежность характеризует способность ИС сохранять уровень определяющих параметров X1, X2, …, Xn в границах соответствующих допусков (a1, b1), (a2, b2), …, (an, bn) в течение требуемого времени (0, t) при заданных режимах и условиях эксплуатации. Критерием параметрического отказа служит критерий нарушения условия: ai ≥ Xi > bi, i = 1, 2, …, n, а количественной мерой надежности – вероятность безотказной работы в течение времени (0, t) или функция параметрической надежности: P{ai ≥ Xi > bi}, i = 1, …, n.
На практике обычно используют некоторый обобщенный параметр X(t) и рассматривают способность ИС сохранять уровень обобщенного параметра X(t) в поле допуска в течение требуемого времени. Тогда математическое выражение для критерия отказа: a ≥ X(t) > b. В качестве параметра X(t) могут выступать не только контролируемые параметры ИС, но и параметры, характеризующие процесс деградации элементов конструкции ИС, например, изменение сопротивления токопроводящих дорожек алюминиевой металлизации вследствие электромиграции.
Обозначим вероятность выхода параметра X(t) за границу поля допуска на интервале [t, t + Δt] через PX(t), при условии, что в начальный момент времени t0 параметр находился внутри поля допуска. Если процесс деградации по характеру поведения математического ожидания X(t) является монотонным, то вероятность выполнения условия выхода за пределы допуска можно записать в виде:
dPX(t) = PX(t + dt) – PX(t) = –[PX(a < x < b, t + dt) – PX(a < x < b, t)] =
= –d[FX (x,t)] | ba = dFX(x,t) | ab, где FX(x, t) – закон распределения случайного значения параметра x в момент времени t. В теории случайных процессов FX(x, t) трактуют как функцию распределения вероятностей, которая определяется как вероятность события, заключающегося в том, что наблюдаемая случайная величина меньше или равна допустимому ее значению x, то есть FX(x) = PX(X ≤ x). Если FX(x, t) – нормальный случайный процесс, то:
,
где u = (x – mx(t)) / σx(t); mx(t) – математическое ожидание случайной величины X; sx(t) – среднеквадратичное отклонение.
В данном случае математическое ожидание mx является технологически предопределенным значением параметра, к которому тяготеют ИС в партии. Стабильность математического ожидания для разных партий – свидетельство стабильности технологического процесса производства ИС. Напротив, существенные отличия в парамет-
рическом поведении (нестабильность математического ожидания, существенные различия дисперсии параметра в разных партиях, непостоянство формы распределения) должны настораживать. Это признаки потенциальной ненадежности ИС.
Необходимо отметить сложность получения функций параметрической надежности. На практике функции параметрической надежности аппроксимируют нор-
мальным распределением, которое удачно используется на протяжении многих лет для прогнозирования количественных показателей надежности ИС.
Необходимо отметить, что характеристики параметрической надежности можно получить, не доводя исследуемые ИС до отказа, а выяснив лишь характерную тенденцию временного поведения определяющего параметра. Это дает возможность сократить длительность испытаний ИС на долговечность. Тем не менее, в литературе отсутствуют практические примеры прогнозирования долговечности ИС по параметрическим отказам с использованием функций параметрической надежности.
Для прогнозирования надежностных характеристик ИС по результатам ускоренных испытаний используют вероятностно-физические модели (диффузионные расп-
ределения), которые выгодно отличаются от строго вероятностных моделей тем, что их параметры могут быть вычислены как на основе статистики отказов, так и на основе статистики процесса деградации, а также путем совместного использования обоих типов информации [15, 16].
Элементы теории «диффузионных процессов» (Броуновский процесс) могут быть использованы для создания методов прогнозирования надежностных характеристик ИС. В настоящее время используются четыре вида двухпараметрических вероятностно-физических моделей отказов – альфа-распределение («веерный» процесс), нормальное («сильно перемешанный» гауссовский процесс), DM- (diffusive monotonic, марковский монотонный процесс) и DN- (diffusive non-monotonic, марковский немонотонный процесс) распределения [15, 16].
Прогнозирование долговечности ИС
с помощью теории катастроф
и теории нечетких множеств
Катастрофой называется скачкообразное изменение выходного параметра при плавном изменении управляющего параметра. В сущности, это понятие тождественно понятию отказа в теории надежности применительно к катастрофическим отказам. Теория катастроф базируется на теории особенностей гладких отображений (проектирование на плоскость каких-либо тел, например проектирование сферы на плоскость дает складку) и теории бифуркации динамических систем. При проектировании годографа вектора состояния ИС (считается, что эволюция процесса деградации определяется векторным полем в фазовом пространстве, вектор – параметр состояния ИС) на плоскость, образуются особенности типа складки и сборки – особые (бифуркационные) точки. Проекции точек складок или сборок образуют кривую катастроф [12].
Задача состоит в следующем. Рассматривается система, динамика которой принадлежит к градиентному типу, и делается попытка минимизировать некоторую функцию. В распоряжении имеется k параметров управления α1, α2, …, αk, и при этом выходные параметры системы принимают такие значения X1*, X2*, …, Xn* в состоянии равновесия, что достигается локальная минимизация некоторой функции Ψ(X1, X2, …, Xn; α1, α2, …, αk), где Ψ является потенциальной функцией системы. Значения Xi*, соответствующие состоянию равновесия, зависят от выбора параметров αk. Рассматривая скачкообразные изменения значений Xi*, происходящие в результате плавного изменения параметров αk, приходим к понятию катастрофа. В общем случае функция Ψ неизвестна. Предполагается, что она является потенциалом, описывающим динамику данной системы.
Прогнозирующие модели процесса деградации технических характеристик ИС возможно представить с использованием разных канонических форм функции Ψ(x, α), имеющих полиномиальную структуру, в виде: X13 + αX1 – складка; X14 + α1(X12 / 2) + α2X1 – сборка; (X15 / 5) + α1(X13 / 3) + α2(X12 / 2) + α3X1 – ласточкин хвост. Критические точки X1*(α),X2*(α) и т.д. получаются из решений уравнений ∂Ψ / ∂X1 = 0; ∂Ψ / ∂X2 = 0 и т.д.
Также предлагается использовать при прогнозировании долговечности ИС теорию нечетких множеств в виду большого числа влияющих факторов на долговечность ИС и недостаточно полного (четкого) знания о влиянии этих факторов на долговечность ИС [13].
Прогнозирование деградации
параметров ИС с использованием
теории временных рядов
Процесс деградации параметров ИС может быть эффективно описан с использованием метода Бокса – Дженкинса (АРПСС-модель – авторегрессия проинтегрированного скользящего среднего, или ARIMA-модель), применяемого для анализа и прогнозирования нестационарных случайных процессов и временных рядов различной природы [17, 18]. В настоящее время метод Бокса – Дженкинса считается одним из наиболее эффективных при прогнозировании поведения параметров объектов и динамических систем, описываемых случайными временными рядами. Метод нашел широкое применение для стохастического моделирования интенсивности отказов деталей узлов различного рода радиоэлектронного оборудования, высокочастотных генераторов, используемых в бортовой аппаратуре, для построения прогнозов в рамках теории управления производственными процессами, в эконометрии.
Одним из преимуществ метода Бокса – Дженкинса перед другими методами прогнозирования является то, что он позволяет построить модели временных рядов деградации параметров ИС, которые можно использовать для прогнозирования долговечности ИС по параметрическим отказам. Метод также позволяет прогнозировать долговечность ИС независимо от ее степени интеграции, конструкционных особенностей и учитывать влияние случайной составляющей (белый шум) на процесс деградации параметров. В отличие от методов теории распознавания не требуется предварительного исследования выборки, то есть разделения ИС в пространстве параметров на ряд совокупностей «группы схожести», проведения «обучения» и «распознавания».
Следует отметить, что существует хорошо развитое программное обеспечение для этого метода, способствующее его широкому применению в различных областях науки и техники. Для прогнозирования долговечности ИС возможно использовать статистические пакеты программ (СПП), позволяющие проводить анализ временных рядов, такие как пакет STATGRAPHICS (Statistical Graphics System), разработанный корпорацией STSC (США), пакет Statistica for Windows от Stat Soft, пакет GENSTAT, разработанный в Numerical Algorithms Group, пакет SPSS (Statistical Package for the Social Science) – статистический пакет для общественных наук и др. Разработан ряд отечественных программ анализа и прогнозирования временных рядов, такие как ПАРИС и МАВР, STANP. Модуль СПП SPSS Trends позволяет с использованием ARIMA-моделей решать задачи управления производством и обработки данных, прогнозировать объемы продаж, проводить социальные исследования. Программный продукт Forecast Expert (система построения экономических прогнозов), базирующийся на методе Бокса – Дженкинса, предназначен для прогнозирования любого параметра, в отношении которого имеется должное количество замеров в конкретном промежутке времени. Всего насчитывается более десятка таких статистических пакетов.
Применение метода Бокса – Дженкинса для изучения долговечности ИС есть результат подхода к исследуемой ИС как к «черному ящику». Метод «черного ящика» – кибернетический: объект исследования представляется как некоторая кибернетическая система, которая может быть описана своим функциональным оператором. Такой подход ставит своей целью посредством построения некоторой модели установить изоморфизм не с внутренней структурой и ее функционированием, а с внешними проявлениями ее информативных параметров. Недостатком метода является то, что АРПСС-модель вытекающая из метода Бокса – Дженкинса не основывается на физических представлениях.
Метод Бокса – Дженкинса основывается на том, что гладкий нестационарный временной ряд путем взятия разностей некоторого d-го порядка можно свести к эквивалентному стационарному, то есть к случаю, для которого разработаны методы анализа и прогнозирования. В методе Бокса – Дженкинса нестационарный временной ряд представляется в виде модели авторегрессии проинтегрированного скользящего среднего (АРПСС, ARIMA) в прямой и возвратной форме:
прямая форма: φ(B)∇d Xt = θ(B)at, t = +1, +2, …,
возвратная форма: φ(F)∇d Xt = θ(F)et, t = –1, –2, …,
φ(B) = 1 – φ1B – φ2 B2 – ... – φP BP = 1 – ∑j=1p φj B j = АР-оператор,
θ(B) = 1 – θ1 B – θ2 B2 – ... – θq Bq=1 – ∑i=1q θi Bi = CC-оператор,
где φ(B), φ(F) – операторы авторегрессии (АР); B, F – операторы сдвига; ∇ – оператор разности: ∇1Xt = Xt–Xt–1 = (1 – B)Xt; d – порядок разности, обеспечивающий переход от нестационарного ряда к эквивалентному стационарному; p – порядок авторегрессии; q, θ(B), θ(F) – порядок и операторы проинтегрированного скользящего среднего (СС) соот-
ветственно; at, et – последовательности независимых случайных величин, имеющих одинаковое нормальное распределение (белый шум): E[at] = 0, Var[at] = σat. Значения белого шума at, согласно методу Бокса – Дженкинса, должны быть получены прогнозированием назад с использованием возвратной формы временного ряда.
Смешанный процесс АР и СС представляет собой эффективную модель для описания стационарных и нестационарных временных рядов. АРПСС-модель эффективна, когда временной ряд является суммой независимых составляющих, каждая из которых описывается либо моделью АР, либо моделью СС, но которые непосредственно не наблюдаются. Модель является параметрической. Чтобы подобрать такую модель, нужно оценить по наблюдаемым данным набор параметров АРПСС-модели φ, θ, σa.
Согласно авторегрессионной модели формирования сигнала, он формируется путем пропускания дискретного белого шума через рекурсивный фильтр, в модели скользящего среднего (СС-модель) используется нерекурсивный фильтр, а в АРПСС-модели – фильтр общего вида, содержащий рекурсивную и нерекурсивную ветви. Фильтры АРПСС имеют обратную связь, обладают большой гибкостью, экономичны, то есть малое число параметров в левых и правых частях.
Заключение
В статье кратко рассмотрены различные подходы, которые возможно использовать для прогнозирования долговечности ИС по параметрическим отказам, базирующие-
ся на методах, заимствованных из различных областей науки и техники, таких как статистика и теория вероятностей, теория распознавания образов и теория ка-
тастроф, теория временных рядов, термодинамика и др. С математической точки зрения применение этих методов для прогнозирования долговечности ИС по парамет-
рическим отказам не имеют ограничений.
Следует отметить, что в настоящий момент ни одним из рассмотренных методов прогнозирования долговечности ИС по параметрическим отказам не получены сведения о фактической долговечности ИС в условиях эксплуатации РЭА. Однако в современной работе [19] проведена верификация физико-статистического подхода теории надежности для простейших случаев, показавшая его правомочность.
Обзор литературных источников показывает, что наиболее отработанными с точки зрения математического аппарата и программного обеспечения являются методы, основанные на теории временных рядов и позволяющие по траекториям процесса деградации параметров ИС прогнозировать долговечность ИС.
ЛИТЕРАТУРА
Ушаков И.А. Вероятностные модели надежности
информационно-вычислительных систем. М.:
Радио и связь, 1991. 132 с.
Барлоу Р., Прошан Ф. Математическая теория
надежности / Пер. с англ.; под ред. Б.В. Гнеденко.
М.: Сов. радио, 1969. 537 с.
Барлоу Р., Прошан Ф. Статистическая теория
надежности и испытания на безотказность /
Пер. с англ. И.А. Ушакова. М.: Наука, 1984. 452 с.
Строгонов А.В. Характеристики надежности
современных ПЛИС // ЭЛЕКТРОНИКА: Наука,
Технология, Бизнес. 2019. № 4. С. 52–58.
Соколов В.И. Проблемы микроэлектроники
(1. Диффузия. 2. Дефектообразование. 3. Деградация) // Физика и техника полупроводников. 1995. Т. 29. С. 842.
Горин Б.М., Кив А.Е, Плотникова Л.Г. и др.
Механизмы естественного старения и вынужденной
деградации полупроводниковых приборов // Обзоры
по электронной техники, 1983. 57 с.
Серия «Полупроводниковые приборы».
Вавилов В.С., Горин Б.М., Данилин Н.С. и др.
Радиационные методы в твердотельной электронике.
М.: Радио и связь, 1990. 184 с.
Алексанян И.Т., Черняев Н.В. Обобщенная модель
надежности изделий // Известия вузов. Электроника. 1998. № 1. C. 85–90.
Алексанян И.Т., Черняев Н.В. Выражения для основных количественных показателей надежности в физико-
статистическом подходе // Петербургский журнал
электроники. 1994. № 1. C. 56–58.
Алексанян И.Т., Черняев Н.В.
Метод изучения надежности интегральных микросхем // Микроэлектроника. 1992. Вып. 2. Т. 21. C. 105–111.
Алексанян И.Т., Черняев Н.В. Управление надежностью интегральных микросхем на основе информационной
избыточности // Известия вузов. Электроника. 1998. № 7. C. 62–66.
Малков Я.В. Диагностика деградационных процессов
в полупроводниковых структурах с позиций междисциплинарных теорий // Зарубежная радиоэлектроника. 1996. № 5. C. 77–82.
Дубицкий Л.Г. Предвестники отказов в изделиях
электронной техники. М.: Радио и связь, 1989. 96 с.
Сыноров В.Ф., Пивоварова Р.П. Параметрическая
надежность и физические модели отказов интегральных схем. Воронеж: Изд-во ВГУ, 1983. 152 с.
Погребинский С.Б., Стрельников В.П.
Проектирование и надежность многопроцессорных ЭВМ. М.: Радио и связь, 1988. 168 с.
Стрельников В.П. Вероятностно-физические методы
исследования надежности машин и аппаратуры //
Надежность и контроль качества. 1989. № 9. С. 3–7.
Алексанян И.Т., Кривошапко В.М.,
Романов А.А. Построение модели параметрической
надежности ИС по данным о деградации
их характеристик и имитация испытаний на ЭВМ //
Электронная техника. 1979. Вып. 1. С. 15–19.
Серия «Микроэлектроника».
Алексанян И.Т., Кривошапко В.М. Моделирование
параметрических отказов и изучение надежности
интегральных схем // Электронная техника. 1981. Вып. 4. С. 52–57. Серия «Управление качеством, стандартизация, метрология, испытания».
Кожевников В.С., Матюшкинa И.В., Черняев Н.В.
Анализ основного уравнения физико-статистического подхода теории надежности технических систем //
Математические основы и численные методы
моделирования. 2020. Т. 12. № 4. С. 721–735.
А. Строгонов, д.т.н.
Современные радиоэлектронные системы для ответственных применений рассчитаны на эксплуатацию в активном или дежурном режимах на протяжении 15–20 лет. Продление срока службы этих дорогостоящих систем требует научного обоснования, которое базируется на прогнозировании долговечности используемых в них высоконадежных ИС. В настоящее время в России разрабатывают новый подход к оценке надежности ИС и полупроводниковых приборов с помощью математических методов прогнозирования. Минпромторг России в 2023 году выделил на эти цели более 450 млн руб. Проводимые исследования должны базироваться на новейших российских и мировых достижениях в таких областях, как анализ характеристик надежности, теория вероятностей, математическая статистика и автоматизированные методы проектирования изделий микроэлектроники с учетом эффектов старения. В статье рассмотрены основные методологические подходы, принятые в зарубежной и отечественной инженерных школах в области прогнозирования долговечности ИС по параметрическим отказам.
Введение
Спектр методов прогнозирования долговечности ИС по параметрическим отказам очень широк: методы теории катастроф и нечетких множеств, прогнозирование методами распознавания образов, имитационная физико-статистическая теория надежности, методы, основанные на анализе и прогнозировании рядов деградации контролируемых параметров, методы параметрической теории надежности (статистические и вероятностно-статистические методы прогнозирования), термодинамические методы прогнозирования, энергетически-временная модель деградации параметров ИС, эволюционно-энтропийная модель процесса деградации, феноменологический метод физического анализа и прогнозирования надежности, методы идентификации систем (ARX и ARMAX-модели) и др.
Многие теоретические исследования не были доведены до практического их применения. Почти все методы требуют разработки дополнительного программного обеспечения.
Единственный подход, который широко используется на практике и имеет строгое математическое обоснование, это подход с использованием статистических и вероятностно-статистических методов прогнозирования. Однако, для современных ИС с интенсивностью отказов порядка 10–9 1/ч и менее, характеризующихся преобладанием параметрических отказов над катастрофическими, более эффективным оказался подход с использованием имитационной физико-статистической теории надежности.
Под «старением» в статистической теории надежности понимают явление постепенной утраты одного или нескольких свойств в расчетных условиях эксплуатации [1]. Существует и другое определение термина «старение». Например, в [2, 3] под старением подразумевают ухудшение показателей надежности с течением времени, то есть старение может означать возрастание с течением времени интенсивности отказов.
Старение характеризуется скоростью потери работоспособности, определяемой как внутренними индивидуальными особенностями конструкции ИС, так и внешними условиями эксплуатации. В отличие от старения, процесс постепенной утраты работоспособности в результате неблагоприятного сочетания внешних условий называется износом [2]. Если старение характеризуется какой-то усредненной скоростью потери работоспособности, то износ может быть определен как старение со скоростью, превышающей номинальную. Считается, что износ для ИС не характерен, так как при эксплуатации ИС в радиоэлектронной аппаратуре не наблюдается период старения, характерный для кривой интенсивности отказов ИС, который выявляется только по результатам ускоренных испытаний [4].
Установлено, что деградация является отражением процессов перестройки дефектной структуры ИС. Кинетика деградации определяется условиями и режимом эксплуатации, при этом дефектная структура перестраивается, приобретая свойства, отвечающие измененным внешним условиям [5]. Параметры ИС являются прямым отражением состояния дефектной структуры и ее эволюции. Вносимые в дефектную структуру изменения и соответственно наблюдающиеся изменения параметров при эксплуатации могут быть либо обратимыми, либо необратимыми [5]. Считается, что после изготовления полупроводникового прибора его дефектная структура находится в потенциальной яме на вершине холма, высота которого характеризует степень неравновесности дефектной структуры, а крутизну склонов определяют условия эксплуатации [5].
На рис. 1 условно показан потенциальный рельеф, описывающий жизненный цикл полупроводникового прибора [5]. Срок службы прибора определяется эффективной длиной радиуса-вектора, проведенного из центра потенциальной ямы на уровень основания холма. При этом ориентация радиуса-вектора характеризует особенности эксплуатации полупроводникового прибора в выбранном режиме. Срок службы прибора будет максимальным при хранении в комфортных условиях (вектор a) или эксплуатации в щадящем режиме (вектор b), когда для перехода в равновесное состояние дефектная структура должна преодолеть высокий потенциальный барьер. При эксплуатации в штатном режиме (вектор c) величина потенциального барьера ниже. Это, однако, не означает, что в данном случае срок службы прибора будет меньше.
Для необратимого изменения параметров необходимо накопление определенной концентрации дефектов, присущих выбранному режиму эксплуатации. Происходящая при этом трансформация дефектной структуры может перевести ее в другую потенциальную яму, которых на пути эволюции может быть несколько [5].
Считается, что естественное старение ИС происходит в процессе их нормальной эксплуатации и в условиях хранения при оптимальных параметрах среды. Вынужденная деградация происходит под воздействием дестабилизирующих факторов, таких как радиация, механические и тепловые нагрузки и др. [6, 7].
В работе [6] развита концепция естественного старения ИС, основанная на учете электрон-решеточных взаимодействий, вызванных электронными возбуждениями. Выявлены механизмы, приводящие к старению двух типов: макроскопическое размытие примесных профилей (p-n-переходов) и микроскопическая перестройка дефектно-примесных комплексов.
Получение оценок долговечности ИС по результатам изучения процесса деградации параметров ИС представляется достаточно сложной задачей. Эволюционные процессы, протекающие в дефектной структуре, могут сопровождаться и улучшением параметров ИС. Зачастую многие важные электрические параметры ИС не изменяют своего значения вплоть до отказа, другие же не несут информацию, необходимую для прогнозирования долговечности, а лишь указывают на то, что схема в какой-то момент времени является работоспособной. Поэтому для предсказания долговечности приходится выбирать параметры, лишь косвенно характеризующие естественное старение, или те параметры, контроль которых является обязательным по ТУ.
На механизмы и скорость деградационных процессов оказывают влияние различные случайные возмущения (дефекты), исходно имеющиеся или возникающие в элементах конструкции ИС на различных этапах их жизненного цикла. В зависимости от степени влияния (размера) этих дефектов характер проявления отказов может носить детерминированный, полуслучайный или случайный характер. Если предположить, что дефекты отсутствуют или не оказывают влияния на скорость или механизмы деградационных процессов, то проявления отказов будут носить детерминированный характер, а изучение надежности ИС при известных механизмах отказов может проводиться аналитическими методами [8].
Если размеры дефектов малы, то есть их влияние сказывается только на скорости деградации, а не на механизме деградации, то проявление отказов может быть описано методами теории функций случайных аргументов [8]. Если в изделии имеются «грубые» дефекты, то есть дефекты, изменяющие механизмы и скорость деграда-
ционных процессов, то проявление отказов будет носить чисто случайный характер, и поэтому изучение надежности должно проводиться на основе методов теории вероятности и математической статистики [8].
Согласно работам [8–11] в теории надежности ИС выделяют три методологических подхода – физический, физико-статистический и статистический (вероятностный), которые позволяют прогнозировать количественные показатели надежности, такие как интенсивность отказов, вероятность безотказной работы и др., а также время наступления параметрических отказов в выборке. Однако существуют и другие подходы. Это термодинамический подход [12], подход с использованием теории катастроф и нечетких множеств [13] и др. Из всего многообразия предложенных подходов прогнозирования долговечности ИС физико-статистический подход, основанный на решении физико-статистического уравнения, является наиболее законченным.
Разделяют краткосрочные и долгосрочные прогнозы технического состояния РЭА. Краткосрочное прогнозирование основано на изучении эволюции контролируемых технических характеристик исследуемого объекта. Техническую характеристику (параметр) рассматривают как случайный процесс. Для оценки вероятности выхода параметра за поле допуска используют различные математические методы, основанные на отыскании распределения экстремальных значений, анализе многомерных плотностей распределений значений контролируемого параметра, канонических разложениях случайного процесса, марковской аппроксимации (дискретной модели) поведения контролируемых параметров. При долгосрочном прогнозировании технического состояния сложных РЭА, как правило, рассчитывают только показатели надежности. Это связано с трудностями прогнозирования поведения контролируемых параметров РЭА. Для ИС такой классификации не проводят.
Физико-статистический подход
при прогнозировании долговечности ИС
по параметрическим отказам
В настоящее время в теории надежности уделяется большое внимание поиску выражений для основных количественных показателей (физико-статистический подход [8–11]), сущность которого заключается в предположении о доминирующем влиянии на надежность естественных деградационных процессов.
Пусть работоспособная ИС характеризуется N-мерным вектором технических характеристик X с компонентами X1, X2, …, XN. Компонентами вектора технических характеристик X ИС могут быть, например, входное напряжение высокого уровня (UIH), входное напряжение низкого уровня (UIL), выходное напряжение низкого уровня (UOL) и др. Параметры режимов и условия эксплуатации, деградационные процессы и статистические характеристики дефектов, на которые может оказываться прямое воздействие, рассматриваются как управляющие переменные. Обозначим эти переменные вектором U. Поскольку влияние возмущающих параметров проявляется фактически в виде случайных отклонений технических характеристик изделий от их номинальных значений, то эти параметры учитываются путем придания переменным состояния вероятностно-статистического смысла. Возмущения приводят к тому, что значения переменных состояния становятся случайными величинами. Для описания последних вводится зависящая от времени и управляющих переменных N-мерная функция плотности их совместного распределения f (X, t, U). В качестве целевой функции выберем вероятность безотказной работы P (t). Выражение для целевой функции имеет следующий вид [8]:
P(t,U) = ∫V f(X,t,U)dV,
где V – область допустимых значений технических характеристик изделий с элементом объема dV = dX1dX2…dXN.
Необходимо найти такие значения управляющих переменных U, при которых значения количественных показателей надежности были бы в пределах заранее установленных границ. При заданных времени наработки t и минимально допустимой вероятности безотказной работы α необходимо отыскать область значений управляющих переменных D(U), удовлетворяющую следующему критерию эффективности [9]:
P(t,U) = ∫V f(X,t,U)dV ≥ α для всех U ∈ D(U).
Подынтегральная функция f(X, t, U) является решением физико-статистического уравнения непрерывности для значений технических характеристик изделий [9]:
f(X,t,U) + div(f(X,t,U)C(X,t,U)) = 0
с начальным условием f(X, 0, U) = f0(X, U), где f0(X, U) – функция плотности распределения исходных (в момент времени t = 0) значений технических характеристик X; div – оператор дивергенции N-мерный; C(X, t, U) = dX / dt – вектор скоростей деградации.
В общем случае решение задачи возможно получить на основе априорной информации аналитическими или численными методами. Сформулированная задача относится к классу задач нелинейного программирования. Однако общего решения уравнения для прогнозирования количественных показателей надежности серийно выпускаемых современных БИС нет. Тем не менее, для некоторых частных случаев, когда исходное распределение и скорость деградации можно аппроксимировать достаточно простыми аналитическими функциями, можно получить искомое решение.
В работе [10] показано применение физико-статистического подхода на примере деградации тока потребления КМОП БИС в течение наработки 30 тыс. ч. Исходная плотность распределения значений параметра Iпит имеет вид:
f0(Y) = , 0 ≤ Y ≤ 1,
где Y – вектор технических характеристик, состоящий из параметра Iпит; [0, 1] – границы нормированных значений Y; L – значение параметра распределения, равное 0,27. Скорость деградации этой характеристики аппроксимируется зависимостью вида:
,
где E – случайная, асимптотически нормально распределенная величина с характеристиками M(E) = 0, σ(E) = 0,14. В общем случае для достижения поставленной цели решается уравнение непрерывности [10]:
df / dt + div(fC) = 0,
f(Y, 0) = f0(Y).
Решение этого уравнения методом Коши есть функция плотности [10]:
f(Y,t,E) = .
Вероятность безотказной работы и интенсивность отказов имеют вид [10]:
P(t,E) = 1 – exp(–Z / L) / [1 – exp(–1 / L)],
λ(t,E) = Z2 exp(–Z / L) / [L(1 + E)(t + e)(1 – exp(–1 / L))],
Z = 1 + E / [E + ln(t + e)].
Таким образом, решение уравнения непрерывности позволяет получить важный вывод: скорость деградации является неявной функцией времени, то есть деградация носит автономный характер. При t→∞ λ(t) не зависит от исходного распределения значений технических характеристик f0(Y). Следовательно, влияние исходных статистических распределений на интенсивность отказов весьма значительно при малых временах наработки и существенно уменьшается при больших. Напротив, разброс скоростей деградационных изменений в материалах ИС гораздо сильнее проявляется при больших временах наработки, чем при малых [10].
В работе [8] авторы приводят вероятность безотказной работы и выражение для интенсивности отказов λ с учетом внезапных и постепенных отказов:
P(t) = exp(–∫0tq(t)dt) ∫V0 f0 (Z)(dZ / dX)dX1 dX2 … dXN,
,
где q(t) – имеет смысл интенсивности внезапных отказов, обусловленных грубыми дефектами; V0 – элемент объема технических характеристик ΔV = ΔX1, ΔX2…, ΔXN; Z = (Z1, Z2, …, ZN) – N-мерный вектор частных интегралов, удовлетворяющих условию Z(X, 0) = X системы уравнений в обыкновенных производных: dX / dt = C(X, t); S0 – часть поверхности S, охватывающей объем V0. Интегрирование проводится по поверхности S0, охватывающей объем V0, для которого скалярное произведение CdS положительно, что отражает условие о необратимости отказов. Для высоконадежных изделий предполагается отсутствие «грубых» дефектов и q(t) = 0. Тогда первые члены выражений для P(t) и λ(t) равны 1 и 0 соответственно, а сами выражения дают физико-статистические определения количественных показателей надежности.
Вероятностный и вероятностно-физический методы прогнозирования надежностных характеристик ИС
Долговечность ИС в условиях эксплуатации может быть оценена с помощью функций параметрической надежности [14]. Параметрическая надежность характеризует способность ИС сохранять уровень определяющих параметров X1, X2, …, Xn в границах соответствующих допусков (a1, b1), (a2, b2), …, (an, bn) в течение требуемого времени (0, t) при заданных режимах и условиях эксплуатации. Критерием параметрического отказа служит критерий нарушения условия: ai ≥ Xi > bi, i = 1, 2, …, n, а количественной мерой надежности – вероятность безотказной работы в течение времени (0, t) или функция параметрической надежности: P{ai ≥ Xi > bi}, i = 1, …, n.
На практике обычно используют некоторый обобщенный параметр X(t) и рассматривают способность ИС сохранять уровень обобщенного параметра X(t) в поле допуска в течение требуемого времени. Тогда математическое выражение для критерия отказа: a ≥ X(t) > b. В качестве параметра X(t) могут выступать не только контролируемые параметры ИС, но и параметры, характеризующие процесс деградации элементов конструкции ИС, например, изменение сопротивления токопроводящих дорожек алюминиевой металлизации вследствие электромиграции.
Обозначим вероятность выхода параметра X(t) за границу поля допуска на интервале [t, t + Δt] через PX(t), при условии, что в начальный момент времени t0 параметр находился внутри поля допуска. Если процесс деградации по характеру поведения математического ожидания X(t) является монотонным, то вероятность выполнения условия выхода за пределы допуска можно записать в виде:
dPX(t) = PX(t + dt) – PX(t) = –[PX(a < x < b, t + dt) – PX(a < x < b, t)] =
= –d[FX (x,t)] | ba = dFX(x,t) | ab, где FX(x, t) – закон распределения случайного значения параметра x в момент времени t. В теории случайных процессов FX(x, t) трактуют как функцию распределения вероятностей, которая определяется как вероятность события, заключающегося в том, что наблюдаемая случайная величина меньше или равна допустимому ее значению x, то есть FX(x) = PX(X ≤ x). Если FX(x, t) – нормальный случайный процесс, то:
,
где u = (x – mx(t)) / σx(t); mx(t) – математическое ожидание случайной величины X; sx(t) – среднеквадратичное отклонение.
В данном случае математическое ожидание mx является технологически предопределенным значением параметра, к которому тяготеют ИС в партии. Стабильность математического ожидания для разных партий – свидетельство стабильности технологического процесса производства ИС. Напротив, существенные отличия в парамет-
рическом поведении (нестабильность математического ожидания, существенные различия дисперсии параметра в разных партиях, непостоянство формы распределения) должны настораживать. Это признаки потенциальной ненадежности ИС.
Необходимо отметить сложность получения функций параметрической надежности. На практике функции параметрической надежности аппроксимируют нор-
мальным распределением, которое удачно используется на протяжении многих лет для прогнозирования количественных показателей надежности ИС.
Необходимо отметить, что характеристики параметрической надежности можно получить, не доводя исследуемые ИС до отказа, а выяснив лишь характерную тенденцию временного поведения определяющего параметра. Это дает возможность сократить длительность испытаний ИС на долговечность. Тем не менее, в литературе отсутствуют практические примеры прогнозирования долговечности ИС по параметрическим отказам с использованием функций параметрической надежности.
Для прогнозирования надежностных характеристик ИС по результатам ускоренных испытаний используют вероятностно-физические модели (диффузионные расп-
ределения), которые выгодно отличаются от строго вероятностных моделей тем, что их параметры могут быть вычислены как на основе статистики отказов, так и на основе статистики процесса деградации, а также путем совместного использования обоих типов информации [15, 16].
Элементы теории «диффузионных процессов» (Броуновский процесс) могут быть использованы для создания методов прогнозирования надежностных характеристик ИС. В настоящее время используются четыре вида двухпараметрических вероятностно-физических моделей отказов – альфа-распределение («веерный» процесс), нормальное («сильно перемешанный» гауссовский процесс), DM- (diffusive monotonic, марковский монотонный процесс) и DN- (diffusive non-monotonic, марковский немонотонный процесс) распределения [15, 16].
Прогнозирование долговечности ИС
с помощью теории катастроф
и теории нечетких множеств
Катастрофой называется скачкообразное изменение выходного параметра при плавном изменении управляющего параметра. В сущности, это понятие тождественно понятию отказа в теории надежности применительно к катастрофическим отказам. Теория катастроф базируется на теории особенностей гладких отображений (проектирование на плоскость каких-либо тел, например проектирование сферы на плоскость дает складку) и теории бифуркации динамических систем. При проектировании годографа вектора состояния ИС (считается, что эволюция процесса деградации определяется векторным полем в фазовом пространстве, вектор – параметр состояния ИС) на плоскость, образуются особенности типа складки и сборки – особые (бифуркационные) точки. Проекции точек складок или сборок образуют кривую катастроф [12].
Задача состоит в следующем. Рассматривается система, динамика которой принадлежит к градиентному типу, и делается попытка минимизировать некоторую функцию. В распоряжении имеется k параметров управления α1, α2, …, αk, и при этом выходные параметры системы принимают такие значения X1*, X2*, …, Xn* в состоянии равновесия, что достигается локальная минимизация некоторой функции Ψ(X1, X2, …, Xn; α1, α2, …, αk), где Ψ является потенциальной функцией системы. Значения Xi*, соответствующие состоянию равновесия, зависят от выбора параметров αk. Рассматривая скачкообразные изменения значений Xi*, происходящие в результате плавного изменения параметров αk, приходим к понятию катастрофа. В общем случае функция Ψ неизвестна. Предполагается, что она является потенциалом, описывающим динамику данной системы.
Прогнозирующие модели процесса деградации технических характеристик ИС возможно представить с использованием разных канонических форм функции Ψ(x, α), имеющих полиномиальную структуру, в виде: X13 + αX1 – складка; X14 + α1(X12 / 2) + α2X1 – сборка; (X15 / 5) + α1(X13 / 3) + α2(X12 / 2) + α3X1 – ласточкин хвост. Критические точки X1*(α),X2*(α) и т.д. получаются из решений уравнений ∂Ψ / ∂X1 = 0; ∂Ψ / ∂X2 = 0 и т.д.
Также предлагается использовать при прогнозировании долговечности ИС теорию нечетких множеств в виду большого числа влияющих факторов на долговечность ИС и недостаточно полного (четкого) знания о влиянии этих факторов на долговечность ИС [13].
Прогнозирование деградации
параметров ИС с использованием
теории временных рядов
Процесс деградации параметров ИС может быть эффективно описан с использованием метода Бокса – Дженкинса (АРПСС-модель – авторегрессия проинтегрированного скользящего среднего, или ARIMA-модель), применяемого для анализа и прогнозирования нестационарных случайных процессов и временных рядов различной природы [17, 18]. В настоящее время метод Бокса – Дженкинса считается одним из наиболее эффективных при прогнозировании поведения параметров объектов и динамических систем, описываемых случайными временными рядами. Метод нашел широкое применение для стохастического моделирования интенсивности отказов деталей узлов различного рода радиоэлектронного оборудования, высокочастотных генераторов, используемых в бортовой аппаратуре, для построения прогнозов в рамках теории управления производственными процессами, в эконометрии.
Одним из преимуществ метода Бокса – Дженкинса перед другими методами прогнозирования является то, что он позволяет построить модели временных рядов деградации параметров ИС, которые можно использовать для прогнозирования долговечности ИС по параметрическим отказам. Метод также позволяет прогнозировать долговечность ИС независимо от ее степени интеграции, конструкционных особенностей и учитывать влияние случайной составляющей (белый шум) на процесс деградации параметров. В отличие от методов теории распознавания не требуется предварительного исследования выборки, то есть разделения ИС в пространстве параметров на ряд совокупностей «группы схожести», проведения «обучения» и «распознавания».
Следует отметить, что существует хорошо развитое программное обеспечение для этого метода, способствующее его широкому применению в различных областях науки и техники. Для прогнозирования долговечности ИС возможно использовать статистические пакеты программ (СПП), позволяющие проводить анализ временных рядов, такие как пакет STATGRAPHICS (Statistical Graphics System), разработанный корпорацией STSC (США), пакет Statistica for Windows от Stat Soft, пакет GENSTAT, разработанный в Numerical Algorithms Group, пакет SPSS (Statistical Package for the Social Science) – статистический пакет для общественных наук и др. Разработан ряд отечественных программ анализа и прогнозирования временных рядов, такие как ПАРИС и МАВР, STANP. Модуль СПП SPSS Trends позволяет с использованием ARIMA-моделей решать задачи управления производством и обработки данных, прогнозировать объемы продаж, проводить социальные исследования. Программный продукт Forecast Expert (система построения экономических прогнозов), базирующийся на методе Бокса – Дженкинса, предназначен для прогнозирования любого параметра, в отношении которого имеется должное количество замеров в конкретном промежутке времени. Всего насчитывается более десятка таких статистических пакетов.
Применение метода Бокса – Дженкинса для изучения долговечности ИС есть результат подхода к исследуемой ИС как к «черному ящику». Метод «черного ящика» – кибернетический: объект исследования представляется как некоторая кибернетическая система, которая может быть описана своим функциональным оператором. Такой подход ставит своей целью посредством построения некоторой модели установить изоморфизм не с внутренней структурой и ее функционированием, а с внешними проявлениями ее информативных параметров. Недостатком метода является то, что АРПСС-модель вытекающая из метода Бокса – Дженкинса не основывается на физических представлениях.
Метод Бокса – Дженкинса основывается на том, что гладкий нестационарный временной ряд путем взятия разностей некоторого d-го порядка можно свести к эквивалентному стационарному, то есть к случаю, для которого разработаны методы анализа и прогнозирования. В методе Бокса – Дженкинса нестационарный временной ряд представляется в виде модели авторегрессии проинтегрированного скользящего среднего (АРПСС, ARIMA) в прямой и возвратной форме:
прямая форма: φ(B)∇d Xt = θ(B)at, t = +1, +2, …,
возвратная форма: φ(F)∇d Xt = θ(F)et, t = –1, –2, …,
φ(B) = 1 – φ1B – φ2 B2 – ... – φP BP = 1 – ∑j=1p φj B j = АР-оператор,
θ(B) = 1 – θ1 B – θ2 B2 – ... – θq Bq=1 – ∑i=1q θi Bi = CC-оператор,
где φ(B), φ(F) – операторы авторегрессии (АР); B, F – операторы сдвига; ∇ – оператор разности: ∇1Xt = Xt–Xt–1 = (1 – B)Xt; d – порядок разности, обеспечивающий переход от нестационарного ряда к эквивалентному стационарному; p – порядок авторегрессии; q, θ(B), θ(F) – порядок и операторы проинтегрированного скользящего среднего (СС) соот-
ветственно; at, et – последовательности независимых случайных величин, имеющих одинаковое нормальное распределение (белый шум): E[at] = 0, Var[at] = σat. Значения белого шума at, согласно методу Бокса – Дженкинса, должны быть получены прогнозированием назад с использованием возвратной формы временного ряда.
Смешанный процесс АР и СС представляет собой эффективную модель для описания стационарных и нестационарных временных рядов. АРПСС-модель эффективна, когда временной ряд является суммой независимых составляющих, каждая из которых описывается либо моделью АР, либо моделью СС, но которые непосредственно не наблюдаются. Модель является параметрической. Чтобы подобрать такую модель, нужно оценить по наблюдаемым данным набор параметров АРПСС-модели φ, θ, σa.
Согласно авторегрессионной модели формирования сигнала, он формируется путем пропускания дискретного белого шума через рекурсивный фильтр, в модели скользящего среднего (СС-модель) используется нерекурсивный фильтр, а в АРПСС-модели – фильтр общего вида, содержащий рекурсивную и нерекурсивную ветви. Фильтры АРПСС имеют обратную связь, обладают большой гибкостью, экономичны, то есть малое число параметров в левых и правых частях.
Заключение
В статье кратко рассмотрены различные подходы, которые возможно использовать для прогнозирования долговечности ИС по параметрическим отказам, базирующие-
ся на методах, заимствованных из различных областей науки и техники, таких как статистика и теория вероятностей, теория распознавания образов и теория ка-
тастроф, теория временных рядов, термодинамика и др. С математической точки зрения применение этих методов для прогнозирования долговечности ИС по парамет-
рическим отказам не имеют ограничений.
Следует отметить, что в настоящий момент ни одним из рассмотренных методов прогнозирования долговечности ИС по параметрическим отказам не получены сведения о фактической долговечности ИС в условиях эксплуатации РЭА. Однако в современной работе [19] проведена верификация физико-статистического подхода теории надежности для простейших случаев, показавшая его правомочность.
Обзор литературных источников показывает, что наиболее отработанными с точки зрения математического аппарата и программного обеспечения являются методы, основанные на теории временных рядов и позволяющие по траекториям процесса деградации параметров ИС прогнозировать долговечность ИС.
ЛИТЕРАТУРА
Ушаков И.А. Вероятностные модели надежности
информационно-вычислительных систем. М.:
Радио и связь, 1991. 132 с.
Барлоу Р., Прошан Ф. Математическая теория
надежности / Пер. с англ.; под ред. Б.В. Гнеденко.
М.: Сов. радио, 1969. 537 с.
Барлоу Р., Прошан Ф. Статистическая теория
надежности и испытания на безотказность /
Пер. с англ. И.А. Ушакова. М.: Наука, 1984. 452 с.
Строгонов А.В. Характеристики надежности
современных ПЛИС // ЭЛЕКТРОНИКА: Наука,
Технология, Бизнес. 2019. № 4. С. 52–58.
Соколов В.И. Проблемы микроэлектроники
(1. Диффузия. 2. Дефектообразование. 3. Деградация) // Физика и техника полупроводников. 1995. Т. 29. С. 842.
Горин Б.М., Кив А.Е, Плотникова Л.Г. и др.
Механизмы естественного старения и вынужденной
деградации полупроводниковых приборов // Обзоры
по электронной техники, 1983. 57 с.
Серия «Полупроводниковые приборы».
Вавилов В.С., Горин Б.М., Данилин Н.С. и др.
Радиационные методы в твердотельной электронике.
М.: Радио и связь, 1990. 184 с.
Алексанян И.Т., Черняев Н.В. Обобщенная модель
надежности изделий // Известия вузов. Электроника. 1998. № 1. C. 85–90.
Алексанян И.Т., Черняев Н.В. Выражения для основных количественных показателей надежности в физико-
статистическом подходе // Петербургский журнал
электроники. 1994. № 1. C. 56–58.
Алексанян И.Т., Черняев Н.В.
Метод изучения надежности интегральных микросхем // Микроэлектроника. 1992. Вып. 2. Т. 21. C. 105–111.
Алексанян И.Т., Черняев Н.В. Управление надежностью интегральных микросхем на основе информационной
избыточности // Известия вузов. Электроника. 1998. № 7. C. 62–66.
Малков Я.В. Диагностика деградационных процессов
в полупроводниковых структурах с позиций междисциплинарных теорий // Зарубежная радиоэлектроника. 1996. № 5. C. 77–82.
Дубицкий Л.Г. Предвестники отказов в изделиях
электронной техники. М.: Радио и связь, 1989. 96 с.
Сыноров В.Ф., Пивоварова Р.П. Параметрическая
надежность и физические модели отказов интегральных схем. Воронеж: Изд-во ВГУ, 1983. 152 с.
Погребинский С.Б., Стрельников В.П.
Проектирование и надежность многопроцессорных ЭВМ. М.: Радио и связь, 1988. 168 с.
Стрельников В.П. Вероятностно-физические методы
исследования надежности машин и аппаратуры //
Надежность и контроль качества. 1989. № 9. С. 3–7.
Алексанян И.Т., Кривошапко В.М.,
Романов А.А. Построение модели параметрической
надежности ИС по данным о деградации
их характеристик и имитация испытаний на ЭВМ //
Электронная техника. 1979. Вып. 1. С. 15–19.
Серия «Микроэлектроника».
Алексанян И.Т., Кривошапко В.М. Моделирование
параметрических отказов и изучение надежности
интегральных схем // Электронная техника. 1981. Вып. 4. С. 52–57. Серия «Управление качеством, стандартизация, метрология, испытания».
Кожевников В.С., Матюшкинa И.В., Черняев Н.В.
Анализ основного уравнения физико-статистического подхода теории надежности технических систем //
Математические основы и численные методы
моделирования. 2020. Т. 12. № 4. С. 721–735.
Отзывы читателей