Математическая модель рамной конструкции из составных упругопластических элементов
Автор: Рочев Анатолий Алексеевич
Журнал: Ученые записки Петрозаводского государственного университета @uchzap-petrsu
Рубрика: Архитектура и строительство
Статья в выпуске: 1 (90), 2008 года.
Бесплатный доступ
Построена математическая модель, позволяющая рассчитать напряженно-деформированное состояние и проверить устойчивость рамной конструкции из составных упругопластических элементов. Переменные состояния модели учитывают изменение геометрических параметров элементов системы по длине, податливость их на изгиб, сдвиг и кручение, нелинейные геометрические и физические факторы, возникающие при работе элементов под нагрузкой.
Рамная конструкция, деформированное состояние, составной упругопластический элемент
Короткий адрес: https://sciup.org/14749366
IDR: 14749366
Текст научной статьи Математическая модель рамной конструкции из составных упругопластических элементов
Исследуется рамная конструкция, включающая в себя составные прямолинейные элементы с ветвями, выполненными из тонкостенных профилей, соединенных между собой решеткой с образованием элемента переменного сечения по длине, имеющего две оси симметрии в поперечном сечении. Элементы рамы имеют начальные несовершенства в виде пространственно изогнутых осей. Исследование устойчивости осуществляется на примере рамной конструкции пролетом l , имеющей произвольный неразветв-ленный ломаный контур и загруженной по длине произвольной поперечной нагрузкой G , приложенной с эксцентриситетами относительно продольной оси. Из своей плоскости рамная система имеет жесткие опоры в узлах и между ними (в местах расположения связей). Узлы системы элементов - жесткие. Сечения элементов в узлах системы не имеют возможности поворота в своей плоскости.
В работе используется приближенное выражение для кривизны оси элементов. Для материа- ла элементов устанавливается произвольная зависимость между деформациями и напряжениями. Влияние разгрузки не учитывается. Используется гипотеза неплоских сечений [1]. Учет деформаций сдвига осуществляется способом, предложенным Ф. Энгессером и С. П. Тмошенко [2]. Не учитывается влияние касательных напряжений на развитие пластических деформаций. Геометрическая неизменяемость поперечного сечения элементов обеспечивается постановкой поперечных диафрагм жесткости. По торцам составных элементов рамной системы предполагается наличие опорных плит, препятствующих депланации поперечного сечения.
Первоначально, аналогично [3], исследование рамной конструкции осуществляется без учета влияния деформаций сдвига. Используется шаговое загружение конструкции. В деформированном состоянии хорда i -го элемента рамы, расположенного между узлами i - 1 и i , повернется на угол ηi относительно первоначального положения. В этом состоянии уравнения равно- весия части силовой системы, отделенной сечением, проходящим через центр узла i перпендикулярно к хорде i -го элемента, будут иметь вид:
R oy o — G y o — P i ( i + 1) cos(t> i + / i) + Q i(i + 1) sin( U i + / i ) = 0 ,
R o z o — G z o — P ( i + 1) sin( U i + / ) + Q i(i + 1) cos( U i + / ) = 0,(1)
RoyoZoi — RoZoyoi — MG + Mi(i+1) = 0 , где Royo и Rozo - проекции главного вектора реакций опоры с координатами Xo = yo = zo = o на координатные оси yo и z0; (Р. и Giyo - проекции главного вектора внешних сил на оси у0 и zo; MG - проекция главного момента внешних сил относительно центра узла i ; P(i+1), Qi(i+1) -силы, действующие соответственно вдоль и поперек хорды i -го элемента со стороны элемента i +1 ; Mi(i+i) - момент, действующий в i -м узле со стороны элемента i +1; и - угол между осью уо и осью i -го элемента недеформированной силовой системы; уoi и zoi - координаты i -го узла деформированной системы.
Координаты узлов деформированной системы элементов Z oi и y oi определяются из выражений, учитывающих измененную вследствие деформаций длину l i элементов системы:
)
k = i — 1
У oi = £ У o , k + l i cos( u i + / i ) , k = 1 (2)
Z oi = zo , k + 6 sin( U i + / ).
Связь между силами и моментами, действующими по разные стороны i -го узла, определится из решения системы уравнений равновесия
P ( i + 1) cos( U i + 1 + / i + 1 ) — Q i( i + 1) sin( U i + 1 + / i + 1 ) —
— P i ( i + 1) cos( U i + / i ) + Q i — 1) sin( U i + / i ) = o,
P i ( i + 1) cos( U i + 1 + / i + 1 ) — Q i( i + 1) cos( U i + 1 + / i + 1 ) — (3)
— P i ( i + 1) sin( и + / i ) — Q i( i — 1) cos( U i + / i ) = o,
M i ( i + 1) = M o— ) .
Разделив длину i -го элемента системы на m равных частей длиной si, изогнутую ось i -го элемента в местной системе координат можно описать интерполяционными многочленами Лагранжа m -й степени xoi = Pm1(zoi) , yoi = Pm2(zoi) . (4)
В j -м сечении i -го элемента системы прогиб будет складываться из
——. >——
——— >——
X oij = X oij + X oij , y oij = y oij + y oij ,
где —oij и —oij - начальные погибы в j -ом сечении i -го элемента системы; Xoj и yoj - прогибы в j -м сечении i -го элемента после приложения к силовой системе внешних сил.
Проведя в соответствии с [1] через точку хорды i -го элемента деформированной системы с координатой Z oi = Z oij цилиндрическое поперечное сечение радиусом r j получим в волокне радиального направления деформацию
E rij = E r o ij + E r1 ij + E r 2 ij + ^ r o ij ,
где E r oij , E r и, и E r 2ij - деформации, возникающие от действия соответственно продольной силы N ij , изгибающих моментов M xj и M yj ; E r m ij - деформация, возникающая при стесненном кручении [4].
С учетом (5) относительные деформации в радиальном направлении определятся по формулам:
E roij = E oij cos в ,
E r 1ij = — yoij r ij sin 2 в / 2 , E r 2 ij = — x o ij x , (7)
Eraij = — ^ij^rij , где Arij - площадь поперечного j -го цилиндрического сечения i -го элемнта ; в - угол между направлением рассматриваемого волокна и плоскостью xoz; torij - секториальная координата точки контура j -го цилиндрического сечения i -го элемента.
Угол ^ ij определяется как сумма приращений углов поворота сечений i -го элемента системы начиная от сечения с координатой z oi = o.
k = j g = k
V ij = ^ ( A M itZ k k ^ ^ "1 tg ), (8)
k = 1 g = 1
где N M itz k k - приращение крутящего момента по длине k -го участка i -го элемента системы относительно оси, проходящей через центр изгиба сечения; V 1 tg - угол закручивания g -й пространственной панели k -го участка элемента от действия единичных крутящих моментов, приложенных по ее торцам, определяемый аналогично [6] с учетом развития пластических деформаций [7]; k - число панелей на k -ом участке.
Выражение для определения V аппроксимируем интерполяционным многочленом k,=+2
V ij = ^ D j + k ^ ij + k (9)
k , =— 2
Условия равновесия для части i -го элемента, отделенной j -м сечением элемента, в проекциях на центральные оси инерции этого сечения x 1 , y 1 и Z 1 параллельные координатным осям x oi , y oi и Z oi с учетом деформированного состояния элемента запишутся в виде:
P i ( i + 1) y oij - M i ( i + 1) + Q i( i + 1) Z oj - M G + M y cos l ij 2 cos V ij - M y cos ^ ij 3 sin V ij - M y’ cos 6 j 1 sin x Oij = 0, P i ( i + 1) X oij - M Ryi + M y cos < ij 2 sin V ij - M y’ cos l ij 3 cos V ij - M iy ’ sin l j 1 = 0,
Q i( i + 1) X oj ± M G - M y ’ sin l 2 + M y ’ sin l 3 - M y cos l j1 cos x Oij = 0, ± Ry, + Q y ’ cos l ij 2 cos V ij ± Q y ’ cos < ij 3 sin V ij ± P, y’ cos l j1 sin X oij = 0, Q ( i + 1) - G ijy - + Q y" ’ cos l ij 2 sin V ij ± Q y ’ cos l ij 3 cos V ij ± P, y’ sin l ij 1 = 0, P i ( i + 1) + Q i/X" ’ sin l ij 2 ± Q y ’ sin l ij 3 ± P y’ cos l ij 1 cos X 'oij = 0,
гдеξij1=yo′ij[1-(xo′ij)2/2],ξij2=xo′ij[1-(ψij)2/2], ξij3 = yo′ij[1-(ψij)2/2], Gijy1 - проекция главного вектора внешних сил, действующих на i -й элемент системы; MiGjx и M iGjz - проекции вектора главного момента внешних сил относительно центра тяжести j -го сечения i -го элемента системы на оси соответственно x1 и z1 ; Pij‰’ -главный вектор эпюры нормальных напряжений j -го сечения i -го элемента системы; Mi‰jx ’ и Mi‰jy ’ - проекции вектора главного момента эпюры нормальных напряжений j -го сечения i -го элемента системы относительно центра тяже- сти этого сечения на главные центральные оси инерции; M i‰jz ’ - главный вектор касательных сил j -го сечения i -го стержня; Qi‰jx ’ и Qi‰jy ’ -проекции главного вектора эпюры касательных напряжений изгиба на главные центральные оси инерции; Rijx1 - главный вектор горизонтальных опорных реакций; M iRjy - проекция главного момента от действия горизонтальных опорных реакций на ось y1 j -го сечения i -го элемента системы.
Внутренние силы j-го сечения в (10) опреде- лятся из
P ij ‰’ = r ij ∫ σ rij cos β dA rij , A rij
M i‰jx’ = r ij ∫ σ rij y cos β dA rij , A rij
M i‰jy’ = ∫ σ rij x cos β dA rij . A rij
Уравнения, раскрывающие статическую неопределимость задачи, будут иметь вид:
y o ′ 1 o = 0, y o ′ km = 0, l = const . (12)
При составлении выражения для l в (12) учитывается то, что длина стержней деформированной системы с учетом деформаций изгиба, осевого растяжения-сжатия равна li li = li + Δli , Δli li
ε oi dx oi . (13)
Используя (1) - (13), решается задача определения напряженно-деформированного состояния рамной конструкции без учета деформаций сдвига.
Учет деформаций сдвига дает следующие углы поворота сечений yo′ ij = y′oij + Qijвyнγˆ1yj , xo′ij = xo′ ij + Qiвjxнγˆ1xj ,
Qij‰y ’ =Qi‰jy’± ∫σrijsinβdArij , (14) Arij где γˆ1xj и γˆ1yj - усредненные углы поперечного сдвига j -го сечения в главных пло скостях инерции на j -м участке системы от действия единичных поперечных сил, определяемые из рассмотрения картины деформации пространственной панели составного элемента с учетом влияния пластических деформаций [7].
Зная y o ′ ij , x o ′ ij и используя формулы чи с ленного анализа [9], можно найти величины y oij и x oij . После подстановки их в (10) вместо величин y oij , y ′ oij , x oij , x o ′ ij и учета влияния деформаций сдвига в выражениях (2) и (13) задача устойчивости рамной конструкции решается методом, изложенным в [8].
Список литературы Математическая модель рамной конструкции из составных упругопластических элементов
- Верховский А. В. Новый способ определения напряжений в деталях сложной формы (метод неплоских сечений)//Тр. Горьковского политехн. института. Вып. 1. Горький, 1951. Т. IХ. 102 с.
- Тимошенко С. П. Устойчивость упругих систем. М.; Л.: ОГИЗ-Гостехиздат, 1946. С. 133.
- Рочев А. А. Математическая модель упругопластического составного элемента несущей конструкции//Труды ПетрГУ. Сер. Прикладная мат. и информатика. Вып. 12. Петрозаводск: Изд-во ПетрГУ, 2007. С. 55-61.
- Борисов М. Д. Расчет на кручение балочных и рамных систем из тонкостенных составных стержней на планках. Л.: Стройиздат,1970. 150 с.
- Трофимов В. И. Исследование устойчивости трехгранных сквозных стержней//Исследования по стальным конструкциям: Тр. ЦНИИСК. М.: Стройиздат, 1962. Вып.13. С. 173.
- Рочев А. А. Исследование устойчивости стальных перфорированных внецентренно сжатых стержней в упругопластической стадии//Металлические конструкции и испытания сооружений: межвузовский тематический сборник трудов. Л.: ЛИСИ, 1977. №1(134). С. 119-123.
- Санжаровский Р. С. Устойчивость элементов строительных конструкций при ползучести. Л.: Изд-во ЛГУ, 1984. 280 с.
- Березин И. С., Жидков Н. П. Методы вычислений: В 2 т. М.: Наука, 1966. Т. 1. 632 с.