Full Text

Введение При исследовании задач, связанных с прогнозированием и повышением надежности и безопасности технических систем, имеющих покрытия, важную роль играют методики диагностики отслоений в покрытиях и построение адекватных моделей, описывающих процесс отслоения. Это связано с практической значимостью идентификации такого рода дефектов, возникающих как на макро-, так и на микроуровне. При этом особо стоит отметить функционально-градиентные [1] покрытия, свойства которых зависят от координат, как наиболее динамично развивающийся вид покрытий. Одной из причин возникновения отслоений, а впоследствии - трещин и полостей между подложкой и покрытием, является наличие полей предварительных напряжений [2, 3] и остаточных деформаций, в связи чем необходим учет данных факторов при моделировании. В современной литературе моделированию неоднородных покрытий посвящено множество работ. Так, в работе [4] на основе метода акустической эмиссии разработана методика оценки поведения функционального покрытия, нанесенного на материал. В частности, методика применялась для керамических покрытий, для которых характерно наличие механизмов накопления повреждений, постепенно приводящих к нарушению работоспособности конструкции. Авторами работы [5] рассмотрена модель пьезоэлектрического покрытия полупространства, исследованы динамические характеристики системы и особенности распространения поверхностных волн. В статье [6] представлен анализ распределения диссипации энергии, передаваемой на покрытие дорожного полотна при проезде автотранспорта, для чего была разработана специальная динамическая модель напряженно-деформированного состояния многослойного полупространства. Проведены эксперименты для подложек различной степени жесткости и капитальности. В работе [7] построена итерационная схема идентификации неизвестных функций, характеризующих параметры Ламе и плотность толстого покрытия, нанесенного на цилиндрический волновод, на основе анализа акустического отклика. Особый интерес представляют работы, в которых рассмотрены аспекты моделирования отслоений в покрытиях. Среди отечественных авторов стоит отметить Р.В. Гольдштейна, К.Б. Устинова. Например, в работе [8] отслаивающееся покрытие моделируется как пластина с упругим защемлением к подложке. Проведен анализ параметров потери устойчивости, определения критического напряжения в случаях изотропного и анизотропного покрытий и подложек. В результате выявлено, что потеря устойчивости для широких диапазонов свойств может описываться одним параметром, зависящим от упругих свойств покрытия и подложки, а также от отношения длины отслоения и толщины покрытия. В статье [9] проведено асимптотическое исследование смещения берегов трещины при нормальном отрыве покрытия от подложки, моделируемых двумя упругими полосами, как вблизи, так и вдали от ее вершины. В работах [10-13] представлены подходы к решению задач об отслоении покрытий, основанные на моделировании отслоившегося участка с помощью одного из вариантов теории пластин, граничные условия для которых ставятся исходя из рассмотрения задачи о контакте полубесконечного отслоения с основанием, решаемой аналитически. В рамках этого подхода представлено решение нескольких задач о полосе, контактирующей с полуплоскостью из другого материала. В статье [14] произведено конечно-элементное моделирование колебаний полосы с преднапряженным покрытием при наличии отслоения. Также стоит отметить статьи зарубежных исследователей, в частности Т. Хатчинсона и соавторов [15-18]. В работе [16] решена задача о потере устойчивости отслоившегося участка покрытия в случае жесткой подложки с прямолинейной поверхностью. Статья [19] посвящена исследованию отслоений, вызванных остаточными напряжениями. Метод исследования заключается в измерении ударной вязкости покрытий при последовательном растяжении и изгибе. На основе балочной теории удается достаточно точно описать процесс распространения отслоения. Также исследовано влияние остаточного напряжения на процесс отслоения. Выделим ряд работ, посвященных индентированию (вдавливанию жесткого штампа) неоднородных покрытий [20-23]. В статье [22] исследовано поведение линейно-упругого полупространства с функционально-градиентным покрытием при внедрении в него индентора (штампа) со сферической формой наконечника, построена схема расчета полей деформаций и контактных напряжений для произвольных законов изменения модуля Юнга для покрытия и подложки. В работе [23] представлено сравнение нескольких моделей деформирования неоднородной упругой полосы под воздействием индентора, учитывающее различный характер неоднородности свойств. Результаты сравниваются конечно-элементными решениями в однородном случае. Настоящее исследование посвящено разработке моделей отслоения предварительно напряженных неоднородных покрытий на примере задачи о колебаниях неоднородной полосы с покрытием. Основной целью работы является анализ влияния механических и геометрических характеристик покрытий, в частности уровня предварительного напряжения, на функции раскрытия в области отслоения, а также функции смещения и амплитудно-частотные характеристики на верхней границе полосы. Такого рода исследование необходимо для выявления наиболее информативных частотных диапазонов зондирования и построения подхода к решению обратных задач [24] идентификации параметров отслоения на основе анализа акустического отклика. 1. Постановка задачи Рассмотрим установившиеся колебания с частотой упругой изотропной неоднородной по толщине полосы, состоящей из основного слоя-подложки толщиной Н и предварительно напряженного (ПН) слоя-покрытия толщиной h, располагающегося в области в плоской постановке. Между основным слоем и слоем-покрытием имеется отслоение размерностью . Нижняя грань слоя S1 жестко защемлена, на части верхней границы S20 приложена нагрузка, определяемая вектором , где (рис. 1). В рамках принятых гипотез у вектора смещений ненулевыми являются компоненты . Рис. 1. Схематичное изображение рассматриваемой полосы с ПН покрытием Fig. 1. The schematic illustration of a coated strip Линеаризованная постановка задачи об установившихся колебаниях предварительно напряженно-деформированного анизотропного тела с плотностью после отделения временного множителя имеет вид [2, 25] (1.1) Здесь величины с верхним индексом «0» относятся к начальному состоянию, величины без индекса - компоненты соответствующих добавочных векторов и тензоров: - компоненты вектора перемещений, - компоненты симметричного тензора напряжений Кирхгофа, - компоненты несимметричного добавочного тензора напряжений Пиолы. Запишем уравнения установившихся колебаний полосы: (1.2) Для подложки без ПН определяющие соотношения для изотропного материала имеют вид: (1.3) В случае наличия одноосного ПН для покрытия определяющие соотношения записываются в форме: (1.4) Подобная модель ПНДС, учитывающая только предварительные напряжения, является довольно распространенной в литературе [26-28]. Запишем граничные условия для жесткого защемления нижней границы полосы и нагружения на верхней границе: (1.5) Также должны выполняться условия отсутствия напряжений в области отслоения: (1.6) В сформулированной задаче требуется определить функции смещения на верхней границе слоя с целью анализа влияния исходных параметров задачи, в частности, уровня предварительных напряжений на них. 2. Вспомогательные задачи Для решения исходной задачи рассмотрим две постановки вспомогательных задач для каждой из полос. Задача 1. Рассматриваются колебания основной полосы-подложки, для которой заданы функции смещения в основании и на границе полос: (2.1) Задача 2. Рассматриваются колебания полосы-покрытия, для которой заданы функции смещения на нижней грани и нагрузка верхней: (2.2) Для решения задач используется преобразование Фурье [29] по переменной : (2.3) Уравнения задачи (1.2), записанные в трансформантах, примут вид: (2.4) После исключения из первого уравнения запишем систему (2.4) в каноническом виде: (2.5) Граничные условия для основного слоя в трансформантах: (2.6) Аналогично можно записать итоговую систему для предварительно напряженного слоя-покрытия в каноническом виде: (2.7) Граничные условия для слоя-покрытия в трансформантах: (2.8) Проведем обезразмеривание задач (2.5), (2.6) и (2.7), (2.8) с помощью введения следующих обозначений: , , , , , , , , , , : (2.9) Граничные условия для основного слоя в безразмерных параметрах: (2.10) Аналогично можно записать уравнения для слоя-покрытия: (2.11) Граничные условия для слоя-покрытия в безразмерных параметрах: (2.12) В дальнейшем знак "^" в записи будем опускать. 3. Методика решения вспомогательных задач Решение каждой из вспомогательных задач будем искать в виде линейных комбинаций решений специальных подзадач. 1. Задача для основного слоя Рассмотрим две задачи, уравнения для которых будут иметь вид (2.9), а граничные условия в первом случае: (3.1) во втором случае: (3.2) Тогда решение первой вспомогательной задачи можно записать в виде (3.2) где - решение первой и второй подзадач соответственно. По найденным находятся соответствующие . 2. Задача для покрытия Для второй вспомогательной задачи ввиду ее граничных условий будем рассматривать три подзадачи, уравнения которых будут иметь вид (2.11), а граничные условия: - для первой: (3.3) - для второй: (3.4) - для третьей: (3.5) Тогда решение второй вспомогательной задачи можно записать в виде , (3.6) где - решение первой, второй и третьей подзадач соответственно. Аналогично определяются и . Численное решение вспомогательных задач осуществлялось с помощью метода пристрелки [30], показавшего достаточно высокую точность. 4. Построение операторных соотношений В результате решения вспомогательных задач в трансформантах определяются функции на границе слоев , и , , соответствующие решениям 1-й и 2-й вспомогательных задач, для которых составляется система 4 уравнений: (4.1) В результате решения системы получаются выражения для и через . По ним вычисляются соответствующие , , выраженные через . Их можно записать в виде Затем проводится обратное преобразование Фурье функций : С использованием свойства прямого и обратного преобразований Фурье, переводящих произведение функций в их свертку [31], из условий (1.6) в безразмерном виде (4.2) строятся интегральные операторные соотношения относительно оригиналов вида: (4.3) Отметим, что ядра полученных интегральных соотношений (4.3) являются гиперсинулярными [32, 33], для их нахождения требуется применение специальных методов [34, 35]. 5. Методика решения интегральных уравнений Для решения интегральных уравнений (4.3) использован метод коллокаций, при котором отрезок разбивается на N отрезков , на каждом из которых (5.1) Тогда для уравнений (4.3) можно записать следующую систему из уравнений, полученную с помощью дискретизации: (5.2) При этом для вычисления интегралов в узловых значениях требуется исследование поведения передаточных функций при . В ходе его проведения было выявлено, что функции являются четными при и , а функции являются нечетными при и , где - некоторые константы. Функции при , а интегралы являются сходящимися. Рассмотрим отдельно вычисление интегралов: (5.3) а также (5.4) Отметим, что и являются сходящимися. Растущую часть интеграла из (5.3) будем рассматривать в главном значении в смысле Коши: Тогда интеграл (5.4) примет вид (5.5) Результаты вычислительных экспериментов по построению функций раскрытия для различных наборов характеристик покрытия и подложки будут представлены далее в соответствующем разделе. 6. Определение функции смещения на верхней границе полосы По найденным из системы (4.3) функциям раскрытия можно определить значения функций смещения на верхней границе слоя-покрытия. Трансформанты функций смещения слоя-покрытия в области отслоения из (4.1) имеют вид (6.1) Также трансформанты функции смещения верхнего слоя можно определить из решения второй вспомогательной задачи: (6.2) Подставляя (5.1) в (5.2) получим выражение Для верхней границы слоя: (6.3) Тогда, обращая преобразование Фурье, получим (6.4) Для вычисления используются те же принципы вычисления интегралов, что и для функций раскрытия . В следующем разделе также будут приведены результаты вычислительных экспериментов по определению функций смещения на верхней границе полосы. 7. Вычислительные эксперименты В данном разделе представлен анализ влияния характеристик покрытия и подложки на функции раскрытия, определяемые из решения интегральных уравнений (4.3), и функции смещения на верхней границе полосы, вычисляемые по формулам (5.4) на основе проведения соответствующих вычислительных экспериментов. В каждом из представленных наборов экспериментов соответствующие функции строились для различных значений характеристик покрытия (отношения параметров Ламе, толщины покрытия, уровня предварительных напряжений), а также для различных способов приложения нагрузки к верхней грани. Набор экспериментов 1 (анализ влияния уровня параметров Ламе покрытия). Представлены результаты расчета функций раскрытия (рис. 2) и функций смещения на верхней грани полосы (рис. 3) для различных соотношений ( ) параметров Ламе покрытия толщиной относительно параметров подложки в статическом случае . Из полученных графиков можно сделать вывод о том, что чем выше значения уровня параметров Ламе покрытия, тем меньше амплитуда функций раскрытия и смещения. Рис. 2. Графики функций раскрытия для различных соотношений параметров Ламе Fig. 2. Graphs of disclosure functions for various ratios of Lame parameters Рис. 3. Графики функций смещения на верхней границе полосы для различных соотношений параметров Ламе Fig. 3. Graphs of the displacement functions at the upper boundary of the strip for various ratios of Lame parameters Рис. 4. Графики функций раскрытия для различных толщин покрытия Fig. 4. Graphs of the disclosure functions for different coating thicknesses Рис. 5. Графики функций смещения на верхней границе полосы для различных толщин покрытия Fig. 5. Graphs of the displacement functions at the upper boundary of the strip for different coating thicknesses Набор экспериментов 2 (анализ влияния толщины покрытия). Представлены результаты расчета функций раскрытия (рис. 4) и функций смещения на верхней грани полосы (рис. 5) для различных толщин покрытия в статическом случае . В данном случае видно, что чем толще покрытие, тем меньше амплитуда функций раскрытия и смещения. Набор экспериментов 3 (анализ влияния уровня ПН). Представлены результаты расчета функций раскрытия (рис. 6) и функций смещения на верхней грани полосы (рис. 7) для различных уровней предварительных напряжений в статическом случае . В данном наборе экспериментов и в следующем рассматривалось сжимающее поле предварительных напряжений в покрытии. Рис. 6. Графики функций раскрытия для различных уровней предварительных напряжений Fig. 6. Graphs of the disclosure functions for various prestressing levels Рис. 7. Графики функций смещения на верхней границе для различных уровней предварительных напряжений Fig. 7. Graphs of the displacement functions at the upper boundary of the strip for various prestressing levels Рис. 8. Графики модуля АЧХ на верхней границе полосы для различных уровней предварительных напряжений Fig. 8. Graphs of the frequency response of at the upper edge of the strip for different prestress levels В данном примере заметим, что чем выше уровень сжимающих ПН в покрытии, тем меньше амплитуда функций раскрытия и смещения. При этом зависимость от уровня ПН, ввиду его малости относительно других характеристик, значительно меньше, чем в предыдущих экспериментах для других величин. Набор экспериментов 4 (анализ влияния уровня ПН на АЧХ). Исследована зависимость модулей АЧХ, измеренных на верхней грани полосы от уровня ( ) поля ПН в покрытии при в окрестности частоты первого толщинного резонанса. Графики АЧХ для различных уровней ПН приведены на рис. 8. Также отметим, что при увеличении уровеня сжимающих ПН в покрытии значения частот толщинных резонансов сдвигаются вправо. Заключение Проведенный анализ влияния механических и геометрических параметров покрытия, в частности уровня предварительных напряжений, на функции раскрытия отслоения и функции смещения и АЧХ на верхней грани полосы показывает, что изменение данных характеристик может служить диагностическим фактором для решения обратных задач идентификации уровня предварительных напряжений, вызвавших отслоение. Также стоит отметить, что наибольшее влияние на АЧХ наблюдается в окрестности частот толщинных резонансов, соответственно, данные частотные диапазоны являются наиболее благоприятными для построения схем идентификации характеристик отслоения, основанных на анализе акустического отклика.

About the authors

I V Bogachev

Southern Federal University

A O Vatulyan

Southern Federal University

References

Statistics

Views

Abstract - 403

PDF (Russian) - 309