CC BY
36
УДК 539.3
ДОСЛІДЖЕННЯ РВНИХ МОДEЛEЙ ДОPОЖНЬОЇ КОНCTPУКЦIЇ ЗА НЕСТАЦІОНАРНОГО НАВАНТАЖЕННЯ
О.В. Воропай, доцент, к. т.н., С.І. Новаляєв, доцент, к.т.н.,
В.Т. Г рншакін, асистент, к.т.н., ХНАДУ
Анотація. При моделюванні було використано теорію тонких пластин Кірхгофа та теорію пружних пластин середньої товщини типу Тимошенка. Розглянуто можливість урахування в ’язкопружних властивостей матеріалу пластини. Отримано аналітичні залежності для параметрів напружено-деформованого стану дорожньої конструкції та проведено числовий експеримент.
Ключові слова: дорожня конструкція, моделювання, пластина, основа, коливання, переміщення, деформація.
ИССЛЕДОВАНИЕ РАЗЛИЧНЫХ МОДЕЛЕЙ ДОРОЖНОЙ КОНСТРУКЦИИ НРИ НЕСТАЦИОНАРНОМ НАГРУЖЕНИИ
А.В. Воропай, доцент, к.т.н., С.И. Новаляев, доцент, к.т.н.,
В.Т. Г ришакин, ассистент, к.т.н., ХНАДУ
Аннотация. При моделировании были использованы теория тонких пластин Кирхгофа и теория упругих пластин средней толщины типа Тимошенко. Рассмотрена возможность учета вязкоупругих свойств материала пластины. Получены аналитические зависимости для параметров напряженно-деформированного состояния дорожной конструкции и проведен численный эксперимент.
Ключевые слова: дорожная конструкция, моделирование, пластина, основание, колебания, перемещения, деформация.
INVESTIGATION OF VARIOUS MODELS OF ROAD STRUCTURES UNDER NON-STATIONARY LOADING
A. Voropay, Associate Professor, Candidate of Technical Science,
S. Povalyaev, Associate Professor, Candidate of Technical Science,
V. Grishakin, Assistant Professor, Candidate of Technical Science, KhNAHU
Abstract. The theory of thin Kirhgoffplates and the Timoshenko’s refined theory of medium thickness plates were used. The possibility of using plate material viscoelastic properties were taken into account. Analytical dependences for road structure stress-strain state parameters were received and a numerical experiment was conducted.
Key words: road construction, modeling, plate, base, vibration, displacement, deformation.
Вступ
Для дослідження напружено-деформованого стану та відновлення навантажень, що діють на дорожню конструкцію, а також ідентифікації її параметрів необхідно, щоб математи-
чна модель достатньо точно відображала нестаціонарні деформаційні процеси, що відбуваються в дорожній конструкції. Особливу увагу цьому необхідно приділити під час розв’язання так званих обернених задач механіки, коли за заданими переміщеннями або
деформаціями дорожнього одягу слід встановити навантаження, що спричинило деформування, або визначити невідомі параметри дорожньої конструкції.
Аналіз публікацій
З огляду на те, що моделювання дорожньої конструкції є актуальною проблемою прикладного характеру, їй присвячено достатньо багато наукових праць, зокрема [1-3]. Проте уточнення моделей і розробка моделей для розв’язання обернених задач є важливим напрямом роботи. У статті [4] було підкреслено, що для розв’язання саме цих задач однією з найбільш перспективних, з погляду простоти за достатньої адекватності опису властивостей дорожньої конструкції, є модель динамічного навантаження пружної ізотропної пластини, що лежить на пружній інерційній основі типу Власова-Леонтьева [5].
де Н - товщина шару основи, що стискається; у - об’ємна вага матеріалу основи; R(x,y,t) -реакція пружної основи.
Відповідно до методики, запропонованої в монографії [5], величина реакції визначається у формі
R( х, у, t) = -2tf А^(х, у, t) +
+ k г • ^(х, у, t) +
т
де tf =
ї? п г = , 0 , [ у2 (z) dz ;
7 4 (1 + Уо )^и ’
00
Е0 =
Ет
гр
V
гр
0 і 2 ’ 0 л 1 - V 1 - V
гр гр
(2)
Мета і постановка задачі
Метою запропонованої роботи є вибір моделі дорожньої конструкції, яка б найбільш адекватно відображала деформаційні процеси, що відбуваються в дорожній конструкції. Розглянуто декілька моделей дорожньої конструкції різного рівня складності, отримано чисельні результати розрахунку параметрів напружено-деформованого стану і проведено їх аналіз. Дорожню конструкцію було представлено у вигляді пластини на пружній основі, а при моделюванні самої пластини було використано 2 різні теорії:
- тонких пластин Кірхгофа з урахуванням в’язкопружних властивостей матеріалу пластини (дорожнього покриття) та без їх урахування;
- теорію пружних пластин середньої товщини типу Тимошенка.
Моделювання основи
Контакт між пластиною й основою приймемо двостороннім. Вертикальні переміщення основи описуються функцією у(г), тобто w0(x,y,z,t)=w(x,y,t)■^(z). При цьому функція розподілу напружень у матеріалі основи, відповідно до матеріалу [5], може мати як лінійний характер, так і нелінійний, наприклад
у( z) =
sh (у(Н - z)) sh (уН )
(1)
Коефіцієнти реакції основи на пластину визначалися за такими залежностями [5]
Е0 У
2т (1 -;2 )^
=5,793-107
0
Па
м
sinh (у*) cosh (у*)
* \ *
У +У
sinh (у*)
tf =
Е0т
8^ (у*) cosh (у*)
у -у
(3)
= 2,844 •107,Па • м;
т/ =Р0
sinh (у*)cosh (у*) - у
**
з КГ
= 1,274-103 —;
м
де у =
н у
т
Модель пластини Кірхгофа
Розглянемо тонку пластину (Кірхгофа), що лежить на пружній одношаровій інерційній основі типу Власова-Леонтьева [5]. Рівняння, що описує процес коливання такої пластини, наведено в роботах [4, 6].
*
Модель, що враховує дисипацію енергії коливань пластини, описується диференційним рівнянням (4)
і д ] д w(x,y,t)
1 + а— I DЛЛw(x, у, t) + ph----------2-----
dt) У У dt2
= P( x, y, t) - R( x, y, t),
(4)
де а - час релаксації деформацій.
Відзначимо, що ця модель коливань прямокутної пластини, з урахуванням в’язко-пружних властивостей, базується на гіпотезі Кельвіна-Фойгта.
Розв’язок рівняння (4) було знайдено у вигляді подвійного ряду Фур’є
/ Ч V- V- / ч • knx . nny ...
w(x,y,t) = Ееwkn(t)sin - sin--------. (5)
k n L m
Розв’язок рівняння (4) відносно прогинів пластини було отримано з використанням перетворення Лапласа. Власні частоти коливань (згідно з моделлю Кірхгофа - одномо-дова апроксимація) можуть бути визначені у вигляді[7]
4
wkn (t) =—(----------Vх
lm I mj + mf I xj P(t) sin (Pk t) sin (ц j y0) e~nkn (t~t) х
sin pdk (t -т) d p 0
:---------k-----------d т, если pл > 0;
pi °kn
P dk
wkn (t) =-
(6)
lm (ml + mf ) xj Р(т) sin (Pk т) sin (ц j yo) e~nkn (t-т) x
shpdk (t -т) p *
dт, если pd < 0,
де pdkn =Vpln-nln; ml =p-h;
т П
Fkn(s):^ P(t )sin я kVot sin ц,У0 ; ц n =—;
Модель пластини середньої товщини
Систему рівнянь, яка описує процес коливання прямокутної пластини середньої товщини, що лежить на основі, наведено в роботі [4]. Розв’язок цієї системи у припущенні нульових початкових умов для шарнірно-обпертої пластини на пружній основі отримано у вигляді[4]
wo (x, У, t) = jo Р(т)KW (x, у, t - z)dт; у x (x, у, t) = |0Р(т) K ^ (x, y, t - x)d т; У у (x, у, t) = j0 POOK^ (x, y, t -z)d т,
(7)
де KW (x, y, t), K(x, y, t), KTy (x, y, t) - ядра відповідних інтегралів Дюамеля (згорток) виду
Cf ■ k nx .
nny
m
K (x,y,t) = ЕЕл^sin l sin
k=1 n=l \n '
x[^iftn®lfkn •sin ®fknt - •sin ft ];
Cibl?, knx .
nny
k=l n=l л kfn l
m
sin to{kJ sin ft
to
l kn
to
2 kn
K - (x, У, t) = Е Е s.. k n cos nny
k=1 n=l Лг.
kn
-sin-----------cos-
l
m
sin ®ifknt sin ft
to
l kn
to
2 kn
22
Тут я 2n =n2(-^ + ^r); я k = n^; ц n = n- ;
22
lm
m
G h 2tf kf
a =-----------; a =--------—; a =
1 ph + mf’ 2 ph + mf’ 3 ph + mf
p-1 p- I
. k n aD
яk =~t ; 2nkn =-------------------
l ml + mf
& =
D(я2+ц2) - 2tf (я2+цП )+kf
ml + mf
Л(п =
(я2,(al + a2 + d) + a3 + b)2 -
-4
^ (яld (al + a2 ) + a3d + a2b ) + a3b
tolkn = -J°.5 [(^n (al + a2 + d) + a3 + b) + ЛЬ! ] ;
X
x
X
4
x
ш
1ІП
ш2іп = «/0.5
д/05 [(ЯЬ! (°1 + °2 + Л) + а3 + Ь)-АІП ] ;
делюванні контактних навантажень може бути враховано, що пляма контакту автомобільної шини з дорожнім полотном має круглу або еліптичну форму (рис. 2).
ш
2ІП
Причому шІ, ю2уы - власні частоти коливань (згідно з моделлю типу Тимошенка -двомодова апроксимація).
Вкажемо, що на основі залежностей (7) нескладно отримати вирази для напружень і деформацій у довільній точці пластини.
Для деформацій аналітичні вирази можна
дУх :
дх
одержати на основі формул: вх = z ■
с>Уу . .
в у = z------з використанням співвідношень
су
для гх та &у із (7), а саме
І
вх (X, у, і) = |Рг (т)К“ (х, у, і - т^т;
0
в у (X, у, І) = | Р2 (т) К в (X, у, І - т^ т,
(8)
ї- ад ад
де К “( х, у, t) = - - ХЕ
2 k=1 П=1
СП • ь •яІ2
А і,
. і л • х . пл • у
Х8ІП----------8ІП----------
т
8іпш(и 8іпш2п
ш
1ІП
ш
2ІП
— ад ад
Кєу (х, у, t) = -- ЕХ
2 І=1 п=1
. Іл • х . пл • у
Х8ІП---------------8ІП-------------
Сі • ь у/
і
т
"■ІП
----ї-------- л
А 7
^ІП
8ІП ш(и 8ІП ш2п
ш
1ІП
ш
2ІП
Слід зазначити, що в ядрах виразів (7) та (8) коефіцієнти С[п є коефіцієнтами розкладання функції зовнішнього навантаження в ряди Фур’є.
У випадку випробувань дорожньої конструкції рухомими навантаженнями (рис. 1, 2) в
коефіцієнтах С[п потрібно замінити параметр, що відповідає координаті х, на У0 ■ і , де V - швидкість руху навантаження. При мо-
щ ш г
Рис. 1. Схема навантаження колесом
Рис. 2. Схема пластини, що лежить на основі, для моделювання дорожньої конструкції
Наприклад, у випадку навантаження дороги колесом (кругла рухома пляма контакту з рівномірним розподілом тиску по області)
С[п можна записати у вигляді
С7 =
ІП
Х8ІП
4
1
2л • г
1 •т р •— + т0 Яіп
П%- у0
) Х
8ІП
т
(9)
Числовий розрахунок
Далі наведено розрахунки напружено-деформованого стану для розглянутої моделі дорожньої конструкції за таких вихідних даних:
постійні дорожнього покриття: щільність р=2250 кг/м3; модуль пружності Е=3,2-109 Па; коефіцієнт Пуассона у=0,3; товщина пластини (дорожнього покриття) й=0,1 м.
Х
Постійні пружної основи: щільність р0= = 2000 кг/м3; модуль пружності Е0=1,53-Ш8 Па; коефіцієнт Пуассона у0=0; глибина основи (нормальні переміщення цих точок вважалися рівними 0) Н=1,5 м.
Довжина пластини (досліджуваної ділянки дороги) /=10 м. Ширина пластини т=3,5 м. Координати лінії траєкторії рухомого навантаження у0= 1,75 м. Координати точки, в якій обчислювалися прогини і деформації пластини, а також точки прикладення нерухомого навантаження х= м, уж=1,75 м (геометричний центр пластини). Швидкість рухомого навантаження У=25 м/с.
Величина рухомого навантаження
Р=32000 Н; при розрахунках з рівномірно розподіленим по колу навантаженням
^=352453,85 Па радіус плями контакту колеса з дорогою г=0,17 м. Число кроків за часом - 500. Розглянутий проміжок часу 7=0,4 с.
На рис. 3 представлено зіставлення прогинів пластини при дії рухомого навантаження на дорожню конструкцію (механічну систему, що складається із пластини і пружної основи), які було отримано для розглянутих моделей.
На рис. 3 крива 1 - розрахована на основі теорії тонких пластин Кірхгофа; крива 2 -розрахована на основі теорії тонких пластин Кірхгофа з урахуванням дисипації (в’язко-пружних властивостей матеріалу пластини -дорожнього покриття); крива 3 - розрахована на основі уточненої теорії пластин середньої товщини типу Тимошенка. Відзначимо, що вертикальна вісь Oz, вздовж якої відкладено прогини пластини, спрямована вниз.
Однією з переваг моделювання дорожньої конструкції у вигляді саме пластини середньої товщини, що лежить на основі, є можливість легко отримати аналітичні вирази для деформацій та напружень.
На рис. 4 представлено результати розрахунку деформацій пластини: крива 1 - поздовжня деформація вх (і) , крива 2 - поперечна деформація ву (і) .
На рис. 5 представлено двовимірну епюру прогину пластини у момент часу і=0,2 с, що відповідає знаходженню рухомого навантаження точно посередині пластини.
Рис. 3. Зіставлення прогинів пластини за різними моделями
Час t, с
Рис. 4. Результати розрахунку деформацій пластини
Рис. 5. Епюра прогину пластини
Висновки
На основі отриманих у роботі аналітико-чисельних залежностей та чисельних розрахунків, виконаних для різних типів динамічних моделей дорожньої конструкції, результати яких у графічному вигляді наведено на рис. 3-5, можна зробити такі висновки:
1. Модель дорожньої конструкції, що складається з тонкої пластини теорії типу Кірхгофа та пружної інерційної основи є більш простою, але дозволяє враховувати дисипацію енергії коливань, тобто в’язкопружні властивості матеріалу покриття без особливих ускладнень.
2. Модель, яка базується на уточненій теорії типу Тимошенка, більш точно моделює нестаціонарні деформаційні процеси для пластин середньої товщини (оскільки є двомодо-вою апроксимацією та дозволяє отримати розв’язки, близькі до тривимірної теорії пружності). Прогин та дотичні переміщення в теорії типу Тимошенка (а тому і деформації) є незалежними змінними. Урахування в’язкопружних властивостей в рамках цієї теорії вносить значні ускладнення математичного характеру без значного уточнення результатів.
3. Числові результати розрахунків для обох моделей, які було розглянуто, добре узгоджуються.
Література
1. Углова Е.В. Моделирование деформиро-
вания нежестких дорожных конструкций при воздействии движущегося транспорта / Е.В. Углова // Изв. вузов. Сер. Строительство. - 2009. - № 3-4. -С.31-35.
2. Матвиенко Ф.В. Моделирование деформа-
ции дорожной конструкции от воздействия транспортного потока / Ф.В. Матвиенко, А.Н. Канищев, В.В. Волков // Научный вестник ВГАСУ. Строительство и архитектура. - 2009. - № 3. -С.67-72.
3. Смирнов А.В. Теоретические и экспери-
ментальные исследования работоспособности нежёстких дорожных одежд: автореф. дис. на соискание ученой степени доктора техн. наук: спец. 05.23.11 «Проектирование и строительство дорог, метрополитенов, аэродромов, мостов и транспортных тоннелей» /
А.В. Смирнов. - М., 1991. - 38 с.
4. Вибір моделі дорожньої конструкції для
розв’язку задач ідентифікації за нестаціонарного навантаження / Є.Г. Янютін, Д.І. Богдан, О.В. Воропай, С.І. Поваляєв // Автомобільний транспорт. - 2010. -№27. - С. 148-152.
5. Власов В.З. Балки, плиты и оболочки на уп-
ругом основании / В.З. Власов, НИ. Леонтьев. - М.: Физматгиз, 1960. - 492 с.
6. Пановко Я.Г. Внутреннее трение при коле-
баниях упругих систем / Я.Г. Пановко. -М.: Наука, 1960. - 186 с.
7. Гринев В.Б. Идентификация подвижного
нагружения, воздействующего на вязкоупругую пластину на упругом основании / В.Б. Гринев, Е.Г. Янютин,
В.Т. Гришакин // Вісник НТУ «ХПІ». -2011. - №13. - С. 69-74.
Рецензент: В.К. Жданюк, професор, д.т.н., ХНАДУ.
Стаття надійшла до редакції 15 грудня 2011 р.
