Научная статья на тему 'ОЦЕНКА АДЕКВАТНОСТИ ТЕПЛОВОЙ МОДЕЛИ КМОП ЦИФРОВЫХ ИНТЕГРАЛЬНЫХ СХЕМ ПО ПЕРЕХОДНЫМ ТЕПЛОВЫМ ХАРАКТЕРИСТИКАМ'

ОЦЕНКА АДЕКВАТНОСТИ ТЕПЛОВОЙ МОДЕЛИ КМОП ЦИФРОВЫХ ИНТЕГРАЛЬНЫХ СХЕМ ПО ПЕРЕХОДНЫМ ТЕПЛОВЫМ ХАРАКТЕРИСТИКАМ Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
22
23
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЦИФРОВЫЕ ИНТЕГРАЛЬНЫЕ СХЕМЫ / ТЕПЛОВАЯ СХЕМА / ТЕПЛОВЫЕ ПАРАМЕТРЫ / ПЕРЕХОДНАЯ ТЕПЛОВАЯ ХАРАКТЕРИСТИКА / ИЗМЕРЕНИЕ / АЛГОРИТМ РАСЧЕТА

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Сергеев Вячеслав Андреевич, Тетенькин Ярослав Геннадьевич

Эффективным диагностическим методом контроля качества цифровых интегральных схем (ЦИС) является измерение их тепловых параметров. Значения тепловых параметров реальных изделий определяются качеством их производства и могут существенно отличаться от расчетных значений. Описан способ измерения переходной тепловой характеристики (ПТХ) ЦИС с использованием температурной зависимости частоты кольцевого генератора, построенного на логических элементах ЦИС. На примере ПТХ ЦИС типа CD4011 показан итерационный алгоритм расчета параметров линейной одномерной тепловой схемы ЦИС в приближении Фостера по значениям ПТХ в точках наименьшей крутизны, которые соответствуют нулям второй производной. Алгоритм апробирован при нахождении нулей второй производной ПТХ путем прямого численного дифференцирования ПТХ и дифференцирования аппроксимирующего ПТХ полинома 9-й степени. Показано, что ПТХ исследуемой ЦИС соответствует трехзвенной тепловой схеме. Рассмотренные варианты расчета дают практически одинаковые значения тепловых параметров тепловой схемы, которые достаточно хорошо идентифицируются с характерными слоями конструкции ЦИС.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Сергеев Вячеслав Андреевич, Тетенькин Ярослав Геннадьевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ASSESSMENT OF ADEQUACY OF THE CMOS LINEAR THERMAL MODEL OF DIGITAL INTEGRATED CIRCUITS ON TRANSIENT THERMAL CHARACTERISTICS

Effective diagnostic method of quality control of the digital integrated circuits (DIC) is measurement of their thermal parameters. Values of thermal parameters of real products are defined by quality of their production and can significantly differ from calculated values. A method for measuring the transient thermal characteristics (TTC) of the digital integrated circuits (DIC) with the temperature dependence of the frequency of the DIC ring oscillator has been described. An iterative algorithm for calculating parameters of the linear one-dimensional Foster thermal circuit uses the values of TTC, which correspond to the zeros of the TTC second derivative, has been shown. The algorithm was tested in determining the zeros of the second derivative of TTC by direct numerical differentiation method and by differentiating the 9-th order polynomial approximating function. It has been shown that TTC of the tested DIC corresponds to the three-stage thermal circuit. Both methods of calculation give almost identical values of thermal parameters of the thermal scheme which well corresponds to different technological layers of the DIC structure.

Текст научной работы на тему «ОЦЕНКА АДЕКВАТНОСТИ ТЕПЛОВОЙ МОДЕЛИ КМОП ЦИФРОВЫХ ИНТЕГРАЛЬНЫХ СХЕМ ПО ПЕРЕХОДНЫМ ТЕПЛОВЫМ ХАРАКТЕРИСТИКАМ»

ЭЛЕМЕНТЫ ИНТЕГРАЛЬНОЙ ЭЛЕКТРОНИКИ INTEGRATED ELECTRONICS ELEMENTS

УДК 621.382.8.017.7 DOI: 10.24151/1561-5405-2017-22-4-350-360

Оценка адекватности тепловой модели КМОП цифровых интегральных схем по переходным тепловым характеристикам

В.А. Сергеев, Я.Г. Тетенькин

Ульяновский филиал Института радиотехники и электроники им. В.А. Котельникова Российской академии наук, г. Ульяновск, Россия Ульяновский государственный технический университет, г. Ульяновск, Россия

sva@ulstu.ru

Эффективным диагностическим методом контроля качества цифровых интегральных схем (ЦИС) является измерение их тепловых параметров. Значения тепловых параметров реальных изделий определяются качеством их производства и могут существенно отличаться от расчетных значений.

Описан способ измерения переходной тепловой характеристики (ПТХ) ЦИС с использованием температурной зависимости частоты кольцевого генератора, построенного на логических элементах ЦИС. На примере ПТХ ЦИС типа СБ4011 показан итерационный алгоритм расчета параметров линейной одномерной тепловой схемы ЦИС в приближении Фостера по значениям ПТХ в точках наименьшей крутизны, которые соответствуют нулям второй производной. Алгоритм апробирован при нахождении нулей второй производной ПТХ путем прямого численного дифференцирования ПТХ и дифференцирования аппроксимирующего ПТХ полинома 9-й степени. Показано, что ПТХ исследуемой ЦИС соответствует трехзвенной тепловой схеме.

Рассмотренные варианты расчета дают практически одинаковые значения тепловых параметров тепловой схемы, которые достаточно хорошо идентифицируются с характерными слоями конструкции ЦИС.

Ключевые слова: цифровые интегральные схемы; тепловая схема; тепловые параметры; переходная тепловая характеристика; измерение; алгоритм расчета.

Для цитирования: Сергеев В.А., Тетенькин Я.Г. Оценка адекватности тепловой модели КМОП цифровых интегральных схем по переходным тепловым характеристикам // Изв. вузов. Электроника. - 2017. - Т. 22. - № 4. - С. 350-360. БО!: 10.214151/1561-5405-2017-22-4-350-360

© В.А. Сергеев, Я.Г. Тетенькин, 2017

Assessment of Adequacy of the CMOS Linear Thermal Model of Digital Integrated Circuits on Transient Thermal Characteristics

V.A. Sergeev, Ya.G. Tetenkin

Ulyanovsk Branch of Kotel'nikov Institute of Radio-Engineering and Electronicsof Russian Academy of Sciences, Ulyanovsk, Russia Ulyanovsk State Technical University, Ulyanovsk, Russia

sva@ulstu.ru

Effective diagnostic method of quality control of the digital integrated circuits (DIC) is measurement of their thermal parameters. Values of thermal parameters of real products are defined by quality of their production and can significantly differ from calculated values.

A method for measuring the transient thermal characteristics (TTC) of the digital integrated circuits (DIC) with the temperature dependence of the frequency of the DIC ring oscillator has been described. An iterative algorithm for calculating parameters of the linear one-dimensional Foster thermal circuit uses the values of TTC, which correspond to the zeros of the TTC second derivative, has been shown. The algorithm was tested in determining the zeros of the second derivative of TTC by direct numerical differentiation method and by differentiating the 9-th order polynomial approximating function. It has been shown that TTC of the tested DIC corresponds to the three-stage thermal circuit.

Both methods of calculation give almost identical values of thermal parameters of the thermal scheme which well corresponds to different technological layers of the DIC structure.

Keywords: digital integrated circuits; thermal equivalent schemes; thermal parameters; transient thermal characteristics; measurement; algorithm of calculation.

For citation: Sergeev V.A., Tetenkin Ya.G. Assessment of Adequacy of CMOS Linear Thermal Model of Digital Integrated Circuits on Transient Thermal Characteristics // Proc. of Universities. Electronics. - 2017. - Vol. 22. - No. 4. - P. 350-360. DOI: 10.214151/1561-5405-2017-22-4-350-360

Введение. Функциональные возможности и физическая надежность цифровых интегральных схем (ЦИС) в эксплуатационных режимах во многом определяются температурой активной области. Тепловые свойства ЦИС закладываются при проектировании. Однако значения тепловых параметров реальных изделий определяются качеством их производства и могут существенно отличаться от расчетных значений. Поэтому эффективным диагностическим методом контроля качества ЦИС является измерение их тепловых параметров [1, 2].

Согласно принципу теплоэлектрической аналогии [1, 2], тепловая схема ЦИС пла-нарной конструкции с односторонним отводом тепла, в которой тепловой поток направлен от активной области на рабочей поверхности чипа к корпусу прибора и далее в окружающую среду, в приближении Фостера представляет собой m последовательно соединенных RC-звеньев (рис.1). Задача контроля тепловых свойств ЦИС сводится к определению значений тепловых сопротивлений RT и теплоемкостей CT или тепловых постоянных времени (tt = RtCt) RC-звеньев тепловой схемы, соответствующих определенным слоям конструкции ЦИС.

Рис.1. Схематическое изображение конструкции ЦИС (а) и соответствующая ей тепловая схема (б) в приближении Фостера (0п - температура перехода, 0о - температура окружающей среды)

Fig.1. A design of DIC schematically (a) it is schematical, and the thermal scheme (b) corresponding to it in Foster's approach (0п - transition temperature, 0о - ambient temperature)

Способы измерения переходной тепловой характеристики ЦИС. Тепловые параметры ЦИС могут быть определены по переходной тепловой характеристике (ПТХ), которая получена по изменению температуры Д0П (t) активной области ЦИС после подачи на нее ступеньки греющей мощности. Для тепловой схемы, состоящей из m .RC-звеньев, зависимость Д0П (t) описывается выражением [3]

Д0п (t) = PTZT (t) = PT ]Г RTi (l - exp (t/It )), (1)

i=i

где PT - тепловая мощность, рассеиваемая ЦИС; ZT (t) = ^ R^ (l - exp (t/xT, )) - переход-

i=i ' '

ное тепловое сопротивление, или ПТХ ЦИС.

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

В зарубежной практике для измерения ПТХ полупроводниковых изделий с р-и-переходами используется способ по стандарту JESD51-14 [3], состоящий в том, что на полупроводниковое изделие подают ступеньку электрической греющей мощности заданной величины. В процессе разогрева изделия в определенные моменты времени tj на короткое время (длительностью до нескольких микросекунд) греющую мощность отключают, через р-и-переход пропускают малый прямой ток и измеряют температу-рочувствительный параметр - напряжение U n(t) нар-и-переходе c известным температурным коэффициентом Кц. Приращение температуры Д0п (tj ) р-и-перехода в момент времени tj рассчитывают по формуле

Д0 п (tj ) =

Up-n (tj) - Up-n (0)

K

и

где Up-n (0) - падение напряжения на р-и-переходе до разогрева изделия; U p-n (tj) -

p-n\ j,

падение напряжения на р-и-переходе в момент времени tj.

Точность измерения ПТХ определяется числом отсчетов на декаду временной оси, и в известном измерителе Т3Б1ег [4], реализующем описанный способ, число отсчетов достигает 20 точек на декаду временной оси.

Недостатком этого способа является значительная погрешность измерения напряжения на p-n-переходе сразу после выключения греющей мощности из-за влияния паразитных переходных электрических процессов при переключении полупроводникового изделия из греющего режима в измерительный режим [5]. Для снижения этой погрешности температурочувствительный параметр необходимо измерять через некоторое время после выключения греющей мощности, за которое электрический переходный процесс в основном завершится. За это время температура p-n-перехода может заметно измениться. Для исключения этой погрешности измеряется кривая охлаждения полупроводникового изделия, т. е. изменение температуры активной области после нагрева до стационарного теплового режима и выключения греющей мощности. При этом время измерения ПТХ увеличивается как минимум в два раза. Следует учитывать, что кривые нагрева и охлаждения полупроводникового изделия не являются строго симметричными. Кроме того, в изделиях со сложной топологией активной области, к которым относятся ЦИС, источники тепла и датчик температуры (^^переход), как правило, пространственно разнесены, что приводит к заниженному значению температуры. Применительно к ЦИС в качестве датчика температуры используется один из защитных диодов по входам логических элементов либо диод, специально сформированный в структуре ЦИС.

Измерение переходной тепловой характеристики ЦИС по температурной зависимости частоты кольцевого генератора. Применительно к КМОП ЦИС предлагается легко автоматизируемый способ [6, 7], лишенный описанных недостатков. Разработан аппаратно-программный комплекс [8] для измерения ПТХ по изменению частоты кольцевого генератора, построенного на логических элементах контролируемой ЦИС. Суть способа состоит в том, что в заданные моменты времени ^ после включения кольцевого генератора измеряются мгновенная мощность, потребляемая ЦИС от источника питания, и частота колебаний кольцевого генератора. Значение ПТХ в момент времени ^ определяется по формуле

2 ^ ) =

/к.г (0) - /к.г )

где /кг (0) и /кг ) - частота генерации кольцевого генератора в моменты времени

^ = 0 и ^ соответственно; ^- температурный коэффициент частоты колебаний кольцевого генератора; Р (г ■) - средняя мощность, потребляемая ЦИС за время от начала нагрева до момента времени

В основе способа лежат эффекты саморазогрева ЦИС рассеиваемой при работе кольцевого генератора мощностью и уменьшения частоты его генерации с ростом температуры. Относительный температурный коэффициент частоты генерации кольцевого генератора определяется температурной зависимостью времени задержки распространения сигнала логических элементов ЦИС, является практически постоянным в диапазоне 0-100 °С и имеет значение порядка 0,2-0,3 %/°С, т.е. сравним с температурными коэффициентами электрических параметров ЦИС [1-4]. При реализации предлагаемого способа не требуется переключать ЦИС из режима нагрева в режим измерения температурочувствительного параметра, ПТХ измеряется в режиме нагрева, а датчики

температуры совмещены с источниками нагрева. Отметим, что способ также дает некоторое усредненное по логическим элементам ЦИС значение температуры, но это значение адекватно отражает тепловые свойства ЦИС в эксплуатационных режимах.

Алгоритм определения тепловых параметров ЦИС по переходной тепловой характеристике. Для определения тепловых параметров полупроводниковых изделий по кривым нагрева (или охлаждения) в [9] предложено использовать так называемую дифференциальную структурную функцию, которая определяется как производная суммарной теплоемкости по суммарному тепловому сопротивлению K = dCY/dRz слоев полупроводниковых изделий в заданный момент времени по мере распространения тепла в конструкции изделия. Алгоритм расчета структурной функции основан на сложных многоступенчатых преобразованиях при переходе от модели Фостера к модели Кауэра. Это приводит к накоплению систематических ошибок и, как следствие, к потере точности [10, 11].

Более простой алгоритм определения тепловых параметров ЦИС по ПТХ рассмотрим на примере ПТХ ЦИС типа СБ4011, измеренной в спокойном воздухе с помощью аппаратно-программного комплекса по изменению частоты кольцевого генератора с 10 отсчетами на каждую (кроме первой) декаду шкалы времени (табл.1).

Таблица 1

Результаты измерения переходной тепловой характеристики ЦИС типа CD4011

по частоте кольцевого генератора

Table 1

Results of measurement of ТТС of the CD4011 DIC on RG frequency

tj, Zj tj, Z(t;), tj, Z(t;), tj, Zj tj, Z(tj)

мс К/Вт мс К/Вт мс К/Вт мс К/Вт мс К/Вт

0,13 0,524 11,81 7,393 228,07 21,97 4308 42,25 80988 65,16

0,29 1,170 14,67 7,839 281,07 23,66 5313 43,34 99861 68,18

0,48 2,013 18,19 8,484 347,1 25,34 6552 44,41 123132 71,28

0,72 2,757 22,53 9,054 428,1 26,95 8080 45,45 151825 74,55

1,01 3,476 27,88 9,798 528,1 28,52 9964 46,56 186825 77,92

1,37 3,922 34,48 10,591 651,1 30,10 12287 47,68 229825 81,45

1,81 4,418 42,62 11,48 803,1 31,56 15152 49,02 282825 84,77

2,35 4,914 52,65 12,48 991,1 33,03 18685 50,33 348825 87,87

3,01 5,261 65,02 13,57 1223 34,47 23041 51,92 429825 90,57

3,83 5,558 80,27 14,73 1509 35,73 28412 53,68 529825 92,80

4,84 5,856 99,07 16,00 1861 37,04 35034 55,59 653825 94,66

6,08 6,203 122,07 17,41 2296 38,26 43199 57,65

7,61 6,550 150,07 18,90 2832 39,80 53267 59,95

9,49 6,947 185,07 20,43 3493 41,11 65681 62,45

Тепловые постоянные времени слоев конструкции исследуемой ЦИС различаются на несколько порядков величины. На графике ПТХ (рис.2), построенном в полулогарифмическом масштабе, отчетливо наблюдаются три участка быстрого и медленного изменения. Участок быстрого роста соответствует накоплению тепла в соответствующем слое конструкции ЦИС (подобно накоплению заряда в конденсаторе), а участок медленного роста - квазистационарному процессу роста градиента температуры на границе слоев. Число слоев в конструкции ЦИС и, соответственно, число слагаемых в

2(t), К/Вт

—' шш

S / /

/

У 2 Л ' \

------ * % % 3 \ f \

ч. \ »»

0.1 .-■ 1........ 1 Э 1 I2 .......1 б3 '•■/ 1 О4 1 $ / Т, мс

*****

Рис.2. Исходная ПТХ ЦИС типа CD4011 (кривая 1) и ее первая (кривая 2) и вторая (кривая 3) производные, полученные прямым численным дифференцированием

по пяти точкам со сглаживанием Fig.2. Transient thermal characteristic of CD4011 DIC (curve 1) and its first (curve 2) and second (curve 3) derivatives, received by direct numerical

differentiation

правой части выражения (1) определяются числом участков быстрого или медленного роста на графике ПТХ. Очевидно, что точки максимальной и минимальной крутизны ПТХ соответствуют максимуму и минимуму ее первой производной или нулям второй производной. Точки перехода через ноль второй производной снизу вверх соответствуют участку ПТХ с минимальной крутизной, и значение полного теплового сопротивления ZТ{tk} в к-й точке минимального роста ПТХ равно сумме тепловых сопротивлений к первых .КС-звеньев тепловой цепи:

к

Z Т К } = X Кт1 .

7=1

Поскольку ^ {^ } = К^ , то для нахождения теплового сопротивления К^ первого

звена тепловой цепи надо определить первую точку минимальной крутизны ПТХ.

Тепловая постоянная времени к-го звена определяется в точках ПТХ, в которых выполняется условие 1 = тТ/с. Эти точки согласно (1) при выполнении условия тТ<<тТк+1

соответствуют моменту времени, при котором тепловое сопротивление к-го слоя достигает значения, равного (1- е~ ) ~ 0,63 от полного значения КТ :

( \ к-1

2т (' = ч ) = X Кт7 + °,63КТк = 2Т {к-1} + 0,63КТк.

7=1

Данные точки легко находятся на графике ПТХ. При этом первая тепловая постоянная времени соответствует точке временной оси, для которой переходное тепловое сопротивление 2Т ) принимает значение, равное 0,63 от теплового сопротивления Кт

первого звена. Для исследуемого образца получено КТ =5,8 К/Вт, т^ = 1,35 мс.

В результате последовательного выполнения итераций для всех точек минимальной крутизны (третья точка найдена экстраполяцией второй производной) получены значения тепловых параметров для всех трех звеньев. Выражение для функции (1) имеет вид (значения тепловых постоянных времени даны в миллисекундах)

1тх (1) = 5,8(1 - ^(-1/1,35)) + 34,6(1 - exp(-t/480))+58(1 - ^(-1/176000)).

Для проверки возможности повышения точности определения тепловых параметров тепловой схемы ЦИС проведен их расчет с использованием полиномиальной регрессии экспериментальной ПТХ. Путем подбора в среде Mathcad найден полином 9-й степени, аппроксимирующий исходную ПТХ с минимальной средней квадратической ошибкой (СКО):

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ft) = -1,15 + 1,84t - 0,324t2 + 0,0395t3 - 0,0029t4 + 0,00013t5 - 3,55-10-6t6 +

+ 5,4-10-8t7 - 4,375-10-10t8 + 1,43-10-12t9 (2)

Вторая производная аппроксимирующего полинома (2) также имеет три точки перехода через ноль снизу вверх, которые несколько отличаются от точек, полученных прямым численным дифференцированием ПТХ (рис.3).

ZQ), К/Вт _______

f "-V

* \ % % \ 1 и // / V

\ \ \ 3 \ % \

....... % \ \

0,1 1 / 1 0 1 О2"........1 ..... 1 О4 1 К Т, мс

Рис.3. Исходная ПТХ (кривая 1) и ее первая (кривая 2) и вторая (кривая 3) производные аппроксимирующего полинома 9-й степени (для наглядности

первая и вторая производные увеличены соответственно в 20 и 100 раз) Fig.3. Initial ТТС (curve 1) and first (curve 2) and second (curve 3) derivatives of the approximating polynom of the 9th degree (for descriptive reasons the first and second derivatives are increased, respectively in 20 and 100 times)

По итерационному алгоритму определения тепловых параметров получено выражение для функции вида (1), но с другими параметрами:

ZГ2(0 = 5,4(1 - ^(-^1,30)) + 37,8х х(1 - ^(-^520)) + 54,6(1 - ^(-^180000)).

Как видно из сравнения тепловых параметров ЦИС, найденных с помощью итерационного алгоритма по точкам наименьшей крутизны ПТХ разными способами обработки исходной ПТХ, их различие (кроме т^) не превышает 10 % (табл.2).

Графики исходной экспериментальной ПТХ и аппроксимирующих ее функций приведены на рис.4. СКО аппроксимации экспериментальной ПТХ функцией (¿) составила 1,48 К/Вт, функцией 2Тг (¿) - 1,42 К/Вт. Таким образом, оба способа обработки экспериментальной ПТХ дают практически одинаковые результаты. Наибольшее различие наблюдается у параметров среднего звена, что указывает на неодномерный характер распространения тепла в пластмассовом корпусе ЦИС.

Таблица 2

Тепловые параметры ЦИС, найденные разными способами обработки исходной ПТХ

Table 2

The DIC thermal parameters found with use various ways of processing of initial TTC

Алгоритм расчета Первое звено Второе звено Третье звено

Rt1 , К/Вт тт , мс Rt2 , К/Вт Тт2, мс Rt3 , К/Вт Тт3, с

Прямое дифференцирование 5,8 1,37 34,6 480 58 176

Полиномиальная регрессия 5,4 1,32 37,8 520 54,6 180

Относительная разность,% 7,4 3,7 9,3 11,5 6,2 5,0

Рис.4. Исходная ПТХ (кривая 1) и аппроксимирующие ее функции ZTl(t) (кривая 2) и ZT2(t) (кривая 3), построенные на основе одномерной тепловой

модели

Fig.4. Graphs of initial ТТС (curve 1), function Zn(t) (curve 2) and function ZT2(t) (curve 3) approximating it constructed on the basisof one-dimensional

thermal model

Первое звено тепловой схемы исследуемой ЦИС хорошо идентифицируется как звено, соответствующее кристаллу. Согласно данным производителя [12], площадь ак-

»-» тт 2

тивной области кристалла этого типа ЦИС £а « 1,2 мм , толщина кристалла d = 400 мкм. Принимая теплопроводность сильнолегированного кремния X = 80 Вт/м и учитывая, что в кольцевой генератор (область нагрева) включены из четырех только три логических эле-

3 d

мента ЦИС, для теплового сопротивления кристалла имеем КТкр =-« 5,6 К/Вт,

что практически равно среднему значению RТ , полученному по ПТХ двумя разными способами обработки. Расчетная тепловая постоянная времени кристалла хТ = d2/а

при температуропроводности кремния а = 9Т0-5 м2/с равна 1,77 мс, что примерно соответствует значениям хТ , полученным по ПТХ.

Точный расчет других звеньев затруднителен, но второе звено можно идентифицировать как соответствующее слою кристаллодержатель - основание корпуса, а третье -участку тепловой цепи основание корпуса - среда. Полное тепловое сопротивление переход - среда аналога ЦИС типа CD4011 микросхем типа 561ЛА7, согласно [13], не превышает 120 К/Вт, что соответствует полученным значениям.

Для проверки алгоритма определения тепловых параметров линейных тепловых схем КМОП ЦИС по ПТХ на рис.5 приведены графики исходной ПТХ и ее первой и второй производных, полученных прямым численным дифференцированием и усредненных по результатам измерения девяти образцов микросхем типа CD4011. Все образцы имеют три нуля второй производной ПТХ, и их тепловые схемы содержат три звена. Усредненные тепловые параметры этих схем приведены в табл.3. Наибольший разброс значений наблюдается для тепловых параметров второго звена, что связано с качеством теплоотвода от кристаллодержателя и неодномерностью теплового потока.

1(1), К/Вт

— _

> /

1/

/ 2 3

----- „.У- 1\

____— * : -

0.1 '1...........1 0 1 о2 "i О3 1 О4 1 О5 :: /: Ml

/ :

Рис.5. Исходная ПТХ (кривая 1) и ее первая (кривая 2) и вторая (кривая 3) производные ЦИС типа CD4011, усредненные по девяти образцам Fig.5. ТТС (curve 1) and their first (curve 2) and second (curve 3) derivative CD4011 DIC average on 9 samples

Таблица 3

Усредненные по девяти образцам тепловые параметры трехзвенных тепловых моделей ЦИС типа CD4011 и оценки их выборочных СКО

Table 3

The thermal parameters of the CD4011 and estimates of DIC 3-unit thermal models of their selective mean square deviation average on 9 samples

Параметр Первое звено Второе звено Третье звено

RTl, К/Вт тт , мс RTl, К/Вт Тт2, мс RTs, К/Вт Тт3, с

Усредненный тепловой 5,7 1,35 37,2 545 56,0 180

СКО 0,2 0,05 1,8 30 0,8 4

Отметим, что на графике второй производной ПТХ в диапазоне 1-10 с наблюдается локальный максимум, который у большинства образцов не достигает нуля. Этот максимум свидетельствует о наличии слоя конструкции в теле ЦИС, который может быть связан с кристаллодержателем. Особенности профиля второй производной ПТХ могут служить диагностической информацией при контроле качества.

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

ческих элементах ЦИС, имеет сравнимую с известным способом чувствительность, лишен недостатков, связанных с переходными электрическими процессами, и позволяет получать ПТХ с разрешением, достаточным для ее анализа на различных участках. ПТХ исследованных образцов КМОП ЦИС типа CD4011 имеет три характерных участка медленного роста (полки), что соответствует трехзвенной тепловой схеме, и адекватно описывается функцией в виде суммы трех экспоненциальных функций.

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

Литература

1. Закс Д.И. Параметры теплового режима полупроводниковых микросхем. - М.: Радио и связь, 1983. - 128 с.

2. Чернышов А.А., Тюхин А. А. Контроль тепловых характеристик интегральных схем // Зарубежная радиоэлектроника. - 1983. - № 5. - С. 90-95.

3. IC thermal measurement method - electrical test method (single semiconductor device) EIA/JEDEC JESD51-1 standard. - URL: www.jedec.org/download/search/jesd51-1.pdf (дата обращения: 8.02.2017).

4. T3Ster - Thermal Transient Tester. - URL: www.mentor.com/micred (дата обращения: 8.02.2017).

5. Сергеев В.А., Юдин В.В. Измерение тепловых параметров полупроводниковых изделий с применением амплитудно-импульсной модуляции греющей мощности // Метрология. - 2010. - № 4. - С. 37-47.

6. Сергеев В.А., Тетенькин Я. Г., Юдин В.В. Способ определения теплового сопротивления переход-корпус цифровых интегральных микросхем // Патент России № 2569922, 2015. Бюл. № 34.

7. Сергеев В. А., Тетенькин Я. Г. Определение тепловых параметров цифровых микросхем по температурным зависимостям времени задержки сигнала // Автоматизация процессов управления. - 2015. -№3. - С. 89-97.

8. Сергеев В. А., Тетенькин Я. Г. Алгоритм определения тепловых параметров цифровых интегральных схем по переходной тепловой характеристике // Автоматизация процессов управления. -2016. - № 1. - С. 112-119.

9. Szekely V., Tran van Bien. Fine structure of heat flow path in semiconductor devices: а measurement and identification method // Solid-State Electronics. - 1988. - Vol. 31. - No. 9. - Р. 1363-1368.

10. Вяхирев В. Измерение тепловых характеристик полупроводниковых электронных компонентов // Технологии в электронной промышленности. - 2013. - № 3. - С. 90-92.

11. Schweitzer D., Раре Н., Œеп L. Transient measurement of the junction-to-case thermal resistance using structure functions: chances and limits // Proc. of 24th IEEE SEMI-THERM Symposium. - San Jose, California, USA. - 2008. - Р. 193-199.

12. CD4011 Datasheet(PDF) - Texas Instruments. -

URL: http://pdf1.alldatasheet.com/datasheet-pdf/view/26846/TI/CD4011.html (дата обращения: 8.02.2017).

13. Горлов М. И., Николаева Е.П. Расчет надежности интегральных микросхем: методические указания к выполнению индивидуальных заданий. - Воронеж: ВГТУ, 2006. - 31 с.

Поступила после доработки 13.02.2017 г.; принята к публикации 25.04.2017 г.

Сергеев Вячеслав Андреевич - доктор технических наук, доцент, директор Ульяновского филиала Института радиотехники и электроники им. В.А. Котельникова РАН, заведующий базовой кафедрой радиотехники, опто- и наноэлектроники Ульяновского государственного технического университета при УФИРЭ им. В.А. Котельникова РАН (Россия, 432011, г. Ульяновск, ул. Гончарова, д. 48/2), sva@ulstu.ru

Тетенькин Ярослав Геннадьевич - ведущий инженер Ульяновского филиала Института радиотехники и электроники им. В.А. Котельникова РАН, соискатель базовой кафедры радиотехники, опто- и наноэлектроники Ульяновского государственного технического университета при УФИРЭ им. В.А. Котельникова РАН (Россия, 432011, г. Ульяновск, ул. Гончарова, д. 48/2), a732041@yandex.ru

References

1. Zaks D.I. Parametry teplovogo rezhima poluprovodnikovykh mikroskhem [Parameters of the thermal mode of semiconductor integrated circuits]. Moscow, Radio i svyaz', 1983. 128 p. (In Russian).

2. Chernyshov A.A., Tyukhin A. A. Kontrol' teplovykh kharakteristik integral'nykh skhem [Control of thermal characteristics of integrated circuits]. Zarubezhnaya radioehlektronika, 1983, no. 5, pp. 90-95. (In Russian).

3. IC Thermal Measurement Method - Electrical Test Method (Single Semiconductor Device) EIA/JEDEC JESD51-1 standard. Available at: www.jedec.org/download/search/jesd51-1.pdf (accessed: 10.02.2017).

4. T3Ster - Thermal Transient Tester. Available at: www.mentor.com/micred (accessed: 10.02.2017).

5. Sergeev V.A., Yudin V.V. Izmerenie teplovykh parametrov poluprovodnikovykh izdelij s primeneniem amplitudno-impul'snoj modulyatsii greyushhej moshhnosti [Measurement of thermal parameters of semiconductor products with application of pulse-amplitude modulation of the heating power]. Metrologiya - Metrologia, 2010, no. 4, pp. 37-47. (In Russian).

6. Sergeev V.A., Teten'kin YA. G., Yudin V.V. Sposob opredeleniya teplovogo soprotivleniya perekhod-korpus tsifrovykh integral'nykh mikroskhem [Method of determination of thermal resistance junction-case of digital integrated circuits]. Patent 2569922 RF, MPK G01R31/28, 2015. (In Russian).

7. Sergeev V. A., Teten'kin YA. G. Opredelenie teplovykh parametrov tsifrovykh mikroskhem po temperaturnym zavisimostyam vremeni zaderzhki signala [Determination of thermal parameters of digital integrated circuits on temperature dependences of time of a delay of a signal]. Avtomatizaciaprocessov upravlenia -Avtomatizatsiyaprotsessov upravleniya, 2015, no. 3, pp. 89-97. (In Russian).

8. Sergeev V.A., Teten'kin Ya.G. Algoritm opredeleniya teplovykh parametrov tsifrovykh integral'nykh skhem po perekhodnoj teplovoj kharakteristike [Algorithm of determination of thermal parameters of digital integrated circuits according to the transitional thermal characteristic]. Avtomatizatsiya protsessov upravleniya -Automation of control processes, 2016, no. 1, pp.112-119. (In Russian).

9. Szekely V., Tran van Bien. Fine structure of heat flow path in semiconductor devices: a measurement and identification method. Solid-State Electronics, 1988, vol. 31, no. 9, pp. 1363-1368.

10. Vyakhirev V. Izmerenie teplovykh kharakteristik poluprovodnikovykh ehlektronnykh komponentov [Measurement of thermal characteristics of semiconductor electronic components. Tekhnologii v ehlektronnoj promyshlennosti - Technologies in Electronic Industry, 2013, no. 3, pp. 90-92. (In Russian).

11. Schweitzer D., Pape H., Chen L. Transient measurement of the junction-to-case thermal resistance using structure functions: chances and limits. Proc. of 24th IEEE SEMI-THERM Symposium. San Jose, California, USA, 2008, pp. 193-199.

12. CD4011 Datasheet(PDF) - Texas Instruments. Available at: http://pdf1.alldatasheet.com/datasheet-pdf/view/26846/TI/CD4011.html (accessed: 10.02.2017).

13. Gorlov M.I., Nikolaeva E.P. Raschet nadezhnosti integral'nykh mikroskhem : metodicheskie ukazaniya k vypolneniyu individual'nykh zadanij [Calculation of reliability of integrated chips: methodical instructions to performance of individual tasks]. Voronezh, VGTU, 2006. 31 p. (In Russian).

Submitted 13.02.2017; accepted 25.04.2017.

Sergeev Viacheslav A. - doctor of engineering, associate professor; head of Ulyanovsk Branch of Kotel'nikov Institute of Radioengineering and Electronics of Russian Academy of Sciences (Kotel'nikov UBIRE RAS); head of the Radioengineering, Opto- and Nanolectronics Department, Ulyanovsk State Technical University at Kotel'nikov UBIRE RAS (Russia, 432071, Ulyanovsk, Goncharov st., 48/2,), sva@ulstu.ru.

Tetenkin Iaroslav G. - lead engineer, Ulyanovsk Branch of Kotel'nikov Institute of Radioengineering and Electronics of Russian Academy of Sciences (Russia, 432071, Ulyanovsk, Goncharov st., 48/2), ufire@mv.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.