CC BY
55
УДК 539.3
А.Ф. ВЕРЛАНЬ, д-р техн. наук, Ин-т проблем моделирования в энергетике им. Г.Е. Пухова НАН Украины, Киев,
Б.А. ХУДАЯРОВ, д-р техн. наук, зав. каф., Ташкентский институт ирригации и мелиорации, Ташкент,
Э.Ф. ФАЙЗИБОЕВ, проф., Ташкентский институт ирригации и мелиорации, Ташкент,
З.У. ЮЛДАШЕВ, ст. преп., Ташкентский институт ирригации и мелиорации, Ташкент
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ФЛАТТЕРА ВЯЗКОУПРУГИХ ОРТОТРОПНЫХ ПЛАСТИН В СВЕРХЗВУКОВОМ ПОТОКЕ ГАЗА
В работе рассматривается нелинейный флаттер вязкоупругих ортотропных пластин в сверхзвуковом потоке газа. Уравнения движения получены на основе наиболее общей теории пластин Фёппля-Кармана в перемещениях, для которых справедливы гипотезы Кирхгоффа. Аэродинамическое давление учитывается согласно поршневой теории А.А. Ильюшина. При помощи метода Бубнова - Галеркина, основанного на многочленной аппроксимации перемещений, задача сведена к исследованию системы обыкновенных интегро-дифференциальных уравнений (ИДУ), где независимой переменной является время. Решение ИДУ находится численным методом, основанным на использовании квадратурных формул. Разработаны методика и алгоритм численного решения интегро-дифференциальных уравнений. Определены критические скорости флаттера ортотропных пластин в потоке газа. Табл.: 1. Библиогр.: 9 назв.
Ключевые слова: вязкоупругая ортотропная пластина, алгоритм, интегро-
дифференциальное уравнение, флаттер, пластина.
Постановка проблемы и анализ литературы. Большинство конструктивных элементов летательного аппарата, покрытых тонкой обшивкой, можно схематизировать либо в виде прямоугольных, либо в виде цилиндрических панелей или оболочек. Различие в геометрии этих элементов не является чисто формальным признаком и в значительной степени обуславливает специфику их аэроупругого поведения. В силу этого задачи моделирования панельного флаттера можно разбить на две большие группы. Первая включает задачи, относящиеся к плоским пластинам, вторая содержит задачи, связанные с цилиндрическими панелями и оболочками. Математическому моделированию флаттера упругих элементов второй группы посвящено большое количество теоретических исследований. Однако моделированию задач флаттера вязкоупругих пластин, цилиндрических и пологих оболочек посвящено сравнительно небольшое количество работ. Последнее объясняется специфическими аналитическими трудностями исследования пластин, ISSN 2079-0031 Вестник НТУ "ХПИ", 2012, № 62 (968)
8
цилиндрических и пологих оболочек. Потребности новой техники, в особенности бурно развивающейся авиационной промышленности за последние пятьдесят лет, определили необходимость развития нелинейной теории пластин и оболочек в сверхзвуковом потоке газа.
В настоящее время задача учета наследственных эффектов деформируемых материалов представляет большой теоретический и прикладной интерес. Ее решение представляет собой эффективное приложение теории вязкоупругости к реальным процессам. Поэтому методы и проблемы теории наследственной упругости привлекают большое внимание исследователей.
Математическое моделирование динамических задач о колебаниях и устойчивости вязкоупругих систем является также весьма актуальным в связи с тем, что, с одной стороны, расширяются возможности использования в авиационной промышленности и других отраслях машиностроения материалов, обладающих ярко выраженными вязкоупругими свойствами, а с другой стороны, при использовании наследственных моделей [1 - 4] для описания внутреннего демпфирования материала уравнения колебаний упругих систем записываются в такой же форме, как и для вязкоупругих систем. Часто при рассмотрении упругих систем внутреннее трение материала учитывается с помощью модели Фойгта, хотя известно, что даже в системах с конечным числом степеней свободы, большим единицы, она приводит к некорректным результатам, поскольку для большинства материалов внутреннее трение фактически не зависит или, по крайней мере, слабо зависит от скорости колебаний в достаточно широком частотном диапазоне. В этом смысле более предпочтительной является модель, отображающая наследственные свойства [5].
Задача нелинейного флаттера ортотропной прямоугольной пластинки, обтекаемой сверхзвуковым потоком газа была впервые рассмотрена в работе [6], где для аэродинамических сил использовалась приближенная формула "поршневой теории". В работе [7] задача решается в нелинейной постановке для трехслойной ортотропной пластинки. Позднее [8] эта задача рассматривалась для упругих ортотропных панелей, а динамический процесс рассмотрен без учета распространения упругих волн, однако в вязкоупругой постановке задача не исследована.
В связи с этим целью данной статьи и является теоретическое исследование нелинейного флаттера вязкоупругих ортотропных пластин, обтекаемых в потоке газа.
Постановка задачи. Рассмотрим прямоугольную вязкоупругую ортотропную пластинку со сторонами а и b, которая обтекается с одной стороны сверхзвуковым потоком газа со скоростью V. Аэродинамическое давление учитываем по поршневой теории А.А. Ильюшина [2].
Нелинейные уравнения движения вязкоупругих ортотропных пластин выпишем в виде:
(
б:
2,.Л
du dw д w
+
dx2 dx dx
(
+ Б
12
d2 V
dw d w
dxdy dy dxdy
Í я2
+ Б
d u d V dw d w dw d w
-----T" +-----------I----------------I-------------
dy 2 dxdy dy dxdy dx dy 2
du n Фту = °
dt2
b:
d2V dw d2w
dy 2 dy dy 2
+ Б2!
2
d2u
+
dw d2w
+ Б
2
d2u
d 2 V dw d 2 w dw d 2 w
dxdy dx 2 dx dxdy dy dx
dxdy dx dxdy
P
d2V
= °,
dt2
(1)
(2)
— \B 12
* d W 11
dx
d Idw dx | dx
d Idw dy I dy
4 +|8B + БІ2 + B2*ij
d4w , n*
dx2 dy2 22
d 4 w|
* du * dv
Bii^T + B12^T dx dy
* du * dV
B21^T + B22 _ dx
dy
+ dwB*\dU + dV |l_
dy ^dv dx )j
^•f* +dV ÏL
dx \ dy dx H
d2w Б dw BV dw B1V2 f dw|2 _Q
dx ) ,
+ p—— +------------------+-+ -
" h dt h dx
dt
h
(3)
t t : R*(p(t) = JR (t -т)ф(т)й?т, R*(p(t) = JR (t -х)ф(х)й?х ; R(t -x).
где В* = Ъ - (1-Я*); (/', ] = 1, 2); В* = Ъ(1-R*); Я*-, R* - интегральные г г
операторы: Я*ф(г) =.я ^ *)ф(*)^*, я-Т(.) .я-ч
о о
Я- (г - х) - ядра релаксации, определяемые при рассмотрении одноосной
задачи; р - плотность материала; к - толщина оболочки.
Граничные условия будут иметь вид:
w = 0; V = 0; N = 0; Мх = 0 при х = 0, х = а; w = 0; и = 0; Ыу = 0; Му = 0 при у = 0, у = Ь.
*
+
+
+
+
+
2
4
Этим условиям удовлетворим, представляя искомые функции
и(х, у, I), v(x, у, ^ и w(x, у, 0 в виде
N M
пкх . тку
u
(X, у, t) => > ипт(ґ)С08---------в1П
" а Ь
п=1 m=1
N M
' П
V (Ґ )sln- _____
а Ь
, . ППХ тпу
V(X y, t) = ^ X (t--------------------------------с0^ ~Г~~ , (4)
п=1 ш=1
N М
, ппх шпу
*(X y, t) = X X *пт (t) 81п-------------81п
11 а Ь
п=1 ш=1
Исследуем случай пластинки, шарнирно опертой по контуру, совершающей колебания с образованиям пит полуволн вдоль координатных линий х и у. Решение уравнений (1) - (3),
удовлетворяющее граничным условиям задачи, будем искать в виде (4). Подставляя (4) в систему (1) - (3) и применяя метод Бубнова-Галеркина, получим систему ИДУ. Введя в ИДУ следующие безразмерные величины х у и V
---, ,-----, ----, -------, —, и сохраняя при этом прежние
а Ь Ь Ь ї а
обозначения, запишем
и и + (^С1 1 + а11С )ик1 + (а12С1 2 + аізс ~)ук1 +
■Л М Г 1 (5)
+ X X [ (а16С1 і + а17С )\кіптіг “ (а18С1 2 + а19С )^2кІптіг \ *пт*іг = 0,
п,і=1 т,г=1
^ кі + (а20С22 + а21С ^укІ + (а22С2 і + а23С )икІ +
^ М (6)
+ X X [ (а26С2 2 + а27С )^3кіптіг ~ (а28С2 і + а29С )^4кіптіг \ *пт*іг = 0,
* м +М ы + (а30СП + +а31Сі*2 + а32С2*2 + а33С* )кЫ + Х X Ь2МптігСП™пт™іг +
п,і=1 ш,г=1
N М (
+ X X Кпт\^^1кіптігС11 + ^2кіптігС21 + ^3кіптігС )иіг + (^4кІптітС12 + п,і=1 т,г=1
+ ^5кіптітС22 + ^6кіпшігС )vir + (^7кіптітС11 + ^8КкушігС12 + ^9кіпшігС22)Кіг }+
Г N К+1 2 N М
+ Х М I 2^М * Х^пЛ + 7 М * X X Гк!пті^птЩГ
V п=1 4 п,і=1 т,г=1
п,і=1 т,г=1
где
a10 =п2к 2 A; au = я21 2X2g; a12 =n2kl XA; a13 =n2kl Xg;
a16 = nni2A/Xj; a17 =nnr2X2g/Xt; a18 =nmirX?A/Xt;
a19 =nmirX2g / X1; a20 =n2X2l2/A; a21 =n2k2 g; a22 = %2Xkl / A;
a23 =n2kl Xg; a26 = rcX3r 2m /(AX1); a27 =nmi 2Xg / X1;
a28 = nXnir /(AX1); a29 = nnirXg / X1;
я 4 к 4 A X4 я 4 A X4 я 4 g 1 X4 я414
------■ - -------■ a„ =------------■ - ---------
30 = ; a31 = ; a33 = ; a32 . , 2 ;
12X] 6X1 3X1 A 12X1
d1klnmir = X1niA%(nA7 klnmir + iA 5klnmir); d5klnmir
X 2П1(rA6klnmir
mA7klnmir )/(AX1); d2klnmir = ß1nAkx (iA5klnmir nA Iklnmir ); d3klnmir = rX2A^-(imA6klnmir + nrA5klnm,r - 2nrA8klnmir) / X1; d8klnm,r = ^m™ X x (Xr A 6klnmir - mA 7 k ln mir ) /(AX1 ); d9klnmir mX (X1rkyA6klnmir mß1kxAIklnmir )/ A; d6klnmir = ( - nmrgяXЛ8klnm,r (X + 1) + ^'^g^ A6Unmir + nr ^g^ A5klnmir ) / Xb d7klnmir = ( - nmІgяXA8klnmir (X + 1) + mi VXA6kl«mir + nirKgXA5klnmir ) / X^ d4klnmir = nß1ky A(iA5k/nmir nA7klnmir ) + mX2kx (rß1A 6klnmir mA7klnmir / X1)/A;
^* _ 1 - R1 1 ^ 2 (1 - R1 2). Г1 * _Н-1(1 - R2 1).
C11 = , ; C12 = , ; C21 = , ;
1 -1^Ц 2 1 -l^ 2 1 -l^ 2
1 - R* IE П
C2*2 = 1 R12 ; C* = 1 -R*; A = РЦ g = -=^=.
1 -H-1^ 2 V E2 л/е1 e2
Требуется найти критическую скорость флаттера F^, из решения системы ИДУ (5) - (7), удовлетворяющую начальным условиям:
unm(0) u0nm, unm(0) u0m , vnm(0')~v0nm, ^nm(0) ^0m wnm(0) w0nm,,
^nm(0) = 4m •
Для решения задачи нелинейного флаттера вязкоупругих ортотропных пластинок, описываемой системой ИДУ (5) - (7), используем численный метод, основанный на применении квадратурных формул [9, 10].
Интегрируя систему (5) - (7) два раза по t, запишем ее в интегральной форме и с помощью рационального преобразования исключим слабо-сингулярние особенности интегрального оператора R . Затем, полагая t = t, t, = iAt, i = 1, 2, ... (At = const) и заменяя интегралы ISSN 2079-0031 Вестник НТУ "ХПИ", 2012, № 62 (968)
12
квадратурными формулами трапеций для вычисления игк1 = ик1(^), Vikl = vkl(tг) и wikl = wkl(ti), получим следующие рекуррентные формулы для ядра Колтунова - Ржаницына:
ирМ = и он + и° (
р-1 I А
А11
+ и ом!р -£ Ау (!р - )]аю --------11 £ В<111 ехр(-Ри/^ )и
5 )uj-skt
+
а
11
М12
А ч ' а -2 А V '
-М1------£ В5 еХР(-Р?5 )и --5« +— vjkl - ~ £ В 5 еХР(-Р?5 V
5=0
М-12
У ум ' £ Б 5 ’
а "
5=0
(
+ а12М 2
У ум '
- ~ £ В5 2) еХР(-Р 12 >5 >
12 5 =0
££
г
+ А,.
А1к1пт1га16 (
а17 М12
" А 2Мпт/га18М 2
М12
7-пт™-гг - — £ В5П) еХР(-Р11 ¡5)
а11 5=0 - А £ В5 еХР(-Р )
а 5=0
иг-—£ В512) еХР( "Р12 >5 )™
(8)
1-5Пт --51
42 5=0
5 / ^-5пт ^~51Г
А х~'
“ £ В5 еХР("Р г5 )^--5Пт"--51Г
Урк? У0к? + У ок1 >р -£ А (<р - О а201 ------------22 £ В?1) еХР(-Р 22г5 >/-5« 1 +
М12
+------21 а21 --£В5 еХР(-Р>5 )У--5к1 1+------------231 и-к?-£В5 еХР(-Рг5 )и--5к1 1 +
М121 а 5=0
+ а22М 11 и-к1 - — £В521) еХР(-Р21г5 >--*1 1 +
££
+ А,
М12
А -
™-пт™Лг - — £ В522) еХР("Р22 <5 --5пт™ -
а22 5=0
А -
^-пт^-Г - — £ В5 еХР(-Р<5 >-^пт™--51Г
5=0
А 4Мптгга28М1
- —21 £ В521) еХР(-р21 >5 >
п ..
21 5=0
- А4
М12
™ ™ ..----------------> В
у -пт у-гг *
А
----£ В5 еХР(-Р ¡5
п ■>
М12 “ 0,
+
+
+
+
^-5пт J-5гг
а
А
5=0
-=0
5=0
5=0
А 3к1птгга26
+
п.г=1т.г=1
а
27
™ ™ ■■ у -пт у-гг
а
29
5=0
(
урМ '
М12
1 + АрЫ х
™0к? +( ™0к? + МXМ,0Ы )гр-£ А--=0
р-1
— АМ„
М12
— +1 *2
N М
рI 2ХгМ*£укп™-п? +—— М* £ £
4
Л ] А11^
+ «зс| -----“£Я,(11) ехр(-р11 ts)^.-,4г 1 + 0^2 \ w
п,г=1 т,г=1 .
-к1
'А"2 £ Вв<12) ехр(-рп tí] 1 +
•*12 в=0
А,-
+ а„\ - А22 £в<22)е>ф(-р22 ^ ;>]« |+ ~ ( Wyk^--------£Вв еХР(-Р ts >]-вИ | +
32 \ уу]к1
а
^22 в=0
N М
+ £ £ М2Ь2
п,г=1 т,г=1
М12 V а в=0
- А12£в(12) ехр(-Р19 г ^. | +
в т V А-^12 в/ ]-впт ]-в^ |
п,г =1т,г =1
А,
+ £ £ ^1к1пт[г \ и-гг--------------11 £ ВВ‘ ^ еХР (-Р11 tв )и]-в„- | +
А
+ М^2к 1п,»г I и]г - ~ еХР(-Р21 ^ )и]-вгг 1+ -3к= \ ]-------------£ Д
М12
А
(10)
X ^С-^ ^--в* | + М2 ^ 4кМт{г | V]iГ - ~еХР(-Р12 ^ >]-в/г | +
а12 в=0
X ехр(-р tв )V]^ | + М2^!
22 в >г] -siг
А
^КМптг \ >]---------12£ -^В1 2> еХр (-Р1 2 ^ )
а,
12 в=0
+ ^,г - ~£ В(11) >]-..г | + ¿9Шт,г[^,г - —£ БВ22) еХР( Р22 ^ )w
где р = 1, 2, ...;к = 1, Ж; / = 1, М .
На основе этого алгоритма создан пакет прикладных компьютерных программ. Результаты вычислений представлены в таблице.
V
N
+
п=1
12 в=0
11 в=0
X
21 в=0
в=0
22 в=0
11 в=0
22 в=0
Таблица
Зависимость критической скорости флаттера от физико-механических и геометрических параметров пластинки
Д g Ц2 X XI материал панели Ур
5,0129 4,8482 3,0976 2,2761 2,1467 1,1088 0,0226 0,1549 0,1398 0,1356 0,3 0,26 0,02 0,15 0,13 0,012 0,017 0,32 0,01963 0,023 2,5 50 вязкоупругий 980 892 615 463 440
5,0129 1,1088 0,3 0,012 2,1 2,3 2,6 2,8 50 вязкоупругий 532 734 1124 1433
5,0129 1,1088 0,3 0,012 2,5 40 45 50 вязкоупругий 1022 725 532
5,0129 1,1088 0,3 0,012 2,5 50 вязкоупругий 985 1062 1341
2,2761 0,1398 0,15 0,01963 2,5 50 упругий вязкоупругий 625 463
Проведены расчеты для вязкоупругих ортотропных пластин применительно к пяти маркам стеклопластика: Д = 5,0129; g = 1,1088;
Ц1 = 0,3; ц2 = 0,012 (марка 1); (Д= —1; g = —=■; Еь Е2 -
\Е2 V — ,—2
линейные модули упругости; О - модуль сдвига; ць ц2 - коэффициенты Пуассона); Д = 4,8482; g = 0,0226; щ = 0,26; ц2 = 0,017 (марка 2); Д = 3,0976; g = 0,1549; ц = 0,02; ц2 = 0,32 (марка 3); Д = 2,2761; g = 0,1398; ц = 0,15; ц3 = 0,01963 (марка 4); Д = 2,1467; g = 0,1356; щ = 0,13; ц2 = 0,023 (марка 5). Результаты расчета для этих марок стеклопластика представлены в первых пяти строках таблицы. Видно, что критическая скорость (при постоянных X и X]) сильно зависит от свойств материала.
Изучено влияние вязкоупругих свойств пластинки для марки 4. Критическая скорость для упругой пластинки составляет 625 м/с, а для вязкоупругой пластинки 463 м/с. Разница между ними составляет 25,9%.
Исследовалось влияние параметра X (отношение длины стороны а к Ь) на критическую скорость флаттера вязкоупругих ортотропных пластин, материал которых относится к стеклопластикам марки 1. С ростом параметра X увеличивается протяженность пластины в направлении течения и происходит сближение удлиненных краев
пластины. Последнее способствует повышению относительной жесткости и росту критической скорости, которое можно проследить по табл. При X = 2,1 эта скорость равна 532 м/с. С увеличением параметра X до 2,8, критическая скорость возрастает в 2,69 раза.
Далее исследовано влияние параметра относительной толщины Х1 (отношение стороны а к толщине h пластинки) на критическую скорость флаттера. Для Х1= 40, 45 и 50 значения соответственно равны 1022, 725 и 532, то есть с уменьшением толщины стенок пластинки (с увеличением параметра Х1) критическая скорость флаттера уменьшается.
Выводы. На основании полученных результатов можно заключить, что учет вязкоупругих свойств материала пластинки приводит к уменьшению критической скорости флаттера Vkp, с которой начинается явление флаттера.
Отметим также, что при скорости потока меньшей, чем Vкр, влияние вязкоупругого свойства материала уменьшает амплитуду и частоту колебаний. Если же скорость потока превышает Vkp, то вязкоупругое свойство материала оказывает уже дестабилизирующее влияние.
Список литературы: 1. Григолюк Э.И. Нелинейное деформирование тонкостенных конструкций // Э.И. Григолюк, В.И. Мамай. - М.: Наука: Физматлит, 1997. - 272 с. 2. Ильюшин А.А. Закон плоских сечений в аэродинамике больших сверхзвуковых скоростей // А.А. Ильюшин // Прикладная математика и механика. - 1956. - Т. ХХ. - Вып. 6. - С. 733755. 3. Кочнева Л.Ф. Внутреннее трение в твердых телах при колебаниях // Л.Ф. Кочнева. -М.: Наука, 1979. - 96 с. 4. Сорокин Е.С. Об учете упругих несовершенств материалов методами наследственной упругости // Е.С. Сорокин, Г.Б. Муравский. — Строит. механика и расчет сооружений. - 1975. - № 4. - С. 41—46. 5. Потапов В.Д. Исследование динамической устойчивости вязко-упругих систем с помощью показателей Ляпунова // В.Д. Потапов. — Известия АН. Механика твердого тела. - 2000. - № 6. - С. 82-90. 6. Амбарцумян С.А. Об устойчивости ортотропных пластинок, обтекаемых сверхзвуковым потоком газа // С.А. Амбарцумян, Ж.Е. Багдасарян. — Известия АН СССР, ОТН, Механика и машиностроение. - 1961. - № 4. - С. 91-96. 7. Багдасарян Ж.Е. Об устойчивости трехслойной ортотропной пластинки в сверхзвуковом потоке газа // Ж.Е. Багдасарян. — Известия АН Арм. ССР, Физико-математические науки. - 1961. - XIV. - № 5. - С. 21-29.
8. Jehad F.A. Nonlinear Flutter of Orthotropic Composite panel Under Aerodynamic Heating // F.A. Jehad, R.A. Ibrahim, F.G. Ronald. — J. AIAA. - 1993. - Vol. 31. - № 8. - P. 1256-1268.
9. Бадалов Ф.Б. О некоторых методах решения систем интегродифференциальных уравнений, встречающихся в задачах вязкоупругости // Ф.Б. Бадалов, Х. Эшматов, М. Юсупов. — Прикладная математика и механика. - 1987. - Т. 51. - № 5. - С. 867-871.
10.Верлань А.Ф. Численное решение нелинейных задач динамики вязкоупругих систем //А.Ф. Верлань, Х. Эшматов, Б. Худаяров, Ш.П. Бобоназаров. — Электронное моделирование. - 2004. - Т. 26. - № 3. - С. 3-14.
УДК 539.3
Комп’ютерне моделювання флаттеру в’язкопружних ортотропньїх пластин в надзвуковому потоці газу / Верлань А.Ф., Худаяров Б.А., Файзибоев Е.Ф., Юлдашев З.У. // Вісник НТУ "ХПІ". Серія: Інформатика та моделювання. - Харків: НТУ "ХПІ". - 2012. - № 62 (968). - С. 8 - 17.
У роботі розглядається нелінійний флаттер в'язкопружних ортотропньїх пластин у надзвуковому потоці газу. Рівняння руху отримані на основі найбільш спільної теорії пластин Феппля-кармана в переміщеннях, для яких справедливі гіпотези Кірхгофа. Аеродинамічний тиск враховується згідно поршневої теорії А.А. Іл‘юшина. За допомогою методу Бубнова-Галеркіну, заснованого на багаточленній апроксимації переміщень, завдання зведене до дослідження системи звичайних інтегро-дифференційних рівнянь (ІДУ), де незалежною змінною є час. Рішення ІДУ знаходиться чисельним методом, заснованим на використанні квадратурних формул. Розроблені методика і алгоритм чисельного рішення інтегро-дифференційних рівнянь. Визначені критичні швидкості флаттеру ортотропньїх пластин у потоці газу. Табл.: 1. Бібліогр.: 10 назв.
Ключові слова: в'язкопружні ортотропні пластини, алгоритм, інтегро-дифференційне рівняння, флаттер, пластина.
UDC 539.3
Computer modeling of flutter of viscoelastic orthotropic plates in supersonic flow / Verlan A.F., Khudayarov B.A., Fayziboyev E.F., Yuldashev Z.U. // Herald of the National Technical University "KhPI". Subject issue: Information Science and Modeling. - Kharkov: NTU "KhPI". - 2012. - №. 62 (968). - P. 8 - 17.
This paper considers the nonlinear flutter of viscoelastic orthotropic plates in supersonic flow. The equations of motion derived from the most general theory of plates Foppl-Karman in displacements, which satisfy the Kirchhoff hypothesis. Aerodynamic pressure is accounted for under the piston theory A.A. Ilyushin. Using the Bubnov - Galerkin method based on polynomial approximation of displacements, the problem is reduced to the study of a system of ordinary integro-differential equations (IDE), where the independent variable is time. The IDE is a numerical method based on quadrature formulas. The technique and algorithm for the numerical solution of integro-differential equations. The critical velocity of flutter of orthotropic plates in a gas stream.Computional experiment in task about flutter of viscoelastic plate, streamlined under arbitrary corner, in stream of gas. Tabl.: 1. Refs.: 10 titles.
Keyworts: viscoelastic orthotropic plates, algorithm, integro-differential equations, flutter, plates.
Поступила в редакцию 06.08.2012
