Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

THE MATHEMATICAL MODEL OF STRESSES INITIATED BY LASER PULSE

Isaeva A.S. 1 Ryndin E.A. 2, 1 Ryzhuk R.V. 3
1 Southern Federal University
2 Southern Scientific Center of Russian Academy of Sciences
3 National Research Nuclear University MEPhI
Triboluminescent materials application is an important area of nondestructive testing methods development. The mathematical model that sets a correlation between stresses on the triboluminescent material and parameters of initiating laser pulse for the purpose of nondestructive testing is proposed. The heat equation and thermoelasticity equation are used in this model. Both equations are solved numerically with using finite differences method. The thermoelasticity equation was solved in two dimensions, heat equation was solved in three dimensions and then converted into a two dimension temperature distribution. Nonstationary of temperature and stress tensor components for thin film triboluminescent material (ZnS), irradiated by laser pulse are obtained. Solutions are got both for uniform and Gauss distributions of power density of laser ray. The mathematical model could be useful in field of defect monitoring system development.
model
method
triboluminescence
nondestructive testing
1. Banishev A.F., Bolshukhin V.A., Azarov A.D. Fizika i khimiya obrabotki materialov. 2008. no 3, pp. 24–27.
2. Glushkov A.A., Prostakishin A.S., Slyusareva E.A., Sizykh A.G. Journal of Siberian Federal University. Mathematics & Physics. 2012. no 5(1), pp. 57–62.
3. Kovalenko A.D. Osnovy termouprugosti [The fundamentals of termoelasticity]. Kiev, Naukova Dumka, 1970, pp. 43–48.
4. Korshunov V.V., Senchukov F.D., Shmurak S.Z. Pisma v ZHETF (JETP Letters) 1971. t. 13, pp. 408–412.
5. Makarova N.Y. Tezisy Vserossiyskoy nauchno-prakticheskoy konferencii Aktualnye problemy aviacii i kosmonavtiki (The actual problems of aviation and astronautics) Krasnoyarsk 2005, pp. 65–66.
6. Makarova N.Y., Sparzhakin A.G., Kornilov P.P., Klimenko Y.S., Skornyakov R.A. Tezisy Vserossiyskoy nauchno-prakticheskoy konferencii Aktualnye problemy aviacii i kosmonavtiki (The actual problems of aviation and astronautics) Krasnoyarsk 2005, pp. 67–68.
7. Morozova N.K., Karetnikov I.A., Blinov V.V., Gavrischuk E.M. Fizika I tekhnika poluprovodnikov. 2001. Т. 35, no. 1, pp. 25–33.
8. Tatmyshevsky K.V. Mikrosistemnaya tekhnika. 2004. no 12, pp. 4–10.
9. Tatmyshevsky K.V. Nauchnye osnovy rascheta I proektirovaniya mehanolyuminescentnyh chuvstvitelnykh elementov datchikov impulsnogo davleniya. Vladimir, 2009, pp. 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.