Расчёты на устойчивость
Равновесное состояние упругой системы называется устойчивым, если оно мало изменяется при малых возмущениях. Явление потери устойчивости. Определение величины критической силы для стержня, теряющего устойчивость в упругой стадии, по формуле Эйлера.
Рубрика | Физика и энергетика |
Вид | реферат |
Язык | русский |
Дата добавления | 08.01.2009 |
Размер файла | 37,6 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
2
2
Московский государственный технический
Университет им. Н.Э. Баумана
Калужский филиал
РАСЧЕТЫ НА УСТОЙЧИВОСТЬ
Равновесное состояние упругой системы называется устойчивым, если оно мало изменяется при малых возмущениях. Если остановится на случае прямолинейных, достаточно гибких и центрально загруженных стержней, то явление потери устойчивости состоит в следующем. При силах, не превышающих некоторой величины, называемой критической силой, прямолинейное равновесное состояние является единственным и устойчивым. Однако, если сила больше критической, то прямолинейное равновесное становится неустойчивым и стержень переходит в криволинейное равновесное состояние - изгибается. Происходит бифуркация (раздвоение) равновесных форм. Величина критической силы для стержня, теряющего устойчивость в упругой стадии, определяется по формуле Эйлера:
2EIx
Fcr = , (1)
(l)2
где:
Е - модуль упругости;
Ix - минимальный момент инерции поперечного сечения стержня;
- коэффициент приведения длины, зависящий от закрепления стержня;
l - длина стержня.
Формула Эйлера может использоваться в том случае, если потеря устойчивости происходит в упругой стадии, т.е. если критическое напряжение не превосходит предела пропорциональности:
cr= Fcr /A = 2E / 2 pr , (2)
____
где: = l / ix - гибкость стержня; ix = Ix / A радиус инерции поперечного сечения; А - площадь поперечного сечения.
В 1899 г. русским инженером Ф.С.Ясинским был предложен способ расчёта сжатых стержней на устойчивость, состоящий в том, что расчёт на устойчивость заменяется расчётом на обыкновенное сжатие, но допускаемые напряжения при этом полагаются переменными, зависящими от гибкости:
cradm = adm ; = cr / yc , (3)
здесь = коэффициент снижения допускаемого напряжения; adm - допускаемое напряжение на сжатие.
За пределами применимости формулы Эйлера, т.е. для малых значений гибкости, величины коэффициента рассчитываются с учётом возникновения упруго- пластических деформаций. Разумеется, что зависит не только от гибкости, но и от свойств материала. Для наиболее употребительных материалов составлены таблицы. Приведём такую таблицу для Ст.3, материала наиболее часто используемого для сжатых элементов конструкций.
Табл. 1
0 1.00 110 0.52
10 0.99 120 0.45
20 0.96 130 0.40
30 0.94 140 0.36
40 0.92 150 0.32
50 0.89 160 0.29
60 0.86 170 0.26
70 0.81 180 0.23
80 0.75 190 0.21
90 0.69 200 0.19
100 0.60 --- ---
Для промежуточных значений соответствующие им значения определяются путём линейной интерполяции.
Для стержней, имеющих гибкость больше 200 (редко встречающийся в практике случай), коэффициент снижения допускаемых напряжений может быть определён по формуле (3) с учётом (2).
Возможны два типа задач.
1)Задано поперечное сечение сжатого стержня - требуется найти усилие, которое может быть допущено, исходя из условия устойчивости: Nadm= = admA.
Пример 1. Определить нагрузку, которую можно допустить на ферму, исходя из устойчивости поясов АВ и СD (Рис.1). Поперечное сечение поясов выполнено из двух уголков 63; материал Ст.3, adm= 160 МПа.
F 10
B
x
1м
А С D
Рис.1
Площадь поперечного сечения А = 26,13 = 12,26см2 ( ГОСТ 8509-72). Осью, относительно которой момент инерции минимален, является ось x. Очевидно, что радиус инерции сечения относительно оси x , будет равняться радиусу инерции одного уголка ix = 1,94см (по сортаменту). Т.к. узлы фермы считаются шарнирными, то коэффициент приведения длины .Приведенная длина l = 1 2м = 200см. Гибкость = 200 / 1,94 103,1.
По таблице 1 имеем: = 0.60 - ( 0, 08/10),1 = 0,575.
Нормальная сила, которую можно допустить на стержень АВ равняется:
Nadm= admA = 0,57516кН/см212,26см2 = 112,8 кН.
Свяжем между собой силу F и усилие N.
y Fy = 0 F/2 - NSin() = 0; F = N;
N Fadm = Nadm = 112,8kH.
S
x
F/2 Рис.2
2) Второй тип задач: задана сила F - требуется подобрать размеры поперечного сечения. Можно записать: A = F /(adm ), но зависит от гибкости, следовательно, от радиуса инерции, т.е. опять от размеров поперечного сечения. Таким образом, круг замкнулся. Задача может быть решена методом попыток и сводится, по сути, к последовательности задач первого типа.
Пример 2. Подобрать размеры квадратного поперечного сечения для сжатой стойки (Рис.3). Сила F = 80кН. Материал Ст.3 с adm= 160 МПа.
Разберёмся с геометрическими харак-
теристиками сечения: А = а2, Ix = a4/12;
______ ___
F ix = Ix/A = (1/ 12 )a 0,289a.
Зададимся некоторым средним значени-
a l=0,5м ем коэффициента снижения допускаемого
напряжения 1 = 0,5. Тогда:
по табл.1 = 110 + 10(0,02/0,07) = 112,9.
Коэффициент для данного случая зак-
Рис.3 репления равняется 2.
Радиус инерции сечения ix = l / = 2,9 = 0,8857 см; сторона квадрата а = ix/0,289 = 0,8857/0,289 = 3,06см; площадь сечения А = а2 = 3,062 = = 9,36 см2. На данное сечение можно допустиь усилие:
1adm A = 0,516кН/см29,36см2 = 74,9 кН 80кН = F, т.е. сечение мало.
Пусть 2 = 0.52; = 110; ix = 100/110 = 0,909см; а = 0,909/0,289 =
= 3,15см; А = 3,152 = 9,92 см2; сечение воспринимает:
2admA = 0,52169,92 = 82,5кН, что отличается лишь на 3%.
Обычно считается, что результат достигнут, если сила, которую воспринимает сечение, отличается от силы, действующей на стержень, не более чем на 5% в ту или другую сторону:
0,95F admA 1,05F.
Округляя до более технологичного размера, примем а = 32мм. В последнем примере данных методических указаний мы покажем другой подход к организации попыток подбора, при котором образуется некоторый сходящийся итерационный процесс.
Энергетический способ определения критических сил. Изложенные выше подходы, применимы тогда, когда условия закрепления стержня и способы приложения нагрузки простейшие [1]. В более сложных случаях интегрирование дифференциального уравнения изогнутой оси стержня достаточно громоздко и целесообразно воспользоваться приближённым энергетическим способом. Рассмотрим стержень центрально сжатый силой F (Рис.4). Стержень на рисунке условно показан шарнирно опёртым, но вопрос о граничных условиях оставим пока открытым.
y
v
F F x
z dz Fп
l
Рис.4
Пусть сила F меньше эйлеровой критической силы. Если приложить к стержню некоторую малую поперечную силу Fп , то стержень изогнётся, но будет находится в устойчивом равновесном состоянии. Сжимающая сила F совершит при этом работу на перемещении , которое можно найти следующим образом. Укорочение малого элемента длиной dz будет равно:
d = dz dzCos = dz(1 Cos) = 2dzSin2( ?dz2 .
Учтём, что угол поворота равен первой производной от прогиба: v , тогда перемещение точки приложения силы найдётся:
l
? (v)2dz .
0
Потенциальная энергия изогнутого стержня выразится
l M2dz l
U = = ? EIx(v2dz , здесь учтено, что M = EIxv
0 2EI x 0
Изменение полной энергии стержня при малом изгибе будет складываться из потенциальной энергии деформации и изменения потенциала внешних сил на перемещении .
Э = U F
Если Э 0, то состояние стержня устойчиво. Если же Э 0, т.е. F U, то сила производит работу большую, чем может накопиться в стержне в виде энергии упругой деформации. Избыточная работа идёт на сообщение кинетической энергии, стержень приходит в движение и прогибается дальше, т.е. состояние его неустойчиво. Очевидно, что критическому состоянию соответствует случай
Fcr = U , или
l l
Fcr? (v)2dz = ? EIx(v)2dz , откуда получим:
0 0
l
EIx(v)2dz
0
Fcr = (4)
l
(v2dz
0
Для получения величины критической силы необходимо задаться формой изогнутой оси v = v(z), удовлетворяющей граничным условиям задачи. С математической точки зрения (4) является функционалом, т.е.отображением из множества функций определённого класса (дважды дифференцируемых и удовлетворяющих граничным условиям) в множество действительных чисел.
Пример 3. Найдём критическую силу для стержня, шарнирно опёртого
по обоим концам (Рис.5). Точное решение в этом
F случае известно: Fcr = 2EIx / l2 9,8696EIx / l2.
Если задаться функцией v = CSin(z/l), то
получим точное решение. Допустим, что мы
l этого не знаемем и попробуем аппроксимировать
z изогнутую ось полимиальной функцией. Возьмём
для начала полином второй степени:
Рис.5
v = Az2 + Bz + C (5)
Запишем граничные условия: 1) при z = 0: v = 0; 2) при z = l: v = 0. Подставляя в (5), получим:
С = 0; Al2 + Bl = 0 B = Al.
Дифференцируя (5) и учитывая полученные выражения для коэффициентов, имеем:
v = A(2z - l); v = 2A.
Подставив в (4), имеем:
l
4A2 dz
0
Fcr = = 12EIx / l2.
l
A2 (2z - l )2dz
0
Мы видим, что полученное значение критической силы отличается от точного решения, более чем на 20%.
Известно, что чем выше степень аппроксимирующего полинома, тем выше точность решения. Аппроксимируем иэогнутую ось стержня полиномом четвёртой степени:
v = Az4 + Bz3 + Cz2 + Dz + E (6)
при z = 0: v = 0;
при z = 0: M = EIxv v = 0;
при z = l : v = 0;
при z = l : M = 0 v = 0.
Возьмём производные от (6):
v = 4Az3 + 3Bz2 + 2Cz + D;
v = 12Az2 + 6Bz + 2C.
Реализуем граничные условия, получив при этом систему из четырёх алгебраических уравнений.
Е = 0; 2) С = 0; 3) Al4 + Bl3 + Dl = 0; 4) 12Al2 + 6Bz =0 B = 2Al , подставляя это в предыдущее уравнение, имеем: D = Al3.
Подставляя это в выражения для производных, получим:
v = A(4z3 6lz2 + l3); v = 12A(z2 - lz).
Подставив в (4),будем иметь:
l
144A2EIx (z2 - lz)2dz
0
Fcr = = 168EIx / (17l2) 9,8824EIx / l2 .
l
A2 (4z3 - 6lz2 + l3)2dz
0
Как видим, полученное решение практически совпадает с точным. Обратим внимание на тот факт, что приближённые решения всегда дают завышенные значения критических сил. Это происходит по той причине, что в приближённом решении стержень система с бесконечным числом степеней свободы, заменяется более жёсткой системой с конечным числом степеней свободы.
Пример 4. Найти критическую нагрузку для стержня, показанного на
Рис.6. В этом случае значение коэффициента
приведения длины неизвестно и нет возмож-
F ности непосредственно использовать для вы-
числения критической силы формулу (1). При-
меним энергетический способ. Для аппрок-
l симации изогнутой оси стержня используем
выражение (6). Граничные условия будут
выглядеть:
при z = 0: v = 0;
при z = 0: v = 0 (в заделке угол поворота
z 2l равен нулю);
3) при z = l: v = 0;
4) при z = 3l: M = 0 v = 0 ( на верхнем
конце стержня изгибающий момент равен
Рис.6 нулю);
Выражения для производных см. стр.8. Воспользуемся граничными условиями:
Е = 0; 2) D = 0; 3) 4l2A +2lB + C = 0; 4) 54l2A + 9lB +
+ C = 0, решая относительно А систему из двух последних уравнений: 9lB + C = 54l2A,
2lB + C = 4lA ,
получим: B = ( 50/7)Al; C = (72/7)Al.
Первая и вторая производные от прогиба запишутся:
v = 2A[2z3 - (75/7)lz2 + (72/7)l2z], v = 12A[z2 - (25/7)lz + (12/7)l2].
Подставив в (4), вычислим критическую силу
3l
144A2EIx [z2 - (25/7)zl + (12/7)l2]2dz
0
Fcr = 1,09EIx / l2 .
3l
4A2 [2z3 - (75/7)z2l + (72/7)zl2]2dz
0
Пример 5. Для сжатой стойки, показанной на Рис.7, используя энер-
гетический способ определить коэффи-
z циент приведения длины . Подобрать
размеры поперечного сечения стойки. Си-
ла F =200кН. Материал Ст.3: adm=160 МПа.
D F d/D = = 0,8; l = 3м. Использовать методи-
ку расчёта по коэффициенту снижения до-
d l пускаемых напряжений.
0,6l Решение
z Запишем граничные условия:
при z =0: v = 0;
при z = 0: v = 0;
при z = l : v = 0;
. при z = l : v = 0
Рис.7 Подставим это в аппроксимирующий полином (6) и во вторую производную от него. В результате получим те же самые выражения для v и v, что и в примере 3. Различие будет состоять лишь в том, что в выражении для критической силы интеграл, стоящий в знаменателе придётся брать в пределах от 0 до 0,6l. Это легко понять, если вспомнить, что знаменатель в формуле (4) представляет собой удвоенное перемещение точки приложения силы, а оно зависит от укорочения части стержня, лежащей ниже сечения, в котором приложена сила.
l
144A2EIx (z2 - lz)2dz
0
Fcr = 19,53EIx / l2 .
0,6l
A2 (4z3 - 6lz2 + l3)2dz
0
Найдём коэффициент приведения длины. Для этого представим выражение для критической силы:
Fcr = 19,53EIx / l2 = 2EIx / (2/ 19,53)l2.
Сопоставляя полученный результат с формулой Эйлера (1), получим:
2 = 2 / 19,53 0,711.
Приведенная длина стержня l = 0,711300 213см.
Подберём размеры поперечного сечения. Определим геометрические характеристики. Площадь сечения:
A = D2(1 2)/4 0,785D2(1 - 0,82) 0,283D2.
Момент инерции:
Ix = D4 (1 4) / 64.
Радиус инерции:
_____ ______ _______
ix = Ix / A = 0,25 D 1 + 2 = 0,25D 1 + 0,82 0,32D.
Для подбора размеров сечения используем следующий сходящийся алгоритм: зададимся некоторым средним значением коффициента :
1 = 0,5, тогда площадь А = F / 1adm = 200kH/ (0,516kH/ см2 =25см2;
_______ _______
диаметр сечения равен D = A/0,283 = 25/0,283 = 9,4см;
радиус инерции ix = 0,32D = 0,329,4 = 3,01см;
гибкость = l / ix = 213/3,01 = 70,8;
коэффициент приведенной длины, соответствующий этой гибкости
1 = 0,81 - (0,06/10)0,8 = 0,8052 (табл.1);
следующее значение примем равным среднему арифметическому из двух предыдущих:
2 = ( + 1)/2 = (0,5 + 0,8052)/2 0,653 и повторим расчёт:
__________
А = 200/(0,65316) = 19,14см2; D = 19,14/0,283 = 8,22см; ix = 0,328,22
= 2,63см; = 213/2,63 = 81; 2 = 0.75 (0,06/10) = 0,744;
= (0,653 + 0,744)/2 = 0,7; А = 200/(0,716) = 17,86см2;
__________
D = 17,86/0,283 = 7,94см; ix = 0,327,94 = 2,54см; = 213/2,54 = 83,9
= 0,75 - (0,06/10)3,9 = 0,727; разница между двумя соседними значениями составляет:
100(0,727 - 0,7)/ 0.7 = 3,9% 5%.
Расчёт можно считать законченным. Требуемый размер D = 7,94см, но из конструктивных соображений примем D = 80мм.
Список литературы
Феодосьев В.И. Сопротивление материалов. - М.: Наука. - 1986. 512с.
Работнов Ю.Н. Механика деформируемого твердого тела. - М.: Наука. - 1979. - 744с.
Подобные документы
Понятие равновесного состояния, его виды. Пределы применимости формулы Эйлера. Влияние условий закрепления концов стержня на величину критической силы. Понятие коэффициента запаса на устойчивость. Энергетический способ определения критических сил.
курс лекций [888,8 K], добавлен 23.04.2009Определяющие соотношения модели нелинейно упругой среды, вычисление компонент тензора напряжений. Определение автомодельного движения. Сведение модельных соотношений к системе дифференциальных уравнений. Краевая задача разгрузки нелинейно упругой среды.
курсовая работа [384,1 K], добавлен 30.01.2013Понятие устойчивости применительно к электрической системе. Определение взаимных и собственных проводимостей при различных системах возбуждения, определение коэффициента запаса статической устойчивости. Расчёт динамической устойчивости данной системы.
курсовая работа [403,9 K], добавлен 26.01.2011Построение векторных диаграмм неявнополюсного и явнополюсного генераторов. Запас статической устойчивости простейшей электрической системы, а также меры по её повышению. Критерии статической устойчивости. Внутренняя реактивная мощность генератора.
контрольная работа [287,7 K], добавлен 19.08.2014Уравнения для поперечных компонент смещения плазмы, минимизация функционал Крускала-Обермана потенциальной энергии МГД-возмущения. Невозмущенное состояние, потенциальная энергия возмущения. Преобразование кинетического слагаемого, условие устойчивости.
реферат [567,9 K], добавлен 22.07.2011Расчет и анализ установившихся режимов схемы электроэнергетической системы (ЭЭС). Оценка статической устойчивости ЭЭС. Определение запаса статической устойчивости послеаварийного режима системы. Отключение сетевого элемента при коротком замыкании.
курсовая работа [563,4 K], добавлен 11.09.2015Расчет параметров схемы замещения в относительных единицах. Определение электродвижущей силы генератора и соответствующих им фазовых углов. Расчет статической устойчивости электрической системы. Зависимость реактивной мощности от угла электропередачи.
курсовая работа [941,9 K], добавлен 04.05.2014Расчет электромеханических переходных процессов в системе электропередачи. Предельное снижение напряжения на шинах асинхронного двигателя. Оценка статической и динамической устойчивости системы. Аварийный и послеаварийный режимы при коротком замыкании.
курсовая работа [1,2 M], добавлен 10.12.2014Описание принципа действия системы автоматического регулирования (САР) для стабилизация значения давления газа в резервуаре. Составление структурной схемы с передаточными функциями. Определение запасов устойчивости системы по различным критериям.
дипломная работа [4,6 M], добавлен 22.10.2012Повышение динамического качества станков с помощью возмущений подшипников качения. Колебания при отсутствии вынуждающей силы и сил вязкого сопротивления. Незатухающие гармонические вынужденные колебания. Нарастание амплитуды во времени при резонансе.
реферат [236,6 K], добавлен 24.06.2011