Научная статья на тему 'Алгоритм смешанного метода для определения прогибов в деревянных балках при учёте длительного загружения'

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

CC BY
117
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
FINITE ELEMENT OF THE MIXED METHOD / VARIABLE MODULUS OF ELASTICITY / DEFLECTION OF WOODEN BEAMS / LONG-TERM LOADING / КОНЕЧНЫЙ ЭЛЕМЕНТ СМЕШАННОГО МЕТОДА / ПЕРЕМЕННЫЙ МОДУЛЬ УПРУГОСТИ / ПРОГИБЫ ДЕРЕВЯННЫХ БАЛОК / ДЛИТЕЛЬНОЕ НАГРУЖЕНИЕ

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Раков В. Д., Тухфатуллин Б. А., Путеева Л. Е.

В статье изложен алгоритм определения прогибов в деревянных балках прямоугольного поперечного сечения с учётом длительного загружения. Алгоритм расчёта основан на предложенном Д.К. Арлениновым подходе задания модуля упругости древесины в зависимости от уровня расчётных напряжений в табличной форме. Для решения поставленной задачи используется конечный элемент смешанного метода с переменным по длине элемента и высоте поперечного сечения модулем упругости. Приведены примеры расчёта шарнирно опертых балок при различных условиях загружения. Полученные численные результаты подтверждаются экспериментальными данными, приведёнными в работах других авторов.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по строительству и архитектуре , автор научной работы — Раков В. Д., Тухфатуллин Б. А., Путеева Л. Е.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Algorithm of the mixed method for determining the deflections in wooden beams, taking into account the long loading

TThe article presents an algorithm for determining deflections in wooden beams of rectangular cross-section, taking into account the long load. The calculation algorithm is based on the proposed approach of D. Arleninov to specify the modulus of elasticity of wood, depending on the level of design stresses in tabular form. The aim of the research is to develop a software oriented algorithm for calculating deflections of wooden beams under long-term loading. The object of study a wooden beam of rectangular cross-section, located in a single load; the modulus of elasticity of the material varies both in length and height of the beam, depending on the level of design stresses. The algorithm uses the FEM of a mixed method with a variable length modulus of elasticity. The stiffness characteristics of the cross sections are refined during the calculation by a two-stage algorithm. Examples of the calculation of simply supported beams under various conditions of loading. The results of comparison with the experimental data of other authors indicate the correct operation of the proposed algorithm; the discrepancy between the numerical results and the experimental results is in the range from 0,31% to 11,6 %. This algorithm can also be used to clarify deflections and internal forces in flat statically indeterminate rod systems, whose modulus of elasticity varies depending on the acting stresses

Текст научной работы на тему «Алгоритм смешанного метода для определения прогибов в деревянных балках при учёте длительного загружения»

Алгоритм смешанного метода для определения прогибов в деревянных балках при учёте длительного загружения

В.Д. Раков, Б.А. Тухфатуллин, Л.Е. Путеева Томский государственный архитектурно-строительный университет

Аннотация: В статье изложен алгоритм определения прогибов в деревянных балках прямоугольного поперечного сечения с учётом длительного загружения. Алгоритм расчёта основан на предложенном Д.К. Арлениновым подходе задания модуля упругости древесины в зависимости от уровня расчётных напряжений в табличной форме. Для решения поставленной задачи используется конечный элемент смешанного метода с переменным по длине элемента и высоте поперечного сечения модулем упругости. Приведены примеры расчёта шарнирно опертых балок при различных условиях загружения. Полученные численные результаты подтверждаются экспериментальными данными, приведёнными в работах других авторов.

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

Введение

В настоящее время в мировой практике в качестве строительного материала широко используется древесина как экологичный и возобновляемый природный материал [1]. Для конструкций, выполненных из дерева, наряду с требованиями прочности и устойчивости, обеспечение условия жёсткости также является немаловажным.

Исследования [2, 3] посвящены разработке инженерной методики определения прогибов деревянных балок при длительном загружении с учётом ползучести материала, опирающейся на данные экспериментов ряда авторов [4 - 7]. Предложенная методика заключается в определении прогибов деревянных балок с учётом того обстоятельства, что модуль упругости древесины является переменной величиной, зависящей от уровня расчётных напряжений, возникающих в конструкции при длительном нагружении. Закон изменения модуля упругости в зависимости от уровня напряжений по данным профессора Д.К. Арленинова [3] представлен в табл. 1.

Таблица № 1

Значение модуля упругости в зависимости от уровня нормальных напряжений

а, МПа До 5 6 7 8 9 10 11 12 Свыше 13

E, МПа 10000 9200 8300 7400 6500 5600 5000 4400 3900

В предложенной методике [3] конструкция разбивается на ряд участков с постоянным модулем упругости, принимаемым в соответствии с данными табл. 1, по максимальному в пределах участка напряжению. Таким образом, результаты расчёта будут напрямую зависеть от количества назначенных участков постоянной жёсткости. Следует отметить, что предложенная инженерная методика не учитывает влияния изменения уровня напряжений по высоте сечения на величину модуля упругости; используются данные только о фибровых напряжениях. На наш взгляд это объясняется нацеленностью методики [3] на применение программных комплексов SCAD, ЛИРА САПР и др., в которых возможности задания переменного по высоте сечения модуля упругости в стержневой системе ограничены.

В предыдущих исследованиях [8, 9] был предложен конечный элемент (КЭ) смешанного метода с переменным по длине модулем упругости, закон изменения которого задан в явном виде или таблично. Для вычисления элементов матрицы откликов [10], связывающей между собой реакции в отброшенных связях с перемещениями в добавленных связях, используется процедура численного интегрирования. Модификация алгоритма, предложенная в [11], позволяет вычислять жесткостные характеристики сечения в соответствии с эпюрой нормальных напряжений.

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

N

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

Материалы и методы

Для вычисления модуля упругости с учётом данных табл. 1 необходимо вычислять величины нормальных напряжений как по длине КЭ, так и по высоте поперечного сечения. В связи с этим предлагается следующий двухэтапный алгоритм. Каждый КЭ разбивается по длине на пх участков (рис. 1, а), а поперечное сечение в пределах участка разбивается по

высоте на п слоёв (рис. 1, б).

а

б

У

и

1

1 / %ау

х

пУ

Рис. 1. - Схемы разбиения КЭ: а) по длине; б) по высоте сечения

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

( \ М (х) (1)

^ У ) . (1)

Ъ

N

В формуле (1) усилие М(х) принимается без учёта знака. По найденным напряжениям определяется модуль упругости Е(у) для очередного слоя в соответствии с данными табл. 1. Полученные значения Е(у) подставляются в формулы для вычисления продольной Е(х)А и изгибной Е(х)Jz жёсткостей в сечении с координатой х:

Е (х) А = ]ГЕ (у)Ъёу, (2)

i=1

E ( х) Jz = £ E (y)

ЬШ3 ^ 2

+ bdy ■ y; 12 г

(3)

где b - ширина поперечного сечения; h - высота поперечного сечения; h

dy = — - высота слоя (рис. 1, б). Результаты вычислений жесткостей E (х) A

ny

и E(x)Jz по формулам (2), (3) используются для формирования матрицы откликов и вектора грузовых коэффициентов и расчёта системы на втором этапе. Осевую жёсткость E(х)A требуется вычислять из-за особенностей алгоритма смешанного метода [8, 9].

Результаты и обсуждение Пример № 1. Рассматривается шарнирно опертая балка [4] пролётом i = 6 м прямоугольного поперечного сечения с шириной b = 12 см и высотой h = 33 см, загруженная равномерно распределённой нагрузкой q = 4,8 кН/м (рис. 2, а). Геометрические характеристики сечения (момент инерции Jz и момент сопротивления WZ) определены по формулам:

= 12-33^ = 35937 см4; w = 12^ = 2178см3. z 12 6

Изгибающий момент Mz max и нормальные напряжения a max составят:

л о /г2 OIA 1 Г\3

Mz max = -3—— = 21,6 кН ■ м; amax =--- = 9,9 ■ 106 Па = 9,9 МПа.

z ,max 8 max 2178 10~ 6

i=1

а

б

в

.....д.....

н

> <

£ £

3 3

2

2

Рис. 2. - Расчётные схемы: а) пример № 1; б) примеры № 2, 3; в) пример № 4

По результатам испытания двух балок в [4] установлены величины упругого /0(э) и длительного /(э) прогибов (время испытания г = 367 суток):

/0(э) = 23,0 мм; /0(э2) = 22,5 мм; /Д3-1 = 30,2 мм; = 28,8 мм. Результаты вычисления

прогибов по предложенному алгоритму смешанного метода составили: /0(мкэ) = 22,3 мм; /(мкэ) = 26,67 мм.

Пример № 2. Однопролётная шарнирно опертая балка [3, 12] прямоугольного поперечного сечения загружена двумя сосредоточенными силами ^, приложенными в третях пролёта (рис. 2, б). Исходные данные для расчёта: £ = 2,9 м; Ь = 10 см; И = 20,3 см; ^ = 8,82 кН. Вычисленные величины: Jz = 6971 см4; = 687 см3; Ыг,тах = 8,53 кН• м; атах = 12,4 МПа. Эксперимент

проводился в течение 8 месяцев; прогибы в середине пролёта составили: /0(э) = 11,3 мм; /(э) = 17,5 мм. По предложенному алгоритму вычислены прогибы:

/0(мкэ) = 10,8 мм; /(мкэ) = 15,9 мм.

Пример № 3. Для той же схемы балки, что и в примере № 2 приняты следующие исходные данные для расчёта: £ = 1,46 м; Ь = 10,0 см; И = 2,7 см; ^ = 0,317 кН [6]. Результаты расчёта: Jz = 16,4 см4; = 12,15 см3;

Ыг тах = 0,158 кН • м; атах = 13 МПа. Время проведения эксперимента - 6 месяцев; прогибы в середине пролёта составили: /0(э) = 15,6 мм; /(э) = 32,0 мм. Прогибы, вычисленные по предложенному алгоритму:

/0(мкэ) = 21,2 мм; ¡((мкэ) = 31,9 мм.

Пример № 4. Однопролётная шарнирно опертая балка [5] прямоугольного поперечного сечения Ь = 9,4 см и высотой к = 1,91 см загружена в середине пролёта I = 1,5 м сосредоточенной силой Г = 0,196 кН (рис. 2, в). Кратковременный прогиб составил: /0(э) = 26,1 мм, а после 8 месяцев испытаний увеличился до / {(э) = 33,4 мм. Результаты, полученные с использованием смешанного метода: /0(мкэ) = 24,9 мм; (мкэ) = 32,4 мм.

На рис. 3 представлены графики изменения вычисленной величины прогиба в зависимости от принятого по высоте сечения количества слоёв. Для примеров № 1 - № 4 количество слоёв принималось равным 2, 4, 6, 8, 10 и 12. Как следует из приведенных графиков, достаточным будет назначать число слоёв от 10 до 12.

Рис. 3. - Результаты вычисления прогибов в балках в зависимости от назначенного по высоте сечения числа слоёв: а) пример № 1; б) пример № 2, в) пример № 3; г) пример № 4

Сводка результатов и сравнение с инженерной методикой [3, 12] приведена в табл. 2.

Таблица № 2

Результаты вычисления прогибов в середине пролётов балок

Номер примера Данные эксперимента Инженерная методика [3, 12] Смешанный метод

/о(э), мм /э), мм / т ), мм Дмкэ) f 0 % мм /(мкэ), мм

1 23 30,2 31,5 (-4,3 %) 22,3 (3,04 %) 26,7 (11,6 %)

22,5 28,8 31,5 (-9,4 %) 22,3 (0,89 %) 26,7 (7,3 %)

2 11,3 17,5 19,2 (-9,7 %) 10,8 (4,4 %) 15,9 (9,1 %)

3 15,6 32,0 — 21,2 (-35,9 %) 31,9 (0,31 %)

4 26,1 33,4 - 24,9 (4,6 %) 32,4 (3,3 %)

Заключение

В статье предложен алгоритм вычисления прогибов в деревянных балках постоянного по длине поперечного сечения, работающих в условиях длительного загружения. Алгоритм использует КЭ смешанного метода с переменным по длине модулем упругости. Жесткостные характеристики поперечных сечений уточняются в ходе расчёта по двухэтапному алгоритму. Результаты сравнения с экспериментальными данными других авторов свидетельствуют о корректной работе предложенного алгоритма; расхождение численных результатов с результатами эксперимента находится в диапазоне от 0,31 % до 11,6 %. Данный алгоритм также может быть использован для уточнения прогибов и внутренних усилий в плоских статически неопределимых стержневых системах, модуль упругости которых меняется в зависимости от действующих напряжений.

Работа выполнена при поддержке государственного задания Министерства образования и науки РФ (номер проекта 7.8899.2017/8.9).

Литература

1. Christoforo A.L. Numerical evaluation of the modulus of longitudinal elasticity in structural round timber elements of the Eucalyptus genus / Christoforo A.L. [and others]. // Eng. Agric. Vol. 31. № 5. Jaboticabal Sept./Oct., 2011. pp. 110.

2. Арленинов Д.К., Арленинов П.Д. Переменный модуль упругости древесины // Вестник МГСУ, 2011. № 1. С. 150-152.

3. Арленинов Д.К. Инженерный метод учёта ползучести древесины при определении прогибов // Вестник МГСУ, 2011. № 1. С. 147-149.

4. Рощина С.И. Прочность и деформативность клеёных армированных деревянных конструкций при длительном действии нагрузки: автореферат дис. ...докт. техн. наук. Москва, 2009. 38 с.

5. Арленинов Д.К., Линьков Н.В. Предложения по методике испытаний для оценки ползучести древесины // Известия КГАСУ, 2017. № 4 (42). С. 136141.

6. Арленинов Д.К., Беккер Д.А. Влияние уровня напряжений на ползучесть древесины при изгибе // Известия вузов. Лесной журнал, 2015. № 6. С. 128-137.

7. Фурсов В.В., Пурязданхах М., Бидаков А.Н. Сравнительный анализ результатов теоретических и экспериментальных исследований натурной арки из клеёной древесины // Инженерный вестник Дона, 2014. № 2. URL: ivdon.ru/ru/magazine/archive/n2y2014/2395

8. Тухфатуллин Б.А., Путеева Л.Е., Раков В.Д. Разработка и апробация конечного элемента смешанного метода для расчёта стержневых систем с элементами переменной жёсткости // Инвестиции, строительство, недвижимость как материальный базис модернизации и инновационного развития экономики: материалы VIII Международной научно-практической конференции, 13-15 марта 2018 г.: в 2 ч. Ч. 1. Томск: Изд-во Том. гос. архит.-строит. ун-та, 2018. C. 456-466.

9. Тухфатуллин Б.А., Путеева Л.Е., Красина Ф.А. Расчёт плоских стержневых систем с переменным по длине модулем упругости методом конечных элементов в форме смешанного метода // Инженерный вестник Дона, 2018. № 3. URL: ivdon.ru/ru/magazine/archive/n3y2018/5111.

10. Ignatyev A.V., Ignatyev V.A., Onischenko E.V. Analysis of Systems with Unilateral Constraints through the Finite Element Method in the Form of a Classical Mixed Method // Procedia Engineering. Vol. 150. 2016. pp. 1754— 1759.

11. Тухфатуллин Б.А., Путеева Л.Е., Раков В. Д., Подшивалов И.И. Расчёт деревянной рамы с использованием конечного элемента смешанного метода // Информационно-вычислительные технологии и их приложения: сборник статей XXII Международной научно-технической конференции, 27-28 августа 2018 г. / МНИЦ ПГАУ. Пенза: РИО ПГАУ, 2018. C. 109-113.

12. Арленинов Д.К. Расчётная оценка прогибов деревянных балок при длительной нагрузке // ПГС, 2012. № 11. С. 40-41.

References

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1. Christoforo A.L. Eng. Agric. Vol. 31. № 5. Jaboticabal Sept. Oct. 2011. pp. 110.

2. Arleninov D.K., Arleninov P.D. Vestnik MGSU, 2011. № 1. pp. 150-152.

3. Arleninov D.K. Vestnik MGSU, 2011. № 1. pp. 147-149.

4. Roshhina S.I. Prochnost' i deformativnost' klejonyh armirovannyh derevjannyh konstrukcij pri dlitel'nom dejstvii nagruzki [Durability and deformability of glued reinforced wooden structures under prolonged load action]: avtoreferat dis. ...dokt. tehn. nauk. Moskva, 2009. 38 p.

5. Arleninov D.K., Lin'kov N.V. Izvestija KGASU, 2017. № 4 (42). pp. 136-141.

6. Arleninov D.K., Bekker D.A. Izvestija vuzov. Lesnoj zhurnal, 2014. №6, pp. 17.

7. Fursov V.V., Purjazdanhah M., Bidakov A.N. Inzhenernyj vestnik Dona (Rus), 2014. № 2. URL: ivdon.ru/ru/magazine/archive/n2y2014/2395

8. Tukhfatullin B.A., Puteeva L.E., Rakov V.D. Investicii, stroitel'stvo, nedvizhimost' kak material'nyj bazis modernizacii i innovacionnogo razvitija jekonomiki: materialy VIII Mezhdunarodnoj nauchno-prakticheskoj konferencii, 13-15 marta 2018 g.: v 2 ch. Ch. 1. Tomsk: Izd-vo Tom. gos. arhit.-stroit. un-ta, 2018. pp. 456-466.

9. Tukhfatullin B.A., Puteeva L.E., Krasina F.A. Inzhenernyj vestnik Dona (Rus), 2018, № 3. URL: ivdon.ru/ru/magazine/archive/n3y2018/5111.

10. Ignatyev A.V., Ignatyev V.A., Onischenko E.V. Procedia Engineering. Vol. 150. 2016. pp. 1754-1759.

11. Tukhfatullin B.A., Puteeva L.E., Rakov V.D., Podshivalov I.I. Informacionno-vychislitel'nye tehnologii i ih prilozhenija: sbornik statej XXII Mezhdunarodnoj nauchno-tehnicheskoj konferencii, 27-28 avgusta 2018 g. MNIC PGAU. Penza: RIO PGAU, 2018. pp. 109-113.

12. Arleninov D.K. PGS, 2012. № 11. pp. 40-41.

i Надоели баннеры? Вы всегда можете отключить рекламу.