dira)1
SHSgJ
УДК 624.04
РАСЧЕТ ДЕРЕВЯННОЙ БАЛКИ ПРЯМОУГОЛЬНОГО СЕЧЕНИЯ С УЧЕТОМ ПОЛЗУЧЕСТИ
Чепурненко А. С., Ульянская В. В.
Донской государственный технический университет, Ростов-на-Дону, Российская Федерация
[email protected] [email protected]
В статье приводится вывод разрешающих уравнений для расчета деревянной балки с учетом ползучести. Приведен пример расчета шарнирно-опертой по концам балки прямоугольного поперечного сечения, загруженной равномерно-распределенной нагрузкой. Величина прогиба найдена по методу конечных разностей. Разработана программа для расчета в пакете МЛТЬЛБ с возможностью изменения исходных данных и выводом графика зависимости перемещений от времени. Приведено сравнение максимальной величины прогиба с аналитическим решением. Установлено, что в процессе ползучести напряжения остаются постоянными.
Ключевые слова: ползучесть, метод конечных разностей, уравнение Максвелла-Томпсона, деревянная балка.
UDC 624.04
CREEP CALCULATION OF RECTANGULAR SECTION WOOD BEAM
Chepurnenko A. S., Ulianskaya V. V.
Don State Technical University, Rostov-on-Don, Russian Federtion
anton_chepurnenk@mail .ru [email protected]
The article provides the derivation of the resolving equations for creep calculation of the wood beam. An example of hinged on edges and loaded with a uniformly distributed load rectangular transverse section beam is presented. The deflection has been found by finite differences method. The program has been developed in MATLAB environment with the opportunity of initial data change and derivation of dependency graph of displacement on time. The comparison of deflection maximum with the analytical solution is given. It has been found that stresses remain constant during creep process.
Key-words: creep, finite differences Maxwell-Thompson equation, wood beam.
method,
Введение. Древесина является нелинейным материалом при действии как кратковременных, так и длительных нагрузок, и может быть исследована по времени в двух направлениях: длительной прочности и ползучести. Большинство исследований посвящено первому направлению, в то время как в Еврокодах определение деформаций ползучести является обязательным [1,2]. В современной литературе приводятся некоторые частные задачи расчета с учетом ползучести деревянных конструкций. Так, расчет сжатых деревянных стоек на устойчивость при ползучести рассматривается в работе [3]. В монографии Н. Х. Арутюняна [4] приводятся методы расчета на ползучесть железобетонных балок с возможностью применения к другим материалам, в том числе и к дереву. Однако теория Н. Х. Арутюняна связана с решением интегральных уравнений, что зачастую сопряжено с математическими трудностями. В настоящей статье приводится относительно простой метод расчета деревянных балок с учетом ползучести. При этом закон ползучести может быть произвольным: как линейным, так и нелинейным.
Вывод разрешающих уравнений. Схема деревянной балки прямоугольного сечения с размерами ЬхИ приведена на рис. 1.
Рис. 1. Деревянная балка прямоугольного сечения
Полная деформация в произвольный момент времени определяется как сумма упругой деформации и деформации ползучести е*:
о *
е = —+ е Е
Из гипотезы плоских сечений полная деформация определяется следующим образом:
ё2 ^
(1)
е = -2-
ёх2
Подставим это выражение в (1), предварительно выразив напряжение о через деформации:
ё2 ^ ёх2
Изгибающий момент связан с напряжением интегральным соотношением:
М = -1огёЛ
и ГУ *Л
о = Е(-г—^--е )
(2) (3)
Подставим (2) в (3):
Раскрыв скобки, получим:
М = — Е (- ^-е*) гёЛ; л ёх
М = Е| г2ёЛ ^ + Е| е*гёЛ. л ёх л
Введем обозначение М* = Е| е*гёЛ . Величина | г2ёЛ представляет осевой момент инерции
I = Ькъ /12. Элемент площади ёЛ можно заменить произведением Ь • ёг , тогда получим:
ы2
М = ЕЬ | е гёг
(4)
-ь/2
Окончательно разрешающее уравнение принимает вид:
М = Е1ё-^- + М \ ёх
Методика расчета. Интеграл (4) вычисляется численно по формуле трапеций:
Л = Аг | ^ + ЪУ- + — I I 2 7=1 ' 2 )
где Аг — шаг интегрирования; у0 и уп - начальное и конечное значения функции; у/ -промежуточные значения функции.
Величина прогиба была найдена методом конечных разностей. Для этого балка длиной I была разбита на определенное число интервалов п по длине величиной Ах . Вторая производная в /-ом узле определяется по формуле:
ё2 ы _ Ы/+1 - + Ы/-1
ёх2
Ах2
Решение сводится к системе линейных алгебраических уравнений:
[ 4[ »]=[ г] •
где [А] — матрица коэффициентов при размерностью (п+7; п+7).
"1 0 0 ... 0 0 0
(5)
А =
а2 Ь2 С2 0 а3 Ь3
000 000
0 0 0 0 0 0
ап Ьп Сп 0 0 1
где а = А7; Ь' = А7; с 4?.
[ы] — вектор неизвестных узловых перемещений; [ы] = [
, ^2, ^3,..., Ып+1
]; [ В ] — вектор
свободных членов
[ В] = [0,/1, /2, /з,..,0],
где / = М , Мг = ^ • х-<(0 < х < I) . ' Е1 ' 2 2 4 '
Перемещения найдем, выразив их из формулы (5):
[ ы] = [ А]-1 [ г ].
В качестве закона связи деформаций ползучести и напряжений будем использовать уравнение Максвелла-Томпсона [5]:
дв „ тт да
—Еп + Нв = а + п—,
ы ы
где Н — длительный модуль деформации (Н = 104 МПа), п — время релаксации (п =18 сут), Е —
мгновенный модуль деформации (Е=1,48-104 МПа).
Подставив в закон ползучести вместо в выражение (1), получим:
,1 до дв* Н тт * да
(--+ —)Еп +—а + Н в =а + п— ;
Е & дГ Е Ы
да „ дв* Н тт * да
п--+ Еп--1--а + Н в =а + п —;
дг дг Е дг
Еп дв* = а(1 - Н) - Н в*; дг Е
^ = а(1 - Н) - Н в*]; дг пЕ Е
Для вычисления деформаций ползучести вводится равномерная сетка по времени. Величина е* в момент времени t + At вычисляется на основе метода Эйлера [6]:
■ч+А
* де = е* +-At.
t дt
Аналогичный подход для решения задач ползучести применяется в работах [7-10].
Результаты и обсуждение. Расчет был выполнен в программном комплексе МАТЬАВ. Полученный график зависимости максимального прогиба от времени приведен на рис. 2:
Рис. 2. Зависимость максимального прогиба от времени
Из представленного графика видно, что деформации резко возрастают в начальный период времени, а затем приходят к конечному значению.
Формула для определения максимального прогиба для шарнирно-опертой балки при действии равномерно-распределенной нагрузки выглядит следующим образом:
5д14
н = ■ тах 384Е1
(6)
Для определения прогиба в конце процесса ползучести достаточно заменить мгновенный модуль упругости Е на длительный Н.
При подстановке заданных размеров балки получим:
5 • 2 • 64
н = -
тах
= 0,0506 м.
384-107 • 0,1-0,23 /12 На графике значение нтах совпадает с вычисленным по формуле (6).
Отметим также, что при выбранном законе ползучести напряжения во времени не меняются. Библиографический список
1. Ягнюк, Б. Н. Теоретические основы проектирования деревянных конструкций по нормам европейского союза — Еврокоду 5 / Б. Н. Ягнюк. — Москва-Берлин: Директ-Медиа, 2015. — 140с.
2. ББ БК 1995-1-1:2004+А2:2014/ Еврокод 5. Проектирование деревянных конструкций.
3. Вареник, К. А. Расчет центрально-сжатый деревянных элементов с учетом ползучести: дис. ... канд. тех. Наук / К. А. Вареник. — Великий Новгород, 2015. — 167 с.
4. Арутюнян, Н. Х. Некоторые вопросы теории ползучести / Н. Х. Арутюнян. — 1952. —
324 с.
5. Клименко, Е. С. Устойчивость сжатых неоднородных стержней с учётом физической нелинейности материала: дис. ... канд. тех. наук / Е. С. Клименко. — Ростов-на-Дону, 2011. — 112 с.
6. Турчак, Л. И. Основы численных методов. / Л. И. Турчак, П. В. Плотников. — Москва : Физматлит, 2003. — 304 с.
7. Аваков, А. А. Напряженно-деформированное состояние железобетонной арки с учетом нелинейной ползучести бетона / А. А. Аваков, А. С. Чепурненко, С. Б. Языев // Научно-технический вестник Поволжья. — 2015. — № 1. — С. 27-30.
8. Юхнов, И. В. Продольный изгиб гибкой железобетонной стойки при нелинейной ползучести [электронный ресурс] / Современные проблемы науки и образования // И. В. Юхнов [и др.]. — 2014. — № 5. — Режим доступа : https://www.science-education.ru/ru/article/view?id=14705.
9. Чепурненко, А. С. Расчет внецентренно сжатого железобетонного стержня на ползучесть при различных законах деформирования / А. С. Чепурненко [и др.] // Научное обозрение. — 2014. — № 8. — С. 935-940.
10. Никора, Н. И. Устойчивость полимерного стержня при ползучести с учетом дискретного спектра времен релаксации полимера / Н. И. Никора [и др.] // Научное обозрение. — 2016. — № 4. — С. 40-43.