Моделирование динамического поведения усиливаемых стержневых систем

УДК 539.3

Ф. С. Хайруллин

МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКОГО ПОВЕДЕНИЯ УСИЛИВАЕМЫХ СТЕРЖНЕВЫХ СИСТЕМ

Ключевые слова: стержневая система, усиление конструкции, вариационный метод, вынужденные колебания,

аппроксимирующая функция.

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

Keywords: rods system, structural reinforcement, variational methods, forced oscillations, approximating functions.

Variational method is proposed for solving problems offorced oscillations of rod systems that are pre-exposed to static loads before amplification, and after amplification. At each stage of solving the problem of deformation of the system was taken into account, caused loading in the previous stages. Numerical results are given.

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

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

динамическим нагрузкам. Для данных конструкций определяется напряженно-деформированное

состояние в случае вынужденных колебаний.

Решение задачи разбивается на три этапа:

1) Определение напряженно-деформированного состояния первоначальной (ослабленной) стержневой системы от действия статических ремонтных нагрузок.

2) Определение напряженно-деформированного состояния усиленной стержневой системы от действия дополнительных статических нагрузок с учетом ремонтных напряжений.

3) Определение напряженно-деформированного состояния усиленной стержневой системы от действия динамических нагрузок с учетом ремонтных и дополнительных статических напряжений.

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

Для решения задачи на первых двух этапах используется вариационный принцип Лагранжа [4]. На первом этапе вариационное уравнение

записывается так:

5пр-5'WР = 0,

где П - потенциальная энергия деформации первоначальной стержневой системы при действии ремонтных нагрузок; 5'МР - вариация работы

ремонтных нагрузок.

С учетом принятых гипотез

5Пр = | |(стр15бр1 + стр25ер2 + стр35ер3 )САС,

ЬЛр

5^" =| (qp5up + qp5up + qp5up ^ +

L

/

+ М/25ф2Р/ + М/35фр/).

Здесь ст/1, ст/2, ст/3, е/1, е/2, е/3 - ремонтные напряжения и деформации в произвольной точке стержня с координатами 5, у, г; ир, ир, из,

фр, фр, фр - перемещения и углы поворота оси стержня, вызванные ремонтными нагрузками; д/, д/, д/ - компоненты интенсивности

распределенной нагрузки; з - компоненты

сосредоточенной силы, М/,М/2,М/3 - компоненты

сосредоточенного момента, действующих в точке с координатой 5,-; 5 - длина дуги оси стержня; у, г -главные центральные оси инерции поперечного сечения стержня; Ар - длина и площадь

поперечного сечения стержня.

На втором этапе после усиления конструкции вариационное уравнение имеет вид:

5Пc-5 ' Wc =0,

где

5ПС =Ц[(а Р +ai1 )5в?1 +(ap2 +а^2 )5е?2

LA

+

+ кэ )seÎ3 \lAdi, SWc = j[(q? + q? )>uC +(q2p + q2C +

L

+ + q3)su3 d + z F + FC)suf + /

+ F/S + F2 jSuC + (f? + F3 Цс'|+

+ z[(M? + М^фС +(м?2 + Mc2 >ф? + + (M,?3 + мсэ )>фЭ |.

Здесь индексом «с» отмечены величины, относящиеся ко второму этапу, когда в результате усиления площадь поперечного сечения стержня стала равна A и на нее действуют дополнительные статические нагрузки. При этом оси y, z могут не быть главными центральными осями инерции измененного поперечного сечения стержня. Также при расчетах учитывается, что в результате действия нагрузок на первом этапе ось стержня подверглась деформации и изменились её геометрические характеристики.

Для решения задачи на третьем этапе используется вариационный принцип

Остроградского-Гамильтона [4], на основании которого истинное движение системы за отрезок времени [t0, t1] удовлетворяет вариационному уравнению:

+SKд -S'Wд)dt = 0,

10

где

snd = jj [(ст ? + ст^ + ct|1 )Sedi +(ст ?2 + ст^ +

LA

+ ст?2 )se?2 +(ст?э +стСэ +стдэ )ss?3 |dAdi, sk д=jj p^sù«1+u;dsu;d+u;dsu;d )dAdi+

LA

+ 2 Z mk (uÎSuÎ + u2kSu2k + uIkSuIk ),

2 k=1

SWд = j [(q? + q1 + qd )su| +(q? + qc2 +qI )sul +

L

+(q?+q3 + ql )su| |di+z[(f,? + fc + FI )su? + /

+ (f?+FC2+F2 }sud2i + (f§ + FC3 + F3 )s u|' |+

+z[(M? + МС1 + М,д)>фд'' + (M?2 + МС2 +

+m д )5ф2+(м?3+мс3+M «3 ^ |,

Kд - кинетическая энергия стержневой системы; «д «д *д

Ui , U2 , и3 - перемещения точки стержня с координатами s, y, z ; u j - частная производная от перемещения u j по времени t; p - плотность материала стержня; индексом « д » отмечены величины, относящиеся к третьему этапу, когда на стержень действуют динамические нагрузки.

На основании принятых гипотез можно получить следующие формулы, определяющие

перемещения и углы поворота точек поперечного сечения стержня:

и;

ф; (5, у, г, г ) = ф5, г) ф2 (5, у, г, г ) = ф2 (5, г), фЭ (5, у, г, г) = фэ (5, г) где и; = и , и2, и3, ф1, ф2, фЭ| - вектор компонентов перемещения и;, и2, иЭ и углов поворота ф1, ф2, фЭ произвольной точки стержня;

и = , и2, из, Ф1, Ф2, Фэ} - вектор компонентов перемещения и углов поворота оси стержня.

Деформации стержня определяются через деформации оси стержня по формулам 8ц = 8 + г -к 2 - у -Кз ,

е12 = У3 - г К1, 813 = У 2 + у К1. (5) Деформации оси стержня находятся из соотношений Клебша:

8 = иг5 - кз и2 + к2 из,

У2 = -из,э + к2 и1 -к1 и2 -Ф2,

У3 = и2^ - к из + кз и -Фз,

К1 = Ф1,5 - кз Ф2 + к2 Фз, К2 = Ф2,я - к1 Фз + кз Ф1, Кз = Фз,я - к2 Ф1 + к Ф2, где 8, у2, уЭ - продольная деформация и углы сдвига; К1, К2, К3 - изменения кручения и кривизн (деформации кручения и изгиба) оси стержня; к , к2, кз - кручение и кривизны оси стержня,

На основании обобщенного закона Гука и введенных гипотез получаются следующие физические соотношения, связывающие напряжения и деформации:

0*11 = Е 8ц, 012 = 6 812 , (7) °13 = 6 81Э,

где Е, 6 - модули упругости и сдвига стержня.

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

Для решения задачи стержневая система разбивается на элементы X, (рис.1), оси которых

являются гладкими кривыми.

Введем для всей конструкции глобальную систему координат у, г и для каждого стержня X I локальные системы координат 5/, у, г, где координатная линия 5, направлена вдоль оси стержня, линии у ¡, г , совпадают с главными

центральными осями инерции поперечного сечения первоначального не усиленного стержня.

Рис. 1

Для каждого стержня X, за искомые неизвестные принимается вектор перемещений

и, = ,и3, у1 >Ф2'Фэ} , заданный в глобальной системе координат у, у, ф. На каждом этапе этот вектор обозначается соответствующим индексом. Решение ищется в виде [5,6]: Ф м

и,

т=1

где вт(()=' 1, вт 2, вт э, вт 4, вт 5, вт 6} -

вектор неизвестных постоянных на первых двух этапах и вектор неизвестных функций от времени t

на третьем этапе; р, = -

-, 0 <р, < 1; 5,12

в,- 2 - 5,1

- координаты начала и конца оси стержня X .

Функции формы /т (р,-) записываются так:

/1(Р1)= 1 -Р1, /2 (Р1) = Р1, _

/т

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

Если определить значения искомой функции в узловых точках, то в соответствии с

функциями формы (9) имеем Ц (Эц^) = В{ (),

и, (Э/2, t) = В2 ((). Следовательно, коэффициенты

В1 ((), В'2 (() определяют значения вектора перемещений соответственно в начальной и конечной точках оси стержня X . Это

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

положить В{ = В2, то будет выполнено условие стыковки начала , -го стержня с концом } -го стержня. Чтобы удовлетворить, например, граничному условию жесткого защемления на конце

, -го стержня, т.е. условию и,(Эц^) = 0, следует положить В{ (() = 0 .

Компоненты перемещения стержня X в

локальной системе координат связаны с компонентами перемещения в глобальной системе координат соотношениями:

и, = [с, ] ц,

где

[с, ] =

[со ] 0

[с

Со ] =

ю 11 ю с12 Г ю с13

ю 21 ю с22 ю с23

ю э1 ю сэ2 ^ ю с33

[с о ] - матрица направляющих косинусов локальной системы координат в глобальной системе координат.

Подставляя для каждого стержня X

аппроксимирующие функции (8) в равенства (10), определяя затем перемещения, деформации и напряжения стержней через перемещения осей стержней по формулам (4) - (7), подставляя полученные выражения для каждого этапа в вариационные уравнения (1) - (3), после соответствующих преобразований получаются системы уравнений относительно неизвестных параметров В'т ^) т = 1, М, \ = 1,/ :

[рр ¿л }={рр },

рС В }={РС },

[рд В}- [мд ]Вд }=^д},

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

уравнений.

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

Матричные уравнения

используется метод Ньюмарка [7].

В качестве примера приведены результаты расчета усиленной рамной конструкции, показанной на рис.2,а. До усиления поперечные сечения всех стержней рамы одинаковые (двутавр № 18) и действует первоначальная (ремонтная) нагрузка Fp = 10 кН. Усиление производится путем увеличения поперечного сечения левой стойки в два раза с помощью такого же стержня двутаврового сечения (рис.2,б). После этого на раму действует сначала статическая нагрузка Fc = 15 кН, а затем постоянная динамическая нагрузка Fд = 5 кН .

0

а)

600 см

77777

f

б)

77777

Рис. 2

На рис.3 приводится графики изменения по времени максимального напряжения для усиленной рамы (линии 1) и для усиленной рамы (линии 2), у которой не учитываются ремонтные напряжения, предполагается, что вся нагрузка действует после усиления: Рр = 0, Рс = 25 кН .

250

200

150

1 Г~

---- -----

-- --- — ^Сх --_ ----

-А»

0.25

0,5 0,75

Рис. 3

t. сек

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

Рс = 20 кН , Р6 = 0 . Статическое решение совпадает с решением, приведенным в работе [8].

Результаты показывают, что учет напряжений, действующих в конструкции до усиления, приводит к увеличению напряжений, появляющихся в конструкции после усиления.

Литература

3

Бондаренко С.В. Усиление железобетонных конструкций при реконструкции зданий / С.В. Бондаренко, Р.С. Санжаровский. - М.: Стройиздат, 1990. - 352 с.

Будин А.Я. Усиление портовых сооружений / А.Я. Будин, М.В. Чекренева. - Л.: Транспорт, 1983. - 178 с.

Ребров И.С. Усиление стержневых металлических конструкций/ И.С. Ребров. - Л.: Стройиздат, 1988. - 288 с.

4. Абовский Н.П. Вариационные принципы теории упругости и теории оболочек / Н.П.Абовский, Н.П.Андреев, А.П.Деруга - М.: Наука, 1978. - 288 с.

5. Серазутдинов М. Н. Метод расчета криволинейных стержней / М.Н.Серазутдинов, Ф.С.Хайруллин // Известия вузов. Строительство и архитектура. - 1991. -№5. - С. 104 - 108.

6. Хайруллин Ф.С. Моделирование динамики стержневых систем с движущимися массами. / Ф.С. Хайруллин, Х.А. Абрагим, О.И. Алхалили // Вестник Казанского технологического университета. № 8. Казань. 2013 г. С.133-136.

7. Бате К. Численные методы анализа и метод конечных элементов / Бате К., Вилсон Е. - М.: Стройиздат, 1982. -448 с.5.

8. Серазутдинов М.Н. Особенности расчета усиливаемых упругопластических рамных систем с учетом продольных сил/ Серазутдинов М. Н., Убайдуллоев М.Н. // Вестник Казанского технологического университета. № 9. Казань. 2010 г.- С.523-529.

© Ф. С. Хайруллин - д-р физ.-мат. наук, доц., проф. каф. теоретической механики и сопротивления материалов КНИТУ, tmsm@kstu.ru.