150815 (732819), страница 2
Текст из файла (страница 2)
Подставив в (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с.