ПРИБОРОСТРОЕНИЕ INSTRUMENTATION ENGINEERING
https://doi.org/10.21122/2227-1031-2023-22-6-477-486 УДК 621.923
Обеспечение равномерности съема припуска при магнитно-абразивной финишной обработке полупроводниковых пластин за счет управления режимами процесса обработки
Докт. техн. наук Д. А. Степаненко1), асп. Е. С. Еромин1)
^Белорусский национальный технический университет (Минск, Республика Беларусь)
© Белорусский национальный технический университет, 2023 Belarusian National Technical University, 2023
Реферат. В статье описана методика решения прямых и обратных задач моделирования процесса магнитно-абразивной финишной обработки (МАФО) полупроводниковых пластин. Решение прямой задачи позволяет производить расчет функции съема припуска при заданных параметрах обработки, а решение обратной задачи - расчет значений параметров обработки, требуемых для реализации заданной функции съема припуска. Прямая задача решается на основе уравнения Престона, обычно используемого для описания скорости съема припуска при полировании оптических деталей. Обратная задача рассматривается в матричной формулировке, а ее решение в смысле наименьших квадратов определяется с помощью обобщенной обратной матрицы Мура-Пенроуза. На основе решения прямой задачи с постоянными значениями кинематических и магнитных параметров обработки показано, что МАФО при постоянных значениях параметров не обеспечивает равномерности съема припуска. На основе численных примеров показано, что близкий к равномерному съем припуска может быть обеспечен за счет управления магнитными параметрами обработки, закон изменения которых определяется путем решения обратной задачи. Показано, что гладкость решения исходной плохо обусловленной обратной задачи может быть повышена путем ее регуляризации, по Тихонову, что, в свою очередь, технически упрощает управление параметрами обработки.
Ключевые слова: магнитно-абразивная финишная обработка; полупроводниковые пластины; равномерность съема припуска; уравнение Престона; плохо обусловленные обратные задачи; регуляризация, по Тихонову
Для цитирования: Степаненко, Д. А. Обеспечение равномерности съема припуска при магнитно-абразивной финишной обработке полупроводниковых пластин за счет управления режимами процесса обработки Д. А. Степа-ненко, Е. С. Еромин // Наука и техника. 2023. Т. 22, № 6. С. 477-486. https://doi.org/10.21122/2227-1031-2023-22-6-477-486
Ensuring Uniform Material Removal with Magnetic Abrasive Finishing of Semiconductor Wafers by Mode Control of Machining Parameters
D. A. Stepanenko, E. S. Eromin
''Belarusian National Technical University (Minsk, Republic of Belarus)
Abstract. The paper describes a solution technique for direct and inverse problems of modeling the process of magnetic abrasive finishing (MAF) of semiconductor wafers. Solution of the direct problem enables calculation of allowance removal
Адрес для переписки
Степаненко Дмитрий Александрович Белорусский национальный технический университет ул. Я. Коласа, 22,
220013, г. Минск, Республика Беларусь Тел.: +375 17 293-91-01 stepd@tut.by
■ Наука
итехника. Т. 22, № 6 (2023)
Address for correspondence
Stepanenko Dmitry A.
Belarusian National Technical University
22, Ya. Kolasa str.,
220013, Minsk, Republic of Belarus
Tel.: +375 17 293-91-01
stepd@tut.by
function for the prescribed machining parameters, and solution of the inverse problem enables determination of machining parameters required for realization of the prescribed allowance removal function. The direct problem is solved by means of Preston equation, which is usually used for description of material removal rate when polishing optical parts. The inverse problem is considered in a matrix formulation, and its least squares solution is determined by means of generalized inverse Moore-Penrose matrix. Based on the solution of the direct problem with constant values of kinematic and magnetic machining parameters shows that MAF with constant values of machining parameters does not ensure uniformity of material removal. On the basis of numerical examples it is shown that close to uniform material removal can be ensured by control of magnetic machining parameters with the law of variation determined from solution of the inverse problem. It is demonstrated that the smoothness of solution of the initial ill-conditioned inverse problem can be improved by means of Tikhonov's regularization, which in turn technically simplifies control of machining parameters.
Keywords: magnetic abrasive finishing, semiconductor wafers, uniformity of material removal, Preston equation, ill-conditioned inverse problems, Tikhonov's regularization
For citation: Stepanenko D. A., Eromin E. S. (2023) Ensuring Uniform Material Removal with Magnetic Abrasive Finishing of Semiconductor Wafers by Mode Control of Machining Parameters. Science and Technique. 22 (6), 477-486. https://doi.org/10.21122/2227-1031-2023-22-6-477-486 (in Russian)
Введение
Магнитно-абразивная финишная обработка (МАФО) представляет собой способ финишной обработки поверхностей, основанный на создании силы резания путем воздействия магнитного поля на частицы ферро-абразивного (магнитно-абразивного) порошка (ФАП), в качестве которого могут применяться композиционные материалы на основе технического железа (магнитная фаза) и электрокорунда, алмаза или карбидов вольфрама, молибдена, ниобия и титана (абразивная фаза) [1]. Под действием магнитного поля ФАП может приобретать условно связанное состояние с формированием так называемой абразивной щетки, жесткость которой может регулироваться путем изменения характеристик магнитного поля. Таким образом, путем регулировки режимов МАФО можно в зависимости от требований к характеристикам обработанной поверхности реализовать процесс обработки, близкий к обработке свободным или связанным абразивом. При МАФО количество режущих элементов на единицу площади обрабатываемой поверхности значительно выше, чем при традиционном шлифовании, что позволяет в отличие от шлифования обеспечить высокую производительность обработки при малых давлениях на обрабатываемую поверхность, что, в свою очередь, приводит к существенному уменьшению глубины дефектного слоя. Одним из перспективных направлений применения МАФО является обработка пластин из кремния и других полупроводниковых материалов [1, с. 188; 2].
Традиционно финишная обработка кремниевых пластин производится методом химико-
механического полирования (ХМП), основанным на использовании абразивной суспензии, содержащей сильный окислитель, например пероксид водорода [3]. Очевидным недостатком такой технологии является необходимость использования вредных для здоровья человека и окружающей среды химических реактивов.
Важная характеристика поверхностей кремниевых пластин - их шероховатость: среднеарифметическое значение шероховатости Sa по поверхности (двумерное обобщение известного параметра Ra), достижимое с помощью ХМП, составляет 0,4 нм [4]. Как показывают исследования, шероховатость поверхности раздела полупроводниковой подложки и диэлектрического затвора оказывает существенное влияние на подвижность носителей заряда в электронных устройствах типа МДП-транзисторов и соответственно на их работоспособность [5].
Вторая важная характеристика кремниевых пластин, которую требуется обеспечить при их обработке, - общее изменение (отклонение) толщины по пластине TTV (Total Thickness Variation) [6]. Требования к этому параметру становятся все более жесткими с уменьшением размеров элементов микросхем. Предельно достижимый размер элементов микросхемы при фотолитографии обратно пропорционален числовой апертуре используемого объектива, а глубина резкости обратно пропорциональна квадрату числовой апертуры [6]. Как следствие, миниатюризация элементов микросхемы сопровождается значительным уменьшением глубины резкости, что, в свою очередь, приводит к повышенным требованиям к плоскостности пластин, так как отклонения от плоскостно-
Наука
итехника. Т. 22, № 6 (2023)
сти не должны превышать глубины резкости. Для сверхплоских кремниевых пластин параметр ТТУ не должен превышать 1 мкм. МАФО способна обеспечивать требуемые для кремниевых пластин значения шероховатости (0,7-2,0 нм по параметру Яа) и минимальную глубину дефектного слоя [1], однако обеспечение равномерного съема припуска (и, как следствие, минимальных отклонений от плоскостности) требует управления скоростью съема припуска и/или продолжительностью контакта инструмента с различными участками поверхности заготовки в зависимости от положения инструмента относительно заготовки. Данная статья посвящена исследованию возможности повышения равномерности съема припуска при МАФО пластин из полупроводниковых материалов за счет управления режимами процесса обработки.
Прямая задача моделирования
процесса МАФО
При моделировании съема припуска будем считать, что скорость съема припуска может быть рассчитана по формуле Престона [7, 8]:
h(p, Ф) = ф, t)dt,
(2)
RR(p, ф, t) = = к • Р(Р, Ф, Рс (t), Фc (t))v(p, Ф, pc (t), фс (t)),
(1)
где к - коэффициент, зависящий от условий обработки; р(р, ф, рс(0, фс(0) - давление инструмента на обрабатываемую поверхность в точке с координатами (р, ф) в момент времени I, соответствующий расположению оси вращения инструмента в точке (рс, фс); у(р, ф, рс(0, фс(0) -скорость движения точки (р, ф) относительно инструмента.
Для упрощения модели считаем коэффициент к постоянным и не будем учитывать его в дальнейшем анализе. При моделировании рассматривается кинематика, представляющая собой сочетание вращательного движения заготовки с угловой скоростью вращательного движения инструмента (полюсного наконечника магнитной системы) с угловой скоростью ю и радиального поступательного движения инструмента со скоростью \г.
Функция съема припуска может быть определена путем интегрирования уравнения (1)
■ Наука
итехника. Т. 22, № 6 (2023)
где Т - общая продолжительность процесса обработки.
Прямая задача моделирования процесса МАФО состоит в определении функции (2) при известной скорости (1) съема припуска и легко решается путем численного интегрирования. Для решения прямой задачи в двумерной постановке введем на поверхности обрабатываемой пластины регулярную сетку точек, лежащих на концентрических окружностях с радиусами р{ = /Ар, где / = 0...^г, Ар = Я/Мг, Я - радиус пластины. В случае равномерного окружного распределения точек с угловым шагом Аф = 2л/Ыа они будут иметь угловые координаты ф; = /Аф, где / = 0...Ыа -1. Если рс - текущее значение радиальной координаты оси вращения инструмента, а фс - текущее значение ее угловой координаты, то точка р (рг-, ф/) поверхности пластины принадлежит области проекции полюсного наконечника при выполнении следующих условий:
pc -r < p, < pc+r;
(3)
Фс - arcsin(r/pc ) < ф; < Фс + arcsin(r/pc ); (4)
p2 + p2 - 2pcp, cos^c - Фj) < r2,
(5)
где г - радиус полюсного наконечника; 2а1гат(г/рс) - угол охвата зоны обработки.
Условие (5) является достаточным, однако его непосредственная проверка для всего массива точек требует значительных затрат времени на вычисления. В связи с этим вводятся два избыточных условия (3)-(4), позволяющих сузить множество индексов ((, /), для которых требуется проверка условия (5). В случае рс(0 > г проверка (5) проводится для значений индексов, лежащих в следующих диапазонах:
ceil((pc - r)/Ap) <, <
< min(floor((pc + r )/Ap), Nr );
(6)
/1 < / < /2 при условии /2 >/1;
>
(0 < j < j2) u ( ji < j < Na-1) при условии j2 < j1 ;
j = if(ceil(mod(9c - arcsin(r/pc), 2я)/Аф) = = Na ,0, ceil(mod(фc - arcsin(r/pc), 2я)/Аф));
j2 = floor(mod(фc + arcsin(r/pc), 2я)/Аф),
где ceil - оператор округления до ближайшего большего целого значения; floor - оператор округления до ближайшего меньшего целого значения; mod(a, b) - остаток от деления числа a на число b; функция if(cond, vall, val2) принимает значение vall при cond = 1 (истинное логическое условие) и val2 при cond = 0 (ложное логическое условие).
В формуле (6) при определении верхней границы изменения индекса i значение Nr принимается в случае, когда pc(t) > R - r, то есть полюсный наконечник имеет неполное перекрытие с поверхностью заготовки. Случай j2 < j1 возникает при пересечении зоной обработки радиуса пластины ф = 0.
Функция mod приводит углы к диапазону [0, 2п), а функции ceil и floor дают диапазоны значений [1, Na] и [0, Na - 1]. Так как по определению значения индекса j должны лежать в диапазоне [0, Na - 1], то в случае ceil(mod^c - arcsin(r/pc), 2я)/Аф) = Na нижняя граница изменения индекса j принимается равной нулю, что описывается с помощью функции if.
При pc(t) < r условие (5) проверяется для всех значений индекса j, так как понятие угла охвата зоны обработки теряет в этом случае смысл и для индексов i, удовлетворяющих условию
0 < i < floor((pc + r)/Ар).
При моделировании рассматривались следующие варианты начальных условий:
pc (0) = R - r, фc (0) = я, p c (T) = r; (7) pc (0) = R, фc (0) = я, p c(T) = r; (8)
pc(0) = R, фc(0) = я, pc(T) = 0. (9)
Начальное условие фД0) = п удобно использовать, так как в этом случае угол фc - arcsin(r/pc) все время остается положительным.
Пластину будем условно считать неподвижной, а полюсный наконечник - совершающим радиальное движение подачи с линейной скоростью vr в сочетании с вращательным дви-
жением относительно оси симметрии пластины с угловой скоростью В этом случае:
Pc (,) = Pc (0) - vrt;
Фc (t) = фc (0) + Qt;
T = (pc(0)-Pc(T))/Уг .
При моделировании будет использоваться дискретная временная переменная:
= к At; к = 0... N; А = Т/Н,.
Давление р будет считаться постоянным в области проекции полюсного наконечника, однако в случае необходимости в рассматриваемую модель может быть без труда введена истинная функция распределения давления, определенная расчетным или экспериментальным путем. Скорость V движения точки Р поверхности заготовки относительно инструмента определялась по формуле
ур=7 ^ + у12,
где V! - составляющая скорости, перпендикулярная радиус-вектору оси вращения инструмента; Уп - составляющая скорости, параллельная радиус-вектору оси вращения инструмента.
Составляющие скорости рассчитывались по формулам:
• 008(6) + Орг ео8(фс (,к) - Ф;);
VI = юЛ ^ 1 -0082 (6) + + ^Рг Ф - 0°52(Фс (,к ) -ф1) + V,
где Л = д/р^ ) + Р2 - 2Рс (,к )Рг 008(Фс (,к ) - Ф] ) -расстояние между точкой поверхности заготовки и осью вращения инструмента; со8(6) = = (Рг 008(Фс(А)-Ф])-Рс(,к))/Л - косинус угла
между радиус-вектором оси вращения инструмента и вектором, соединяющим рассматриваемую точку поверхности заготовки с осью вращения инструмента.
Наука
итехника. Т. 22, № 6 (2023)
Для расчета функции (2) вводилась матрица {Ь2} съема припуска с нулевыми начальными значениями элементов. Если в момент времени tk точка PiJ- находится в области обработки, то значение соответствующего ей элемента матрицы {Ь2} увеличивается на величину vP(tk)At. Для определения истинных значений съема припуска полученную таким образом матрицу нужно умножить на коэффициент k, характеризующий условия обработки, и давление p в зоне обработке, однако для качественного анализа возникающих при обработке отклонений формы можно принять k = 1 и p = 1.
Результаты расчета при начальных условиях (7)-(9) приведены на рис. 1.
Расчет производился при следующих исходных данных: R = 150 мм; радиус полюсного наконечника r = 10 мм; Q = 1 об/с; ю = 10 об/с; vr = 1 мм/с; N = 150; N = 360; At = 0,56 мс.
Все рассмотренные случаи характеризуются значительной неравномерностью съема припуска. В случае начальных условий (7) наблюдается «краевой эффект», состоящий в снижении величины съема припуска до нуля в центре
а
1,5х103.-1-
50 100 150 Радиальная координата, мм
с
и на периферии заготовки, не пересекаемых инструментом в процессе обработки. В случае начальных условий (8) краевой эффект не периферии заготовки становится менее выраженным, так как инструмент пересекает ее в процессе обработки. Аналогичным образом при начальных условиях (9) снижается краевой эффект в центре заготовки. Кривые, аналогичные представленным на рис. 1, были ранее получены при исследовании имеющего схожую кинематику процесса напыления покрытий на поверхность вращающегося диска [9]. Таким образом, для исключения возникновения краевого эффекта инструмент должен пересекать при обработке центр и периферию заготовки. Что касается неравномерности съема припуска вне зоны краевого эффекта, она сохраняется во всех рассмотренных случаях и для ее компенсации требуется управление режимами обработки, которое может быть реализовано несколькими способами, вытекающими из уравнения (1), за счет изменения:
1) кинематических параметров процесса обработки (скорости V).
2) давления р.
Ь
50 100 150
Радиальная координата, мм
■ ; ■ <D
/1 20 1 1 130
w о s -& S - ! -
\ -
о
0 1 1
50 100 150
Радиальная координата, мм
Рис. 1. Результаты расчета функции съема припуска: а - для начального условия (7); b - для (8); с - для (9)
Fig. 1. Calculation results of material removal function: a - for initial condition (7); b - for (8); с - for (9)
Наука
итехника. Т. 22, № 6 (2023)
Предпочтение следует отдать второму способу, так как давление р влияет на скорость съема припуска линейным образом (в отличие от кинематических параметров ю и у), хотя существуют исследования, посвященные управлению съемом припуска при МАФО за счет изменения кинематических параметров [10]. Так как давление при МАФО определяется характеристиками магнитного поля [11], технически возможны следующие варианты его регулирования:
1) изменение индукции магнитного поля за счет регулирования силы тока, протекающего через обмотку электромагнита (для магнитных систем на основе электромагнитных индукторов);
2) варьирование величины воздушного зазора между полюсным наконечником и поверхностью заготовки (для магнитных систем на основе электромагнитных индукторов и постоянных магнитов);
3) изменение скважности импульсов магнитного поля при фиксированной частоте их повторения (широтно-импульсная модуляция) (для магнитных систем на основе электромагнитных индукторов при обработке в импульсных магнитных полях).
Обратная задача
моделирования процесса МАФО
Рассматриваемая в данном разделе обратная задача состоит в определении закона регулирования давления р при МАФО, обеспечивающего заданное (близкое к равномерному) распределение величины съема припуска по поверхности заготовки. Для вывода описывающих обратную задачу уравнений рассмотрим одномерную функцию съема припуска
Н10 (р) = Л(Р, 0)
и представим ее в виде суммы функций съема за один оборот заготовки (пооборотных функций съема):
ЕН+1 г * (к+1) х—\Ы+1
^ [(к) ДД(Р,0, *)Л = Хн=1/к (р), (10)
где N - число полных оборотов заготовки за время Т; * (к) = 2%(к -1)/О при к < N + 1; *(N + 2) = Т .
Функции /к(р) рассчитываются путем численного интегрирования функции (1) скорости съема припуска на интервале (*(к); *(к + 1)).
Как было показано в предыдущем разделе, при постоянных режимах обработки функция h1D(p) является непостоянной, что соответствует возникновению при обработке отклонений формы. Для минимизации этих отклонений необходимо, чтобы функция h1D(p) принимала близкое к постоянному значение, что может быть достигнуто заменой суммирования в выражении (10) взвешенным суммированием:
hD (p) = ZN=+1 wh (p), (11)
где весовые коэффициенты wk должны быть подобраны таким образом, чтобы обеспечить минимальное отклонение функции h1D(p) от заданного значения.
С технической точки зрения, введение весовых коэффициентов в выражении (11) соответствует пропорциональному изменению давления p при обработке. Таким образом, обратная задача моделирования процесса МАФО может быть сведена к определению весовых коэффициентов wk в выражении (11), обеспечивающих минимальное отклонение одномерной функции h1D(p) съема припуска от заданного постоянного значения. Для решения этой задачи представим выражение (11) в дискретной форме
hopt(i)=Е I+1 wkhi, k, (12)
где h^ = h0pt(pi); \к = hk (p); р. = (i-1)R (Nr-1), i = 1... Nr; hopt(p) - требуемое (оптимальное) распределение одномерной функции h1D(p) съема припуска.
Уравнение (12) может быть записано в матричной форме
{h}w = h opt, (13)
где {h} - матрица пооборотного съема припуска, столбцы которой соответствуют дискрети-зированным значениям пооборотных функций съема hk(p); w - вектор весовых коэффициентов; hopt - вектор оптимального распределения съема припуска.
Как было показано ранее (рис. 1), при перемещении инструмента из положения pc (0) = R - r в положение pc(T) = r функция h1D(p), рассчитанная для постоянных режимов обработки, характеризуется наличием зон краевого эффекта р > R - 2r и р < 2r, в которых величина съема припуска плавно снижается до нуля. Поэтому
Наука
итехника. Т. 22, № 6 (2023)
оптимальное распределение /^(р) съема припуска должно задаваться в виде трапецеидальной функции
Мр) =
/0 • > Я - 2г, (Я - Р)/2г, < 2г, р/2г , 1)).
Так как матрица {Ь} в (13) в общем случае является прямоугольной, для решения этого уравнения требуется использование обобщенной обратной (псевдообратной) матрицы Му-ра-Пенроуза:
w,
= min(H )/Hk < 1;
{h}+= ({h}T {h})-1{h}T ;
w = {h}+ h
opt'
(14)
где индекс T обозначает транспонирование матрицы.
Численный анализ показывает, что матрица {h}T{h} и соответственно рассматриваемая обратная задача моделирования процесса МАФО являются плохо обусловленными, что характерно для многих обратных задач. Как будет показано далее, представление решения обратной задачи в виде (14) дает достаточно точное (в смысле величины отклонения || {h}w -hopt ||,
где |||| - евклидова (12) норма) решение, однако распределение весовых коэффициентов wk является в этом случае недостаточно гладким, что делает затруднительным или невозможным соответствующее управление режимами обработки. Сглаживание решения и улучшение обусловленности задачи могут быть достигнуты путем регуляризации, по Тихонову, которая в простейшем варианте описывается уравнением [8]
w = ({h}T {h} + M)-1{h}T hopt,
где X > 0 - параметр регуляризации; I - единичная матрица.
Выбор параметра X является компромиссной задачей, так как его увеличение приводит к повышению гладкости решения, однако при этом происходит потеря точности (увеличение отклонения ||{h}w - hopt ||).
В работах [9, 12], посвященных равномерности нанесения покрытий на поверхность вращающихся деталей, использовалось следующее интуитивное (не имеющее строгого математического обоснования, однако демонстрирующее практически приемлемые результаты) определение весовых коэффициентов:
■ Наука
итехника. Т. 22, № 6 (2023)
Hk = max({h}<k > ),
(15)
(16)
где индекс <к> соответствует извлечению к-го столбца матрицы.
Технически подобное определение весовых коэффициентов соответствует снижению давления р при обработке обратно пропорционально максимальному значению пооборотного съема припуска, то есть на оборотах, соответствующих максимальному съему припуска будет прикладываться минимальное давление. Альтернативным вариантом является определение весовых коэффициентов через интегральное (или, что равносильно, среднее) значение пооборотного съема припуска. В этом случае
Hk =Е S
Н, к-
(17)
Оба описанных определения весовых коэффициентов будут использоваться в дальнейшем для оценки достоверности результатов, полученных путем решения регуляризованной обратной задачи.
На рис. 2 представлены графики изменения весовых коэффициентов, рассчитанных по формулам (15), (16) и (15), (17) для тех же исходных данных, что и приведенное выше решение прямой задачи.
1,0
M
о H 0,8
s и и
к M s я 0,б
я
а
к СО •е- 0,4
о
•Л, 0,2
i
;2 i \
i \ » % - \ i
50 100 Номер оборота
150
Рис. 2. Графики изменения весовых коэффициентов: 1 - расчет по формулам (15), (16); 2 - по (15), (17) Fig. 2. Plots of variation in weight coefficients:
1 - calculation by formulas (15), (16);
2 - calculation by formulas (15), (17)
На рис. 3 представлен график оптимального распределения съема припуска и графики съема припуска, рассчитанные по формуле (11) с учетом приведенных на рис. 2 значений весовых коэффициентов.
50 100 150
Радиальная координата, мм
400
(D g
О
Е? s а а
(D tfl О
200
100
150
Радиальная координата, мм
Рис. 3. Графики: а - оптимального распределения съема; b - рассчитанные по формуле (11); 1 - расчет весовых коэффициентов по (15), (16); 2 - то же по (15), (17) Fig. 3. Plots of material removal: a - optimal distribution of removal; b - plots calculated by formula (11); 1 - calculation of weight coefficients by formulas (15), (16); 2 - calculation of weight coefficients by formulas (15), (17)
Как видно из графиков, распределение съема припуска становится значительно более равномерным по сравнению с распределениями, представленными на рис. 1, однако все же достаточно сильно отклоняется от оптимального, особенно в случае расчета весовых коэффициентов по (15), (16).
На рис. 4 представлены графики изменения весовых коэффициентов, рассчитанных путем решения нерегуляризованной и регуляризован-ной обратных задач.
На этих же графиках для сравнения приведено нормированное распределение весовых коэффициентов, рассчитанных по (15), (17). Нормировка производилась путем умножения на постоянный коэффициент, минимизирую-
8x10
«
о
H
« и
s (D
« S
<D я
S s
л и го -е
-е
m §
6x10
4x10
2x10
-2x10
50 100 Номер оборота
150
щий среднеквадратическое отклонение между коэффициентами, рассчитанными двумя способами. Как видно из графиков, общая тенденция изменения весовых коэффициентов соответствует расчету по (15), (17), однако для изменения коэффициентов, рассчитанных путем решения обратных задач, характерна более высокая степень колебательности. Как предсказывает теория, решение регуляризованной обратной задачи является более гладким по сравнению с решением нерегуляризованной задачи.
На рис. 5 представлены графики съема припуска, рассчитанные по формуле (11) с учетом значений весовых коэффициентов, полученных путем решения обратных задач.
Ь
5x10
4x10 "
3x10
2x10 ""
1x10 "
50 100 Номер оборота
150
Рис. 4. Графики изменения весовых коэффициентов: а - расчет путем решения нерегуляризованной обратной задачи; b - то же путем решения регуляризованной обратной задачи; 1 - расчет по формулам (15), (17); 2 - решение обратной задачи
Fig. 4. Plots of variation in weight coefficients: a - calculation using solution of non-regularized inverse problem; b - calculation using solution of regularized inverse problem; 1 - calculation by formulas (15), (17); 2 - solution of inverse problem
Наука
итехника. Т. 22, № 6 (2023)
b
a
a
1,5
nS w
о
1,0
& 0,5
<D fi О
_L
0 50 100
Радиальная координата, мм
1,05
7- 0,95,
50 100
150
1,5
1,05
S
о
1,0
& 0,5
<D tfl О
0 50 100 150 50 100 150
Радиальная координата, мм
Рис. 5. Графики съема припуска: а - расчет весовых коэффициентов путем решения нерегуляризованной обратной задачи; b - то же регуляризованной обратной задачи
Fig. 5. Plots of material removal: a - calculation of weight coefficients using solution of non-regularized inverse problem; b - calculation of weight coefficients using solution of regularized inverse problem
a
b
В обоих случаях (для нерегуляризованной и регуляризованной обратных задач) расчетные распределения съема припуска достаточно близки к требуемому оптимальному распределению. Как предсказывает теория, для решения регуляризованной задачи характерно более сильное отклонение расчетного распределения от требуемого оптимального.
На практике важной является не только задача о равномерном съеме припуска, но и более общая задача о съеме припуска по определенному закону, так как заготовка поступает на операцию МАФО с отклонениями от формы, возникшими на предшествующих операциях обработки, и соответственно возникает проблема коррекции этих отклонений (корректирующей обработки со съемом припуска по определенному закону) [10, 13]. Решение этой задачи будет рассмотрено в дальнейших исследованиях. Корректирующая обработка может производиться как методом МАФО [10], так и с помощью других методов, таких как локальное жидкостное травление [13], однако последний
■ Наука
итехника. Т. 22, № 6 (2023)
способ требует применения химически агрессивных сред.
ВЫВОДЫ
1. Разработана методика решения прямых и обратных задач моделирования процесса МАФО полупроводниковых пластин, позволяющих рассчитывать функцию съема припуска при заданных параметрах обработки (прямая задача) и требуемые для реализации заданной функции съема припуска значения параметров обработки (обратная задача).
2. На основе анализа прямой задачи с постоянными значениями кинематических и магнитных параметров обработки показано, что обработка при постоянных значениях параметров не обеспечивает равномерности съема припуска.
3. На основе анализа обратной задачи показано, что управление магнитными параметрами обработки позволяет с достаточно высокой точностью реализовать заданную функцию съема припуска и, как частный случай, близкий к равномерному съем припуска.
4. Путем сравнительного анализа решения исходной обратной задачи и решения обратной регуляризованной, по Тихонову, задачи показано, что регуляризация позволяет повысить гладкость решения и соответственно упростить процесс управления параметрами обработки.
ЛИТЕРАТУРА
1. Хомич, Н. С. Магнитно-абразивная обработка изделий / Н. С. Хомич. Минск: БНТУ, 2006. 217 с.
2. Mosavat, M. Nano-Finishing of the Monocrystalline Silicon Wafer Using Magnetic Abrasive Finishing Process / M. Mosavat // Applied Optics. 2019. Vol. 58, № 13. P. 3447-3453. https://doi.org/10.1364/A0.58.003447.
3. The use of Potassium Peroxidisulphate and Oxone® as Oxidizers for the Chemical Mechanical Polishing of Silicon wafers / A. Pineiro [et al.] // Wear. 2013. Vol. 303, № 1-2. P. 446-450. https://doi.org/10.1016/j.wear.2013. 03.030.
4. An Efficient Approach for Atomic-Scale Polishing of Single-Crystal Silicon Via Plasma-Based Atom-Selective Etching / Z. Fang [et al.] // International Journal of Machine Tools and Manufacture. 2020. Vol. 159, Part B. Article 103649. https://doi.org/ 10.1016/j.ijmachtools. 2020. 103649.
5. Influence of Silicon Wafer Surface Roughness on Semiconductor Device Characteristics / K. Mori [et al.] // Japanese Journal of Applied Physics. 2020. Vol. 59. Article SMMB06. https://doi.org/10.35848/1347-4065/ab918c.
6. Silicon Wafer Thickness Variation Measurements Using the National Institute of Standards and Technology Infrared Interferometer / T. L. Schmitz [et al.] // Optical Engineering. 2003. Vol. 42, N 8. P. 2281-2290. https://doi. org/10.1117/1.1589757.
7. Моделирование кинематики процесса магнитно-абразивного полирования кремниевых пластин / Н. С. Хомич [и др.] // Вестник БНТУ. 2009. № 1. С. 32-38.
8. Bouvier, C. Investigation of Polishing Algorithms and Removal Processes for a Deterministic Subaperture Polisher: PhD Thesis / C. Bouvier. - University of Rochester, 2007. 272 p.
9. Stepanenko, D. A. Modeling of Spraying with Time-Dependent Material Feed Rate / D. A. Stepanenko // Applied Mathematical Modelling. 2007. Vol. 31, № 11. P. 2564-2576. https://doi.org/10.1016/j.apm.2006.10.005.
10. Zhang, Y. Study of Corrective Abrasive Finishing for Plane Surfaces Using Magnetic Abrasive Finishing Processes / Y. Zhang, Y. Zou // Nanotechnology and Precision Engineering. 2021. Vol. 4, № 3. Article 033001. https://doi.org/10.1063/10.0004961.
11. Modeling of Material Removal in Magnetic Finishing Based on Maxwell's Stress Tensor Theory and its Experimental Validation / X. He [et al.] // Journal of Materials Processing Technology. 2023. Vol. 312. Article 117808. https://doi.org/10. 1016/j.jmatprotec.2022.117808.
12. Киселев, М. Г. Исследование равномерности нанесения покрытий на примере напыления абразивной суспензии на поверхность ограночного диска / М. Г. Киселев, Д. А. Степаненко // Вестник БНТУ. 2006. № 3. С. 51-56.
13. Improvement of Thickness Uniformity of Bulk Silicon Wafer by Numerically Controlled Local Wet Etching / M. Nagano [et al.] // Journal of Crystal Growth. 2009. Vol. 311, N 8. P. 2560-2563. https://doi.org/10.1016/). jcrysgro.2009.01.131.
Поступила 10.07.2023 Подписана в печать 12.09.2023 Опубликована онлайн 30.11.2023
REFERENCES
1. Khomich N. S. (2006) Magnetic Abrasive Treatment of products. Minsk, Belarusian National Technical University. 217 (in Russian).
2. Mosavat M. (2019) Nano-Finishing of the Monocrystalline Silicon Wafer Using Magnetic Abrasive Finishing Process. Applied Optics, 58 (13), 3447-3453. https://doi.org/ I0.1364/A0.58.003447.
3. Pineiro A., Black A., Medina J. C., Dieguez E., Parra V. (2013) The use of Potassium Peroxidisulphate and Oxone® as Oxidizers for the Chemical Mechanical Polishing of Silicon Wafers. Wear, 303 (1-2), 446-450. https://doi.org/ 10.1016/j.wear.2013.03.030.
4. Fang Z., Zhang Y., Li R., Liang Y., Deng H. (2020) An Efficient Approach for Atomic-Scale Polishing of Single-Crystal Silicon Via Plasma-Based Atom-Selective Etching. International Journal of Machine Tools and Manufacture, 159, Part B, 103649. https://doi.org/10.1016/). ijmachtools.2020.103649.
5. Mori K., Samata S., Mitsugi N., Teramoto A., Kuroda R., Suwa T., Hashimoto K., Sugawa S. (2020) Influence of Silicon Wafer Surface Roughness on Semiconductor Device Characteristics. Japanese Journal of Applied Physics, 59, Article SMMB06. https://doi.org/10.35848/1347-4065/ab918c.
6. Schmitz T. L., Davies A., Evans C. J., Parks R. E. (2003) Silicon Wafer Thickness Variation Measurements Using the National Institute of Standards and Technology Infrared Interferometer. Optical Engineering, 42 (8), 2281-2290. https://doi.org/10.1117/L1589757.
7. Fedortsev R. V., Khomich N. S., Lougovik A. I., Kor-zun A. E., Kukhto P. V. (2009) Kinematics Modeling of Magnetic and Abrasive Silicon Wafer Polishing. Vest-nik BNTU, (1), 32-38 (in Russian).
8. Bouvier C. (2007) Investigation of Polishing Algorithms and Removal Processes for a Deterministic Subaperture Polisher: PhD Thesis. University of Rochester. 272.
9. Stepanenko D. A. (2007) Modeling Of Spraying With Time-Dependent Material Feed Rate. Applied Mathemati-calModelling, 31 (11), 2564-2576. https://doi.org/10.1016/). apm.2006.10.005.
10. Zhang Y., Zou Y. (2021) Study of Corrective Abrasive Finishing for Plane Surfaces Using Magnetic Abrasive Finishing Processes. Nanotechnology and Precision Engineering, 4 (3), Article 033001. https://doi.org/10.1063/10.0004961.
11. He X., Jin H., Zhou C., Gao C., Zhang G., Shiju E. (2023) Modeling of Material Removal in Magnetic Finishing Based on Maxwell's Stress Tensor Theory and its Experimental Validation. Journal of Materials Processing Technology, 312, Article 117808. https://doi.org/10.1016/). jmatprotec.2022.117808.
12. Kiselev M. G., Stepanenko D. A. (2006) Study of Uniformity of Coating Application Using the Example of Spraying an Abrasive Suspension Onto the Surface of a Cutting Disc. Vestnik BNTU, (3), 51-56 (in Russian).
13. Nagano M., Mitani T., Ueda K., Zettsu N., Yamamura K. (2009) Improvement of Thickness Uniformity of Bulk Silicon Wafer by Numerically Controlled Local Wet Etching Journal of Crystal Growth, 311 (8), 2560-2563. https://doi.org/10.1016/j.jcrysgro.2009.01.131.
Received: 10.07.2023 Accepted: 12.09.2023 Published online: 30.11.2023
Наука
итехника. Т. 22, № 6 (2023)