Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ МЕХАНИЧЕСКИХ НАПРЯЖЕНИЙ, ИНИЦИИРОВАННЫХ ЛАЗЕРНЫМ ИМПУЛЬСОМ

Исаева А.С. 1 Рындин, Е.А. 2, 1 Рыжук Р.В. 3
1 ФГАОУ ВПО «Южный федеральный университет»
2 ФГБУН «Южный научный центр Российской академии наук»
3 3ФГАОУ ВПО «Национальный исследовательский ядерный университет МИФИ»
Применение механолюминесцентных материалов является важным направлением разработки методов неразрушающего контроля. Предложена математическая модель, связывающая величины механических напряжений в механолюминесцентном материале с параметрами инициирующего лазерного импульса для целей неразрушающего контроля. В математической модели используются уравнения теплопроводности и термоупругости, которые решаются численно с применением метода конечных разностей. Получены нестационарные распределения по координатам температур и компонент тензора напряжений в тонкой пленке механолюминофора ZnS, при длительности лазерных импульсов 1 мс для равномерного и гауссовского распределений плотности мощности по площади сечения лазерного луча. Поскольку экспериментально механолюминесценция, как правило, генерируется лазерными импульсами, разработанная модель позволяет достичь более эффективного применения экспериментальных данных для разработки систем мониторинга дефектов на основе явления механолюминесценции.
модель
метод
механолюминесценция
неразрушающий контроль
1. Банишев А.Ф., Большухин В.А., Азаров А.Д. Механолюминесценция мелкодисперсного порошкообразного твердого раствора SrAl2O4:(Eu2 + ,Dy3 +), возбуждаемая воздействием лазерных импульсов // Физика и химия обработки материалов. – 2008. – № 3. – С. 24–27.
2. Глушков А.А., Простакишин А.С., Слюсарева Е.А., Сизых А.Г. Влияние пространственного распределения интенсивности лазерного излучения на кинетику фотообесцвечивания красителя в полимерной матрице // Journal of Siberian Federal University. Mathematics & Physics. – 2012. – № 5(1). – С. 57–62.
3. Коваленко А.Д. Основы термоупругости. – Киев: Наукова Думка, 1970. – C. 43–48.
4. Коршунов В.В., Сенчуков Ф.Д., Шмурак С.З. Исследование временных характеристик деформационной люминесценции // Письма в ЖЭТФ. – 1971. – Т. 13. – С. 408–412.
5. Макарова Н.Ю. Интеллектуальные композиционные панели со встроенными механолюминесцентными датчиками давления // Актуальные проблемы авиации и космонавтики: тезисы Всеросийской научн.-практ, конференции. – Красноярск, 2005. – С. 65–66.
6. Система мониторинга напряженно-деформированного состояния крыльев самолета / Н.Ю. Макарова, А.Г. Спажакин, П.П. Корнилов, Ю.С. Клименко, Р.А. Скорняков // Актуальные проблемы авиации и космонавтики: тезисы Всероссийской научн.-практ. конференции. – Красноярск, 2005. – С. 67–68.
7. Исследование центров люминесценции, обязанных присутствию меди и цинка в ZnSe / Н.К. Морозова, И.А. Каретников, В.В. Блинов, Е.М. Гаврищук // Физика и техника полупроводников. – 2001. – Т. 35, Вып. 1. – С. 25–33.
8. Татмышевский К.В. Механолюминесцентные (светогенерационные) сенсорные элементы для современных информационно-измерительных технологий // Микросистемная техника. – 2004. – № 12. – С. 4–10.
9. Татмышевский К.В. Научные основы расчета и проектирования механолюминесцентных чувствительных элементов датчиков импульсного давления: дис. ... д-ра тех. наук. – Владимир, 2009. – С. 46.

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

Многообещающими в этом отношении являются механолюминесцентные материалы. Анализ литературы по данной тематике показывает недостаточное число прикладных исследований механолюминесцентных материалов с целью разработки методов НК, в то время как в работах
[5, 6, 8, 9] обозначается перспективность данного направления.

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

Математическая модель
механических напряжений

Приведенная в работе [8] математическая модель механолюминесцентного внутрицентрового излучения цинксульфидных порошковых люминофоров не всегда может быть применена при разработке систем мониторинга, поскольку содержит большое количество несправочных величин и связывает поток излучения при механолюминесценции только с деформационными напряжениями в материале. В то же время экспериментально механолюминесценция, как правило, генерируется путем воздействия на материал лазерного импульса малой длительности (10–6–10–3 с) и радиуса [1, 4, 7] и, следовательно, величина напряжения в материале не следует явно из экспериментальных данных. Таким образом, разработка математической модели, связывающей величины механических напряжений в материале с параметрами инициирующего лазерного импульса, является актуальной задачей.

Для решения этой задачи в данной работе предложена модель, включающая нестационарное уравнение теплопроводности, описывающее распределение температуры в объекте исследования под воздействием лазерного импульса, и уравнения термоупругости [3], позволяющие определить механические напряжения, возникающие вследствие теплового воздействия:

Eqn21.wmf

 (1) 

 (2)

 (3)

 (4)

где ρ – плотность материала; c – удельная теплоемкость; λq – коэффициент теплопроводности; T – абсолютная температура; t – время; x, y, z – координаты; f(x, y, z) – функция плотности мощности источников тепла; E – модуль Юнга; αТ – коэффициент линейного теплового расширения; σxx, σxy, σyy – компоненты тензора напряжения;
Т0 – начальная температура тела.

Уравнения (1)–(4) решались для граничных условий:

Eqn22.wmf (5)

Eqn23.wmf (6)

Eqn24.wmf (7)

Eqn25.wmf (8)

где D – диаметр лазерного пучка в месте его падения на пластину; ЕL – энергия лазерного импульса; τ – длительность лазерного импульса.

Лазерный импульс воздействует на малую площадь поверхности тела S = πD2/4 (рис. 1).

рис_35.wmf

Рис. 1. Область воздействия лазерного импульса

В случае равномерного распределения плотности мощности по площади сечения лазерного луча задаются граничные условия второго рода (7).

Распределение интенсивности лазерного импульса в пучке, как правило, соответствует нормальному закону распределения Гаусса [2]. Для лучшей корреляции с экспериментальными данными следует задать для области воздействия лазерного импульса следующие граничные условия второго рода:

Eqn26.wmf (9)

В общем случае задачи тепропроводности и термоупругости решаются для трех измерений, но с учетом того, что объектом исследования является тонкая пленка или пластина, целесообразно решать задачу термоупругости для двух измерений (x, y), а уравнение теплопроводности для трех измерений с последующим преобразованием в двумерное распределение температуры в соответствии с выражением:

Eqn27.wmf (10)

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

Система (1)–(10) решалась численно с использованием метода конечных разностей.

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

• для изучения процесса механолюминесценции существенное значение имеют напряженные состояния на малом отрезке времени, недостаточном для релаксации материала, и, следовательно, потери тепла за счет конвективного теплообмена не будут заметно влиять на погрешность решения;

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

По результатам проведенных вычислений непосредственное решение системы уравнений (1)–(10) приводит к неадекватным результатам вне зависимости от метода дискретизации уравнений первого порядка (3) и (4). С целью решения данной проблемы уравнения (3) и (4) системы (1)–(10) были приведены посредством дополнительного дифференцирования к виду:

Eqn28.wmf (11)

Eqn29.wmf (12)

Численное решение системы уравнений (1), (2), (5)–(12) выполнялось для исходных данных, соответствующих сульфиду цинка, обладающему механолюминесцентными свойствами.

Результаты моделирования

рис_36.tif

Рис. 2. Распределение по координатам температуры и напряжений на поверхности пластины
для равномерного распределения плотности мощности по площади сечения лазерного луча:
а – распределение температуры в момент прекращения действия лазерного импульса;
б, г, е – распределение напряжений σxx, σxy, σxz в пластине в момент прекращения действия лазерного импульса соответственно; в, д, ж – распределение напряжений в σxx, σxy, σxz пластине спустя 1 мс после прекращения действия лазерного импульса соответственно

рис_37.tif

Рис. 3. Распределение по координатам температуры и напряжений на поверхности пластины для гауссовского распределения плотности мощности по площади сечения лазерного луча:
а – распределение температуры в момент прекращения действия лазерного импульса;
б, г, е – распределение напряжений σxx, σxy, σxz в пластине в момент прекращения действия лазерного импульса соответственно; в, д, ж – распределение напряжений в σxx, σxy, σxz пластине спустя 1 мс после прекращения действия лазерного импульса соответственно

На рис. 2–3 приведены распределения по координатам температур и компонент тензора напряжений σxx, σxy, σyy в различные моменты времени, полученные из решения системы (1), (2), (5)–(12) при длительности лазерных импульсов 1 мс для равномерного и гауссовского распределений плотности мощности по площади сечения лазерного луча.

Заключение

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

Работа выполнена при финансовой поддержке Министерства образования и науки Российской Федерации (Гос. контракт № 16.740.11.0425 от 03.12.2010, гос. соглашение № 14.А18.21.0126) в рамках ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009–2013 годы.

Рецензенты:

Агеев О.А., д.т.н., профессор, директор НОЦ «Нанотехнологии»;

Жорник А.И., д.ф.-м.н., профессор кафедры теоретической, общей физики и технологии ФГБОУ ВПО ТГПИ.

Работа поступила в редакцию 26.10.2012.


Библиографическая ссылка

Исаева А.С., Рындин, Е.А., Рындин, Е.А., Рыжук Р.В. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ МЕХАНИЧЕСКИХ НАПРЯЖЕНИЙ, ИНИЦИИРОВАННЫХ ЛАЗЕРНЫМ ИМПУЛЬСОМ // Фундаментальные исследования. – 2012. – № 11-3. – С. 609-614;
URL: http://www.fundamental-research.ru/ru/article/view?id=30582 (дата обращения: 17.01.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074