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

1846

.pdf
Скачиваний:
0
Добавлен:
16.06.2024
Размер:
3.28 Mб
Скачать

наиболее опасной точке xm , ym , z имеем плоское напряжённое состояние

z 0, zx 0, zy 0

и напряжения x xm , ym , z

и y xm , ym , z

являются главными нормальными напряжениями.

Поскольку M x 0 и

M y 0 , то знак у напряжений x

и y будет

определяться знаком

координаты z , то есть на верхней грани пластины z

h

оба напряжения

 

 

 

2

 

будут отрицательными, на нижней z

h

 

– положительными. Учитывая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

значения моментов M x и M y

в точке xm , ym , получаем: x 2

и y 1 ,

поскольку, как известно, 1 2

3 и для пластин 3

0 .

 

 

1

В

соответствии

с теорией

прочности

Губера

– Мизеса

Генки

 

 

2

2

 

2

 

3

2

3

2

R ,

 

толщина

пластины

будет

 

 

 

 

2

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определяться по формуле

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

6

 

M x xm , ym 2 M x xm , ym M y xm , ym M y xm , ym 2 ,

 

 

 

 

 

 

R

 

и составит:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h 4,073 10 2

м 4,073 см.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.Определяем распределённые реакции. На основании соотношений

(34)получаем:

1)

0 x a,

y b и

y b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vx 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 3Y1 y

 

 

2

dY1 y

 

 

x

 

 

 

 

Vy x D

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

dy

3

 

 

 

 

dy

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 3Y3 y

 

 

 

3 2

dY3 y

 

 

3 x

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

sin

 

 

.

 

 

 

 

dy

3

 

 

 

 

 

 

dy

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

x 0

и x a, b

 

y b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

y D

 

3

 

 

 

 

y 2

d Y1 y

cos x

 

 

V

 

Y

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

1

 

 

 

a

 

 

dy

2

 

 

 

 

 

 

a

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

y

2 3

d Y3 y

cos 3 x ;

 

 

 

D

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

3

 

 

 

 

 

 

a

 

 

dy

2

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vy 0.

181

Вычисления показывают, что распределённые реакции вдоль контура пластины направлены вверх и одинаковы на противоположных сторонах пластины.

8. Определяем сосредоточенные реакции в углах пластины:

R x0 , y0 2H

 

 

 

 

 

 

 

 

 

 

 

 

2D 1

dY1

y0

 

 

x

2D 1

3 dY3 y0

 

3 x

 

 

 

 

cos

0

 

 

cos

0 .

a dy

 

a dy

 

 

 

a

 

 

 

a

В итоге получаем:

 

 

 

 

 

 

 

 

 

 

 

 

R 1,6 10 10

0 H;

R 1,6 10 10

0 H;

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

R 1,6 10 10

0 H;

R 1,6 10 10

0 H,

 

 

3

 

 

 

 

 

4

 

 

 

 

 

 

то есть сосредоточенные реакции в углах пластины отсутствуют.

Эпюры распределённых реакций по контуру пластины изображены на рис. 96.

 

-6,22

 

3,78

-6,22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,62

-3,49

 

-3,49

 

 

2,62

3,53

 

 

 

 

 

 

2,62

3,78

-6,22

 

3,53

 

-6,22

 

2,62

 

 

 

 

-3,49

 

-3,49

 

 

 

 

 

 

 

Рис. 96

Численные расчёты показывают, что распределённая реакция вдоль жёсткой заделки практически равна нулю.

182

Глава 5. МЕТОДЫ НА ОСНОВЕ ВАРИАЦИОННЫХ ПРИНЦИПОВ

§19. Вариационная формулировка задач теории упругости

В самом общем случае любая механическая задача может быть сформулирована с двух точек зрения:

1.С точки зрения движения механической системы и действующих на неё сил.

2.С точки зрения энергетического баланса движущейся механической системы.

Первая точка зрения приводит к математической формулировке задачи

ввиде дифференциальных уравнений. Вторая позиция математически формирует задачу в виде некоторых интегральных соотношений или функционалов. Оба подхода эквивалентны и приводят к одному и тому же результату. Разница заключается в трудоёмкости решения задачи.

Понятие о функционале

Пусть дан некоторый класс (множество) функций K y x , где х

независимая переменная или совокупность нескольких независимых переменных x x1, x2 , ..., xn .

Говорят, что переменная величина I I y x есть функционал от функции y y x , если каждой функции y x K по заданному закону

ставится в соответствие определённое число I .

Замечание. Функция каждому числу ставит в соответствие другое число; функционал каждой функции ставит в соответствие определённое число; производная каждой функции ставит в соответствие другую функцию.

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

П р и м е р 1 .

1) Площадь криволинейной трапеции, как известно, определяется

b

соотношением S y x dx , то есть является функционалом. Требуется: для

a

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

183

b

dy x

2

2) Длина дуги вычисляется по формуле L

dx , то есть

1

dx

 

a

 

 

 

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

или минимальна.

Рис. 97

Уравнения теории упругости в матричной форме

Введём следующие векторы-столбцы: u u,v, w T ;

 

 

 

 

 

x , y , z , xy , yz , zx T ;

 

ε

 

 

 

 

x , y , z , xy , yz , zx T ;

(1)

σ

 

 

Fx , Fy , Fz T ;

 

F

 

 

pnx , pny , pnz T .

 

P

 

При этом основные группы соотношений теории упругости можно, очевидно, записать в следующем виде.

1. Геометрические соотношения:

 

 

 

 

 

AT

 

 

.

 

 

 

 

 

 

(2)

 

 

ε

u

 

 

 

 

 

 

Здесь матрица A имеет следующую структуру:

 

 

 

 

0

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

0

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

y

 

 

 

x

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

y

 

 

 

 

 

 

 

 

 

 

 

x

184

2. Уравнения равновесия:

 

A

 

 

 

 

0.

(3)

σ

F

3. Физические уравнения.

 

a. В форме закона Гука:

 

 

 

С

 

.

(4)

ε

σ

Здесь матрица упругой податливости С имеет следующую структуру:

 

 

 

1

 

 

0

0

0

 

 

 

 

 

1

 

0

0

0

 

 

 

 

 

С

1

 

 

 

1

0

0

0

 

 

 

0

0

0

2 1

0

0

.

 

 

E

 

 

 

 

0

0

0

0

2 1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

0

 

 

 

 

 

2 1

b. В форме Ламе:

σ D ε.

Здесь матрица жёсткости D равна:

K

0

 

4 G

K

0

 

2 G

K

0

 

2 G

 

 

3

0

 

 

3

0

 

 

3

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 G

 

 

 

4 G

 

 

 

2 G

K

 

 

K

 

 

K

 

 

 

0

 

3

0

 

0

 

3

0

 

0

 

3

0

 

 

 

2 G

 

 

 

2 G

 

 

 

4 G

D K

 

 

K

 

 

K

 

 

 

0

 

3

0

 

0

 

3

0

 

0

 

3

0

 

 

0

 

 

 

0

 

 

 

0

 

 

 

0

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

(5)

0

0

0

 

 

 

 

 

 

0

0

0

 

 

 

 

 

0

0

0

.

 

 

 

 

G

0

0

 

0

G0

 

 

0

0

 

0

0

G

 

 

 

0

 

4. Уравнения совместности деформации:

 

B ε 0 .

(6)

185

Здесь матрица B определяется следующим образом:

 

 

 

 

2

 

 

 

 

2

 

 

 

0

 

 

 

2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2

 

 

 

x

2

 

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

2

 

 

 

 

2

 

0

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

 

y2

 

 

 

 

y z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

0

 

 

 

 

 

 

2

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

2

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y z

 

 

 

 

 

 

z x

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

0

 

 

 

2

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z x

 

 

y z

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

2

 

 

 

0

 

 

 

 

0

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

z

2

 

z x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5. Условия на поверхности:

0

 

 

 

 

0

 

 

 

 

2

2 .z xx y

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y z

 

 

L

 

 

 

.

 

 

 

(7)

 

 

σ

P

 

 

 

Здесь матрица L имеет следующую структуру:

l

0

0

 

m

0

n

 

 

0 m

0

 

l

n 0

 

;

L

 

 

 

0

0

 

n

0

m

l

 

 

 

 

 

 

l, m, n направляющие косинусы нормали к рассматриваемой площади

поверхности тела.

При записи уравнений теории упругости в матричной форме решение прямой задачи теории упругости в перемещениях сводится к интегрированию системы дифференциальных уравнений вида:

 

 

T

 

 

 

 

 

 

 

A

 

 

F 0.

(8)

D A

 

u

Решение прямой задачи теории упругости в напряжениях сводится к интегрированию системы дифференциальных уравнений следующего вида:

B C σ 0,

(9)

A σ F 0.

Решения системы уравнений (8) и (9) должны удовлетворять условиям на поверхности в перемещениях и напряжениях (7).

186

Вариационный принцип Лагранжа

Полная потенциальная энергия деформированного тела "Э" состоит из потенциальной энергии деформации тела (потенциал внутренних сил) "U"

и энергии внешних сил (потенциал внешних сил) "П":

 

Э U П .

(10)

Будем считать, что в начальном, недеформированном, состоянии полная потенциальная энергия деформированного тела Э0 0 . Следовательно,

полная энергия деформированного тела представляет собой изменение энергии внутренних и внешних сил при переходе тела из начального в деформированное состояние.

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

Удельный потенциал внутренних сил определяется, как известно, соотношением

U0

 

1

 

T

 

.

(11)

σ

ε

 

 

2

 

 

 

 

 

Полный потенциал внутренних сил будет равен:

U U0 dV

1

 

 

T

 

dV .

(12)

σ

ε

V

2

V

 

Внешние силы, определяемые вектором p, и объёмные силы, задава-

емые вектором F, на перемещениях, соответствующих переводу тела из деформированного в недеформированное состояние, совершают работу, равную:

П

 

T

 

 

dS

 

T

 

dV .

(13)

p

u

F

u

S

 

 

V

 

Отметим, что при этом работа внутренних сил будет положительной, а работа внешних и объёмных сил – отрицательной. То есть соотношение (10) получает вид:

Э U0

dV

 

T

 

 

dS

 

T

 

dV .

(14)

p

u

F

u

V

S

 

 

V

 

187

Используя закон упругости (5), геометрические соотношения (2), соотношение (11), формулу (14) можно записать в терминах функций перемещений:

Э 1

AT

 

T

D AT

 

 

dV

 

T

 

 

dS

 

T

 

dV . (15)

u

u

p

u

F

u

2

 

 

 

 

S

 

 

V

V

 

 

 

 

 

Итак, полная потенциальная энергия деформированного тела вполне определяется заданием функций перемещений u, v, w, зависящих, вообще

говоря, от пространственных координат.

Принцип Лагранжа формулируется так: для истинных перемещений u, v, w функционал полной энергии деформиро-

ванного тела Э имеет экстремальное (стационарное) значение, то есть его первая вариация равна нулю.

Согласно этому принципу, если u истинные перемещения точек тела,

при которых имеет место равновесие, то работа внешних и внутренних сил на произвольном бесконечно малом изменении перемещений

u u, v, w T , допускаемом связями тела, должна быть равна нулю, то есть

A Aвнутр Aвнешн 0.

(16)

Так как работа внутренних сил пропорциональна изменению потенциала внутренних сил, а работа внешних сил пропорциональна изменению

потенциала внешних сил, то есть Aвнутр U ; Aвнешн П, то на основании (16) получаем: A Э.

Итак, для истинных перемещений u изменение полной потенциальной энергии деформированного тела Э, вызванное вариациями перемещений

uT , должно быть равно нулю:

 

Э Э u u Э u 0 .

(17)

Соотношение (17) представляют в виде ряда, в котором каждое

слагаемое зависит от соответствующей степени u :

 

Э Э1 u Э2 u2 ... 0 .

(18)

Здесь Э Э1 u первая вариация функционала энергии;

2Э Э2 u2 вторая вариация функционала энергии.

188

Устремляя в (18) вариации перемещений u к нулю и отбрасывая все слагаемые, кроме первого, как бесконечно малые более высокого порядка малости, приходим к равенству:

Э 0 .

(19)

Равенство нулю первой вариации функционала (19) является необходимым условием локального экстремума этого функционала, то есть в локальной зоне изменения функций-аргументов, функционал с точностью до бесконечно малых первого порядка сохраняет неизменное (стационарное) значение. В соответствии с теоремой Лежен Дирихле:

 

при

Э 0

и

2Э 0

полная энергия деформированного тела

минимальна и его равновесие устойчиво;

 

при

Э 0

и

2Э 0

эта энергия максимальна и равновесие

системы внутренних и внешних сил неустойчиво;

 

при

Э 0

и

2Э 0

энергия стационарна, а тело находится в

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

Итак, решая задачу теории упругости с использованием принципа Лагранжа, следует найти такие функции перемещений u, v, w, зависящие, во-

обще говоря, от пространственных координат, при которых первая вариация полной потенциальной энергии деформированного тела будет равна нулю:

Э 0 .

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

В принципе Лагранжа варьируемыми функциями служат перемещения

u, а соответствующий функционал (15) определяет полную потенциальную энергию деформированного тела.

В принципе Кастильяно варьируемыми функциями являются напряже-

ния σ, а соответствующий функционал определяет дополнительную энергию деформируемого тела:

Эk

 

T

 

 

dV

 

TS2

 

S2 dS2 .

(20)

ε

σ

u

p

V

 

 

S2

 

Здесь S2 часть поверхности тела, на которой заданы перемещения.

Принцип Кастильяно говорит о том, что истинные напряжения сообщают дополнительной энергии деформируемого тела стационарное значение.

189

В принципе Рейсснера в качестве варьируемых функций приняты функции и перемещений u, и напряжений σ. Функционал Рейсснера имеет вид:

Эр

1

 

 

T

C

 

dV

 

T A

 

 

 

 

dV

 

σ

σ

u

σ

F

 

 

2

V

 

 

 

 

 

 

V

(21)

 

 

T L

 

 

 

S1 dS

 

TS2 L

 

 

dS.

u

σ

p

u

σ

 

 

S1

 

 

 

 

 

 

 

 

S2

 

Здесь S1 часть поверхности тела, на которой задана нагрузка.

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

В принципе Ху и Вашицу в качестве независимо варьируемых функций

приняты и перемещения u, и напряжения σ, и деформации ε Ху – Вашицу можно записать в следующей форме:

ЭХ В

1

 

 

 

T D

 

dV

 

 

T AT

 

 

 

 

 

dV

ε

ε

σ

u

ε

 

2

V

 

 

 

 

 

 

V

 

 

T

 

 

dV

 

T

 

dS

 

 

 

S2 T L

 

 

u

F

u

p

u

u

σ

 

 

V

 

 

 

 

S1

 

 

 

 

 

S2

. Функционал

(22)

dS.

Принцип Ху – Вашицу говорит о том, что для истинных перемещений, напряжений и деформаций первая вариация функционала Ху – Вашицу равна нулю.

Решение задач теории упругости на основе сформулированных принципов сводится к отысканию экстремума соответствующего функционала.

§20. Метод Ритца

В 1909 году Вальтер Ритц, швейцарский физик-теоретик и математик, предложил метод решения общей задачи об экстремуме функционала. Метод предусматривает выбор пробной функции, которая должна минимизировать определённый функционал, в виде суперпозиций известных функций, которые удовлетворяют граничным условиям. При этом задача

сводится к поиску неизвестных коэффициентов суперпозиции.

 

Пусть для функции

x, y, z

,

составляющей

оператор

R , задан

функционал

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

,

 

,

2

 

dV .

(1)

I R x, y, z, ,

x

y

z

x

2

,...

V

 

 

 

 

 

 

 

 

 

190

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]