УДК 539.3
Б.А. ХУДАЯРОВ, д-р техн. наук, зав. каф., Ташкентский институт
ирригации и мелиорации, Ташкент,
Р. АБДИКАРИМОВ, доц., Ташкентский финансовый институт,
Ташкент
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЗАДАЧИ ДИНАМИКИ
ВЯЗКОУПРУГИХ СИСТЕМ С ПЕРЕМЕННОЙ ЖЕСТКОСТЬЮ
В работе на основе интегральных моделей разработаны и развиты расчетные математические модели для исследования нелинейных колебаний вязкоупругих пластин, круговой цилиндрической панели и цилиндрической оболочки с переменной жесткостью. Ил.: 1. Библиогр.: 8 назв.
Ключевые слова: математическая модель, переменная жесткость, вязкоупругие системы, нелинейные колебания, задача динамики.
Постановка проблемы и анализ литературы. Пластины переменной толщины широко используются в качестве подпорных стенок железобетонных зданий, прямоугольного элемента резервуаров, контрфорсных плотин, бетонного покрытия для укрепления взлетной полосы аэропорта, фундаментных плит сооружений и несущих элементов метрополитена, элементов подземных сооружений. Также они широко внедрены и в кораблестроение, судостроение, самолетостроение, и многих других отраслях техники и промышленности. В связи с интенсивным развитием промышленности, большое развитие приобрела механика композитных материалов. Интерес к проблемам деформирования и прочности пластин, панелей и оболочек из
композиционного материала связан с тем, что они представляют собой основные несущие элементы конструкций. Использование новых
композиционных материалов в инженерной практике, а также проектирование и создание прочных, легких и надежных конструкций требует совершенствования механических моделей деформируемых тел и разработки более совершенных математических моделей их расчета с учетом реальных свойств конструкционных материалов и их геометрии [1 - 4]. Поэтому разработка эффективных алгоритмов, которые
используются для решения нелинейных задач динамической
устойчивости элементов тонкостенных конструкций из композиционных материалов, является актуальной задачей.
В последнее время большое внимание уделяется изучению динамики существенно нелинейных вязкоупругих механических систем.
© Б.А. Худаяров, Р. Абдикаримов, 2013
Начало теоретическим исследованиям нелинейных колебаний упругих оболочек положено Э.И. Григолюком [5] и Е. Рейсснером [6], которые разработали основную нелинейную модель, использовав приближенные вариационные методы. В качестве разрешающих уравнений в данных работах использовались геометрически нелинейные уравнения теории тонких оболочек, а динамический прогиб и функции напряжений в срединной поверхности аппроксимировались одночленными выражениями. Для определения обобщенного перемещения получено нелинейное дифференциальное уравнение, приближенное решение которого определялось с помощью метода Бубнова-Галеркина и метода возмущений. Последующие исследования нелинейных собственных и вынужденных колебаний упругих систем выполнены многими авторами. Наиболее существенные результаты получены С.А. Амбарцумяном, Н.Х. Арутюняном, Г.Е. Багдасаряном, Ф.Б. Бадаловым, А.Е. Богдановичем, В.В. Болотиным, А.С. Вольмиром, И.И. Воровичем, К.З. Галимовым, В.Ц. Гнуни, Я.М. Григоренко, А.Н. Гузем, А.А. Ильюшиным, В.С. Калининым, В.Г. Карнауховым, Г. Каудерером, Я.Ф. Каюком, И.Г. Кельдибековым, И.А. Кийко, Н.А. Кильчевским, М.С. Корнишиным, А.С. Кравчуком, В.Д. Кубенко, В.А. Крыско, Г.В. Мишенковым, Х.М. Муштари, Ю.В. Петровым, Г.С. Писаренко, Е.И. Шемякиным, А.Н. Филатговым, Л.Г. Донеллом, Дж. Навинским, Е. Рейснером и др.
Целью данной работы является разработка комплекса математических моделей в форме систем интегро-дифференциальных уравнений с сингулярными ядрами, численных алгоритмов и компьютерных программ для решения широкого класса задач нелинейных колебаний и динамической устойчивости вязкоупругих тонкостенных элементов конструкций с переменной жесткостью.
Постановка задачи. Построим математическую модель задачи о нелинейном колебании вязкоупругой изотропной пластины и оболочек переменной толщины (рис.) в геометрически нелинейной постановке по кинематической гипотезе Кирхгофа-Лява.
Рис. Вязкоупругая изотропная пластина и оболочка переменной толщины
Интегральную модель Больцмана-Вольтерра, которая характеризует закон изменения между напряжениями стх, сту, ст^ и деформациями
ех, еу, е ху в срединной поверхности, запишем в интегральном виде [7]
-Е2 ( -Г* )(е х + Це у); ст у = 7^2 (! -Г* )(е у + Це у)
1 -ц2 4 " ' 1 -ц
ст = ,Е (1 -Г* )е, (1)
^ 2(1 + ц) ’ ху
где Е - модуль упругости; ц - коэффициент Пуассона; Г* -
интегральный оператор; ех, еу, еху - компоненты конечной
деформации.
Связь между деформациями в срединной поверхности е х, еу, у^ и перемещениями и, V, w по направлениям х, у, 2 имеют вид [8]:
ди 1 ( дw I2 дv 1 ( дw I
е х =----+ —I — I ; е у =-+ — I— I ;
дх 2 ^ дх1 у ду у 2 [ ду) ^
ди дv дw дw
ду дх дх ду
Изгибающие и крутящие моменты элемента оболочки примем в виде [6]:
М = -О (1 -Г* )(д^ + цд^|, х о у, Н = -О (1 -ц)1-Г*) [дх ду )
. (3)
дхсу
Ек 3 (х, у)
где О =------------^-переменная цилиндрическая жесткость оболочки.
12(1 -ц 2)
Уравнения движения элемента вязкоупругой изотропной оболочки имеют вид [8]:
ЭЫх 5Ы д 2и . дМ дМу д \
+------- + Рг -рк Г- = 0 ; ---- +----- + п -ок
дх ду дг дх ду у дг
д2М дМ л2Н , ,т , ,т д (,т дw ,т дw
‘у ~ д Н , , г , -.г д I ^ г д^ ^ т дw I /Л\
+-----о—+ 2-----------+ кхМх + у +--1 Nх---------------+ N — I + (4)
дх2 ду2 дхду х х у у дх | х дх ^ ду)
д (АГ дw дw I д2 w
+------1 Nху---+ N' — I + а - рк —— = 0 .
ду [ ху дх у ду) а р дг2
Здесь Nx, Ny и N - усилия, отнесенные к единице длины сечения
оболочки:
Nx = стхк , х 0 У, Nxy = Ухук , (5)
рх, р^ и а - интенсивность заданных внешних нагрузок, приложенных
к элементу по направлениям х, у и г соответственно.
Подставляя (3) и (5) в (4), получим систему нелинейных интегро-дифференциальных уравнений с переменными коэффициентами вида:
('-г* )
де-
йт
+ ц-
де у дх
1-ц_ду51Л
2 ду
дк I .
+йх (е х+це у;+
1 -ц дк
2 ду
ху
1 -ц 2
Е
Рх -Рк
1 -ц д и Е д2
= о;
+цдк+1-ц£Уху
ду ду 2 дх
дк ( \
+ д^(е у+цех )^^-у
1 -ц дк
~Г дху ху
1 2 1 2 ^2 1 -ц ; 1-ц д V _
+ ——Ру-рк—с - = о ;
Е
Е дг
2
(6)
('-Г* )|
ОУ4 w + 2 — — V 2 w + 2— — V2 w + У2 ОУ2 w -
дх дх
-(1 -ц)
д2О д2w „ дд2лм д2О д2w
- 2-------+■
дх2 ду2 дхду дхду ду2 дх2
ду ду
Ек
1-ц
2 [(кх +цку )ех + (цкх + ку )еу ]}+
дк
ду
ду ('-Г*)(еу +-х)+(-Г*Ьху
+а -рк
_д^ дг2
= о.
С учетом (2) система уравнений (6) относительно перемещений и, V, V примет вид:
(. „Л Г, д2и I, \дw 1 -цд2ы 1 + ц д2v дw д2^
(1 -Г ЦИ—--(кг+цк„)— + —-—- + — ----------+------- +
| дх2 дх 2 ду2 2 дхду дх дх2
1-цдw дw 1 + цйw дw + —-------т- + -
2 дх ду2 2 ду дхду
дк
дх
^ + ц^. -(кх+цку ) „ +
[дх ду
1 ( дw
21 дх ) 2 I ду
1 -ц дк (ди дv дw дwI! (1 -ц2)рк д
2 ду I ду дх дх ду
+ — + — —!> + Рх— = 0;
Е
дг2
1л т-*\ I > д2v \дw 1 -ц д2v 1 + ц д2и дw д2^
(1 -Г Цк —- - (к у+цкх)— + —-—- +— ----------------------------------------+------------ +
ду у ду 2 дх2 2 дхду ду ду
к
+
+
+
2
2
+
1 -ц dw d2 w 1 + Hdw d2 w 2 dv dx2 2 dx dxdy
dh
"dv
^ + ц^-(kv+Hkx ) w +
dv dx v y x>
1 ( dw I ц ( dw I2
2 I|dv J 2 I dx J
1 -ц dh ( du dv dw dw 11 (l -ц2 )ph d 2v
2 dx I|dv dx dx dv Jj Py 77 я'2
(l-Г* ) h3
^d4 w d 4w +d 4 w
dx4 dx 2dy2 dv4
4
[
+ 3 (l-Г*)
J
( dh I2 ,2d2h 2h I — I +h2—-
[dx J dx
6 (l - г* )h
dh
3_. A
dx
+3 fl -г* )
d w d w +
dx3 dxdv2
,(l-г*) )
2hI
6 (l - г* )h
h2
ddh
dv
(
d w d w +
E dt2
d2w d2w
—т+ц—T
dx2 dv2
,3... ^
dv3 dx2dy
d 2h_ ,,2
dv2
22 d w d w
9 + H" 9
dv2 dx 2
+
6 (1 - ц) (l-г* )
dh dh 2 d2h
2h—— + h
dx dv dxdv
d 2 w
dxdv
(7)
-12 (l-Г ) h [(kx +ykv )^+ ^Hkx + ky )~ -ik!2 +k2 + 2Hkxky ) w +
kx +Hkv ( dw I kv +Hkx Í dw
dx
dv
-12dw (l-Г* )jh
1 -Hdu 1 +h d v dw dw 1 -ц dw d w 1 + ц dw d w
2 dv2 2 dxdv dx dx2 2 dx dv2 2 dv dxdv
dh
dx
du
dx
dv L , \ 1 (dwI2 ц Í dw
+ ^-fkx +Hkv j w+t|^-| +d^-
dv
2 [ dx J 2 I dv
1 -ц dh (du dv dw dw 11
2 dv [dv dx dx dv J J
,-dw/, *\, du dv/, , \ 1 (dw I ц Í dw
12ТТІ1 -Г )h ^T + H^-\kx + Hky)w+t!^I ^IlT
dx
dv dv
2 [dx J 2 [dv
-12dw (l -Г* )jh
d v (k t^dw 1 -Hd v l + ц d u
dv2 y x dy 2 dx2 2 dxdv
dw d2w 1 -ц dw d2w 1 + ц dw d2w
dv dv2 2 dv dx2 2 dx dxdy
dh
dv
dv du
------+ц----------
dv dx
г , \ 1 ( dw I ц ( dw I2
- Iky + Hkx ) w +— I — I +— I — I
Vy x 2 [dyJ 2 [dxJ
1 -ц dh (du dv dw dw 11 2 dx I^ dy dx dx dy J J
+
2
+
2
+
2
2
2
2
+
2
+
du dv dw dw dy dx dx dy
dx dy J + E dt2 E q.
Заметим, что система нелинейных интегро -дифференциальных уравнений движения вязкоупругой оболочки (7) является достаточно общей, из которой в частном случае можно получить уравнения движения вязкоупругих пластин, пологих и цилиндрических оболочек постоянной и переменной толщины.
Таким образом, математическая модель в геометрически нелинейной постановке задач динамики вязкоупругих изотропных пластин, панелей и оболочек с переменной жесткостью описывается, ИДУ в частных производных вида (7), при соответствующих начальных и граничных условиях.
Выводы. На основе интегральных моделей разработаны обобщенные математические модели нелинейных задачи динамики вязкоупругих изотропных пластин, панелей и оболочек с переменной жесткостью. Разработаны и развиты математические модели задач о нелинейных колебаниях и динамической устойчивости пластин, круговой цилиндрической панели, пологой оболочки двоякой кривизны и цилиндрических оболочек с симметричной и несимметричной толщиной.
Список литературы. 1. Верлань А.Ф. Математическое моделирование флаттера летательных аппаратов. Монография // А.Ф. Верлань, Б.А. Худаяров. - Из-во "Palmarium Academic Publishing". - Saarbrucken, Germany, 2012. - 310 c. 2. Верлань А.Ф. Компьютерное моделирование флаттера вязкоупругих ортотропных пластин в сверхзвуковом потоке газа /А.Ф. Верлань, Б.А. Худаяров, Э.Ф. Файзибоев, З.У. Юлдашев // Вестник НТУ "ХПИ". -2012. - № 62 (968). - С. 8-17. 3. Верлань А.Ф. Численное моделирование нелинейшх задач динамики вязкоупругих систем с переменной жесткостью / А.Ф. Верлань, Р.А. Абдикаримов, Х. Эшматов // Электронное моделирование. - 2010. - Т. 32. - N° 2. - С. 3-14. 4. ЭшматовХ. Клмпьютерное моделирование задачи о нелинейном флаттере вязкоупругой пластины из композиционого материала с сосредоточенными массами / Х. Эшматов, Д.А. Ходжаев, Б.Х. Эшматов, О.Р. Кучаров // Электронное моделирование. - 2010. - Т. 32. - № 5. - С. 310. 5. Григолюк Э.И. Устойчивость оболочек / Э.И. Григолюк, В.В. Кабанов. - М.: Наука, 1978. - 360 с. 6. Reissner E. On transverse vibration of thin shallow shells / E. Reissne // J. Quart. Appl. Math. - 1955. - Vol. 13. - № 2. - P. 169-180. 7. ИльюшинА.А. Основы математической теории термовязкоупругости // А.А. Ильюшин, Б.Е. Победря. - М.: Наука, 1970. - 280 с. 8. Вольмир А.С. Нелинейная динамика пластинок и оболочек / А.С. Вольмир. - М.: Наука. -1972. - 432 с.
Поступила в редакцию 08.08.2013
УДК 539.3
Математична модель задачі динаміки в'язкопружних систем із змінною жорсткістю / Худаяров Б.А., Абдикаримов Р. // Вісник НТУ "ХПІ". Серія: Інформатика та моделювання. - Харків: НТУ "ХПІ". - 2013. - № 19 (992). - С. 160 - 166.
У роботі на основі інтегральних моделей розроблені і розвинені розрахункові математичні моделі для дослідження нелінійних коливань в'язкопружних пластин, кругової циліндрової панелі і циліндрової оболонки із змінною жорсткістю. Іл.: 1. Бібліогр.: 8 назв.
Ключові слова: математична модель, змінна жорсткість, в'язкопружні системи, нелінійні коливання, задачі динаміки.
UDC 539.3
Mathematical model problem of dynamics viscoelastic systems with variable rigidity / Khudayarov B.A., Abdikarimov R. // Herald of the National Technical University "KhPI". Subject issue: Information Science and Modelling. - Kharkov: NTU "KhPI". - 2013. - N° 19 (992). - P. 160 - 166.
In the work on the basis of integral models developed mathematical model for research of nonlinear vibrations of viscoelastic plates, circular cylindrical panel and cylindrical shell with variable rigidity. Figs.: 1. Refs.: 8 titles.
Keywords: mathematical model, variable rigidity, viscoelastic systems, nonlinear oscillations, problem of dynamics.