eng
Выпуск #2/2023
О. Почтарь, А. Почтарь
ОЦЕНКА СТАБИЛЬНОСТИ И ВОСПРОИЗВОДИМОСТИ ТЕХНОЛОГИЧЕСКОГО ПРОЦЕССА МЕТАЛЛИЗАЦИИ ЧЕРЕЗ СВОБОДНУЮ МАСКУ
ОЦЕНКА СТАБИЛЬНОСТИ И ВОСПРОИЗВОДИМОСТИ ТЕХНОЛОГИЧЕСКОГО ПРОЦЕССА МЕТАЛЛИЗАЦИИ ЧЕРЕЗ СВОБОДНУЮ МАСКУ
Просмотры: 691
DOI: 10.22184/1992-4178.2023.223.2.72.76
В статье на основе статистической обработки экспериментальных данных представлена оценка стабильности и воспроизводимости технологического процесса двухслойной металлизации кристаллов ниобата лития через свободную металлическую маску методом электронно-лучевого испарения.
В статье на основе статистической обработки экспериментальных данных представлена оценка стабильности и воспроизводимости технологического процесса двухслойной металлизации кристаллов ниобата лития через свободную металлическую маску методом электронно-лучевого испарения.
Теги: electron beam evaporation process stability and reproducibility shewhart control charts statistical process control two-layer crystal metallization двухслойная металлизация кристаллов контрольные карты шухарта стабильность и воспроизводимость технологического процесса статистический контроль технологических процессов электронно-лучевое испарение
Оценка стабильности и воспроизводимости
технологического процесса металлизации через свободную маску
О. Почтарь 1, А. Почтарь 2
Статистические методы контроля технологических процессов – эффективный инструмент отработки новых технологий, повышения качества выпускаемой продукции, снижения производственных затрат. В статье представлена оценка стабильности и воспроизводимости технологического процесса двухслойной металлизации кристаллов ниобата лития через свободную металлическую маску методом электронно-лучевого испарения, предложены рекомендации по отладке данного технологического процесса. Оценка выполнена на основе статистической обработки экспериментальных данных (ширина зазора и величина подпыления) с применением контрольных карт Шухарта.
Статистика в технологии
Как известно, контроль качества выпускаемой продукции может осуществляться различными методами, которые можно разделить на две большие группы: контроль качества изготавливаемых изделий и контроль технологических процессов.
При контроле качества изготавливаемых изделий вероятны дефекты изготовления, которые могут приводить к ранним отказам и необходимости переделки изделия. Однако при увеличении объемов производства такой подход перестает быть эффективным, и целесообразным становится переход к контролю технологических процессов. Такой метод контроля продукции позволяет минимизировать вероятность появления дефектов изготовления.
Инструментом для контроля технологических процессов являются статистические методы, в частности применение контрольных карт Шухарта, которые позволяют оценить стабильность и воспроизводимость технологических процессов, обнаружить и устранить особые причины изменения свойств продукции, приводящие к выпуску изделия, несоответствующего конструкторским требованиям.
Отработка любой технологии включает в себя проведение серии экспериментов по подбору режимов и необходимой оснастки. При разработке технологии металлизации кристаллов подбираются режимы металлизации и соответствующая оснастка – свободная металлическая маска. Режимы металлизации влияют на толщину и молекулярный состав покрытия, а также на зародышеобразование (поведение осаждаемых атомов на поверхности кристалла) [1]. Конфигурация оснастки определяет топологию металлизации.
Оснастка разрабатывается индивидуально под каждый тип кристалла и должна обеспечивать изготовление продукции, соответствующей требованиям конструкторской документации. Также необходимо, чтобы оснастка, при прочих равных условиях технологического процесса, обеспечивала стабильный и воспроизводимый процесс металлизации кристаллов. Для проверки оснастки на соответствие указанным требованиям возможны два подхода:
последовательное изготовление оснастки с различными геометрическими размерами и экспериментальное определение соответствия топологии металлизации конструкторской документации при использовании каждого типа оснастки;
применение статистических методов анализа процесса (в частности, использование контрольных карт Шухарта).
Последовательное изготовление оснастки нецелесообразно ввиду неэффективного расхода рабочего времени специалистов.
Применение статистических методов, а именно контрольных карт Шухарта, позволяет оценить стабильность процесса, а также обнаружить и устранить особые причины изменения свойств продукции, приводящие к получению изделия, несоответствующего конструкторским требованиям [2]. Это дает возможность не только подобрать оптимальные размеры оснастки, но и выявить зависимость влияния геометрических размеров оснастки на размеры элементов топологии металлизации. Определение такой зависимости обеспечит возможность перехода от контроля технологического процесса к его управлению, то есть позволит в дальнейшем изготавливать кристаллы с иной топологией металлизации без проведения длительных экспериментов.
Описание исследуемого технологического процесса
Металлизированный кристалл (рис. 1) представляет собой кристалл ниобата лития (LiNbO3) цилиндрической формы с нанесенными на одну торцевую сторону субмикронными слоями ванадия и никеля через свободную металлическую маску.
Свободная маска необходима для формирования металлических контактов заданных размеров, разделенных зазором, который обеспечивается перемычкой маски. Процесс металлизации выполнялся методом электронно-лучевого испарения (рис. 2).
Ширина зазора между металлическими контактами является основной контролируемой величиной, которая обеспечивается размерами элементов в маске. Ширина зазора измеряется по всей длине зазора. Так как от ширины зазора зависит работоспособность изделия, необходимо оценить стабильность технологического процесса металлизации кристаллов, а также выявить факторы, влияющие на эту величину.
Результаты статистического контроля
Параметром технологического процесса металлизации, влияние которого на стабильность процесса оценивается в работе, является ширина перемычки свободной металлической маски.
В работе использовалось два вида масок с одинаковой толщиной, но с разным значением ширины перемычки. В процессе отработки технологии было обнаружено, что на готовых кристаллах наблюдается уменьшение ширины зазора относительно ширины перемычки, что объясняется процессом подпыления материалов под свободную металлическую маску.
Таким образом, влияние ширины перемычки маски на стабильность и воспроизводимость технологического процесса металлизации кристаллов первично оценивается через оценку величины подпыления. Для этого были построены контрольные карты Шухарта для величины подпыления на одну сторону зазора (рис. 3), где UNPLx – верхняя контрольная граница; CLx – центральная линия (математическое ожидание); LNPLx – нижняя контрольная граница.
Каждый результат представляет собой один отдельный запуск процесса металлизации, поэтому краткосрочная вариация величин подпыления для последовательности индивидуальных значений была представлена как вариация от одной точки к другой [3].
Модуль скользящего размаха – это модуль разницы между двумя последующими значениями измеряемого параметра. Он отражает скорость изменения процесса, то есть является производной функцией процесса. Для построения контрольных карт Шухарта не нужно подтверждение нормальности распределения, так как метод контрольных карт не чувствителен к этому [3]. Но без точного определения типа распределения данных нельзя однозначно определить закон, по которому вычисляется производная функции процесса, поэтому для корректного анализа стабильности процесса необходимо отдельно производить оценку модуля скользящего размаха (рис. 4).
На основании экспериментальных данных, отраженных на контрольных картах Шухарта (см. рис. 3 и 4), были рассчитаны значения средней величины подпыления на одну сторону зазора ψ и средний скользящий размах: ψ = 7,214 мкм; mR = 1,038 мкм.
Были рассчитаны величины пределов контрольных карт Шухарта [2]:
для величины подпыления на одну сторону зазора: UNPLx = 9,976 мкм; CLx = 7,214 мкм; LNPLx = 4,452 мкм;
для скользящего размаха величины подпыления на одну сторону зазора: UCLr = 3,394 мкм; CLr = 1,038 мкм; LCLr = 0,000 мкм (поскольку отрицательного значения физически быть не может).
Расстояние между UNPLx и CLx для величины подпыления на одну сторону зазора и расстояние между UCLr и CLr для скользящего размаха принимаются за величину 3σ (σ – среднее квадратическое отклонение) [4]. Исходя из этого были рассчитаны величины 2σ и σ. Данные величины отражены на контрольных картах Шухарта.
При анализе контрольных карт Шухарта индивидуальных значений и скользящего размаха было выявлено, что 13 из 14 точек (92,8%) на графике подпыления на одну сторону зазора лежат в пределах ±σ, а также 14 из 14 точек (100%) лежат в пределах ±2σ, а соответственно и в пределах ±3σ.
Полученные данные укладываются в эмпирические правила определения стабильности процесса [3]. Процесс остается стабильным при замене свободной металлической маски: значения величин зазоров не выходят за пределы статистических границ ±3σ. Данный факт подтверждает, что ширина перемычки не оказывает влияния на стабильность процесса и величину подпыления на одну сторону зазора. Однако ширина перемычки маски влияет на фактическую величину получаемого зазора.
Для оценки фактической получаемой ширины зазора были построены контрольные карты Шухарта для индивидуальных значений ширины зазора (рис. 5) и рассчитаны статистические пределы по двум различным маскам:
для маски с шириной зазора 251 мкм: UNPLx = 243,08 мкм; CLx = 236,62 мкм; LNPLx = 230,16 мкм;
для маски с шириной зазора 185 мкм: UNPLx = 175,28 мкм; CLx = 170,50 мкм; LNPLx = 165,71 мкм.
Ширина зазора, установленная конструкторской документацией, составляет: b = 200–240 мкм.
По данным рис. 5 видно, что при металлизации через маску № 2 контрольные пределы и полученные результаты измерений находятся ниже значения, установленного конструкторской документацией, что говорит о невозможности использовать такую маску для металлизации кристаллов. При металлизации через маску № 1 наблюдается превышение значением UNPLx = 243,08 мкм максимального значения ширины зазора, установленного конструкторской документацией. Такой результат говорит о том, что имеется вероятность несоответствия изделия техническим требованиям при сохранении стабильности технологического процесса.
Для анализа воспроизводимости процессов металлизации кристаллов были рассчитаны индексы воспроизводимости Ср для данных по двум маскам: Ср(маска 1) = 3,09 и Ср(маска 2) = 4,11. Так как процессы, согласно [4], имеют естественный допуск ±3σ, и величина Ср для обоих случаев больше 1,33 (коэффициент воспроизводимости процесса), то можно говорить о стабильности и воспроизводимости процесса. Однако, очевидно, что процесс нуждается в центровке (смещение величины математического ожидания к центру).
Проведенный статистический анализ позволяет рассчитать требуемые геометрические размеры маски для обеспечения заданной конструкторской документацией ширины зазора, без экспериментального подбора конфигурации оснастки.
Для смещения статистических границ процесса металлизации кристаллов внутрь конструкторских границ необходимо смещение математического ожидания (CLx) до 220 мкм (см. рис. 5). При возможном отклонении ширины зазора в диапазоне ±3σ процесс, находящийся в стабильном состоянии, будет обеспечивать сохранение ширины зазора согласно допуску по конструкторской документации b = 220 ± 20 мкм.
Установлено, что величина подпыления на одну сторону зазора составляет 7,2 мкм (математическое ожидание для величины подпыления на одну сторону зазора). Исходя из полученных данных с помощью соотношения (1) рассчитывается ширина перемычки металлической маски, которая обеспечит зазор 220 мкм:
w = b + (2 · ψ). (1)
Подставим численные значения в формулу (1):
w = 220 + (2 · 7,2) = 234,4 мкм.
Для удобства изготовления металлической маски полученное значение ширины перемычки было округлено до w = 240 мкм.
Заключение
В работе методами статистической обработки результатов технологического процесса была доказана стабильность и воспроизводимость процесса металлизации кристаллов.
На основе полученных данных вынесены рекомендации по отладке технологического процесса: для получения стабильного процесса металлизации кристаллов с шириной получаемого зазора, укладывающегося в допуски, установленные конструкторской документацией b = 220 ± 20 мкм, необходимо использовать свободную металлическую маску толщиной 0,12 мм и шириной перемычки w = 240 мкм.
Для дальнейшего повышения однородности процесса и уменьшения величины среднего квадратического отклонения необходимо рассмотреть факторы, влияющие на величину подпыления. Такими факторами являются: температура кристалла в процессе металлизации и температура испаряемых материалов в тигле, которую можно варьировать путем изменения тока луча электронно-лучевой пушки. Уменьшив разброс значений величины подпыления, появится возможность изготовления кристаллов с шириной зазора, близкой к величине математического ожидания.
В результате проведенной работы получена замкнутая модель производственного процесса металлизации кристалла из LiNbO3 слоями ванадий-никель через свободную металлическую маску методом электронно-лучевого испарения с формированием заданной топологии.
ЛИТЕРАТУРА
Иванов А., Смирнов Б. Электронно-лучевое напыление: технология и оборудование // НАНОИНДУСТРИЯ. 2012. № 6 (36). С. 28–34.
Локтев Д. А. Статистическое управление производственными процессами – ключ к успеху современного промышленного предприятия // Известия МГТУ «МАМИ». 2014. № 1(19), т. 2. С. 128–136.
Уилер Д., Чамберс Д. Статистическое управление процессами. Оптимизация бизнеса с использованием контрольных карт Шухарта / Пер. с англ., 2‑е изд. М.: Альпина Паблишер, 2016. 582 с.
ГОСТ Р ИСО 7870-2-2015. Статистические методы. Контрольные карты Шухарта. Ч. 2. М., 2015. 46 с.
технологического процесса металлизации через свободную маску
О. Почтарь 1, А. Почтарь 2
Статистические методы контроля технологических процессов – эффективный инструмент отработки новых технологий, повышения качества выпускаемой продукции, снижения производственных затрат. В статье представлена оценка стабильности и воспроизводимости технологического процесса двухслойной металлизации кристаллов ниобата лития через свободную металлическую маску методом электронно-лучевого испарения, предложены рекомендации по отладке данного технологического процесса. Оценка выполнена на основе статистической обработки экспериментальных данных (ширина зазора и величина подпыления) с применением контрольных карт Шухарта.
Статистика в технологии
Как известно, контроль качества выпускаемой продукции может осуществляться различными методами, которые можно разделить на две большие группы: контроль качества изготавливаемых изделий и контроль технологических процессов.
При контроле качества изготавливаемых изделий вероятны дефекты изготовления, которые могут приводить к ранним отказам и необходимости переделки изделия. Однако при увеличении объемов производства такой подход перестает быть эффективным, и целесообразным становится переход к контролю технологических процессов. Такой метод контроля продукции позволяет минимизировать вероятность появления дефектов изготовления.
Инструментом для контроля технологических процессов являются статистические методы, в частности применение контрольных карт Шухарта, которые позволяют оценить стабильность и воспроизводимость технологических процессов, обнаружить и устранить особые причины изменения свойств продукции, приводящие к выпуску изделия, несоответствующего конструкторским требованиям.
Отработка любой технологии включает в себя проведение серии экспериментов по подбору режимов и необходимой оснастки. При разработке технологии металлизации кристаллов подбираются режимы металлизации и соответствующая оснастка – свободная металлическая маска. Режимы металлизации влияют на толщину и молекулярный состав покрытия, а также на зародышеобразование (поведение осаждаемых атомов на поверхности кристалла) [1]. Конфигурация оснастки определяет топологию металлизации.
Оснастка разрабатывается индивидуально под каждый тип кристалла и должна обеспечивать изготовление продукции, соответствующей требованиям конструкторской документации. Также необходимо, чтобы оснастка, при прочих равных условиях технологического процесса, обеспечивала стабильный и воспроизводимый процесс металлизации кристаллов. Для проверки оснастки на соответствие указанным требованиям возможны два подхода:
последовательное изготовление оснастки с различными геометрическими размерами и экспериментальное определение соответствия топологии металлизации конструкторской документации при использовании каждого типа оснастки;
применение статистических методов анализа процесса (в частности, использование контрольных карт Шухарта).
Последовательное изготовление оснастки нецелесообразно ввиду неэффективного расхода рабочего времени специалистов.
Применение статистических методов, а именно контрольных карт Шухарта, позволяет оценить стабильность процесса, а также обнаружить и устранить особые причины изменения свойств продукции, приводящие к получению изделия, несоответствующего конструкторским требованиям [2]. Это дает возможность не только подобрать оптимальные размеры оснастки, но и выявить зависимость влияния геометрических размеров оснастки на размеры элементов топологии металлизации. Определение такой зависимости обеспечит возможность перехода от контроля технологического процесса к его управлению, то есть позволит в дальнейшем изготавливать кристаллы с иной топологией металлизации без проведения длительных экспериментов.
Описание исследуемого технологического процесса
Металлизированный кристалл (рис. 1) представляет собой кристалл ниобата лития (LiNbO3) цилиндрической формы с нанесенными на одну торцевую сторону субмикронными слоями ванадия и никеля через свободную металлическую маску.
Свободная маска необходима для формирования металлических контактов заданных размеров, разделенных зазором, который обеспечивается перемычкой маски. Процесс металлизации выполнялся методом электронно-лучевого испарения (рис. 2).
Ширина зазора между металлическими контактами является основной контролируемой величиной, которая обеспечивается размерами элементов в маске. Ширина зазора измеряется по всей длине зазора. Так как от ширины зазора зависит работоспособность изделия, необходимо оценить стабильность технологического процесса металлизации кристаллов, а также выявить факторы, влияющие на эту величину.
Результаты статистического контроля
Параметром технологического процесса металлизации, влияние которого на стабильность процесса оценивается в работе, является ширина перемычки свободной металлической маски.
В работе использовалось два вида масок с одинаковой толщиной, но с разным значением ширины перемычки. В процессе отработки технологии было обнаружено, что на готовых кристаллах наблюдается уменьшение ширины зазора относительно ширины перемычки, что объясняется процессом подпыления материалов под свободную металлическую маску.
Таким образом, влияние ширины перемычки маски на стабильность и воспроизводимость технологического процесса металлизации кристаллов первично оценивается через оценку величины подпыления. Для этого были построены контрольные карты Шухарта для величины подпыления на одну сторону зазора (рис. 3), где UNPLx – верхняя контрольная граница; CLx – центральная линия (математическое ожидание); LNPLx – нижняя контрольная граница.
Каждый результат представляет собой один отдельный запуск процесса металлизации, поэтому краткосрочная вариация величин подпыления для последовательности индивидуальных значений была представлена как вариация от одной точки к другой [3].
Модуль скользящего размаха – это модуль разницы между двумя последующими значениями измеряемого параметра. Он отражает скорость изменения процесса, то есть является производной функцией процесса. Для построения контрольных карт Шухарта не нужно подтверждение нормальности распределения, так как метод контрольных карт не чувствителен к этому [3]. Но без точного определения типа распределения данных нельзя однозначно определить закон, по которому вычисляется производная функции процесса, поэтому для корректного анализа стабильности процесса необходимо отдельно производить оценку модуля скользящего размаха (рис. 4).
На основании экспериментальных данных, отраженных на контрольных картах Шухарта (см. рис. 3 и 4), были рассчитаны значения средней величины подпыления на одну сторону зазора ψ и средний скользящий размах: ψ = 7,214 мкм; mR = 1,038 мкм.
Были рассчитаны величины пределов контрольных карт Шухарта [2]:
для величины подпыления на одну сторону зазора: UNPLx = 9,976 мкм; CLx = 7,214 мкм; LNPLx = 4,452 мкм;
для скользящего размаха величины подпыления на одну сторону зазора: UCLr = 3,394 мкм; CLr = 1,038 мкм; LCLr = 0,000 мкм (поскольку отрицательного значения физически быть не может).
Расстояние между UNPLx и CLx для величины подпыления на одну сторону зазора и расстояние между UCLr и CLr для скользящего размаха принимаются за величину 3σ (σ – среднее квадратическое отклонение) [4]. Исходя из этого были рассчитаны величины 2σ и σ. Данные величины отражены на контрольных картах Шухарта.
При анализе контрольных карт Шухарта индивидуальных значений и скользящего размаха было выявлено, что 13 из 14 точек (92,8%) на графике подпыления на одну сторону зазора лежат в пределах ±σ, а также 14 из 14 точек (100%) лежат в пределах ±2σ, а соответственно и в пределах ±3σ.
Полученные данные укладываются в эмпирические правила определения стабильности процесса [3]. Процесс остается стабильным при замене свободной металлической маски: значения величин зазоров не выходят за пределы статистических границ ±3σ. Данный факт подтверждает, что ширина перемычки не оказывает влияния на стабильность процесса и величину подпыления на одну сторону зазора. Однако ширина перемычки маски влияет на фактическую величину получаемого зазора.
Для оценки фактической получаемой ширины зазора были построены контрольные карты Шухарта для индивидуальных значений ширины зазора (рис. 5) и рассчитаны статистические пределы по двум различным маскам:
для маски с шириной зазора 251 мкм: UNPLx = 243,08 мкм; CLx = 236,62 мкм; LNPLx = 230,16 мкм;
для маски с шириной зазора 185 мкм: UNPLx = 175,28 мкм; CLx = 170,50 мкм; LNPLx = 165,71 мкм.
Ширина зазора, установленная конструкторской документацией, составляет: b = 200–240 мкм.
По данным рис. 5 видно, что при металлизации через маску № 2 контрольные пределы и полученные результаты измерений находятся ниже значения, установленного конструкторской документацией, что говорит о невозможности использовать такую маску для металлизации кристаллов. При металлизации через маску № 1 наблюдается превышение значением UNPLx = 243,08 мкм максимального значения ширины зазора, установленного конструкторской документацией. Такой результат говорит о том, что имеется вероятность несоответствия изделия техническим требованиям при сохранении стабильности технологического процесса.
Для анализа воспроизводимости процессов металлизации кристаллов были рассчитаны индексы воспроизводимости Ср для данных по двум маскам: Ср(маска 1) = 3,09 и Ср(маска 2) = 4,11. Так как процессы, согласно [4], имеют естественный допуск ±3σ, и величина Ср для обоих случаев больше 1,33 (коэффициент воспроизводимости процесса), то можно говорить о стабильности и воспроизводимости процесса. Однако, очевидно, что процесс нуждается в центровке (смещение величины математического ожидания к центру).
Проведенный статистический анализ позволяет рассчитать требуемые геометрические размеры маски для обеспечения заданной конструкторской документацией ширины зазора, без экспериментального подбора конфигурации оснастки.
Для смещения статистических границ процесса металлизации кристаллов внутрь конструкторских границ необходимо смещение математического ожидания (CLx) до 220 мкм (см. рис. 5). При возможном отклонении ширины зазора в диапазоне ±3σ процесс, находящийся в стабильном состоянии, будет обеспечивать сохранение ширины зазора согласно допуску по конструкторской документации b = 220 ± 20 мкм.
Установлено, что величина подпыления на одну сторону зазора составляет 7,2 мкм (математическое ожидание для величины подпыления на одну сторону зазора). Исходя из полученных данных с помощью соотношения (1) рассчитывается ширина перемычки металлической маски, которая обеспечит зазор 220 мкм:
w = b + (2 · ψ). (1)
Подставим численные значения в формулу (1):
w = 220 + (2 · 7,2) = 234,4 мкм.
Для удобства изготовления металлической маски полученное значение ширины перемычки было округлено до w = 240 мкм.
Заключение
В работе методами статистической обработки результатов технологического процесса была доказана стабильность и воспроизводимость процесса металлизации кристаллов.
На основе полученных данных вынесены рекомендации по отладке технологического процесса: для получения стабильного процесса металлизации кристаллов с шириной получаемого зазора, укладывающегося в допуски, установленные конструкторской документацией b = 220 ± 20 мкм, необходимо использовать свободную металлическую маску толщиной 0,12 мм и шириной перемычки w = 240 мкм.
Для дальнейшего повышения однородности процесса и уменьшения величины среднего квадратического отклонения необходимо рассмотреть факторы, влияющие на величину подпыления. Такими факторами являются: температура кристалла в процессе металлизации и температура испаряемых материалов в тигле, которую можно варьировать путем изменения тока луча электронно-лучевой пушки. Уменьшив разброс значений величины подпыления, появится возможность изготовления кристаллов с шириной зазора, близкой к величине математического ожидания.
В результате проведенной работы получена замкнутая модель производственного процесса металлизации кристалла из LiNbO3 слоями ванадий-никель через свободную металлическую маску методом электронно-лучевого испарения с формированием заданной топологии.
ЛИТЕРАТУРА
Иванов А., Смирнов Б. Электронно-лучевое напыление: технология и оборудование // НАНОИНДУСТРИЯ. 2012. № 6 (36). С. 28–34.
Локтев Д. А. Статистическое управление производственными процессами – ключ к успеху современного промышленного предприятия // Известия МГТУ «МАМИ». 2014. № 1(19), т. 2. С. 128–136.
Уилер Д., Чамберс Д. Статистическое управление процессами. Оптимизация бизнеса с использованием контрольных карт Шухарта / Пер. с англ., 2‑е изд. М.: Альпина Паблишер, 2016. 582 с.
ГОСТ Р ИСО 7870-2-2015. Статистические методы. Контрольные карты Шухарта. Ч. 2. М., 2015. 46 с.
Отзывы читателей
