Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Основы механики твердого деформируемого тела

.pdf
Скачиваний:
128
Добавлен:
14.05.2015
Размер:
26.55 Mб
Скачать

ГЛАВА 5. МНОГОСТЕРЖНЕВЫЕ СИСТЕМЫ

5.1. Введение. Итак, анализ работоспособности сжатых стержней сводится к проверкам их прочности и устойчивости по формулам (4.10), а для того, чтобы вычисления по второй из этих формул были возможны, необходимо знать коэффициент µj приведенной длины рассматриваемого элемента. И в самом деле, множитель ϕ в соотношении (4.10b) зависит от гибкости λj стержня, а последнюю, согласно равенству (4.1), можно установить только тогда, когда известна величина µj . Эта величина определяется формулой (3.12), в которую входит сила Niкр(Pкр), так что становится совершенно очевидной необходимость отыскания критической нагрузки на рассматриваемую конструкцию в целом.

Если в зависимостях (3.1) и (3.5) модуль упругости E заменить на касательный модуль Et или приведенный модуль Er, то критическая сила для каждого из стержней, представленных в таблице 3.1, уменьшится, но коэффициенты их приведенных длин останутся прежними. Это означает, что величины (3.12) можно найти по критической силе, определенной на упругой стадии деформирования конструкции, а затем применить полученные результаты и при пластической работе материала.

5.2. Использование метода перемещений. При определении критической нагрузки на линейно деформируемую многостержневую конструкцию метод перемещений применяется наиболее часто. Пусть конструкция загружена так, что в исходном состоянии ее элементы испытывают только осевую деформацию (рис. 5.1a). Данная модель системы используется при расчете рам, воспринимающих вертикальную нагрузку. Последняя сводится в узлы точно так же, как и любое силовое воздействие на ферму. При смещении наложенных связей (рис. 5.1b, c) силы kiP порождают дополнительные изгибающие моменты в тех стержнях, которые искривляются при воздействиях Zj = 1. Например, к изгибающему моменту M 1, который возникает в стойке AB только от воздействия Z1 = 1, добавится усилие k1P u(x), так что в конечном счете в этой стойке появится дополнительный изгибающий момент k1P u(x)Z1. Таким образом, мера влияния узловых сил на напряженнодеформированное состояние конструкции определяется величиной перемещений ее узлов, а потому эти силы должны быть отнесены к единичным (а не к грузовому) состояниям основной системы. Но тогда M0 0 и все свобод-

692

Часть VI

ные члены Ri0 канонических уравнений метода перемещений обращаются в

нуль:

n

rikZk = 0, i = 1, 2, . . . , n.

(5.1)

k=1

Единичные реакции rik находят статическим способом по эпюрам ”Mk”, которые строят в основной системе при помощи стандартных табличных решений (см. таблицу 5.1 на с. 693), получаемых при помощи формул (3.5) или (3.5a). Пусть, например, верхний защемленный торец сжатого стержня, изображенного на рис. 5.2, поворачивается на угол, равный единице. Так как u0 = 0, u0 = 1, то первые две из формул (3.5) принимают вид:

 

1

 

 

M0

 

 

V0

 

 

 

u =

 

 

sin px +

 

 

(cos px − 1) +

 

 

(sin px − px),

 

p

p2EI

p3EI

 

u = cos px

 

M0

sin px +

V0

(cos px

 

1).

 

pEI

 

 

 

 

 

 

 

p2EI

 

 

С учетом условий закрепления второго конца стержня отсюда следует:

M0

(cos pl−1)+

V0

 

(sin pl−pl) = sin pl,

pEI

p2EI

 

 

M0

sin pl−

 

V0

(cos pl−1) = cos pl.

 

pEI

p2EI

Из этих двух неоднородных уравнений находят константы M0 и V0, после чего для определения усилий в любом сечении стержня останется воспользоваться двумя последними равенствами из соотношений (3.5). Результаты такого рода решений для стержней двух типов приведены в таблице 5.1. Совершенно аналогично, но уже при помощи формул (3.5a), получают усилия в растянуто-изогнутых стержнях.

Глава 5

693

Таким образом, коэффициенты системы (5.1) зависят от параметра P внешнего воздействия, а потому критическое значение этого параметра можно найти как наименьший корень трансцендентного уравнения

det (rik) = 0,

(5.2)

которое представляет собой условие существования нетривиального решения однородной системы линейных алгебраических уравнений, каковыми и являются уравнения (5.1).

Таблица 5.1

 

Сжатие

 

 

 

γ =

 

 

0, 5t3

 

 

 

2tg(t/2) − t

 

β =

 

γ

t − sin t

 

 

 

 

 

 

t2

sin t

 

 

 

2γ +t2

 

 

=

 

 

t3

α+β =

 

 

 

 

,

γ

 

 

 

 

4

 

tgt−t

 

 

 

 

 

 

α = γ + t2

Растяжение

0, 5t3

γ = t − 2th(t/2)

γsht − t

β= t2 sht

α+β =

2γ −t2

,

 

=

t3

γ

4

t − tht

 

 

 

 

α = γ + t2

l = J ,

t = l

EI

EI

 

 

P

Завершается расчет вычислением коэффициентов приведенной длины сжатых элементов по формуле (3.12). Для стержня с номером j, сжатого

µ2j lj2

694

Часть VI

силой Pj = kj P , будет (см. табл. 5.1)

t2j = kj P lj2/EIj .

При P = Pкр эта формула дает

t2j EIj kj Pкр = lj2 .

С другой стороны, для отдельного стержня (формула (3.10))

kj Pкр = π2EIj .

Из сравнения двух последних зависимостей следует, что

µj = π/tj .

Ничего, кроме примеров, больше не требуется.

На рис. 5.3a изображена стойка, критическая нагрузка на которую была найдена в п. 3.4. Рядом, т. е. на рис. 5.3b, приведены основная система метода перемещений и эпюра усилия M1. При ее построении использовалась таблица 5.1, согласно которой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2 tg t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2 tg t

 

 

 

 

1

=

 

1 + t12

=

 

1

 

 

 

1

 

,

 

 

2 =

 

 

2 + t22 =

 

2

2

.

 

α

γ

α

γ

 

tg t1 − t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg t2 − t2

Нижний индекс отвечает номеру участка стойки. Таким образом,

 

 

r11 =

EI

 

 

+

EI

 

 

=

EI

 

 

2t12 tg t1

+

 

t22 tg t2

.

 

 

 

 

 

 

 

 

 

α1

 

 

α2

 

 

 

 

 

 

 

 

 

 

l

2l

 

2l

tg t1 − t1

tg t2 − t2

Пусть t ≡ pl = l

 

. Тогда t1 = t, t2 = 2t и условие (5.2), сводящееся к

P/EI

уравнению r

11

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0, примет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg t

 

 

 

2tg 2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

= 0

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

tg t − t

tg 2t − 2t

 

 

 

 

 

 

 

 

 

 

 

 

 

3tg t · tg 2t − 2t(tg t + tg 2t) = 0.

 

 

Поскольку

 

 

 

 

 

 

 

(5.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin 3t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg t + tg 2t =

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos t · cos 2t

 

 

 

 

 

Глава 5

695

то равенство (5.3) сводится к уравнению (3.16a), а потому решение задачи будет таким же, что было получено еще в п. 3.4.

Следующий пример. На рис. 5.4a изображена портальная рама, находящаяся под воздействием равномерно распределенной нагрузки. Последняя сводится к узловым силам так, как это показано на рис. 5.4b. Портальная рама с вертикальными узловыми силами теряет устойчивость по обратно симметричной форме. И в самом деле, критическая нагрузка на раму будет тем больше, чем выше погонная жесткость ригеля. При kEI/l →∞ эта сила максимальна и в случае обратносимметричной формы потери устойчивости она будет такой же, как и для стойки, имеющей с одной стороны неподвижное защемление и подвижную заделку с другой. Критическая сила для такой стойки известна (см. рис. 5.5a). При нулевой изгибной жесткости ригеля критическая сила рамы минимальна. Если форма потери устойчивости симметрична, то эта сила совпадает с той, что имеет стержень, показанный на рис. 5.5b. Но тогда

min P

сим

=

π2EI

> max Pкробр.сим =

π2EI

,

кр

(0, 7h)2

 

 

 

 

h2

 

 

 

 

что и доказывает сделанное выше утверждение об обратной симметрии формы потери устойчивости рассматриваемой рамы.

Из сказанного следует, что основную систему необходимо назначить так, чтобы она отвечала обратносимметричной форме потери устойчивости (рис. 5.4c). Эпюры ”M1” и ”M2” строятся при

помощи таблицы 5.1. По ним обычным для метода перемещений образом определяются реакции

(η = kh/l)

r11 = 2(α1 + 3η), r12 = 2(α + β)1/h, r22 = 2γ1/h2,

после чего записывается условие (5.2), которое элементарными преобразованиями приводится к виду:

(α1 + 3η)γ1 (α + β)12 = 0.

(5.4)

696

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Часть VI

В формуле (5.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0, 5t3

 

 

 

 

 

 

2γ +t2

 

γ1

=

 

 

 

 

 

, (α+β)1

=

 

1

 

 

 

,

 

 

 

2tg(t/2)

−t

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

t2 sin t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

α

= (α + β)

 

 

γ1

(t−sin t)

;

t = h

 

 

P

.

(5.5)

 

 

 

 

 

 

Уравнение (5.3) решается численно в интервале 0 < t < π. В таблице 5.2 приведены результаты решения задачи при k = 1 и трех различных отношениях размеров h и l.

Достоинством описанного варианта применения метода перемещений является его простота. Но есть у него и недостатки. Во-первых, этот вариант требует модификации в случае, когда внешнее воздействие и после приведения его к узловым силам вызывает в исходном состоянии системы изгиб (рис. 5.6). Во-вторых, при сведении нагрузки в узлы меняется характер деформирования системы. Например, при нагрузке, показанной на рис. 5.4a, продольные силы возникают во всех трех стержнях рамы, а при узловом воздействии (рис. 5.4b) в ригеле осевой силы не будет. А поскольку заранее неизвестно, как отмеченное обстоятельство отразится на значении критической силы, то лучше выполнять расчеты на устойчивость по схеме, излагаемой ниже (см. также заключительный пункт предыдущей главы). Данная схема хороша еще тем, что позволяет решать вопрос об устойчивости конструкций, исходное безмоментное состояние которых невозможно даже после приведения заданной нагрузки к узловой (рис. 5.6).

 

 

 

 

 

 

 

 

 

 

Таблица 5.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходные

 

k = 1,

h = const

 

 

 

 

 

 

 

 

 

 

 

 

 

данные

l = h (η = 1)

 

l = 4h (η = 1/4)

l = 8h (η = 1/8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значение

 

 

 

 

 

 

 

 

 

 

 

 

2,7165

 

2,1746

 

1,9397

 

 

 

Pкр

 

 

tкр = h

 

 

 

 

 

 

 

 

 

EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчет начинается с определения продольных сил во всех стержнях конструкции от заданного (непреобразованного) воздействия. Пусть это силы

Глава 5

697

N (1), N (2), . . . , N (C) (C – число стержней). Очевидно,

N (j) = N (j)P,

где P – параметр нагрузки, N (j) – значение продольной силы в стержне j от

силы P = 1. После того, как числа N (j) будут найдены, вычисления ведутся обычным образом.

Найти продольные силы в раме, изображенной на рис. 5.4a, проще всего методом перемещений (одно неизвестное, если учитывать симметрию конструкции). Результаты расчета представлены на рис. 5.7a. Расчетная схема (см. рис. 5.7b) для решения задачи устойчивости определяется эпюрой продольных сил, при этом

 

 

P = ql/2,

η = kh/2l,

 

(1)

 

 

(2)

 

l

N

= 1, N

≡ k2 =

 

 

 

.

 

 

2h(2 + η)

При обратносимметричной форме потери устойчивости основная система и эпюры усилий M1 и M2 остаются такими же, что и в предыдущем примере (рис. 5.4c), но только теперь значение изгибающего момента на левом конце ригеля при смещении Z1 = 1 будет равно числу 2ηα2, а не 6η. Это приводит к иному значению реакции r11:

 

 

 

 

 

 

 

 

t2 tg t

 

 

 

 

 

 

 

 

 

 

2

 

r11

= 2(α1 + ηα2), α2

=

 

,

tg t2

− t2

 

 

 

 

 

 

 

 

 

и к новому виду характеристического уравнения задачи:

 

(α1 + ηα

2)γ1 (α + β)12 = 0.

(5.4a)

Величины γ1, α1 и (α + β)1 по-прежнему определяются формулами (5.5), но,

кроме того,

 

 

 

 

 

 

 

 

 

 

 

 

t2

2

 

 

 

= ρt;

ρ = 2h

 

 

 

.

kEI

 

k

 

 

l

 

 

k2P

 

l

 

 

k2

698

Часть VI

Однако самое существенное отличие между системами, помещенными на рис. 5.4b и 5.7b, заключается в том, что первая из них всегда теряет устойчивость по обратносимметричной форме, а вторая в зависимости от отношения размеров h и l имеет обратносимметричную либо симметричную форму потери устойчивости. Стало быть, при анализе устойчивости портальной рамы со всеми тремя сжатыми элементами ограничиваться изучением лишь

ееобратносимметричного деформированного состояния нельзя.

Всимметричной задаче (nпер = 1) условие r11 = 0 существования невырожденного решения (рис. 5.8) принимает вид:

 

α1 + η(α − β)2 = 0,

 

 

(5.6)

где

 

 

 

 

 

 

(α − β)2 =

t2

1

− t2

1 + cos t

2

.

2tg (t2/2) − t2

2 sin t2

 

Уравнения (5.4a) и (5.6) решаются численно. Результаты решения при k = 1 и трех различных отношениях h/l приведены в таблице 5.3 (см. также табл. 5.2). Критические

значения tкр = h Pкр/EI параметра t в таблице 5.3 набраны жирным шрифтом, а значения t, отвечающие неустойчивым (нереализуемым) отклоненным формам равновесия системы, заключены в скобки. В третьей строке снизу указана разница ∆ между критическими силами, найденными с учетом и без учета (т. е. при k2 = 0) сжатия ригеля. Из таблицы видно, что только от исходных данных зависит, будут ли вычисления при k2 = 0 и k2 = 0 приводить к практически одинаковым результатам (случай η = 1), к весьма далеким результатам (в случае η = 1/4) или даже к результатам, не имеющим между собой ничего общего, как при η = 1/8, когда различны даже формы потери устойчивости.

Следует также обратить внимание на характер изменения коэффициентов приведенной длины элементов рамы (2 последние строки табл. 5.3), вычисляемых по формуле (3.12). С ростом отношения l/h величина µ2 уменьшается, что говорит о повышении жесткости закрепления ригеля. Так оно и есть на самом деле, ибо ригель представляет собой упругозащемленный по торцам стержень, причем роль упругих заделок берут на себя стойки. Жесткость этих заделок тем больше, чем больше отношение длины l ригеля к высоте h стоек. Однако гибкость λ = µ2l/r ригеля для всех трех помещенных в таблице 5.3 конструкций одинакова, так как уменьшение числа µ с ростом величины l/h компенсируется увеличением размера l. Понижением погонной жесткости ригеля с возрастанием величины l/h, а значит, и

Глава 5

699

уменьшением его (ригеля) влияния на деформируемость стоек, объясняется увеличение с ростом l/h коэффициента µ1.

 

 

 

 

 

Таблица 5.3

 

 

 

 

 

 

 

 

 

 

 

 

 

Схема

 

 

 

 

 

рамы

 

 

 

 

 

EI = const

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

ql/2

2ql

4ql

 

 

tкр (k2

= 0)

2,7165

2,1746

1,9397

 

 

Табл. 5.2

 

 

 

 

 

 

 

 

 

 

 

 

 

tкр (k2 = 0)

2,7096

1,6313

(1,1940)

 

 

Обр. сим.

 

 

 

 

 

tкр (k2

= 0)

(4,9828)

(2,8534)

1,0758

 

 

Симметр.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,51%

77,7%

225%

 

 

 

 

 

 

 

 

 

µ1

 

1,159

1,925

2,920

 

 

 

 

 

 

 

 

 

µ2

 

4,016

1,021

0,532

 

 

 

 

 

 

 

 

5.3. Использование метода сил. Этот метод при решении задач устойчивости применяют редко, но познакомиться с методом сил имеет смысл по двум причинам. Во-первых, надо понять, чем объясняется сравнительно малая его популярность. Во-вторых, на его основе проще обсудить некоторые особенности, которые могут встретиться при исследовании критических состояний многостержневых конструкций.

Как и при расчете методом перемещений, предварительно находят продольные силы во всех стержнях, после чего пролетную нагрузку отбрасывают (рис. 5.9a). Затем выбирается основная система метода сил. Делается это так, чтобы основная система в исходном состоянии под действием только осевых сил kj P не изгибалась (рис. 5.9b). Если же грузовое состояние основной системы безмоментно, то свободные члены канонических уравнений

700

Часть VI

равны нулю: Ui0 = 0, а потому эти уравнения становятся однородными:

n

uikXk = 0, i = 1, 2, . . . , n.

(5.7)

k=1

Что же касается сил k1P, . . . , kCP , то они дадут о себе знать при вычислении перемещений uik, ибо воздействия Xk = 1 приведут к отклонению основной системы от исходного положения равновесия и, как следствие, – к дополнительному изгибу ее стержней силами kj P . Величины uik находят при помощи интегралов (M i, Mk) либо (Mi, M k), где M i – изгибающие моменты, порождаемые в основной системе силой Xi = 1, приложенной по направлению искомого перемещения, а Mk – изгибающие моменты, обусловленные одновременным приложением силы Xk = 1 и нагрузок kj P (т. е. воздействием, приводящим к перемещению uik). Если эпюры M i получают обычным для метода сил образом, то усилия Mk находят при помощи формул (3.5).

Дальнейшее ясно. Поскольку uik = uik(P ), надо, уходя от тривиального решения системы линейных алгебраических уравнений (5.7), приравнять к нулю ее определитель:

det(uij ) = 0,

(5.8)

а затем найти наименьший корень P данного трансцендентного уравнения.

В качестве примера можно повторно решить задачу о критической нагрузке на раму, приведенную на рис. 5.4b. Поскольку она теряет устойчи-