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

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

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

Глава 5

471

Это и понятно: эпюры усилий M 1 и M 2, показанные на рис. 5.21c, взаимно ортогональны, т. е. их произведение равно нулю:

u12 = (M 1, M 2) = 0.

Использование ортогональных эпюр приводит к тому, что система канонических уравнений распадается на независимые уравнения. В данном примере это уравнения

u11X1 + U10 = 0, u22X2 + U20 = 0,

из которых сразу же могут быть получены значения основных неизвестных X1 и X2. Объем вычислений минимален, тогда как их точность максимально возможная, ибо отсутствует погрешность, связанная с исключением неизвестных при решении системы алгебраических уравнений.

Следующий пример. Портальная рама, изображенная на рис. 5.22a, обладает симметрией (считается, что стойки рамы имеют одинаковую изгибную жесткость). Симметричными являются и основная система A, показанная на рис. 5.22b, и первые две ее единичные эпюры изгибающих моментов. Что же касается эпюры "M 3", то она относится к классу обратносимметричных эпюр. Поэтому u13 = u23 = 0, но u12 = 0. Такие эпюры называют частично ортогональными. При частично ортогональных единичных состояниях основной системы канонические уравнения полностью не распадаются. В данном примере отделяется только третье уравнение:

472

 

 

 

Часть IV

DA =

u21

u22

0

.

 

u11

u12

0

 

 

 

 

 

 

 

 

 

 

 

0

0

u33

 

Полученный результат не так уж и плох, но в данной задаче можно добиться и большего´. Для этого внутренние силы взаимодействия в месте разреза рамы надо привести к силам X1, X2 и X3, показанным на рис. 5.22c (о. с. B). В такой основной системе второе и третье единичные состояния будут такими же, что и в о. с. A, а усилие M 1 станет ортогональным и к

усилию M 2. Теперь

u12 = u13 = u23 = 0

и система канонических уравнений распадается на три самостоятельных уравнения.

Аналогичный прием используется и при расчете бесшарнирных арок (рис. 5.23a). Бесшарнирные арки часто применяются в большепролетных мостовых сооружениях. Чтобы уменьшить их собственный вес, поперечное сечение конструкции делают переменным, уменьшая его размеры от опор к замку. Мостовые арки – дорогостоящие конструкции, их проектируют весьма тщательно. В частности, при выполнении прочностного расчета учитывают влияние на перемещения не только изгибающего момента, но и поперечных и продольных сил. Если принять во внимание еще и нелинейность контура арки, то станет ясным, что расчет данной конструкции многотруден и по возможности его надо упрощать. Для этого как раз и используется ортогонализация всех единичных состояний основной системы.

Полная ортогонализация достигается назначением основной системы, показанной на рис. 5.23b, и специальным подбором расстояния a, на которое от центра тяжести замкового сечения переносится сила X2. По единичным эпюрам, приведенным на рис. 5.23b, видно, что

 

1 = 1,

 

1 =

 

 

1 = 0,

 

2

= f −a−y,

 

2

= sin ϕ,

 

2

= cos ϕ,

M

Q

N

M

Q

N

 

 

 

 

 

 

3 = l/2−x,

 

 

3 = cos ϕ,

 

 

3 = sin ϕ.

(5.29)

 

 

 

 

M

 

Q

 

N

 

Глава 5

 

 

 

 

 

 

473

Тогда

 

 

 

 

 

 

 

u11 = (M 1, M 1) =

l

 

EI

 

 

 

 

 

 

 

 

 

ds

 

 

 

 

0

 

 

 

или (ds = dx/ cos ϕ)

 

 

 

 

 

 

 

u11 =

l

 

 

 

 

 

 

dx

.

(5.30)

 

EI cos ϕ

0

 

 

 

 

 

 

 

Размер a выбирается таким, чтобы первое и второе единичные состояния арки стали ортогональными, ибо равенства u13 = 0 и u23 = 0 выполняются при любом значении величины a. Так как (см. формулы (5.29) и (5.30))

u

 

=

l

f −a−y

dx = (f

a)u

11

l

y dx

,

 

 

EI cos ϕ

 

EI cos ϕ

 

12

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

то условие u12 = 0 ортогональности состояний 1 и 2 основной системы дает

1

 

l

 

y dx

a = f −

 

 

 

.

u11

EI cos ϕ

 

 

0

 

 

 

При таком размере a система канонических уравнений решаемой задачи распадается на изолированные равенства, из которых следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U10

 

 

 

 

 

 

 

 

U20

 

 

 

 

 

 

 

 

 

U30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X1 =

 

 

 

 

, X2 =

 

 

 

 

, X3 =

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u11

u22

u33

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

=

 

l

 

 

dx

 

,

 

 

u

 

 

 

=

l

 

(f

 

 

a−y)2

+µ

sin 2ϕ

+

cos 2ϕ

 

 

dx

,

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

22

0

 

 

 

EI

 

 

 

 

 

 

 

 

GF

 

 

 

 

 

 

 

EF

 

cos ϕ

 

 

 

EI cos ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

=

 

l

 

 

(l/2−x)2

 

+µ

cos 2ϕ

+

sin

2ϕ

 

 

 

dx

,

U

 

 

=

 

l

M0dx

,

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

33

 

 

 

EI

 

 

 

 

 

 

 

 

GF

 

 

 

 

 

EF

 

cos ϕ

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI cos ϕ

 

 

 

 

 

 

 

 

 

U =

l

 

(f −a−y)M0

 

 

 

µ

Q0tg ϕ

 

N0

 

 

dx ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

0

 

 

 

EI cos ϕ

 

 

 

 

 

 

 

GF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

=

 

l

 

(l/2−x)M0

µ

Q0

+

N0tg ϕ

 

dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

GF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

EI cos ϕ

 

 

 

 

 

 

EF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

474 Часть IV

5.7.3. Использование симметрии конструкции. Совокупность единичных состояний метода сил называют базисом этого метода. Базисы, о которых шла речь в п. 5.7.1–5.7.2, именуются локальным и ортогональным соответственно. Построение последнего из них в общем случае связано с большим объемом вычислений, значительно превышающим тот, что необходим для решения задачи с первым попавшимся базисом, а потому к ортогонализации единичных эпюр прибегают лишь тогда, когда никакими другими средствами обеспечить требуемую точность расчета не удается. Однако частичную ортогонализацию базиса, основанную на взаимной ортогональности клас-

сов симметричных и обратносимметричных состояний основной системы, следует применять для любой симметричной конструкции.

Если конструкция симметрична, то любое воздействие на нее, в том числе и единичное, можно разбить на симметричную и обратносимметричную составляющие. На рис. 5.24b указаны такие составляющие для неравных в общем случае сил R1 и R2, приведенных на рис. 5.24a. Это составляющие X1 = (R1 + R2)/2 и X2 = (R1 − R2)/2. Используя данный прием, можно все основные неизвестные метода сил разбить на группы симметричных и обратносимметричных воздействий, вследствие чего отвечающие им единичные состояния основной системы будут взаимно ортогональными.

Рама, изображенная на рис. 5.25a, имеет четыре лишние связи. Разбиение реакций этих связей на симметричные X1, X2 и обратносимметричные

X3, X4 составляющие приведет к равенствам u13 = u14 = u23 = u24 = 0. В итоге система канонических уравнений четвертого порядка разобьется на

две подсистемы второго порядка каждая. Для конструкции, изображен-

ной на рис. 5.25b, возможно такое представление основных неизвестных, при котором первые четыре единичные состояния будут симметричными, а три последние – обратносимметричными. В результате вместо системы семи разрешающих уравнений получатся одна система четвертого порядка и одна система третьего порядка.

Глава 5

475

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

5.8. Линии влияния усилий. В статически неопределимых конструкциях нельзя построить линии влияния усилий, опираясь только на уравнения равновесия. Имеется еще одно усложнение по сравнению со статически определимым случаем: теперь все участки искомых линий влияния имеют нелинейное очертание.

Усилие SA, возникающее в сечении A при действии вертикальной единичной силы, положение которой на конструкции характеризуется параметром t, может быть найдено по формуле

SA(t) = S0A(t) + S1AX1(t) + · · · + SnAXn(t).

Через S0A(t) обозначено усилие, порождаемое воздействием P (t) = 1 в основной системе метода сил, а S1A, . . . , SnA – значения названного усилия от воздействий X1 = 1, . . . , Xn = 1 соответственно. Графики левой и правой частей любого равенства одинаковы, т. е.

л. в. SA = л. в. S0A +

 

1A л. в. X1 + · · · +

 

nA л. в. Xn.

(5.31)

S

S

Таким образом, для построения линии влияния усилия в статически неопределимой конструкции необходимо иметь линии влияния всех основных неизвестных. Составляющая л. в. S0A искомой линии влияния может быть получена по правилам, изложенным в главе 4.

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

DX + B = 0,

(5.32)

где D – матрица коэффициентов (единичных перемещений), а

X =

..

 

,

B =

..

 

 

 

X.

1

 

 

 

U.10

 

 

Xn

 

Un0

 

 

 

 

 

 

 

 

 

– столбцы реакций лишних связей и свободных членов соответственно. Так как положение подвижной силы характеризуется параметром t, то зависят от этого параметра и перемещения U10, . . . , Un0. Следовательно, столбец B является функциональным. При вычислениях функциональные матрицы так или иначе заменяются на числовые. Например, загружаемую часть осевого контура конструкции можно разбить на участки и определять перемещения Ui0 при приложении единичной силы на границах участков. Пусть точек

476

Часть IV

приложения этой силы будет m. Вычисленные значения перемещений вносятся в матрицу B , которая, в отличие от массива B, состоит не из одного функционального столбца, а из m численных столбцов. Для удовлетворительного описания линий влияния каждый загружаемый стержень конструкции требуется разбить на 6–10 частей, поэтому матрица B и определяемая ею матрица

X = −D1B

(5.32a)

имеют большие размеры.

Другой способ численного представления зависимости Ui0(t) состоит в разложении перемещения Ui0(t) в ряд по так называемым координатным

функциям fk(t):

 

 

k

 

Ui0(t) =

akfk(t).

(5.33)

 

=1,2,...

 

В матрицу B помещаются коэффициенты ak представления (5.33). При вычислениях удерживают конечное число слагаемых. Такой подход к решению задачи целесообразен тогда, когда число p сравнительно невелико, например, для конструкции, состоящей из призматических стержней. На каждом призматическом стержне перемещение (5.33) представляет собой конечный полином не выше четвертого порядка.

В качестве примера можно рассмотреть задачу построения линий влияния усилий в неразрезной балке (рис. 5.26). В первую очередь должны быть получены линии влияния основных неизвестных X1, X2 и X3. Единичная сила последовательно размещается во всех четырех пролетах балки. Каждый раз ее положение отмечается параметром t, отсчитываемым от левой опоры соответствующего пролета. При выбранной основной системе свободные члены канонических уравнений определяются формулами (см. также рис. 5.27, на котором индекс j отвечает номеру загружаемого пролета):

(j)

Uj−1,0

Uj(0j)

=

 

1

 

 

1

 

t(lj − t)

l

j ·

2lj − t

=

 

t(lj − t)(2lj − t)

,

EI 2

 

lj

 

3lj

6ljEI

 

 

·

 

 

 

 

=

 

 

1

 

1

 

t(lj − t)

l

j ·

 

lj + t

=

t(lj − t)(lj + t)

.

 

 

 

 

 

lj

 

 

 

 

 

EI 2

 

·

 

 

3lj

 

 

6ljEI

 

Для того, чтобы сделать дальнейшие записи более компактными, следует перейти к безразмерной координате τ = t/lj. В этом случае

Глава 5

477

(j)

l2

 

 

(j)

 

l2

 

 

j

 

 

 

j

 

 

Uj−1,0 =

 

[2τ (1

−τ )−τ 2(1 − τ )], Uj0

=

 

[τ (1−τ )+τ 2

(1−τ )].

6EI

6EI

Пусть lj2/6EI = bj

и

 

 

 

 

 

Тогда

 

f1(τ ) = τ (1 − τ ), f2(τ ) = τ 2(1 − τ ).

(5.34)

 

 

 

 

 

 

 

 

(j)

= bj(2f1

(j)

= bj(f1 + f2).

 

 

Uj−1,0

f2), Uj0

 

 

Так как перемещений U00(1) и U40(4) при n = 3 быть не может, то матрица B , содержащая множители при функциях f1 и f2 на каждом из четырех пролетах балки, имеет вид:

B =

 

1

b2

b2

2b3

b3

 

 

 

.

 

b

 

b1 2b2

b2

b3

b3

2b4

 

b4

 

 

 

 

 

 

 

 

 

 

 

Пусть не только изгибные жесткости EI стержней, но и их длины lj будут одинаковыми. Тогда b1 = · · · = b4 = = l2/(6EI) и

B =

l2

 

1

1

1

1

2

1

 

 

 

. (5.35)

 

6EI

 

 

2

1

1

1

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При lj = l (см. рис. 5.26)

u11 = (M 1, M 1) = 4l/(6EI), u12 = (M 1, M 2) = l/(6EI), u13 = 0,

u22 = u33 = u11, u23 = u12.

Стало быть,

 

 

 

1

 

D1 = 56l

 

 

 

 

4 .

D = 6EI

1

4

,

4

16

 

l

 

4

1

0

 

 

6EI

15

4

1

 

 

 

0

1

4

 

 

 

 

1

 

4

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При матрице B , определяемой равенством (5.35), вычисления по формуле (5.32a) дают:

X = 56

 

 

4

 

4

8

20

 

28

 

20

8

4

.

(5.36)

 

l

 

 

15

 

15

 

26

19

 

7

5

2

1

 

 

 

 

 

 

1

 

 

1

2

5

 

7

 

19

30

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

478 Часть IV

Теперь можно в явном виде записать уравнения линий влияния всех основных неизвестных. По данным первой строки матрицы X получают уравнения всех четырех участков л. в. X (номеру участка отвечает верхний индекс у символа X):

(1)

 

=

 

 

 

15l

+ f2),

X1

 

 

 

 

(f1

 

56

(2)

 

 

l

 

 

 

 

 

 

 

 

X1 =

 

 

 

 

(26f1

+ 19f2),

 

56

(3)

 

 

 

 

 

l

(5.37)

 

 

 

 

 

 

X1

 

=

 

 

 

(7f1 5f2),

 

 

56

(4)

=

 

 

l

+ f2).

X1

 

(2f1

56

Уравнения л. в. X2 и л. в. X3 получают по второй и третьей строкам матрицы X .

По формулам (5.34) видно, что функции f1(τ ) и f2(τ ) по концам каждого стержня, т. е. при τ = 0 и τ = 1, обращаются в нуль. Поэтому при построении графиков зависимостей (5.37)

необходимо вычислять ординаты лишь в промежуточных точках каждого пролета. При выполнении таких вычислений удобно пользоваться заранее подготовленными таблицами (табл. 5.1). Линии влияния основных неизвестных рассматриваемой балки приведены на рис. 5.28a.

Имея линии влияния основных неизвестных, можно приступать к построению линий влияния любых усилий. Так, линия влияния изгибающего момента MA, действующего в середине второго пролета балки, определяется равенством (см. формулу (5.31) и эпюры "M i" на рис. 5.26):

1

л. в. MA = л. в. M0A + 2 (л. в. X1 + л. в. X2).

Линия влияния изгибающего момента в сечении A основной системы изображена на рис. 5.28b. Добавление к этому графику полусуммы л. в. X1 и л. в. X2 приводит к искомой линии влияния (рис. 5.28c). Можно указать и аналитические записи всех ветвей л. в. MA. Для этого надо воспользоваться первой из формул (4.5) и данными, содержащимися в двух первых строках матрицы (5.36):

(1)

=

11l

(f1

+f2),

(2)

=

 

l

[1−τ −(12τ )θ(1,

2τ )]

l

(34f1

+f2),

MA

 

MA

 

 

 

 

 

112

2

112

 

 

 

 

 

3l

 

 

 

 

 

 

 

3l

 

 

 

 

(5.38)

 

 

 

(3)

 

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

MA

=

 

 

(7f1 +5f2), MA

=

 

(2f1 +f2).

 

 

 

 

112

112

 

 

ωA(j)
ak

Глава 5

 

 

 

 

 

 

 

 

479

 

 

 

 

 

 

 

 

 

 

Таблица 5.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fi

 

 

Безразмерная абсцисса τ

 

 

 

 

 

 

 

 

рассматриваемой точки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1/8

1/4

 

3/8

1/2

5/8

 

3/4

 

7/8

 

 

 

 

 

 

 

 

 

 

 

 

 

f1

7/64

3/16

 

15/64

1/4

15/64

 

3/16

 

7/64

 

 

f2

7/512

3/64

 

45/512

1/8

75/512

 

9/64

 

49/512

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если нагрузки вертикальны, то по готовой линии влияния усилия можно найти его значение, воспользовавшись формулами (4.10)–(4.13). При учете равномерно распределенной нагрузки необходимо знать площадь загружаемого участка линии влияния. Пусть нагрузка занимает весь пролет балки, а уравнение л. в. SA на данном пролете имеет вид:

SA(j) = a(kj)fk(τ ). (5.39)

k=1,2,...

Искомая площадь ωA(j) равна интегралу от этой функции:

= (j) fk(τ )dt, dt = lj dτ.

k=1,2,... lj

Линии влияния усилий в неразрезной балке определяются двумя функциями (5.34), поэтому

 

 

1

 

 

 

 

1

 

ωA(j) = lj

 

a1(j)

τ (1 − τ )+ a2(j)

τ 2(1 − τ )

 

 

0

 

 

 

 

 

0

 

или

 

ωA

 

=

12

2a1

+ a2

.

(5.40)

 

 

 

 

 

(j)

 

 

lj

(j)

(j)

 

 

Нетрудно найти и производную от функции (5.39) по параметру t, которую нужно знать при загружении линий влияния сосредоточенными моментами. Так как dt = lj, то

 

dSA(j)

= tg α(j) =

1

 

 

a(j)

dfk

.

 

 

 

 

k

 

 

 

dt

 

 

 

lj

k

 

 

 

 

 

 

 

 

=1,2,...

 

 

 

Для неразрезных балок

 

a1

(1 2τ ) + a2

τ (3 2τ ) .

 

tgα(j) = lj

(5.41)

 

 

1

(j)

 

 

 

(j)

 

 

 

 

480

 

 

 

 

 

 

 

 

 

 

Часть IV

При τ = 0 и τ = 1 формула (5.41) дает

 

 

 

 

 

 

tg α

(j)

(0) =

a(j)

, tg α

(j)

(1) =

a(j) + a(j)

.

(5.41a)

 

lj

lj

2

 

 

 

1

 

 

1

 

 

Пусть в рассматриваемой балке требуется найти изгибающий момент от нагрузки, приведенной на рис. 5.29. Уравнения ветвей л. в. MA известны: они даются формулами (5.38). На первом пролете

a(1)1 = a(1)2 = 11l/112

и, согласно равенству (5.40),

(1)

=

11l2

ωA

 

.

448

Сосредоточенный момент приложен к концу четвертого пролета. В этом месте a(4)1 = 6l/112, a(4)2 = 3l/112 и в соответствии со второй из формул (5.41a)

tg α(4)(1) = 1123 .

Силе P = 2ql отвечает на л. в. MA ордината yA = 155l/896. Таким образом,

= q ·

MA = A(1) + 2ql · yA + 3ql2 · tg

α(4)(1) =

 

448

 

+ 2ql ·

896

+ 3ql2 · −

112

 

=

112 ql2.

(5.42)

 

 

11l2

 

 

155l

 

3

 

 

27

 

 

В заключение несколько слов о кинематическом способе построения линий влияния усилий в статически неопределимых конструкциях. Решение задачи начинается с устранения связи, воспринимающей искомое усилие SA. Далее по направлению удаленной связи задаются перемещением, которое приводит к изгибу конструкции. Изогнутые оси загружаемых стержней образуют линию, отличающуюся от л. в. SA только множителем. Полученная модель линии влияния состоит из нелинейных участков, а потому переход от м. л. в. SA к л. в. SA требует вычисления ординат в большом числе точек. Трудоемкость расчета здесь не меньше, чем при использовании аналитического метода. Но кинематический способ удобен для определения вида искомых линий влияния, а знать их форму полезно не только в целях зрительного контроля результатов вычислений. Ведь для того, чтобы указать опасное положение временной равномерно распределенной нагрузки на конструкции, зачастую достаточно одних лишь моделей линий влияния. Более подробно об этом говорится в следующем пункте.