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

Размещено на http://www.allbest.ru/

Л.А. Барагунова, М.М. Шогенова

Кабардино-Балкарский государственный университет, г. Нальчик

Аннотация

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

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

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

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

На рисунке 1 изображён однородный упругий стержень переменного сечения, на который действует переменная распределённая нагрузка q(x) и сосредоточенные силы Fi, действующие вдоль оси. Внешние нагрузки считаются «мёртвыми», т. е. не меняются ни по величине, ни по направлению при деформировании стержня. Так же заданы функции изменения переменной жёсткости EJ(x) и распределённой нагрузки q(х).

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

, x [0, l], B(x) = EJ(x) (1)

В правой части уравнения (1) изгибающий момент определяется из условия равновесия отсечённой верхней части стержня. Например, по рис. 2

. (2)

При учёте (2) уравнение (1) становится интегро-дифференциальным

, x [0, l] (3)

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

 x [0, l] (4)

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

Уравнение (4) имеет очевидное тривиальное решение, что соответствует простому сжатию, что в данном случае исследуется. Задача состоит в том, чтобы определить критические силы, которым могут соответствовать ненулевые решения, т.е. искривлённые положения равновесия.

Отыскание аналитических решений в таких и других более сложных задачах практически невозможно. Выход состоит в использовании численных методов. Используя метод конечных разностей и разобьём длину стержня на n одинаковых отрезков с шагом h = l/n, c номерами узловых точек i = 1, 2 , … , n, n + 1.

Вместо непрерывной функции непрерывного аргумента v(x) введём сеточную функцию Тогда производные в уравнении (4) можно представлять приближённо в виде конечноразностных соотношений

в силу чего оно примет вид

(5)

При этом к левой части (4) процедура замены второй производной конечноразностным соотношением применена дважды. Коэффициенты уравнения имеют значения

Здесь

Ni = N(xi).

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

Верхний конец (х = 0) свободен и к нему приложена сила F. Поэтому:

1) Изгибающий момент равен нулю.

(6)

2) Поперечная сила определяется с помощью функции прогибов

(7)

По рис. 2 с учётом малости угла поворота концевого сечения

(8)

Приравнивая правые части (7) и (8), получим

Конечноразностная аппроксимация производных даёт уравнение

(9)

Где b = -3B1+4B2-B3, c21=2b - 5B1-3Fh2, c22= -5b + 18B1+ 4Fh2,

c23 = 4b - 24B1 - Fh2, c24 = -b + 14B1, c25 = -3B1.

Правый конец (х = l) защемлён, поэтому угол поворота и перемещение равны нулю

(10)

4) (11)

Уравнения (5), (6), (9)-(11) образуют линейную однородную алгебраическую систему

(12)

где y = {y1, y2, …, yn+1} - вектор, компонентами которого являются отклонения стержня, С-квадратная матрица порядка n + 1

Нулевые элементы матрицы здесь не выписаны.

Критические значения параметров нагрузок, т.е. собственные значения определяются из уравнения

(13)

которое является условием существования нетривиального решения уравнения (12). Решение уравнения аналитическими методами не представляется возможным, поэтому используем приближённые решения с помощью численных и графических методов. Численные методы приводят к громоздким компьютерным программам и реализуются с трудом. Графический метод является более предпочтительным. Данный метод основан на возможностях моментальной визуализации графика левой части уравнения (13) с помощью современных компьютерных программ MatLab, MathCad и т.д. Для этого в координатной системе F - det(C) на экране строится соответствующая кривая. Точки пересечения ею оси F и определяют значения критических сил. Имеется возможность увеличения любых фрагментов рисунка, что легко позволяет достигать высокой степени точности, и отличает компьютерный графический способ от обычного ручного способа небольшой точности.

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

l = Е = 1, J = 1, q

При этом получены ожидаемые первые критические силы Fк = 1, 9, 25 с высокой степенью точности, что позволяет уверенно переходить к задачам более общего характера.

Пример. Рассмотрим стальной стержень, имеющий поперечное сечение с переменным моментом инерции

J(x) = 30(1+20x/l) см4.

К стержню приложена сосредоточенная сила F и неравномерно распределённая нагрузка

q(x) = (l - x) Н/м,

ассоциированная с сосредоточенной силой так, что будет отыскиваться критическое значение только одного параметра F. Прочие исходные данные примем следующими: l = 3 м, Е = 200 ГПа.

Результат счёта, выданный на экран монитора, показан на рис. 3 в виде графика. После увеличения рисунка легко читаются первые элементы спектра собственных значений Fк = {168,9 849,2} кН.

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

Литература

1. Алфутов Н.А. Основы расчёта на устойчивость упругих систем. - М.: Машиностроение, 1978. - 312 с.

2. Вольмир А.С. Устойчивость деформируемых систем. М.: Наука. 1967. - 984 с.

3. Варвак П.М., Варвак Л.П. Метод сеток в задачах расчёта строительных конструкций. М.: Стройиздат. 1977. -154 с.

4. Культербаев Х.П., Барагунова Л.А. О реализации проблемы собственных значений сжато-растянутого стержня на компьютере. Компьютерные технологии в строительстве: Материалы Всероссийской научно-технической конференции. ДГТУ. - Махачкала: Алеф (ИП Овчинников), 2012. С. 90-94.

5. Культербаев Х.П. О структурировании пространства параметров сжато-растянутого стержня по механическому состоянию // Известия вузов. Северо-Кавказский регион. Технические науки. 2009. № 3. С.85-88.

6. Литвинов С.В., Языев Б.М., Бескопыльный А.Н., Ананьев И.В. Расчёт на устойчивость стержней из ЭДТ-10 при начальной погиби стержня в виде S-образной кривой. // Инженерный вестник Дона, 2012, №1.

7. Барагунова Л.А. Устойчивость предварительно сжимаемой арматуры в железобетонных балках. // Инженерный вестник Дона, 2016, №4.

8. Sargin M. Stress-strain relationships for concrete and the analysis of structural concrete sections. SM study, №4, Solid Mechanical Division, University of Waterloo. Ontario, Canada. - 1970. p. 167.

Размещено на Allbest.ur