Размещено на http://www.allbest.ru/

Содержание

• Введение
• 1. Распределение температуры по толщине пластины кремния, нагреваемой в режиме теплового потока
• 2. Описание работы программы
• 3. Исходный код программы
• 4. Выполнение индивидуального задания
• Заключение
• Список использованных источников

Введение

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

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

Учет распределения температуры по пластине очень важен, так как при изменении её на 50 градусов химические процессы могут изменить свою скорость в 10-100 раз, и часть пластины соответственно может не соответствовать режимам техпроцесса.

1. Распределение температуры по толщине пластины кремния, нагреваемой в режиме теплового потока

Диффузионное перераспределение тепла из слоя поглощения энергии начинает проявляться в полупроводниковых материалах при нагреве импульсами излучения длительностью более 10^-6 с. В нестационарных тепловых процессах перенос тепла описывается обобщённым уравнением Фурье, которое для одномерного случая имеет вид:

(1)

где с (Т) - удельная теплоёмкость;

с - объёмная плотность материала;

К (Т) - коэффициент теплопроводности.

Тепловые потери учитываются в граничных условиях:

(2)

Решение данного уравнения с записанными граничными условиями возможно только численными методами, для упрощения которых обычно используют преобразование Кирхгофа. Наиболее распространенные - конечноразностные методы:

метод Кранка-Никольсона (нелинейность исходного уравнения устраняется в процессе счета итерациями);

метод численного преобразования Лапласа-Карсона;

метод интегральных преобразований.

Аналитическое решение данного уравнения возможно при незначительных упрощениях, предполагающих адиабатические граничные условия, и при пренебрежении температурными и спектральными зависимостями параметров полупроводникового материала. Одно из них для импульса длительностью t_p имеет вид:

(3)

где x=K/ (сc), щ=р/t_p, 0В щtВр.

Аналитические решения, справедливые, как правило, в приповерхностной области полупроводниковой пластины для коротких длительностей воздействия (t_p<10^-4 с), по точности не удовлетворяют требованиям описания многих практически важных случаев, поэтому наибольшее распространение получили численные решения уравнения теплопроводности.

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

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

программа физический процесс график

Рисунок 1 - Профили относительного изменения температуры по глубине пластины кремния (380 мкм), нагреваемой импульсами излучения рубинового лазера (сплошные линии), ксеноновой лампы (штриховые) различной длительности с экспозиционной энергией 5 Дж/см²

Представленные зависимости относительной температуры (T/T_пов) по толщине кремневой пластины (380 мкм) позволяют сделать вывод о наличии значительных температурных градиентов по толщине при длительностях импульса 10^-5 - 10^-4 с. С увеличением длительности до 10^-3 - 10^-2 с они быстро уменьшаются и при экспозиции более 10^-2 с температурное поле становится практически равномерным по толщине полупроводниковой пластины. При длительностях лучистого нагрева 10^-2 с и более в полупроводниковых пластинах устанавливается квазиравномерное распределение температуры по толщине и её тепловой режим определяется балансом проводимой мощности, в котором доминирующую роль играют потери на тепловое излучение. С учётом других потерь энергии уравнение теплового баланса имеет вид:

(4)

где

(5)

Традиционно полагали, что параметры в этом уравнении не зависят от температуры и все потери, кроме радиационных, пренебрежимо малы. С этими допущениями получено аналитическое выражение для оценки максимальной температуры и времени достижения температуры T:

(6)

Более корректное описание изменения полупроводниковой пластины (в режиме теплового потока) требует учёта температурных и спектральных зависимостей параметров полупроводникового материала. При этом аналитическое решение невозможно и проводится численное интегрирование уравнения теплового баланса.

Расчёты тепловых полей при секундной термообработке излучением индуцируют незначительные градиенты температуры по толщине пластины, оцениваемые соотношением:

?T/d = уе_eT⁴/K (7)

где K - значение коэффициента теплопроводности при температуре обработки.

Проведённые расчёты для кремниевых пластин толщиной 300 мкм - 500 мкм дают даже вблизи температуры плавления разницу температур облучаемой и необлучаемой поверхности в пределах 7 - 10 градусов. Это значительно меньше возникающих при такой термообработке радиальных перепадов температуры.

Радиальные градиенты могут быть рассчитаны в приближении тонкой пластины (d<<a, a - радиус пластины) с использованием стационарного уравнения теплопроводности:

(8)

где Q = 2уе_eT⁴_e - интегральная плоскость мощности, выделенная в объёме пластины; Т_{е -} установившееся значение температуры в центре пластины (r=0); r - радиальная координата.

При этом предполагается однородность поверхностного облучения пластины и наличия лишь радиационного теплоотвода с её поверхностей.

С учётом дополнительной излучающей границы, которой является боковая поверхность пластины, граничные условия представляются в виде:

(9)

Решением уравнения теплопроводности с учетом этих граничных условий является:

(10) где (11)

I₀ и I₁ - функции Бесселя нулевого и первого порядка соответственно, которые для аргумента F??1 приближенно определяются как

(12)

Член со знаком минус в правой части решения учитывает наличие дополнительной излучающей границы.

Радиальный градиент температуры достигает максимального значения на краю пластины, а область его локализации в основном определяется её диаметром и распространяется на расстояние 0,2 - 0,3А от края пластины.

Рассмотренное описание вертикальных и горизонтальных температурных градиентов в нагреваемых излучением полупроводниковых пластинах позволяет рассчитывать связанные с ними термоупругие напряжения.

2. Описание работы программы

В данном курсовом проекте программа была разработана на языке программирования Delphi.

Для расчета графика зависимости необходимо задать толщину пластины, мощность, время в нижней части программы и нажать клавишу "РАСЧЕТ”. При этом появится график. Если изменить значение "ВРЕМЯ" и произвести новый расчет, то новый график будет отображен рядом со старым и цвет линии графика будет другим.

При нажатии на кнопку "ОЧИСТИТЬ”, производится очистка графика от предыдущих расчётов. Если график не очистить, то следующий график выводится на панель с сохранением графиков предыдущих расчётов, что удобно для наблюдения изменений расчетов.

Для выхода из программы нажимаем кнопку "ВЫХОД”.

Информация о разработчике указана вверху программы.

Общий вид главного окна программы представлен на рисунке 2.

Рисунок 2 - Общий вид главного окна программы

3. Исходный код программы

unit Unit1 interface, uses	Название модуля
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms, Dialogs, StdCtrls, ExtCtrls, TeeProcs, TeEngine, Chart, Series;	Библиотеки
type TForm1 = class (TForm) Panel1: TPanel; Chart1: TChart; Edit1: TEdit; Edit2: TEdit; Edit4: TEdit; Button1: TButton; Button2: TButton; Button3: TButton; Label1: TLabel; Label2: TLabel; Label3: TLabel; Label4: TLabel; Label5: TLabel; Label6: TLabel; ComboBox1: TComboBox; Image1: TImage; Label7: TLabel; Label8: TLabel; procedure Button1Click (Sender: TObject); procedure Button2Click (Sender: TObject); procedure Button3Click (Sender: TObject); private { Private declarations } Public { Public declarations } end;	Описание классов программы
const Ct=204; p=2330; ee=0.95; E=1.9E-11; v=0.25; av=15E6; Ts=300; g=5.68E-8;	Описание констант
var Form1: TForm1; T: extended; diam: extended; dtol: extended; W: extended; tp: extended;	Описание переменных
Implementation {$R *. dfm} function I1 (x: extended): extended; begin
I1: =exp (x) /sqrt (6.28*x); end;	Функция Бесселя
Procedure TForm1. Button1Click (Sender: TObject); var qr,qq: TLineSeries; Qrr,Qqq: extended; a,r: extended; K,Y,x: extended; Tmax: extended; begin
dtol: =strtofloat (Form1. Edit1. Text) *1E-6; diam: =strtofloat (Form1. Edit2. Text) *1E-3; W: =strtofloat (Form1. Edit4. Text) *1E4; tp: =strtofloat (Form1.comboBox1. Text) *1E-3; Tmax: =sqrt (sqrt (W*0.67/ (2*g*ee))); x: = ( (8*tp*Tmax*Tmax*Tmax*g*ee) / (Ct*p*dtol)) + +ln ( (Tmax+Ts) / (Tmax-Ts)) +2*arctan (Ts/Tmax) - tp*5.3; T: =Tmax* (exp (x) - 1) / (exp (x) +1); qr: =TLineSeries. Create (Form1. Chart1. Owner); qq: =TLineSeries. Create (Form1. Chart1. Owner);	Ввод толщины пластины Ввод диаметра пластины Ввод мощности Ввод времени Аналитическое выражение для оценки максимальной температуры и времени достижения температуры T
Form1. Chart1. AddSeries (qr); Form1. Chart1. AddSeries (qq);
qr. Title: ='Qrr T='+inttostr (round (T)); qq. Title: ='Qqq T='+inttostr (round (T));
a: =diam/2; K: =802.99/T; Y: =sqrt (k*dtol/ (8*g*ee*T*T*T)); r: =0.9*a; while r<=a do begin
Qrr: =-1E3* (av*E*g*ee*T*T*T*T*Y*Y) * ( (1/a) - (1/r) *I1 (r/Y) /I1 (a/Y)) / (K* (1-v)); Qqq: =1E3* (av*E*g*ee*T*T*T*T*Y*Y) * ( (1/a) + ( (1/r) *I1 (r/Y) /I1 (a/Y)) - ( (1/Y) *I1 (r/Y) /I1 (a/Y))) / (K* (1-v));	Уравнение теплопроводности
qr. AddXY (r/a,Qrr); qq. AddXY (r/a,Qqq); r: =r+a/300; end; end; procedure TForm1. Button2Click (Sender: TObject); begin
while Form1. Chart1. SeriesCount<>0 do Form1. Chart1. Series [0]. Destroy; end;	Очистка графика
procedure TForm1. Button3Click (Sender: TObject); begin close; end; end.	выход

4. Выполнение индивидуального задания

Выполнение тестового расчет в соответствии с данными исходного задания (кремниевая пластина, нагреваемая в режиме теплового потока).

Рисунок 3 - Тестовый расчет для d = 380 мкм, W = 500 Вт/см², tp =10^-3с

На рисунке 3 видна зависимость распределения температуры по толщине пластины кремния при заданной длительности нагрева. Градиент температуры невелик.

Рисунок 4 - Тестовый расчет для d = 380 мкм, W = 500 Вт/см², tp =10^-4с

Рисунок 5 - Тестовый расчет для d = 380 мкм, W = 500 Вт/см², tp =10^-5с

Рисунок 6 - Графики сравнения тестовых расчетов

На основе данных, полученных при тестовом расчете и графиков, изображенных на рисунке 6, можно сделать вывод о том, что при увеличении времени нагрева пластины, градиент распределения температуры уменьшается.

Заключение

В настоящее время проектирование процесса производства интегральных схем невозможно представить без автоматизации, поэтому практически все процессы рассчитываются с помощью персонального компьютера. Конкретная программа была разработана на языке Delphi.

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

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

Основываясь на данных, полученных при тестовом расчете в ходе выполнения данного курсового проекта и представленных зависимостей относительной температуры по толщине кремневой пластины, можно сделать вывод о наличии значительных температурных градиентов по толщине при длительностях импульса 10^-5 - 10^-4 с. С увеличением длительности до 10^-3 - 10^-2 с они уменьшаются. При длительностях лучистого нагрева 10^-2 с и более в полупроводниковых пластинах устанавливается квазиравномерное распределение температуры по толщине, и её тепловой режим определяется балансом проводимой мощности, в котором доминирующую роль играют потери на тепловое излучение.

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

[1] Бубенников, А.Н. Моделирование интегральных микротехнологий, приборов и схем / А.Н. Бубенников - М.; Высшая школа, 1989. - 320 с.

[2] Борисенко, В.Е. Твердофазные процессы в полупроводниках при импульсном нагреве / В.Е. Борисенко - Мн.; Наука и техника, 1991. - 260 с.

[3] Абрамов, И.И. Моделирование элементов интегральных схем: курс лекций / И.И. Абрамов - Мн.; БГУ, 1999. - 148 с.

Размещено на Allbest.ru