CC BY
42
rZa28_thesis.htm (дата обращения: 20.11.2009).
2. Волгин В.В. Склад: организация, управление, логистика. М.: Издат.-торг. корпор. «Дашков и К», 2005.
3. Прикладные нечеткие системы; под ред. Т. Тэрано, К. Асам, М. Сугено. М.: Мир, 1993.
4. Синюк В.Г., Пивненко Е.В. Применение нечеткой истинности в логическом выводе, основанном на активации сети
продукций // Интеллектуальные системы (AIS'07); Интеллектуальные САПР (CAD-2007): тр. Междунар. науч.-технич. конф. В 4-х т. М.: Физматлит, 2007. Т. 2. С. 221-223.
5. Штовба С.Д. Введение в теорию нечетких множеств и нечеткую логику // Консультационный центр Matlab компании Softline, 2006. URL: http://matlab.exponenta.ru/fuzzylogic/book1/ index.php (дата обращения: 15.12.2009).
УДК 681
ПОДСИСТЕМА УЧЕТА ВЛИЯНИЯ ОДИНОЧНЫХ СОБЫТИИ НА ХАРАКТЕРИСТИКИ МИКРОСХЕМ
Ю.К. Фортинский, к.т.н. (Воронежская государственная лесотехническая академия, [email protected])
В статье рассматриваются программное обеспечение и модели расчета реакции микросхем на воздействие отдельных ядерных частиц. Для решения задач моделирования на различных уровнях адаптированы известные средства функционально-логического моделирования и генерации тестов, схемотехнического, топологического проектирования.
Ключевые слова: автоматизация проектирования, радиация, отдельные ядерные частицы, схемотехническое и функционально-логическое моделирование, генерация тестов.
Расчет поглощенной дозы (DOS)
Расчет заряда в диэлектрике (QD)
Расчет параметров транзистора (PR)
Для проектирования микросхем, стойких к одиночным событиям, разработано проблемно-ориентированное программное обеспечение, позволяющее учесть радиационные эффекты от одиночных ядерных частиц (ОЯЧ).
Оно является одной из подсистем программного обеспечения, с помощью которого рассчитывается стойкость изделий к факторам космического происхождения. На рисунке 1 приведена структура подсистемы расчета стойкости сверхбольших интегральных схем (СБИС) к факторам космического пространства.
Основное внимание при разработке данного программного обеспечения уделялось созданию математических моделей, позволяющих учесть воздействие ОЯЧ, а для решения задач собственно моделирования на различных уровнях адаптированы известные средства функционально-логического моделирования и генерации тестов, схемотехнического, топологического проектирования.
При адаптации указанных средств проведена их модификация с целью повышения эффективности и в первую очередь снижения вычислительных затрат для обеспечения решения задач моделирования базовых элементов в реальном времени.
Алгоритм процесса моделирования сбора заряда предложено свести к решению упрощенной задачи, основанной на следующих предположениях.
Во-первых, чувствительный объем представляется в виде цилиндра, по оси которого проходит трек тяжелых заряженных частиц. Вдоль оси ци-
Графический редактор (GR)
Расчет параметров неисправных элементов (RN)
Интерфейс (IN)
Помощь (HP)
Монитор (MN)
Анализ и выбор
структуры с минимальными размерами активных областей (УМ)
Расчет заряда при возникновении одиночных событий (Q_o)
Расчет величины и формы тока ионизации при возникновении одиночных событий (1_о)
Расчет порогового заряда одиночных событий (Qp_o)
Расчет частоты одиночных событий (N_o)
СУБД
Библиотека элементов на функциональ-но-логичсском уровне
Библиотека элементов на схемотехническом уровне
Библиотека элементов
расчета деградации электропараметров
Рис. 1. Структура программного обеспечения расчета радиационных воздействий космического характера
Т
линдра направлена координатная ось z. Заряд, генерированный в треке, собирается коллекторной областью, которой является верхняя торцевая поверхность цилиндра. В начальный момент задается радиальное распределение носителей заряда в треке, а дальше рассматривается эволюция трека с учетом объемной рекомбинации носителей заряда.
Во-вторых, предполагается отсутствие внешнего поля. Таким образом, рассматривается только заряд, собранный посредством диффузии. Учет дрейфовой составляющей, как это описывалось выше, может быть осуществлен, например, путем удвоения диффузионного тока неосновных носителей.
В-третьих, объемная рекомбинация определяется двумя механизмами: Оже-рекомбиницией и рекомбинацией Шокли-Рида-Холла. Общая скорость объемной рекомбинации Я будет складываться из скоростей этих двух видов рекомбинации:
(1)
Для скорости Оже-рекомбинации Кд^ег можно записать как
дп Эр
- и р-ип 1-
(2)
К» =—=—=С п2р-п2р„ +
Аи8ег а а е 0 0
+С„ р2п-р02п0 ,
где п, р - полные концентрации носителей; п0, р0 - начальные равновесные концентрации носителей.
Значения Оже-коэффициентов Се и Сь для кремния обычно составляют порядка 2-1СГ31 см6/с.
Для скорости рекомбинации Шокли-Рида-Холла ^ин можно записать как
=-;-—-!—-- , (3)
где т ,=
%,(п+п,)+*„/(Р+Р,) 1 1
■ тр' =
- время жизни
электронов и дырок (при низком уровне инжек-ции); у,||п, у,||р - тепловая скорость электронов и дьфок соответственно; су„, стр - сечение захвата электронов и дырок соответственно.
Для низкого уровня инжекции (когда доминирует механизм рекомбинации Шокли-Рида-Холла) можно приблизительно считать, что скорость рекомбинации определяется выражением
п"по Р-Ро
к.
(4)
где тр/=тп,=т.
В-четвертых, предполагается наличие аксиальной симметрии относительно оси z (направление начального трека). Если использовать цилиндрическую систему координат, то данное предположение позволяет превратить решаемую задачу из трехмерной в двухмерную относительно независимых координат z и г (рис. 2). Плотность носи-
телей заряда не будет зависеть от полярного угла в силу наличия аксиальной симметрии относительно оси z.
Трек тяжелых заряженных частиц располагается по оси цилиндра (ось ¿); верхняя торцевая поверхность цилиндра (2=0) - коллекторная область для сбора носителей заряда из трека тяжелых заряженных частиц.
Рис. 2. Модельное представление чувствительного объема в виде цилиндра радиусом Я и высотой 20
При этих предположениях уравнения диффузии для неравновесных носителей заряда, генерированных в треке тяжелых заряженных частиц, можно записать в виде
Эр
—— 1.
Я р
(5)
а "
(6)
где Бр, Бп - коэффициенты диффузии электронов и дырок.
Общий ток, характеризующий значение собранного заряда, определяется диффузионным током неосновных носителей заряда, значение которого умножается на коэффициент порядка 2 для учета дрейфовой компоненты. Известно, что обычно в качестве чувствительной области, отвечающей за возникновение одиночных обратимых сбоев при воздействии тяжелых заряженных частиц, выступает область пространственного заряда стокового р-п-перехода закрытого п-канального транзистора. Если считать, что уровень легирования п+-стоковой области значительно выше, чем в р-подложке, то область пространственного заряда будет практически полностью сосредоточена в подложке. Следовательно, в качестве неосновных носителей заряда можно рассматривать электроны, изменение концентрации которых описывается уравнением (6).
Следует учесть, что при высоком уровне ин-жекции коэффициенты диффузии для электронов и дырок в уравнениях (5) и (6) будут равны коэффициенту амбиполярной диффузии. В качестве
границы между высоким и низким уровнями ин-жекции можно принять равенство концентраций избыточных неосновных носителей (электронов) и основных носителей заряда (дырок): п-п0=ро. Если концентрация избыточных электронов превышает равновесную концентрацию дырок, реализуется состояние с высоким уровнем инжекции, иначе - с низким.
В качестве граничных условий можно принять:
Эр
дг <3р 5/
Эр дг
= 0:
=0; дг
= 0
п| =0; дп
дь
Эп дг
0:
дп дг
= 0 .
(7)
(8)
(9)
(10) (11)
(12)
Начальное распределение плотности носителей заряда можно найти из начального распределения поглощенной энергии в треке тяжелых заряженных частиц, которое в общем случае не является однородным (в некоторых работах указывается, что распределение поглощенной энергии в треке пропорционально 1/г2, где г - расстояние от трека). Однако во многих случаях можно пренебречь неоднородным распределением поглощенной энергии в треке, поскольку начальный радиус трека намного меньше размеров области сбора заряда, для которой решается задача. Типичное значение начального радиуса трека составляет порядка 0,1 мкм, отсюда площадь поперечного сечения трека ~0,03 мкм2. Как уже указывалось, чувствительной областью для эффектов одиночных сбоев является область пространственного заряда стокового перехода закрытого п-канального транзистора. Размеры стоковой области можно определить, если известна топология базового элемента СБИС. Рассмотрим топологию с проектными нормами 350 нм. В соответствии с правилами проектирования транзисторы бывают трех типов: А-типа, Н-типа и с параллельно соединенными стоками.
Если считать, что напряжение на стоке закрытого п-канального транзистора в КМОП-паре равно напряжению питания, и взять его равным 5 В, то ширина области пространственного заряда составит порядка 0,8-2,6 мкм в зависимости от уровня легирования р-подложки (данные оценки получены для концентрации акцепторной примеси в р-подложке 1016 и 1015 см3). При этом площадь чувствительной области, определяемая площадью области пространственного заряда стокового пе-
рехода п-канального транзистора, составит 48112 мкм2. Как видно, данные значения существенно превосходят площадь поперечного сечения исходного трека тяжелых заряженных частиц. Это позволяет пренебречь неоднородностью распределения плотности носителей заряда по толщине трека.
Если считать, что длина пробега тяжелых заряженных частиц значительно превышает глубину области сбора заряда, в результате в пределах этой области можно пренебречь изменением линейной потери энергии (ЛПЭ) по длине трека (такая ситуация реализуется при длине пробега тяжелых заряженных частиц порядка 30-40 мкм и более), то концентрацию электронно -дырочных пар пе.ь в треке можно оценить следующим образом:
АЕ АЕ (13)
Б;АУ
е^г; Аг
где АЕ - энергия, поглощенная на участке трека длиной Аг; г, - энергия ионизации (для 81 8^3,6 эВ); АУ - объем участка трека длиной Аг.
Если линейные потери энергии тяжелых заряженных частиц Ь считать приблизительно постоянными, то для них можно записать
1 ДЕ 1 АЕ
р (1г р Аг
(14)
С учетом (14) выражение (13) принимает вид: рЬ
8!71Г1
Начальные концентрации электронов и дырок в треке принимаются равными концентрации электронно-дырочных пар, определяемой формулой (15).
Таким образом, алгоритм моделирования сбора заряда из трека тяжелых заряженных частиц в случае нормального падения тяжелых заряженных частиц по центру области сбора заряда проводится следующим образом.
1. Решается краевая задача для распределений концентрации носителей заряда в чувствительном объеме
Эр
т
дг2 дг2
-Се(п2р-п02р0)-
(16)
-С„(р2п-р02п0) дп
Р"Ро
Ы
-=Б
VII
удх2 + дг2 J
-Се(п2р-п02р0)-
-С„(р2п-р02п0)-п| =0:
Ъ=0
др Эр дп
г=0 ^ 2=г„ ~Эг
Эр =Эр дп Эп
Эг г=0 дг ' дг г=0 ' дг
=0:
(17)
(18)
(19)
0
г=К
р(г,г4=0)=п(г,г4=0) =
рЬ
£:71Г
, г< г(:
О, г>г4;
Б
В,
Оо,п-п0<р0
2Б Б
Ор+О,
,п-п0>р0
Оп,п-п0<р0;
2Б „Б„
Ор+О,
-,п-п0>р0
(20)
(21)
(22)
Данная задача решается численно для различных значений ЛПЭ падающих частиц, которые определяют начальную концентрацию электронно-дырочных пар в треке.
2. Вычисляется диффузионный ток неосновных носителей заряда (электронов) через верхнюю торцевую поверхность цилиндрического чувствительного объема:
г5п(г=0)
(23)
В полярной системе координат с учетом аксиальной симметрии относительно оси z выражение (23) можно преобразовать:
¡•9п(г=0;г^)
дг
■5п(/=0;г#)
в
1п И
о о ил 1 Г. КГ5П(2=0;г;0
г(1г(1ф=
г (1г=
г(1г.
(24)
3. Общий ток неосновных носителей, учитывающий дрейфовую компоненту, определяется путем удваивания диффузионного тока:
кгЭп(г=0;г;1)
I.
дг
-г (||- .
(25)
4. Собранный заряд определяется путем интегрирования выражения (25):
I к
0(1)= 4га,Ом
<3п(/=0;г;1')
Г(1Г .
(26)
Таким образом, предлагаемая методика позволяет рассчитать для заданного значения ЛПЭ падающих частиц форму импульса ионизационного тока (точнее, его диффузионной составляющей) и зависимость собранного заряда от времени после попадания частицы. При вычислении собранного заряда необходимо корректно выбрать верхний предел интегрирования в интеграле по времени. В частности, можно взять такое значение 1 в качестве верхнего предела интегрирования, при котором остаток составит не более 10 % от общего значения заряда.
После того как получена зависимость собранного заряда от значения ЛПЭ тяжелых заряжен-
ных частиц, может быть определено пороговое значение ЛПЭ тяжелых заряженных частиц для одиночных обратимых сбоев. В качестве порогового выбирается такое значение ЛПЭ, при котором собранный заряд равен заряду переключения логического состояния КМОП-триггера, для которого проводится моделирование. Величина заряда переключения обычно известна, как правило, она рассчитывается на этапе разработки конструкции и технологии СБИС.
Для предсказания частоты и вероятности возникновения одиночных сбоев при воздействии тяжелых заряженных частиц с заданными спектрально-энергетическими характеристиками, кроме порогового значения ЛПЭ Ь0, необходимо также знать значение сечения насыщения ст0. Если известны эти два параметра, можно воспользоваться, например, рекомендациями РД В 319.03.38-2000 для определения значений остальных параметров функции Вейбулла, аппроксимирующей зависимость сечения сбоев от ЛПЭ, и рассчитать частоту сбоев по формуле
ьт„
уюп= | а(Ь)ф| п (Ь)с1Ь. (27)
Ь„
где фшп(Ь) - дифференциальный ЛПЭ-спектр тяжелых заряженных частиц, см 2 с '(МэВ-см2/мг) 1: ст(Ь) - зависимость сечения событий от ЛПЭ тяжелых заряженных частиц, которая обычно аппроксимируется функцией Вейбулла.
В качестве а0 можно принять площадь торцевой поверхности цилиндрического чувствительного объема, для которого проводилось моделирование сбора заряда:
сто = яЫ2. (28)
Можно также воспользоваться известным соотношением, характерным для статических ОЗУ, между площадью запоминающей ячейки 8ЗЯ и топологической нормой проектирования /т
8 [мкм2] ~ 300(/т [мкм])2 (29)
и считать сечение насыщения равным площади запоминающей ячейки. Такой подход будет давать завышенное значение ст0, следовательно, завышенное значение частоты возникновения событий, однако относительно обеспечения стойкости изделий к эффектам одиночного события это является вполне оправданным.
Следует также отметить, что при вычислении частоты возникновения событий по формуле (27) берется сечение сбоев, приходящееся на корпус. Для этого сечение одной запоминающейся ячейки умножается на информационную емкость (общее количество бит) моделируемой СБИС. После этого проводится предварительное формирование результатов расчета. Оно заключается во включении в БД результатов расчета четырех составляющих: порогового значения энергии и массивов тока ионизации, частоты сбоев, сечения в зависимости от энергии тяжелых заряженных частиц.
8
и
0
Отметим, что оценка стойкости СБИС к воздействию тяжелых заряженных частиц в полном виде включает в себя определение оценки времени потери работоспособности СБИС и пороговой энергии наступления тиристорного эффекта, если он наблюдается. Однако очень часто проектировщику аппаратуры достаточно знать пороговое значение энергии и зависимости сечения и частоты сбоев от энергии частиц при условии, что ти-ристорный эффект отсутствует. Поэтому, как правило, расчет на этом заканчивается. Необходимо только оценить его адекватность и, если точность расчета вписывается в допустимые нормы, сформировать результаты для окончательного вывода.
Разработанная подсистема позволяет проводить расчет параметров внешней среды (радиационного воздействия), изменения параметров моделей стандартных элементов при воздействии статических видов радиации космического характера - электронного и протонного излучений, характерных для космического пространства.
Программное обеспечение радиационного моделирования включает модули:
• монитора - центрального звена средств моделирования, который обеспечивает взаимодействие всех программных модулей, а также выход на операционную систему;
• расчета накопленной дозы, определяющего характеристики внешней среды;
• расчета заряда в подзатворном диэлектрике, который накапливается при воздействии статических факторов космического пространства;
• расчета параметров моделей транзистора для схемотехнического моделирования с использованием ¿р/се-подобных программ.
Полученные значения параметров транзисторов соответствуют их деградации при воздействии данной дозы радиации. Данные об элементах заносятся в БД. Они относятся к так называемым неисправным элементам. Для моделирования предусмотрены три несправных элемента, соответствующих значению дозы уровней начального, среднего и конечного изменений параметров.
Эти уровни соответствуют следующим величинам дозы: начальной, когда изменения параметров только начинаются; конечной, при которой изменения выходят на норму технического задания или технических условий, и средней, расположенной в средней между ними точке.
Численные значения накопленной дозы вы-считываются в зависимости от вида воздействия и его ионизационной составляющей. Другими словами, данный модуль позволяет свести всю сложную совокупность внешнего воздействия к понятиям дозы.
Комплекс программ расчета стойкости к электронному и протонному излучениям обеспечивает полное представление о стойкости изделия к факторам космического пространства.
Графический редактор обеспечивает ввод информации и вывод результатов расчета в наглядной форме. Программный модуль интерфейса пользователя предназначен для взаимодействия человека и ЭВМ.
Предусмотрена обучающая подсистема.
Литература
Ачкасов В.Н. Проектирование микроэлектронных компонентов космического назначения: монография. Воронеж: ВГУ, 2005. 270 с.
УДК 004.81:159.942.52
МОДЕЛИРОВАНИЕ ЭМОЦИОНАЛЬНОГО СОСТОЯНИЯ ЧЕЛОВЕКА НА ОСНОВЕ ГИБРИДНЫХ МЕТОДОВ
В.Л. Розалиев, к.т.н.; А.В. Заболеева-Зотова, д.т.н.
(Волгоградский государственный технический университет, rozaliev_v@mail-ru, [email protected])
Объектом рассмотрения в настоящей работе являются эмоции человека. Основные задачи статьи - это определение основных понятий, характеристик, функций из теории эмоций и построение математической модели эмоции. Предложенная модель не только интересна теоретически, но и имеет практическую направленность.
Ключевые слова: эмоции человека, иммунные системы, нечеткий вывод, математические модели, акустические параметры речи.
Вопрос рассмотрения эмоций и практического применения знаний о них все чаще встает в самых разных сферах нашей жизни [1].
В современной теории эмоций положительные и отрицательные эмоции трактуются как противо-
положные, и человек в каждый момент может находиться либо в радостном, либо в подавленном эмоциональном состоянии [2]. Корреляция между положительными и отрицательными либо обратная, либо близка к нулю. Зависимость между ними
