Основы механики твердого деформируемого тела
.pdf
Глава 5 |
461 |
силовым воздействием. Перемещение же от силового воздействия вычисляется по формуле Мора:
UAX = (M, |
|
Aсил). |
(5.18) |
M |
Здесь M – изгибающий момент от сил X1, . . . , Xn. Он определяется равенством (5.16). Изгибающий момент M Aсил отыскивается в основной системе от воздействия PA = 1, прикладываемого по направлению искомого перемещения.
Второе слагаемое в формуле (5.17) представляет собой перемещение основной системы, возникающее от заданных осадок Uj ее опор (рис. 5.16d). Его находят из геометрических рассмотрений либо по формуле, аналогичной равенству (5.13):
|
|
m |
UAсил = |
|
j |
|
− |
jA j |
|
|
=1 |
Реакция rjA связи j основной системы вызывается воздействием PA = 1. Таким образом (см. формулы (5.17)–(5.19)),
UA = (M, |
|
|
Aсил) + UAсил |
|
(5.20) |
|||
M |
|
|||||||
или |
|
m |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
Aсил) |
|
j |
|
|
|
UA = (M, |
M |
− |
r U . |
(5.20a) |
||||
|
|
|
|
|
jA |
j |
|
|
|
|
|
|
|
|
=1 |
|
|
На равенстве (5.20a) основана кинематическая проверка эпюры "M ". Поскольку перемещение по направлению устраняемой связи должно равняться
i |
U = (M, M |
сил) |
m |
r |
U . |
||
заданному значению U |
, то, согласно формуле (5.20a), |
||||||
|
|
|
|
i − |
j |
|
|
|
i |
ji |
j |
||||
|
|
=1 |
|||||
|
|
|
|
|
|
|
|
Второе слагаемое в правой части последнего равенства равно перемещению Ui0 (см. зависимость (5.13)), а усилие M силi есть не что иное, как единичный изгибающий момент M i. Стало быть,
i = ( |
|
|
i) + |
|
i0 |
, i = 1, 2, . . . , n. |
U |
M, M |
|
U |
|
||
Все эти n равенств можно сложить:
n |
n |
n |
||
|
i |
|
|
|
Ui = |
(M, M i) + |
Ui0. |
||
i=1 |
=1 |
|
|
i=1 |
Интеграл от суммы равен сумме интегралов от отдельных слагаемых, поэтому
n
(M, M i) = (M, M (S))
i=1
462 |
|
|
|
|
Часть IV |
и окончательно |
n |
|
|
||
|
|
|
|
|
|
|
|
(S)) = |
i |
|
|
(M, |
M |
|
(Ui − Ui0). |
(5.21) |
|
|
|
|
=1 |
|
|
Если при вычислении произведения (M, M (S)) используются коэффициенты податливости αj = EI0/EIj, то равенство (5.21) заменяется соотноше-
нием |
|
|
||
|
|
n |
|
|
|
|
i |
|
|
(M, |
M |
(S)) = EI0 |
(Ui − Ui0). |
(5.21a) |
=1 |
|
|
||
Завершает рассказ о расчете статически неопределимых конструкций на кинематическое воздействие следующий пример. На рис. 5.17a изображена балка, правая опора которой испытывает поворот на угол ϕ. Требуется найти вызываемые этим воздействием усилия и прогиб в середине пролета. Балка имеет только одну лишнюю связь, а потому и условие совместности перемещений только одно:
u |
11X1 |
+ U10 |
U |
. |
(5.22) |
|
= 1 |
|
|
В основной системе, показанной на рис. 5.17b, связей, испытывающих смещения, нет, так что U10 = 0. По рис. 5.17a, b видно, что
|
= ( |
|
|
|
|
|
) = |
l |
, |
|
= |
|
|
u11 |
M |
, M |
U |
ϕ. |
|||||||||
|
3EI |
||||||||||||
1 |
1 |
|
|
1 |
|
||||||||
В этом случае из уравнения (5.22) следует
U |
/u |
|
= |
3EIϕ |
. |
|
|
||||
X1 = 1 |
|
11 |
l |
||
Эпю ра изгибаю щих моментов M = M1X1 приведена на рис. 5.17c. Эпюру поперечных сил изображать нет смысла, ибо она имеет вид прямоугольника с
ординатой 3EIϕ/l2.
Глава 5 |
463 |
Для определения перемещения UA надо приложить в основной системе вертикальную силу PA = 1 и построить эпюру "M силA ". Результаты указанных действий представлены на рис. 5.17d. В основной системе нет связи, получившей смещение, поэтому UAсил = 0. Стало быть (см. формулу (5.20)),
|
|
сил |
|
1 |
|
1 |
l |
|
1 |
|
3EI |
|
3 |
|
|||
|
|
|
|
|
|
|
|||||||||||
UA = (M, M A |
) = |
|
· |
|
|
l · |
|
· |
|
· |
|
|
ϕ = |
|
lϕ. |
||
EI |
2 |
4 |
2 |
|
l |
16 |
|||||||||||
Интересно оценить величину угла поворота ϕ, при котором наибольшие нормальные напряжения в балке становятся равными допускаемым. Пусть поперечное сечение бруса имеет горизонтальную ось симметрии и потому
W = 2I/h, где h – высота сечения. Тогда |
|
|
|
||||
max σx = |
max M |
= |
3EIϕ/l |
= |
3h |
||
|
|
|
|
Eϕ. |
|||
W |
2I/h |
2l |
|||||
При max σx = [σ] отсюда следует
ϕ= 2 [σ] l .
3 E h
Для металлических балок [σ]/E ≈10−3, l/h ≈ 10, так что ϕ ≈ 10−2 радиан. Другими словами, предельно допустим поворот фундамента балки, примерно, на 40 минут. Это говорит о том, что кинематическое воздействие на конструкцию может быть опасным. Имеется немало примеров того, как здания и сооружения получали большие повреждения и даже разрушались в результате неравномерных осадок фундаментов. Усилия в статически неопределимых конструкциях, возникающие при кинематическом воздействии, пропорциональны изгибным жесткостям их элементов, что видно и по рис. 5.17c. Вот почему обычно жесткость сооружений стремятся понизить, устраивая, например, противоосадочные швы.
5.6. Расчет на тепловое воздействие. Задача состоит в отыскании усилий и перемещений в статически неопределимой конструкции при изменении температуры окружающей среды. Канонические уравнения метода сил этой задачи записываются в форме (5.2), но только теперь свободные члены имеют смысл перемещений основной системы, вызываемых температурным, а не силовым, воздействием. Если температура меняется по поперечным сечениям стержней линейно и находится в пределах, позволяющих считать физические свойства материала неизменными, то перемещения Ui0 могут быть найдены по формулам, приведенным в п. II.9.8. В частности, при призматических стержнях с поперечными сечениями, высоты hj которых делятся центральными осями ozj пополам,
C |
|
|
|
|
|
|||
|
|
|
∆Tjωij |
|
|
|||
Ui0 = j=1 αj N ijTср.jlj + |
. |
(5.23) |
||||||
hj |
||||||||
464 |
Часть IV |
Через αj и lj обозначены коэффициент теплового расширения материала, из которого изготовлен стержень с номером j, и его длина. Средняя температура крайних волокон j-го стержня и разница между их температурами обозначаются через Tср.j и ∆Tj соответственно. Наконец, N ij и ωij – единичная продольная сила и площадь эпюры изгибающих моментов M i на стержне j. Из сказанного следует, что при расчете конструкции на изменение температуры в основной системе обязательно должны быть найдены единичные продольные силы, т. е. кроме эпюр "M i" нужно строить и эпюры "N i". Поскольку в формулу (5.23) не входят изгибные жесткости стержней, то в случае, когда при вычислении единичных перемещений uij используются коэффициенты податливости, канонические уравнения записываются следующим образом:
u11X1 + u12X2 + · · · + u1nXn + EI0U10 = 0, |
|
|
|||
u21X1 + u22X2 + |
· · · |
+ u2nXn + EI0U20 = 0, |
|
(5.24) |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
un1X1 + un2X·2·+· |
|
+· · ·unnXn· ·+· |
|
|
|
|
EI0Un0 = 0. |
|
|
||
|
· · · |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для построения эпюры искомых изгибающих моментов используется та же формула (5.16), что и при расчете на кинематическое воздействие: ведь и изменение температуры не вызывает усилий в статически определимой основной системе, т. е. M0 = 0. Да и формула для перемещений конструкции внешне не отличается от равенства (5.20), предназначенного для вычисления смещений фиксированных сечений сооружения при осадке его опор:
UA = (M, |
|
Aсил) + UAсил. |
(5.25) |
M |
И в самом деле, к этому результату можно прийти при помощи тех же рассуждений, что привели в предыдущем пункте к равенству (5.20). Но теперь перемещение UAсил порождается изменением температуры, а не осадкой опор, и для его определения должна использоваться формула
C |
|
|
|
|
|
|
|
|
|
сил |
Tср.jlj + |
∆TjωAjсил |
. |
(5.26) |
|||
|
|
|||||||
|
|
|||||||
UAсил = j=1 αj N Aj |
hj |
|||||||
|
|
|
|
|
|
|
|
|
Через N силAj и ωсилAj обозначены продольная сила и площадь эпюры изгибающих моментов M силA , которые построены на стержне j от силы PA = 1. Это воздействие прикладывается по направлению искомого перемещения в любой основной системе метода сил.
Глава 5 |
465 |
Для кинематической проверки эпюры изгибающих моментов используется равенство
|
|
|
C |
|
(M, |
|
(S)) = − |
j |
(5.27) |
M |
Ui0, |
|||
|
|
|
=1 |
|
которое вытекает из формулы (5.25). Для обоснования сказанного можно использовать те же доводы, что приводились при выводе соотношения (5.21). Если интеграл (M, M (S)) вычисляется с привлечением коэффициентов податливостей стержней конструкции, то формулу (5.27) нужно изменить:
|
|
C |
|
|
|
j |
|
(M, |
M |
(S)) = −EI0 Ui0. |
(5.27a) |
=1 |
|
||
Последовательность вычислений при расчете на тепловое воздействие демонстрирует следующий пример. На рис. 5.18a изображена рама, которая в свое время была рассчитана на силовое воздействие (см. п. 5.4). Теперь конструкция подвергается воздействию, состоящему в указанном на рис. 5.18a изменении температуры окружающей эту конструкцию среды. Коэффициенты канонических уравнений были найдены в п. 5.4.3, так что для формирования системы (5.24) достаточно вычислить лишь свободные члены Ui0. Для этого в той же самой основной системе, которая использовалась раньше, надо построить единичные эпюры продольных сил. Для удобства дальнейших вычислений на рис. 5.18a и e заново изображены основная система и единичные эпюры изгибающих моментов (см. также рис. 5.9–5.11). Кроме того, на рис. 5.18e приведены дополнительно построенные эпюры "N i"и "N (S)".
466 |
Часть IV |
Пусть конструкция удовлетворяет всем требованиям, позволяющим при вычислении перемещений Ui0 пользоваться формулой (5.23). Для более удобного обращения к этой формуле целесообразно представить температурное воздействие в виде суперпозиции средних температур Tср крайних волокон стержней и перепада ∆T температур названных волокон, как это показано на рис. 5.18c, d. Разность температур записывается со стороны тех волокон, где температура выше. Пусть материал всех стержней конструкции одинаков, т. е. αj = α = const, а высоты поперечных сечений элементов рамы в 20 раз меньше их длин.
Техника использования формулы (5.26) такова. Сначала записываются все слагаемые, связанные с осевой деформацией стержней (в нижеприводимых записях они заключены в первую пару квадратных скобок), а затем
– слагаемые, сопряженные с деформацией изгиба. Запись слагаемых осуществляется в порядке нумерации стержней. Значения средних температур и продольных сил берутся со своими знаками. Знаки слагаемых, связанных с изгибной деформацией, устанавливаются следующим образом. Если численное значение температуры на схеме воздействия ∆T и ординаты эпюры "M i" отложены с одной и той же стороны стержня с номером j, то соответствующий член в формуле (5.26) принимается положительным. Остальное ясно из приводимых ниже подсчетов.
U10 |
= α [0+1 · (−6) · 2l+0+0]+ l/20 |
2 l · l+ |
2l/20 |
· 2l · l+0+ l/20 |
2 l · l |
= 448αl, |
||||||||
|
4 |
1 |
|
|
|
20 |
|
|
2 |
|
1 |
|
||
U20 |
= α [−1 · 2 · l+0+0+1 · 3 · l]+ 0+ 2l/20 2 |
· 2l · 2l+0+ l/20 l · 2l |
= 481αl, |
|||||||||||
|
|
20 |
|
|
1 |
|
|
2 |
|
|
|
|||
U30 = α [0+0+0+1 · 3 · l]+ 0+0+ 2l/20 |
2 · 2l · l− l/20 · l · 2l = 283αl; |
|||||||||||||
|
|
18 |
|
|
1 |
|
|
|
2 |
|
|
|
|
|
3
Ui0 = (448+481+283)αl = 1212αl.
i=1
Для проверки вычислений нужны эпюры усилий M (S) и N (S), приведенные на рис. 5.18e. С их помощью находят суммарное перемещение
U (S) = i=1 αj[N j(S)Tср.jlj +(∆Tj/hj)ωj(S)] = α [−1 |
· 2 · l+1 · (−6) · 2l+0+2 · 3 · l]+ |
||||||||||||||||
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ l/20 2 |
· l · l+ 2l/20 |
2 (l+3l) · 2l+ |
2l/20 2 |
· |
2l · 2l+ l/20 |
2 l · l |
= 1212αl. |
||||||||||
4 |
1 |
20 |
|
1 |
|
|
|
18 |
1 |
|
2 |
|
1 |
|
|
||
Глава 5 |
467 |
Так как U (S) = U10 +U20 +U30, то свободные члены канонических уравнений вычислены правильно.
Далее коэффициенты uij, установленные в п. 5.4.3, и только что вычисленные перемещения Ui0 подставляются в систему (5.24). После деления обеих частей получившихся равенств на 4l3 и использования обозначения
B = αEI0/l2
разрешающие уравнения задачи примут вид:
|
|
|
6 X1 + 2 X2 − X3 + 448B = 0, |
|
||||||||||||
3 |
7 |
|
|
3 |
|
|
|
|
|
|
|
|||||
|
|
|
14 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X1 + |
|
|
|
|
X2 |
|
|
4X3 + 481B = 0, |
|
||||
2 |
|
|
|
3 |
|
− |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
X1 |
− |
4X2 + 14 X3 + 283B = 0. |
|
||||||||||||
|
|
|
|
|
|
|
3 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
И этой системе, и суммарному уравнению |
||||||||||||||||
5 |
X1 |
+ |
13 |
X2 |
− |
1 |
X3 + 1212B = 0 |
|||||||||
|
|
|
|
|
|
|
||||||||||
|
|
3 |
6 |
3 |
||||||||||||
удовлетворяет решение |
|
|
|
|
||||||||||||
X1 = −82956B/386, |
|
|
X2 = −182514B/386, |
|||||||||||||
X3 = −197625B/386.
Искомая эпюра изгибающих моментов и все ее составляющие, которые используются при построении этой эпюры по фор-
муле M = M 1X1+M 2X2+M 3X3, представлены на рис. 5.19. Кинематическая проверка полученной эпюры ради краткости не приводится. Она сводится к вычислениям по формуле (5.20a), т. е. к выполнению уже знакомых операций. Нет нужды тратить время и на построение эпюр "Q" и "N ".
Любопытен следующий результат найденного решения: растянутыми оказались менее нагретые волокна стержней рассчитанной рамы (см. эпюру "M " на рис. 5.19 и схему температурного воздействия ∆T на рис. 5.18c). Данный эффект наблюдается при расчете на изменение температуры всех статически неопределимых стержневых конструкций. Исключения возможны лишь на небольших участках осевого контура (в рассматриваемом примере – на нижней части стержня с номером 4). Все дело в том, что элементы статически неопределимой стержневой конструкции деформируются не только от изменения температуры окружающей среды, но и от вызванных таким изменением реакций лишних связей. Влияние последних на харак-
468 |
Часть IV |
тер деформирования оказывается более значительным. Нечто подобное наблюдается при повышении температуры жидкости или газа, заключенных в сферической или цилиндрической емкости: более растянутыми оказываются менее нагретые внешние слои сосуда.
Ординаты эпюр усилий, возникающих в статически неопределимой конструкции при изменении температуры, как и в случае кинематического воздействия, прямо пропорциональны изгибной жесткости EI. Стало быть, чем выше жесткость конструкции, тем опаснее будет ее состояние при одном и том же тепловом воздействии. В рассматриваемой раме наибольшие напряжения изгиба возникают в правом торцевом сечении стержня с номером 2. При h = 2l/20 = l/10, W = 4I/0, 5h = 80I/l
max σ |
|
= |
74664A |
= |
74664 · 12αEI/193l |
≈ |
58αE. |
|
x |
W |
80I/l |
||||||
|
|
|
|
Для стали α = 1, 3 · 10−5 1/град, E = 2, 1 · 106 кГ/см2 и
max σx = 1583, 4 кГ/см2.
Эти напряжения почти равны допускаемым, а ведь при их вычислении не учитывалась продольная сила. Это означает, что прочность рассматриваемой конструкции будет исчерпана при одном лишь температурном воздействии. То, что изменение температуры может привести к большим напряжениям и деформациям, инженерам известно. Чтобы уменьшить неблагоприятные последствия температурного воздействия на конструкцию, прибегают к теплоизоляции ее элементов, устройству температурных компенсаторов на трубопроводах и температурных швов в зданиях и сооружениях.
5.7. О рациональной основной системе. Как говорилось в п. 5.4.1, существует множество основных систем, пригодных для расчета конкретной статически неопределимой конструкции методом сил, а значит, существует и проблема выбора лучшей основной системы. Данное понятие относительно. Для начинающего изучать методы прочностных расчетов лучше та основная система, в которой проще строятся эпюры усилий, т. е. консольная, а опытный специалист отдаст предпочтение основной системе, обеспечивающей высокую точность вычислений при наименьшем их объеме. Такие основные системы называют рациональными. О рациональных основных системах и пойдет речь в настоящем параграфе.
5.7.1. Локальные эпюры усилий. На рис. 5.20a изображена неразрезная балка, имеющая четыре лишние связи. Для расчета балки предлагается использовать либо основную систему A (рис. 5.20b), либо основную систему B (рис. 5.20c). Под каждой основной системой показаны единичные и грузо-
Глава 5 |
469 |
вая эпюры изгибающих моментов, причем единичные эпюры, компактности ради, построены на одной общей базе, т. е. совмещены.
Можно показать, что второй вариант во всех отношениях предпочтительнее. Так, при построении эпюр в основной системе A требуется вычислить 16 ординат на эпюрах усилий M 1–M 4 и еще 7 ординат на эпюре "M0" (см. рис. 5.20b), чему должно предшествовать определение опорных реакций. В основной системе B любые эпюры строятся гораздо проще, а вычислить (по известным формулам) придется всего две ординаты эпюры "M0". Тот, кто выберет о. с. A, затратит много труда на вычисление коэффициентов и свободных членов канонических уравнений. Ведь любой из интегралов (M i, M k) и (M i, M0) должен браться на всех пяти участках балки, т. е. в суммах (5.5a) не выпадет ни одно из слагаемых. Кроме того, ни один из 10 вычисляемых интегралов (учитывается, что uij = uji) не обращается в нуль. Если же выбрать о. с. B, то при вычислении величин uii и Ui0 в формулах (5.5a) ненулевыми останутся только два слагаемых, а при подсчете коэффициентов uik – одно слагаемое для случая k = i+1 и ни одного слагаемого для k > i+1. Единичные эпюры "M i" и "M k" при k > i+1 оказываются ортогональными, ибо они расположены на разных участках основной системы B.
Из сказанного выше ясно, что матрицы DA и DB коэффициентов разрешающих уравнений задачи, соответствующие сопоставляемым вариантам основной системы, имеют вид:
|
u21 |
u22 |
u23 |
u24 |
|
|
|
u21 |
u22 |
u23 |
0 |
|
|||||||
|
|
u11 |
u12 |
u13 |
u14 |
|
|
|
|
u11 |
u12 |
0 |
0 |
. |
|||||
DA = |
u31 |
u32 |
u33 |
u34 |
, |
DB = |
0 |
u32 |
u33 |
u34 |
|||||||||
|
u |
|
u |
|
u |
u |
|
|
|
|
|
0 |
0 |
u |
|
u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
41 |
|
42 |
43 |
|
44 |
|
|
|
|
|
|
|
43 |
|
44 |
|
470 |
Часть IV |
Матрица DB является трехленточной, т. е. разреженной, тогда как матрица DA заполнена целиком. Объем вычислений при решении системы алгебраических уравнений с ленточной матрицей существенно меньше, чем с заполненной. Не менее важно и следующее отличие между матрицами DA и DB. Судя по эпюрам "M i", изображенным на рис. 20b, интегралы (M i, M k) близки друг к другу по величине, а потому близки по величине диагональные uii и внедиагональные (побочные) uik элементы матрицы DA. Наоборот, диагональные элементы матрицы DB, по крайней мере, вдвое больше ее побочных элементов. О таких матрицах говорят, что они имеют диагональное преобладание. Из вычислительной математики известно, что решение системы линейных алгебраических уравнений с матрицей, которой присуще диагональное преобладание, может быть получено с высокой точностью, тогда как при отсутствии такового точность решения системы гарантировать нельзя. Качественное объяснение этому факту дается в п. 8.3.
Наконец, о построении эпюры изгибающих моментов по формуле
M = M0 + |
|
1X1 + · · · + |
|
nXn. |
(5.28) |
M |
M |
При основной системе A надо все ординаты эпюр "M i" в точках, где расположены опоры балки, умножить на полученные значения основных неизвестных, а затем сложить результаты друг с другом и с соответствующими ординатами эпюры "M0". При выборе о. с. B никаких арифметических действий, связанных с использованием формулы (5.28), выполнять не приходится (операция 1 · Xi, конечно же, не в счет).
Преимущества основной системы B на всех этапах расчета конструкции объясняются тем, что единичные и грузовые состояния основной системы локальны. Локальными называют усилия и перемещения, которые распространяются не на весь осевой контур конструкции. Как будет показано в главе 8, локальные единичные и грузовые эпюры метода сил могут быть получены сравнительно просто для любых стержневых конструкций.
5.7.2. Ортогональные единичные эпюры. Имеется еще один класс единичных эпюр, использование которых ведет к повышению точности расчета, а в ряде случаев – и к уменьшению объема вычислений. На рис. 5.21a изображена дважды статически неопределимая рама, а на рис. 5.21b, c приведены две ее основные системы. Эти основные системы отличаются друг от друга только способом разложения реакции шарнирной опоры на составляющие X1 и X2. Но если более привычный способ приводит к каноническим уравнениям с заполненной матрицей, то при выборе основной системы B матрица канонических уравнений будет диагональной:
DA = |
u11 |
u12 |
, DB = |
0 u |
|
. |
|
|
u |
u |
|
u11 |
0 |
|
|
|
21 |
22 |
|
|
|
22 |
|
