О методах расчета сжато-изгибаемых конструкций из клееной древесины Текст научной статьи по специальности «Строительство и архитектура»

Интернет-журнал «Науковедение» ISSN 2223-5167 http ://naukovedenie.ru/ Том 8, №3 (2016) http ://naukovedenie. ru/index.php?p=vol8-3 URL статьи: http://naukovedenie.ru/PDF/56TVN316.pdf Статья опубликована 30.05.2016. Ссылка для цитирования этой статьи:

Гранкин К.В., Шмидт А.Б. О методах расчета сжато-изгибаемых конструкций из клееной древесины // Интернет-журнал «НАУКОВЕДЕНИЕ» Том 8, №3 (2016) http://naukovedenie.ru/PDF/56TVN316.pdf (доступ свободный). Загл. с экрана. Яз. рус., англ.

УДК 624.044.2

Гранкин Кирилл Владимирович

ФГБОУ ВПО «Санкт-Петербургский государственный архитектурно-строительный университет»

Россия, Санкт-Петербург1 Аспирант E-mail: gran1256@mail.ru

Шмидт Александр Борисович

ФГБОУ ВПО «Санкт-Петербургский государственный архитектурно-строительный университет»

Россия, Санкт-Петербург Заведующий секции «Деревянных конструкций» Кандидат технических наук, доцент E-mail: ukftpz@lan.spbgasu.ru

О методах расчета сжато-изгибаемых конструкций

из клееной древесины

Аннотация. В данной статье авторами рассмотрены существующие методики расчета сжато-изгибаемых прямолинейных стержней из клееной древесины. Для каждого из рассматриваемых методов расчета сжато-изгибаемых конструкций, в том числе по действующим нормативным документам на территории Российской Федерации (СП 64.13330.2011) и Европы (EuroCode 5), выявлены недостатки и установлены рамки применяемости каждого из методов расчета. Авторами статьи выполнено сравнение результатов, полученных при решении задачи с помощью имеющихся аналитических способов, а также с применение современных программных комплексов, построенных на методе конечных элементов, таких как SCAD Office и Lira Soft. Также в данной статье авторы выполняют численный эксперимент по всем рассматриваемым методикам, для прямолинейной сжато-изгибаемой балки пролетом 7 метров. В результате численного эксперимента строится график зависимости изменения коэффициента учета деформированного состояния, а также график зависимости изменения максимальных краевых нормальных напряжений от гибкости рассматриваемого элемента. В конце статьи авторами данной статьи приводятся выводы на основании сравнения полученных данных в результате проведенного численного эксперимента, для различных методов расчета сжато-изгибаемых элементов.

Ключевые слова: клееные деревянные конструкции; сжато-изгибаемые конструкции; геометрическая нелинейность; формула Максвелла-Мора; метод перемещений; метод начальных параметров

1 190005, Россия, Санкт-Петербург, 2-я Красноармейская ул., д. 4 1

На сегодняшний день имеется множество научных трудов по расчету сжатых и сжато-изгибаемых прямолинейных стержней, работающих в упругой стадии, а также за ее пределами, с учетом геометрической нелинейности, ползучести, а также с учетом анизотропии свойств древесины.

Современные нормы предлагают выполнять расчеты таких конструкций на основании методов решений, разработанных для прямолинейного стержня, опертого на двух опорах с синусоидальной нагрузкой (см. рис. 1). В основе этих формул кривая изгиба представлена в виде быстросходящегося тригонометрического ряда:

.. . пх г 2пх г . 3пх

у = f1smY +f2sm—f3sm— +

(1)

Рисунок 1. Геометрическая интерпретация разложения в тригонометрический ряд

^т(пжх/1).

В реальных же конструкциях равномерно распределенная нагрузка практически всегда распределена не по закону синуса и такое допущение может приводить к неправильным результатам.

Известно, что для конструкций подверженных сжатию с изгибом, принцип суперпозиции не применим, так как от продольной силы появляется дополнительный изгибающий момент, что в свою очередь приводит к дополнительным деформациям (см. рис. 2).

Рисунок 2. Геометрическая схема деформаций прямого стержня при продольно-поперечном изгибе

Согласно первой методики расчета сжато-изгибаемых элементов, предложенной в тридцатых годах прошлого столетия, краевые напряжения находятся по формуле Ф.С. Ясинского:

О =— + - = -(— + -);

фЛ щ Р \фт р)

(2)

где <рт - коэффициент центрального продольного изгиба в плоскости действия момента.

Неточность данной формулы обуславливается тем, что в ней воздействие момента накладывается на напряжения потери устойчивости центрального продольного изгиба. На эту неточность обращали внимание авторы [1], [2]. При этом, решение данной задачи в [1] и [2] предлагается выполнять с учетом дополнительного изгибающего момента полученного от

сжимающей нагрузки в деформированном состоянии системы. Краевые напряжения в таком случае находятся по формуле сложного сопротивления:

а = Е + м + т+о

А Ш Ш у '

где: Го - прогиб от поперечной нагрузки; Гц - дополнительный прогиб от изгибающего момента.

Формула (3) и ее решение выглядит в [1], [2] по-разному. Так в [1] она представлена в следующем виде (вторая методика расчета сжато-изгибаемых конструкций):

а = 1 + м+мг

А Ш Ш у '

п _М

где Г = „„ N

таким образом (3) приобретает вид:

№ , м

° = т +--(5)

Величину 1 —тКр = % преобразуем, заменив Мэкр = РбрКсфэ, где

3100

<Рэ=— (6),

Таким образом, окончательно получаем формулу следующего вида:

№ , М

в = — +—т; (7)

N

где - - -

* = (1--—) (8)

ъ V фЯсРбр) 4 '

При этом тот факт, что при Х<55 коэффициент ф>1, никаким образом не комментируется с точки зрения устойчивости как физического явления (см. рис. 3). На этот факт также обращали внимание Е.Н. Серов [6], В.З. Клименко [7].

Рисунок 3. График изменения коэффициента продольного изгиба от гибкости

Третья методика расчета сжато-изгибаемых конструкций изложена в [2]. Согласно ей, максимальные краевые напряжения находятся по формуле сложного сопротивления в следующем виде:

N М АМ

а = — +---1--;

А Ш Ш

(9)

где АМ = - увеличение момента от прогиба стержня, где прогиб от действия продольной силы и момента, определяется по формуле:

То

г =

(1-

№

М12

(10)

Здесь: /0 - прогиб от действия момента /0 = —,

аЕ1

Рэ - критическая сила Эйлера, равная

а - коэффициент, зависящий от вида приложения нагрузки.

Основным недостатком данной формулы является то, что она не является универсальной при расчетах сжато-изгибаемых конструкций, ввиду сложности использования коэффициента а, который как известно напрямую зависит от вида приложения нагрузки и представлен в [2] только для прямых стержней.

Четвертый метод расчета был предложен Е.Н. Серовым [3], который предложил использовать дифференциальное уравнение изогнутой оси для прямолинейного стержня в виде:

п й2Мх

" . 7 7 Мх ,

у + к у = — + ■ ,

у у Е1 вР йх2

(11)

где ; ^ = ^

Е1 йх2 йх

э

Том 8, №3 (май - июнь 2016)

publishing@naukovedenie.ru

п - коэффициент формы сечения.

Решение было получено путем непосредственного интегрирования дифференциального уравнения прогибов (11) методом начальных параметров. Такой метод позволяет получить уравнения прогибов, изгибающих моментов и поперечных сил с учетом деформированного состояния конструкции.

Основным недостатком данной методики является то, что она не учитывает возможную потерю устойчивости стержня, в тех случаях, когда напряжения от изгиба значительно меньше напряжений от сжатия.

Еще одна методика расчета сжато-изгибаемых стержней, рассматривается в [4], она основана на определении деформированного состояния системы по формуле Максвелла-Мора.

В общем случае для прямолинейного стержня перемещения ьтой осевой точки находятся из следующей выражений:

горизонтальные перемещения:

где: Mui, Nui, Qui - внутренние усилия от единичной горизонтальной силы, приложенной в i-той точке;

Mvi, Nvi, Qvi - внутренние усилия от единичной вертикальной силы, приложенной в i-той точке;

Mqj, Nqj, Qqj - внутренние усилия от j-той рабочей нагрузки, приложенной произвольно,

д - безразмерная постоянная, зависящая от формы поперечного сечения.

Напряженно-деформированное состояние определяется следующим образом. После приложения нагрузки определяются внутренние усилия М, N Q и вычисляются перемещения УМ^УШ^УО! и иМ11,иШ1,и011 на 1-ом этапе деформирования

После необходимого количества итераций, находят деформированную ось конструкции от приложенной нагрузки и внутренние усилия с учетом деформированного состояния конструкции Мд, Кд, Од.

Краевые напряжения в этом случае определяются:

j=lL EiJi j—1L Ei±i j=lL Gi±i

вертикальные перемещения:

i

(i2)

Интернет-журнал «НАУКОВЕДЕНИЕ» Том 8, №3 (май - июнь 2016)

http://naukovedenie.ru publishing@naukovedenie.ru

Эта методика может быть использована при расчете большепролетных непрямолинейных конструкций, таких как арки, рамы и т.п. Так как при расчете таких конструкций отношение напряжений от изгиба к напряжениям от сжатия как правило всегда больше 0.1 [5], а также такие конструкции не подвержены потери устойчивости первого рода, что подтверждается испытаниями и теоретическими выкладками в работах Е.Н. Серова [6], А.Ю. Зобачевой [9].

Г.В. Свенцицкий [5] при изучении напряженно-деформированного состояния сжато-изгибаемых стержней, отмечал, что при малых гибкостях теория краевых напряжений дает преуменьшенную несущую способность, а при больших гибкостях получаются недостоверные результаты, так как еще до достижения краевыми напряжениями сжатия величины расчётного сопротивления стержень теряет устойчивость. Г.В. Свенцицкий [5] решил задачу устойчивости прямого стержня, рассматривая предельное равновесие между моментом внешних сил Ma = P(e + f) и моментом внутренних напряжений M, = EJ / р . Несмотря на удовлетворительные результаты, этот метод, вследствие большой трудоемкости практического применения не нашел.

Так же следует отметить, что совершенствование метода расчета сжато-изгибаемых деревянных элементов строительных конструкций описано в [10], [12]. Распределение напряжений по поперечному сечению определяется распределением деформаций в сечении и зависимостью g-8. Принимается, что деформации распределяются по линейному закону и ось нулевых деформаций совпадает с осью нулевых напряжений.

Для аппроксимации работы древесины при сжатии используется кубическая парабола

вида:

ас = А1Е — А2£3

Связь между напряжениями и деформациями в растянутой зоне сечения соответствуют закону Гука:

Ор = ЕрЕ

Главный вектор Pвн и главный момент Mвн эпюры нормальных напряжений наиболее нагруженного сечения определяют из соотношения:

rh rh h

Рвн=Ь I a(x)dx; Мвн = b I a(x)(x — —)dx, Jo Jo 2

h

где разность x — — - расстояние от рассматриваемого волокна до центра тяжести сечения.

Для нахождения условия критического состояния используется критерий, предложенный Р.С. Санжаровским [11].

SM = ЗМвн; 5Рвн = 0.

Основным недостатком этой методики является допущение о деформированном состоянии конструкции по закону синуса, что в свою очередь не позволяет применять данный метод расчета как универсальный для любого типа конструкций.

Расчет сжато-изгибаемых конструкций, так же можно проводить с помощью современных программных комплексов, таких как Lira, SCAD, основанных на МКЭ. Для учета деформированного состояния, решение ведется с учетом геометрической нелинейности. В SCAD для учета геометрической нелинейности необходимо использовать тип элемента 301 - 310. Расчет ведется в нелинейной постановке. При моделировании нелинейных нагрузок

предлагается три метода: шаговый, шаговый с учетом невязки и шагово-итерационный. Для учета геометрической нелинейности в Lira, необходимо использовать тип элемента 309. Решение ведется с учетом моделирования нелинейных нагрузок, следующими методами решения: простой шаговый, автоматический выбор шага с поиском новых форм равновесия, автоматический выбор шага для геометрически и физически нелинейных задач. В результате такого расчета получают внутренние усилия с учетом деформированного состояния конструкции Мд, Мд, Од. Краевые напряжения находятся по (12).

Применение такого метода, позволяет решать задачи любого очертания геометрической оси конструкции, но при решении задач с переменной жесткостью по длине конструкции, необходимо увеличивать число конечных элементов, либо производить решение с помощью пластинчатых или объемных элементов, что в значительной степени увеличивает временные затраты на решение задачи.

В EuroCode5 расчет сжато-изгибаемых конструкций в плоскости изгиба выполняют по следующей формуле:

°c.0.d °m.y.d < ^

кс.у • ^га.О^

где: кс.у. - коэффициент продольного изгиба [8], ^.ол, - расчетная несущая

способность сжатию и изгибу, ос.ол, Ош.у.а - максимальные напряжения от сжатия и изгиба.

Для сравнительной оценки вышеописанных методик, проведем количественное сравнение максимальных краевых напряжений, а также коэффициентов учета деформированного состояния при различных гибкостях стержня. Расчетная длинна стержня 7м, сжимающая сила N=3.2 тс, и равномерно-распределенная нагрузка д = 130 кгс/м. (см. рис. 4).

Рисунок 4. Расчетная схема сжато-изгибаемого стержня

Для сравнительного анализа, взяты 7 деревянных балок прямоугольного сечения шириной сечения 16см, и с разной высотой сечения, кратной 33 мм, для моделирования различных гибкостей стержня. Что соответствует следующему соотношению изгибных и сжимающих напряжений ош/осж=7.5/1; 5.5/0.8; 4.2/0.7; 3.3/0.7; 2.7/0.6; 2.2/0.5; 1.9/0.5.

Результаты расчетов приведены в таблице 1, 2, а также на графиках зависимости %-Х,

о-Х.

Таблица 1

Коэффициент по разным методикам

Гибкость ^снип По СНиП ^снип LIRA ^1=M/M1 SCAD ^2=M/M2 НС. Стрелецкий ^3=M/M3 По методу перемещений ^4=M/M4 По методу начальных параметров Серов Е.Н. ^5=M/M5

122 0.656 0.805 0.826 0.846 0.804 0.792

105 0.783 0.877 0.887 0.903 0.873 0.868

92 0.855 0.917 0.924 0.935 0.915 0.911

82 0.898 0.942 0.946 0.954 0.939 0.937

Гибкость ^снип По СНиП &снип LIRA &1=M/M1 SCAD &2=M/M2 НС. Стрелецкий &3=M/M3 По методу перемещений &4=M/M4 По методу начальных параметров Серов Е.Н. &5=M/M5

73 0.926 0.957 0.961 0.967 0.955 0.954

67 0.944 0.968 0.970 0.975 0.966 0.965

61 0.957 0.975 0.977 0.981 0.973 0.973

Таблица 2

Разница коэффициента £ от нормативного по СП

Гибкость ^снип &снип & 100% &снип 100% &3 &снип &3 100% &снип 100% &снип 100%

122 18.52% 20.57% 22.46% 18.41% 17.20%

105 10.71% 11.75% 13.29% 10.31% 9.83%

92 6.80% 7.43% 8.56% 6.56% 6.15%

82 4.65% 5.08% 5.87% 4.37% 4.14%

73 3.29% 3.60% 4.24% 3.04% 2.93%

67 2.48% 2.71% 3.18% 3.18% 2.19%

61 1.88% 2.06% 2.45% 2.45% 1.69%

122 105 92 82 73 67 61 Придельная гибкость, к

Рисунок 5. График зависимости f-1

= Q.

я я ЕВ К a S S й

= я - -

I 3 g li

= I

£ S

3 s

I u

£

=

S u И я

14.0

12.0

10.0

B.O

6.0

4.0

2.0

0.0

1

; Як \\\

X w

\ \\

\\

\ \

\

\

1 ч \ N 1

х\ \

122 105 92 82 73 Придельная гибкость,}.

67

61

^^Максимальные нормальные напряжения, SCAD:Mna

^^Максимальные нормальные напряжения, СНИП: Мл а

W Максимальные нормальные напряжения, Lira, Мл а

Ж Максимальные нормальные напряжения, Стрелецкий, Мпа

-♦ Максимальные нормальные напряжения, по Максвелпу-Mopv, Мпа ■ Максимальные нормальные напряжения, Формула Ясинского, Мпа ^-Максимальные нормальные напряжения, по EuroCode, Мпа

Максимальные нормальные напряжения, Метод начальных

параметров, Мпа

Рисунок 6. График зависимости о-Х Проанализировав данные графики можно сделать следующий вывод:

1.

2.

3.

При гибкости стержней Х<100, максимальные нормальные напряжения в стержнях при всех предложенных методиках имеют приблизительно одинаковый результат, за исключением расчетов по БигоСоёе 5 и расчету с помощью формулы Ф.С. Ясинского. Это связано с тем, что в этих методиках воздействие изгибающего момента накладывается на напряжения потери устойчивости центрального продольного изгиба.

При гибкости стержней Х>100, некоторые методики дают завышенные значения напряжений (до 22%).

Большинство имеющихся на сегодняшний день методик, позволяющих с большой точностью выполнять расчеты сжато-изгибаемых конструкций, в том числе и описанная в СП 64.13330.2011, позволяют рассчитывать только прямолинейные сжато-изгибаемые конструкции. Для непрямолинейных конструкций решение по существующим методикам приводит к завышенным результатам, что приводит к увеличению материалоемкости и как следствие снижению конкурентоспособности на рынке строительных материалов.

ЛИТЕРАТУРА

1. Карлсен, Г.Г. Деревянные конструкции / Г.Г. Карлсен, В.В. Большаков, М.Е. Каган, Г.В. Свенцицкий. - М.: СИ, 1961. - 642 с.

2. Стрелецкий, Н.С. Основы металлических конструкций. Часть 1 / Н.С. Стрелецкий. - М.: СИ, 1940. - 844 с.

3. Светозарова Е.И., Серов Е.Н., Душечкин С.А., Попов В.Д. Опыт проведения испытания клееных рам пролетом 18 м. - В кн.: Повышение надежности и долговечности строительных конструкций: Межвуз. темат. сб. тр. Л.: ЛИСИ, 1972. с. 97-114.

4. Шмидт, А.Б. К расчету криволинейных сжато-изгибаемых трехшарнирных систем из клееной древесины и фанеры по деформированной схеме.-Известия вузов. Строительство, 1999 г.

5. Свенцицкий Г.В. Устойчивость внецентренно сжатых деревянных стержней // Исследование прочности и устойчивости деревянных стержней: Сб. ЦНИПС. -М.: Стройиздат 1940. - с. 14-55.

6. Серов Е.Н. Рациональное использование анизотропии прочности материалов в клееных деревянных конструкциях массового изготовления: дис. ... д-ра техн. наук / Е.Н. Серов; Л.: ЛИСИ, 1989. - 521 с.

7. Клименко В.З. Феноменологический подход к расчету сжато-изгибаемых деревянных элементов / В.З. Клименко // Строительная механика и расчет сооружений. - 2011. - №1. - С. 7-11.

8. Порто Д. Руководство для проектировщиков к Еврокоду 5: Проектирование деревянных конструкций // Москва:МГСУ. 2013. C. 308.

9. Зобачева, А.Ю. Совершенствование расчетных моделей и экспериментальное обоснование методики оценки устойчивости большепролетных арочных конструкций из клееной древесины: дис. ... канд. техн. наук: Перм. нац. исслед. политехн. ун-т. - Пермь, 2012. - 111 с.

10. А.С. Вареник Устойчивость внецентренно-сжатых деревянных элементов при кратковременном нагружении // Вестник Новгородского Государственного Университета. 2013. №75 Т.1 с. 7-12.

11. Санжаровский Р.С. Устойчивость элементов строительных конструкций при ползучести. Л.: Изд-во ЛГУ, 1984. 216 с.

12. Шеховцов А.С. К вопросу устойчивости внецентренно-сжатых деревянных стержневых элементов сетчатых куполов // Промышленное и гражданское строительство. - 2007. - №3. - с. 49-50.

Grankin Kirill Vladimirovich

Saint Petersburg state university of architecture and civil engineering, Russia, Saint Petersburg

E-mail: gran1256@mail.ru

Aleksand Borisovich »Schmidt

Saint Petersburg state university of architecture and civil engineering, Russia, Saint Petersburg

E-mail: ukftpz@lan.spbgasu.ru

The calculation methods compressed-bent constructions

from glulam timber

Abstract. The authors of the article consider the existing methods calculation of linear glued laminated timber compressed-bent bars. For each of these methods calculation of compressed-bent bars including existing standard methods of the Russian Federation (Cn 64.13330.2011) and of the European Union (EuroCode 5) the errors and framework for the use methods have been found. The authors of the article have carried out comparative analysis of results of existing analytic methods and results got from solving the task using finite element analysis with modern software such as SCAD Office and Lira Soft. Also in the article the authors make numerical experiment for all being considered methods for linear compressed-bent bar with 7 m width of span. Based on results of the numerical experiment deformed status accounting coefficient- flexibility and maximums of marginal direct stress - flexibility charts are plotted. There are conclusions in the end of the article based on comparing results got from the spent numerical experiment for different methods calculation of compressed-bent bars.

Keywords: glulam timber constructions; compressed - bent constructions; geometric nonlinearity; formula Maxwell - Mora's; displacement method; method of initial parameters

REFERENCES

1. Karlsen, G.G. Derevyannye konstruktsii / G.G. Karlsen, V.V. Bol'shakov, M.E. Kagan, G.V. Sventsitskiy. - M.: SI, 1961. - 642 s.

2. Streletskiy, N.S. Osnovy metallicheskikh konstruktsiy. Chast' 1 / N.S. Streletskiy. -M.: SI, 1940. - 844 s.

3. Svetozarova E.I., Serov E.N., Dushechkin S.A., Popov V.D. Opyt provedeniya ispytaniya kleenykh ram proletom 18 m. - V kn.: Povyshenie nadezhnosti i dolgovechnosti stroitel'nykh konstruktsiy: Mezhvuz. temat. sb. tr. L.: LISI, 1972. s. 97-114.

4. Shmidt, A.B. K raschetu krivolineynykh szhato-izgibaemykh trekhsharnirnykh sistem iz kleenoy drevesiny i fanery po deformirovannoy skheme.-Izvestiya vuzov. Stroitel'stvo, 1999 g.

5. Sventsitskiy G.V. Ustoychivost' vnetsentrenno szhatykh derevyannykh sterzhney // Issledovanie prochnosti i ustoychivosti derevyannykh sterzhney: Sb. TsNIPS. - M.: Stroyizdat 1940. - s. 14-55.

6. Serov E.N. Ratsional'noe ispol'zovanie anizotropii prochnosti materialov v kleenykh derevyannykh konstruktsiyakh massovogo izgotovleniya: dis. ... d-ra tekhn. nauk / E.N. Serov; L.: LISI, 1989. - 521 s.

7. Klimenko V.Z. Fenomenologicheskiy podkhod k raschetu szhato-izgibaemykh derevyannykh elementov / V.Z. Klimenko // Stroitel'naya mekhanika i raschet sooruzheniy. - 2011. - №1. - S. 7-11.

8. Porto D. Rukovodstvo dlya proektirovshchikov k Evrokodu 5: Proektirovanie derevyannykh konstruktsiy // Moskva:MGSU. 2013. C. 308.

9. Zobacheva, A.Yu. Sovershenstvovanie raschetnykh modeley i eksperimental'noe obosnovanie metodiki otsenki ustoychivosti bol'sheproletnykh arochnykh konstruktsiy iz kleenoy drevesiny: dis. ... kand. tekhn. nauk: Perm. nats. issled. politekhn. un-t. -Perm', 2012. - 111 s.

10. A.S. Varenik Ustoychivost' vnetsentrenno-szhatykh derevyannykh elementov pri kratkovremennom nagruzhenii // Vestnik Novgorodskogo Gosudarstvennogo Universiteta. 2013. №75 T.1 s. 7-12.

11. Sanzharovskiy R.S. Ustoychivost' elementov stroitel'nykh konstruktsiy pri polzuchesti. L.: Izd-vo LGU, 1984. 216 s.

12. Shekhovtsov A.S. K voprosu ustoychivosti vnetsentrenno-szhatykh derevyannykh sterzhnevykh elementov setchatykh kupolov // Promyshlennoe i grazhdanskoe stroitel'stvo. - 2007. - №3. - s. 49-50.