10. Kovalenko, V. I., Marynin, V. G. (1998). Obladnannja dlja doslidzhennja erozii' pokryttiv pry mikro udarnomu dijanni. Voprosy atomnoj nauky y tehnyky, 5, 83-89.
11. Marinin, V. G., Kovalenko, V. I. (2013). Cavitation resistance of vacuum-arc coverings deposited on substrate of different orientation. Technology audit and production reserves, 4/2 (12), 15-18. Available at: http://journals.uran.ua/tarp/ article/view/16243/13672
12. Pogrebnjak, A. D., Shpak, A. P., Azarenkov, N. A., Beresnev, V. M. (2009). Struktura i svojstva tverdyh i sverhtverdyh nanokompozitnyh pokrytij. Uspehi fizicheskih nauk, 179 (1), 35-64. doi: 10.3367/ufnr.0179.200901b.0035
13. Musiel, J., Vlcek, J., Jezek, V., Benda, M., Kole-ga, M., Boomsma, R. (1995). Production of Ti films with controlled texture. Surface and Coatings Technology, 76-77, 274279. doi: 10.1016/0257-8972(95)02582-0
Дата надходження рукопису 18.04.2016
Андреев Анатолий Афанасьевич, доктор технических наук, старший научный сотрудник, ведущий научный сотрудник отдела, Отдел интенсивных вакуумно-плазменных технологий, Институт физики твердого тела, материаловедения и технологий, Национальный научный центр «Харьковский физико-технический институт», ул. Академическая, 1, г. Харьков, Украина, 61108 E-mail: aandreev@kipt.kharkov.ua
Коваленко Владимир Иванович, научный сотрудник отдела, Отдел интенсивных вакуумно-плазменных технологий, Институт физики твердого тела, материаловедения и технологий, Национальный научный центр «Харьковский физико-технический институт», ул. Академическая, 1, г. Харьков, Украина, 61108
E-mail: kartmazov@kipt.kharkov.ua
Маринин Владимир Григорьевич, кандидат физико-математических наук, доцент, старший научный сотрудник, Отдел интенсивных вакуумно-плазменных технологий, Институт физики твердого тела, материаловедения и технологий, Национальный научный центр «Харьковский физико-технический институт», ул. Академическая, 1, г. Харьков, Украина, 61108 E-mail: marinin@kipt.kharkov.ua
Мартыненко Людмила Ивановна, заместитель начальника отдела, Отдел интенсивных вакуумно-плазменных технологий, Институт физики твердого тела, материаловедения и технологий, Национальный научный центр «Харьковский физико-технический институт», ул. Академическая, 1, г. Харьков, Украина, 61108
Столбовой Вячеслав Александрович, кандидат технических наук, старший научный сотрудник отдела, Отдел интенсивных вакуумно-плазменных технологий, Институт физики твердого тела, материаловедения и технологий, Национальный научный центр «Харьковский физико-технический институт», ул. Академическая, 1, г. Харьков, Украина, 61108
УДК 519.64.539.3
DOI: 10.15587/2313-8416.2016.69828
МЕТОД ЭНЕРГЕТИЧЕСКОГО МОДЕЛИРОВАНИЯ В ПРИМЕНЕНИИ К ДИФРАКЦИИ СТАЦИОНАРНЫХ ПЛОСКИХ УПРУГИХ ВОЛН
© А. М. Назаренко, В. О. Мищенко
Энергетический метод использует теорему взаимности работ и преобразование возникающих гиперсингулярных интегралов для получения системы уравнений, состоящей из сингулярных интегральных уравнений первого и второго рода. Алгоритмы численного решения получаются методом дискретных особенностей. Точность вычислений контролируется за счёт алгоритмической диверсности. Показано приложение к численному моделированию дифракции стационарных плоских волн на конечных и периодических системах неоднородностей разных типов
Ключевые слова: дифракция упругих волн, неоднородности, метод энергетического моделирования, плоская деформация, напряжения
The energy method uses reciprocity works theorem and transformation of occurring hypersingular integrals for a system of equations consisting of singular integral equations of the first and second kind. Algorithms for the numerical solution are obtained by the method of discrete singularities. Precision of calculations is controlled due to algorithmic diversity. Application to the numerical modeling of diffraction of stationary plane waves on finite and periodic systems of different inhomogeneities is shown
Keywords: diffraction of elastic waves, inhomogeneities, energy modeling method, plane strain, stresses
1. Введение
В практике проектирования конструкций, работающих в условиях динамического, как правило, гармонически изменяющегося нагружения, необходимо учитывать наличие в материале различных концентраторов напряжений (трещин, включений, отверстий). Часто при этом физические методы исследования готовых конструкций или проведение экспериментального исследования образцов невозможны или непозволительно дороги. В последнее время всё более востребованным становится численное (компьютерное) моделирование, которое тем привлекательнее для инженеров и конструкторов, чем больше альтернативных вариантов постановки задачи может решаться используемым методом и чем проще перенастраивать программную систему для проведения вычислительных экспериментов.
Поэтому актуально создание и компьютерная реализация таких методов математического моделирования процессов дифракции упругих волн на неодно-родностях в твёрдом теле, которые универсально применимы к разным типам волновых процессов и разным типам неоднородностей.
2. Анализ публикаций и постановка проблемы
Хорошо известны следующие методы математического моделирования дифракции упругих волн на различного типа неоднородностях: метод разделения переменных [1], методы конечных разностей [2], методы конечных элементов [3], асимптотические методы [4]. Возможности их применения в некоторых важных для практики случаях имеют ограничения. Основные результаты в последние годы получены на основе СИУ [5]. При этом после сведения задачи к СИУ, её дискретную модель можно получить методом механических квадратур или методами дискретных особенностей (МДО) [6-8]. Замечательно, что компьютерное моделирование дифракционных процессов на основе МДО может быть ускорено на порядок даже на персональных компьютерах за счёт исследования их архитектурных особенностей [9].
Существующие методы обычно используют те или иные интегральные представления функций, выбранных для решения конкретной дифракционной задачи. В случае плоской деформации распространёнными методами исследования дифрагированного волнового поля являются метод разрывных решений [10] и метод специальных интегральных представлений волновых потенциалов [11, 12]. Однако в этих методах способы построения интегральных представлений разрешающих функций, характеризующих возмущенное волновое поле, зависят от типа неодно-родностей. Возникают также специфические трудности, если одновременно учитываются неоднородности разных типов.
Вследствие указанных обстоятельств ряд задач дифракции плоских волн в упругом теле с неоднородными включениями разных видов в 2D постановке, в частности, периодические задачи, оставались вовсе или недостаточно изученными другими учё-
ными с точки зрения распределения напряжений на границах неоднородностей. Благодаря единому энергетическому методу моделирования мы имеем возможность выполнять исследование, не выделяя отдельные типы неоднородностей и волн вплоть до момента записи конкретной системы граничных сингулярных интегральных уравнений. Потом, при записи таких уравнений данный метод позволяет прозрачно учесть специфику краевых условий на различных неоднородностях. Т. е., исследование систем неоднородностей разного типа в задачах дифракции плоских волн на цилиндрических неоднородностях не только актуально, но найденным нами методом выполнимо.
3. Исходная краевая задача и задачи исследования
С целью моделирования дифракционного процесса рассмотрим в неограниченной среде систему бесконечных вдоль оси Х3 цилиндрических включений с ограниченными замкнутыми или разомкнутыми образующими (рис. 1).
Рис. 1. Система неоднородностей в неограниченной упругой среде: а - жесткое включение; б - трещина; в - полость; г - упругое включение
Пусть из бесконечности на такую систему набегает монохроматическая гармоническая волна (зависимость от времени г~т') расширения-сжатия (Р - случай):
U(0'> = 0, U2(0) = Txe
■ a)v X 2
, Xi(1) =®/c(1), Tx = const (1)
„(1)
или сдвига (БУ- случай):
ц(0) =г2e-'* \ Uf = 0,
П
(1) =©/c(1),
г2 = const,
(2)
где с[' и с2 - скорости Р- и БУ- волн в матрице.
При взаимодействии приходящей волны с цилиндрами возникают отраженные и проходящие внутрь цилиндров (если они являются упругими включениями) волны двух типов: продольные и поперечные, причем другие типы волн не образуются. Пусть и[к > и и(к} - амплитуды перемещений отра-
женного (£=1) и проникающего внутрь упругих включений (£=2) полей соответственно. Тогда общее
г
поле перемещений и = (и1, и2) равно
г
(Ц(1) + Ц(0), и® + —о))
г
в матрице и (Ц(2), и22)) внутри упругих включений.
Компоненты вектора и удовлетворяют уравнениям движения [19]:
/Ди + (1 + и)ЯГасШуи + рс2и + / = Л[и] + f = 0, (3)
где 1 и и - постоянные Ламе; р - плотность среды;
г
со - частота колебаний; / = (/, /) - объемная сила.
Амплитудные значения напряжений связаны с амплитудами перемещений с помощью формул
„д(и. + ¡и.) д(Ц -Ц )ч стп +а22 = 2(1 + /)(-^ + ( ' _ 2) ) , (4)
5z
eg7
д(Ц - /Ц2) ez '
СТ22 - СТ11 - 2'СТ12 = —4^
д (Ц + iü2)
dz
Пусть Ь - некоторая кривая на плоскости Х1Х2. Обозначим через ^ и амплитуды тангенциальной и нормальной компонент вектора напряжений $ на Ь . Тогда в произвольной точке кривой д0 е Ь эти напряжения выражаются через компоненты тензора напряжений:
2г(^ + г) = (стп +ст22) в'Фо + (ст22 - стп -2гст12) е-фо, (5)
—2г (^ - ) = (СТ! + ст22) е-* + (^ - ^ + 2СТ2) е'Фо.
Здесь ф0 - угол положительной касательной к Ь в точке д0 е Ь с осью Х1.
На границе замкнутых неоднородностей представляют интерес также распределения амплитуд напряжений
СТ = ^'51пфо - ^2С05фо, СТ =СТ" +^22 -СТ
= Sicos Фо + S2 sin Фо
(6)
связанных с координатами 50 и п0 на Ь .
Сформулируем теперь граничные условия на контурах неоднородностей. Пусть д0=^0+1 щ -
точка кривой Ь , в которой мы будем удовлетворять граничные условия.
1. На границе раздела двух сред (матрица -упругое включение) моделируем условия контакта типа склейки - непрерывность при переходе через границу Ь области включения вектора перемещений
г г
и = (и , и2) и вектора напряжений Б = (^ , ) :
(и-1> + иГ ) Г ) г ^
№+¿г))т=1,2 (7)
Соотношения (7) следует понимать как равенства граничных значений соответствующих величин при стремлении к точке д0=£0+ гщ е Ь точки наблюдения 2 = х1 + 1х2, которая при этом всё время должна оставаться вне (область Д) или внутри (область Д) границы Ь.
2. На границе Ь неподвижного включения перемещения обращаются в ноль:
(u(1) + U(0)) = 0, z еД, m = 1,2.
(8)
3. На контуре Ь жесткого подвижного включения граничные условия вытекают из уравнений движения абсолютно жесткого тела:
(u« + ui0)) = в -0ш%, (Ц1 + u(0)) = B2 +
(9)
Здесь предполагается, что включение с границей Ь перемещается и поворачивается вместе с матрицей; В и В - амплитуды поступательного движения, ав - амплитуда жесткого поворота этого включения.
В случае (9) на основании второго закона Ньютона для поступательного движения включения имеем:
| ^Ио = -дВт, д = ®2р £ , т = ',2. (10)
Ь
Уравнение, описывающие вращательное движение жесткого включения, имеет вид:
(Щ - Х2о)- ^2 (^о - ХЮ)) С?о = -с2■7в С , (11)
где ре - плотность, - площадь, - момент инерции включения Д относительно произвольной точки Л(х10, Хо).
Дополнительные соотношения (10), (11) необходимо использовать для определения неизвестных постоянных В , В и С в (9).
4. На границе Ь полости обращаются в нуль компоненты вектора напряжений:
( s? + S? )
= 0, z e D
fo
iL, m = 1,2.
(12)
5. В случае криволинейного жесткого включения граничные условия на L имеют вид (9), а неизвестные постоянные, характеризующие поступательные и вращательное движения данного включения определяются из условия равенства нулю результирующих сил и момента сил, действующих на L . Берега каждой криволинейной трещины-разреза предполагаются свободными от сил (12), а необходимые дополнительные условия для однозначной разрешимости СИУ первого рода будут вытекать из условия равенства нулю скачков перемещений в вершинах этой трещины.
Однозначная разрешимость краевых задач для уравнений (3) для отраженного и проникающего полей при любом наборе конечной системы неоднородностей обеспечивается условиями подобия отраженного поля расходящейся цилиндрической волне (условие на бесконечности типа Зоммерфельда) и локальной принадлежности Соболевскому классу W22.
3. Методы исследования
Энергетический метод моделирования, который не сразу достиг методической зрелости, фактически сложился в период с 1990 по 2015 г. в ряде работ, начиная с [13], и к настоящему времени уже позволил построить вычислительные методы для широких классов задач (например, [13-19]). Отличительной чертой этого метода является построение интегральных представлений перемещений дифрагированного волнового поля с помощью теоремы взаимности работ [20] для двух различных состояний механической системы с тем, чтобы выразить указанные перемещения в интегральной форме через кинематические и силовые величины. Таким способом удаётся получить интегральные представления перемещений, обладающие желательными свойствами.
Теорема взаимности работ для систем уравнений (3) представляет собой аналог формулы Грина в теории уравнения Пуассона-Лапласа. Она применяется нами для получения сингулярной системы граничных интегральных уравнений по следующей схеме. Введём над вектор-функциями перемещений дифференциальную операцию
!
А[И] = (Л1[^1], Л^]) .
Тогда для замкнутого контура Ь = дD, лежащего в области выполнения (3) для пары перемещений и, V, соответствующих полям объёмных сил Г, g, по обобщенной теореме Стокса имеем
(ш , И = -/ (Л1 V] и +Л 2 [г2] и2)ау =
D
= -[(VI-Л1 +и2-Л2 [иг])ск +
где 3 (V, И) - билинейная форма, представляющая собой интеграл по контуру от произведений компонент перемещений на их частные производные.
Подынтегральное выражение можно привести к виду произведений компонента перемещений на компоненты вектора напряжений Б (5). Выбор в качестве V фундаментального (параметризованного координатами особенности) решения, «подправленного» для выполнения граничных условий на границах неоднородностей (специальная функция Грина конкретной задачи), превращает (13) в интегральное представление для и. Его следует дважды применить, выбирая контуры для границ областей так, как показано для примера на рис. 2 (ABCDEFGHIA и AIHGFEDCBA). Точнее, вначале необходимо отступить от берегов разрезов в глубь областей на е .
Ясно, что для одной из областей (в которую особенность не попала) интеграл по замкнутому контуру обнулится, но вычитание полученных выражений приведёт к новому интегральному представлению, в котором после предельного перехода Я ^ да можно избавиться от интегралов по окружности радиуса Я . Однако затруднение, которого старались избегнуть многие исследователи данной предметной области, состоит в том, что интегралы, содержащие производные функции Грина, не имеют конечных пределов при натягивании контуров на границы неоднородностей. Однако, следуя известному приёму Адамара, можно добавлять и вычитать выражения, которые при предельном переходе е ^ 0 дадут т. н. конечные части интегралов, а тогда окажется, что все «мешающие» члены представления уходят. Более того, удаётся обосновать возможность «перебросить» производные с функции Грина на перемещение, с результатом, как при использовании обобщенной теоремы Стокса.
Рис. 2. Два замкнутых контура с общими участками и дуговыми участками радиуса Я , имея в виду устремить потом Я ^ да
В результате для неоднородностей рассматриваемых типов дифрагированное поле перемещений представимо в виде (предполагая суммирование по повторяющемуся индексу):
+3l (V, U) = (V, f) + 3l (V, U),
(13)
um (P) = i ff)(s) Gmn (P, Q)ds +
L(1)
(i fn (s) Gmn (p, Q)ds - j (s) smn (P, Q)ds),
D
X = 5 ], & = и ], к, т = 1,2. (14)
Здесь Ьт , Ь(2), Ь3 и Ь(4) - совокупности границ упругих включений, полостей, жестких включений и криволинейных трещин соответственно; 81к -символ Кронекера; перемещения вычисляются в точке области Р( х, х2), а интегрирование проводится по точкам Qп); [ии ] - скачки амплитуд, перемещений на трещине совокупности Ь(4), [5и ] - скачки амплитуд напряжений на криволинейном жестком включении совокупности Ь3) (п=1,2); О® - амплитуда перемещения вдоль оси Хт при действии гармонической сосредоточенной силы, приложенной в точке Q и направленной вдоль оси Хп (п-е состояние); - компоненты вектора напряжений Б, соответствующие перемещениям О®.
Удовлетворение граничных условий по перемещениям на соответствующей совокупности контуров приводит к интегральным уравнениям с логарифмическими особенностями. С целью получения СИУ 1-го рода указанные уравнения дифференцировались по дуговой координате. Необходимые дополнительные условия для однозначной разрешимости СИУ 1-го рода вытекают из эквивалентности последних и соответствующих интегральных уравнений с логарифмическими особенностями в случае упругих включений; из (10), (11) в случае замкнутых жестких включений и из равенства нулю главного вектора сил и момента сил, действующих на криволинейном жестком включении совокупности Удовлетворение граничных условий по напряжениям приводит к СИУ 2-го рода в случае упругих включений и полостей и к СИУ 1 -го рода в случае криволинейной трещины совокупности Ь(4).
В качестве примера рассмотрим 2^-перио-дическую решетку, состоящую из однородных полостей или жестких включений (их границы входят в совокупность Ь(2) или Ь3)). В этом случае в интегральном представлении (14) О® - элементы периодической матрицы Грина. Далее верхний индекс опускаем и запишем выражения для следующих комбинаций элементов матрицы Грина (суммирование по п=1,2):
G11 + G22 = hj
» О+А^^-А cos«
i=0 Kni
(X -)
При указанном в (15) выборе знаков для Яи
поле, порождаемое периодической системой источников, носит характер расходящихся волн, что соответствует условиям излучения на бесконечности.
Присутствующий в (15) ряд, соответствующий Оц + О22, в точке приложения сосредоточенного источника (Р = Q) расходится (общий член ряда ведет себя как 1/1). Для выделения логарифмической особенности и ускорения сходимости функциональных рядов-остатков выделим в комбинациях (15) статические члены, соответствующие значению ю=0. Для них получены следующие выражения:
G
101 + G02 = ho (jj cos « (Xi -))- — \x2 - i \)
i=1 ai
G11 - G00 ± 2iG10 = -
-2ho\ X2 -i\ X(j e
,(X -#)-
H— + isign(x2 - rj) x ''sina, (xj), (16)
Здесь И, к0, % - константы, зависящие от упругих постоянных среды.
Фигурирующие в (16) функциональные ряды можно просуммировать в явном виде. Находим
G11 + G22 = -
2ho —d
Re(ln2sin-
л( z
(z-С)
G101 - G22 + 2iG102 = h0 ' (X2 - i) Ctg -
2d ^ 2d
G° - G22 - 2/G12 =
= -h' (*2 -r) cts~
jz -С) 2d '
(17)
1 + £„
Анализ полученных равенств показывает, что в точке Р = Q амплитуды О°, О202 и их сумма О0 + О°2 имеют логарифмическую особенность. Компоненты О°2, О0Х и комбинации О°° - О02 + 2г'О°2 являются непрерывными функциями. Представляя в формулах (15)
G = G0 + (G - G0 )
mn mn v mn mn
G11 - G22 ± 2iG12 =
=h (jj (-1)-1 « -K Ax-r vcos« (x -))
x
1
+2sign(x2 -r) xjj (-1)n
« e 2 sina
(x-)),
л / 2 Г
«i = , K = - > ai
K = Ц« - ^ Гк < ai, к = U.
(15)
выделим там статические члены (16), которые суммируются в явном виде согласно (17), что приводит к следующим равенствам (суммирование по п=1, 2):
о„+Оп=^ (Iе Ке(1„251П ^)+
2 а п 2 а
у е'гп\*2-н %
+ У--- -ЦХ2 -Н +
2 -Г12) 2'
bj +—e «'h-i) x cos«i(x1 -I) ),
1=1 Г2 -/1 i Ki a
i=1
л
i=1
G11 - G22 + 2iG12 =
_ с (Л + ц)
2d
x2 -Л
g -с)
2i 2i
2d
z -С)
2d
(-1)" re
п.. „ir„ix2-m
ne
x2- m
2i (Г -tf)
+
£ ((-1)n
i=1
-I x - m e
AI -а,2 -m
2 2 Г2 -r
i A
'l)cosa( x -£)-
+
^/g"(x2 -m) - 2а
i 1=1 Г22 - Г12
A x -m I I -a x -m
xe nil 2 1 - x - me " 2 1
)sinai (x -#)). (18)
Если в точке р = 2 ряд, фигурирующий в выражении для + С22 в (15), расходился (общий член ряда ведет себя как 1/1), то соответствующий этой сумме функциональный ряд в (18) сходится равномерно и абсолютно: при Р Ф 2 в силу присутствия затухающих экспонент, а при P=Q общий член ряда ведет себя как 1/13. Таким образом, сумма + С22 в точке приложения сосредоточенного источника имеет логарифмическую особенность, соответствующую статическому случаю. Комбинации - С22 + 2гОп в точке P=Q являются непрерывными функциями.
Теперь удобно вместо представлений (14), которые мы рассматриваем в случае совокупности Ь(2) или Ь(3), перейти к интегральным представлениям для комбинаций и1(1) + г и2(1) . Для них находим:
Ц(1) + iU® =
1 x J (p (s) (G® + G22)) + P2 (s) (G® - G22) + 2/G® )) ds,
2 X
U!(1) - iU® =1 x 1 2 2
<J(A (s)(G« -G22) -2/G«) + P2 (s)(Gr + G22)))ds, (19)
Р' ^ = /1 С^О + г /2 Р2 (я) = /1 00 - г /2 (Я).
Здесь Ь - замкнутый контур центральной полости или жесткого включения соответствующей периодической решетки.
Подстановка интегральных представлений (19) в соответствующие граничные условия сводит рассматриваемые краевые задачи к системе двух СИУ 2-го рода в случае полостей и к системе двух СИУ 1-го рода в случае жестких включений. Необходимые дополнительные условия для СИУ 1-го рода вытекают из условий (10), (11).
Численная реализация построенных алгоритмов осуществлялась методом дискретных особенно-
стей (МДО) и механических квадратур [19] ради алгоритмической диверсности. На проведенных тестах оба продемонстрировали равноценную сходимость в зависимости от количества узлов колокации N на контурах (обычно точность вычислений 10-3 достигалась при N=55), причем оба требуют практически одинаковых вычислительных ресурсов (память, быстродействие). Однако для МДО имеется возможность использования программного компонента, обеспечивающего ускоренное выполнение на современном персональном компьютере [9]. Поэтому в качестве основного метода рекомендуем МДО, а при необходимости контролировать точность даже ценой многократного замедления работы, подключать метод механических квадратур.
Пусть неограниченная упругая среда, содержит бесконечную 2^-периодическую систему цилиндров эллиптического сечения. Параметрическое уравнение эллипса, находящегося в основном периоде, задавалось в виде:
£= аа тД щ=- Ъ соаД, о <Д< 2л. (20)
На границе центрального жесткого эллиптического включения и центральной эллиптической полости проводились вычисления максимальных безразмерных напряжений стп, стт (рис. 3) и ст! (рис. 4) соответственно в зависимости от 2С/1, которые получаются делением соответствующих амплитуд контурных напряжений на максимальное напряжение с в падающей волне. Очевидно, а = т1у(>(1+ 2и) в случае набегания на цилиндры продольной волны (1) и ст = т2у() и, если из бесконечности излучается поперечная волна (2). Длина падающей волны 1 = 2л/ у[11 в Р-случае и 1 = 2л/^ в 5У-случае. Выбирались следующие входные параметры:
а / С = о, 5; Ъ / а = о, 5; р / р = 2, о .
Кривые 1, 2 и 3 соответствуют значениям коэффициента Пуассона среды ц = о,1; о, 25 и 0,4 соответственно.
Вычисления показывают, что при дифракции плоской гармонической волны на периодической системе жестких цилиндрических включений и полостей наблюдается резкое увеличение максимальных контурных напряжений возле точек скольжения, что соответствуют значениям 2С/1 = к и 2С/1 = к (к =1, 2, ...). В низкочастотной области с увеличением коэффициента Пуассона ц наблюдается уменьшение максимальных напряжений. Кроме того, в области длинных волн с увеличением отношения периода решетки к длине набегающей волны значения максимальных контурных напряжений увеличиваются в случае периодической решетки, составленной из эллиптических жестких включений, и уменьшаются в случае периодической решетки, составленной из эллиптических полостей.
+
2
L
L
а
б
Рис. 3. Зависимость от максимальных напряжений на контуре эллиптического жесткого включения в
случае Р-волны: а - для; б - для .та х
0 0.5
Рис. 4. Зависимость от максимальных напряжений на контуре эллиптической полости
(в случае Р-волны)
5. Результаты исследований
Схема энергетического метода реализована и доведена до дискретных и компьютерных моделей дифракционных процессов на примере новых задач с периодическими системами неоднородностей. Показано, что в частном случае периодичной системы жестких арматурных прутьев решение базируется на системе двух граничных СИУ 1-го рода, тогда как линейно периодическая система параллельных полых каналов приводит к решению системы СИУ второго рода. Численное моделирование на основе МДО с диверсным контролем точности методом механических квадратур позволило выявить практически важные закономерности в распределении напряжений и наличие их скачков в зависимости от типа неоднородности, волнового диапазона и отношения периода решетки к длине набегающей волны. Показано, что максимальным напряжением на границах включений можно управлять выбором коэффициента Пуассона.
6. Выводы и перспективы дальнейших исследований
Для ряда задач дифракции стационарных плоских волн в упругом теле с неоднородными включениями разных видов (в 2D постановке) на основе энергетического метода моделирования разработаны методы вычисления напряжений, возникающих на границах включений в процессе их взаимодействия с волной. В совокупности с результатами работ этого направления прежних лет, полученные результаты позволяют говорить о создании единой методики решения 2D задач дифракции упругих волн, приводящей к системам граничных уравнений, включающих СИУ 1-го рода. В рассмотренной постановке исследуемые задачи впервые получили решение единообразным методом, что облегчило также программную реализацию и проведение численных экспериментов.
В качестве примера приложения результатов в работе исследовано влияние линейно-периодических
неоднородностей (жесткой арматуры или пустых каналов) на величину максимальных напряжений на их границах.
Литература
1. Гузь, А. Н. Дифрация упругих волн [Текст] / А. Н. Гузь, В. Д. Кубенко, М. А. Черевко. - К.: Наук. думка, 1978. - 307 с.
2. Гладкий, А. В. Численно-аналитические методы исследования волновых процессов [Текст] / А. В. Гладкий, И. В. Сергиенко, В. В. Скопецкий. - К.: Наук. думка, 2001. - 453 с.
3. Зенкевич, О. Конечные элементы и аппроксимация [Текст] / О. Зенкевич, К. Морган. - М.: Мир, 1986. - 318 с.
4. Подильчук, Ю. Н. Применение лучевых методов в задачах распространения и рассеивания волн [Текст] / Ю. Н. Подильчук, Ю. К. Рубцов // Прикл. механика. -1996. - Т. 32, № 12. - С. 3-27.
5. Nazarchuk, Z. T. Singular integral equations in diffraction theory [Text] / Z. T. Nazarchuk // National Academy of Sciences of Ukraine, Karpenko Physico-Mechanical Institute. -1994. - 210 p.
6. Belotsrkovskii, S. M. Method of Discrete Vortises [Text] / S. M. Belotserkovskii, I. K. Lifanov // Boca Raton, CRC Press. - 1993. - 452 p.
7. Gandel, Yu. V. Boundary-value problems for the Helmholz equation and their discrete mathematical models [Text] / Yu. V. Gandel // Journal of Mathematical Sciences. -2010. - Vol. 171, Issue 1. - P. 74-88. doi: 10.1007/s10958-010-0127-3
8. Гандель, Ю. В. Математические модели двумерных задач дифракции: сингулярные интегральные уравнения и численные методы дискретных особенностей [Текст] / Ю. В. Гандель, В. Д. Душкин. - Х.: АВВ МВСУ, 2012. - 544 с.
9. Mishchenko, V. O. Accelerating the computation of the discrete currents method by modification takes into account the architectural features of a modern PCs [Text] / V. O. Mi-shchenko, B. V. Patochkin // Вгсник ХНУ шет В. Н. Ка-разша. - 2015. - № 1056. - С. 129-139.
10. Попов, В. Г. Взаимодействие плоских упругих волн с системами радиальных дефектов [Текст] / В. Г. Попов // Мех. твердого тела. - 1999. - № 4. - С. 118-129.
11. Фильштинский, Л. А. Дифракция упругих волн на трещинах, отверстиях и включениях в изотропной бреде [Текст] / Л. А. Фильштинский // Мех. твердого тела. -1991. - № 4. - С. 119-127.
12. Shibahara, M. Application of the integral equation method to the elastodynamic boundary-value problems [Text] / M. Shibahara, Y. Taniquchi // Bull JSME. -1983. - Vol. 26, Issue 222. - P. 2054-2059. doi: 10.1299/jsme1958.26.2054
13. Nazarenko, A. M. Diffraction of shear waves on cylindrical inclusions and cavities in an elastic half space [Text] / A. M. Nazarenko // Strength of Materials. - 1990. - Vol. 22, Issue 11. - P. 1669-1674. doi: 10.1007/bf00767156
14. Nazarenko, A. M. Dynamical problems of longitudinal shear of half space with nonhomogeneities [Text] / A. M. Nazarenko // Journal of Mathematical Science. - 1994. -Vol. 70, Issue 5. - P. 1985-1990. doi: 10.1007/bf02110825
15. Nazarenko, A. M. A problem of dynamic bending of a sheet with a curvilinear crack [Text] / A. M. Nazarenko, V. I. Ostrik // Soviet Material Science. - 1991. - Vol. 25, Issue 5. - P. 499-502. doi: 10.1007/bf00731941
16. Nazarenko, A. M. Plane problem of diffraction of elastic harmonic waves on periodic curvilinear inserts [Text] / A. M. Nazarenko, O. M. Lozhkin // Materials Science. - 2007. - Vol. 43, Issue 2. - P. 249-255. doi: 10.1007/s11003-007-0028-x
17. Назаренко, А. М. Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде [Текст] / А. М. Назаренко,
Б. Е. Панченко // Проблемы программирования. - Киев, 2010. - № 2-3. - С. 604-610.
18. Назаренко, А. М. Моделирование взаимодействия гармонических P- и SV- волн с периодической системой криволинейных разрезов [Текст] / А. М. Назаренко // Вюник ХНУ iменi В. Н. Каразша. - 2013. - № 1063. - C. 125-135.
19. Назаренко, А. М. Вычислительные методы в задачах дифракции упругих волн на системах неоднородно-стей на базе сингулярных интегральных уравнений [Текст] / А. М. Назаренко. - Сумы: СумДУ, 2015. - 220 с.
20. Новацкий, В. Теория упругости [Текст] / В. Но-вацкий. - М.: Мир, 1975. - 872 с.
References
1. Guz, A. N., Kubenko, V. D., Cherevko, M. A. (1978). Diffraction of elastic waves. Kyiv: Nauk. dumka, 307.
2. Gladkij, A. V., Sergienko I. V., Skopetsky V. V. (2001). Numerical and analytical methods of research of wave processes. Kyiv: Nauk. dumka, 453.
3. Zenkevich, O., Morgan, K. (1986). Finite Elements and Approximation. Moscow: Mir, 318.
4. Podilchuk, Yu. N., Rubtsov, Yu. K. (1996). The use of radiological methods in problems of propagation and scattering of waves. Prikl. Mekhanika, 32 (12), 3-27.
5. Nazarchuk, Z. T. (1994). Singular integral equations in diffraction theory. National Academy of Sciences of Ukraine, KarpenkoPhysico-Mechanical Institute, 210.
6. Belotserkovskii, S. M., Lifanov, I. K. (1993). Method of Discrete Vortises. Boca Raton, CRC Press, 452.
7. Gandel, Yu. V. (2010). Boundary-value problems for the Helmholtz equation and their discrete mathematical models. Journal of Mathematical Sciences, 171 (1), 74-88. doi: 10.1007/ s10958-010-0127-3
8. Gandel, Yu. V., Dushkin, V. D. (2012). Singular integral equations and numerical methods of discrete singularities. Kharkiv: ABB MVSU, 544.
9. Mishchenko, V. O., Patochkin, B. V. (2015). Accelerating the computation of the discrete currents method by modification takes into account the architectural features of a modern PCs. Visnyk KhNU imeni V. N. Karazina, 1056, 129-139.
10. Popov, V. G. (1999). Interaction of plane elastic waves with longitudinal defects. Mekh. tvyordogo tela, 4, 118-129.
11. Filshtinsky, L. A. (1991). The diffraction of elastic waves by cracks, holes and inclusions in isotropic medium. Mekh. tvyordogotela, 4, 119-127.
12. Shibahara, M., Taniguchi, Y. (1983). Application of the Integral Equation Method to the Elastodynamic Boundary-value Problems. Bulletin of JSME, 26 (222), 2054-2059. doi: 10.1299/jsme1958.26.2054
13. Nazarenko, A. M. (1990). Diffraction of shear waves on cylindrical inclusions and cavities in an elastic half space. Strength of Materials, 22 (11), 1669-1674. doi: 10.1007/bf00767156
14. Nazarenko, A. M. (1994). Dynamical problems on longitudinal shear of a half-space with nonhomogeneities. Journal of Mathematical Sciences, 70 (5), 1985-1990. doi: 10.1007/ bf02110825
15. Nazarenko, A. M., Ostrik, V. I. (1990). A problem of dynamic bending of a sheet with a curvilinear crack. Soviet Materials Science, 25 (5), 499-502. doi: 10.1007/bf00731941
16. Nazarenko, O. M., Lozhkin, O. M. (2007). Plane problem of diffraction of elastic harmonic waves on periodic curvilinear inserts. Materials Science, 43 (2), 249-255. doi: 10.1007/s11003-007-0028-x
17. Nazarenko, A. M., Panchenko, B. E. (2010). The scheme of parallel computations in problems of diffraction of shear waves in the hole system in an infinite isotropic medium. Problems of programming. Kyiv, 2-3, 604-610.
18. Nazarenko, A. M. (2013). Modeling the interaction of harmonic P- and SV-waves with a periodic system of curvilinear sections based on the method of singular integral equations. Visnyk KhNU imeni V. N. Karazina, 1063, 125-135.
19. Nazarenko, A. M. (2015). Computational methods in problems of diffraction of elastic waves by inhomogeneities systems based on singular integral equations. Sumy: SumDU, 220.
20. Nowacki, W. (1975). Theory of Elasticity. Moscow: Mir, 872.
Дата надходження рукопису 12.04.2016
Мищенко Виктор Олегович, доктор технических наук, профессор, заместитель заведующего кафедры, Кафедра моделирования систем и технологий, Харьковский национальный университет имени В. Н. Ка-разина, пл. Свободы, 4, г. Харьков, Украина, 61022 E-mail: mischenko@univer.kharkov.ua
Назаренко Александр Максимович, кандидат физико-математичних наук, доцент, кафедра моделирования сложных систем, Сумской государственный университет, ул. Римского-Корсакова, 2, г. Сумы, Украина, 40007
E-mail: aleksandr-nazarenko54@mail.ru
УДК 622.765:542.61:546.571 DOI: 10.15587/2313-8416.2016.69476
ВИДАЛЕННЯ С1НТЕТИЧНИХ БАРВНИК1В З1 СТ1ЧНИХ ВОД © Т. I. Обушенко, Н. М. Толстопалова, I. М. Астрелш
Проаналгзовано сучасний стан проблеми забруднення cmi4Hux вод барвниками, розглянуто iснуючi мето-ди вилучення i деструкцИ барвниюв. В якостi альтернативной запропонована технологiя флотоекстракцИ. До^джено закономiрностi видалення барвника бенгальського рожевого з модельних розчинiв за до-помогою метода флотоекстракцИ. На основi проведених експериментiв встановлено вплив рН середо-вища, тривалостi процесу, мольного спiввiдношення реагуючих речовин та визначено рацiональнi умови проведення флотоекстракцИ
Ключовi слова: барвник, флотоекстракщя, поверхнево-активнi речовини, стiчнi води, бенгальський ро-жевий, гексадецилтриметиламотю бромiд
Current state of wastewater pollution caused by dyestuff was analyzed. Traditional techniques for dyes elimination and destruction were covered. Solvent sublation technology was proposed as an alternative method. The main characteristics of removal of rose bengal dye from wastewater by solvent sublation were investigated. The influence ofpH, time, molar ratio of reactants was studied and the reasonable conditions of the process were determined
Keywords: dye, solvent sublation, surface-active substance, wastewater, rose bengal, hexadecyltrimethylammo-nium bromide
1. Вступ
У зв'язку з обмеженою шльшстю запас1в пр1сно! води на Земл1 (запас доступно! прюно! води на планет становить всього 5-6 тис. м3 на душу населення) проблема охорони пдросфери невпинно загострюеться, хоча для И виршення людство при-кладае чималих зусиль. Рют мют, бурхливий розви-ток промисловосл, штенсиф1кац1я сшьського госпо-дарства, значне розширення площ1 земель, що зро-шуються, покращення культурно-побутових умов 1 ряд шших фактор1в - ва щ чинники ще б1льше ускладнюють проблеми забезпечення прюною водою, незважаючи на значний вклад у захист пдросфери, за рахунок розроблення та впровадження су-часних мехашчних, х1м1чних та бюлопчних техно-логш очищення стоив. Негативний внесок створю-
ють оргашчш сполуки, яш особливо небезпечш як забруднювач1 навколишнього середовища унаслщок комплексного впливу 1 непрогнозованосп насладив. Шюдлива д1я оргашчних речовин, що потрапляють у водойми, посилюеться за рахунок кумулятивного ефекту (прогресуюче збшьшення вмюту шквдливих сполук у кожнш наступн1й ланц1 троф1чного лан-цюга). Серед таких оргашчних забруднювач1в ч1ль-не м1сце займають барвники, як1 широко застосову-ють у р1зних галузях промисловост1 та в побут1. Су-часн1 синтетичш барвники з1 складною х1м1чною структурою не шдлягають б1ох1м1чн1й деструкцН у водних системах. 1х концентрац1я, наприклад, у ст1чних водах текстильних виробництв в залежност1 в1д типу барвника коливаеться в межах в1д 5 до 40 мг/дм3, що багаторазово перевищуе гранично