ЗАДАНИЕ НА КУРСОВОЙ ПРОЕКТ

по проектированию и конструированию полупроводниковых приборов и элементов ИМС

1. Тема: Определение параметров модели биполярного транзистора в программе OrCAD 9.2.

2. Срок представления проекта к защите:

3. Исходные данные для научного исследования: научно-техническая информация

4. Содержание пояснительной записки курсового проекта:

   1. Модель биполярного транзистора в программе схемотехнического анализа Pspice. Статические и динамические параметры

   2. Уравнения, описывающие электрические характеристики транзистора

   3. Методы экстракции статических параметров модели биполярного транзистора из результатов измерения характеристик и параметров

   4. Методы экстракции динамических параметров модели биполярного транзистора из результатов измерения динамических характеристик и параметров

5. Перечень графического материала:

   1. Эквивалентная схема модели биполярного транзистора

Руководитель проекта, к.ф.-м.н., зав. кафедрой_________________

Задание принял к исполнению_______________________________

Содержание

Введение

1. Модель биполярного транзистора (БТ) в программе схемотехнического анализа PSpice

2. Модель биполярного транзистора в статическом режиме работы

2.1 Уравнения, описывающие статические характеристики транзистора

2.2 Параметры, описывающие статическую модель БТ

2.3 Методы экстракции статических параметров модели БТ из результатов измерения статических характеристик и параметров

3. Динамическая модель БТ в PSpice

3.1 Уравнения, описывающие электрические характеристики БТ в динамическом режиме

3.2 Параметры модели БТ в динамическом режиме

3.3 Методы экстракции динамических параметров модели БТ из результатов измерения динамических характеристик и параметров

4. Зависимость параметров модели БТ от температуры и площади

Заключение

Список использованных источников

Приложение А

Реферат

Пояснительная записка содержит 50 листов, 18 рисунков, 5 источников, 1 приложение

Перечень ключевых слов: биполярный транзистор, модель, эквивалентная схема, параметр модели, вольтамперная характеристика, ток, напряжение

Объект разработки: модель биполярного транзистора в программе схемотехнического анализа OrCAD 9.2

Цель работы: изучение и анализ характеристик и параметров модели биполярного транзистора

Методы разработки: проводился анализ научно-технической литературы

Полученные результаты: в процессе работы получено объяснение уравнений математической модели биполярного транзистора, ее параметров

Степень внедрения: не внедрено.

Эффективность: не рассчитывалась.

Область применения: применение работы в качестве методического материала, для изучения модели диода в PSpice

Основные конструктивные и технико-эксплуатационные характеристики: нет

Введение

В настоящее время машинные методы все шире используются при разработке радиоэлектронной аппаратуры. Особенно большое эти методы имеют при проектировании интегральных схем (ИС), что обусловлено их сложностью, наличием паразитных связей между компонентами.

Точность машинного расчета характеристик любой ИС практически полностью определяется точностью используемых математических моделей элементов схемы. Биполярный транзистор – один из самых распространенных приборов в радиоэлектронике, соответственно модели биполярного транзистора (БТ) есть во всех современных моделирующих программах.

При моделировании различных схем на БТ важно знать его основные параметры, способы их измерения и определения, зависимости этих параметров от различных факторов. Поэтому целью данного курсового проекта стало изучение принципов построения БТ на примере известной программы схемотехнического моделирования PSpice, которая является основой системы OrCAD 9.2.

1. Модель биполярного транзистора в программе схемотехнического анализа PSpice

Основу системы OrCAD 9.2 составляет программа PSpice, которая является наиболее известной модификацией программы схемотехнического моделирования SPICE (Simulation Program with Integrated Circuit Emphasis – моделирующая программа, ориентированная на интегральные схемы). В качестве модели биполярного транзистора (БТ) в PSpice используется один из вариантов модели Гуммеля-Пуна (Г-П), эквивалентная схема которого приведена на рисунке 1. В основе этой модели лежит передаточная эквивалентная схема, представленная на рисунке диодами I_BE1/BF, I_BС1/BR и генератором тока (I_BE1-I_BC1)/Q_B. Эти элементы отражают передачу тока основных носителей в прямом (от эмиттера к коллектору) и обратном направлении [2].

Рисунок 1 – Схема замещения биполярного n-p-n-транзистора: модель Гуммеля-Пуна [5]

Поскольку токи I_BE1 и I_BC1 отражают токи неосновных носителей, то иногда их называют идеальными составляющими, так как в идеальном транзисторе только эти токи вносят полезный вклад.

Диоды с токами I_BE2 и I_BC2 отражают паразитные токи основных носителей. Токи основных носителей влияют на режимные зависимости коэффициентов передачи токов транзистора. Эти диоды включены в эквивалентную схему рисунка 1, для того чтобы моделировать рост коэффициентов передачи токов транзистора a и b в области низких уровней инжекции [2].

Рассмотрим назначение остальных элементов эквивалентной схемы транзистора. Как показано на рисунке 1, БТ моделируется как внутренний транзистор с омическими сопротивлениями, включенными последовательно с коллектором (с), базой (b) и эмиттером (е). RE, RC, – объемные сопротивления коллектора и эмиттера. Сопротивление RB – нелинейное сопротивление базы.

Для расчета транзистора в динамическом режиме необходимо учитывать емкости p-n-переходов. Поэтому в эквивалентную схему параллельно p-n-переходам необходимо включить конденсаторы. Емкость эмиттерного перехода, представленная в схеме конденсатором C_JBE, состоит из двух составляющих: барьерной и диффузионной. Емкость коллекторного перехода, как видно из рисунка 1, разделена на две составляющие: одна из них C_JBC относится к активной части коллекторного перехода, лежащей под эмиттером, и подключена параллельно переходу коллектор – база (К-Б); другая C_BX относится к пассивной части перехода К-Б, лежащей вне эмиттера, и подключена непосредственно к базовому выводу. Диффузионная емкость коллекторного перехода целиком отнесена к емкости его активной части.

Емкость перехода коллектор – подложка (К-П) в модели содержит только одну составляющую – барьерную емкость C_JS. Вывод подложки S подключается по-разному в зависимости от конструкции БТ. Обычно цепочка, моделирующая переход К-П, подключена к коллектору. Однако для горизонтальных p-n-p транзисторов, описываемых в PSpice моделью с именем LPNP, подложка образует p-n-переход с базой транзистора и соответственно цепочка, моделирующая этот переход, подключается к базе.

Элементы I_EPI – управляемый источник тока, и две нелинейных емкости, заряды которых на рисунке 1 обозначены Q₀ и Q_W, характеризуют режим квазинасыщения. Этот режим характеризуется прямым смещением перехода внутренняя база – коллектор, в то время как переход наружная база – коллектор остается смещенным в обратном направлении. В результате искажается начальный участок выходных характеристик транзистора вблизи области насыщения [2].

Модель Г-П (рисунок 1) автоматически упрощается до более простой модели Эберса-Молла (Э-М), если опустить некоторые параметры. Эквивалентная схема этой модели для n-p-n-структуры изображена на рисунке 2. Параметры полной математической модели БТ в PSpice приведены в таблице приложения А [5].

Рисунок 2 – Схема замещения биполярного n-p-n-транзистора: передаточная модель Э-М

Электрические модели приборов делятся на статические и динамические. Статические электрические модели отражают стационарное состояние прибора при неизменных внешних управляющих воздействиях и не учитывают переходные процессы и характеристики в приборе. Динамические модели дополнительно отражают переходные процессы, электрические характеристики, изменяющиеся во времени, определяющие поведение прибора при изменении во времени управляющих воздействий.

Обычно моделированию предшествует математический анализ. Под математической моделью понимается система уравнений (или математическое описание другого вида), позволяющая определить с требуемой точностью необходимые характеристики компонента в различных условиях работы. Например, статическая модель транзистора описывается системой уравнений, связывающих токи, напряжения на выводах прибора с его параметрами при работе в статическом режиме [4].

Поэтому, для того чтобы разобраться в принципах построения модели БТ в программе PSpice рассмотрим основные уравнения, описывающие электрические характеристики транзистора. В следующих разделах будут рассмотрены статическая и динамическая PSpice модели БТ и параметры, описывающие их, зависимость параметров моделей от температуры и площади. Малосигнальная модель, а также модель, учитывающая источники шумов в рамках данной работы не рассматриваются.

2. Модель БТ в статическом режиме работы

2.1 Уравнения, описывающие статические характеристики транзистора

Отправной точкой для анализа служит "передаточный вариант" уравнений Эберса-Молла. БТ может рассматриваться как взаимодействующая пара p-n-переходов. Модель Э-М основана на суперпозиции нормального и инверсного транзисторов, работающих в активном режиме. P-n-переходы представляют в виде диодов. Часть тока диода передается через базу транзистора и собирается другим электродом. Этот ток учитывается генератором тока (ГТ), включенным в эквивалентную схему. Модель Э-М связывает токи на выводах транзистора с напряжением на p-n-переходах, поэтому она удобна для автоматического анализа электронных схем [4].

В идеальной передаточной модели в качестве токов связи используются токи, коллектируемые p-n-переходами. Эти токи в эквивалентной схеме моделируются ГТ I_CC и I_EC. В Spice программах два опорных ГТ заменяют на один источник тока I_CT, включенный между эмиттером и коллектором. I_CT определяется следующим выражением [1]:

(1)

где I_S – ток насыщения, kT/q = V_T – тепловой потенциал, V_BE – напряжение перехода Э-Б, V_BC – напряжение перехода К-Б.

Эквивалентная схема передаточной модели Э-М приведена на рисунке 3.

Рисунок 3 - Передаточная модель Э-М в Spice [1]

Тогда, согласно [1] токи, текущие через, диоды станут равными:

,(2)

(3)

где b_F и b_R, соответственно, коэффициенты передачи по току в схеме с общим эмиттером в нормальном и инверсном режимах.

Выходные токи будут равны

(4)

Эта модель рекомендуется, когда БТ работает как переключатель на постоянном токе или в определенно узком диапазоне напряжений.

Включение трех постоянных резисторов (r_E, r_B, r_C) улучшает характеристику модели по постоянному току. Они представляют омические сопротивления транзистора между его активной областью и его К, Э, и Б выводами, соответственно. Эти резисторы включены в модель, как показано на рисунке 4.

Рисунок 4 – Статическая модель Э-М в PSpice [1]

Данная модель учитывает модуляцию ширины базовой области (эффект Эрли). Как показано в [3], ширина обедненной области обратносмещенного p-n-перехода зависит от напряжения. В БТ изменение напряжения смещения коллекторного перехода вызывает изменение ширины области его объемного заряда (ООЗ) и, следовательно, ширины квазинейтральной области базы. Эти изменения представляют собой источник ряда физических эффектов, усложняющих анализ характеристик транзистора в режиме линейного усилителя.

В общем, эффект Эрли – есть зависимость тока I_S а, следовательно, и тока коллектора I_C, коэффициента передачи b_F как функции от V_BC. На рисунке 5 представлен график зависимости I_C от V_CE, иллюстрирующий эффект модуляции ширины базы.

пунктирные линии – идеальная модель;

сплошные линии – реальная модель (с учетом эффекта Эрли)

Рисунок 5 – График зависимости I_C от V_CE [3]

Анализируя работу БТ в линейной области, сначала определяется ширина базы относительно эффекта Эрли, а затем параметры, зависящие от ширины базы W_B. Результаты анализа будут записаны следующим образом [1]:

(5)

(6)

(7)

где V_А – дополнительный параметр – напряжение Эрли – определяется как

(8)

и - I_S(0), b_F(0) есть значения параметров при V_BC = 0.

Учитывая эффект модуляции ширины базы, выражение (1) следует изменить, согласно (6). Выражение для I_CT станет следующим:

(9)

А выходные токи модели будут описаны следующими уравнениями:

(10)

где GMIN - минимальная проводимость включается параллельно с каждым переходом, которая автоматически добавляется для того чтобы осуществить сходимость (ее значение по умолчанию равно 10^-12 См).

Модель Э-М не учитывает ряд физических эффектов, наблюдаемых в реальных структурах.

Два наиболее важных из них – эффекты низкого и высокого уровней инжекции.

Чтобы учесть эти эффекты второго порядка, в PSpice была реализована модель Г-П.

В модели используется обобщенное соотношение управления зарядом, которое позволяет выразить ток, передаваемый от эмиттера к коллектору через напряжения на p-n-переходах и общий заряд в базе Q_B [4].

В результате, модель Э-М можно модифицировать, описав следующие эффекты второго порядка [3]:

1 Рекомбинацию в ООЗ перехода Э-Б при малых напряжениях смещения V_BE;

2 Снижение коэффициента усиления по току, наблюдаемое при больших токах;

3 Полное описание эффекта Эрли на ток связи между Э и К.

На рисунке 6 представлен типичный график зависимости b_F от I_C.

Рисунок 7 – График зависимости b_F от I_C при V_BC = 0[1]

Как видно из рисунка 6, выделяют три области зависимости b_F от I_C: I – область слабых токов, II – область средних токов, III – область больших токов.

В области II модель Э-М имеет силу, и токи даются следующими выражениями (при V_BC = 0):

, (11)

,

где b_F в этой области принимает максимальное значение и определяется параметром b_FМ.

Снижение b_F при слабых токах вызвано дополнительными компонентами I_B, которые до сих пор не учитывались.

Для нормальной активной области при V_BC = 0, существует три дополнительных компонента, которые обусловлены следующими явлениями:

1 Рекомбинация носителей на поверхности;

2 Рекомбинация носителей в ООЗ перехода Э-Б;

3 Формирование каналов на поверхности перехода Б-Э.

Сумма трех составляющих равняется составному току, который необходимо прибавить к току базы соотношения (11) [1]:

(12)

где n_EL – коэффициент неидеальности перехода Б-Э.

Доминирующей составляющей обычно является рекомбинация в ООЗ перехода Э-Б.

Поэтому, при V_BC = 0, ток базы приближенно равен

(13)

где член I_S.composite в формуле (12) был заменен на C₂I_S(0); то есть соотношение было приведено к I_S(0).

С₂ – коэффициент, определяющий ток насыщения утечки перехода Б-Э.

Если переход Б-К смещен в прямом направлении, то выражение для I_B тогда становится равным

(14)

где n_CL - коэффициент неидеальности перехода Б-К, и C₄ – коэффициент, определяющий ток насыщения утечки перехода Б-Э.

Выражение для I_C тогда становится равным:

(15)

Дополнительные составляющие тока базы I_B включены в схему модели посредством двух неидеальных диодов, как показано на рисунке 7.

Рисунок 7 – Статическая модель Гуммеля-Пуна[1]

Для того чтобы учесть эффект сильных токов и эффект Эрли, достаточно модифицировать параметр I_S, который является множителем в выражении (1) для тока связи между эмиттером и коллектором. Согласно [1], I_S может быть определено как

(16)

где D_n – эффективный коэффициент диффузии в базе, n_i – собственная концентрация, A_J - одномерная площадь поперечного сечения, p₀(x) – равновесная концентрация дырок в нейтральной базе и областях обеднения.

Новое выражение для I_S получено из рассмотрения уравнений для плотности тока в БТ. Это новое выражение неотъемлемо включает эффекты модуляции ширины базы и высокого уровня инжекции.

В модели Г-П полный заряд основных носителей в нейтральной базе выражается символом Q_B и определяется следующим выражением:

(17)

На рисунке 8 приведен базовый профиль для одномерного n-p-n-транзистора.

Рисунок 8 – Базовый профиль для n-p-n-транзистора []

В модели Г-П заряд Q_B представлен составляющими, которые зависят от напряжения смещения и легко рассчитываются. В модели имеется встроенный заряд - заряд основных носителей в базе при нулевом смещении Q_B0, определяемый как

,(18)

где N_A(x) – концентрация акцепторов в базе.

Приведенный заряд основных носителей в базе q_B определен как

(19)

Помимо члена соотношения (18), в модель входят накопленные заряды эмиттерного C_JEV_BE и коллекторного C_JCV_BC переходов плюс заряд, связанный с прямой t_BFI_CC и инверсной t_BRI_EC инжекцией неосновных носителей базы. В результате их сложения выражение для приведенного заряда базы имеет вид:

(20)

где t_BF и t_BR - двухзарядные управляющие постоянные времени.

Уравнение (20) может быть приведено к следующей форме:

(21)

где и - напряжения Эрли для инверсного режима и эквивалентного напряжения Эрли для нормального режима, совпадают с определением, данным в выражении (8); и - смотрите выражение (1).

Заметим, что q_е и q_с моделируют эффект модуляции ширины базы, в то время как q_BF и q_BR моделируют эффекты высокого уровня инжекции.

С помощью введенных нормированных переменных уравнение (21) можно переписать в следующем виде:

(22)

где

(23)

Слагаемое q₁ моделирует эффекты модуляции ширины базы. Если эффект Эрли можно не учитывать, то q₁ стремится к единице. Слагаемое q₂ заряда базы учитывает избыточный заряд основных носителей базы, который получен из инжектированных неосновных носителей, следовательно, q₂ моделирует эффекты высокого уровня инжекции. Если эффекты высокого уровня инжекции не влияют на работу транзистора, то q₂ будет малым.

Выражение (22) выражается квадратным уравнением для q_B, конечное решение которого может быть написано как [3]

,(24)

отрицательное решение опущено, т.к. q_b положительно.

(25)

где - ток начала спада зависимости b_F от тока коллектора в нормальном режиме; - ток начала спада зависимости b_R от тока эмиттера в инверсном режиме.

Для высокого уровня инжекции, когда эффект Эрли действует значительно меньше по сравнению с эффектами высокого уровня инжекции, имеем q₂>q₁. В этой ситуации нормированный заряд базы описывается асимптотическим выражением для высоких напряжений смещения

(26)

Следовательно, формула для коллекторного тока будет иметь вид

(27)

Данная модель также принимает во внимание зависимость сопротивления базы от тока (шнурование тока). Здесь объемное сопротивление базы описывается сопротивлением базы r_BB’ между внешним и внутренним выводами базы.

Полная модель Г-П, которая учитывает описанные выше эффекты, показана на рисунке 9.

Рисунок 9 – Статическая модель БТ в PSpice

Согласно вышеприведенным соотношениям характеристики БТ по постоянному току могут быть заданы следующими уравнениями [1]:

Нормальная активная область

(28)

Инверсная область

(29)

Область насыщения

(30)

Область отсечки

(31)

2.2 Параметры, описывающие статическую модель БТ

Параметры, требуемые для модели БТ могут быть заданы в опции .МОDEL. Для описания статической модели БТ необходимо задать следующие параметры [1]:

ISТок насыщения при температуре 27°С (I_S);

BFМаксимальный коэффициент передачи тока в нормальном режиме в схеме с ОЭ (без учета токов утечки) (b_F);

BRМаксимальный коэффициент передачи тока в инверсном режиме в схеме с ОЭ (без учета токов утечки) (b_R);

NEКоэффициент неидеальности перехода Б-Э (n_EL);

NCКоэффициент неидеальности перехода К-Б (n_CL);

VAFНапряжение Эрли в нормальном режиме (V_А);

VARНапряжение Эрли в инверсном режиме (V_B);

IKFТок начала спада зависимости BF от тока коллектора в нормальном режиме (I_KF);

IKRТок начала спада зависимости BR от тока эмиттера в инверсном режиме (I_KR);

RBСопротивление базы при нулевом смещении перехода Б-Э (r_B);

RBMМинимальное сопротивление базы при больших токах (r_BM);

IRBТок, при котором сопротивление базы падает на 50% к его минимальному значению (I_rB);

ISEТок насыщения утечки перехода Б-Э (I_SE = C₂I_S);

ISCТок насыщения утечки перехода Б-К (I_SC = C₄I_S);

NFКоэффициент неидеальности в нормальном режиме (n_F);

NRКоэффициент неидеальности в инверсном режиме (n_R);

NKКоэффициент, определяющий множитель q_b (n_KF);

ISSОбратный ток p-n-перехода подложки (I_SSUB);

NSКоэффициент неидеальности перехода подложки (n_S).

Обозначения, используемые в тексте, указаны в круглых скобках.

Параметр n_F моделирует прямой ток коллектора при слабых токах - это есть показатель экспоненты тока коллектора, который определяет наклон характеристики I_C от V_BE в логарифмическом масштабе; n_R имеет тот же самый смысл в инверсном режиме.

Четырьмя параметрами модели: C2, n_EL (для b_F) и C4, n_KL (для b_R) описывается спад b при низких токах.

Модуляция ширины базы. Этот эффект описан двумя параметрами, V_А и V_B, которые отождествлены в компоненты q_C и q_E [см. выражение (21)].

Должно быть отмечено, что т.к. V_А, V_B, I_KF и I_KR не могут быть равны нулю, PSpice интерпретирует нулевые значения для этих параметров как бесконечные.

Включены три омических сопротивления RB, RC и RE, где RB может быть зависимым от сильных токов. Сопротивление r_C уменьшает наклон кривых в области насыщения для низких напряжений К-Э. Сопротивление коллектора может ограничивать токопередаточную способность БТ, и это также воздействует на максимальную рабочую частоту при больших токах.

Эмиттер - наиболее сильно легированная область в транзисторе. По этой причине, доминирующий компонент сопротивления эмиттера r_E обычно является сопротивление контакта (порядка единиц Ом). r_E, которым часто пренебрегают, необходимо обычно принимать маленьким, постоянным значением.

Объемное сопротивление базы r_BB’ между внешним и внутренним выводами базы состоит из двух отдельных сопротивлений [1]. Внешнее постоянное сопротивление r_B (внешнее сопротивление базы) состоит из сопротивления контакта и поверхностного сопротивления внешней области базы. Сопротивление внутренней области базы r_BM характеризует сопротивление активной области базы, являющейся частью базы, находящейся непосредственно под эмиттером. Это сопротивление является функцией тока базы. Зависимость этого сопротивления по току устройства возникает в результате отличного от нуля удельного сопротивления базовой области.

Можно показать, что полное сопротивление базы может быть выражено как [5]

(32)

где r_BM - минимальное сопротивление базы, которое имеет место при больших токах; r_B - сопротивление базы при нулевом смещении (маленькие токи базы), и z – переменная относительно удельного сопротивления базы, теплового потенциала и внутренней (эффективной) длины базы.

Чтобы сократить сложность вычисления в расчете z, используется метод приближения, отображающий cos z в соответствии с первыми двумя слагаемыми ряда Макларена. Значение z из этого приближения есть

(33)

где I_rB - ток, при котором сопротивление базы падает на 50 % к его минимальному значению.

Зависимость тока от сопротивления базы смоделирована в PSpice следующим образом:

(34)

где z находится по формуле (33).

В формулах (14) и (15) вместо C₂I_S, C₄I_S были введены два новых параметра: обратный ток насыщения Э перехода I_BЕ (IBЕ) и обратный ток насыщения К перехода I_BC (IBC), оба с размерностью [А]. Если и IBЕ, и IBC определены в опции .MODEL, PSpice использует их вместо IS, вычисляет соответствующие компоненты перехода в уравнении (14), (15).

Переход подложки с диодом, соединенным или с коллектором или с базой, в зависимости от того, горизонтальный или вертикальный БТ, моделируется посредством параметра SUBS. Горизонтальная геометрия подразумевается когда параметр модели SUBS = - 1, и вертикальная геометрия когда SUBS = +1. Ток подложки – от подложки к коллектору - для вертикального БТ, и от подложки к базе - для горизонтального БТ определяется следующим образом [1]:

Вертикальный БТ

(35)

Горизонтальный БТ

(36)

В приведенные выше уравнения были введены два новых параметра: I_SSUB (ISS) – обратный ток p-n-перехода подложки; и n_S (NS) - коэффициент неидеальности перехода подложки.

Окончательное уравнение заряда принимает форму []

(37)

где новый параметр n_KF (NK) – коэффициент, определяющий множитель q_b (по умолчанию равен 0,5).

2.3 Методы экстракции статических параметров модели из результатов измерения характеристик и параметров

Ток насыщения I_S в модели Э-М экстраполируется отрезком прямой тока зависимости lnI_C от V_BE в прямой области и lnI_E от V_BC в обратной области, как показано на рисунке 10. График зависимости рисунка 10 показывает смысл параметров I_S и _F.

Рисунок 10 –График зависимости lnI_C (lnI_B)от напряжения V_BE (V_BC)

Напряжение Эрли V_A может быть получено непосредственно из зависимости I_C от V_CE. Наклон этих характеристик в нормальной активной области g₀ получается из уравнений (4) и (9), путём ослабления незначительный второго члена, и тогда дифференцируя относительно V_BC (V_BE принята постоянной), получим

(38)

Геометрический смысл соотношения (38) показывает, что V_A получается из отрезка прямой экстраполирующей наклон относительно оси V_CE (как показано на кривой рисунка 5). Например, наклон (50 кОм)^-1 при I_C(0) = 1 мА дает, из уравнения (38), V_A = 50 В.

Для определения параметра _F необходимо построить график зависимости lnI_C и lnI_B как функции V_BE, как показано на рисунке 11. Так как вертикальная ось логарифмическая, b_F получается непосредственно из графика как расстояние между кривыми I_C и I_B.

Рисунок 11 – График зависимости lnI_C и lnI_B от V_BE при V_BC = 0 [1]

График зависимости lnI_B от qV_BE/kT при V_BC = 0, приведенный на рисунке 12, иллюстрирует два компонента тока I_B: идеальная компонента с наклоном 1 и неидеальная с наклоном равным 1/n_EL. Экстраполяция этих прямолинейных участков на ось у дает значения C₂I_S(0) и I_S(0)/b_FM(0).

Подобный график lnI_B как функции от V_BC для инверсного режима работы дает значения для параметров модели C₄ и n_CL. Типичное значение для C₂ (и C₄) составляет 10³, а типовое значение для n_EL (и n_CL) - 2.

Из характеристики lnI_C от V_BE в двух экстремумах – высокого и низкого уровня инжекции – можно экспериментально определить параметр I_KF.

Асимптота низкого тока дается следующим уравнением (для q_e»q_c»0):

(39)

Рисунок 12 – График зависимости lnI_C и lnI_B от qV_BE/kT при V_BC = 0

Асимптота больших токов дается следующим соотношением

(40)

Пересечение двух асимптот определяет ток излома I_KF и соответствующее ему напряжение излома V_KF.

Из соотношения (40) следует, что при высоком уровне инжекции

,(41)

в то время как из формулы (39) следует, что при низком уровне инжекции

(42)

Из решений уравнений (41) и (42) получается

(43)

Аналогично можно получить I_KR, если построить график lnI_E от V_BC.

Дополнительный параметр модели Г-П может быть получен из поведения асимптоты в результате увеличения тока короткого замыкания I_L (см. рисунок 6). Усиление по току постоянно при значении _FМ для I_C > I_L; уменьшается с наклоном 1-n_EL при I_C < I_L; и уменьшается с наклоном минус единица для токов I_C > I_KF. I_L дается приближенным соотношением:

(44)

3. Динамическая модель БТ в PSpice

3.1 Уравнения, описывающие электрические характеристики БТ в динамическом режиме

Рассмотрим эффекты накопления заряда в устройстве на примере модели Э-М. Накопление заряда в БТ моделируется включением трех типов конденсаторов: двух нелинейных конденсаторов, учитывающих барьерные емкости p-n-переходов, двух нелинейных конденсаторов, учитывающих диффузионные емкости переходов, и постоянного конденсатора перехода подложки [1].

Заряд, связанный с подвижными носителями в БТ, моделируется диффузионными емкостями. Этот заряд разделен на две составляющие: один связан с опорным источником коллекторного тока I_CC, и другой с опорным источником тока эмиттера I_EC. Каждый компонент отображается конденсатором.

Чтобы вычислить диффузионную емкость, связанную с I_CC, необходимо рассмотреть общее число подвижных зарядов, связанных с этим током. Поэтому примем, что переход Б-Э прямо смещен и V_BC=0.

Для упрощенного одномерного случая постоянно легированной базы, незначительной рекомбинации в базе, и низкого уровня инжекции в БТ (см. рисунок 13), сумма подвижных зарядов Q_DE, связанная с I_CC, может быть записана как сумма отдельных неосновных зарядов:

,(45)

где Q_E – заряд неосновных носителей, запасенный в эмиттерной области, Q_JE – заряд неосновных носителей в обедненной области перехода Э-Б, связанный с I_CC (обычно принимают равным нулю), Q_BF – заряд неосновных носителей, накопленный в нейтральной области базы, Q_JC- заряд неосновных подвижных носителей в обедненной области перехода К-Б, связанный с I_CC.

Рисунок 13 – Поперечный разрез n⁺-p-n^- -транзистора, показывающий размещение компонент заряда [1]

Чтобы определить диффузионную емкость, необходимо рассмотреть только одну составляющую. Из уравнения (45) полный заряд подвижных носителей, связанный с I_CC может также быть выражен как [1]

,(46)

где t_E – время задержки эмиттера; t_EB - время пролета через ООЗ перехода Э-Б; t_BF - время пролета базы; t_CB - время пролета ООЗ перехода Б-К; и t_F - общее прямое время пролета (принятое здесь постоянным), которое представляет среднее время для неосновных носителей, необходимое для того чтобы диффундировать через нейтральную область базы из эмиттера к коллектору. t_ЕB стремится к нулю.

Подобный анализ общего заряда подвижных носителей, связанного с I_EC приводит к

(47)

где Q_C – заряд неосновных подвижных носителей, накопленный в нейтральной области коллектора; Q_JC – заряд неосновных носителей в ООЗ перехода К-Б, связанный с I_EC; Q_BR – заряд неосновных носителей, накопленный в нейтральной области базы; Q_JE - заряд неосновных носителей в ООЗ Э-Б, связанный с I_EC. Если заряд Q_JC принять равным 0, тогда из уравнения (47) следует

,(48)

где t_C – время задержки коллектора, t_BR - обратное время пролета Б, и t_R - полное обратное время пролета (принят постоянным). t_СB стремится к нулю.

Два заряда Q_DE и Q_DC моделируются двумя нелинейными конденсаторами

(49)

как показано на рисунке 14.

Рисунок 14 – Модель Эберса-Молла для большого сигнала

С ростом приращений на переходах неподвижные заряды Q_JE и Q_JC, накопленные в обедненных областях БТ, могут быть смоделированы двумя конденсаторами – называемыми барьерными емкостями. Эти емкости, обозначенные C_JE для перехода Б-Э и C_JC для коллекторного перехода, включены в модель, как это показано на рисунке 14. Каждая емкость перехода - нелинейная функция от напряжения на выводах перехода, с которым соединена.

В [1] показано, что обе эти зависимости имеют следующий вид:

(50)

где C_J(0) – барьерная емкость при нулевом смещении, V – приложенное напряжение, f - контактная разность потенциалов перехода, m – показатель плавности перехода.

Для эмиттерного перехода и коллекторного переходов коэффициенты плавности равны по умолчанию m_E = m_C = 0,33.

Чтобы получить неподвижные заряды Q_JE и Q_JC, необходимо проинтегрировать барьерные емкости по их напряжению, то есть

(51)

На рисунке 15 показаны три кривые зависимости барьерной емкости как функции напряжения [1].

кривая (а) соответствует выражению (50)

кривая (b) показывает конечное изменение барьерной емкости

кривая (с) описывается соотношением (52)

Рисунок 15 – График изменения барьерной емкости с напряжением

Кривая (с) рисунка 15 представляет прямолинейное приближение, сделанное в соответствии с обычными компьютерными программами для V> f/2. Уравнение для этой прямой линии, полученной для соответствующего наклона в f/2, определяется как

для V ³ f/2(52)

В этом приближении удается избежать бесконечной емкости. Но оно не столь точно как кривая Чавла-Гуммеля (b), однако, приемлемо потому, что под прямым смещением диффузионные емкости, доминируют и неотъемлемо включают эффект заряда подвижных носителей в обедненных областях [1].

В Spice используется прямолинейная аппроксимация для C_J подобная линии (с) рисунка 15. Уравнение (52) заменено следующим общим соотношением:

для V ³ 0(53)

Помимо C_JE и C_JC, при проектировании интегральных схем должна быть принята во внимание еще одна емкость: емкость подложки C_JS.

Хотя фактически это барьерная емкость в области с изменяющимся потенциалом эпитаксиальный слой – подложка, здесь она смоделирована как конденсатор с постоянным номиналом.

Это представление адекватно для большинства случаев, так как переход ЭС – подложка смещен в обратном направлении в целях изоляции.

Установив основные соотношения эффектов накопления заряда, покажем, как реализована модель Э-М для большого сигнала в Spice. Компоненты накопленных зарядов Q_BE = Q_DE + Q_JE и Q_BC = Q_DC + Q_JC моделируются конденсаторами C_BE и C_BC, включенными в эквивалентную схему модели так, как показано на рисунке 16.

Рисунок 16 – Модель большого сигнала Эберса-Молла в Spice2

Компоненты заряда накопления представлены в PSpice следующими зависимыми от напряжения уравнениями емкости [1]:

,(54)

,(55)

,(56)

где для эмиттерного перехода

,(57)

для коллекторного перехода

,(58)

где FC – коэффициент нелинейности барьерных емкостей прямосмещенных переходов, принимающий значения от 0 до 1. Коэффициенты плавности переходов хотя и включены в соотношения (54) – (58), фактически не учитываются в модели Э-М

В PSpice схема модели большого сигнала Г-П идентична схеме, приведенной на рисунке 16. Зависимые от напряжения емкости, определены соотношениями (54) - (58), причем здесь учитываются коэффициенты плавности m_E, m_C и m_S (обычно, они изменяются между 0,33 и 0,5), а I_EC и I_CC рассматриваются как функции I_SS и q_B. Кроме того, модель большого сигнала Г-П, учитывает три дополнительных эффекта: распределенная емкость перехода Б-К, модуляция времени переноса заряда t_F, и распределенные явления в области базы (стадия избытка).

3.2 Параметры модели БТ в динамическом режиме

Для описания модели реального БТ, работающего на большом сигнале, необходимо задать следующие параметры модели [5]:

CJEЕмкость эмиттерного перехода при нулевом смещении (C_JE);

CJCЕмкость коллекторного перехода при нулевом смещении (C_JC);

CJSЕмкость перехода коллектор - подложка при нулевом смещении (C_JS);

VJEКонтактная разность потенциалов перехода Б-Э (f_E);

VJCКонтактная разность потенциалов перехода Б-К (f_C);

VJSКонтактная разность потенциалов перехода коллектор - подложка (f_S);

TFВремя переноса заряда через базу в нормальном режиме (t_F);

TRВремя переноса заряда через базу в инверсном режиме (t_R);

FCКоэффициент нелинейности барьерных емкостей прямосмещенных переходов (FC);

MJEКоэффициент плавности эмиттерного перехода (m_E);

MJCКоэффициент плавности коллекторного перехода (m_C);

MJSКоэффициент плавности перехода коллектор – подложка (m_S);

ITFТок, характеризующий зависимость TF от тока коллектора при больших токах (I_tF);

PTFДополнительный фазовый сдвиг на граничной частоте БТ f_Т=1/(2pTF) (P_tF);

VTFНапряжение, характеризующее зависимость TF от смещения база-коллектор (V_tF);

XCJCКоэффициент расщепления емкости база-коллектор CJC (X_CJC);

XCJC2Коэффициент расщепления емкости база-коллектор CJC (X_CJC2);

XTFКоэффициент, определяющий зависимость TF от смещения база-коллектор (X_tF)

Обозначения, используемые в тексте, обозначаются в круглых скобках.

Коэффициенты плавности перехода в модели Эберса-Молла устанавливаются по умолчанию равными 0,33.

FC – принимает значения от 0 до 1 и используется для вычисления напряжения (FC ´ f_E и FC ´ f_C) в области прямого смещения, вне которой, емкость смоделирована линейной экстраполяцией. Это сделано, для того чтобы предотвратить бесконечные емкости при V = f_Е и при V = f_С, и следовательно, гарантировать непрерывную функцию для емкостей и производных. По умолчанию FC установлена в PSpice равной 0,5.

Емкость перехода К-Б разделена на две составляющие: одна из них относится к активной части коллекторного перехода (емкость X_CJCC_JC включена между внутренним выводом базы и коллектором), другая к пассивной части коллектора (C_JC(1 - X_CJC) - емкость от внешней базы до коллектора). В качестве параметра в модели задается величина X_CJC, которая равна отношению барьерной емкости активной части перехода К-Б к полной барьерной емкости. X_CJC изменяется между 0 и 1, а по умолчанию задается равным 1, то есть емкость пассивной части вообще не учитывается. При задании параметра XCJC необходимо учитывать, что он определяется не только геометрическими размерами активной и пассивной частей перехода, но и различием удельных барьерных емкостей донной и боковой частей перехода [2].Этот параметр обычно важен только на СВЧ.

Компоненты накопленного заряда представлены в PSpice следующими уравнениями емкости зависимыми от напряжения:

,(59)

где F₁, F₂, и F₃ определяются выражением (58).

На рисунке 17 изображена полная модель БТ на большом сигнале с дополнением эффекта распределенной емкости перехода Б-К.

Рисунок 17 – Модель большого сигнала Г-П в PSpice [1]

Диффузионная емкость пропорциональна среднему времени пролета t_F носителей заряда через базу в прямом направлении и дифференциальной проводимости .

Параметр t_F учитывает зависимость времени пролета от уровня инжекции и от напряжений на переходах, то есть учитываются эффекты модуляции ширины базы. В модели принята следующая аппроксимация режимных зависимостей t_F [2]:

,(60)

где x = .

В этом выражении сомножитель 3x²-2x³ при ITF>0 отражает рост t_F при повышении уровня инжекции, что характерно для эффектов Кирка. Впрочем, по умолчанию ITF = 0 и, следовательно, этот эффект не описывается. Экспоненциальный сомножитель описывает спад t_F с ростом запирающего напряжения на коллекторном переходе, что связано с эффектом Эрли. Но по умолчанию VTF = ¥, и этот эффект не учитывается. Кроме того, XTF = 0, и если не задать XTF > 0, то никакие режимные зависимости t_F не учитываются. В этом случае t_F = TF, где TF – параметр модели [2].

Таким образом, прямой и обратный переносимый заряд можно смоделировать заданием параметров t_F, X_tF, V_tF , I_tF и t_R.

3.3 Методы экстракции динамических параметров модели БТ из результатов измерения динамических характеристик и параметров

Характер изменения t_F от I_C показан на рисунке 18. Изменение t_F при больших токах К обычно определяется эмпирическим уравнением, полученным из произведения полосы усиления f_T и тока I_C, при различных напряжениях К-Э V_CE.

Рисунок 18 – График зависимость t_F от lnI_C

В области средних токов, f_T находиться по его пиковому значению и почти постоянно; при этом время переноса заряда - время, необходимое н.з. для того, чтобы пересечь область базы и ООЗ коллекторного перехода. Диффузионная емкость перехода Э-Б увеличивается с током, уменьшая рост дифференциальной проводимости, имеющий результатом определенный предел для f_T. Таким образом, идеальный максимум t_F определяется из выражения [1]:

(61)

На больших токах, f_T и к тому же t_F становятся функцией I_C и V_CE и перестают быть постоянными (см. рисунок 18). Физические эффекты типа эффекта Кирка увеличивают время переноса и уменьшают f_T. Эти эффекты смоделированы следующей эмпирической функцией [1]:

(62)

t_F умножают на ATF в уравнениях заряда. Постоянная 1,44 просто дает интерпретацию V_tF как значение V_BC, где экспонента равняется 1/2. X_tF управляет полным спадом из-за f_T ; V_tF преобладает над изменением по f_T относительно V_CE; I_tF доминирует над изменением по f_T относительно тока.

Как видно из рисунка 18, I_tF может быть получен путем экстраполяции прямой до пересечения с осью lnI_C.

Можно показать [1], что

(63)

В области слабых токов или высоких V_CE (АТF = 1), выражение (63) сокращается до

(64)

При больших I_C, таких, что I_CC/(I_CC + I_tF) » 1 и средних V_CE , формула (63) сократиться до

(65)

Таким образом, асимптота сильных токов для данных V_BC определена параметрами X_tF и V_tF. Аналогично, асимптотическая зависимость выражения (63) в экстремумах от V_BC будет

, для I_C®¥, V_BC®0(66)

когда ATF » 1+X_tF , при I_C ® ¥, V_BC ® 0(67)

Таким образом, максимально возможный спад в f_T управляется параметром X_tF. В PSpice, этот эффект выражен следующим выражением для заряда и эквивалентной емкости:

, (68)

,

где t_FF - модулированное время переноса заряда через базу, данное выражением:

,(69)

и t_F - идеальное время переноса заряда в активном режиме. t_F экстраполируется на ось y зависимости t_FF от lnI_C.

4. Зависимость параметров модели БТ от температуры и площади

Температурные зависимости параметров элементов эквивалентной схемы БТ устанавливается с помощью следующих выражений [1].

Здесь могут устанавливаться несколько температурных уравнений для РSpice параметров модели БТ, которые можно выбрать, установив параметры TLEV и TLEVC в опции .MODEL. В последующем, мы будем рассматривать только уравнения, выбранные с TLEV.

Температурная зависимость ширины ЗЗ Еg (ЕG) следует из выражения

Температурная зависимость b_F (BF) определяются уравнением

Температурная зависимость I_S (IS) моделируется формулой

I_BЕ (IBE) и I_BC (IBC) определены

Температурная зависимость I_SSUB (ISS) определена как

Зависимости параметров I_KF (IKF), I_KR (IKR) и I_rB (IRB) от температуры представлены следующим образом:

где ТIKF1, ТIKR1, TIRB1 и TIKF2, TIKR2, TIRB2 температурные коэффициенты первого и второго порядка для соответствующих параметров, соответственно.

Следующие параметры определены для случая, когда соответствующие температурные коэффициенты определяются независимо от значения TLEV

Наконец, сопротивления, как функция температуры независимо от значения TLEV, определены следующим образом []:

В вышеупомянутых уравнениях, коэффициенты, заканчивающиеся на 1 - температурные коэффициенты первого порядка, заканчивающиеся на 2 – температурные коэффициенты второго порядка для соответствующего параметра.

Скалярный коэффициент AREA позволяет учесть параллельное соединение однотипных транзисторов, для чего в приведенной выше модели БТ изменяются следующие параметры [5]. На параметр AREA необходимо умножить все токи, емкости и заряды, а все сопротивления поделить на AREA. AREAВ и AREAC масштабный размер области базы и области коллектора. AREAВ или AREAC используются для вычисления, и выбирается в зависимости от вертикальной или горизонтальной геометрии (задание параметра модели SUBS). Для вертикальной геометрии AREAВ – масштабный коэффициент (коэффициент пересчета) для IBC, ISC и CJC. Для горизонтальной геометрии масштабный коэффициент AREAC.

Значение AREA указывается в задании на моделирование при включении транзистора в схему, по умолчанию AREA = 1.

Заключение

В результате проведенной работы изучена PSpice модель БТ и параметры для ее описания. В данном проекте были получены основные соотношения для расчета некоторых параметров модели транзистора, зависимости этих параметров от температуры и конструкции, рассмотрены методы экстракции параметров модели из экспериментальных характеристик.

Анализ PSpice модели БТ показал, что наряду с достоинствами этой модели есть и существенные недостатки. В целом модель биполярного транзистора в PSpice может с высокой точностью и в широком диапазоне напряжений, токов и частот описывать характеристики реальных приборов. Но для этого параметры модели должны быть тщательно идентифицированы по достоверным экспериментальным данным. Для идентификации может использоваться входящая в OrCAD 9.2 программа Model Editor. А получение достоверных исходных данных, особенно на высоких частотах, требует применения очень точной измерительной аппаратуры. Поэтому рядовой пользователь обычно не может квалифицированно идентифицировать параметры модели. А использование значений параметров по умолчанию, как было показано выше, не может обеспечить приемлемой точности расчетов.

Автоматическое проектирование ИС распространяется все шире и становится практически единственным инструментом в этой области. Поэтому знание основ модели необходимо для проектировщиков любого уровня.

Список использованных источников

1 Massobrio G., Antognetti P. Semiconductor Device Modeling with SPICE. Second Edition. McGraw-Hill, Inc. 1988. – 479 p.

2 Архангельский А.Я. PSpice и Design Center. В 2-х ч. Часть 1. Схемотехническое моделирование. Модели элементов. Макромоделирование. Учебное пособие. - М.: МИФИ, 1996. - 236 с.

3 Маллер Р., Кейминс Т. Элементы интегральных схем: Пер. с англ. – М.: Мир, 1989. – 630 с., ил.

4 Носов Ю.Р. и др. Математические модели элементов интегральной электроники. - М.: Сов. Радио. 1976. – 304 с.

5 Разевиг В.Д. Система сквозного проектирования электронных устройств Design Lab 8.0. – М.: Солон 1999. – 698 с.

Приложение А

Таблица А - Параметры модели биполярного транзистора

Обозначение параметра	Параметр	Разме-рность	Значение по умолчанию
AF	Показатель степени, определяющий зависимость спектральной плотности фликкер-шума от тока через переход		1
BF	Максимальный коэффициент передачи тока в нормальном режиме в схеме с ОЭ (без учета токов утечки)		100
BR	Максимальный коэффициент передачи тока в инверсном режиме в схеме с ОЭ		1
CJC	Емкость коллекторного перехода при нулевом смещении	Ф	0
CJE	Емкость эмиттерного перехода при нулевом смещении	пФ	0
CJS (CCS)	Емкость коллектор-подложка при нулевом смещении	Ф	0
EG	Ширина запрещенной зоны	эВ	1,11
FC	Коэффициент нелинейности барьерных емкостей прямосмещенных переходов		0,5
GAMMA	Коэффициент легирования эпитаксиальной области		10-11
IKF (IK)*	Ток начала спада зависимости BF от тока коллектора в нормальном режиме	А	∞
IKR*	Ток начала спада зависимости BR от тока эмиттера в инверсном режиме	А	∞
IRB*	Ток базы, при котором сопротивление базы уменьшается на 50% полного перепада между RB и RBM	А	∞
IS	Ток насыщения при температуре 27°С	А	10-16
ISC (C4)*	Ток насыщения утечки перехода база-коллектор	А	0
ISE (C2)*	Ток насыщения утечки перехода база-эмиттер	А	0
ISS	Обратный ток p-n-перехода подложки	А	0
ITF	Ток, характеризующий зависимость TF от тока коллектора при больших токах	А	0
KF	Коэффициент, определяющий спектральную плотность фликкер-шума		0
MJC (МС)	Коэффициент, учитывающий плавность коллекторного перехода		0,33
MJE (ME)	Коэффициент, учитывающий плавность эмиттерного перехода		0,33
MJS (MS)	Коэффициент, учитывающий плавность перехода коллектор-подложка		0
NC*	Коэффициент неидеальности коллекторного перехода		1,5
NE*	Коэффициент неидеальности перехода база-эмиттер		1,5
NF	Коэффициент не идеальности в нормальном режиме		1
NK	Коэффициент, определяющий множитель Qb		0,5
NR	Коэффициент неидеальности в инверсном режиме		1
NS	Коэффициент неидеальности перехода подложки		1
PTF	Дополнительный фазовый сдвиг на граничной частоте транзистора fГР=1/(2Πtf)	градус	0
QCO	Множитель, определяющий заряд в эпитаксиальной области	Кл	0
RB	Объемное сопротивление базы (максимальное) при нулевом смещении перехода база-эмиттер	Ом	0
RBM*	Минимальное сопротивление базы при больших токах	Ом	RB
RC	Объемное сопротивление коллектора	Ом	0
RCO	Сопротивление эпитаксиальной области	Ом	0
RE	Объемное сопротивление эмиттера	Ом	0
TF	Время переноса заряда через базу в нормальном режиме	с	0
TR	Время переноса заряда через базу в инверсном режиме	с	0
TRB1	Линейный температурный коэффициент RB	0C-1	0
TRB2	Квадратичный температурный коэффициент RB	0C-2	0
TRC1	Линейный температурный коэффициент RC	0C-1	0
TRC2	Квадратичный температурный коэффициент RC	0C-2	0
TRE1	Линейный температурный коэффициент RE	0C-1	0
TRE2	Квадратичный температурный коэффициент RE	0C-2	0
TRM1	Линейный температурный коэффициент RBM	0C-1	0
TRM2	Квадратичный температурный коэффициент RBM	0C-2	0
T_ABS	Абсолютная температура	0C
T_MEASURED	Температура измерений	0C
T_REL_GLOBAL	Относительная температура	0C
T_REL_LOCAL	Разность между температурой транзистора и модели-прототипа	0C
VAF (VA)*	Напряжение Эрли в нормальном режиме	В	∞
VAR (VB)*	Напряжение Эрли в инверсном режиме	В	∞
VJC (PC)	Контактная разность потенциалов перехода база-коллектор	В	0,75
VJE (PE)	Контактная разность потенциалов перехода база-эмиттер	В	0,75
VJS (PS)	Контактная разность потенциалов перехода коллектор-подложка	В	0,75
VO	Напряжение, определяющее перегиб зависимости тока эпитаксиальной области	В	10
VTF	Напряжение, характеризующее зависимость TF от смещения база-коллектор	В	∞
XCJC	Коэффициент расщепления емкости база-коллектор CJC		1
XCJC2	Коэффициент расщепления емкости база-коллектор CJC		1
ХТВ	Температурный коэффициент BF и BR		0
XTF	Коэффициент, определяющий зависимость TF от смещения база-коллектор		0
ХТI (РТ)	Температурный коэффициент IS		3
* Только для модели Гуммеля-Пуна