Динамика конструкций и сооружений
СВОБОДНОЕ КОЛЕБАНИЕ НЕОДНОРОДНОГО РАЗНОМОДУЛЬНОГО СТЕРЖНЯ, ЛЕЖАЩЕГО НА ДВУХКОНСТАНТНОМ ОСНОВАНИИ
НАТИГ С. РЗАЕВ
Институт Математики и Механики НАН Азербайджана, ул. Б Вагабзаде 9, Баку, AZ1141, Азербайджан, natiq.rzayev.1984@list.ru
В работе рассматривается задача о собственном колебании стержня из неоднородного по толщине и длине разномодульного материала и лежащего на двухкон-стантном основании типа Пастернака П.Л.
Уравнение движения является сложным дифференциальным уравнением с частными производными четвертого порядка. Решение задачи строится с помощью метода разделения переменных и метода ортогонализации Бубнова-Галеркина. При конкретных значениях характерных параметров проведен расчет, и результаты представлены в виде таблиц и графиками зависимости между круговой частотой и параметрами неоднородности. Результаты расчетов показывают, что учет неоднородности модуля упругости и плотности существенным образом влияет на величины круговых частот.
КЛЮЧЕВЫЕ СЛОВА: растяжение, сжатие, изгиб, кручение, эластичность, круговая частота, колебание.
В работе исследуется задача свободного поперечного колебания стержня прямоугольного поперечного сечения, изготовленного из разномодульного неоднородного материала, лежащего на двухконстантном основании Пастернака [1]. Решение задачи осуществляется путем применения метода разделения переменных Бубнова- Галеркина. В первом приближении найдена формула для определения значения круговой частоты. При конкретных значениях характерных параметров проведен численный анализ, и результаты представлены в виде таблиц и графиков.
Как известно, в строительстве современных инженерных сооружений, в машиностроении и во многих других отраслях широко используются конструкции из неоднородного разномодульного материала. Сюда можно отнести, в основном, композитные неоднородные материалы, наполненные полимерами [1,6]. Предполагается, что модуль упругости соответственно при растяжении и сжатии и плотность являются непрерывными функциями координат длины х и высоты z:
E + = Eo f (x )л (z); E - = e0 f (x )л (z); p = po¥i (x У 2 (z). (1)
Здесь E0+, E-, p0 - соответствуют однородному изотропному материалу, f (x) со своими производными до второго порядка и y1 (x)- являются непрерывными функциями. Отметим, что специфическим фактором в данном случае является то, что граница растянутых и сжатых областей существенным образом зависит от функции f2 (z).
Распределение напряжений по сечению определяется с использованием следующих соотношений:
ff = E+fi (x)f2 (z)(l -pz); ff > 0; ff = E f (xf (z)(l -pz); ff < 0. (2)
Предполагается, что уравнение нейтральной линии определяется следующим соотношениям:
l - z-p = 0, (3)
где I и значения деформации и кривизны центральной лини, г0 - граница нейтральной линии.
Сюда необходимо добавить условие отсутствия и уравнение момента:
+ h
^аёг = ёг + ёг = 0
+щ
М = |гаёг = |га +ёг + |го~ёг = 0 . -h h1 h2 Здесь Н1 и Н2 - размеры растянутых и сжатых зон.
Учитывая формулы (1) и (4), получим
+ h
(4)
(5)
| аёг = | /2 (г% - г р)ёг + а| f2 (г% - г р)ёг = 0 , [а = Е+ ^;
| /2 (г )гёг + а | г/2 (г ё
отсюда находим:
I = р~
| /2 (г )ёг + а\ /2 (г )ёг
(6)
Теперь определим изгибающий момент:
| Е+/2 (г + г (р)гёг + | Е - /2 (г + г (р)гёг
М=/■(•)
(7)
или же
М
ЬЕ+/1 (х)
= I
| г/2 (г + а | г/2 (г
+ Р
| г / (г )ёг + а | г / (г )ёг
щ
Подставляя сюда выражение (6), получим:
7 V
| г/2 (г )ёг + а | г/2 (г )ёг
М
ЬЕ+/1 (х)
| /2 (г )ёг + а\ /2 (г )ёг
+ | г2 / (г )ёг + а| г2 / (г )ёг
где, как уже указывалось ранее, Щ и к2 - размеры по высоте растянутых и сжатых зон.
С учетом (7) и сопротивления двухконстантного основания уравнение движения записывается в следующем виде:
а2
дх2
/ (х )
дЖ
дх2
— пг — д Ж _ / ча ж + КЖ - К + Р¥(х )—2- = 0.
д 2Ж
^ дх2
дt2
(8)
Здесь К1 и К2 - характеристики основания, определяемые с помощью экспериментов. Приняты следующие обозначения:
К = К (М0К)-1; К2 = (М0К)-1; Р0 = тЩ ]>2 (г)ёг; Р = Р0 (М0К)-1.
Уравнение (8) можно записать в следующем виде:
2
2
2
2
2
2
Э \д ^ (г- V ( \\д
+ 2 /1+ (/1 - к 2 /1 ))—г~+
~ 4 ^ 1 V / ~ 3 V 1 2 л 1 V // ~ 2
дх дх дх
+ к • /!-1 (х)Ж + р¥2 (х/Г1 (х)^ = 0 . (9)
дt
Как видно, уравнение (9) является сложным и поэтому при решении следует применить приближенные аналитические методы. В данном случае будут использованы методы, которые являются эффективными и апробированными, а именно, метод разделения переменных и метод Бубнова - Галеркина.
В первом случае решение будем искать в следующем виде:
Ж(х) = V (х)вш М, (10)
где V (х) должно удовлетворять краевым условиям, а — круговая частота. Подставляя (10) в (9), получим:
а V ) Л У ( - —1 ( )) ё V -
—+2/;{х) •—+(//—к2/11 (х))—-ах ах ах
Решение уравнения (11) будем искать с помощью метода Бубнова- Галеркина, а функцию V (х) примем в следующем виде:
+ 2//(х>^ + (/-— К2/— (х))^ + К • /1—1 (х)У — рм2щ(х)/1—1 (х) = 0. ( 11)
а,
1=1
У(х) = 2 аА (х). (12)
г=1
Здесь а^ — неизвестные постоянные и каждый член ряда А (х) должен удовлетворять соответствующим краевым условиям.
Функция ошибки в данном случае записывается в следующем виде:
/1 (х) + 2 /1(х) "А + (/1'(х) — К2/— (х 0 + ах ах ах
+ К/Г1 (хУ — Ых)/Г1 (х)(2 ]* 0 . (13)
Условие ортогонализации в данном случае имеет вид:
I
= 0 k = 1,2,... (14)
0
В произвольном приближении а2 определяется из системы уравнений
(14). Эта система относительно а2 является линейно однородным уравнением. Для существования нетривиального решения главный определитель данной системы должен обращаться в нуль:
||(2|| = 0. (15)
При раскрытии (15) получается алгебраическое уравнение «-ой степени относительно а2. Хотя при помощи современной компьютерной техники определение а2 при любом приближении не вызывает особого труда, будем определять основной тон частоты. Это соответствует первому приближению:
I
|^(х АА = 0, (16)
/
/, (х) ^ + 2 /;(х) + /(х)— К2 /— (х)) +
ах ах ах
+к • /г1(х)А1 —ММх)Г1(хА Аа=0. (17)
Отсюда находим:
ю2 =
¡
/, (x) Л-Л. + 2f;(x) d-6 + (/;« - K /- rf6 + K, • /г' «6,
dx
dx
dx
6,dx
(18)
P¡v(x)• /- (x,
x К dx
В случае если стержень является неоднородным только по модулю упругости ц/(х ) = 1 и формула (18) принимает следующий вид:
со} =
/1 (x)^ + 2/,'(x)d6, + fi'(x) - K2/,-1 (x))^ + K, • /1-1 (xК dx dx dx
6 , dx
(19)
^ / - (x 6,2 dx
Из выражений (18) и (19) получим следующие соотношения:
со} =
f 2 с
V Ю J
i
i /- (x 6,2 (x d 0_
l
i /1-1 (x )y(x662 (x)dx
(20)
Для анализа рассмотрим случай: 6, (x)= sin nx; /, (x) = 1 + sx; у, (x) = 1 + jx; /2 (x) = 1 + sz; y/2 (z) = 1 + jz . (21)
Здесь x = x • l_1; z = z • h_1; se [0,1]; p e [0,1].
С учетом (21) формула (20) принимает следующий вид:
(„ V 1
с2 =
ю
VJ
i 1
J_ + sx
- sin2 nxdx
i 1
1 + sin2 nxdx.
+ sx
(22)
0 / 0 Для простоты анализа примем следующие обозначения:
1 • 2 - 1 _ . 2 -
, fSin nx , , ex Sin nx ,
b = -dx; b = -dx;
1 J 1 + sx 2 •> 10 0
+ £х •'1 + £х
0
тогда формула (22) записывается в следующем виде:
ю2 =■
b,
1
b1 + Jb2
1 Ь2
1 + J — b
> 0.
л Ь2 Ь2 1 Ь2 1 Отсюда получим: 1 + ц— > 0 т.е: ц— > -1; —>--.
Ь1 Ь1 Ь1 М
Продолжая анализ результата, примем что сопротивление внешней среды отсутствует. Тогда формула (18) принимает следующий вид:
i
/1 (x) ^ + 2 /;(*) + fHx) ^
dx
dx
dx2
6,dx
,, 2 0
®0 =-}-.
Р\у(х)■ /- (х№ёх
0
Из формул (18) и (23) можно получить следующее соотношение:
со,2 =
С С 2
с
Vсо J
\
= 1 + -
KJ- (xЪ - K2/- (x)^
dx
6Adx
i i
/ (x) ^ + 2 f(x) ^ + ^
dx 4 dx 3 dx 2
UA l/fA. UA
(24)
dxdx
Учитывая (21), получим:
с2=1 + !
— sin ш — 1 (ж
K1
1 + sx
+ K
1 + sx V l
sin 7tdx
f / \2 A1 • 2
— — ( ж 1 TSin 77 K1 + K 2
(1 + sx)| 7 1 - 2s( 7 I cos 77
= 1 + -
If 7
+ sx
Sin 77
ж 141
l J 2
с2 = 1 + 2| —
f
K + K2I7J b,.
V J
2
(25)
Расчет можно провести и для случаев:
/ (р) = 1 + seр; (р) = 1 + /ер; /2 (р) = 1 - s cos 27p; \у2 (р) = 1 - и cos 27р . (26)
Результаты проведенных численных расчетов, представлены в виде таблиц и графика:
1-2 1-2
= fsn_7Pdp; b2 = f Psin 7 dp.
1 J 1 + sp 2 J 1
0 ^ 0
+ sp
Таблица 1
Таблица 2
2
е b1 Ь2 е и II 5
0 0,5 0,25 0 0,8
0,25 0,445 0,219 0,25 0,802
0,50 0,402 0,196 0,50 0,804
0,75 0,367 0,177 0,75 0,805
Рис. 1. График зависимости безразмерной частоты от параметра неоднородности
Л и т е р а т у р а
1. Толоконников Л.А. О связи между напряжениями и деформациями в разномо-дульных изотропных средах// Инж. журнал. МТТ. - 1968. - № 6. - С. 108-110.
2. Панферов В.М. Теория упругости и деформационная теория пластичности твердых тел с разными свойствами на сжатие и кручение. - ДАН. СССР. - 1968. - Т. 180. -№ 1. - С. 41-44.
3. Москвитин В.В. Сопротивление вязкоупругих материалов. - М., 1972. - 320 с.
4. Быков Д.Л. Основные уравнения и теоремы для одной модели физически нелинейной среды// Инж. журнал. МТТ. - 1966. - №4. - С. 58-64.
5. Шапиро Г.С. О деформациях тел, обладающими различным сопротивлением растяжению и сжатию// Инженерный журнал. Механика твердого тела. - 1966. - № 2. -С. 123-125.
6. Ломакин Е.В., Работнов Ю.Н. Соотношения теории упругости для изотропного разномодульного тела// Известия АН СССР. МТТ. - 1978. - № 6. - С. 29-34.
7. Пастернак П.Л. Основы нового метода расчета фундаментов на упругом основании при помощи двух коэффициентов постели. - М.: Стройиздат. - 1954. - 89 с.
8. Gadjiev V.D., Rzayev N.S. Sollitions of nonhomegenous viscoelastic foundation// Translation of NAS of Azerbaijan, 2013. - Vol. XXXII, № 4. - Pp-133-138.
9. Новацкий В. Динамика сооружений. - М, 1963. - 376 с.
10. Маркин А.А., Соколова М. Ю., Христич Д.В. Проблемы прочности, пластичности и устойчивости в механике деформируемого твердого тела// Материалы Международного научного симпозиума, Тула: Тул. ГУ. - С. 97-98.
References
1. Tolokonnikov, L.A. (1968). O svyazi mezhdu napryazheniyami i deformaziyami v raznomodul'nyh izotropnyh telah, Inzhenerniy Zhurnal. Mehanika Tverdogo Tela, №6, pp.108-110.
2. Panferov, V.M. (1968). Teoriya uprugosti i deformatzionnaya teoriya plastichnosti tvyordyh tel s raznymi svojstvami na szhatie i kruchenie, DAN SSSR, Vol. 180, №1, .41-44.
3. Moskvitin, V.V. (1972). Soprotivlenie VyazkouprugihMaterialov, Moscow, 320 p.
4. Bykov, D.L. (1966). Osnovnie uravneniya i teoremi dlya odnoj modeli fisicheski nelineynoj sredi, Inzhenerniy Zhurnal. Mehanika Tverdogo Tela, №4, pp. 58-64.
5. Shapiro, G.S. (1966). On deformation of bodies possessing different resistance to tensile and compression, Inzhenerniy Zhurnal. Mehanika Tverdogo Tela, № 2, pp. 123-125.
6. Lomakin, E.V., Rabotnov, Yu.N. (1978). Sootnosheniya teorii uprugosti dlya izotropnogo rasnomodul'nogo tela, Izv. AN SSSR. Mehanika Tverdogo Tela, № 6, pp. 29-34.
7. Pasternak, P.L. (1954). Osnovi Novogo Metoda Rascheta Fundamentov na Uprugom Osnovanii pri Pomoschi Dvuh Koeffitzientov Posteli, Мoscow, Strojizdat, 89 p.
8. Gadjiev, V.D., Rzayev, N.S. (2013). Sollitions of nonhomegenous viscoelastic foundation, Translation of NAS of Azerbaijan, Vol. XXXII, № 4. pp. 133-138.
9. Novatzkiy, V. (1963). Dinamika Sooruzheniy, Moscow, 376 p.
10. Markin, A.A., Sokolova, M.Yu., Hristich, D.V. Problemi prochnosti, plastichnosti, i ustojchivosti v mehanike deformiruemogo tverdogo tela, Materialy Mezhdunarodnogo yauchnogo simposiuma, Tula, TulGU, pp. 97-98.
OSCILLATION OF A HOMOGENEOUS DIFFERENT MODULUS BAR LYING ON TWO CONSTANT FOUNDATION
Natig S. Rzayev Institute of Mathematics and Mechanics of NAS of Azerbaijan, AZ1141, 9, B. Vahabzade str., Baku, Azerbaijan
in the paper, a problem of natural vibration of a bar inhomogeneous in thickness and length made of different modulus material and lying on two Pasternak -type foundation is considered.
The equation of motion is the fourth order complex partial differential equation. The solution of the problem is constructed by means of method of operation of variables and of the Bubnov - Galerkin orthogonalization method. At specific values of characteristic parameters, the calculation in conducted and the results are represented in the form of tables and graphs of dependence between annular frequency and inhomogeneity parameters The results of calculations show that account of inhomogeneity of elasticity modulus density and also mediums resistance essentially influences' on the value of annular frequency.
Key words: tension, compression, bending, tersion, elastic, frequently, vibration.