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

Избранные разделы теории упругости. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
852 Кб
Скачать

Уравнение эллипса, ограничивающего пластинку:

 

x2

+

y2

1 = 0.

 

 

(5.1)

 

a2

 

 

 

 

 

 

b2

 

 

 

 

 

 

 

 

Зададим функцию прогибов в виде

 

 

 

 

 

 

 

 

 

 

 

x2

y2

2

(5.2)

w(x, y) =C

 

 

+

 

 

 

1

,

 

2

b

2

 

 

 

a

 

 

 

 

 

 

 

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

Покажем, что функция прогибов в виде (5.2) удовлетворяет граничным условиям, соответствующим заделке. Действительно, поскольку в выражение (5.2) входит как сомножитель левая часть уравнения эллипса, то сразу видно, что на контуре w = 0. Кроме того, первые производные функции w(x, y) имеют вид:

w

 

4Cx

 

2

 

 

y

2

 

 

 

=

 

x

 

+

 

 

 

1 ,

 

x

a2

 

 

b2

 

 

a2

 

 

 

 

w

 

4Cy2

 

2

 

 

 

2

 

 

 

=

 

x

 

+

y

 

1 .

(5.3)

y

 

2

2

 

b

a

 

 

 

b

 

 

 

 

В оба выражения опять-таки входит в качестве сомножителя левая часть уравнения эллипса (5.1). Следовательно, везде по контуру

wx = 0 и wy = 0.

Итак, функция прогибов в виде (5.2) при всех значениях параметра С автоматически удовлетворяет условиям заделки.

Определим теперь значение константы С, при котором удовлетворяется дифференциальное уравнение (4.31).

101

w

 

4Cx

 

 

x

2

 

 

y

2

 

 

 

 

 

4C

 

x

3

+ xy

2

 

=

 

 

+

 

1

=

 

 

 

x ;

x

a2

 

 

b2

a2

 

 

 

 

a2

 

 

 

 

 

 

 

a2

b2

 

 

 

 

 

2

w

 

4C

 

3x

2

 

y

2

 

 

 

 

 

 

 

 

 

 

 

=

 

 

+

 

 

1 ;

 

 

 

 

 

x2

 

a2

a2

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3w

=

4C 6x

 

 

 

 

 

 

 

 

 

 

 

 

 

x

3

 

a

2

 

 

2

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4w

 

=

24C

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x4

 

 

a4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4w

 

=

24C

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y4

 

 

b4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3w

=

4C

2y

 

 

 

 

2

 

 

2

;

x

2

y

a

 

 

 

 

b

 

 

3w = 8C . x2y2 a2b2

Подставим полученные выражения в уравнение (4.31)

24aC4 +2a82Cb2 + 24bC4 = Dq .

Уравнение удовлетворяется, если

 

C =

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

2

 

 

 

 

3

 

 

 

 

 

8D

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

Поэтому

 

 

 

 

 

4

 

 

 

2

b

2

b

4

 

 

 

 

 

 

a

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

y2

2

w =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

1 .

 

3

 

2

 

 

 

 

3

 

 

 

 

 

b2

 

 

 

 

 

 

a2

 

 

 

 

 

 

8D

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

2

b

2

b

4

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.4)

(5.5)

(5.6)

Полученное решение является точным в рамках сделанных предположений (см. § 1 главы IV). Теперь мы можем вычислить все интересующие нас величины (см. § 9 главы IV).

Наибольший прогиб будет в центре эллипса, в точке (0, 0). Он численно равен С. Теперь становится ясен физический смысл коэффициента С – величина максимального прогиба пластинки.

102

Найдем моменты в направлении оси Ox:

 

 

 

 

 

 

 

 

 

 

 

 

 

M x

 

 

 

 

 

 

 

 

 

2w

2w

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −D

x

2

y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

3x2

 

 

 

y2

 

 

 

 

 

 

 

4

 

 

x2

 

3y

2

 

 

 

 

 

 

 

= −D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

1

=

 

 

 

 

 

 

 

3

 

 

 

 

 

 

2

 

 

 

 

 

 

3

 

 

a2

 

a2

b2

b2

 

 

a2

 

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8D

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

2

b

2

 

b

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

3x2

 

 

 

 

y2

 

 

 

 

 

 

4

 

x2

 

 

 

 

 

3y2

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

+

 

 

 

 

1

 

 

 

 

 

 

 

 

+

 

 

 

2

1 .

 

 

 

(5.7)

 

3

 

 

 

 

 

 

2

 

 

 

 

 

 

 

3

 

 

a

2

 

 

a

 

b

2

b

2

 

 

a

2

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

+

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

a

2

b

2

 

b

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В центре пластинки, в точке (0, 0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

1

 

ν

M x = −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−ν

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

.

 

 

 

 

3

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

a

2

 

 

b

2

 

 

 

3

 

 

 

 

 

 

2

 

 

 

3

 

 

2

b

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

+

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

a

2

b

2

 

b

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

2

b

2

 

 

b

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.8)

Если пластинка круглая, то a = b = r. Тогда максимальный прогиб в центре круга равен

C =

qr 4

.

(5.9)

64D

 

 

 

Изгибающий момент в центре пластинки

M x =

qr 4

1

+

ν

=

(1+ν)qr 2

.

(5.10)

 

 

 

 

 

 

 

 

 

16

 

2

r

2

16

 

r

 

 

 

 

 

 

 

§2. Расчет прямоугольной плиты (пластины)

спомощью тригонометрических рядов

Получить решение в замкнутой форме, как в предыдущем параграфе, удается лишь в редких случаях. В частности, для прямоугольной пластинки получить такое решение не удается. Его приходится искать в виде бесконечного ряда.

103

Рассмотрим шарнирно опертую по всему контуру прямоугольную пластинку с равномерно распределенной по всей пластинке нагрузкой q (рис. 5.1).

Начало координат расположим в углу пластинки. Размер пластинки в направлении оси Ox равен a, а в направлении оси Oy – b.

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

0

 

 

 

 

 

 

x

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 5.2

Для шарнирно опертой по контуру пластинки должны выполняться следующие граничные условия:

w

w

 

 

 

2w

= 0;

 

 

x =0

=

x

2

 

 

 

 

 

x =a

 

 

 

x =0

 

 

 

 

 

 

 

 

x =a

 

 

 

 

 

2w

= 0.

(5.11)

 

y =0

=

y

2

 

 

 

y =b

 

 

 

y =0

 

 

 

 

 

 

 

 

y =b

 

 

Этим условиям удовлетворяют члены двойного тригонометрического ряда

∞ ∞

mπx

nπy

 

w(x, y) = ∑∑Amn sin

 

sin

 

.

(5.12)

a

b

m=1 n=1

 

 

 

 

 

Двойной ряд (5.12) представляет собой бесконечную сумму слагаемых, каждое из которых само является бесконечным рядом.

104

Для нахождения коэффициентов Amn, зависящих от порядковых номеров членов ряда, используем основное уравнение (4.31). Сначала вычислим производные:

 

w

 

 

∞ ∞

mπ

 

 

 

 

 

mπx

 

 

nπy

 

 

 

 

 

 

 

= ∑∑Amn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

cos

 

 

 

 

 

 

sin

 

 

 

 

 

;

 

 

 

 

 

 

 

 

a

 

 

 

 

a

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2w

 

 

∞ ∞

 

m2π2

 

mπx

 

nπy

 

 

x

2

 

= −∑∑Amn

 

 

 

 

2

 

 

sin

 

 

 

 

sin

 

 

 

 

 

 

 

 

;

 

 

 

 

 

a

 

 

 

a

 

 

b

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3w

 

 

∞ ∞

 

m3π3

 

 

 

 

mπx

 

nπy

 

 

x

3

 

= −∑∑Amn

 

 

 

 

3

 

 

cos

 

 

 

 

sin

 

 

 

 

 

 

 

 

;

 

 

 

 

a

 

 

 

a

 

 

 

b

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4w

 

 

∞ ∞

m4 π4

 

 

 

 

mπx

nπy

 

 

 

x

4

 

= ∑∑Amn

 

 

4

 

 

 

sin

 

 

 

sin

 

 

 

 

 

 

 

.

 

 

 

a

 

 

 

 

a

 

 

b

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4w

∞ ∞

 

n4 π4

 

 

 

mπx

nπy

 

 

y

4

= ∑∑Amn

 

 

 

4

 

 

sin

 

 

sin

 

 

 

 

 

 

.

 

 

 

b

 

 

 

a

 

 

b

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И наконец

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3w

 

 

∞ ∞

 

 

m2nπ3

 

 

mπx

 

 

 

 

nπy

 

 

 

 

 

= −∑∑Amn

 

 

 

 

 

 

 

 

 

sin

 

 

cos

 

 

 

 

 

 

 

;

2

y

 

 

 

 

2

b

 

 

 

 

b

 

x

 

m=1 n=1

 

 

 

a

 

 

 

a

 

 

 

 

 

 

 

4w

 

 

∞ ∞

 

m2n2π4

 

 

 

mπx

 

 

 

 

nπy

 

 

 

 

 

 

= ∑∑Amn

 

 

 

 

2

 

 

2

 

 

sin

 

 

 

sin

 

 

 

 

 

 

.

2

y

2

 

 

a

b

 

 

 

 

 

 

 

b

 

x

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

Подставим полученные выражения в основное уравнение (4.31):

 

 

∞ ∞

 

 

 

m4 π4

 

 

mπx

 

 

nπy

 

 

 

 

 

 

 

 

 

 

 

 

∑∑Amn

 

 

 

 

 

 

 

 

sin

 

 

sin

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

a

4

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

∞ ∞

 

 

 

m2n2π4

 

 

mπx

nπy

 

 

 

 

 

 

 

 

+2∑∑Amn

 

 

 

2

b

2

sin

a

sin

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1 n=1

 

 

 

a

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

∞ ∞

 

 

n4 π4

 

 

 

 

mπx

 

nπy

 

 

q

 

 

 

 

 

 

 

 

+∑∑Amn

 

 

 

 

 

 

 

sin

 

 

 

sin

 

 

 

=

 

 

 

 

 

.

 

 

 

 

 

 

 

 

b

4

 

 

 

a

 

 

 

 

D

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

После элементарных преобразований получим

 

 

 

 

 

 

 

4

∞ ∞

m4

 

 

m2n2

 

n4

 

mπx

 

nπy

 

q

 

π

 

∑∑Amn

 

 

4 +

 

2

 

 

2

 

 

2 +

 

4 sin

 

 

 

sin

 

 

 

 

 

=

 

(5.13)

 

a

 

a

b

b

a

 

b

D

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

105

или

 

4

∞ ∞

m2

n2 2

mπx

nπy

 

q

 

 

π

 

∑∑Amn

 

2 +

 

2

 

sin

 

sin

 

 

=

 

.

(5.14)

 

a

b

a

b

D

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

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

sin(mψ)

 

π

 

 

=

,

(5.15)

m

4

m=1

 

 

 

где индекс m принимает только нечетные значения: m =1, 3, 5, … Соотношение (5.15) выполняется для любых значений ψ в преде-

лах от 0 до 2π. В нашем случае ее можно переписать так

1

 

mπx

π

 

 

 

sin

 

 

=

4

 

 

 

 

 

 

m=1 m

 

 

a

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

1

 

mπx

 

1 =

 

 

 

sin

 

 

 

.

(5.16)

 

 

 

 

 

 

 

πm=1 m

 

 

a

 

 

 

Мы получили разложение единицы в ряд по направлению х. Точ-

но такое же разложение можно получить и в направлении у:

 

 

 

4

1

nπy

 

1

=

 

sin

 

.

(5.17)

 

 

 

 

π n=1 n

b

 

 

Перемножая (5.16) и (5.17), получаем разложение единицы в двой-

ной ряд:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

2

 

 

1

 

 

mπx

 

nπy

 

 

 

1 =

 

 

 

 

sin

 

 

sin

 

 

 

.

 

 

 

 

 

a

b

 

 

 

 

 

 

 

π

 

m=1 mn

 

 

 

 

 

 

 

 

Очевидно, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

q

4 2

 

1

 

mπx

 

nπy

 

 

 

 

=

 

 

 

 

 

 

sin

 

 

sin

 

 

 

.

(5.18)

 

D

 

 

 

 

 

 

 

b

 

 

D

π

m=1 mn

 

a

 

 

 

 

 

 

106

Здесь m и n принимают только нечетные значения: m = 1, 3, 5, …; n = 1, 3, 5, …

Подставим выражение (5.18) в уравнение (5.14)

 

 

 

 

 

4

 

∞ ∞

 

 

m2

 

n2 2

 

mπx

 

nπy

 

 

 

 

 

 

π

 

∑∑Amn

 

2 +

 

 

2

sin

 

sin

 

 

 

 

=

 

 

 

 

 

 

a

b

a

 

 

b

 

 

 

 

 

 

 

m=1 n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16q

∞ ∞

1

 

 

 

 

mπx

 

nπy

 

 

 

 

 

 

 

 

 

 

 

=

 

2

∑∑

 

 

 

 

 

sin

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

b

 

 

 

 

 

 

 

 

 

 

 

 

π

D m=1 n=1 mn

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∞ ∞

 

 

4

m2

n2 2

 

16q

 

 

 

 

1

 

 

mπx

 

nπy

 

 

∑∑

π

 

Amn

 

 

2 +

 

2

 

 

2

 

 

 

 

 

 

 

sin

 

sin

 

 

 

 

=0.

(5.19)

 

a

b

 

D

 

 

 

 

a

 

b

 

m=1 n=1

 

 

 

 

 

 

 

 

π

 

 

mn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

4

 

m2

 

 

n2

2

16q

 

1

= 0.

 

 

 

 

 

 

 

π

Amn

 

+

b2

 

π2D

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16q

 

 

 

 

 

 

 

16q

 

 

 

 

 

1

 

 

 

 

 

 

Amn =

 

 

 

 

 

 

 

 

 

=

6

 

 

 

 

 

 

 

 

 

 

 

.

(5.20)

 

m2

n2

 

2

D

 

m2

n2

 

2

 

 

 

 

 

π

 

 

 

 

 

 

 

 

π6Dmn

a

2 +

b

2

 

 

 

 

 

 

 

 

 

 

mn

a

2 +

b

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как в разложении (5.18) имеются только члены с нечетными значениями m и n, то все коэффициенты Amn с четными индексами равны нулю, и остаются только коэффициенты Amn с нечетными индексами. Во второй части равенства в скобку отдельно взят постоянный коэффициент, не зависящий от индексов. Он может быть вычислен один раз, и в дальнейшем это число будет просто умножаться на переменную часть формулы.

Возвращаясь к выражению (5.12) для функции прогибов w(x, y), подставим в него значения коэффициентов (5.20):

107

m=∞ n=∞
∑ ∑
m=∞ n=∞
∑ ∑

 

 

 

 

 

mπx

 

nπy

 

 

 

16q

m=∞ n=∞ sin

 

sin

 

 

 

 

 

a

b

 

 

∑ ∑

 

 

 

 

 

 

w(x, y) =

6

 

 

 

 

 

 

 

 

 

.

(5.21)

 

m2

n2

 

2

 

 

π

D

m=1

n=1

 

 

 

 

 

 

 

 

 

 

mn

2 +

b

2

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

Напомним, что здесь m =1,3,5..., n =1,3,5, ....

Имея выражение для w(x, y), можно найти моменты и поперечные силы. В частности, изгибающие моменты по формулам (4.17) и первой из (4.22):

=164qπ

=16π4q

m=1

m=1

 

 

 

 

 

2w

2w

=

 

 

 

 

M x = −D

x 2

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m2

 

 

 

n2

 

 

 

 

 

 

 

 

 

 

 

a2 b2

 

 

 

mπx

nπy

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

sin

 

.

 

m2

 

 

n2

 

2

 

a

b

n=1

+

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

2

 

 

b

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2w

2w

=

 

 

 

 

M y = −D

y2

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m2

 

n2

 

 

 

 

 

 

 

 

 

 

 

ν a2 + b2

 

 

 

mπx

nπy

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

sin

 

.

 

m2

 

 

n2

 

2

 

a

b

n=1

+

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

2

 

 

b

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Mx

в середине квадратной плиты со стороной а (а = b), т. е. в точке

с координатами x = 0,5a; y = 0,5a. Коэффициент Пуассона ν = 0,25.

 

В этом случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16qa2

m=∞ n=∞

 

m2 +0,25n2

 

 

 

π

 

π

 

M x =

 

 

4

 

∑ ∑

 

 

 

sin m

sin n

.

 

π

mn(m2 +n2 )

2

 

 

 

 

m=1 n=1

 

 

 

2

 

2

 

Обозначим для упрощения записи

 

 

 

 

 

 

 

 

 

 

m2

+0,25n2

 

 

 

 

π

 

π

 

 

R =

 

 

 

 

 

 

 

; S = sin m

sin n

,

 

 

mn(m

2

2

)

2

 

 

 

 

 

+n

 

 

 

 

2

 

2

 

108

так что

M x = 16qa4

2

 

∑ ∑R S .

 

 

 

m=∞ n=∞

 

π

 

 

m=1 n=1

Примем сначала m = 1 и найдем сумму, давая значения n = 1, 3, 5, …

Тогда коэффициенты R будут последовательно равны:

R11 = 12 +0,25 12 = 0,3125;

1 1(12 +12 )2

R13 = 12 +0,25 32 = 0,0108; 1 3(12 +32 )2

. . . . . . . . . . . . . . . . . . . . . .

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

Затем примем m = 3 и вновь проделаем то же самое и т. д. Будем продолжать эту процедуру до тех пор, пока первый же коэффициент Rm1 не будет иметь четыре нуля после запятой. На этом вычисления остановим и просуммируем полученные суммы.

Запись будем вести в табличной форме.

 

m =1;

 

π

= +1

m =3;

 

 

π

m =5;

 

 

π

 

sin m

2

 

sin m

2

= −1

sin m

2

= +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

R

 

π

RS

R

 

 

π

 

RS

R

 

 

π

 

RS

sin n

2

 

 

sin n

2

 

 

sin n

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,3125

+1

 

 

+0,3125

0,0308

 

+1

 

 

–0,0308

0,0075

 

+1

 

 

+0,0075

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

0,0108

–1

 

 

–0,0108

0,0039

 

–1

 

 

+0,0039

0,0016

 

–1

 

 

–0,0016

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

1,0021

+1

 

 

+1,0021

0,0009

 

+1

 

 

–0,0009

0,0005

 

+1

 

 

+0,0005

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

0,0008

–1

 

 

–0,0008

0,0003

 

–1

 

 

+0,0003

0,0002

 

–1

 

 

–0,0002

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

0,0004

+1

 

 

+0,0004

0,0001

 

+1

 

 

–0,0001

0,0001

 

+1

 

 

+0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

0,0002

–1

 

 

–0,0002

0,0001

 

–1

 

 

+0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

109

13

0,0001

+1

 

 

 

 

 

+0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

0,0001

–1

 

 

 

 

 

–0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

17

0,0001

+1

 

 

 

 

 

+0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

+0,3033

 

 

 

 

 

 

 

–0,0275

 

 

 

 

 

+0,0063

 

m = 7;

 

 

π

 

 

m = 9;

 

 

 

π

m =11;

 

 

π

 

sin m

 

 

= −1

 

 

sin m

 

= +1

sin m

2

= −1

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

R

 

π

 

 

RS

 

R

 

 

 

π

RS

R

 

 

π

 

 

 

sin n

2

 

 

 

 

 

sin n

2

 

 

sin n

2

 

 

RS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,0028

+1

 

 

 

 

 

–0,0028

 

0,0013

 

 

+1

 

 

+0,0013

0,0007

 

+1

 

 

 

–0,0007

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

0,0007

–1

 

 

 

 

 

+0,0007

 

0,0004

 

 

–1

 

 

–0,0004

0,0002

 

–1

 

 

 

+0,0002

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

0,0003

+1

 

 

 

 

 

–0,0003

 

0,0002

 

 

+1

 

 

+0,0002

0,0001

 

+1

 

 

 

–0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

0,0001

–1

 

 

 

 

 

+0,0001

 

0,0001

 

 

–1

 

 

–0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

–0,0023

 

 

 

 

 

 

+0,0010

 

 

 

 

 

 

–0,0006

 

m =13;

 

 

π

 

= +1

 

m =15;

 

 

π

 

 

 

 

π

 

sin m

2

 

 

sin m

= −1

m =17; sin m

 

= +1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

π

 

 

 

 

 

 

 

 

π

 

 

 

 

π

 

 

 

n

R

sin n

2

 

 

 

RS

 

R

 

sin n

2

 

RS

R

 

sin n

2

 

 

RS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,0005

+1

 

 

 

 

 

+0,0005

 

0,0003

 

 

+1

 

 

–0,0003

0,0002

 

+1

 

 

 

+0,0002

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

0,0001

–1

 

 

 

 

 

–0,0001

 

0,0001

 

 

–1

 

 

+0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

+0,0004

 

 

 

 

 

 

 

–0,0002

 

 

 

 

 

 

+0,0002

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m =19;

 

 

π

 

= −1

 

m = 21;

 

 

π

 

 

 

 

 

 

 

 

sin m

2

 

 

sin m

= +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

n

R

sin n

2

 

 

 

RS

 

R

 

sin n

2

 

RS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,0002

+1

 

 

 

 

 

–0,0002

 

0,0001

 

 

+1

 

 

+0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

–0,0002

 

 

 

 

 

 

+0,0001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Просуммируем все результаты:

∑∑= +0,2805.

110

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