Моделирование крутильных колебаний вязкоупругого круглого стержня, вращающегося с постоянной угловой скоростью

УДК 539.3

Моделирование крутильных колебаний вязкоупругого круглого стержня, вращающегося с постоянной угловой скоростью

© Х.Х. Худойназаров1, А. Абдирашидов2, Ш.М. Буркутбоев1

:ТГТУ, Ташкент, 100095 Узбекистан 2СамГУ, Самарканд, 140100 Узбекистан

Выведены общие и приближенные уравнения крутильных колебаний круглого вязко-упругого стержня, вращающегося с постоянной угловой скоростью вокруг оси симметрии. Разработан алгоритм, позволяющий определить напряженно-деформированное состояние этого стержня. На основе полученных приближенных уравнений колебаний численно решена задача о его крутильных колебаниях. Проведен сопоставительный анализ результатов, полученных для экспоненциального и слабосингулярного ядер оператора вязкоупругости. Даны оценки влияния вращения на колебания стержня.

Ключевые слова: стержень, крутильные колебания, перемещение, напряжения, вращение, угловая скорость, вязкоупругость.

Введение. Изучению задач о колебаниях вращающихся стержней и оболочек с учетом их вязкоупругих свойств посвящено небольшое количество работ. Динамическое поведение балки, вращающейся с постоянной угловой скоростью вокруг своей продольной оси, исследовано в работе [1] для свободной, защемленной, свободно опертой балки, а также для других видов закрепления ее торцов. Задача о колебаниях стержня, вращающегося вокруг своей оси, с учетом геометрической нелинейности, обусловленной применением неподвижных в продольном направлении шарнирных опор, решена в [2]. Свободные и вынужденные колебания вращающейся слоистой цилиндрической оболочки исследованы в работах [3-5].

Численное исследование влияния реологических параметров на характер колебаний наследственно-деформируемых систем имеет большое значение. В работе [6] показано, что дифференциальная зависимость между усилиями и деформациями, имеющая место при решении динамических задач вязкоупругости, приводит к определенной погрешности, особенно в начальный момент времени. Исследования нелинейной устойчивости и колебания вязкоупругой балки, движущейся в осевом направлении, выполнены на основе моделей Кельвина и Максвелла в работе [7]. Проведено сравнение численных решений полных нелинейных и линеаризованных уравнений колеба-

ния. Работа [8] посвящена решению спектральной задачи, связанной с описанием малых поперечных колебаний однородного вязкоупругого стержня. В работе [9] предложена методика расчета эффективных вязкоупругих характеристик композиционных материалов при установившихся циклических колебаниях, которая основана на методе асимптотического осреднения периодических структур и конечно-элементном решении локальных задач вязкоупругости на ячейке периодичности композитов.

Построение основных соотношений стержневой теории заключается в приведении трехмерной по пространственным координатам задачи к двумерной или одномерной. Исследователи при этом стремятся вывести уточненные дифференциальные уравнения колебаний, учитывающие те или иные физические, механические или геометрические факторы [10, 11].

Постановка задачи. В цилиндрической системе координат (г, 9, 2) рассмотрим однородный и изотропный круглый вязкоупругий стержень радиусом Го, вращающийся вокруг оси симметрии с постоянной угловой скоростью О. Считается, что механические свойства стержня описываются в рамках линейной теории вязкоупругости, в которой связь между напряжениями и деформациями задана соотношениями Больцмана — Вольтерра.

Задача о крутильных колебаниях круглого цилиндрического стержня может быть рассмотрена отдельно от задачи продольно-радиальных его колебаний [10, 11]. В силу осесимметричности задачи компоненты полей напряжений и перемещений не зависят от угловой координаты 9. При этом отличны от нуля только перемещение и9 (г, г, г) и напряжения а^ (г , г, г), а ^ (г, г, г) [12]:

где у — компонента векторного потенциала поперечных волн;

функция, к которой применяется оператор); ц — модуль сдвига; Я(г) — ядро оператора вязкоупругости.

Предполагается, что оператор М обратимый.

При учете центробежной силы, вызванной вращением стержня, уравнение его движения относительно ненулевых компонент напряжения имеет вид [5]

do

rQ

дг dz

до zQ. + 2 ürQ

= р

dt2

-Q2Uq

где р — плотность материала стержня.

Предполагается, что крутильные колебания стержня вызваны напряжением (г, I) на его поверхности, т. е. граничное условие

задачи при г = г имеет вид

°rQ (^t) = frQ (^t)-

Начальные условия считаются нулевыми.

Вывод уравнений колебаний. Подставив выражения перемещения и напряжений (1) в уравнение движения (2) получим следующее интегродифференциальное уравнение:

,/d2y 1 dy d2y^

м —2-+—-+—2-

^ dr2 r dr dz2 )

= P

fd 2у „2

HF

Л

У

Искомую функцию в формуле (4) представим в виде [13] y(t, z ) = Ц У(0) (p, k )ept+kzdpdk,

h i

1 '2

где у(t, z) — непрерывная функция в области t > 0 и z > 0; у(0) (p, k) — образ при ее двухмерном преобразовании (Re p > 0, -ж arg p < я; Rek > 0, -ж arg k < я); 'i и '2 — неограниченные открытые контуры в комплексной плоскости.

Подставляя (5) в уравнение (4), получим обыкновенное дифференциальное уравнение Бесселя. Общее решение этого уравнения, ограниченное при r = 0, имеет вид

у (0)( p, k) = BI0( yr).

Здесь 10 — модифицированная функция Бесселя нулевого порядка; В — постоянная интегрирования;

y2 =рМ0-1 (p2-Q2)- k 2,

где М0 =ц[1 -^(p)].

Представив перемещение и9 в виде (5), выразив полученное

таким образом перемещение и90) через общее решение (6) и разложив в нем функцию Бесселя в ряд по степеням радиальной координаты, будем иметь

ч2и+1

и 90)

п = 0

|!(и +1) !

где и(0) = - -у2В — главная часть перемещения и90). Введем операторы Xп (£) по формулам

Хп = Л у2пер+кгёрСк, п = 0, 1, 2,....

1 '2

Тогда, представив напряжение аг9 также в виде (5), из граничного условия (3) получим

\2 п + 2

21

(г/ 2 )2

п = 0

I! (п + 2) !

х п+1и = М-1

1т9

где М 1 — оператор, обратный оператору М.

Если учесть, что параметр у имеет вид (7), то нетрудно увидеть, что операторы Хп в переменных (2, г) имеют вид [10]

Хп =

рМ

-1

( а2

^-О2

дл2

дг2

п = 1, 2, ....

С учетом вида введенных интегродифференциальных операторов

Хп (п = 0,1,2,...) полученное уравнение (10) является

интегродифференциальным уравнением бесконечно высокого порядка относительно главной части крутильного перемещения. Поэтому назовем его общим уравнением крутильных колебаний вязкоупругого круглого стержня, вращающегося с постоянной угловой скоростью вокруг своей оси.

Отметим, что уравнение (10) содержит производные любого порядка по координате 2 и по времени г. Кроме того, правая часть этого уравнения учитывает внешнее усилие, приложенное к

поверхности стержня, и отражает зависимость от оператора вязко-вязкоупругости, если он известен для конкретной среды.

Формулы для перемещения и напряжений. Применив преобразование (5) к выражению преобразованного перемещения (8), с учетом (9) получим

Uе

и=0 п!(п +1)!

Аналогично для ненулевых компонент напряжений будем иметь

» (rl2)2п+2

а* (r, z, t) = 2M £ Л, n+U;

n=0 n !(n + 2)! , ч

2 ,

(rl2)2n+1 Я а ze (r, z, t) = 2M £ Xn fu.

n=0 n !(n + 1)! dz

Формулы (12) и (13) позволяют с заданной точностью по радиальной координате и времени определить перемещение и напряжения в точках произвольного сечения стержня через решение уравнения (10).

Приближенные уравнения колебаний. Как уже отмечалось, уравнение (10) имеет очень высокий порядок и, естественно, в таком виде не может быть применено в инженерных расчетах. Отсюда следует, что необходимо ограничить количество его членов, т. е. ограничиться нулевым (п = 0), первым (п = 1), вторым (п = 2) и последующими приближениями.

Предполагая, что для случая вращающегося стержня также выполнены условия относительно области применимости «усеченного» таким образом уравнения, аналогичные полученным в работах [10, 11], можно найти различные приближенные уравнения из уравнения (10). Так, ограничиваясь в нем нулевым приближением, имеем уравнение второго порядка

..2 (

Г0_

4

Я2

Л

at2

Я

,2 ^

v^ / у

az2

pM"1 -J-Q2 U = M"1 [/re],

которое при отсутствии вращения в точности совпадает с уравнением, полученным профессором И.Г. Филипповым [10].

Первое приближение уравнения (10) дает уравнение четвертого порядка

(

48

р2 м -2_ 2р2М -202 Н дУ4 Н дУ2

р2м-2о4 -

- 2рМ-1

4

„2 Г

дУ 2 дг 2 Г л2

2РМ-1"2 £

рм

-1

дг-о 2

дУ2

дг2

и = м-1 [у^]

которое является уточненным по отношению к (14) и учитывает в своей структуре, кроме угловой скорости вращения, деформацию поперечного сдвига и инерцию вращения.

При ограничении нулевым, первым или другими приближениями уравнения колебаний (10) следует также ограничиться соответствующим приближением в формулах для перемещения (12) и напряжений (13): • при п = 0

Цв= ги;

а гв =-

-2 Г д2и _ д2и

Л

- м-

дг2

а гв= гм ди;

дг

при п =1

ив = ги+—Ш; агв= м 16

Г г 2

г\2 >

хи+—х2и

4 48

а гв = м

Г ди г\ диЛ

г-+—X-

дг 8 дг у

Если найдены решения уравнения (14) или (15), то по ним легко определить напряженно-деформированное состояние стержня посредством формул (16) или (17).

Краевая задача. Исследуем влияние вращения и вязкоупругих свойств материала на напряженно-деформированное состояние стержня на основе выведенных уравнений колебания. Допустим, что поверхность стержня свободна от внешних нагрузок. Тогда, полагая /гв = 0 в правой части уравнения (14), учитывая вид оператора вязкоупругости и вводя безразмерные переменные по формулам

У = У

* 10. ъ'

г = г*г0; 0 = 0*Ъ; и = и

'о

где Ъ = -у/р, / р — скорость распространения поперечных волн в материале стержня, получим

г

d2U d2U í чd2U

j R

dt2 dz2 0 ' dz

Для дальнейших расчетов для оператора вязкоупругости примем слабосингулярное ядро А.Р. Ржаницына [14]:

R (t) = At a-1e-pt,

где A, а, p — реологические параметры, удовлетворяющие условиям A > 0, 0 <а< 1, p > 0.

Рассмотрим круглый стержень длиной l, подвергнутый воздействию синусоидальной внешней нагрузки на одном конце. Будем считать, что другой его конец закреплен. Тогда перемещения U (z, t) точек стержня удовлетворяют следующим граничным и начальным условиям:

U(0, t) = A sin — при z = 0; U(l, t) = 0 при z = l; (19) t1

U

dt

где ^ — время действия нагрузки.

Численный алгоритм решения. Для решения краевой задачи выразим интегродифференциальное уравнение (18) с помощью неявной схемы Кранка — Никольсона [15], при этом на каждом временном шаге интеграл в уравнении (18) вычисляется по квадратурной формуле прямоугольников

U++1 + трГ + nU-1 = f, 0 < i < N +1; n > 0.

тП+1

Tn+1

Здесь

m =

r ±+^ h2 + A2 2

; n = 1;

fi = -2A

1 m-1

h

n-1

h2 ( -2Un-1 )+¿ (Un-1 -2Un-1 +Un-1 ) +

[

-Q2Uf-1 - A—у У

o i 1^2 ¿—i

,-P(n-k )h

A2 k=

k=0

где ^ — шаг по времени; А — шаг по координате. Граничные условия (19) имеют вид

К X

иои = А МП ^ при г = 0; Щ+1 = 0 при г = I,

ч

а начальные условия (20) —

ТТ1 - ТТо

и0 = 1 1 = 0 при X = 0.

к

Система алгебраических уравнений (21)-(23) решена с помощью метода прогонки.

Результаты численных расчетов. Для расчетов приняты следующие значения безразмерных параметров: г0 = 1; I = 10. Полученные результаты представлены на рис. 1-6.

и/А 1,0 0,8 0,6 0,4 0,2 О -0,2

0 2 4 6 8 10 X

Рис. 1. Зависимости перемещения от времени в фиксированных (г = 2, 5, 8) сечениях вязкоупругого стержня с экспоненциальным ядром (А = 0,1; а = 1) при различных значениях параметра Р :

1 — упругий стержень; 2 — Р = 1; 3 — Р = 0,5; 4 — Р = 0,1

11/А

гг = 2 -г=8

/ * Г

// Л / Л / 'г ПС "

/ \ ? Л у Л ' 11

/ \ / // к

1 г^ ч^ г, ; ___

?г = 2

/ Л / ' V / 1 1/ 2

Л> / \д /

// \\ / V < / /

// \ У 4  \ Ч //У Г\ ' / V '• ' \

У/- / V» >

1,0 0,8 0,6 0,4 0,2 О -0,2

0 2 4 6 8 10 Г

Рис. 2. Зависимости перемещения от времени в фиксированных (г = 2, 5, 8) сечениях вязкоупругого стержня со слабосингулярным ядром (А = 0,1; а = 1) при различных значениях параметра Р :

1 — упругий стержень; 2 — Р = 1; 3 — Р = 0,5; 4 — Р = 0,1

U/A 1,0 0,8 0,6 0,4 0,2 0 -0,2

yZ = 2 = 5 1z = 8

/ /

5- * —к.' vk / -¿A / ji

4~ Jr/ \w / > \v / / 7 ь \ / У \\ - Vv V

I? ///' V //; YW/ / / у, ч7 ''/ > / ^' ^ \\ XV \ \ 4 ^

Sr. \w, Vi-JA. rtf"" \ ч

0 2 4 6 8 10 t

Рис. 3. Зависимости перемещения от времени в фиксированных (z = 2, 5, 8) сечениях вязкоупругого стержня со слабосингулярным ядром (A = 0,1; р = 0,1) для различных значений параметра а:

1 — упругий стержень; 2 — а = 0,5; 3 — а = 0,25; 4 — а = 0,1

«WW 0,! 0,4 0

-0,4 -0,! -1,2

6

Рис. 4. Изменения напряжения с течением времени в фиксированных (z = 2, 5, 8) сечениях стержней:

1 — упругого

U/A 1,4 1,2 1,0 0,8 0,6 0,4 0,2 О

-0,2

0 2 4 6 8 10 /

Рис. 5. Зависимости перемещения точек различных сечений упругого стержня от времени при различной угловой скорости:

1 — Q = 0; 2 — Q = 0,2; 3 — Q = 0,3

** ' ч

z = 2 / z= 5 / wz=8

/ • / / ' V/

Г \\ \ - // >

> \\ \ \ 4 У \A \ 4

V3 / \\ \ \ v / \ 4 \ Л  ' -

V / \ \ , /._-

V J* i" \

и/А 1,0 0,8 0,6 0,4 0,2 0 -0,2

0 2 4 6 8 10 Г

Рис. 6. Зависимости перемещения точек различных сечений вязкоупругого стержня (А = 0,1; а = 0,25; в = 0,5) от времени при различной

угловой скорости:

1 — О = 0; 2 — О = 0,2; 3 — О = 0,3

Анализ полученных результатов. Зависимости перемещения и вязкоупругого стержня от времени (см. рис. 1) для различных значений параметра Р экспоненциального ядра (а = 1) показывают снижение амплитуды перемещения при уменьшении р. В сечении г = = 8 максимальное значение перемещения снижается на 33 %.

Кривые изменения перемещения во времени при различных значениях параметров р и а в вязкоупругом стержне со слабосингулярным ядром (см. рис. 2 и 3) показывают, что и в этом случае уменьшение значений параметра р приводит к уменьшению амплитуды перемещения. При близких к нулю значениях показателя сингулярности а наблюдается запаздывание волн, сопровождаемое снижением амплитуды. Максимальное значение амплитуды перемещения при р = 0,1 и а = 0,1 снижается на 51 % по сравнению с максимальным значением в случае упругого стержня (А = 0). Отсюда следует, что влияние параметра а на амплитуду перемещения существеннее, чем влияние параметра Р (см. рис. 3). Это подтверждает полученный в работе [1] вывод о предпочтительности применения слабосингулярных ядер вязкоупругости.

С удалением от торца наблюдается затухание во времени вязкоупругих напряжений агв (см. рис. 4) в фиксированных сечениях (г = 2, 5, 8) (А = 0,1; а = 0,25; в = 0,5) стержня. Например, максимальное значение напряжения в сечении г = 2 достигается при I = 4, 8, а в сечениях г = 5 и г = 8 — при I = 8 и I = 11,2 соответственно. При этом максимальные значения вязкоупругих напряжений в сечениях г = 2 и г = 8 различаются более чем на 52 %.

В случае упругого стержня, вращающегося вокруг своей оси, из зависимостей изменения перемещения во времени при различных значениях угловой скорости О (см. рис. 5) следует, что с увеличением

угловой скорости возрастает амплитуда перемещения. В сечении z = 8 при угловой скорости Q = 0,3 перемещение увеличивается в 1,5 раза и форма колебаний изменяется.

Возрастание амплитуды перемещения с увеличением угловой скорости наблюдается и в случае учета вязкоупругих свойств (см. рис. 6). При Q = 0,3 амплитуда перемещения в вязкоупругом стержне со слабосингулярным ядром уменьшается в 1,12 раза, а в случае упругого стержня — в 1,5 раза.

Выводы. В работе предложена новая математическая модель крутильных колебаний вязкоупругого круглого стержня, вращающегося с постоянной угловой скоростью вокруг своей оси. Исследовано влияние на колебания стержня его вязкоупругих характеристик и вращения. Установлено, что учет вязкоупругости приводит к уменьшению амплитуды колебаний и запаздыванию волны, а учет вращения — к увеличению амплитуды колебаний. Создан алгоритм, позволяющий при проведении численных экспериментов варьировать параметры вязкоупругости, угловую скорость и внешние силы, действующие на стержень.

ЛИТЕРАТУРА

[1] Bauer H.F. Vibration of a rotating uniform beam. Part I: Orientation in the axis of rotation. Journal of Sound and Vibration, 1980, vol. 72, no. 2, pp. 177-189.

[2] Муницын А.И. Пространственные изгибные колебания стержня, вращающегося вокруг своей оси. Математическое и компьютерное моделирование машин и систем, 2008, вып. 3, с. 64-67.

[3] Gong S.W., Lam K.Y. Rotating multilayered cylindrical shells to impact loading. AIAA Journal, 2002, vol. 41, no. 1, Technical notes, рр. 139-142.

[4] Ng T.Y., Lam K.Y. Vibration and critical speed of a rotating cylindrical shell subjected to axial loading. Applied Acoustics, 1999, no. 56, pp. 273-282.

[5] Rand O., Stavsky Y. Free Vibrations of Spinning Composite Cylindrical Shells. International Journal Solids and Structures, 1991, vol. 28, no. 7, pp. 831-843.

[6] Бадалов Ф.Б., Абдукаримов А., Худояров Б.А. Численное исследование влияния реологических параметров на характер колебаний наследственно-деформируемых систем. Вычислительные технологии, 2007, т. 12, № 4, с. 17-26.

[7] Marynowski K. Non-linear dynamic analysis of an axially moving viscoelastic beam. Journal of theoretical and applied mechanics, 2002, vol. 40, no. 2, pp. 465-481.

[8] Горохова И.В. Малые поперечные колебания вязкоупругого стержня.

Математические заметки, 2011, т. 89, вып. 6, с. 825-832.

[9] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Конечно-элементное моделирование эффективных вязкоупругих свойств однонаправленных композиционных материалов. Математическое моделирование и численные методы, 2014, № 2, с. 28-48.

[10] Филиппов И.Г., Чебан В.Г. Математическая теория колебаний упругих и вязкоупругих пластин и стержней. Кишинев, Штиинца, 1988, 190 с.

[11]Худойназаров Х.Х. Нестационарное взаимодействие цилиндрических оболочек и стержней с деформируемой средой. Ташкент, Изд-во им. Абу Али ибн Сино, 2003, 326 с.

[12] Гузь А.Н., Кубенко В.Д., Черевко М.А. Дифракция упругих волн. Киев, Наукова думка, 1978, 308 с.

[13] Снеддон И. Преобразование Фурье. Москва, Иностранная литература, 1955, 667 с.

[14] Колтунов М.А. Ползучесть и релаксация. Москва, Высшая школа, 1976, 276 с.

[15] Формалев В.Ф., Ревизников Д.Л. Численные методы. Москва, Физматлит, 2004, 400 с.

Статья поступила в редакцию 21.03.2016

Ссылку на эту статью просим оформлять следующим образом:

Худойназаров Х.Х., Абдирашидов А., Буркутбоев Ш.М. Моделирование крутильных колебаний вязкоупругого круглого стержня, вращающегося с постоянной угловой скоростью. Математическое моделирование и численные методы, 2016, № 1 (9), с. 38-51.

Худойназаров Хайрулла Худойназарович родился в 1950 г., окончил Самаркандский государственный университет в 1972 г. Д-р техн. наук, профессор ТГТУ. Автор более 100 научных работ в области механики деформируемого твердого тела. e-mail: khayrullakhudoynazarov@gmail.com

Абдирашидов Аблакул родился в 1958 г., окончил Самаркандский государственный университет в 1980 г. Канд. физ.-мат. наук, доцент СамГУ. Автор более 100 научных работ в области механики деформируемого твердого тела. e-mail: abdira@samdu.uz

Буркутбоев Шерзод Мамасолиевич родился в 1977 г., окончил Самаркандский государственный университет в 1999 г. Старший научный сотрудник-соискатель ТГТУ. Автор 21 научной работы. e-mail: sherzodb@rambler.ru

Torsional vibrations of the viscoelastic round bar rotating with the constant angular velocity

© Kh.Kh. Khudoynazarov1, A. Abdirashidov2, Sh.M. Burkutboyev1

1Tashkent State Technical University, Tashkent, 100095, Uzbekistan 2Samarkand State University, Samarkand, 140100, Uzbekistan

The purpose of this article is to deduce general and approximate equations for the torsional vibration of the viscoelastic round bar rotating around the symmetry axis with the constant angular velocity. Within the research we develop the algorithm allowing to define the bar deflected mode. The received approximate equations enabled to numerically solve the problem of the bar torsional vibrations. Moreover, we carry out a comparative analysis of the results obtained for exponential and weakly singular kernels of the viscoelastic operator. As a result, we estimate the rotation influence on the bar vibrations.

Keywords: bar, torsional vibrations, displacement, stress, rotation, angular velocity, viscoelasticity.

REFERENCES

[1] Bauer H.F. Journal of Sound and Vibration, 1980, vol. 72, no. 2, pp. 177-189.

[2] Munitsyn A.I. Matematicheskoye i kompyuternoe modelirovanie mashin i sistem — Mathematical and Computer Modeling of Machines and Systems, 2008, vol. 3, pp. 64-67.

[3] Gong S.W., Lam K.Y. AIAA Journal, 2002, vol. 41, no. 1, Technical notes, pp. 139-142.

[4] Ng T.Y., Lam K.Y. Applied Acoustics, 1999, no. 56, pp. 273-282.

[5] Rand O., Stavsky Y. International Journal Solids and Structures, 1991, vol. 28, no. 7, pp. 831-843.

[6] Badalov F.B., Abdukarimov A., Xudoyarov B.A. Vychislitelniye tekhnologii — Computing Technologies, 2007, vol. 12, no. 4, pp. 17-26.

[7] Marynowski K. Journal of theoretical and applied mechanics, 2002, vol. 40, no. 2, pp. 465-481.

[8] Gorokhova I.V. Matematicheskiye zametki — Mathematical notes, 2011, vol. 89, no. 6, pp. 825-832.

[9] Dimitriyenko Yu.I., Gubareva E.A., Sborschikov S.V. Matematicheskoye modelirovaniye i chislenniye metody — Mathematical Modeling and Computational Methods, 2014, no. 2, pp. 28-48.

[10]Filippov I.G., Cheban V. G. Matematicheskaya teoriya kolebaniy uprugikh i vyazkouprugikh plastin i sterzhney [Mathematical theory of vibrations of elastic and viscoelastic plates and bars]. Kishinev, Shtiintsa Publ., 1988, 190 p.

[11]Khudoynazarov Kh.Kh. Nestatsionarnoye vzaimodeystviye tsylindricheskikh obolochek i sterzhney s deformiruyemoy sredoy [Non-stationary interaction of cylindrical shells and bars with deformable medium]. Tashkent, Abu Ali ibn Sino Publ., 2003, 326 p.

[12]Guz A.N., Kubenko V.D., Cherevko M.A. Difraktsiya uprugikh voln [Diffraction of elastic waves]. Kiyev, Naukova Dumka Publ., 1973, 308 p.

[13] Sneddon I. Preobrazovaniye Furye [Fourier transforms]. Moscow, Inostrannaya Literatura Publ., 1955, 667 p.

[14]Koltunov M.A. Polzuchest i relaksatsiya [Creep and relaxation], Moscow, Visshaya Shkola Publ., 1976, 276 p.

[15]Formalyov V.F., Reviznikov D.L. Chislennyye metody [Numerical methods]. Moscow, Fizmatlit Publ., 2004, 400 p.

Khudoynazarov Kh. Kh. (b. 1950) graduated from the Faculty of Mathematics and Mechanics, Samarkand State University in 1972. Dr. Sci. (Eng.), Professor of Tashkent State Technical University. Author of more than 100 scientific works in the field of deformable solid mechanics. e-mail: khayrullakhudoynazarov@gmail.com

Abdirashidov A. (b.1958) graduated from the Faculty of Applied Mathematics of Samarkand State University in 1980. Cand. Sci. (Phys.-Math), Assoc. Professor of the Department of Applied Mathematics, Samarkand State University. Author of more than 100 scientific works in the field of deformable solid mechanics and computational mechanics. e-mail: abdira@samdu.uz

Burkutboyev Sh.M. (b.1977) graduated from the Faculty of Mathematics and Mechanics, Samarkand State University in 1999. Senior Research Scientist, Tashkent State Technical University. Author of 21 scientific works. e-mail: sherzodb@rambler.ru