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

книги из ГПНТБ / Киселев В.А. Расчет пластин

.pdf
Скачиваний:
41
Добавлен:
23.10.2023
Размер:
4.3 Mб
Скачать

Вместо функций Рг

(у), F2

(у),

F3

(у) и FA

(у) основного случая

•в выражениях

(3.82) — (3.85)

будут

функции:

 

 

/•3COSrj(/

.

Г l cos г3

у

 

 

 

— 2

—2

 

 

— 2

—2

 

 

 

 

гі—лз

 

 

гі— гз

 

 

о

/ , л _

r i s i n g

 

 

,

nsin/-8 i/

 

 

 

^(гі —лз)

 

 

г3 (п—гз)

 

S3

(У) = * - 2

(cos /-! у cos г3

у);

 

 

п — гъ

 

 

 

 

 

 

s4

(у) =

sin гх у

 

 

,

sin г3

у

 

- /-а - 2 \

 

-

/ - 2

- 2 \

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

5-й случай

с с т < 0 ;

Я,*, >

0;

А,,"„ = о с ш

> 0

 

Корни характеристического

уравнения:

 

гs = —/-4 = V~ am—i YKi — a?„ = ßm ym i.

Этот случай сводится к первому — см. (2.108) — (2.109).

6-й случай а 2 , > 0 ; \ т

—ат>0

Такой случай возможен для изотропной пластины без упругого основания, когда Нх = Ну = 0. Заменяющие функции были опре­ делены ранее, см. (2.176) — (2.179).

Все полученные функции для различных корней характеристи­ ческого уравнения, удовлетворяющие единичной матрице и имею­ щие производные по табл. 4, входящие в выражения (3.83) — {3.95), могут быть обозначены единообразно через F± (у), F2 (у),

F3 (у) и FA (у) и сведены в общую таблицу.

Если на пластину действует вибрационная нагрузка

q (х, у, t) = q (х, у) (A sin Qt + В cos Ѳ^),

то к дифференциальному уравнению (3.56) в левую его часть надо

добавить инерционные силы т*

где т* — масса пластины и

приближенно присоединяемого основания на единицу площади,

111

учитывающей влияние основания на колебания пластин. Диф­ ференциальное уравнение без учета сил сопротивления будет

1 од:4

3

дх2 ду2

2

ai/4

дх2

ѵ

ду*

 

4- с. ш) — Со

d-w , d2w \ ,

3 2 ш

 

 

 

1

Ч- m

 

 

 

1

- ^дл:2

ду* )

 

dt2

 

 

=

—q(x,

y){AsmQt

+ BcosBt).

(3.99>

Здесь w (x, y, i) — прогиб пластины, изменяющийся во времени. Для установившихся колебаний положим

 

w (x, у,

f) =

w (х, у) (A sin

Qi +

В cos

Ѳ/).

(3.100)

Подставляя

(3.100)

в (3.99),

получим

 

 

 

 

1

дх*

3 дх2

diß

2 ду*

дх"-

J

diß

 

 

+ c \

w - c J ^

+ ^ \ =

-q(x,y),

 

(3.101)

 

 

 

дх2

dtf-

 

 

 

 

где

 

 

c î = d — / я * Ѳ .

 

 

 

(3.102)

 

 

 

 

 

 

Это дифференциальное уравнение по форме совпадает с исход­ ным дифференциальным уравнением (3.56). Поэтому все ранее полученные результаты на основе уравнения (3.56) и все ранее рассмотренные случаи загружения в примерах изотропных пластин,, при соответствующей небольшой их обработке, применимы и для вибрационной нагрузки с заменой во всех выражениях сх на новое значение с* = сг — /п*Ѳ2 .

В заключение приведем краткие соображения о свободных коле­ баниях и устойчивости пластин. Частоты собственных колебанийпластин со, а равно и критические нагрузки Нх и Ну могут быть найдены из определителя уравнений, написанных для вычисления неизвестных начальных параметров, при отсутствии поперечных нагрузок, если его (определитель) приравнять нулю. Тогда частоту &

действующей

вибрационной

нагрузки

следует считать неизвест­

ной частотой

собственных колебаний

со или нагрузки

НхиНу

критическими,

подлежащими

определению. При этом

множители

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

112

§ 21. Расчет прямоугольных пластин, по-разному опертых по всем четырем сторонам

Рассматриваются пластины, имеющие шарнирные опирания и заделки. Метод решения проследим на пластине с защем-

.ленной стороной х = 0 (рис. 54, а). Решение проводим по методу сил. За основную систему принимаем пластину, шарнирно опертую по всем сторонам (рис. 54, б). Неизвестные,опорные моменты по краю X — 0 раскладываем в ряд

 

с о

 

 

Мх(у) =

2

CÄsinüM .

(3-103)

-

п = 1

b

 

Каноническое уравнение перемещений по освобожденному от за­

щемления краю

пластины

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ( ^ - )

 

= 0 .

.(3.104)

 

 

дх /от Mxj

 

\

дх

/от q

 

 

 

Используем готовые решения для пластины, шарнирно опертой

по четырем сторонам в двойных

 

рядах

(см. 3.57) — (3.61):

 

 

 

о о

с о

 

 

 

 

тлх

 

.

пли

 

/

\

V I

V I

 

 

 

 

 

t

v i

v

 

 

 

 

sin

Ь*

 

w(x,

у =

. у

у

 

 

a„,„sm

 

а

 

 

(х,

у)=

=

1 Zn =

 

1

 

~^~

 

 

 

 

 

m- 2

 

s m

s

 

 

 

 

 

 

 

a

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

_

o o

с о

 

 

 

 

 

 

 

 

 

 

 

aw

v i

v

 

mn

тлх

 

пли

 

— =

У

>

атп

a

cos

a

sin — J ~

 

дх

m=1n=1

 

 

 

 

 

 

b

При X = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с о

с о

 

 

 

 

 

 

 

 

 

i r =

2 2

a

™ ~ s

i n

4 -

( З Л 0 5 )

Здесь

 

m =

I n = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fl»»

= -

ï ™ r ' L

по

(3.60),

где определяют: c7mn по (3.59),

 

ЬПп

по (3.61).

 

От неизвестных моментов по (2.83) имеем

 

 

 

 

атп

= 2птС%:Ьтпа\

 

 

 

(3.106)

Подставляя (3.105) в уравнение (3.104) с учетом (3.106), получим

с о

с о

 

 

 

 

 

 

 

 

 

 

 

 

 

- — • —— sin —— +

 

У

>

am 7 l (от q) — sin —— = 0.

т=

i п=1

 

 

 

m = 1

п = 1

 

 

 

 

 

113

'y)

s)

Рис. 54

Рис. 55

Умножим это равенство на sin ^~

dy (где п конкретное, а не

текущее) и проинтегрируем от нуля до Ь:

СО

с о

 

/гад \ . «я« ,

J

2 ,

1 n =

2

^

-

г

8 Ш

— — 1 sin — — ау-

 

0

\ m =

1

 

 

 

 

 

 

 

 

^

I

с о

о о

 

 

 

 

 

,

 

 

+

 

2

2

a »»<0 T »>

sin

- ^ j s i n Ä

А/ = 0.

0 \ m = l

n = l

 

 

 

 

 

 

 

 

Учитывая

свойство

ортогональности, см. (2.82),

будем иметь

 

 

 

2я/пСп

mit

6 .

 

 

тп

b

п

 

от=

1

 

 

 

 

 

m = 1

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tnamn

(от с) ï

m a

 

(3.107}

 

 

 

 

 

m = 1

 

 

 

m = 1

 

 

В частном

случае

от

равномерной

нагрузки

было получено

(см. стр. 31).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

4<7

(1—cos тп) (1—cos пп)

 

 

 

 

 

 

я 2

тп

 

 

Ьтп

 

 

 

114

Подставим в (3.107)

Сп* '•

 

4q (I—cosmst)

( 1 — cos пп)

, 2п

 

 

 

 

т= 1

 

П2

nbmn

 

а2

2А

 

 

 

 

 

 

 

 

 

а

т~1

и

т п

 

 

 

 

 

 

 

 

 

 

После преобразований

 

 

 

 

 

 

 

 

 

 

 

2qa2

(1 — cos пп)

^

[U — cos т я ) : 6 m n ]

 

 

 

г*

_

 

 

 

 

 

 

 

 

 

 

 

 

 

ПП3

2

№-Ьтп)

 

 

 

 

 

 

 

 

 

т= 1

 

 

 

 

 

 

Учитывая,

далее, что при m = 2, 4, ... и

п. = 2,

4,

... (1 —

— cos mn) =

0 и (1 — cos пп)

=

0, окончательно получим,

С* —

Sqa2

о о

 

j

 

ОО

 

 

 

(3.108)

 

2

 

&„,„

2

 

 

 

 

 

 

о,

о

 

 

 

 

 

 

. m = 1,

о

ш =

1, о,

 

 

 

После определения

С,* опорные

моменты получим по (3.103).

Аналогично проводится расчет и на другие нагрузки.

 

 

Если пластина имеет несколько заделок, то надо составить

столько же канонических уравнений,

которые в этом случае будут

сложнее уравнения (3.104) за счет

взаимного

влияния

опорных

моментов на углы поворота. Естественно, что расчет таких

пластин

значительно сложнее рассмотренной. Можно расчет проводить и последовательно, приняв, например, пластину с одним защемлением в качестве основной системы для расчета пластины с двумя защем­ лениями, и применить далее, что было изложено здесь.

§ 22. Понятие о расчете пластин,

лежащих на упругом полупространстве

Просадка полупространства от сосредоточенной силы Р определяется по формуле (1.17)) (рис. 55, а):

к,* = р ; Е*Г,

(3.109)

где Е* = пЕ0 : (1 [A2,); ßQ и ^ о — соответственно

модуль упру­

гости и коэффициент Пуассона основания.

 

 

Обозначим, как и ранее, через г (х, у) давление от пластины

на

основание в точке k с координатами х и у.

 

 

Рассмотрим некоторую точку основания с координатами х =

\ и

у = т] в пределах пластины (рис. 55, б). Бесконечно малая сила на площадке d£,dr\ в этой точке будет

dP = г (g, r\)dldr\.

115

Бесконечно малое перемещение от нее в точке k основания с ко­ ординатами X и у

dw* (х,у) =

' в '

.

 

£ * V ( * - ! ) 2 + ù / - ' n ) 2

 

Полное перемещение точки k основания от всей неизвестной на­ грузки г(х, у)

{ж, у) =

і

.

(3.110)

J

J

JS*V

(у-т,)»

Это перемещение основания должно быть равно прогибу пла­ стины, определяемому из дифференциального уравнения

D l

+

2D3-*ÜL

+

 

Dn-^_Я;сÜÜL_

 

1

ox*

3 А-3ОГ/2

 

ду*

x

дх2

 

 

-НУ-^~-=~ІЯ(Х,

 

У) +

Г ( Х ,

y)].

(3.111)

В это уравнение также входит неизвестная реакция упругого основания.

Таким образом, расчет пластины, лежащей на упругом полу­ пространстве, сводится к совместному решению уравнений (3.110) и (3.111). Решение этих уравнений представляет значительные ма­ тематические трудности. С подробностями расчета читатель может ознакомиться в монографии М. И. Горбунова-Посадова [16].

Г л а в а 4

ПРИБЛИЖЕННЫЕ МЕТОДЫ РАСЧЕТА ПРЯМОУГОЛЬНЫХ ОРТОТРОПНЫХ И ИЗОТРОПНЫХ ПЛАСТИН, ЛЕЖАЩИХ НА УПРУГОМ ОСНОВАНИИ

СДВУМЯ КОЭФФИЦИЕНТАМИ ПОСТЕЛИ

§23. Метод Бубнова — Галеркина

Дифференциальное уравнение изгиба ортотропной пла­ стины на упругом основании с двумя коэффициентами постели в об­ щем случае при вибрационной или статической нагрузке запишем по '(3.101) в таком виде:

ох

-Ну~- ду2-

+ c\w-c0y*w

 

+ q(x, у)=0,

 

(4.1)

где

 

 

 

 

 

 

 

 

I = D

, — + 2 D 8

— + A» —

(4.2)

1

дх*

 

дх*ду°-

2

ду1

 

ѵ

'

линейный оператор.

116

Будем искать приближенное решение уравнения (4.1) в виде

ряда функций фг (х, у) с неопределенными

коэффициентами at.

п

 

У ) 4 » 2 ЩЧі(х, у)-

(4.3)

І= \

Произвольные функции этого ряда должны удовлетворять гео­ метрическим и статическим граничным условиям пластины, а коэф­ фициенты at подбираются по некоторым соображениям хорошего его приближения к действительно упругой поверхности до (х, у). Кроме того, функции фг (х, у) должны быть «похожими» на ожидае­ мую упругую поверхность пластины. Так, например, в симметрич­ ных случаях относительно двух или одной оси симметрии эти функции также должны быть соответственно симметричными. Иначе включаемая в общее выражение произвольной функции кососимметричная часть будет только искажать симметричную картину пере­ мещений, если соответствующие коэффициенты аг отличны от нуля.

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

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

Подставляя приближенное решение (4.3) в дифференциальное уравнение (4.1), получим левую часть ее, отличную от нуля, которая- и определяет ошибку решения

лп.

-fei 2 в | ф | — c 2

V a 2

0гФг +

<7(*. У)фО.

(4.4)-

J = I

;=ï

 

 

 

По методу Бубнова — Галеркина так

подберем

коэффициенты'

а, в (4.3), чтобы эта ошибка

(по условию минимума

ошибки в сред­

нем) была ортогональна к любой

функции фй (х, у)

из этого

ряда.

Для записи ортогональности

умножим (4.4) на функцию срА (х, у)

и составим интеграл по всей площади пластины, а результат при­ равняем нулю:

1= I

/=

I

<= !

 

л

 

п

 

 

 

+ с:\ 2

a, q^—c,

2 аі

Ѵа Фг+<7(*,

y)\<phdxdy=zO,

(4.5)

где k = 1, 2,

п.

 

 

 

 

117

Представим уравнение (4.5) в развернутом виде, учитывая, что оператор L линейный:

 

ab

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£

Ф і - Л я

-

^

-

 

 

-

^ -

+

d\ Ф і

- с ,

Ѵ Ф

і ] фА dxdy

+

 

 

 

 

 

 

 

ѵ

ду2

 

 

 

 

 

 

 

о о

 

а Ь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ff* ф 2

 

2

ф 2

 

 

 

 

 

+ "'Я

 

 

 

 

 

 

 

 

 

 

 

 

 

d*2

 

 

 

 

 

 

 

 

 

 

 

О О

 

 

 

 

Я

 

 

 

 

 

 

 

 

 

Х ф ь ^ у +

.-. +

а,,

 

 

 

дх-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n о

 

 

a

b

 

 

 

 

+ С ! Ф Й

- С

2 Ѵ2 ФЙ

]<pkdxdy+...+an

 

Jj*[Дф - Н х д - ^ - - Н у

^ +

 

 

 

 

 

 

 

 

 

 

 

0 0

 

 

 

"

 

4- с*Ф„—с2 Ѵ2Ф„1 Ф й с Ы # 4 - ^q(x,

y)yhdxdy=0.

'(4.6)

 

 

 

 

 

 

 

 

 

 

 

о 0

 

 

 

 

Это же

уравнение

в

канонической

форме

 

 

 

 

 

 

 

ai + ch2

a2

+ ...+chn

cn+Chp

= 0

 

(4.7)

при

k — 1,

2, ...

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь по

(4.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a ь

L<pt-H\

^

-

 

-

H

v - ^ Ï L +

C f Ф г

- с 2

Va Ф

| ] Ф й

(4.8)

- =

яо о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

b

 

 

 

 

 

 

(4.9)

 

 

 

 

 

Chp

= \\a{^

 

y)q>hdxdy.

 

 

 

 

 

 

 

 

 

 

о 0

 

 

 

 

 

 

 

Покажем, что если функции фг в (4.3) удовлетворяют всем гра­ ничным условиям пластины, то коэффициенты канонических урав­ нений (4.7) обладают взаимностью, т. е. см = сі1{. Для этого рас­ смотрим две функции фг и Ф л , удовлетворяющие, при соответст­ вующих им нагрузках qt и qh, дифференциальному уравнению (4.1):

L

^ - H

x ^ — H v ^

+

с 1 Ф г _ С 2 Ѵ 2 Ф г

= -<?г (*,У);

(4.10)

bth-Hx

^

Ж - Н ^ + с

^ ъ

- с , Ѵ а Ф , =

- д й ( х , у).

(4.11)

Теперь по (4.8) а учетом (4.10), (4.11) может записать:

и b

 

 

сы = \\ ЯІ(Х,

y)yhdxdy\

(4.12)

a b

 

 

et* = \ \ qh(x,

y)fptdxdy.

(4.13)

Выражение (4.12) можно толковать как работу нагрузки qt на перемещениях ц>к, а выражение (4.13), наоборот, как работу нагруз­ ки qh на перемещениях ср,-. Эти работы по теореме взаимности работ должны быть равны, а потому

chi = Cih-

(4.14)

Если в частном случае выражение (4.3) дает точно упругую поверхность пластины w (х, у), то (4.4) не будет содержать ошибки и канонические уравнения (4.7) дадут нам точные значения коэф­ фициентов at. Покажем это на примере шарнирно опертой по всем четырем сторонам пластины. Здесь (4.3) точно:

*»(х> У)= 2

2

a m n S i n - ^ s i n Ä .

(4.15)

m =

1 п =

I

 

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

а„ =

amn;

cpft =

sin

тлх

.

пли

sin —f- ;

 

 

 

 

a

 

о

а г =

ат »п»;

г =

.

m* лх

.

п*лу

5іп

 

sin

По формуле (4.8)

о о L

.

пг пх

 

 

тлх , плд j j

sin

 

sin

b J

sin—?-dxdy.

\ а 2

 

 

a

 

 

 

 

В

силу

 

ортогональности

тригонометрических

 

функций

ви = 0.

По

той

 

же

формуле

(4.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, „

/

о

о

, „

/

ц п \ 2 ,

 

,

. ' т

 

з т зт

,n

nяjt

\")\T -

mnx

 

 

m i t \ 2

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

2

2

 

2

2

2

Sin

 

м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X sm —— \ sin

 

sin —— dxdy.

 

 

 

 

 

 

 

 

 

 

 

b

)

 

a

 

 

 

b

 

 

 

 

 

 

 

 

 

После

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ab

Гг-. /

mn

, О Г

1

/• mn

\ 2

 

/ л л \ 2

 

 

 

/

rm\*

,

 

 

с

» * - т И ~ )

+

2

М ~ ) ( T ) + D ' ( T ) +

 

-

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

,4.16)

Канонические

уравнения

(4.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при

 

k = 1,

2, ...

 

 

 

 

 

 

 

 

Отсюда

получаем

 

 

 

 

a fr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mux

 

 

 

nit;/

 

 

 

 

 

 

 

 

 

 

 

 

 

г f

q (x,

у)

sin

sin

 

 

 

 

 

 

 

 

 

 

п.

 

 

— \ \

 

 

— — dxdy

 

 

 

oft

= a m n =

-

=

J J

 

 

 

 

 

a

 

b

 

(4.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это выражение совпадает

G ранее полученным значением коэф­

фициента

am 7 i

(3.59) — (3.61),

и

решение

(4.15)

будет точным.

 

Если же в (4.15) ограничиться конечным числом членов, то и при

точных

коэффициентах атп

решение будет уже не точным. Иногда,

вообще, в качестве первого приближения задают всего одну функ­

цию

ф! (х,

у) и тогда имеют

одно

каноническое уравнение

(4.7):

 

 

с п О і

+ С

= 0.

(4.18)

В

этом

случае особенно важно,

чтобы функция

cpj (х, у)

была

как можно больше «похожа» на ожидаемую упругую поверхность пластины w (х, у ) , иначе результаты даже для первого приближения могут оказаться ненадежными.

Все приведенные ранее в этом параграфе выражения могут быть применены и для расчета изотропных пластин, если в них положить

Eh3

Di = D 2 = Ds = 1 2 ( 1 «) и оператор L = ѵ а Ѵа-

Покажем применение уравнения (4.18) на простом примере изот­ ропной шарнирно опертой по четырем сторонам пластины о равно-

'120

Соседние файлы в папке книги из ГПНТБ