Выпуск #8/2023
С. Ванцов, О. Хомутская, Е. Лийн
ВЛИЯНИЕ КОНСТРУКТИВНЫХ ПАРАМЕТРОВ НА ПЛОСКУЮ ДЕФОРМАЦИЮ ПЕЧАТНЫХ ПЛАТ
ВЛИЯНИЕ КОНСТРУКТИВНЫХ ПАРАМЕТРОВ НА ПЛОСКУЮ ДЕФОРМАЦИЮ ПЕЧАТНЫХ ПЛАТ
Просмотры: 459
DOI: 10.22184/1992-4178.2023.229.8.108.112
Описываются построение математической модели плоской деформации слоев печатной платы и экспериментальная проверка ее правильности. Полученные экспериментальные результаты качественно полностью совпадают с результатами расчета; количественное расхождение не превышает 20%, что можно считать удовлетворительным на данном этапе уточнения коэффициентов модели с учетом случайных факторов.
Описываются построение математической модели плоской деформации слоев печатной платы и экспериментальная проверка ее правильности. Полученные экспериментальные результаты качественно полностью совпадают с результатами расчета; количественное расхождение не превышает 20%, что можно считать удовлетворительным на данном этапе уточнения коэффициентов модели с учетом случайных факторов.
Теги: drilling etching fiberglass laminate plane strain plane strain model pressing printed circuit boards модель плоской деформации печатные платы плоская деформация прессование сверление стеклотекстолит травление
Влияние конструктивных параметров на плоскую деформацию печатных плат
С. Ванцов, к. т. н., О. Хомутская, к. т. н., Е. Лийн
В настоящее время общим местом является обсуждение актуальности исследований в области технологии изготовления различного рода печатных плат. Увеличение плотности монтажа и повышение класса точности печатных плат приводит к увеличению требований к точности позиционирования элементов печатной платы (ПП). Особенно важным это требование становится в производстве многослойных плат, в которых несовмещение слоев может привести к неисправимому производственному браку. Последнее связано не столько с погрешностями выполнения отдельных технологических операций, сколько с характеристическими свойствами материала диэлектрического основания слоев печатной платы.
Традиционным, наиболее часто используемым, материалом диэлектрического основания является стеклотекстолит, представляющий собой слоистую структуру из стеклоткани и различного вида смол, чаще всего эпоксидной смолы. Этим определяются и физико-механические характеристики таких материалов. Наиболее существенными из них являются анизотропия свойств по объему материала и зависимость твердости от температуры. Образно говоря, контактные площадки печатной платы в силу этих свойств расположены на льдинах во время ледохода, случайно перемещающихся под влиянием изменения температуры как окружающего воздуха, так и воды, в которой льдины находятся. Такое сравнение тем более наглядно, поскольку пластические материалы, используемые в диэлектрическом основании, по сути представляют собой жидкости с достаточно высокой температурой текучести. Взаиморасположение «льдин» и их форма в плоскости «поверхности воды» определяются связанностью отдельных участков поверхности печатной платы за счет структуры плетения стеклоткани и топологии печатного рисунка, а по объему – с температурным воздействием на толщину «льдины».
Из сказанного следует, что, с учетом первоначально фиксированного общего ледяного покрова поверхности, величина перемещений участков печатной платы друг относительно друга зависит от:
случайного сочетания факторов, удерживающих ледяной пласт от раскалывания;
состояния среды, в которой происходит раскол участков льдин;
внешних воздействий, влияющих на скорость их расползания.
Для печатных плат случайное сочетание факторов зависит от физико-механических свойств материала печатной платы.
Состояние среды зависит от параметров технологических операций, производимых с заготовками в процессе изготовления ПП.
Внешние воздействия в значительной мере определяются состоянием плоских поверхностей слоя печатной платы, а именно количеством меди на отдельных фрагментах поверхности слоя, которое зависит от топологии платы.
Изначально топология ПП задается конструктором исходя из целесообразности размещения дискретных элементов на плате и требований по построению эффективной системы передачи сигналов внутри реализуемой электрической схемы. Отсюда возникает задача предсказания изменения расположения отдельных участков печатной платы в ходе осуществления каждой технологической операции и компенсации этих изменений на этапе проектирования и конструирования печатной платы.
На деформацию многослойной печатной платы в целом оказывают влияние следующие технологические операции:
Сверление базовых и переходных отверстий;
Травление слоев для получения рисунка печатной платы;
Прессование пакета слоев платы.
Влияние этих операций на величину плоской деформации печатной платы определяется следующими действиями:
Процесс сверления заключается в удалении части армирующего поверхностного слоя меди и части тела диэлектрического основания, что в целом ослабляет жесткость слоя платы и, отчасти, способствует появлению релаксирующих усилий за счет снижения степени внутренних напряжений материала основания.
Процесс травления имеет целью получение печатного рисунка и заключается в удалении значительной части армирующего слоя меди, что в еще большей степени приводит к явлениям снижения жесткости слоя и увеличению релаксирующих воздействий.
Процесс прессования слоев оказывает наибольшее влияние на деформацию слоя с уже ослабленной на предыдущих операциях жесткостью за счет внешних воздействий – силового и теплового. Последнее в наибольшей степени проявляется с учетом физико-механических свойств слоистых пластиков на базе термореактивных полимеров. На рис. 1 показан типичный вид изменения линейных размеров пластиков при тепловом воздействии с учетом возникающих в их материале фазовых переходов.
Визуально процесс плоской деформации слоя печатной платы можно показать, используя в качестве «маячков» величины смещения отдельных участков изменения положений центров сквозных отверстий на слое (рис. 2).
Начальным полагается расположение центра отверстия до проведения каждой технологической операции, конечным – после ее проведения. Величина смещения рассчитывается как векторная сумма и для каждого отверстия.
Исходя из этой предпосылки, величину смещения центра произвольного исходного отверстия после проведения технологических операций можно рассчитать, используя следующую математическую модель [1]:
,
, (1)
где – уточненный вектор деформации в квадрате разбиения по оси абсцисс;
– уточненный вектор деформации в квадрате разбиения по оси ординат;
– исходный вектор деформации в квадрате разбиения по оси абсцисс;
– исходный вектор деформации в квадрате разбиения по оси ординат;
N – количество квадратов разбиения в строке;
i, j – индексы квадрата разбиения, деформация которого определяется;
l, k – индексы промежуточных квадратов разбиения, влияющих на деформацию квадрата с индексами i и j;
– поправочный коэффициент изменения деформации из-за топологии слоя и базового материала.
Эффективность построенной модели была экспериментально подтверждена при условии регулярного расположения исходных отверстий на поверхности основания платы [2]. Модель позволяет рассчитывать математическое ожидание величины плоской деформации при осуществлении различных технологических операций.
Однако при произвольном количестве и расположении отверстий на слое существенное влияние на результаты расчета приобретает значение поправочного коэффициента К. Его значение определяется преимущественно количеством просверленных отверстий в расчетном квадрате разбиения модели. С увеличением количества отверстий уменьшается количество базового материала, что естественным образом ведет к уменьшению жесткости основания платы и создает предпосылку к увеличению степени деформации слоя при проведении технологических операций.
В качестве параметров, характеризующих потенциал воздействия отверстий на деформацию на том или ином участке ПП, можно выделить:
расстояние между отверстиями;
диаметр отверстий.
Поскольку на реальных платах отверстия в основном располагаются не матрицей, возникает необходимость перехода к оценке деформации при случайном расположении отверстий (рис. 3).
В этом случае целесообразным является введение интегральной оценки, оказывающей влияние на величину весового коэффициента в модели (1), учитывающей случайное расположение отверстий. В качестве такой оценки можно использовать удельную плотность сверления:
(2)
где Soi – площадь i-го отверстия, n – количество отверстий, Sпп – площадь рассматриваемого участка печатной платы (рис. 4).
Для учета такого конструктивного параметра платы, как толщина диэлектрического слоя, можно также воспользоваться удельной объемной плотностью сверления
(3)
где – объем удаляемого материала при сверлении i-ого отверстия; n – количество отверстий на исследуемом участке слоя платы; Vпп – объем материала на исследуемом участке; – глубина просверливаемого отверстия, в том числе и глухого; hпп – толщина слоя платы.
Введенные оценки удельной плотности сверления позволяют сделать заключение о том, что величина весового коэффициента будет лежать в пределах от 0 до 1. При отсутствии отверстий в данном i-ом квадрате разбиения K = 0 и деформации на этом участке отсутствуют, то есть Li, j = Li ‑ 1, j − 1. В предельном случае, при отсутствии материала на данном участке разбиения K = 1; такая ситуация может возникнуть в случае сложного внешнего контура границ слоя.
Предложенный подход позволяет экспериментально выявить зависимости деформации участков слоя печатной платы, а именно – зависимости математического ожидания от диаметра сквозного отверстия M(D) от расстояния между сквозными отверстиями M(l) и от толщины печатной платы M(h), а также среднеквадратического отклонения от тех же величин: σ(D), σ(l) и σ(h).
Принцип выявления таких зависимостей базируется на анализе наложений двух изображений – анализируемого и эталонного. В качестве реперных точек для получения расчетных данных при определении величины деформации используются центры сквозных отверстий. Разница в смещении реперных точек на эталонном и анализируемом изображении позволяет определить в целом величину плоской деформации слоя, возникшей в результате осуществления определенной технологической операции.
Эксперимент включает в себя следующие этапы:
Получение изображения поверхности слоя до осуществления технологической операции.
Получение изображения поверхности слоя после осуществления технологической операции.
Определение центров отверстий на изображениях.
Приведение центров отверстий на сравниваемых изображениях к единой системе координат.
Расчет величин смещения центров отверстий.
Построение регрессионной модели деформации слоя.
Расчет величины смещения центров отверстий на эталонном и анализируемом изображениях производится по следующей формуле:
(4)
где – величина в пикселях смещения центра i-го отверстия на анализируемом изображении по сравнению с эталонным; , – координаты центра i-го отверстия на эталонном изображении; , – координаты положения i-го отверстия на анализируемом изображении.
Результат расчета величины смещения сквозного отверстия переводится из пикселей в миллиметры по следующей формуле:
(5)
где li – величина смещения центра i-го отверстия в мм; DPI (dots per inch) – разрешение изображения в точках на дюйм.
Эксперимент по оценке влияния технологической операции сверления отверстий в печатной плате был проведен на 180 заготовках с более чем 100 тысячами успешно распознанных на изображениях отверстий.
Результаты оценки математического ожидания величины деформации в зависимости от расположения отверстий и их диаметра приведены на рис. 5.
На основании этой диаграммы можно сделать выводы:
при увеличении диаметра сверла математическое ожидание модуля длины вектора смещения отверстий увеличивается, плата деформируется сильнее;
при увеличении расстояния между отверстиями математическое ожидание модуля длины вектора смещения отверстий уменьшается, плата деформируется меньше.
Экспериментально полученные значения СКО значений деформации представляют собой практически линейные ниспадающие зависимости от расстояния между центрами отверстий и возрастающие от диаметров отверстий.
В результате полнофакторного эксперимента была построена регрессионная модель зависимости величины деформации диэлектрического основания от диаметров просверливаемых отверстий и расстояния между отверстиями на исследуемом квадрате разбиения поверхности.
. (6)
Зачастую в квадрате разбиения имеются отверстия разного диаметра и с разными расстояниями между ними. В этом случае удобным является использование коэффициента удельной плотности сверления (2) и коэффициента удельной объемной плотности сверления (3).
Характерная зависимость математического ожидания деформации от толщины диэлектрического слоя представлена на рис. 6.
Полученные экспериментальные результаты качественно полностью совпадают с результатами расчета, что позволяет сделать заключение о правильности построенной модели плоскостной деформации слоев печатной платы. Количественное расхождение не превышает 20%, что можно считать удовлетворительным на данном этапе уточнения коэффициентов модели с учетом случайных факторов.
Работа выполнена в рамках гранта Президента Российской Федерации для государственной поддержки молодых российских ученых – кандидатов наук (МК‑582.2022.4).
ЛИТЕРАТУРА
Никитин С., Поздняков К., Хомутская О. Оценка деформации печатных плат // ЭЛЕКТРОНИКА: Наука, Технология, Бизнес. 2019. № 5. С. 144–150.
Куликов Н., Хомутская О., Ванцов С. Цифровой метод автоматизированной оценки деформации печатной платы // ЭЛЕКТРОНИКА: Наука, Технология, Бизнес. 2018. № 2. С. 186–191.
С. Ванцов, к. т. н., О. Хомутская, к. т. н., Е. Лийн
В настоящее время общим местом является обсуждение актуальности исследований в области технологии изготовления различного рода печатных плат. Увеличение плотности монтажа и повышение класса точности печатных плат приводит к увеличению требований к точности позиционирования элементов печатной платы (ПП). Особенно важным это требование становится в производстве многослойных плат, в которых несовмещение слоев может привести к неисправимому производственному браку. Последнее связано не столько с погрешностями выполнения отдельных технологических операций, сколько с характеристическими свойствами материала диэлектрического основания слоев печатной платы.
Традиционным, наиболее часто используемым, материалом диэлектрического основания является стеклотекстолит, представляющий собой слоистую структуру из стеклоткани и различного вида смол, чаще всего эпоксидной смолы. Этим определяются и физико-механические характеристики таких материалов. Наиболее существенными из них являются анизотропия свойств по объему материала и зависимость твердости от температуры. Образно говоря, контактные площадки печатной платы в силу этих свойств расположены на льдинах во время ледохода, случайно перемещающихся под влиянием изменения температуры как окружающего воздуха, так и воды, в которой льдины находятся. Такое сравнение тем более наглядно, поскольку пластические материалы, используемые в диэлектрическом основании, по сути представляют собой жидкости с достаточно высокой температурой текучести. Взаиморасположение «льдин» и их форма в плоскости «поверхности воды» определяются связанностью отдельных участков поверхности печатной платы за счет структуры плетения стеклоткани и топологии печатного рисунка, а по объему – с температурным воздействием на толщину «льдины».
Из сказанного следует, что, с учетом первоначально фиксированного общего ледяного покрова поверхности, величина перемещений участков печатной платы друг относительно друга зависит от:
случайного сочетания факторов, удерживающих ледяной пласт от раскалывания;
состояния среды, в которой происходит раскол участков льдин;
внешних воздействий, влияющих на скорость их расползания.
Для печатных плат случайное сочетание факторов зависит от физико-механических свойств материала печатной платы.
Состояние среды зависит от параметров технологических операций, производимых с заготовками в процессе изготовления ПП.
Внешние воздействия в значительной мере определяются состоянием плоских поверхностей слоя печатной платы, а именно количеством меди на отдельных фрагментах поверхности слоя, которое зависит от топологии платы.
Изначально топология ПП задается конструктором исходя из целесообразности размещения дискретных элементов на плате и требований по построению эффективной системы передачи сигналов внутри реализуемой электрической схемы. Отсюда возникает задача предсказания изменения расположения отдельных участков печатной платы в ходе осуществления каждой технологической операции и компенсации этих изменений на этапе проектирования и конструирования печатной платы.
На деформацию многослойной печатной платы в целом оказывают влияние следующие технологические операции:
Сверление базовых и переходных отверстий;
Травление слоев для получения рисунка печатной платы;
Прессование пакета слоев платы.
Влияние этих операций на величину плоской деформации печатной платы определяется следующими действиями:
Процесс сверления заключается в удалении части армирующего поверхностного слоя меди и части тела диэлектрического основания, что в целом ослабляет жесткость слоя платы и, отчасти, способствует появлению релаксирующих усилий за счет снижения степени внутренних напряжений материала основания.
Процесс травления имеет целью получение печатного рисунка и заключается в удалении значительной части армирующего слоя меди, что в еще большей степени приводит к явлениям снижения жесткости слоя и увеличению релаксирующих воздействий.
Процесс прессования слоев оказывает наибольшее влияние на деформацию слоя с уже ослабленной на предыдущих операциях жесткостью за счет внешних воздействий – силового и теплового. Последнее в наибольшей степени проявляется с учетом физико-механических свойств слоистых пластиков на базе термореактивных полимеров. На рис. 1 показан типичный вид изменения линейных размеров пластиков при тепловом воздействии с учетом возникающих в их материале фазовых переходов.
Визуально процесс плоской деформации слоя печатной платы можно показать, используя в качестве «маячков» величины смещения отдельных участков изменения положений центров сквозных отверстий на слое (рис. 2).
Начальным полагается расположение центра отверстия до проведения каждой технологической операции, конечным – после ее проведения. Величина смещения рассчитывается как векторная сумма и для каждого отверстия.
Исходя из этой предпосылки, величину смещения центра произвольного исходного отверстия после проведения технологических операций можно рассчитать, используя следующую математическую модель [1]:
,
, (1)
где – уточненный вектор деформации в квадрате разбиения по оси абсцисс;
– уточненный вектор деформации в квадрате разбиения по оси ординат;
– исходный вектор деформации в квадрате разбиения по оси абсцисс;
– исходный вектор деформации в квадрате разбиения по оси ординат;
N – количество квадратов разбиения в строке;
i, j – индексы квадрата разбиения, деформация которого определяется;
l, k – индексы промежуточных квадратов разбиения, влияющих на деформацию квадрата с индексами i и j;
– поправочный коэффициент изменения деформации из-за топологии слоя и базового материала.
Эффективность построенной модели была экспериментально подтверждена при условии регулярного расположения исходных отверстий на поверхности основания платы [2]. Модель позволяет рассчитывать математическое ожидание величины плоской деформации при осуществлении различных технологических операций.
Однако при произвольном количестве и расположении отверстий на слое существенное влияние на результаты расчета приобретает значение поправочного коэффициента К. Его значение определяется преимущественно количеством просверленных отверстий в расчетном квадрате разбиения модели. С увеличением количества отверстий уменьшается количество базового материала, что естественным образом ведет к уменьшению жесткости основания платы и создает предпосылку к увеличению степени деформации слоя при проведении технологических операций.
В качестве параметров, характеризующих потенциал воздействия отверстий на деформацию на том или ином участке ПП, можно выделить:
расстояние между отверстиями;
диаметр отверстий.
Поскольку на реальных платах отверстия в основном располагаются не матрицей, возникает необходимость перехода к оценке деформации при случайном расположении отверстий (рис. 3).
В этом случае целесообразным является введение интегральной оценки, оказывающей влияние на величину весового коэффициента в модели (1), учитывающей случайное расположение отверстий. В качестве такой оценки можно использовать удельную плотность сверления:
(2)
где Soi – площадь i-го отверстия, n – количество отверстий, Sпп – площадь рассматриваемого участка печатной платы (рис. 4).
Для учета такого конструктивного параметра платы, как толщина диэлектрического слоя, можно также воспользоваться удельной объемной плотностью сверления
(3)
где – объем удаляемого материала при сверлении i-ого отверстия; n – количество отверстий на исследуемом участке слоя платы; Vпп – объем материала на исследуемом участке; – глубина просверливаемого отверстия, в том числе и глухого; hпп – толщина слоя платы.
Введенные оценки удельной плотности сверления позволяют сделать заключение о том, что величина весового коэффициента будет лежать в пределах от 0 до 1. При отсутствии отверстий в данном i-ом квадрате разбиения K = 0 и деформации на этом участке отсутствуют, то есть Li, j = Li ‑ 1, j − 1. В предельном случае, при отсутствии материала на данном участке разбиения K = 1; такая ситуация может возникнуть в случае сложного внешнего контура границ слоя.
Предложенный подход позволяет экспериментально выявить зависимости деформации участков слоя печатной платы, а именно – зависимости математического ожидания от диаметра сквозного отверстия M(D) от расстояния между сквозными отверстиями M(l) и от толщины печатной платы M(h), а также среднеквадратического отклонения от тех же величин: σ(D), σ(l) и σ(h).
Принцип выявления таких зависимостей базируется на анализе наложений двух изображений – анализируемого и эталонного. В качестве реперных точек для получения расчетных данных при определении величины деформации используются центры сквозных отверстий. Разница в смещении реперных точек на эталонном и анализируемом изображении позволяет определить в целом величину плоской деформации слоя, возникшей в результате осуществления определенной технологической операции.
Эксперимент включает в себя следующие этапы:
Получение изображения поверхности слоя до осуществления технологической операции.
Получение изображения поверхности слоя после осуществления технологической операции.
Определение центров отверстий на изображениях.
Приведение центров отверстий на сравниваемых изображениях к единой системе координат.
Расчет величин смещения центров отверстий.
Построение регрессионной модели деформации слоя.
Расчет величины смещения центров отверстий на эталонном и анализируемом изображениях производится по следующей формуле:
(4)
где – величина в пикселях смещения центра i-го отверстия на анализируемом изображении по сравнению с эталонным; , – координаты центра i-го отверстия на эталонном изображении; , – координаты положения i-го отверстия на анализируемом изображении.
Результат расчета величины смещения сквозного отверстия переводится из пикселей в миллиметры по следующей формуле:
(5)
где li – величина смещения центра i-го отверстия в мм; DPI (dots per inch) – разрешение изображения в точках на дюйм.
Эксперимент по оценке влияния технологической операции сверления отверстий в печатной плате был проведен на 180 заготовках с более чем 100 тысячами успешно распознанных на изображениях отверстий.
Результаты оценки математического ожидания величины деформации в зависимости от расположения отверстий и их диаметра приведены на рис. 5.
На основании этой диаграммы можно сделать выводы:
при увеличении диаметра сверла математическое ожидание модуля длины вектора смещения отверстий увеличивается, плата деформируется сильнее;
при увеличении расстояния между отверстиями математическое ожидание модуля длины вектора смещения отверстий уменьшается, плата деформируется меньше.
Экспериментально полученные значения СКО значений деформации представляют собой практически линейные ниспадающие зависимости от расстояния между центрами отверстий и возрастающие от диаметров отверстий.
В результате полнофакторного эксперимента была построена регрессионная модель зависимости величины деформации диэлектрического основания от диаметров просверливаемых отверстий и расстояния между отверстиями на исследуемом квадрате разбиения поверхности.
. (6)
Зачастую в квадрате разбиения имеются отверстия разного диаметра и с разными расстояниями между ними. В этом случае удобным является использование коэффициента удельной плотности сверления (2) и коэффициента удельной объемной плотности сверления (3).
Характерная зависимость математического ожидания деформации от толщины диэлектрического слоя представлена на рис. 6.
Полученные экспериментальные результаты качественно полностью совпадают с результатами расчета, что позволяет сделать заключение о правильности построенной модели плоскостной деформации слоев печатной платы. Количественное расхождение не превышает 20%, что можно считать удовлетворительным на данном этапе уточнения коэффициентов модели с учетом случайных факторов.
Работа выполнена в рамках гранта Президента Российской Федерации для государственной поддержки молодых российских ученых – кандидатов наук (МК‑582.2022.4).
ЛИТЕРАТУРА
Никитин С., Поздняков К., Хомутская О. Оценка деформации печатных плат // ЭЛЕКТРОНИКА: Наука, Технология, Бизнес. 2019. № 5. С. 144–150.
Куликов Н., Хомутская О., Ванцов С. Цифровой метод автоматизированной оценки деформации печатной платы // ЭЛЕКТРОНИКА: Наука, Технология, Бизнес. 2018. № 2. С. 186–191.
Отзывы читателей