Беляев А.Ю., Изюмов Р.И., Свистков А.Л. Влияние поверхностного заряда на формирование рельефа карбонизированного слоя на поверхности полимера при ионно-плазменной обработке // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2023. - № 2. С. 88-97. DOI: 10.15593/perm.mech/2023.2.08
Beliaev A.Yu., Izyumov R.I., Svistkov A.L. Influence of surface charge on the formation of the carbonized layer relief on the polymer surface during ion-plasma treatment. PNRPUMechanics Bulletin, 2023, no. 2, pp. 88-97. DOI: 10.15593/perm.mech/2023.2.08
пермскии политех
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 2,2023 PNRPU MECHANICS BULLETIN
http s://ered.p stu .ru/index.php/mechanics/index
Научная статья
Б01: 10.15593/регш.шесЬ/2023.2.08 УДК 539.3
ВЛИЯНИЕ ПОВЕРХНОСТНОГО ЗАРЯДА НА ФОРМИРОВАНИЕ РЕЛЬЕФА КАРБОНИЗИРОВАННОГО СЛОЯ НА ПОВЕРХНОСТИ ПОЛИМЕРА ПРИ ИОННО-ПЛАЗМЕННОЙ ОБРАБОТКЕ
А.Ю. Беляев1, Р.И. Изюмов2, А.Л. Свистков1
1 Институт механики сплошных сред Уральского отделения РАН -
филиал Пермского федерального исследовательского центра Уральского отделения РАН, Пермь, Россия
2Институт проблем механики им. А.Ю. Ишлинского РАН, Москва, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 15 октября 2022 г. Одобрена: 06 марта 2023 г. Принята к публикации: 30 апреля 2023 г.
Ключевые слова:
ионно-плазменная обработка, углеродный нанослой, численное моделирование, метод конечных элементов, упругость, устойчивость.
Рассматривается гипотеза происхождения волнообразного рельефа на поверхности полиуретана, обработанного плазмой. Выдвинуто предположение, что напряжения и деформации возникают из-за накопленного в приповерхностном слое одноименного заряда. Предложена методика расчета напряжений при условии равномерного распределения заряда. Получены определяющие уравнения упругой среды с распределенным зарядом, исходя из закона сохранения энергии и термодинамического неравенства. В полученных уравнениях тензор напряжений Коши содержит слагаемое, зависящее от плотности распределения заряда. Проведен расчет, показывающий зависимость величины напряжений от энергии и накопленной дозы. Расчет по предложенной модели показывает, что напряжения, возникающие в материале при некоторых дозах обработки достаточно высокие, чтобы вызвать изменение формы поверхности. Показано, что потеря устойчивости и, как следствие, возникновение волн на поверхности материала, характерно для низкомодульных полимеров. Проведено сравнение результатов расчета с реальными снимками поверхностей образцов после обработки, полученных с помощью оптического и атомно-силового микроскопов. Сделаны выводы о состоятельности выдвинутой гипотезы.
© ПНИПУ
© Беляев Антон Юрьевич - м.н.с., е-mail: belyaev@icmm.ru Изюмов Роман Игоревич - м.н.с., е-mail: izumov@icmm.ru Свистков Александр Львович - д.ф.-м.н., зав. лаб., е-mail: svistkov@icmm.ru
Anton Yu. Beliaev - Junior Researcher, е-mail: belyaev@icmm.ru Roman I. Izyumov - Junior Researcher, е-mail: izumov@icmm.ru Aleksandr L. Svistkov - Doctor of Physical and Mathematical Sciences, Head of С
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
Department, е-mail: svistkov@icmm.ru
INFLUENCE OF SURFACE CHARGE ON THE FORMATION OF THE CARBONIZED LAYER RELIEF ON THE POLYMER SURFACE DURING ION-PLASMA TREATMENT
A.Yu. Beliaev1, R.I. Izyumov2, A.L. Svistkov1
Institute of Continuous Media Mechanics UB RAS, Perm, Russian Federation Ishlinsky Institute for Problems in Mechanics RAS, Moscow, Russian Federation
ARTICLE INFO ABSTRACT
The origin hypothesis of the wavy relief on the surface of plasma-treated polyurethane is considered in this paper. It has been suggested that stresses and strains appear due to the eponymous charge accumulated in the near-surface layer. A technique for calculating stresses under the condition of a uniform charge distribution is proposed. The constitutive equations of an elastic medium with a distributed charge, based on the law of conservation of energy and thermodynamic inequality, are obtained. The Cauchy stress tensor contains a term that depends on the charge distribution density in the resulting equations. A calculation, showing the dependence of the magnitude of stresses on energy and fluence of implanted ions, has been carried out. According to the proposed model, calculations show that the stresses in the material are high enough to cause a change in the shape of the surface at certain treatment regimes. It is shown that the loss of stability and, as a consequence, the appearance of waves on the surface of the material is typical for low-modulus polymers. The calculation results are compared with real images, obtained using optical and atomic force microscopes, of the samples surfaces after the treatment. Conclusions about the viability of the proposed hypothesis are drawn.
© PNRPU
Received: 15 October 2022 Approved: 06 March 2023 Accepted for publication: 30 April 2023
Keywords:
ion-plasma treatment, carbon nanolayer, numerical simulation, finite element method, elasticity, stability.
Введение
Ионно-плазменная обработка полимеров является перспективным направлением в создании нанокомпо-зитных материалов. Такая обработка позволяет создавать на поверхности полимеров специфический слой с особыми физико-механическими свойствами [1-5]. Исследования с помощью инфракрасной и рентгеновской спектроскопии показывают высокое содержание углерода в слое. Получение таких слоев интересно с практической точки зрения, особый интерес вызывает использование таких материалов в медицинских целях [6-8]. Во-первых, высокое содержание углерода говорит о возможности получения электропроводящих слоев. Во-вторых, слой имеет достаточно высокую жесткость, в сравнении с исходным материалом [9]. В-третьих, слой имеет уникальные биофизические свойства, что коренным образом меняет поведение посаженных на обработанный материал биологических объектов (клетки, молекулы протеина, бактерии) [10-14].
Много вопросов возникает вокруг рельефа полученного слоя. Образовавшийся слой принимает волнистую форму [15; 16]. Амплитуда и ширина волн зависят от энергии обработки, накопленной дозы ионов и, судя по всему, от жесткости исходного материала. От полученного рельефа зависит трещиностойкость слоя и шероховатость поверхности [17; 18]. Оба фактора имеют значение при взаимодействии с биологическими объектами. В связи с этим возникает фундаментальная задача: объяснить причины образования волнообразного рельефа и научиться контролировать этот процесс для получения подходящих физико-механических свойств полученного материала.
Ряд работ посвящен расчетам, задачей в которых является установить границы появления волн, обладающих критическими значениями амплитуд и длин [19-21]. В других работах показано, что на микроуровне полиуретан содержит жесткие и мягкие блоки [22; 23], которые могут являться фактором формы потери устойчивости. При этом, как правило, фундаментальный вопрос о причинах появления напряжений, следствием которых является потеря устойчивости в материале, обычно не рассматривается.
Некоторые предыдущие попытки объяснить данный феномен давали совсем неубедительные результаты. Например, попытка объяснить феномен с позиции классической термоупругости или давления ионов на материал в процессе обработки. Другие попытки казались более перспективными. Рассмотрение процесса возможно с точки зрения неравновесной термодинамики с релаксирующими потоками тепла [24]. В работе показано, что в процессе обработки в приповерхностном слое возникают напряжения, которые могут быть причиной потери устойчивости. Однако тщательная численная проверка показала, что вклад напряжений вряд ли может считаться определяющим [27]. Подобный подход использовался и другими исследователями при рассмотрении схожих процессов, проходящих при малых временах и на малых масштабах [25; 26].
В данном исследовании осуществлен анализ разумности «электрической» гипотезы возникновения волнообразного рельефа. Так как ионы, которыми обрабатывается поверхность, имеют одноименный заряд, можно предположить, что при определенной накопленной дозе плотность заряда примет критическое значение. В ряде работ наличие заряда установлено экспериментально
[28; 29]. Под критическим значением понимается такое значение плотности заряда, когда силы электростатического отталкивания будут способны вызвать значительные напряжения и деформации материала в приповерхностном слое и произойдет потеря устойчивости. Теоретические основы использования теории для моделирования изгиба балки с учетом эффектов электроупругости изложены в работе [30]. В предложенной нами работе использован несколько иной подход.
В работе получены зависимости для сил, возникающих между равномерно распределенными ионами, и напряжений, возникающих вследствие этого. Проведена оценка возникающих напряжений. Показано, что результаты расчетов коррелируют с реальной картиной, полученной на оптическом и атомно-силовом микроскопе.
1. Электрическая гипотеза формирования волнообразной поверхности
Облучение материала происходит ионами азота. Ионы - это частицы, имеющие электрический заряд, соответственно при увеличении дозы обработки происходит накопление заряда в приповерхностном слое. Опишем влияние заряда на появление деформаций в материале.
Рассмотрим ровную поверхность, на которой расположены частицы с одинаковыми электрическими зарядами. Поскольку все заряды равны, то их расположение образует правильную треугольную решетку (рис. 1, а). Энергия взаимодействия Ж двух одинаковых зарядов q вычисляется по формуле Кулона
В рассматриваемом примере при N=1016 ион/см2 и И =80 нм расстояние между зарядами примет значение I = 1,02 нм. В реальности эта величина может быть больше, поскольку часть заряженных частиц может диффундировать глубже в полиуретановый материал.
W =
kq2
где г - расстояние между зарядами. Коэффициент к имеет значение 8,988-109 Н м2/Кл2. На рис. 1, а, штриховкой выделена шестиугольная область, которую имеет смысл рассматривать как окрестность одного электрического заряда. Этот заряд показан на иллюстрации в виде черного кружка. Между геометрическими размерами /, /ь /2, показанными на рисунке, имеется простая связь
/ = 2/2, /2 /!
2 2 2 ' 1 л/3
Ионы с определенной энергией проникают в полиуретан на некоторую глубину И. Таким образом, заряды распределены в материале в виде множества слоев. Можно считать, что места расположения зарядов находятся в узлах гексагональной решетки. Иллюстрация такой решетки, собранной для наглядности из касающихся друг друга шаров, показана на рис. 1, Ь.
Таким образом, выражение для / будет иметь вид
а b
Рис. 1. Выделенная область около одного заряда, используемая для нахождения плотности энергии взаимодействия электрических зарядов на плоской поверхности (а), гексагональная плотноупакованная решетка (b)
Fig. 1. Illustration of a selected region near a single charge, used to find the interaction energy density of electric charges on a flat surface (a), a hexagonal densely packed lattice (b)
2. Определяющие уравнения упругой среды с распределенным зарядом
Рассмотрим углеродный нанослой как очень тонкий слой около поверхности материала. Поскольку мы полагаем, что слой очень тонкий, то он должен напоминать оболочку, в которой важную роль в уравнениях должны играть касательные вектора к поверхности образца. Используем следующее упрощающее предположение, связанное с особенностями деформирования материала. Будем считать, что составляющая работы электрических зарядов при движении по направлению нормали при изменении толщины нанослоя, существенно меньше, чем составляющая работы, совершаемая зарядами при деформировании нанослоя в касательной плоскости. Поэтому этой составляющей можно пренебречь. В этом разделе будет предложена математическая модель, в которой электрические заряды приведут к появлению растягивающих двумерных напряжений в приповерхностном нанослое. Эти напряжения будут стремиться растянуть материал в касательной плоскости. Все приведенные ниже уравнения будут записаны в актуальной конфигурации.
В качестве первой исходной посылки для получения определяющих уравнений является формулировка первого начала термодинамики для деформируемого континуума
dj[pe + 2pv• v IdV =
l =
4h 1з 3N
= jp rdV + jT> • vdS + jTchn • OvdS - jn • qdS,
r
'y
где V, S, n - произвольно выделенный объем среды в углеродном слое, ограничивающая его поверхность и внешняя нормаль, t - текущий момент времени, e - массовая плотность внутренней энергии, р - плотность массы континуума, v - скорость точки континуума, q -тепловой поток, Tch - электрическая составляющая тензора напряжений Коши, возникающая в результате наличия заряженных частиц в деформируемой среде, Tel -упругая составляющая тензора напряжений Коши, определяющая упругую реакцию деформирования среды. В уравнении присутствует тензор O. Это ортопроектор. Он имеет следующие свойства.
O=OT = O2. (1)
Ортопроектор имеет вид
O = т1 ® т1 + т 2 ® т 2.
Символами т1 и т2 обозначены касательные к плоскости поверхности векторы. Это ортогональные векторы единичной длины. В предлагаемой модели использована посылка, что скорость совершения работы, вызванная электрическими силами, не зависит от движения точек нанослоя по направлению нормали к поверхности. Учитываются только проекции скоростей точек на касательную плоскость. Именно в касательной плоскости совершается работа электрическими зарядами. Поэтому в законе сохранения энергии появилось выражение Ov при учете скорости совершения работы электрическими зарядами.
Опуская промежуточные выкладки, остановимся на получаемых выводах. Справедливость второго закона термодинамики в любой инерциальной системе отсчета (инвариантность к преобразованиям Галилея) возможна только в случае, когда выполняется закон сохранения массы
Х1Х2 Х3 р = р0 = const,
уравнение движения
pv = div(Tel + OTch)
и более простая формулировка закона сохранения энергии, в которую не входит кинетическая энергия точек среды
pe = ( + OTch )• (D+W) - div q. (2)
Символами обозначены X1, X2, X3 - кратности удлинений материала, D и W - симметричная (тензор скоростей деформации) и антисимметричная часть в разложении градиента деформации grad v. Точка между
тензорами T и D означает скалярное умножение векторов девятимерного векторного пространства, которое образует множество тензоров второго ранга. Точка над символами означает взятие материальной производной. Это скорость изменения величины в выбранной точке континуума.
Второе начало термодинамики для деформируемого континуума формулируется в виде неравенства Клау-зиуса - Дюгема
ё
±\9sdv >-i ^ ds
Ht J J ft
ёг
где 5 - массовая плотность энтропии, 9 - температура. В дифференциальной форме оно принимает вид ограничение
ps > - divI ^^ I.
(3)
Из равенства (2) и ограничения (3) нетрудно получить важное термодинамическое неравенство, которое должно автоматически удовлетворяться в рамках предлагаемой математической модели
( + OTch )•( + W )-р f -ps9 - > 0,
(4)
где / = е - 9 5 - массовая плотность свободной энергии материала.
Нас интересуют упругие изотропные среды. Используем для формулировки массовой плотности свободной энергии следующее выражение
г
/ = Д (е,^^^) + ^О • Dёг.
'о
Величина с является функцией параметров состояния среды, которая зависит в том числе от плотности зарядов в нанослое. Неравенство (4) принимает в этом случае вид
(
3 5f т-pVL^i n. ®n.
\
D +
+ (OTCh-ap O) D + (Tel + OTch )• W -
(5)
s +-
f 59
grad 9^q
9
> 0.
Использование следующих определяющих уравнений является достаточным условием выполнения термодинамического ограничения в любом процессе в рассматриваемом материале. Упругая составляющая тензора напряжений Коши задается обычной для теории упругости формулой
3 Ы
е1 ^ ' а. • •
где п1,п2,п3 - ортонормированная тройка собственных векторов левого тензора растяжений. ТА - электрическая составляющая тензора напряжений Коши, имеет вид
ТА = стрО.
Поскольку тензор О является ортопроектором, выполняется равенство (1). То есть второе слагаемое в неравенстве (5) всегда равно нулю.
0ТсЬ-оО=о02 -оО = 0.
Третье слагаемое в неравенстве (5) равно нулю, так как является скалярным умножением симметричного на антисимметричный тензор. Энтропию материала следует вычислять, используя массовую плотность свободной энергии.
s = —
ее
При этом всегда должно выполняться ограничение §гаё < 0.
Это означает, что тепло всегда должно передаваться из горячих областей в холодные. В заключение отметим, что тензор напряжений Коши в рассматриваемой математической модели имеет вид:
3 д/
Т = Т + 0ТсЬ =р^.®п. +ор(т, ®т + т2 ®т2).
¿=1 дк,
Уравнение теплопроводности можно получить из равенства (2).
Для решения задач необходимо задать конкретный вид массовой плотности свободной энергии и электрической составляющей тензора напряжений Коши. Для эластомерного материала можно использовать потенциал:
fq = ^(V + V + V - 3).
Сложнее задать электрическую составляющую тензора напряжений Коши. Но ее легко получить в простом частном случае, когда все заряды располагаются на равном расстоянии друг от друга на некоторой поверхности и кривизна этой поверхности существенно больше расстояния между зарядами. В этом случае электрическая составляющая тензора напряжений Коши определяется выражением
Т =о(т1 ®т1 + т2 ®т2),
где а - величина, зависящая от плотности зарядов на поверхности, т1 и т2 - ортогональные касательные векторы единичной длины к плоскости поверхности. Это означает, что при решении задач на компьютере можно использовать стандартые программы, в которых тензор, отвечающий за механическое поведение материала, имеет вид:
3 д/
Т =-о(т1 ®т1 +т2 ®т2) + р^.®п.
¿=1 дКг
Первое слагаемое в правой части выражения дает деформацию материала вследствие распределенного заряда.
Рассмотрим конечную область материала. Вырежем тонкий слой образца около его границы с воздухом и заменим отброшенный материал силами, которые заме-
няют влияние удаленной части образца (рис. 2). На рисунке действие упругих сил и действие сил поверхностного натяжения показано маленькими черными стрелками. Красными стрелками показаны силы, которые создают заряды, располагающиеся на удаленной части материала. Поскольку система находится в равновесии, то их можно сосчитать как силы им противодействующие, взятые с противоположным знаком. Это силы, возникающие со стороны зарядов в выделенной части материала. Положение этих зарядов отмечено на рисунке черными кружками. Действуют они на крайние заряды, подписанные буквами А и В.
id.—
/////////////////// 'А • • • • • • в^*^
F
У///////////////// / г
Рис. 2. Иллюстрация, поясняющая геометрию выбранного участка материала и действующие силы, которые используются при построении математической модели
Fig. 2. Illustration explaining the geometry of the selected material section and the acting forces used to build the mathematical model
Абсолютную величину интересующей нас силы нетрудно определить. Вычисляется она по формуле
¥=¥ 14
при условии, что выделенная полоса имеет большую длину и число зарядов в ней можно считать очень большим. Значение суммы заряда известно.
V 1 ^ ^=~6~.
Следовательно,
F =
л2 kq2
Силу можно рассматривать как действующую на участок поверхности шириной Толщина
одного слоя электрических зарядов при их расположении в узлах гексагональной решетки, как уже было отмечено ранее, равна 0, /. То есть можно считать, что сила ¥ действует на сечение углеродного слоя площадью БР = 0,5>/3/■ 0,5>/э / = 0,75/2. Поэтому можно утверждать, что электрические заряды вызывают в карбонизированном слое появление сжимающих напряжений вдоль границы образца, определяемых формулой
=¥ = 2л2 kq2 °=^" 9/4 .
Ранее уже отмечалось, величина q равна заряду одного электрона
1,602-10 Кл, коэффициент к имеет значение 8,988-109 Н м2/Кл2. При накопленной дозе
N=1016 ион/см2 и кинетической энергии ионов 20 кэВ расстояние между зарядами равно I = 1,02 нм. В этом случае с = 4,67-108 Н/м2 = 467 МПа. Это очень большие напряжения, сжимающие тонкий углеродный слой толщиной 80 нм вдоль поверхности образца. В результате потери устойчивости углеродного слоя на массивном полиуретановом образце должна появиться волнообразная поверхность.
3. Численный эксперимент
Для осуществления вычислительных экспериментов необходимо знать упругие свойства материалов, задать геометрию рассматриваемой системы и сформулировать граничные условия. Модуль упругости материала был определен из стандартных механических испытаний и составил 10 МПа. Для нахождения толщины возникающего углеродного слоя использована программа 8Я1М. В ней рассматривается глубина проникновения частиц в материал. Она принимает во внимание химическую формулу материала и влетающих в него частиц. Безусловно, при таком подходе возникнет некоторая неточность, связанная с тем, что при инно-плазменной обработке меняется сам материал. В начале процесса это полиуретан. На завершении процесса это углеродный слой. Тем не менее такой анализ позволяет сделать необходимые оценки. Вычислительное моделирование с помощью программы 8Я1М, осуществленное в прошлые годы, позволило установить, что ионы азота с кинетической энергией 10 кеВ проникают в полиуретан на глубину 40 нм, а с энергией 20 кэВ - на глубину 80 нм.
Сложной проблемой является определение упругих свойств карбонизированного слоя. Для этого были использованы кривые индентирования полиуретанового материала с углеродным слоем (рис. 3).
При этом за основу брался наклон кривых при достаточной глубине внедрения зонда. На рисунке эти области показаны черным цветом. Такой прием позволил использовать для анализа те глубины внедрения зонда, при которых роль поверхностных эффектов в формировании силы противодействия зонду незначительна, и можно провести достаточно точный расчет силы противодействия внедрению зонда методом конечного элемента. Не потребуется в этом случае знать константы, характеризующие влияние поверхностных эффектов на особенности индентирования материала. При большом погружении зонда основную роль в противодействии внедрению играют в основном только упругие силы. Проведенные эксперименты позволили установить среднее значение начального модуля. Как видно из рис. 3, разбросы кривых индентирования достаточно большие. Это обусловлено тем, что полученный слой неоднороден: некоторые области могут быть недостаточно карбонизированы, при этом области с высокой степенью карбонизации содержат углерод в разных формах. Модуль упругости слоя, принятый для расчета составляет 1 ГПа.
14 12
100 -50 0 50 100 150
5, нм
a
100 -50 0 50 100 150
S, нм
b
Рис. 3. Кривые индентирования полиуретанового материала с карбонизированным слоем, которые использовались для нахождения упругого модуля слоя. Полиуретан с обработкой ионами азота с кинетической энергией 10 кэВ и дозой обработки: а - 1015 ион/см2; b - 1016 ион/см2
Fig. 3. Indentation curves of polyurethane material with the carbonized layer used to find the elastic modulus of the layer. Polyurethane treated with nitrogen ions with kinetic energy of 10 keV and fluence: a - 1015 ions/cm2; b - 1016 ions/cm2
I/; = 40 нм
а
1Л = 80 нм
- - ; \__I' ' ' J ~ -i-I 1 M KM _ ,
\ I 7\ I .
b
Рис. 4. Геометрия волнообразной поверхности углеродного слоя с толщиной: а - 40 нм; b - 80 нм
Fig. 4. Wave surface geometry of carbon layer with thickness of (a) 40 nm and (b) 80 nm
Расчеты выполнены методом конечного элемента в линейно упругой постановке. Решалась задача плос-кодеформированного состояния подложки с тонким
слоем на поверхности. В слое задавалась деформация и на разных шагах деформирования рассчитывалась сила, действующая в слое.
Результаты вычисления для разных толщин карбонизированного слоя показаны на рис. 4. Также проведены расчеты для разных модулей подложки и слоя.
4. Обсуждение результатов
Картина, полученная в результате численного эксперимента, качественно схожа с реальной. Наблюдается появление волн, размеры которых сопоставимы с размерами волн, измеренных с помощью АСМ (характерные длины волн в несколько микрон и высоты в десятки нанометров). К сожалению, нельзя сказать о хорошем численном совпадении численного эксперимента и реальной картины. Длины и амплитуды волн отличаются на 10-50 %, в зависимости от условий счета.
15 -10 -5 0 5 10 15
а
15 -10 -5 0 5 10 15
b
Рис. 5. Рельеф полиуретана после ионно-плазменной обработки снятый с помощью АСМ. Размеры области и перепады высот приведены в мкм. Накопленная доза 1016 ионов/см2 при энергии ионов 10 кэВ (а) и 20 кэВ (b)
Fig. 5. Relief of polyurethane after ion-plasma treatment taken by AFM. Area dimensions and height differences are given in ^m. Fluence is 1016 ions/cm2 at ion energies of 10 keV (a) and 20 keV (b)
Как можно заметить на рис. 4, а, в численном эксперименте наблюдается несколько мод потери устойчивости. Сначала проявляется мелкая мода, по мере увеличе-
ния деформации появляется вторая, более крупная. В численном эксперименте деформация определяет накопленную дозу ионов, так как чем выше доза, тем выше напряжения в слое. На рис. 5 приведены снимки поверхности обработанных полиуретанов с разной накопленной дозой. Можно заметить, что картина качественно схожа с результатами численного моделирования.
Расхождения численного эксперимента и картины, наблюдаемой для реальных материалов, обусловлено рядом упрощающих предположений, принятых в расчетах. Во-первых, углеродный материал не способен деформироваться на большую величину. Он разрушается при конечных деформациях. Иначе ведет себя полиуретан. Он может удлиниться при растяжении в несколько раз. При этом полиуретан размягчается, и в нем растут большие остаточные деформации. То есть под действием зарядов должен деформироваться только материал пока он еще не стал карбонизированным слоем. Следовательно, количество зарядов в процессе искривления поверхности образца может быть меньше, чем принятое исходя из накопленной дозы.
Во-вторых, по мере накопления зарядов около поверхности образца будет усиливаться электрическое поле, которое начнет тормозить движение ионов. При этом процесс торможения будет нарастать, ионы будут проникать в полиуретан на меньшую глубину. Указанное явление будет усиливаться в процессе обработки образца.
В-третьих, результаты численного моделирования показывают регулярную структуру складок, в то время как реальная картина достаточно хаотична, и геометрические характеристики (амплитуда и ширина) реальных волн имеют некоторое распределение. Связано это с тем, что реальная потеря устойчивости происходит на плоскости, и волны, идущие под разными углами, натыкаются друг на друга, оказывая влияние на конечную форму. Расчет в трехмерной постановке, в свою очередь, показывает картину, качественно схожую с реальной.
Рис. 6. Поле перемещений в трехмерной реализации численного эксперимента
Fig. 6. The displacement field in the three-dimensional numerical experiment implementation
Расчет для жестких подложек (модуль упругости 30 МПа) показал слабое изменение поверхности и отсутствие ярко выраженных волн, что тоже соответствует наблюдаемой картине на обработанных полиуретанах более жестких композиций.
Заключение
Предложена математическая модель, позволяющая связать накопленную дозу (количество ионов, попавших на единицу поверхности образца), толщину и модуль карбонизированного слоя с амплитудой и частотой волн, возникающих на границе материала. Дано обоснование накопления зарядов около поверхности образца в результате проникновения в него ионов азота в качестве физического объяснения причины появления волн.
Полученный в математической модели прогноз искривления карбонизированного слоя на полимерном образце находится в согласии с наблюдаемыми экспериментами. Возрастание жесткости карбонизированного слоя при увеличении накопленной дозы (в определенном диапазоне изменения) приводит к увеличению амплитуды волн и расстояния между вершинами волн.
Библиографический список
1. Oehr C. Plasma surface modification of polymers for biomédical use // Nucl. Instrum. Meth. B. - 2003. - Vol. 208. -P. 40-7.
2. Synthesis of biocompatible surfaces by nanotechnology methods / A.P. Alekhin, G.M. Boleiko, S.A. Gudkova, A.M. Mar-keev, A.A. Sigarev, V.F. Toknova, A.G. Kirilenko, R.V. Lapshin, E.N. Kozlov, D.V. Tetyukhin // Nanotechnologies in Russia 2010. -Vol. 5. - P. 696-708.
3. Surface characterization of plasma treated polymers for applications as biocompatible carriers / P. Slepicka, N.S. Kasal-kova, E. Stranska, L. Bacakova, V. Svorcik // eXPRESS Polymer Letters. - 2013. - Vol. 7. - P. 535-545. DOI: 10.3144/expresspolymlett.2013.50
4. Горошкина У.В., Лапуть О.А., Курзина И.А. Влияние ионно-плазменной обработки на поверхностные физико-химические свойства материалов на основе полилактида и гидроксиапатита // Вестник Томского государственного университета. Химия. - 2020. - № 19. - С. 45-50. DOI: 10.17223/24135542/19/5
5. Nageswaran G., Jothi L., Jagannathan S. Chapter 4 -Plasma Assisted Polymer Modifications // Non-Thermal Plasma Technology for Polymeric Materials. - 2019. - P. 95-127. DOI: 10.1016/B978-0-12-813152-7.00004-4
6. Plasma and Nanomaterials: Fabrication and Biomedical Applications / N.K. Kaushik, N. Kaushik, N.N. Linh, B. Ghimire, A. Pengkit, J. Sornsakdanuphap, S.-J. Lee, E.H. Choi // Nanomaterials. - 2019. - № 9. - Р. 98. DOI: 10.3390/nano9010098
7. Plasmas for medicine / T. Von Woedtke, S. Reuter, K. Masur, K.D. Weltmann // Phys. Rep. - 2013. - No. 530. - P. 291320. DOI: 10.1016/j.physrep.2013.05.005
8. Патент № 2721280 C1 Российская Федерация, МПК A61L 33/02, A61L 27/22, A61L 31/08. Способ создания анти-тромбогенного покрытия на полимерных сосудистых протезах: № 2019119336: заявл. 19.06.2019: опубл. 18.05.2020 / И.Н. Шардаков, В.С. Чудинов; заявитель Общество с ограниченной ответственностью «ИМБИОКОМ».
9. Kondyurin A., Volodin P., Weber J. Plasma Immersion Ion Implantation of Pebax Polymer // Nuclear Instruments and Methods B. - 2006. - Vol. 251, no. 2. - P. 407-412. DOI: 10.1016/j.nimb.2006.06.026
Аналогичная картина наблюдается при увеличении кинетической энергии ионов. Отсутствие волн на жестких подложках также находится в согласии с экспериментом.
При большом накоплении зарядов на поверхности образца происходит его пробой. Этот факт говорит о том, что электрическим явлениям необходимо уделять серьезное внимание.
Понимание причины появления данного явления дает инструмент для практического использования. Например, в результате увеличения размеров волн при деформировании материала будет происходить выравнивание углеродного слоя, и только при значительных деформациях в нем начнут появляться трещины, что положительно скажется на взаимодействии биологических объектов (клетки, молекулы протеина) с материалом.
10. Improvement of polyurethane surface biocompatibility by plasma and ion beam techniques / V. Melnig, N. Apetroaei, N. Dumitrascu, Y. Suzuki, V.J. Tura // Optoelectron. Adv. Mater. - 2005. - Vol. 7. - P. 2521.
11. Free radical functionalization of surfaces to prevent adverse responses to biomedical devices / M.M.M. Bilek, D.V. Bax, A. Kondyurin, Y.B. Yin, N.J. Nosworthy, K.P. Fisher // Proceedings of the National Academy of Sciences of the United States of America. - 2011. - Vol. 108. - P. 14405-14410.
12. Studying the Influence of Surface Topography on Bacterial Adhesion using Spatially Organized Microtopographic Surface Patterns / D. Perera-Costa, J.M. Bruque, M.L. González-Martín, A.C. Gómez-García, V. Vadillo-Rodríguez // Langmuir. - 2014. -Vol. 30. - P. 4633-41.
13. Елинсон В.М., Щур П.А. Технология формирования антимикробных фторуглеродных покрытий с углеродным подслоем на поверхности полимеров // Наноиндустрия. -2021. - Т. 14, № S6. - С. 211-220. DOI: 10.22184/1993-8578.2021.14.6s.211.220
14. Елинсон В.М., Щур П.А., Кравчук Е.Д. Параметры рельефа антиадгезионных антимикробных покрытий, сформированных методами ионно-плазменной технологии на поверхности ПЭТФ // Наноиндустрия. - 2020. - Т. 13, № S2. -С. 228-235. DOI: 10.22184/1993-8578.2020.13.2s.228.235
15. Control of Nanotexture and Wetting Properties of Polydimethylsiloxane from Very Hydrophobic to Super-Hydro-phobic by Plasma Processing / K. Tsougeni, A. Tserepi, G. Bou-lousis, V. Constantoudis, E. Gogolides // Plasma Processes and Polymers. - 2007. - Vol. 4, no. 4. - P. 398-405.
16. Croll A.B., Crosby A.J. Pattern driven localization in thin diblock copolymer films // Macromolecules. - 2012. -Vol. 45. - P. 4001-6.
17. Structural-mechanical and antibacterial properties of a soft elastic polyurethane surface after plasma immersion N2+ implantation / I.A. Morozov, A.S. Mamaev, I.V. Osorgina, L.M. Lemkina, V.P. Korobov, A.Yu. Belyaev, S.E. Porozova, M.G. Sherban // Materials Science & Engineering C. - 2016. -Vol. 62. - P. 242-248.
18. Evaluation of the deformation strength properties of a fluorinated polymer layer on the surface of poly(ethylene tereph-thalate) films / A.L. Volynskii, V.G. Nazarov, L.M. Yarysheva, V.P. Stolyarov, S.V. Tsypysheva, D.V. Bagrov, E.G. Rukhlya,
A.Yu. Yarysheva, N.F. Bakeev // Doklady Physical Chemistry. -2012. - Vol. 442. - P. 5-7. DOI: 10.1134/S0012501612010034
19. Nikravesh S., Ryu D., Shen Y.-L. Instabilities of Thin Films on a Compliant Substrate: Direct Numerical Simulations from Surface Wrinkling to Global Buckling // Scientific Reports. -2020. - Vol. 10. - P. 5728. DOI: 10.1038/s41598-020-62600-z
20. Nikravesh S., Ryu D., Shen Y.-L. Surface Instability of Composite Tin Films on Compliant Substrates: Direct Simulation Approach // Frontiers in Materials. - 2019. - Vol. 6. - P. 214. DOI: 10.3389/fmats.2019.00214
21. Folding wrinkles of a thin stiff layer on a soft substrate / J.-Y. Sun, S. Xia, M.-W. Moon, K. Hwan Oh, K.-S. Kim // Proc. R. Soc. - 2011. - A. 468932. - P. 953. DOI: 10.1098/rspa.2011.0567
22. Surface and subsurface AFM study of carbon-implanted polyurethane / I.A. Morozov, A.S. Kamenetskikh, A.Y. Beliaev, R.I. Izumov, M.G. Scherban, D.M. Kiselkov // Plasma Processes Polym. - 2022. - Vol. 19. - P. 2100156. DOI: 10.1002/ppap.202100156
23. Yilgör I., Yilgör E., Wilkes G.L. Critical parameters in designing segmented polyurethanes and their effect on morphology and properties: A comprehensive review // Polymer. - 2015. -Vol. 58. - P. 1-36. DOI: 10.1016/j.polymer.2014.12.014.
24. Müller I., Ruggeri T. Extended Thermodynamics V. 37. Springer Tracts in Natural Philosophy. - Springer, Berlin, 1993.
25. Lurie S.A., Belov P.A. On the nature of the relaxation time, the Maxwell-Cattaneo and Fourier law in the thermodynam-
References
1. Oehr C. Plasma surface modification of polymers for biomedical use // Nucl. Instrum. Meth. B. - 2003. - Vol. 208. -P. 40-7.
2. Synthesis of biocompatible surfaces by nanotechnology methods / A.P. Alekhin, G.M. Boleiko, S.A. Gudkova, A.M. Mar-keev, A.A. Sigarev, V.F. Toknova, A.G. Kirilenko, R.V. Lapshin, E.N. Kozlov, D.V. Tetyukhin // Nanotechnologies in Russia 2010. -Vol. 5. - P. 696-708.
3. Surface characterization of plasma treated polymers for applications as biocompatible carriers / P. Slepicka, N.S. Kasal-kova, E. Stranska, L. Bacakova, V. Svorcik // eXPRESS Polymer Letters. - 2013. - Vol. 7. - P. 535-545. DOI: 10.3144/expresspolymlett.2013.50
4. U.V. Goroshkina, O.A. Laput', I.A. Kurzina Vliyaniye ionno-plazmennoy obrabotki na poverkhnostnyye fiziko-khimicheskiye svoystva materialov na osnove polilaktida i gidroksiapatita. Vestnik Tomskogo gosudarstvennogo universiteta. Khimiya, 2020, no. 19, pp. 45-50. DOI: 10.17223/24135542/19/5
5. Nageswaran G., Jothi L., Jagannathan S. Chapter 4 -Plasma Assisted Polymer Modifications // Non-Thermal Plasma Technology for Polymeric Materials. - 2019. - P. 95-127. DOI: 10.1016/B978-0-12-813152-7.00004-4
6. Kaushik, N.K.; Kaushik, N.; Linh, N.N.; Ghimire, B.; Pengkit, A.; Sornsakdanuphap, J.; Lee, S.-J.; Choi, E.H. Plasma and Nanomaterials: Fabrication and Biomedical Applications. Nanomaterials, 2019, no. 9, p. 98. DOI: 10.3390/nano9010098
7. Von Woedtke, T.; Reuter, S.; Masur, K.; Weltmann, K.D. Plasmas for medicine. Phys. Rep., 2013, no. 530, pp. 291-320. DOI:10.1016/j.physrep.2013.05.005
8. Patent No 2721280 C1 Russian Federation, IPC A61L 33/02, A61L 27/22, A61L 31/08. Method for creation of antithrombogenic coating on polymeric vascular prostheses : No. 2019119336 : application. 19.06.2019 : publ. 18.05.2020 /
ics of a continuous medium, and the scale effects in thermal conductivity // Continuum Mech. Thermodyn. - 2020. - Vol. 32. -P. 709-728. DOI: 10.1007/s00161-018-0718-7
26. Wave reflection at a free interface in an anisotropic pyroelectric medium with nonclassical thermoelasticity / Aen. Abd-alla, I. Giorgio, L. Galantucci [et al.] // Continuum Mech. Thermodyn. - 2016. - Vol. 28. - P. 67-84. DOI: 10.1007/s00161-014-0400-7
27. Комар Л.А., Свистков А.Л., Беляев А.Ю. Метод решения линеаризированных тепловых задач с учетом явления релаксации теплового потока // Вычислительная механика сплошных сред. - 2022. - № 15(2). - С. 185-192. DOI: 10.7242/1999-6691/2022.15.2.14
28. Kukushkin Surface charge of polymer materials modified by nanostructured fluorocarbon coatings / V.M. Elinson [et al.] // J. Phys.: Conf. Ser. - 2020. - Vol. 1713. - P. 012016. DOI: 10.1088/1742-6596/1713/1/012016
29. Влияние поверхностного заряда электретов на грибо-стойкость фторуглеродных полимерных материалов / В.М. Елин-сон, А.С. Аболенцев, Т.В. Ходырев, П. А. Щур // Вакуумная наука и техника: материалы XXVIII Научно-технической конференции с участием зарубежных специалистов, Крым, Судак, 16-21 сентября 2021 года. - Судак, 2021. - С. 260-267.
30. Lurie S., Solyaev Y. On the formulation of elastic and electroelastic gradient beam theories // Continuum Mech. Thermodyn. - 2019. - Vol. 31. - P. 1601-1613. DOI: 10.1007/s00161-019-00781-3
I.N. Shardakov, V.S. Chudinov ; applicant Imbiokom Limited Liability Company.
9. Патент № 2721280 C1 Российская Федерация, МПК A61L 33/02, A61L 27/22, A61L 31/08. Способ создания анти-тромбогенного покрытия на полимерных сосудистых протезах: № 2019119336: заявл. 19.06.2019: опубл. 18.05.2020 / И.Н. Шардаков, В.С. Чудинов; заявитель Общество с ограниченной ответственностью «ИМБИОКОМ».
10.Kondyurin A., Volodin P., Weber J. Plasma Immersion Ion Implantation of Pebax Polymer // Nuclear Instruments and Methods B. - 2006. - Vol. 251, no. 2. - P. 407-412. DOI: 10.1016/j.nimb.2006.06.026
11. Bilek, M. M. M.; Bax, D. V.; Kondyurin, A.; Yin, Y. B.; Nosworthy, N. J.; Fisher, K. P. et all. Free radical functionalization of surfaces to prevent adverse responses to biomedical devices. Proceedings of the National Academy of Sciences of the United States of America, 2011, vol. 108, pp. 14405-14410.
12. Studying the Influence of Surface Topography on Bacterial Adhesion using Spatially Organized Microtopographic Surface Patterns / D. Perera-Costa, J.M. Bruque, M.L. González-Martín, A.C. Gómez-García, V. Vadillo-Rodríguez // Langmuir. - 2014. -Vol. 30. - P. 4633-41.
13. V.M. Yelinson, P.A. Shchur Tekhnologiya formiro-vaniya antimikrobnykh ftoruglerodnykh pokrytiy s uglerodnym podsloyem na poverkhnosti polimerov [Technology of formation of antimicrobial fluorocarbon coatings with a carbon sublayer on the surface of polymers]. Nanoindustriya, 2021, vol. 14, no. S6, pp. 211-220. DOI 10.22184/1993-8578.2021.14.6s.211.220.
14. V.M. Yelinson, P.A. Shchur, Ye.D. Kravchuk Parametry rel'yefa antiadgezionnykh antimikrobnykh pokrytiy, sformiro-vannykh metodami ionno-plazmennoy tekhnologii na poverkhnosti PETF [Relief parameters of anti-adhesive antimicrobial coatings formed by ion-plasma technology on the surface of PET].
Nanoindustriya, 2020, vol. 13, no. S2, pp. 228-235. DOI 10.22184/1993-8578.2020.13.2s.228.235.
15. Control of Nanotexture and Wetting Properties of Polydimethylsiloxane from Very Hydrophobic to Super-Hydro-phobic by Plasma Processing / K. Tsougeni, A. Tserepi, G. Bou-lousis, V. Constantoudis, E. Gogolides // Plasma Processes and Polymers. - 2007. - Vol. 4, no. 4. - P. 398-405.
16. Croll A.B., Crosby A.J. Pattern driven localization in thin diblock copolymer films // Macromolecules. - 2012. -Vol. 45. - P. 4001-6.
17. Structural-mechanical and antibacterial properties of a soft elastic polyurethane surface after plasma immersion N2+ implantation / I.A. Morozov, A.S. Mamaev, I.V. Osorgina, L.M. Lemkina, V.P. Korobov, A.Yu. Belyaev, S.E. Porozova, M.G. Sherban // Materials Science & Engineering C. - 2016. -Vol. 62. - P. 242-248.
18. Evaluation of the deformation strength properties of a fluorinated polymer layer on the surface of poly(ethylene terephthalate) films / A.L. Volynskii, V.G. Nazarov, L.M. Yarysheva, V.P. Stolyarov, S.V. Tsypysheva, D.V. Bagrov, E.G. Rukhlya, A.Yu. Yarysheva, N.F. Bakeev // Doklady Physical Chemistry. - 2012. - Vol. 442. - P. 5-7. DOI: 10.1134/S0012501612010034
19. Nikravesh S., Ryu D., Shen Y.-L. Instabilities of Thin Films on a Compliant Substrate: Direct Numerical Simulations from Surface Wrinkling to Global Buckling // Scientific Reports. -2020. - Vol. 10. - P. 5728. DOI: 10.1038/s41598-020-62600-z
20. Nikravesh S., Ryu D., Shen Y.-L. Surface Instability of Composite Tin Films on Compliant Substrates: Direct Simulation Approach // Frontiers in Materials. - 2019. - Vol. 6. - P. 214. DOI: 10.3389/fmats.2019.00214
21. Folding wrinkles of a thin stiff layer on a soft substrate / J.-Y. Sun, S. Xia, M.-W. Moon, K. Hwan Oh, K.-S. Kim // Proc. R. Soc. - 2011. - A. 468932. - P. 953. DOI: 10.1098/rspa.2011.0567
22. Morozov, I.A., Kamenetskikh, A.S., Beliaev, A.Y., Izumov, R.I., Scherban, M.G., Kiselkov, D.M. Surface and subsurface AFM study of carbon-implanted polyurethane. Plasma Pro-
cesses Polym, 2022, vol. 19, 2100156. DOI: 10.1002/ppap.202100156
23. Yilgör I., Yilgör E., Wilkes G.L. Critical parameters in designing segmented polyurethanes and their effect on morphology and properties: A comprehensive review // Polymer. - 2015. -Vol. 58. - P. 1-36. DOI: 10.1016/j.polymer.2014.12.014.
24. Müller I., Ruggeri T. Extended Thermodynamics V. 37. Springer Tracts in Natural Philosophy. - Springer, Berlin, 1993.
25. Lurie S.A., Belov P.A. On the nature of the relaxation time, the Maxwell-Cattaneo and Fourier law in the thermodynamics of a continuous medium, and the scale effects in thermal conductivity // Continuum Mech. Thermodyn. - 2020. - Vol. 32. -P. 709-728. DOI: 10.1007/s00161-018-0718-7
26. Wave reflection at a free interface in an anisotropic pyroelectric medium with nonclassical thermoelasticity / Aen. Abd-alla, I. Giorgio, L. Galantucci [et al.] // Continuum Mech. Thermodyn. - 2016. - Vol. 28. - P. 67-84. DOI: 10.1007/s00161-014-0400-7
27. Komar, L.A., Svistkov, A.L., Beliaev, A.Yu. Method for solving linearized thermal problems taking into account the phenomenon of heat flow relaxation. Computational Continuum Mechanicsthis, 2022, no. №15(2), pp. 185-192. DOI: 10.7242/1999-6691/2022.15.2.14
28. V.M. Elinson, P.A. Shchur, D.Yu. Kukushkin Surface charge of polymer materials modified by nanostructured fluoro-carbon coatings. J. Phys.: Conf. Ser., 2020, vol. 1713, 012016. DOI 10.1088/1742-6596/1713/1/012016
29. V.M. Yelinson, A.S. Abolentsev, T.V. Khodyrev, P.A. Shchur Vliyaniye poverkhnostnogo zaryada elektretov na gribostoykost' ftoruglerodnykh polimernykh materialov [The effect of the surface charge of electretes on the fungal resistance of fluo-rocarbon polymer materials]. Vakuumnaya nauka i tekhnika: Materialy XXVIII nauchno-tekhnicheskoy konferentsii s uchastiyem zarubezhnykh spetsialistov, Krym, Sudak, 2021.. 260-267.
30. Lurie S., Solyaev Y. On the formulation of elastic and electroelastic gradient beam theories // Continuum Mech. Thermodyn. - 2019. - Vol. 31. - P. 1601-1613. DOI: 10.1007/s00161-019-00781-3
Финансирование. Работа выполнена при финансовой поддержке Российского научного фонда (проект № 18-19-00574-П). Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов. Вклад авторов равноценен.
Financing. The work was supported by RSF Grant 18-19-00574-P. Conflict of interest. The authors declare no conflict of interest. The contribution of the authors is equivalent.