Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде
Parallel algorithm of numeral solution of theory of elasticity tasks by example of interaction SH-waves with system of apertures with arbitrary normal cross-section is offered.With the help of integral representation of amplitude of shift of reflected wave field the boundary task was reduced to the...
Gespeichert in:
| Datum: | 2026 |
|---|---|
| Hauptverfasser: | , |
| Format: | Artikel |
| Sprache: | Russisch |
| Veröffentlicht: |
PROBLEMS IN PROGRAMMING
2026
|
| Schlagworte: | |
| Online Zugang: | https://pp.isofts.kiev.ua/index.php/ojs1/article/view/959 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Problems in programming |
| Завантажити файл: |  |
Institution
Problems in programming| _version_ | 1866844984113627136 |
|---|---|
| author | Nazarenko, A.M. Panchenko, B.E. |
| author_facet | Nazarenko, A.M. Panchenko, B.E. |
| author_institution_txt_mv | [
{
"author": "A.M. Nazarenko",
"institution": "Sumy State University"
},
{
"author": "B.E. Panchenko",
"institution": "Glushkov Institute of Cybernetics NAS of Ukraine"
}
] |
| author_sort | Nazarenko, A.M. |
| baseUrl_str | https://pp.isofts.kiev.ua/index.php/ojs1/oai |
| collection | OJS |
| datestamp_date | 2026-06-01T06:02:58Z |
| description | Parallel algorithm of numeral solution of theory of elasticity tasks by example of interaction SH-waves with system of apertures with arbitrary normal cross-section is offered.With the help of integral representation of amplitude of shift of reflected wave field the boundary task was reduced to the system of integral, equation which can be numerically solved. Scheme of parallel computing made it possible to research situations with big number of reflecting apertures. New numeric results are resulted.Problems in programming 2010; 2-3: 604-610 |
| first_indexed | 2026-06-02T01:01:54Z |
| format | Article |
| fulltext |
Прикладне програмне забезпечення
© А.М. Назаренко, Б.Е. Панченко, 2010
604 ISSN 1727-4907. Проблеми програмування. 2010. № 2–3. Спеціальний випуск
УДК 004.652, 539.3
СХЕМА ПАРАЛЛЕЛЬНЫХ ВЫЧИСЛЕНИЙ
В ЗАДАЧАХ ДИФРАКЦИИ ВОЛН СДВИГА
НА СИСТЕМЕ ОТВЕРСТИЙ
В БЕСКОНЕЧНОЙ ИЗОТРОПНОЙ СРЕДЕ
А.М. Назаренко, Б.Е. Панченко
Сумский государственный университет, Институт кибернетики имени В.М. Глушкова НАН Украины
Предложен параллельный алгоритм численного решения стационарных задач теории упругости на примере
взаимодействия SH-волн с системой отверстий произвольного поперечного сечения. С помощью интегрального
представления амплитуды перемещения отраженного волнового поля краевая задача сведена к системе
интегральных уравнений, которая решается численно. Схема параллельных вычислений позволила исследовать
ситуации с большим числом отражающих отверстий. Приведены новые численные результаты.
Parallel algorithm of numeral solution of theory of elasticity tasks by example of interaction SH-waves with system of
apertures with arbitrary normal cross-section is offered.With the help of integral representation of amplitude of shift of
reflected wave field the boundary task was reduced to the system of integral, equation which can be numerically solved.
Scheme of parallel computing made it possible to research situations with big number of reflecting apertures. New numeric
results are resulted.
Введение
При оценивании ресурсов конструкций, которые содержат большое число неоднородностей и работают под
действием динамических нагрузок, вызывает интерес анализ взаимодействия волн перемещений и напряжений в
упругой среде с отверстиями или включениями. Поэтому изучение дифракции упругих волн на системах
неоднородностей является актуальным. К тому же такие задачи малоисследованы вследствие необходимости
привлечения больших объемов вычислений и значительных ресурсов цифровой памяти. Учитывая это особое
значение приобретают эффективные параллельные алгоритмы, в основе которых лежат обоснованные
аналитические методы [1].
Среди аналитических методов решения плоских и антиплоских задач теории дифракции следует
отметить метод разложения в ряд по собственным функциям (в случае круговых неоднородностей) [2, 3],
метод конформных отображений (в случае неоднородностей канонической формы) [4, 5], метод
R-функций (в случае неоднородностей сложной геометрической формы) [6, 7], метод интегральных
уравнений (в случае неоднородностей произвольной формы) [8, 9].
Значительный прогресс в теории дифракции достигнут именно благодаря эффективности метода
интегральных уравнений [10–12]. Дополнительные преимущества этого метода заключаются в сокращении
числа пространственных переменных и возможности применения различных эффективных численных методов
решения интегральных уравнений [13, 14]. Кроме того, метод интегральных уравнений обладает большими
возможностями при построении параллельных вычислительных схем. В настоящей работе разрабатывается
параллельный алгоритм численного решения систем интегральных уравнений, возникающих при исследовании
модельной задачи дифракции волн сдвига на системах цилиндрических полостей произвольного поперечного
сечения.
1. Постановка задачи
Рассмотрим в неограниченной изотропной среде m бесконечных вдоль оси OZ полых цилиндров,
поперечные сечения которых ограничены замкнутыми (без общих точек) контурами mjL j ,1, типа
Ляпунова. Пусть L – совокупность указанных контуров и положительное направление выбрано так, что при
движении вдоль L область D остается слева ( рис. 1).
Рис. 1
Прикладне програмне забезпечення
605
Предположим, что источники внешнего поля перемещений W0 размещены в области D. В качестве
такого источника может быть набегающая на цилиндры из бесконечности монохроматическая SH-волна,
нормаль к фронту которой составляет угол с осью OX ( =const),
2 ( cos sin )
0 ,
i x yW e 2
2c
(1)
или гармонический источник интенсивности P, сосредоточенный в точке M0(x0, y0) и порождающий поле
перемещений [8]
(1)
0 0 2
1
( ),
4
P
W H r
i
2 2
0 0( ) ( ) .r x x y y (2)
Здесь c2 – скорость волны сдвига, – частота колебаний, – модуль сдвига, i – мнимая единица
(
2 1i ), )()1( xHn – функция Ханкеля первого рода n-го порядка, зависимость от времени выражается
множеством
tωie .
В результате взаимодействия приходящей волны с отверстиями возникает дифрагированное волновое
поле. Обозначим W1 амплитуду отраженной волны сдвига. Тогда общее поле амплитуд перемещений равно
W= W0+ W1. Неизвестная функция W1 должна удовлетворять однородному уравнению Гельмгольца в области D
с волновым числом 2:
2
1 2 1 0,W W
2 2
2 2
,
x y
(3)
а также условием излучения на бесконечности типа Зоммерфельда [14].
На границе отверстий L нас будут интересовать касательные напряжения ,i t
sz se .i t
nz ne
В случае антиплоской деформации
,s
W
s
,n
W
n
(4)
где s – положительная касательная, n – нормаль в точке i L (рис. 1).
Пусть 0 0 0i – точка L, в которой будем удовлетворять граничным условиям. Так как
L – граница отверстий, то, очевидно,
0 1
0
( ) 0,LW W
n
(5)
где 0n – нормаль к L в точке 0 L .
Таким образом, задача дифракции волны сдвига (1) или (2) на системе отверстий в неограниченной
изотропной среде сводится к решению краевой задачи (3), (5) при выполнении дополнительных условий
излучения на бесконечности.
2. Метод решения
Следуя [8], запишем функцию W1(x,y), характеризующую рассеянную отверстиями волну перемещений в
области D, следующим образом:
1( , ) ( ) ( , , , ) ,
L
W x y f s G x y ds (6)
(1)
0 2
1
( ),
4
G H r
i
,r z ,z x iy .i L
Здесь L – совокупность контуров mjL j ,1, (рис. 1); f(s) – неизвестная функция, удовлетворяющая на
L условию Гельдера.
Интегральное представление (6) удовлетворяет уравнению Гельмгольца (3) в области D и обеспечивает
выполнение условий излучения на бесконечности. Остается выполнить граничное условие (5). Для
осуществления предельного перехода в (5) при 0z L частные производные
0
W
s
и
0
W
n
будем понимать
следующим образом:
0
0 0
0
,
i i
L
z
W W W
e e
s z z
0
00
,0i i
L
z
W W W
i e e
n z z
0 0
0
,
i d
e
ds
0 0 0 .i L (7)
Прикладне програмне забезпечення
606
Воспользуемся также известными соотношениями [8]:
(1) (1)
0 1( ) ( ),
2
iH r e H r
z
(1) (1)
0 1( ) ( ),
2
iH r e H r
z
(1)
1 1
2
( ) ( ),H r H r
i r
,iz re (8)
где 1( )H x непрерывная функция в точке x=0.
Привлечение формулы Сохоцкого – Племеля [14] для вычисления предельных значений интегралов типа
Коши, возникающих при удовлетворении граничного условия (5) с учетом соотношений (6) – (8), приводит
к искомому интегральному уравнению относительно неизвестной функции f(s):
0 0 0
1
( ) ( ) ( , ) ( ),
2
n
L
f s f s E s s ds K s 1, 2,n (9)
(1)
0 2 1 2 0 0 0( , ) ( )sin( ),E s s H r 0
0 0 ,
i
r e
1 0 2 0 0 0( ) ( )sin( ),K s i W s
(1)
2 0 2 1 2 10 10 0( ) ( )sin( ),
P
K s H r
10
0 0 10 ,
i
z r e 0 0 0.z x iy
Здесь функции 1 0( )K s и 2 0( )K s отвечают случаям (1) и (2), соответственно.
Представим ядро 0( , )E s s , учитывая (8), в виде:
0
0 2 1 2 0 0 0
0
2
( , ) Im ( )sin( ).
i
e
E s s H r
i
10)
Нетрудно убедиться [14], что функция 0( , )E s s непрерывна на L. Следовательно, интегральное
уравнение (9) является уравнением Фредгольма второго рода, которое как известно, разрешимо и имеет
единственное решение в классе функций, непрерывных по Гельдеру.
3. Дискретизация задачи
Представим неизвестную плотность f(s) интегрального уравнения (9) как совокупность функций fj(sj),
определенных на контурах mjL j ,1, . Тогда (9) превращается в систему интегральных уравнений
Фредгольма 2-го рода:
0 0 0
1
1
( ) ( ) ( , ) ( ),
2
m
k k j j j k j n k
j L
f s f s E s s ds K s 1, ,k m 1, 2.n (11)
Здесь дуговые координаты js и
k
s относятся к точкам j j j ji L и 0 0 0k k k ki L
соответственно.
Численная реализация интегральных уравнений (11) проводилась методом механических квадратур [14].
Вводилась параметризация контура Lj с помощью соотношений
( ),j j 0 0 0( ),j j 00 , 2 , (12)
причем (0) (2 ).j j Интегральное уравнение, соответствующее контуру Lk, удовлетворялось в узлах вида
(2 1) /l kl n ( 1, )kl n и сводилось к системе линейных алгебраических уравнений относительно
значений функции fj( ) в узлах вида (2 1) /p jp n ( 1, )jp n , где nj – число точек разбиения контура Lj.
Внеинтегральные значения ( )k lf выражались с помощью интерполяционных полиномов Лагранжа через
искомые значения ( ).k pf Для них получены следующие выражения [8]:
1
1
( ) ( 1) ( )cos ,
2
k
k
n
l pn p
k l k p
pk
f f ec
n
(13)
если nk – нечетно, и
1
1
( ) ( 1) ( ) ,
2
k
k
n
l pn p
k l k p
pk
f f ctg
n
(14)
если nk – четно.
Таким образом, при численной реализации системы интегральных уравнений (8) задача сводится к
решению системы линейных алгебраических уравнений с N=n1+n2+…+nm неизвестными.
Прикладне програмне забезпечення
607
4. Схема вычислений
Пусть система интегральных уравнений сведена к системе линейных алгебраических уравнений, все
элементы матрицы которой являются результатом дискретизации контуров. Очевидно, размер матрицы
пропорционален числу отверстий. Для исследования описанного метода при большом числе отверстий, а также
для получения высокоточных результатов с погрешностью вычислений до 10
-6
и проверки сходимости решений
при большом числе точек коллокации потребуются существенные вычислительные ресурсы. Применим
распараллеливание алгоритма. Из (11) следует, что каждый элемент матрицы определяется координатами узлов
дискретизации.
Как показано на рис. 2, данный метод в вычислительном смысле сводится к обходу каждого контура по
точкам коллокации внеинтегральной переменной 0k и одновременному обходу каждого контура по анало-
гичным либо иным узлам переменной интегрирования .k
Рис. 2
Важной особенностью алгоритма такого обхода является то, что результирующая матрица строится в
соответствии с отношением „многие ко многим” и формально является результатом Декартова произведения
этих множеств [15].
Это означает, что все элементы матрицы формально независимы один от другого, что не строго
доказывает возможность применения параллельного вычисления.
Таким образом, переменная 0k формирует строки матрицы СЛАУ, а переменная k – ее столбцы.
Диагональные элементы матрицы соответствуют коэффициентам системы, вычисленным в общих узлах для
0k и k отверстий. Иные коэффициенты вычисляются так, что значения 0k принадлежат множеству точек
коллокации с одних контуров, а значения переменных интегрирования k – с других.
Распараллеливание вычислительного алгоритма может быть осуществлено для основных четырех
процессов – варьирование группами исходных данных, вычисление коэффициентов системы линейных
уравнений, численное решение СЛАУ и вычисление искомых характеристик дифрагированного поля по
полученным после решения СЛАУ функциям. Наиболее емким с точки зрения временных затрат и ресурсов
памяти являются второй и третий процессы. Для параллельного формирования матрицы СЛАУ
вычислительный алгоритм сведется к двум основным процедурам. Первая формирует массив параметрических
координат контуров всех отверстий, который является своеобразной базой данных для формирования клонов.
При этом для обеспечения равномерной загрузки хоста кластерной системы исходный массив дробится
пропорционально числу используемых процессоров.
Так как для данной методики решения краевой задачи основная операция при вычислении каждого
элемента матрицы – это определение разностного аргумента цилиндрических функций Ханкеля, заданного на
множестве значений параметрических координат отверстий, а также вычисление самих этих функций и
коэффициентов при них, то на следующем шаге на каждом клоне хоста запускаются цикловые процедуры
определения указанных коэффициентов. При этом синхронизация каждого процесса не требуется, так как
итоговая матрица собирается по факту завершения последнего.
Вычислительный процесс решения СЛАУ также распараллеливается согласно [16]. Параллельное
вычисление итоговых искомых характеристик осуществляется путем подстановки массивов значений
неизвестных функций ( )k pf в представление (6) аналогично процедурам формирования матрицы СЛАУ.
5. Численные результаты
С целью исследования сходимости построенного алгоритма рассмотрим случай нормального падения
(ψ=π/2) волны сдвига (1) на систему эллиптических отверстий, расположенных вдоль одной линии на
одинаковом расстоянии d один от другого (рис. 3).
Прикладне програмне забезпечення
608
Рис. 3
Предполагаем, что основной эллипс L0 задается уравнением
( ) sin ,a ( ) cos ,b 0 2 . (15)
В этом случае рассматриваемая дифракционная задача обладает свойством симметрии, что позволяет
осуществлять первичное самотестирование получаемых результатов.
В ходе численной реализации вычислялись безразмерные контурные напряжения /s . Точность
вычислений проверялась путем сравнения результатов при различных значениях N. Проводилось также
сравнение полученных результатов с результатами, приведенными в [5, 8] для случая одиночного отверстия.
Применение метода параллельных вычислений, проведенного на кластере «Инпарком-256», позволило
исследовать сходимость до 10
-6
порядка. Обнаружено, что сходимость решения интегрального уравнения
практически не зависит от числа отражателей. При наличии 100 отверстий общее время решения системы
линейных уравнений для матрицы размерностью порядка 10
6
(для одного варианта) с использованием
256 процессоров составляло около 2 часов. Решение аналогичной задачи процедурным методом для 10 отражателей
(использование типовых вычислительных средств с точностью до 10
–3
– 500 точек коллокации) осуществлялось
несколько суток. При симметричном нагружении наблюдается эффект насыщения решетки уже при
9–10 отверстиях (при дальнейшем увеличении числа отверстий контурные напряжения практически не меняются).
В работе проводились вычисления контурных напряжений βσ вдоль контуров центрального L0
и крайнего Lk отверстий (рис. 3) в случае решетки, состоящей из нечетного числа отверстий (p=k). Отсчет угла
ведется от нуля (теневая точка) до (лобовая точка) для центрального отверстия (учитывается симметрия в
случае нормального распределения волны сдвига (1)) и от 0 до 2 для крайнего отверстия (в силу симметрии
распределения напряжений на контурах Lk и L-k одинаковы).
Анализ показывает, что существует принципиальное отличие в распределении βσ при набегании SH-
волны (1) на решетку, состоящую из вытянутых вдоль оси Ox эллипсов (b/a<1), от случая решетки, вытянутой
навстречу волне (b/a>1). Если в первом случае (b/a<1) распределение βσ является моноэкстремальным [5, 8]
(максимальное значение – в точках соскальзывания = /2 и =3 /2), то во втором случае, начиная с некоторого
значения b/a>1, появляется несколько локальных максимумов. Поэтому в работе рассматривается случай
эллипсов, вытянутых навстречу набегающей волне, при этом выбирается значение b/a=2.5.
На рис. 4 показаны распределения вдоль контура центрального отверстия L0 в случае решетки, состоящей
из трех (кривая 1), пяти (кривая 2) и семи (кривая 3) эллипсов, при значениях безразмерного волнового числа γ2a,
равных 1.7 (а) и 2.5 (б). Если в теневой ( = 0) и лобовой ( = ) точках =0, то в зоне соскальзывания
с увеличением γ2a число локальных максимумов также увеличивается, причем растет и максимальное
значение .
Прикладне програмне забезпечення
609
а б
Рис. 4. γ2a = 1.7 – а; γ2a = 2.5 – б
На рис. 5 показано аналогичное распределение вдоль контура крайнего отверстия. Здесь
наблюдаются локальные минимумы в теневой ( = 0) и лобовой ( = ) точках. Число локальных максимумов
в зоне соскальзывания с увеличением γ2a также увеличивается. Максимальное значение растет
с увеличением γ2a.
а б
Рис. 5. γ2a=1.7 – а; γ2a=2.5 – б
Рис. 6, 7 иллюстрируют отклонение максимальных значений βσ на контуре центрального (рис. 6)
и крайнего (рис. 7) отверстий от максимального значения (0) на контуре одиночного включения в зависимости
от параметра решетки d (рис. 3) при b/a=2.5, γ2a=5.0 в случае симметричного нагружения SH-волной (1). Для
этого по аналогии с [12] на данных рисунках представлены величины * (0)max max . Здесь
рассмотрены случаи трех и семи (рис. 6 а, б) отверстий в решетке. Как и следовало ожидать, растягивание
решетки приводит к уменьшению значений
*
βσ и стабилизации пульсаций этих значений на нулевом уровне
при достаточно больших d. Причем с увеличением волнового числа γ2a стабилизация пульсаций разности *
достигается при бóльших значениях параметра решетки d.
а б
Рис. 6. m=3 – а; m=7 – б
Прикладне програмне забезпечення
610
Известно [2], что для установившихся колебаний при сколь угодно больших значениях расстояния
между объектами, вблизи объектов все же наблюдается незначительное влияние решетки. Графики, показаны
на рис. 6, 7, иллюстрируют слабые пульсации значений * вблизи нуля. Этим подтверждается достоверность
используемого алгоритма.
а б
Рис. 7. m=3 – а; m=7 – б
Выводы
Для задачи дифракции сдвиговых волн на системе отверстий некруговой формы параллельные
алгоритмы позволяют значительно сократить время вычислений и более детально проанализировать
характеристики волнового поля. Это очень важно, так как получение точных величин резонансных максимумов
контурных напряжений, а также точных координат дислокации резонансных максимумов, позволит избежать
разрушений конструкций, работающих в условиях динамических нагрузок. Все это требует использования
увеличенного числа вариантов исходных данных. Сочетание же метода интегральных уравнений,
позволяющего на единицу снизить размерность задачи, а также значительная экономия времени вычислений за
счет распараллеливания вычислительных процедур, приводит к существенному увеличению эффективности
предложенного алгоритма.
1. Вертгайм И.И., Терпигуров В.Н. Параллельные технологии вычислений в механике сплошных сред и МДТТ.: Учебное пособие. –
Пермь, 2007. – 84 с.
2. Гузь А.Н., Головач В.Т. О решении основных граничных задач теории установившихся колебаний для бесконечной плоскости с
круговыми отверстиями // Изв. АН СССР. Механика твердого тела. – 1968. – № 2. – С. 58 – 64.
3. Pao Y.H., Mow C.C. The diffraction of elastic waves and dynamic stress concentration. – New York: Crane Rassak Co, 1973. – 689 p.
4. Datta S.K. Diffraction of SH-waves by an elliptic cylinder // Int. J. Solid Struct. – 1974. – 10, № 1. – P. 123 – 133.
5. Гузь А.Н., Кубенко В.Д., Черевко М.А. Дифракция упругих волн. – Киев: Наук. думка, 1978. – 307 с.
6. Рвачев В.Л. Теория R-функций. – Киев: Наук. думка, 1982. – 556 с.
7. Гуляев Ю.В., Кравченко В.Ф., Рвачев В.Л., Сизова Н.Д. Исследование дифракции упругих волн на пластинах, ослабленных двумя
отверстиями произвольной формы // ДАН. Математическая физика. – 1996. – 349, № 2. – С. 175 – 179.
8. Назаренко А.М., Панченко Б.Е. Дифракция волн сдвига на цилиндрических неоднородностях произвольного поперечного сечения //
Динамика и прочность машин. Респ. межвед. научно-техн. сб. – 1991. – Вып. 52. – С. 38 – 45.
9. Фильштинский Л.А. Дифракция упругих волн на трещинах, отверстиях, включениях в изотропной среде // Изв. АН СССР. Механика
твердого тела. – 1991. – № 4. – С. 119 – 127.
10. Назаренко А.М. Дифракция волн сдвига на цилиндрических включениях и полостях в упругом полупространстве // Проблемы
прочности. – 1990. – № 11. – С. 90 – 94.
11. Feng Y.D., Wang Y.S., Zhang Z.M. Time domain BEM analysis of cylinder embedded in soil with frictional slip at the interface //Soil Dyn. Earthquake
Eng. – 2003. – 23, №4. – P. 303-311.
12. Кюркчан А.Г., Скородумова Е.А. Решение трехмерной задачи дифракции волн на группе объектов //Акустический журнал. – 2007. – 53,
№ 1. – С. 5 – 14.
13. Габдулхаев Б.Г. К численному решению интегральных уравнений методом механических квадратур // Изв. вузов. Математика. – 1972.
– № 2. – С. 23 – 39.
14. Панасюк В.В., Саврук М.П., Назарчук З.Т. Метод сингулярных интегральных уравнений в двумерных задачах дифракции. – Киев:
Наук. думка, 1984. – 344 с.
15. Курош А.Г. Общая алгебра. – М.: 1979. – 159 с.
16. Химич А.М., Полянко В.В. Ефективность двумерных блочно-циклических паралельних алгоритмов // Проблеми програмування. – 2008.
– № 3. – С. 145 – 149.
|
| id | pp_isofts_kiev_ua-article-959 |
| institution | Problems in programming |
| keywords_txt_mv | keywords |
| language | Russian |
| last_indexed | 2026-06-02T01:01:54Z |
| publishDate | 2026 |
| publisher | PROBLEMS IN PROGRAMMING |
| record_format | ojs |
| resource_txt_mv | ppisoftskievua/57/d67d822833001e4ea730d8f45f0eca57.pdf |
| spelling | pp_isofts_kiev_ua-article-9592026-06-01T06:02:58Z Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде Nazarenko, A.M. Panchenko, B.E. UDC 004.652, 539.3 УДК 004.652, 539.3 Parallel algorithm of numeral solution of theory of elasticity tasks by example of interaction SH-waves with system of apertures with arbitrary normal cross-section is offered.With the help of integral representation of amplitude of shift of reflected wave field the boundary task was reduced to the system of integral, equation which can be numerically solved. Scheme of parallel computing made it possible to research situations with big number of reflecting apertures. New numeric results are resulted.Problems in programming 2010; 2-3: 604-610 Предложен параллельный алгоритм численного решения стационарных задач теории упругости на примере взаимодействия SH-волн с системой отверстий произвольного поперечного сечения. С помощью интегрального представления амплитуды перемещения отраженного волнового поля краевая задача сведена к системе интегральных уравнений, которая решается численно. Схема параллельных вычислений позволила исследовать ситуации с большим числом отражающих отверстий. Приведены новые численные результаты.Problems in programming 2010; 2-3: 604-610 PROBLEMS IN PROGRAMMING ПРОБЛЕМЫ ПРОГРАММИРОВАНИЯ ПРОБЛЕМИ ПРОГРАМУВАННЯ 2026-06-01 Article Article application/pdf https://pp.isofts.kiev.ua/index.php/ojs1/article/view/959 PROBLEMS IN PROGRAMMING; No 2-3 (2010); 604-610 ПРОБЛЕМЫ ПРОГРАММИРОВАНИЯ; No 2-3 (2010); 604-610 ПРОБЛЕМИ ПРОГРАМУВАННЯ; No 2-3 (2010); 604-610 1727-4907 ru https://pp.isofts.kiev.ua/index.php/ojs1/article/view/959/1027 Copyright (c) 2026 PROBLEMS IN PROGRAMMING |
| spellingShingle | UDC 004.652 539.3 Nazarenko, A.M. Panchenko, B.E. Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| title | Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| title_alt | Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| title_full | Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| title_fullStr | Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| title_full_unstemmed | Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| title_short | Схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| title_sort | схема параллельных вычислений в задачах дифракции волн сдвига на системе отверстий в бесконечной изотропной среде |
| topic | UDC 004.652 539.3 |
| topic_facet | UDC 004.652 539.3 УДК 004.652 539.3 |
| url | https://pp.isofts.kiev.ua/index.php/ojs1/article/view/959 |
| work_keys_str_mv | AT nazarenkoam shemaparallelʹnyhvyčislenijvzadačahdifrakciivolnsdviganasistemeotverstijvbeskonečnojizotropnojsrede AT panchenkobe shemaparallelʹnyhvyčislenijvzadačahdifrakciivolnsdviganasistemeotverstijvbeskonečnojizotropnojsrede |