CC BY
371Научно-практический журнал «Новые исследования в разработке техники и технологий» № 2/2015
Никольский Ю.В.
ПРОГНОЗИРОВАНИЕ НАДЕЖНОСТИ РАДИОЭЛЕКТРОННЫХ УСТРОЙСТВ НА ОСНОВЕ ТЕПЛОФИЗИЧЕСКОГО МОДЕЛИРОВАНИЯ
Современные представления о методах прогнозирования (расчета) надежности как отечественных, так и зарубежных радотехтических устройств (РТУ)отражены в [1]. Там же приведены рассчитанные уровни надежности ЭРИ зарубежного производства. Согласно [1], надежность является функцией многих аргументов (факторов, влияющих на надежность РТУ), в общем случае индивидуальных для каждого класса РТУ. Одним из наиболее важных факторов, определяющих надежность РТУ является температура электрорадиоизделий, а особенно полупроводниковых приборов, так, например, достаточно давно известно, что повышение температуры даже до относительно умеренных значений +40-60°С может приводить [2,3] к увеличению прямых и обратных токов полупроводниковых переходов, повышению проводимости утечки в полупроводниковых приборах и другим эффектам, непосредственно оказывающим влияние на надежность РЭА. Помимо этого, в [1] отмечен факт резкого увеличения в последние годы количества отказов интегральных микросхем из-за дефектов кристалла и корпуса. А дефекты такого рода в значительной степени обусловлены температурной неоднородностью и температурными напряжениями [4].
Для анализа температурных полей в элементах и блоках радиоэлектронной аппаратуры (РЭА) до настоящего времени наиболее часто использовались модели с сосредоточенными параметрами [5] или “нуль-мерные”. Такой подход позволяет существенно упростить процедуру анализа и расчета температурных полей благодаря положенному в основу базовому допущению о том, что как в малоразмерных ЭРИ, так и в достаточно крупных деталях аппаратуры отсутствуют градиенты температуры по всем координатным направлениям. Соответственно температурное поле любого ЭРИ или блока РЭА является однородным и характеризуется некоторой средней температурой, которая зависит только от времени. Поэтому в этом случае достаточно реально может быть отражено тепловое состояние лишь узлов радиоэлектронной аппаратуры, в которых неравномерность температурного поля невелика - порядка единиц градусов. В других случаях обычно остается открытым вопрос о погрешности
23
ISSN 2313-1160
такой интерпретации реальных температурных полей.
Известны также многомерные математические модели, описывающие процессы теплопереноса в РЭА в рамках двумерных и трехмерных постановок [6]. Но в настоящее время отсутствуют такие их реализации, в которых одновременно сочетались бы следующие важные признаки: пространственные области соответствуют реальным узлам и (или) элементам РЭА (например, интегральным схемам); учитываются механизмы конвективного и радиационного теплообмена с внешней средой и режим работы РЭА (нестационарность). Модель, отвечающая всем этим требованиям, требует вычислительных ресурсов, практически недоступных до недавнего времени для исследователей. Важным является и то, что реальная аппаратура космического, авиационного, морского, транспортного базирования в основном работает в нестационарных режимах, обусловленных спецификой ее использования. Поэтому известные аналитические и численно-аналитические методы [6-8], достаточно эффективные при решении сложных стационарных линейных и даже некоторых нелинейных задач теплопроводности не всегда могут применяться при анализе тепловых режимов реальной радиоэлектроники. Целью данной работы является анализ
надежности элементов РТУ с применением современного аппарата теории надежности [1] и методов математического моделирования температурных полей. Достижение поставленной цели планируется путем численного моделирования нестационарных полей температур в типичном элементе (печатном узле с радиоэлементами) по пространственной теплофизической математической модели с учетом реальных механизмов теплообмена с внешней средой и последующего сопоставления показателей надежности типичных ЭРИ, полученных на базе пространственных температурных полей и средних по объему температур ЭРИ.
Постановка задачи.
Решается задача расчета температурного поля в параллелепипеде, размеры которого ограничены по осям х, у и z величинами, равными Lx, Ly и Lz. Параллелепипед включает в себя несколько элементов (также параллелепипедов) с отличающимися теплофизическими характеристиками и размерами. На границах между элементами реализуется идеальный тепловой контакт. Внутри параллелепипеда действует локальный источник тепловыделения заданной интенсивности Q. На границах заданы граничные условия III рода с излучением (смешанный теплообмен).
Геометрия области решения представлена на рисунке 1. Реальная конструкция представляет собой интегральную микросхему (ИС) стабилизатора напряжения, расположенную на двухслойной печатной плате из стеклотекстолита толщиной 1,5 мм.
Представленная на рисунке 1 модель ЭРИ была выбрана не только как одна из наиболее типичных реализаций современных электрорадиоэлементов [6]. По данному объекту удалось получить экспериментально значения температур в
24
Научно-практический журнал «Новые исследования в разработке техники и технологий» № 2/2015
ряде точек рассматриваемой модели при контролируемых условиях теплообмена с окружающей средой в типичных режимах работы. Погрешность определения температуры Т при этом не превышала ±6% для достаточно большого (tmax =180с) времени работы изделия.
Задачавтакойпостановкесводитсякрешениюнелинейногонестационарного уравнения теплопроводности с соответствующими граничными и начальными условиями:
^ ч , ,сТ д ^ ,DT \ д(„ ,дГ } д ( ^ ,DT )
С (х, У, z)P(x, У,z)— = —I Л(х У, z)— I + —I Л(х, y, z)— I + —I Л(х, y, z)— 1 + Q(t, х, y, z)
ot дх ^ дх J ду ^ ду ) дz ^ дz J
t е [о; tmax 1 х е [О; Ьх ] у е [0; Ly ] z е [0; Lz ]
Tt=0 = T0( х у, z);
х = 0; y е[0; Ly ] z е[0; Ь ] дТ -Л— = а(Тв -Т) + snpa(TB -Т4); дх
х = Ьх; у е [°; Ly ]z е [0; Ь] дТ Л— = а(Тв -Т) + *^(TB -Т4); дх
y = 0;х е[0; Lx ] z е[0; L ] дТ -Л =а(Тв -Т) + 8^^ -Т4); ду
у = Ly;х е [0; L1z е [0; Lz] дТ Л— = а(Тв -Т) + 8ПРа(Тв4 -Т4); ду
z = 0; у е [0; Ly ]х е [0; Lx ] дТ -Л— = а(Тв - Т) + 8ПР&Т - Т4); дz
z = Lz; у е [0; Ly ]х е [0; L ] дТ Л— = а(Тв -Т) + 8ПР&Т -Т4); дz
где С - удельная теплоемкость; х, у, z - координаты; р - плотность; Т - температура; t - время; X - коэффициент теплопроводности; Q ~ удельное тепловыделение источника; tmax - время работы; Lx, Ly ,Lz - размеры параллелепипеда по осям х, у и z соответственно; T0 - начальная температура; а -коэффициент теплоотдачи; Тв - температура внешней среды; е пр - приведенный коэффициент черноты тела и окружающей среды; о - постоянная Стефана-Больцмана.
Сформулированная краевая задача решена методом конечных разностей [9]. Разностный аналог исходного нелинейного пространственного дифференциального уравнения решен методом расщепления по координатам. Одномерные разностные уравнения решены методом итераций с применением метода прогонки на каждой итерации.
25
ISSN 2313-1160
Модель с сосредоточенными параметрами
В основу моделей с сосредоточенными параметрами положено допущение о том, что градиенты температуры по всем координатным направлениям отсутствуют [5,6]. Таким образом, параллелепипед заменяется однородным телом, теплофизические характеристики которого вычисляются путем объемного осреднения теплофизических характеристик элементов, составляющих параллелепипед. Для однородного тела запишем уравнение теплового баланса [6]:
dT a TB - T) snM TB - T4) Q
dt Lz Lz Q
t G I0 tmax ] Tt=0 = T0 •
где С - удельная теплоемкость; р - плотность; Т - температура; t - время; а - коэффициент теплоотдачи; Тв - температура внешней среды; Lz - размеры тела по оси z; ; е пр - приведенный коэффициент черноты тела и окружающей среды; о - постоянная Стефана-Больцмана; Q ~ удельное тепловыделение источника; tmax - время работы; Т0 - начальная температура.
Решается это дифференциальное уравнение методом Рунге-Кутта четвертого порядка точности.
Модель расчета надежности ЭРИ.
Согласно [1] современные методы анализа и прогнозирования надежности ЭРИ опираются на определенные математические модели, в основе которых лежит предположение об экспоненциальном распределении числа наработок до отказа и параметр распределения является функцией факторов испытаний. При этом в качестве важнейшего фактора всегда принимается температура Т, а типичной моделью является модель Аррениуса:
_ E
AT) = С • exp(—) (1)
где С-константа, Е -энергия активации, k -постоянная Больцмапа.
Вычисленное по модели с сосредоточенными параметрами значение температуры Т, является среднеобъемной для моделируемого узла РЭА, поэтому значение числа наработок до отказа А, вычисленное по (1), также будет средним по объему.
Пространственная модель позволяет определять значение температуры в тех местах моделируемого узла РЭА (для тех элементов), которые фактически определяют надежность узла РЭА. Поэтому вычисленное с помощью уравнения
26
Научно-практический журнал «Новые исследования в разработке техники и технологий» № 2/2015
(1) по пространственной модели значение числа наработок до отказа X будет реальнее отражать действительное состояние надежности узла РЭА.
Результаты и обсуждение.
Основные численные и экспериментальные исследования проведены при следующих значениях исходных данных и параметров: Lx = 60 мм, Ly = 55 мм, Lz= 6 мм.
Температура окружающей среды Тв= 300 К. Мощность источника 1,6 и 2,1 Вт. Теплофизические характеристики областей 1-6 (см. рисунок 1) взяты из справочных данных [5].
Геометрия области решения: 1) плата; 2) фольга; 3) корпус ИС; 4) теплоотвод ИС; 5) кристалл ИС; 6) выводы ИС;
7) точка измерения температуры термопарой.
Для расчетов по пространственной модели использовалась разностная сетка размерами 301x276x31, шаг по времени 0,06 с.
Коэффициент конвективного теплообмена с внешней средой в обеих моделях рассчитывался по методикам [6] и составлял от 0 до 2,9 Вт/(м2-К).
Для анализа результатов численного эксперимента по пространственной модели выделена температура в фиксированной точке корпуса ИС Тр.
Для проверки точности полученных результатов были проведены эксперименты. В фиксированной точке верхней поверхности корпуса ИС термопарой измерялась температура Тэ соответствующая Тр для численного исследования. Для проверки точности полученных результатов были проведены эксперименты. В фиксированной точке верхней поверхности корпуса ИС термопарой измерялась температура Тэ соответствующая Тр для численного исследования. Результаты расчетов и измерений приведены в таблице 1.
27
ISSN 2313-1160
Таблица 1 - Измеренная и вычисленная по пространственной модели температуры корпуса
Время, с 0 20 40 60 80 100 120 140 160 180
Мощность источника 1 ,6 Вт
Измеренная температура корпуса Тэ, °С 29 38 47 54 59 66 69 74 76 78
Расчетная температура корпуса Тр , °С 29 40 48 55 60 65 69 72 76 78
Мощность источника 2,1 Вт
Измеренная температура корпуса Тэ , °С 27 42 52 61 68 76 81 85 88 92
Расчетная температура корпуса Тр , °С 27 42 53 61 68 74 79 84 88 92
Для проверки точности полученных результатов были проведены эксперименты. В фиксированной точке верхней поверхности корпуса ИС термопарой измерялась температура Тэ соответствующая Тр для численного исследования. Результаты измерений температуры представлены в таблице 1.
Сопоставление результатов численного моделирования по пространственной модели с экспериментальными данными показало, что их расхождение не превышает ±2°С при абсолютных значениях Т до 92°С. Т.е. отклонение теоретических результатов от опытных данных не превышает ±3% или (другими словами) укладывается в диапазон погрешности средств измерения (термопар).
Как показал анализ результатов исследований, с ростом времени разность T = ТС-ТЭ достигает 46°С или 50% от абсолютного значения Тэ. Т.е. даже для такой достаточно простой реализации блока РЭА, работающего в условиях интенсивного теплообмена с внешней средой имеет место существенная неоднородность температурного поля, обусловленная наличием в области решения сильных неоднородностей теплофизических характеристик и источников тепловыделения.
Моделирование показало, что даже при умеренных тепловых нагрузках (относительно малой мощности тепловыделения источника) температурные поля характеризуются большими градиентами [11]. Перепад температуры Т максимально достигает в пределах платы 65°С. Естественно, эти перепады температуры не фиксируются моделью с сосредоточенными параметрами.
Для интегральных микросхем, изготовленных по биполярной и МОП (металл-окисел-полупроводник) технологиям по данным [1] с использованием выражения (1) были вычислены показатели интенсивности отказов для температуры в фиксированной точке корпуса Тр, полученной в многомерной постановке, и средней температуры Тс, полученной по модели с сосредоточенными
28
Научно-практический журнал «Новые исследования в разработке техники и технологий» № 2/2015
параметрами.
Как и следовало ожидать, учитывая экспоненциальныйхарактерзависимости Х(Т), отклонение температуры на 46°С (или на 50% от максимальных значений) приводит к падению показателя надежности для всех типов рассмотренных ИС минимально в 5 раз, а для некоторых (полупроводниковые линейные ИС) в 20 раз и более. Полученный результат позволяет сделать вывод о том, насколько существенны изменения в показателях надежности при переходе к анализу надежности с учетом пространственной неоднородности температурных полей в реальных ЭРИ и деталях РЭА.
Специально исследовалось влияние нестационарности в режиме работы на показатели надежности РЭА. На рисунке 2 приведены зависимости расчетного числа отказов от времени для нестационарного теплового режима работы печатной платы.
Изменение во времени интенсивности отказов интегральной микросхемы, вычисленное для температуры фиксированной точки (1) и средней температуры (2) для источника мощностью 2,1 Вт.
Хорошо видно, что с ростом времени (и соответственно температуры) величина X растет достаточно существенно для выделенной точки пространственной модели, в то время как рост величины для средней температуры обладает существенно меньшей скоростью (как и рост температуры).
Необходимо отметить, что нестационарность теплового режима объекта исследования была обусловлена, как это обычно бывает, не нестационарностью режима работы РЭА, а поглощением тепла за счет теплоемкости материалов печатной платы до выхода на стационарный режим теплообмена с внешней
29
ISSN 2313-1160
средой. Аналогичные закономерности проявляются при работе большинства радиотехнических систем в режимах “включение-выключение”. Так, например, выход на стационарный тепловой режим приборов бортовой аппаратуры модуля полезной нагрузки спутника связи в негерметичном исполнении составляет от 16 до 22 часов [10]. Поэтому, вероятно, существуют такие режимы работы бортовой радиоэлектроники, когда возможно применение стационарных моделей теплопроводности, но для РЭА космического, авиационного, морского, транспортного базирования нестационарные тепловые режимы являются достаточно типичными. Так, например, бортовая электроника космических аппаратов типа спутников связи, теле и радиовещания работает в основном в нестационарных режимах. Для ряда систем бортовой аппаратуры спутников связи и контроля стационарные режимы просто недостижимы.
Полученные результаты в целом хорошо соответствуют основным физическим представлениям как о тепловом состоянии элементов РЭА, так и о надежности радиоэлектронной аппаратуры [1].
На основании проведенных исследований можно сделать вывод, что учет неоднородности температурного поля типичных элементов и блоков радиоэлектронной аппаратуры приводит к существенному (от 5 до 20 и более раз) изменению показателей надежности РЭА, рассчитываемых с использованием современных математических моделей надежности [1]. Это позволяет утверждать о перспективности применения пространственных теплофизических моделей узлов и деталей РЭА на этапе ее конструирования для прогнозирования надежности радиоаппаратуры. Для дальнейшего повышения точности прогностического моделирования необходимо более точное отражение внутренней структуры сложных радиоэлементов, например, интегральных микросхем с большой функциональной насыщенностью и плотностью компоновки (микропроцессоры, большие и сверхбольшие интегральные схемы, оперативные запоминающие устройства).
Литература
1. Борисов А. А,, Горбачева В. М., Карташов Г. Д., Мартынова М. Н., Прыткое С. Ф. Надежность зарубежной элементной базы. // Зарубежная радиоэлектроника. 2000, № 5,с.34-53.
2. Бердичевский Б.Е. Вопросы обеспечения надежности радиоэлектронной аппаратуры при разработке. М.: Советское радио. 1977. 384 с.
30
Научно-практический журнал «Новые исследования в разработке техники и технологий» № 2/2015
3. Маджарова Т.Б., Гаврилова Г.А. Изменение параметров
полупроводниковых линейных интегральных микросхем в зависимости от температуры воздуха //Электронная техника. Серия 8. 1975. Вып.1(31), с.21-32.
4. Алексеев В.П., Кузнецов Г.В., Рудзинский В.П. Теплофизические аспекты надежности элементов радиоэлектронной аппаратуры. Труды второй Российской национальной конференции по теплообмену 1998, Т.1, с.117-120.
5. Дульнев Г.Н. Тепло- и массобмен в радиоэлектронной аппаратуре. М.: Высш.шк., 1984.247с.
6. Дульнев Г.Н., Парфенов В.Г., Сиголов А.В. Методы расчета теплового режима приборов. М.: Радио и связь, 1990.
7. Котляр Я.М., Совершенный В.Д., Стриженов Д.С. Методы и задачи тепломассообмена. М.: Машиностроение, 1987.
8. Карташов Э.М. Аналитические методы в теории теплопроводности твердых тел. М.:Высшая школа, 2001.
9. Самарский А.А. Теория разностных схем. М.: Наука. 1983. 616 с.
10. Ашурков Е.А., Бураков В.А., Козлов А.Г., Кожухов В.П., Корчагин Е.Н., Кузнецов Г.В., Санду С.Ф. Математическое моделирование нестационарных теплофизических процессов в отсеках бортовой аппаратуры космических аппаратов. // Известия вузов.Физика. 1993, №4, с. 119-128.
11. Алексеев В.П., Кузнецов Г.В., Шлома С.В. О теплофизическом моделировании надежности элементов радиоэлектронной аппаратуры // Тезисы докладов ХХVI Сибирского теплофизического семинара. Новосибирск: Изд-во ИТ СО РАН, 2002. с.9-10.
31
