Сопротивление материалов курс лекций
..pdfэнергию деформации с достаточной для инженерных расчётов точностью можно вычислять с учётом только изгибающего момента:
U =∫ |
M x2dz |
. |
(11.7) |
|
|||
l 2EIx |
|
Определение перемещений с помощью теоремы Кастилиано
Для упругих систем, подчиняющихся закону Гука, Кастилиано можно сформулировать следующим образом:
∂U
∂Fi =∆i .
теорему
(11.8)
Здесь U – потенциальная энергия системы; Fi – некоторая обобщённая сила; ∆i – перемещение точки приложения силы по её
направлению, т.е. соответствующее обоб- |
|
||
щённое перемещение. |
|
|
|
Формула |
(11.8) позволяет |
найти |
|
обобщённое перемещение, соответствую- |
|
||
щее любой из приложенных к системе со- |
Рис. 11.4. |
||
средоточенных |
обобщённых сил. |
Напри- |
|
мер, для показанной на рис. 11.4 балки прогиб в точке А определяется следующим образом.
Запишем выражение потенциальной энергии
U =∫ |
M x2dz |
=∫ |
(−Fz)2 dz. |
|
|||
l 2EIx l |
2EIx |
Здесь учтено, что изгибающий момент в произвольном сечении z M x =−Fz .
После интегрирования получаем
U = F 2l3 . 6EIx
151
Взяв частную производную от этого выражения по величине F, находим
v |
|
|
|
|
∂U |
|
Fl3 |
|
A |
=∆ |
F |
= |
∂F |
= |
|
. |
|
|
||||||||
|
|
|
|
3EIx |
Несмотря на очевидную простоту, этот метод обладает существенным недостатком – он позволяет определять перемещения только тех точек, где приложены сосредоточенные нагрузки, причём перемещения должны соответствовать этим нагрузкам. Так, определить угол поворота в сечении А для рассмотренного примера не удаётся, поскольку углу поворота как обобщённому перемещению соответствует сосредоточенный момент, а в нагрузке он отсутствует. Указанные ограничения можно обойти с помощью фиктивных нагрузок. Для этого в интересующей точке конструкции прикладывают обобщённую силу FФ, соответствующую искомому перемещению. Потенциальная энергия системы выражается через заданные нагрузки и эту фиктивную силу, затем полученное выражение дифференцируется по величине FФ, а в конечном результате полагают FФ = 0.
Принцип возможных перемещений. Интеграл Мора
Известный из курса теоретической механики принцип воз-
можных перемещений Лагранжа применительно к деформиру-
емым системам можно сформулировать следующим образом: если система под действием приложенных нагрузок находится в равновесии, то сумма работ внешних и внутренних сил на возможных бесконечно малых перемещениях точек системы равна нулю:
W +W1 =0. |
(11.9) |
Здесь W – работа внешних обобщённых сил Fi на их возможных обобщённых перемещениях ∆Fi, W =ΣFi∆Fi ; W1 – работа внутренних сил.
152
Внутренние силы препятствуют возникновению деформаций от нагрузки, поэтому работу этих сил следует считать отрицательной. Возможными называют перемещения, которые не противоречат наложенным на систему связям, т.е. учитывают условия закрепления конструкции.
Возможные перемещения добавляются к действительным, возникающим от приложенных нагрузок. При этом работа внутренних сил системы на этих дополнительных возможных перемещениях
−W =∫Mкdφ+∫M x dθx +∫M y dθy +
l |
l |
l |
(11.10) |
|
+∫N∆dz +∫kxQx dsx +∫kyQy dsy . |
||||
|
||||
l |
l |
|
|
Силы в процессе бесконечно малого изменения перемещений считаются неизменными, поэтому множитель 1/2 в выражении работ отсутствует.
Для нахождения перемещений в произвольной точке конструкции поступим следующим образом. Рассмотрим два состояния нашей системы: первое состояние, которое возникает от действия приложенных нагрузок, будем называть грузовым, все относящиеся к нему величины – внутренние силовые факторы и соответствующие им обобщённые перемещения пометим индексом «F». Второе со-
стояние будет вызываться фиктивной обобщённой силой F =1 (при нахождении углового перемещения в роли такой силы будет выступать единичный безразмерный момент). Это состояние будем называть единичным. Все величины, относящиеся к единичному состоянию, пометим индексом «1».
Сформулируем принцип возможных перемещений для системы, находящейся в единичном состоянии, т.е. нагруженной фиктивной единичной силой в той точке, где нас интересует перемещение. При этом в качестве возможных перемещений возьмём обобщённые перемещения грузового состояния, т.е. действительные перемещения системы, возникшие от приложенных к ней нагрузок. Эти перемещения хотя и не являются бесконечно малыми, как того требует
153
принцип Лангранжа, но по сравнению с размерами конструкции они малы и удовлетворяют условиям закрепления.
В соответствии с (11.1) работа фиктивной единичной силы будет равна искомому перемещению:
W = |
F |
∆1F =1 ∆1F =∆1F . |
(11.11) |
Обобщённые перемещения грузового состояния системы, входящие в выражение (11.10), имеют вид:
dφ= |
M Fкdz |
; ∆dz = |
NF dz |
; |
|
||||||||
|
|
GIк |
EA |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
dθx = |
M Fхdz |
; dsx = |
QFx ds |
|
; |
(11.12) |
|||||||
|
|
|
|
||||||||||
|
|
|
EIx |
|
|
|
|
|
GA |
|
|
|
|
dθy = |
M Fy dz |
; dsy = |
QFy ds |
. |
|
||||||||
|
|
|
|||||||||||
|
|
|
EIy |
|
|
|
|
|
GA |
|
|
|
Умножая внутренние усилия от единичной нагрузки на соответствующие обобщённые перемещения (11.12) и интегрируя по длине всех стержней конструкции, запишем (11.9) следующимобразом:
∆1F =∫ |
M |
1к |
M |
Fк |
+∫ |
M |
1x |
M |
Fx |
+∫ |
M1y M Fy |
dz + |
|
||||||
|
|
GIк |
|
|
EIx |
|
EI y |
|
|
||||||||||
|
l |
|
|
|
|
l |
|
|
|
l |
|
|
(11.13) |
||||||
|
N N |
|
|
|
|
|
kyQ1xQFx |
|
|
|
kyQ1yQFy |
|
|
||||||
+∫ |
F |
dz +∫ |
dz +∫ |
dz. |
|
||||||||||||||
1 |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
l |
EA |
|
|
l |
|
GA |
|
|
|
|
l |
|
GA |
|
|
|
Полученное выражение называется интегралом Мора. Оно может быть использовано для определения перемещений в любых произвольно нагруженных упругих стержневых системах. Как говорилось ранее, не все интегралы в правой части (11.13) имеют одинаковый порядок величины. При наличии в поперечных сечениях стержней изгибающих и крутящего моментов последние три слагаемые можно не учитывать.
154
В случае простых видов нагружения интеграл Мора записывается проще. Например, для системы, работающей на растяжении и сжатие, перемещение определяется следующим образом:
|
∫ |
EA |
|
|
∆i = |
|
N1i NF |
dz. |
(11.14) |
|
|
l
Выражение (11.14) нужно понимать как сумму интегралов по длине каждого из стержней системы. Здесь ∆i – перемещение точки
приложения единичной фиктивной силы F i = 1 в направлении этой силы; Ni – нормальные силы в стержнях от действия единичной силы; NF – нормальные силы в стержнях системы в грузовом состоянии, т.е. от действия приложенных нагрузок.
При кручении перемещения представляют собой углы закручивания:
φi =∆i =∫ |
M1кM Fк |
dz. |
(11.15) |
|
|||
l |
GIк |
|
В качестве обобщённых фиктивных сил в этом случае высту-
пают единичные скручивающие моменты Мi = 1, поочерёдно прикладываемые в тех сечениях, где требуется найти угловое перемещение.
При плоском изгибе перемещения с достаточной точностью можно найти, учитывая в выражении (11.13) только изгибающий момент:
∆i =∫ |
M1x M Fx |
dz. |
(11.16) |
|
|||
l |
EIx |
|
Индексы «F» в обозначении изгибающих моментов от нагрузки обычно опускают, так же как и индекс «х» для обозначения моментов единичного состояния. Для определения прогиба в интересую-
щей точке прикладывают сосредоточенную силу Fi = 1, для нахожде-
ния угла поворота прикладывают сосредоточенный момент Мi = 1. Если имеется несколько силовых участков, на каждом из которых
155
уравнения изгибающих моментов М1i(z) и (или) Мх(z) имеют различный вид, интегрирование ведётся по каждому из участков отдельно, а результаты алгебраически суммируются:
m |
|
M1i M x |
|
|
||
∆i =∑ |
∫ |
dz j , |
(11.17) |
|||
|
||||||
j=1 |
EI |
x |
|
|||
|
l j |
|
|
где m – число силовых участков; lj – длина j-го силового участка. Необходимо обратить внимание на то, чтобы выделение сило-
вых участков и выбор локальных координат zj на них было одинаковым как в грузовом, так и в единичном состояниях конструкции, т.е. внутренние усилия, стоящие под знаком интеграла, должны быть функциями одной и той же переменной.
Положительный знак результата свидетельствует о том, что перемещение происходит в направлении единичной обобщённой силы.
Пример 11.1. Определение перемещений интегралом Мора. Для шарнирно опёртой балки постоянной жёсткости, находя-
щейся под давлением равномерно распределённой нагрузки интенсивностью q (рис. 11.5, а), определить прогиб в середине пролёта и угол поворота сечения над опорой В.
Решение. Для определения прогиба в сечении С освобождаем балку от действия распределённой нагрузки и прикладываем в этом
сечении фиктивную сосредоточенную силу F = 1. В единичном состоянии балка имеет два силовых участка (рис. 11.5, б), такие же участки необходимо выделить и в грузовом состоянии балки. Запишем по участкам выражения изгибающих моментов:
I. 0≤z1 ≤l2,
M |
|
=R |
z −q |
z2 |
|
ql |
z |
−q |
z2 |
; |
|
1 |
= |
|
1 |
||||||
|
2 |
|
2 |
|||||||
|
x |
A |
1 |
|
2 1 |
|
|
M11 =RA′ z1 = z21 .
156
а
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
в |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 11.5. |
|
|
|
|||||||||||||||
II. 0≤z2 ≤l 2, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
M |
|
|
=R z |
|
|
−q |
|
z |
2 |
|
|
|
ql |
z |
|
|
|
−q |
z2 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
2 |
= |
|
|
|
|
|
|
2 |
; |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
2 |
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
x |
|
|
|
B |
2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
M |
11 |
=R′ z |
2 |
= |
z2 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
B |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Записываем интеграл Мора и вычисляем прогиб: |
|
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
l |
2 |
M11M x |
|
|
|
|
|
l 2 |
M11M x |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
vc =∆1 = ∫ |
dz1 + ∫ |
dz2 = |
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
EIx |
|
|
|
|||||||||||||||||||||||||||||||||||||
= |
|
|
1 |
|
0 |
|
|
EIx |
|
|
|
|
0 |
|
z2 ql z2 −q z22 dz2 = |
||||||||||||||||||||||||||||||
|
|
|
∫ z1 ql z1 −q z12 dz1 |
+ ∫ |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
l 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l 2 |
|
|
|
|
|
|
|
|||||
|
|
|
EIx |
|
0 |
2 2 |
|
|
|
|
|
|
|
2 |
|
0 |
|
|
2 2 |
2 |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
= |
|
|
2 |
|
|
1 |
|
qlz3 − |
|
|
1 |
qz4 |
|
|
l0 |
2 = |
5 |
|
|
ql4 |
. |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
EIx 12 |
|
|
|
|
|
|
|
|
|
16 |
|
|
|
|
|
|
|
|
|
|
384 EIx |
|
|
|
Здесь учтено, что в силу симметрии подынтегральные выражения оказываются одинаковыми. Найденный прогиб совпадает по направлению с приложенной единичной силой.
157
Для нахождения угла поворота вновь рассматриваем вспомогательное состояние конструкции, возникающее под действием сосре-
доточенного безразмерного момента М = 1, приложенного в сечении В (рис. 11.5, в). Поскольку в этом случае, как и в грузовом состоянии, имеется только один силовой участок, интегрировать будем по переменной z1 по всей длине балки.
|
0≤z1 ≤l, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
M |
|
|
|
|
|
ql |
z |
−q |
|
z2 |
; |
|
|
|
|
|
||||||||
|
x |
= |
|
|
|
1 |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
2 |
|
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
M |
12 |
=R′′ z |
= |
z1 |
. |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
A |
1 |
|
|
|
l |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
l |
M12 M x |
|
|
|
|
|
|
1 |
|
|
l |
|
|
|
|
2 |
|
|
|||||||
θB =∆2 =∫ |
dz1 |
= |
|
|
|
∫ |
z1 |
|
|
ql |
z1 |
−q |
z1 |
dz1 |
= |
||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||
0 |
EIx |
|
|
|
|
|
|
EIx 0 l |
|
|
2 |
|
2 |
|
|
= |
1 |
|
1 |
qz3 |
− |
1 |
|
z14 |
|
|
l |
= |
ql3 |
. |
|
||||||||||||||
|
|
|
|
|
|
0 |
|
|||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
||
|
EIx |
6 |
|
|
8 |
|
l |
|
|
|
|
24EIx |
||
|
|
|
|
|
|
Результат получился положительным, следовательно, поворот сечения В происходит в направлении приложенного единичного момента, т.е. против часовой стрелки.
Определение перемещений способом Верещагина
Вычисление интеграла Мора можно упростить, если конструкция состоит из прямолинейных участков постоянной жёсткости. Упрощение основано на том, что внутренние усилия, возникающие от сосредоточенных единичных обобщённых сил, являются линейными функциями продольной координаты z.
Пусть на интервале от а до b требуется взять интеграл от произведения двух функций f1(z) и f (z), причём одна из них, например f1(z) – линейна:
f1 (z)=c+kz. |
(11.18) |
158
Тогда
J =∫b f1 (z) f (z)dz =∫b (c +kz)f (z)dz =c∫b f (z)dz +k ∫b zf (z)dz.
a |
a |
a |
a |
Первый из интегралов может быть представлен графически площадью, ограниченной кривой f(z), или, что то же самое, площадью эпюры этой функции ω (рис. 11.6):
∫b f (z)dz =∫dω=ω.
аω
Рис. 11.6.
Второй из интегралов, входящих в выражение для J, представляет собой статический момент площади ω относительно оси ординат у:
∫b zf (z)dz =∫zdω=ω zцт,
аω
где zцт – абсцисса центра тяжести эпюры f. Тогда
J =c ω+kω zцт =ω(с+kzцт ).
159
Но поскольку
с+kzцт = f1 (zцт ),
окончательно получим
J =ω f1 (zцт ).
Таким образом, интегрирование произведения двух функций, одна из которых линейна, можно заменить произведением площади эпюры функции ƒ на ординату линейной функции ƒ1, взятую под центром тяжести эпюры ƒ. Если линейны обе функции, то безразлично, перемножаются ли площадь первой эпюры с ординатой второй, или наоборот.
Способ Верещагина состоит в применении изложенного подхода к нахождению перемещений с помощью интеграла Мора. При плоском изгибе, например, выражение для перемещений (11.17) примет вид
n |
ωj M cj |
|
|
|
∆i =∑ |
. |
(11.19) |
||
|
||||
j=1 |
(EIx )j |
|
Здесь ωj – площади простейших фигур, на которые разбивается эпюра изгибающих моментов Mx;
M cj – ординаты с эпюры моментов М1i от единичной нагрузки,
взятые под центрами тяжести площадей ωj;
n – число простейших фигур, на которые разбивается эпюра грузовых моментов Мх.
Необходимо учесть, что перед разбиением на простейшие фигуры должны быть выделены участки постоянной жёсткости, на которых эпюра Мх не имеет разрывов, а эпюра моментов от единичной нагрузки линейна, т.е. не имеет разрывов и изломов.
В подавляющем большинстве случаев эпюры внутренних усилий можно разбить на простейшие фигуры трёх типов – прямоугольники, треугольники и сегменты квадратной параболы. Формулы для вычисления их площадей и положений центров тяжести приведено в таблице.
160