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

лекции Каюмова по сопромату

.pdf
Скачиваний:
101
Добавлен:
12.05.2015
Размер:
2.9 Mб
Скачать

Отсюда:

l

M x(T ) M x(q)

2

 

1

l

(T )

(q)

2

 

Jx

 

l

(T )

(q)

 

W ds

 

y

dA

 

dsM x

M x

y

dA

 

 

M x

M x

ds .

EJ 2

EJ 2

EJ

2

0 A

x

 

 

x

0

 

 

A

 

 

x

0

 

 

 

Jx

Запишем закон сохранения энергии:

 

 

 

 

 

U W .

 

 

Отсюда вытекает формула Мора:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

l

 

 

 

 

 

 

T B

 

(T )

(q)

ds .

 

 

 

EJ

M x

M x

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

 

 

 

Перечислим использованные обозначения.

vB

- искомый прогиб в точке B (от рабочих нагрузок);

(16.2)

T =1 – единичная сила, приложенная в интересующем нас направлении

искомого прогиба

vB

в интересующей нас точке В.

M x

- изгибающий момент в фиктивной задаче о приложении к балке

(T )

 

силы Т в точке В.

M x

- изгибающий момент от рабочих нагрузок.

(q)

 

Еще раз напомним, что физический смысл формулы Мора заключается в том, что работа силы Т на искомом перемещении vВ равна работе внутренних сил, вызванных этой силой, на деформациях от внешних сил.

Примечания.

1.Для удобства вычислении обычно принимают, что Т=1.

2.Работой касательных напряжений обычно пренебрегают ввиду ее малости по сравнению с W.

3. При необходимости вычисления угла наклона балки α вместо единичной фиктивной силы Т необходимо прикладывать единичный момент m=1 в интересующей нас точке. Формула Мора примет вид

m

 

 

1

B

EJ

 

 

 

 

 

 

 

 

x

l

0

M

(m)

M

(q)

ds

x

x

 

 

 

.

16.1.1 Методы вычисления интегралов. Формулы трапеций и Симпсона

Для приближенного вычисления интегралов существует много разных методов. Пусть надо найти:

m

I y(x) dx ?

n

101

Для вычисления интеграла Мора часто используют метод Верещагина. Однако более удобными являются приближенные методы.

Рис. 16.5

Разобьем интервал

Формула трапеций

n, m

на малые интервалы (например, на рис.16.5.

их четыре). Поскольку по геометрическому смыслу интеграл представляет собой площадь фигуры mnkl, то С можно вычислить приближенно, представив ее в виде суммы площадей четырех трапеций:

I

y y

 

x1

y

y

x2

 

y

y

 

x3

 

y

y

x4 .

1

2

2

3

3

 

4

4

5

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

2

 

 

 

 

2

 

Формула Симпсона

Формула Симпсона намного точнее формулы трапеций показаться менее удобной). Она имеет вид:

(16.3)

(хотя может

I

a

y

4 y

 

 

 

2

6

1

 

 

 

 

 

 

y

 

в

y

4 y

 

3

 

6

3

4

 

 

 

 

y5

.

(16.4)

При этом в отличие от метода трапеций, отрезок а должен разбиваться

на равные интервалы x1 x2

. Аналогично, должно быть x3

x4

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

формулу:

 

 

 

 

 

I

a в

y1 4 y3 y5 .

 

 

6

 

 

 

 

 

16.2. Вычисление прогибов на основе решения дифференциального уравнения изогнутой оси балки

Прогибы можно находить и другими способами, например, на основе решения дифференциального уравнения изогнутой оси балки. Для вывода этого уравнения, рассмотрим элемент балки (рис.16.6).

102

Рис. 16.6

Ясно, что чем больше M x , тем больше кривизна Эту фразу можно записать в виде:

1

k M x .

 

 

1

изогнутой оси балки.

(16.5)

Выразим кривизну через прогиб. анализа:

1

 

 

 

 

1

 

 

Согласно формулам математического

v

 

 

3 .

v

2

 

2

 

 

 

 

 

 

Рис.16.7

По геометрическому кривой (рис16.7):

Ввиду малости угла

Тогда:

смыслу производная - это тангенс угла наклона

v

 

tg .

 

 

 

 

 

 

 

 

 

можно записать:

 

 

 

 

 

 

 

 

 

 

 

 

 

tg v 1.

 

 

 

 

 

1

 

v

 

 

 

 

 

 

1 v 2 32

v

k M x .

(16.6)

 

 

Очевидно, что k зависит от геометрии сечения и материала балки. Найдем эту зависимость.

Рассмотрим малый элемент балки длины ds (рис. 16.8). После изгиба он превратится в изогнутый элемент. Длина волокна BC, которое проходит через центр тяжести сечения, не изменяется и будет равна ds . Нижнее волокно DH удлиняется на величину, которое обозначим через .

103

Рис.16.8

Вычисляем

, учитывая, что

y

D

 

 

 

DH

h

.

2

 

ds

Согласно определению

.

Используя закон Гука и формулу Навье получаем

 

 

 

DH

 

M

x

 

 

M

x

 

h

 

 

DH ds

 

ds

 

yD ds

 

 

 

 

ds .

 

 

 

EJ x

EJ x

 

 

 

 

 

E

 

2

 

Вычислим

 

теперь по

другому -

через

угол

геометрии известна формула для вычисления длины дуги:

d

(16.7)

(рис.16.8). Из

ds d .

Тогда DH BC ( h / 2)d d

h 2

ds

.

(16.8)

Приравниваем (16.7) и (16.8):

h ds

 

M

x

 

 

 

2

 

EJ

x

 

 

 

h 2

ds

.

Отсюда получаем:

1

 

 

M x

.

 

 

 

 

 

 

 

E Jx

Учитываем, что согласно (16.6):

 

 

 

 

 

 

1

v

 

 

 

 

 

 

 

 

 

Окончательно получаем:

v

M

x

 

 

 

 

 

 

 

E J

x

 

 

 

Это и есть уравнение изогнутой оси балки.

(16.9)

16.2.1. Решение дифференциального уравнения изогнутой оси балки

Если балка имеет постоянную толщину, то есть Jx const , то решение легко записывается в общем виде:

104

1 EJ

1 EJ

x

M x ds C , x

ds M x ds Cs

D

.

(16.10)

(16.11)

Хотя решение получено в общем виде, однако основная трудность заключается в определении Мх и констант C и D, поскольку на разных участках балки M x разные, а значит C и D также разные (в частности, если

балка имеет три участка, то нужно определить 6 констант).

Однако существует способ интегрирования, который сводит все неизвестные только к двум константам (разработан Клебшом)

16.2.2. Правила Клебша

Правила Клебша сводятся к следующему.

1) M x выражаем через внешние силы, которые лежат только слева (или

только справа) от сечения.

2) Если погонная сила q не доходит до правого конца, то ее доводим до этого правого конца и уравновешиваем ее снизу (рис.16.9).

Рис.16.9 3) Если имеется сосредоточенный

изгибающий момент записываем в виде m0

момента mо.

4) Интегрируем, не раскрывая скобок.

момент

mо, то

(s a)

0

,

где а -

 

 

 

его вклад в расстояние до

При выполнении этих условий все константы С на разных участках будут одинаковы. Аналогично будут одинаковы все константы D.

Справедливость этого утверждения доказывается непосредственной проверкой, то есть подстановкой решения в условия стыковки решения на границе участков. Рассмотрим, например, случай, приведенный на рис.16.10.

рис16.10

По правилам Клебша момент M x на участках (I), (II) запишем в виде:

(I): M x P s

(II): M x P s m0 (s a)0

Дифференциальные уравнения на участках:

105

(I)

(II)

v

P s

 

 

 

 

 

EJ

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P s

 

m(s a)

0

v

 

 

EJ

 

EJ

 

 

 

 

x

 

x

 

 

 

 

 

 

 

 

Решение этих уравнений на участках (1), (2) имеет вид:

Участок (I):

Участок (II):

 

P s

2

 

 

 

 

 

P s

3

 

2EJ

 

C,

 

6EJ

Cs

 

x

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

P s

2

 

m(s a)

 

 

 

 

 

C,

2EJ

 

 

EJ

 

 

 

 

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

D

.

P s

 

3

6EJ

x

 

m(s a)2 2EJ x

Cs

D

.

Отсюда видно, что при s = a получим равенство углов наклона и прогибов, вычисленных по разным формулам при любых С и D, т.е. условия гладкости изогнутой оси выполняются. Аналогично проверяются условия гладкости на границе участка, на которой заканчивается погонная сила q.

16.2.3Условия для определения С и D

1)Первый случай . Рассмотрим балку, лежащую на двух опорах

(см. рис.16.11).

l

 

 

1

 

 

l

2

l

 

 

 

 

Рис. 16.11

Рис. 16.12

Из схемы видно, что

 

 

v(l1 ) 0

(16.13)

 

 

 

v(l2 ) 0

 

Таким образом из (16.13), получаем систему уравнений для С и D.

2) Второй случай. Пусть балка

заделана на расстоянии l

(консольная балка, см. рис.16.12).

 

В заделке не может появиться наклона оси, поэтому там не только нет

 

 

прогиба, но и v 0 .

 

Таким образом, из схемы следует, что:

 

v(l) 0

(16.14)

 

v (l) 0

 

Опять получили два уравнения для С и D.

106

Пример вычисления прогиба

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

Ввиду симметричности схемы можно сразу найти реактивные силы – они будут равны ql/2. Тогда изгибающий момент в сечении на расстоянии s

от левой опоры будет равен

 

 

 

l

 

 

s

2

M

 

q

s q

 

x

 

 

2

 

 

 

2

 

 

Первый способ. Использование дифференциального уравнения изогнутой оси балки.

 

 

 

 

 

 

 

 

 

 

 

 

 

EJx v

 

 

 

 

l

 

 

 

 

 

 

s2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

s q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируем 2 раза:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EJ

 

v q

l

 

s2

 

 

q

s3

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EJ

 

v q

l

 

s3

 

q

s4

 

Cs D

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

6

 

 

 

 

 

24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Константы

 

 

 

 

интегрирования

 

 

находим из

условий закрепления:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v 0 0

 

 

 

 

 

0

4

 

q

0

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

C 0

D 0

 

 

 

D 0.

 

 

12

 

24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v l 0

 

 

 

 

 

 

l

4

 

 

ql

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ql

4

 

 

 

q

 

 

 

 

 

Cl 0

 

 

 

 

 

 

 

 

 

 

C

 

 

 

12

24

 

 

 

 

 

 

 

 

 

24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим прогиб в центре балки (при s = l/2):

 

 

l

 

 

4 ql

4

 

 

 

 

 

ql

4

 

 

 

 

 

8 ql

4

1

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

EJ

v

 

 

 

 

 

 

 

 

 

 

 

 

ql

4

 

 

0, 01302ql

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

12 8

 

 

 

 

24 16

 

 

 

 

24

 

 

 

2

 

 

 

 

 

 

384

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Второй способ. Использование интеграла Мора

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

Прогиб в центре балки находим по формуле

 

 

EJ

 

 

v

 

M

 

M

(T )

ds

.

 

 

 

 

 

 

 

x

x

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нарисуем эпюру изгибающих моментов Мх(T) от единичной силы Т=1

(см. рис. 16.13).

Рассмотрим различные приближенные методы интегрирования. 1. Метод трапеций по 2-м участкам.

 

ql2

8 l 4

 

l

 

ql2

8 l 4

 

l

 

ql4

 

 

 

 

 

 

 

 

EJx v (0

 

 

)

 

(

 

 

0)

 

 

 

0,0156ql4 .

 

2

2

 

2

2

64

 

 

 

 

 

 

 

 

 

 

Метод дал ошибку в 17% 2. Метод трапеций по 4-м участкам.

 

1

 

 

 

3

 

 

 

l

l

 

1

 

3

 

 

 

l

 

ql

2

0

 

ql

2

 

 

ql

2

 

 

 

EJx v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

32

 

 

 

8

 

4

 

2

 

32

 

 

 

8

 

8

Метод дал ошибку в 5%.

3. Метод Симпсона по 2-м участкам.

l

 

l

3 3

8

 

4

 

 

 

2

 

 

ql

4

4

128

8

 

 

 

 

 

0, 01367ql

4

 

.

107

 

 

 

 

3

 

 

 

l

 

ql

2

 

l

l

 

5ql

4

EJx v 2

0

4

ql

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32

 

 

 

8

 

8

 

4

6

 

384

.

Таким образом, метод Симпсона в этом примере дает точное решение.

16.2.4. Балки на упругом основании. Закон Винклера

Рассмотрим основание (грунт) Приложим к мерному стержню силу Р (рис.16.14), тогда осадка стержня v , будет тем больше, чем больше сила Р.

Рис.16.14

Запишем это утверждение аналитически:

P

k

v

.

Здесь

k - называется коэффициентом постели. Эта зависимость называется

законом Винклера.

 

Согласно III-му закону Ньютона реакция грунта R=P. Приведем закон

к виду, когда задаѐтся погонная сила q (cм. рис.16.15).

 

q

 

Р

 

a

R r

Рис.16.15

Ее равнодействующая должна быть равна Р, т.е.

q a P ;

 

Отсюда

 

 

 

q a k v ;

q

k

v .

a

 

 

 

108

Обозначая

k

 

 

k

получим

a

 

 

 

 

 

 

q k v .

Погонную реакцию грунта будем обозначать через r (см. Тогда:

(16.15)

рис 16.14).

r

R a

k v

.

(16.16)

16.4. Уравнение изогнутой оси балки на упругом основании

Рассмотрим балку, которая лежит на грунте.

q

r

Рис.16.16

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

Грунт противодействует внешним силам некоторой погонной силой r . Выразим ее через прогиб (s) . Для этого вырежем малый элемент (см. рис.16.16). Сравнивая рис.16.14 и рис.16.16 видим, что элемент балки, фактически представляет собой мерный стержень, для которого реакция определяется по формуле (16.16), т.е. реакция балки с прогибом связана соотношением:

r

k v

.

Далее запишем уравнения равновесия элемента балки и уравнение ее изогнутой оси

Q

y

 

q .

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

M

x

Qy .

 

 

 

 

 

 

 

 

 

 

s

 

 

 

1

 

 

 

v

 

 

 

 

 

2

 

 

 

 

 

 

 

 

s

2

 

 

EJ

 

M x .

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

Из рис.16.16 видно, что на балку действует две погонные силы: r . Тогда получим:

Qy qвнеш r qвнеш k v .s

(16.17)

(16.18)

(16.19)

qвнешн и

(16.20)

109

Здесь знак «-» перед

k

имеет отрицательный знак,

погонной силе q

внешн

.

 

 

'

v

 

 

а

поставлен потому, что осадка v элемента реакция должна быть противоположна

Подставим (16.18) в (16.20):

2 M x qвнешн k v .s2

Подставляя сюда (16.19) получаем искомое уравнение:

 

4

v

 

1

(q

 

 

 

 

 

 

 

 

внеш

s

4

 

EJ

 

 

 

x

 

 

 

 

 

k v)

.

(16.21)

Решение запишем в виде суммы:

v v

v

одн

част .

Простой подстановкой в (16.21) можно проверить, что решение имеет вид:

v

1

C1e

s

sin s C2e

s

cos s C3e

s

sin s C4e

s

cos s B . (16.22)

 

EJ

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

 

 

 

 

4

k

.

 

 

 

 

 

 

 

4EJ

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

Частное решение находим, подставляя

vчаст

= B

в уравнение (16.21):

 

B

 

1

 

 

 

 

4

 

 

 

 

(q

 

k B)

 

 

 

 

 

внеш

 

 

 

 

 

 

 

 

s

4

 

EJ

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

q

 

внеш

B

k

 

.

Остальные константы получают из геометрических соображений (условий закрепления) и условий статики на концах балки.

16.5. Бесконечная балка на упругом основании

Этой моделью можно описать, например,

полотна с автомобилем веса Р. Погонная сила

q

внеш

 

 

 

 

поведение дорожного qвеса представляет собой

погонный вес полотна. Мы можем общее решение представить как сумму

решения 2-х задач: задачи о действии только силы

веса qвеса и задачи о

действии только силы Р.

Здесь

vчаст

B

соответствует случаю

когда,

действует лишь qвеса.

 

 

 

 

 

 

Прогиб vодн , который

содержит,

C1,C2 ,C3 ,C4 ,

соответствует

случаю

q

внеш

 

0

, P 0 . Рассмотрим этот случай. Пусть s – расстояние от силы Р до

сечения. Слева и справа прогиб симметричный, поэтому исследуем прогиб v только справа, то есть, найдем функцию v(s).

Рис.16.17

110