CC BY
16
К РАСЧЕТУ РАМ ИЗ НЕУПРУГИХ СОСТАВНЫХ ЭЛЕМЕНТОВ С УЧЕТОМ ТЕПЛОВОГО ВОЗДЕЙСТВИЯ
А.А. РОЧЕВ, канд. техн. наук, доцент Петрозаводский государственный университет,
185026, Республика Карелия, г. Петрозаводск, пр. Комсомольский, дом 15, кв. 120. Электронный адрес: metalll@bk. ru
В работе рассматривается деформационный расчет сжато-изогнутых рам из неупругих составных элементов, имеющих переменное поперечное сечение и переменную жесткость связей сдвига по длине с учетом влияния теплового воздействия. В основу решения положена теория упругих составных стержней А.Р. Ржаницына. Статическая неопределимость рам раскрыта методом деформаций. Получено выражения для определения эквивалентного модуля деформаций, который учитывает сжимаемость ветвей, деформации сдвига материала ветвей, составляющих стержень, развитие неупругих силовых деформаций в них и влияние температурных деформаций.
КЛЮЧЕВЫЕ СЛОВА: сжато-изогнутая рама, деформационный расчет, тепловое воздействие, эквивалентный модуль деформаций.
Исследуются неупругие рамные системы, включающие в себя составные элементы, имеющие переменное поперечное сечения и переменную жесткость связей на сдвиг по длине элементов.
В работе используются основные положения общей теории упругих составных стержней, разработанной А.Р. Ржаницыным [1]. Для материала ветвей и связей между ветвями устанавливается произвольная зависимость между деформациями и напряжениями. Применяется гипотеза о нелинейно - упругом материале. Раскрытие статической неопределимости рамной системы осуществляется методом деформаций с учетом влияния продольных сил, возникающих в элементах при деформировании рам под действием нагрузки. Предполагается линейное изменение температуры по высоте поперечного сечения элементов рамы.
Для определения перемещений плоской рамы элементы многоконтурного стержня рамы делятся по длине на участки постоянной жесткости (в общем случае неравные) длиной lj между узловыми точками j и j + 1. В узлы элементов рамы вводятся дополнительные моментные 1j и силовые 2j связи, препятствующие их перемещениям. Многоконтурный стержень рамы с введенными связями образует основную систему метода деформаций. Расчет такой рамы производится шаговым методом [2].
Для определения реакций во введенных связях j-й участок элемента рамы рассматривается как жестко защемленный по концам составной стержень длиной lj или стержень, защемленный на одном конце и шарнирно опертый на другом конце. Дифференциальные уравнения для определения сдвигающих усилий
Tв шве и прогибов yji) упругого двухветвевого (n = 2) составного элемента были получены в [1]. При работе за пределом упругости эти уравнения для j -го участка стержня рамы, в котором действует продольная сила Nj), на i -м
шаге загружения, с учетом влияния температурного воздействия будут иметь вид
t?i)=+ ,
(1)
(2)
где М^ - изгибающий момент в составном стержне, лишенном связей сдвига;
С хи (/) - суммарная эквивалентная жесткость ветвей составного элемента на на ]-ом участке элемента рамы при /-ом шаге нагружения равная
п
-<еди(/) _ ^ / т?еди(/)
j = Е (Efi® Jj), (3)
здесь EjU(i) - эквивалентный модуль деформаций для сечений v -й ветви j -го участка элемента рамы на i-м шаге нагружения с учетом влияния температурных деформаций этой ветви; Jxjv - момент инерции поперечного сечения v -й ветви постоянный по длине j-го участка элемента рамы:
Aj = N$/(E%Afl) + st0]l — N«/(E«2Aj2) + stoj2 — c,Mj /C^"(i), (4)
здесь N(1) и N(2 - продольные усилия, возникающие в ветвях 1 и 2, на i -м шаге нагружения; ej и e(0j2 - температурные деформации осевых волокон ветвей 1 и 2; Ají и Aj2 - площадь поперечного сечения ветвей 1 и 2 j -го участка элемента рамы; Cj - расстояние между центрами тяжести ветвей составного элемента; E(jo1 и E(jo2 - секущие модули деформаций осевых волокон ветвей 1 и 2 на i-м шаге нагружения;
Л°) =
(1/( E() Aj v) + c)/(5)
v=1
здесь - коэффициент жесткости продольных связей сдвига на ] -м участке
элемента рамы при /-м шаге нагружения, рассчитанный по [1], с учетом температурных деформаций частей конструкции, образующей эти связи.
Величина Е^(/) в (3) определяется по формуле, которая получена автором данной статьи аналогично формуле, опубликованной в [3] и использованной в [4]:
M(i-1) h (1 — e(i-1) + e(i-1))
pequ(i) =_1V1xjv njv V1 cojv tojv '__(6)
j (AS(i—1) v(i—1)h d'(г—1) + n At(i—1)) T '
(Aej v — h j v hj v Qyjv + nj Atj v ) Txj v
где МX —1 - наибольший изгибающий момент, действующий в сечении V -й ветви ] -го участка элемента рамы на (/ — 1) -м шаге нагружения; h]1) - высота поперечного сечения V -ой ветви ] -го участка элемента рамы; е]и е(;—1 -линейная деформация оси V -й ветви ] -го участка элемента рамы на (/ — 1) -м шаге нагружения соответственно от силового и температурного воздействия;
А е]—1 = е]1) — е]'1; е1]1 и е(— 1 - краевые линейные деформации в волокнах поперечного сечения V -й ветви ] -го участка элемента рамы, возникающие в соответствии с гипотезой плоских сечений на (/ - 1) -м шаге нагружения от силового воздействия; 1 - угол сдвига V -й ветви на ] -м участке элемента рамы от единичной поперечной силы при (/ - 1)-м шаге нагружения с учетом влияния изменения температуры; Q'У]■ 1) - первая производная от поперечной
силы, действующей в и-ой ветви ]-го участка элемента рамы при (/ - 1)-м шаге нагружения; а]- коэффициент теплового линейного расширения для материала
рамы; А'^ 1 - разность температуры краевых волокон V -й ветви составного
элемента в плоскости изгиба у-го участка стержня рамы.
Используя операционное исчисление [5] сведем решение дифференциальных уравнений (1) и (2) к решению алгебраических уравнений ( для случая отсутствия свободы сдвигов в торцах ] -го защемленного по концам участка
стержня рамы, где Т'О') = 0 и Т'^1) = 0,
Т(0 = О + М- Q$$) - N • .у • '31), (7)
У]) = Т(V) + М(М1) - Q «V«, (8)
где Т]'о),Ми Q- соответственно сдвигающая сила, изгибающий момент и поперечная сила на левой опоре ] -го участка стержня рамы (в сечении Zj = 0 ) при отсутствии на нем поперечной нагрузки;
= ек^х), = о^Щ^ х) - МЕ^1^«),
12] = е^ х)/^ - х]/( (1) }),
^ = о )sh(?¿')х)/{Е%(1)}), (9)
V« = О) гш^ШЕ^^-а?), V- = [1 -со8(а«х)]/(Е9Т(%а92(0), у() = [х - 81п(а()х)/а «]/(Е—(3)J„а¡*)),
здесь а-
(0 = ^N(1)/)) . Из (8) для левой опоры, где у) = 0 получаем уравнение
т](1 V) + т(V) + м(>;(1) - QV - № (у;(1) • t« + у]) • ]3°) = о. (10)
Система линейных алгебраических уравнений (7), (8) и (10) в матричном виде будет иметь вид
А1 = В, (11)
где А - квадратная матрица коэффициентов при неизвестных в системе уравнений; В - вектор свободных членов системы уравнений; 2 - вектор неизвестных усилий:
2 = [Т]0\ М«, Q«]. (12)
В системе уравнений (11) матрица А будет включать нижеследующие элементы, полученные из (9) при х = ]
А С1') = ' ') А(1) = t '(1) А(1) = -'')
Л11] ~'о]1 ' Л12] ~ '1 ] ; Л13] ~ '2] '
А О = 'с)V(1') л(0 = '^М1^ + 1/(1) А(1) = '(1)л,(1) -V(1') (13)
А21 ] -/о]Го}1 , А22] — Г1 ] о]1 У1 ]-1 , А23]~'о^о]1 У2]1 , (13)
л(0 = ' ' (1')V(1') + ' (í)V' С ) л(0 = ' V/1) + 1/ '(1) + ' V/1)
А31 ] -'о]1 о]1 о]1 о] , А32] -'1 ] Уо]1 + У1 ] +11 ]Го]1 ,
А(0 = '(1')V(1') - '()V '(1) -V(1') л33 ]- (1 Ауо]1 '2 2 А-
Вектор В для случая смещения правой опоры на 8 = 1 будет включать
следующие элементы:
в (0 = ( у ' (0' (г) у (г)' ' (г)) В И = ^(Г)л/Г) + 1 + у ^ )V(1') (14)
В1) = (у '3-у]1'о}1 ), в2] = V1 ] +1 + ууг '3, (14)
в (о = (у с(1')V(1') + у ' О' а )V 'О + у О' а )V к1)) N(1) + )V ,(1)
в3] = (у]1 13]УоА + у;г '3з^о] + у;?'3)+ К) ^А ,
где К- = а- • А'(1) /к-; А'((1)- разность температур краевых волокон составного элемента; к - высота поперечного сечения составного элемента.
Аналогично с использованием (11) можно найти опорные реакции для других случаев: поворот правой опоры на угол ф = 1, шарнирное опирание j-го стержня на правую опору, отсутствия препятствий для сдвига в шве на опорах составного элемента.
Перемещения узлов стержня рассматриваемой рамы Zjj1 и Z 2j определяются из совместного решения уравнений:
m
2 (RgirZ« + R(j2rZ«) + Rj + Rj = 0, (15)
r=1
m
2 (RM? + Rdj 2rZ 2r) + R(j + R(it = 0, (16)
r=1
где Rj и Rt({j2r - реакции в связях 1 j основной системы метода деформаций от единичного перемещения связей 1r и 2r, определенные с учетом влияние продольных сил на i -м шаге нагружения; Rt(2j1r и Rt(2j2r - тоже, но в связях
2 j; Rj и R2 j - реакции в связях 1j и 2j от силового воздействия на i -м
шаге нагружения; Rj и Ryt - реакции в связях 1j и 2j от теплового воздействия на i-м шаге нагружения; m - число узлов с введенными дополнительными связями. Найденные усилия на опорах j-го стержня при S = 1 и ф = 1 далее используются как реакции в связях R( )1r, R12r, R^^,. и R^)2r в уравнениях (15) и (16).
Деформационный расчет рамы с учетом теплового воздействия осуществляется с использованием эквивалентного модуля деформаций EjU(i), величина которого рассчитывается на каждом следующем шаге расчета по результатам, полученным на предыдущем шаге расчета. Значения параметров, характеризующих НДС рамной системы, найденные на i-м шаге нагружения, могут быть в дальнейшем использованы для проверки прочности и устойчивости сжато-изогнутых рам методом проф. Р.С. Санжаровского [6], а также для решения задач, аналогичных рассмотренным в работе[7].
Л и т е р а т у р а
1. Ржаницын А.Р. Составные стержни и пластинки. - М.: Стройиздат, 1986. 314 с.
2. Биргер И.А. Общие алгоритмы решения задач теории упругости, пластичности и ползучести // Успехи механики деформируемых сред. - М.: Наука , 1975. - С.61 - 73.
3. Рочев А.А. Нелинейная теория расчета сквозных упругопластических статически неопределимых рамных систем // Доклады 58-й конференции профессоров, преподавателей, научных работников, инженеров и аспирантов университета. В 3 ч. Ч. 1. -СПб.: СПбГАСУ, 2001. - С.93 - 94.
4. Рочев А.А. Алгоритм нелинейного расчета круговой составной арки // Ученые записки Петрозаводского государственного университета. Сер. « Естественные и технические науки». - 2010. - №2 (107). - С.25-29.
5. Араманович И.Г. Функции комплексного переменного. Операционное исчисление. Теория устойчивости /И.Г. Араманович, Г.Л. Лунц, Л.Э. Эльсгольц. - М.: Наука, 1968. - 416 с.
6. Санжаровский Р.С. Устойчивость элементов строительных конструкций при ползучести. - Л.: Изд-во ЛГУ, 1984. - 280 с.
7. Chi Kin Iu, Siu Lai Chan A simulation-based large deflection and inelastic analysis of steel frames under fire // Journal of Constructional Steel Research, Volume 60, Issue 10, October 2004. - P. 1495-1524.
References
1. Rzhanitzin, AR (1986). Composite rods and plates. M.: Stroyizdat, 314 p.
2. Birger, LA (1975). Common algorithms for solving problems in the theory of elasticity, plasticity and creep. UspehiMechanikiDeformiraemyh Sred, M.: Nauka, p. 61-73.
3. Rochev, AA (2001). Nonlinear theory of analysis of elastic-and-plastic statically indeterminate frame systems. Proceedings of the 58th Conference of professors, teachers, researchers, engineers and graduate students. Part 1, St. Petersburg: SPGASU, p. 93 - 94.
4. Rochev, AA (2010). Algorithm for non-linear analysis of circular composite arch. Proceedings of Petrozavodsk State University. Ser. "Natural and engineering sciences.", Number 2 (107), p. 25-29.
5. Aramanovich, IG, Luntz, GL, El'sgol'ts LE (1968) Functions of Complex Variable. Operational Calculus. Theory of Stability. M .: Science, 1968. 416 p.
6. Sanzharovskiy, RS (1984). Stability of StructuralELements under Creep. L.: Izd-vo LGU, 280p.
7. Chi Kin Iu, Siu Lai Chan A (2004). Simulation-based large deflection and inelastic analysis of steel frames under fire. Journal of Constructional Steel Research, Vol. 60, Iss. 10, p. 1495-1524.
ON ANALYSIS OF FRAMES MADE OF NON-ELASTrc COMPOSITE ELEMENTS WITH TERMAL ACTION TAKEN INTO ACCOUNT
In this paper, the author solves the problem of calculating the deformation compressed-bent frames of inelastic composite elements having a variable cross-section and variable stiffness relations shift lengths for the effects of heat. The solution is based on the A.R. Rzhanit-sin's theory of elastic composite bars. Static indeterminacy frames revealed by deformation. The expressions for the equivalent modulus of deformation, which takes into account the compressibility of the branches, the branches of the shear deformation of the material, form the backbone of the development of inelastic deformation in them, and the effect of temperature deformations.
KEY WORDS: compressed-bent frame deformation calculation, thermal effects, equivalent modulus of deformation.
Rochev A.A.
Petrozavodskiy Gosudarstvenniy Universitet, Petrozavodsk
