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

b

.pdf
Скачиваний:
22
Добавлен:
17.03.2016
Размер:
1.59 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

 

 

2 w

 

m

 

2 w

l

 

 

 

δ w

d s + ∫∫

2 w

 

 

2δ w

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

y

2

 

 

 

x

 

y

2

 

 

x

2 dA ,

 

 

 

 

 

 

 

 

 

 

 

 

S xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 w

 

2δw

 

 

 

 

 

 

 

2 w

 

 

δ w

 

 

 

 

 

 

 

 

 

 

 

 

 

δ w

 

 

3w

 

 

 

 

 

 

 

 

∫∫

xy

xy dA = xy

 

y

 

l d s ∫∫

 

 

y

 

 

 

∂ ∂

 

2 dA =

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

2 w

 

l

2 w

 

 

 

 

δ w

d s + ∫∫

2 w

 

 

2δ w

dA .

 

 

 

 

 

 

 

 

 

 

 

 

xy

x

2

m

 

y

 

x

2

 

 

y

2

 

Тогда

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 w

 

2δ w

 

 

 

2 w

 

2δ w

 

 

2 w

 

 

2δ w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∫∫

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

dA =

 

 

 

 

 

 

 

 

 

 

 

 

xy

xy

 

x

 

 

y

 

 

 

y

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

w

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

w

 

 

 

 

 

w

 

 

=

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

x

l

 

y

x

 

 

δ

y

+

 

y

m

 

x

 

y

 

 

δ

 

x ds

=

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= S

 

2 w

 

 

 

 

w

 

 

 

2 w

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ x

δ

 

y

τ y δ

 

x ds =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

w

 

 

 

 

w

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

w

 

 

 

 

 

 

 

=

 

 

 

 

δ

 

 

 

 

+l

 

 

 

 

 

 

 

δ

 

 

 

 

m

 

 

 

 

 

 

 

τx

m

n

τ

 

 

 

τy

l

 

n

 

 

τ

 

ds =

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

w

 

 

 

 

 

w

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

w

 

 

 

 

 

w

 

 

 

 

=

 

 

 

 

 

l

 

δ

 

+

 

 

 

 

 

+ m

 

 

δ

 

ds =

 

τ

 

m

x

y

 

 

n

 

τ

l

x

 

y

 

 

τ

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 w

 

 

 

 

 

w

 

 

 

2 w

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

(6.2.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

nτ

 

δ

 

 

τ

 

τ

2

 

δ

 

n

 

ds .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь вариацию потенциальной энергии получим в виде

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

w

 

 

 

2

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w δ wdA +

 

 

w δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δU = D ∫∫

 

 

 

 

 

n

 

 

 

 

n

 

 

 

δw ds +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+(1ν )

 

 

2

w δ w

2

w δ

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ds =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nτ

 

 

 

 

 

τ

 

 

τ

2

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

w

 

 

 

 

 

2

w

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

δ w dA +

 

 

 

 

+ν

 

 

 

 

δ

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= D ∫∫

 

 

 

n

2

 

 

τ

2

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

w

 

 

 

 

 

w

 

 

 

2

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

(1 ν )

 

 

 

 

δ

 

δw ds

 

.

 

 

 

(6.2.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

nτ

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заменяя в контурных интегралах функции прогиба на внутренние усилия в соответствии с формулами (6.1.8), получим

 

4

 

w

 

 

w

δU = D∫∫

 

 

 

 

M nτ δ

 

 

 

w δ w dA + Qn δw M n δ

 

n

τ

ds

A

 

S

 

 

 

 

Вариация работы внешних сил определится формулой

 

 

 

 

 

 

w

 

w

 

 

 

 

δ

+ mnτ δ

 

δT = ∫∫q(x,y) δwdA + p(x, y) δw + mn

n

 

τ

d s

A

 

Sk

 

 

 

 

.

(6.2.5)

.

(6.2.6)

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

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

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

δЭ = ∫∫[D 4 w q(x, y)]δwdA +

 

A

 

 

 

 

 

 

 

+

(Qn p) δwd s (M n + mn ) δ

w d s

 

Sk

Sk

 

 

 

 

n

 

 

 

 

 

 

 

 

 

(M

τ + m

nτ

) δ w d s = 0 .

(6.2.7)

 

 

n

 

τ

 

 

 

Sk

 

 

 

 

 

 

 

 

 

 

 

 

 

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

4 w(x, y)

q(x, y)

= 0 ;

(6.2.8)

D

 

 

 

Qn – p = 0; Mn – mn = 0;

Mnτ – mnτ = 0,

(6.2.9)

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

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

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

6.3. Метод Ритца Тимошенко в задачах изгиба пластинок.

Для удобства выкладок и проведения вычислений перейдем к безразмерным координатам:

 

ξ =

x

;

η =

y

,

λ =

a

,

(6.3.1)

где а, b

 

a

 

b

 

b

 

- характерные размеры пластинки, в частности можно

принимать

а = b.

 

 

 

 

 

 

 

Для прямоугольной пластинки а, b -

размеры пластинки. Область

изменения координат

0 х а,

0 у b переходит для

безразмерных координат в фиксированную область

0 ξ 1,

0 η 1.

Дифференциальные операторы заменяются согласно формулам:

..

=

..

 

ξ

=

 

1 ..

;

 

 

 

..

=

..

 

η =

1

 

..

;

 

 

 

 

ξ

 

 

 

 

 

 

 

 

y

η

 

 

 

 

 

 

 

x

 

 

x a ξ

 

 

 

 

 

 

y b η

 

 

 

 

k ..

=

1 k ..

;

 

k ..

=

1

 

k ..

;

 

 

 

k ..

=

 

1

 

k ..

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xk

ak ξ k

 

yk

bk

 

ηk

 

xm y k m

ambk m

 

ξ m ηk m

 

 

 

 

 

 

 

 

 

 

 

 

dx = adξ ;

 

 

dy = bdη ;

 

2

 

2 ..

 

 

2

..

 

1

 

 

.. =

 

 

+

 

 

 

 

=

 

 

 

 

2

 

 

 

 

2

 

 

 

x

 

 

x

2

 

a

 

 

 

 

 

 

 

 

 

 

(6.3.2)

 

 

 

 

 

 

 

 

 

 

 

 

~

 

dA = dxdy = ab dξ dη = ab dA .

 

2 ..

2 2 ..

 

1 ~ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

2 +λ

 

η

2

 

=

a

2 ..

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

2

 

2

 

2

w

2

 

2

w

2

w

 

 

 

 

 

 

D ~ 2

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U =

 

2λa

2 ∫∫ ( w)

+ 2(1 ν )λ

 

 

 

 

 

 

ξ

2

η

2

 

dA ; (6.3.3)

 

 

 

A

 

 

 

 

ξη

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

~

 

w(ξPi

 

 

 

)+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T =

 

 

 

∫∫q(ξ,η) w(ξ,η)dA + Pi

,ηPi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mx (ξ,η(ξ))w

 

 

 

 

my (ξ,η(ξ))w

~

+ ap(ξ, y(η)) w(ξ,η(ξ))+

 

 

 

 

 

 

 

 

+ λ

 

 

 

 

 

 

 

 

ds .

 

a

 

 

 

ξ

 

 

 

 

a

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.3.4)

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

∞ ∞

 

w(ξ,η)= G ∑∑Amn wmn (ξ,η),

(6.3.5)

m=1n=1

где Аmn - неопределенные коэффициенты (определяются из условия минимума полной энергии деформаций); wmn (ξ,η) -

функции, удовлетворяющие кинематическим граничным условиям опирания пластинки (статические граничные условия могут не выполняться); G – коэффициент, который принимается в процессе решения, для приведения решения к наиболее удобному виду.

При подстановке решения в функционал полной энергии деформаций, функционал становится функцией неопределенных коэффициентов Аmn, которая достигает минимума, если частные производные по всем аргументам равны нулю, т.е.

Э

=

U

T

= 0 или

U

=

T

, k ,l =1,2,3,...

(6.3.6)

A

A

A

A

A

 

 

 

 

 

 

kl

 

kl

 

kl

 

kl

 

kl

 

 

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

 

 

Bklmn Amn = Ckl ,

k ,l =1,2,3,...

(6.3.7)

m=1

Чтобы получить коэффициенты Bklmn и Сkl проведем

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

U

 

 

D

 

 

 

 

 

 

 

~ 2

 

2

 

2

 

2

w

2

2

w

2

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

( w)

+ 2(1 ν )λ

 

 

 

 

 

 

 

 

 

 

 

 

 

dA =

 

=

 

 

2 ∫∫

 

 

 

 

 

 

 

2

 

 

2

Akl

2λa

 

 

 

 

 

 

 

ξ

η

 

 

 

 

A Akl

 

 

 

 

 

 

ξη

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

~

2

 

 

~

2

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

∫∫ (

 

w)

A

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λa2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ (1ν )λ

2

 

 

 

 

2 w

 

 

 

 

 

∂ ∂2 w

 

 

 

2 w

 

 

 

 

 

2 w

 

2 w

 

 

2 w

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η2

 

η2

 

 

 

 

 

 

ξ 2

dA .

 

ξη

 

 

A

 

 

ξη

ξ

 

 

 

A

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Или, подставляя решение (6.3.5), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

=

 

 

D

 

 

G2

 

 

 

 

 

[

2 w 2 w +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

λa2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1∫∫

 

 

 

 

mn

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 w

 

2 w

 

 

 

2 w

 

 

2 w

 

 

2 w

 

 

 

2 w

 

 

~

;

 

 

 

 

 

 

+ (1 ν )

ξη

ξη

 

 

 

ξ

 

2

 

 

 

 

η2

 

η2

 

 

ξ 2

dA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

kl

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

kl

 

 

 

 

mn

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

P

 

 

 

w(ξPi

,ηPi )+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

G∫∫q(ξ,η)w(ξ

,η)dA

 

+

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

λ

 

 

 

ab

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p(ξ

, y(η))

w(ξ,η(ξ))+

 

mx (ξ

,η(ξ))w

 

 

 

my (ξ

,η(ξ))w

 

~

 

 

 

 

+λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+λ

 

 

 

 

 

 

 

 

 

 

 

 

ds .

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

a

2

 

 

 

 

 

 

 

ξ

 

 

 

a

2

 

 

 

 

η

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приравнивая производные потенциальной энергии и работы

внешних сил, сокращая общие множители и положив

 

 

G =

 

q

0

a4

 

,

 

 

 

 

 

D

 

 

 

 

 

q0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

некоторое значение

 

 

 

 

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

 

произвольной

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

~

2

 

 

 

~

2

wkl +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bkl

= ∫∫[

 

 

wmn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 w

 

 

 

2 w

 

2 w

 

 

 

 

 

2 w

 

2 w

 

 

2 w

 

~

;

 

 

 

 

 

 

+ (1 ν ) ξη

ξη

 

 

 

ξ 2

 

 

 

 

η2

 

η2

 

 

ξ

2

dA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

kl

 

 

 

 

mn

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

mn

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q(ξ,η)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

P

 

 

 

 

w(ξPi ,ηPi

)+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ckl = ∫∫

 

 

 

 

 

 

 

 

 

 

w(ξ,η)dA + λ

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q0

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q0 a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

my (ξ,η(ξ))w

 

 

 

 

 

 

 

 

 

 

 

p(ξ, y(η))

w(ξ

,η(ξ))+

mx (ξ,η(ξ))w

 

 

 

 

 

~

 

 

 

+λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

λ

 

 

 

 

 

 

 

 

 

 

 

 

ds .

 

 

 

 

 

 

 

q

 

 

a

 

 

 

 

 

q a

2

 

 

 

 

ξ

 

q

 

a

2

 

 

 

η

 

 

 

 

 

 

S

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.3.8)

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

mn

~

2

~

2

wkl dA .

(6.3.9)

Bkl

= ∫∫

 

wmn

 

A

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

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

 

wmn (x,y)= X m (x) Yn (y),

(6.3.10)

где

X m (x),Yn (y) - функции, удовлетворяющие

условиям

опирания пластинки на краях х = 0, х = а и у = 0, у = b

соответственно, при этом удовлетворяются обычно балочные

условия опирания, а сами функции X m (x),Yn (y)

называют

балочными функциями.

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

Bklmn = Cxmk

Aynl

+νλ2 (Bxmk Byl n + Bxkm Bynl )+ 2(1 ν )Dxmk Dynl , (6.3.11)

mk

1

 

mk

1

 

 

 

′′

 

= X m (x) X k (x)dx ;

= X m (x) X n (x)dx ;

где Ax

 

Bx

 

 

0

 

 

0

 

mk

 

1

′′

mk

1

 

′′

= X m (x) X k (x)dx ;

= X m (x) X k (x)dx .

Cx

Dx

 

 

0

 

 

0

 

Коэффициенты с нижним индексом у получаются с соответствующей заменой функции Х(ξ) на функцию Y(η), dξ на dη.

6.3.Метод КанторовичаВласова в задачах изгиба пластинок.

Решение методом КанторовичаВласова для прямоугольных пластин ищем в виде одинарного ряда

 

w(ξ,η)= G X m (ξ)Ym (η),

(6.4.1)

 

m=1

 

где

Ym (η) - функция, удовлетворяющая граничным

условиям

опирания пластинки на продольных краях η = 0, η = 1

(y = 0, у =

b);

X m (ξ) - неизвестные функции; G - коэффициент, задаваемый

в процессе решения для придания ему более удобной формы записи.

Для нахождения неизвестных функций X m (ξ) используем

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

1 1

4 w q(x, y)

a4

δ w(ξ,η)dξ dη =0.

(6.4.2)

 

 

 

∫∫

 

 

 

 

D

 

 

 

0 0

 

 

 

 

 

 

 

 

Так как варьируются только неизвестные функции

 

 

 

δ w(x, y)= G Ym (η) δ X (ξ),

 

 

 

 

 

 

 

m=1

 

 

то, разделяя интегрирование, получим

 

 

1 1

~

4

 

 

a

4

 

 

 

 

 

 

 

 

 

 

 

∫ ∫

 

w q(x, y)

 

 

 

Yn (η)dη δ X n (ξ)dξ = 0 .

 

D

 

 

 

 

 

 

 

 

 

 

 

0 0

 

 

 

 

 

 

 

 

δ X n (ξ) на

Учитывая независимость

 

вариаций

функций

основании основной леммы вариационного исчисления, получим систему функциональных уравнений метода КанторовичаВласова

 

1

 

~

4

 

a4

 

 

 

 

 

 

 

 

w q(x, y)

 

 

 

Yn (η)dη = 0 ,

n = 1, 2, 3…

(6.4.3)

 

 

 

D

 

0

 

 

 

 

 

 

 

 

 

 

 

Подставляя в соотношения (6.4.3) решение (6.4.1), получим

1

[X m

 

(ξ)Ym (η)+ 2λ

 

X m (ξ)Ym

(η)+ λ

 

X m (ξ)Ym (η)] Yn (η)dη =

G∫ ∑

 

2

4

 

IV

 

 

 

 

 

′′

′′

 

IV

 

0 m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= a4 1 q(ξ,η)Yn (η)dη . D 0

Принимая G = q(ξ,η)a4 и интегрируя ряд почленно, получим q0 D

систему обыкновенных дифференциальных уравнений для

функций

X m (ξ)

 

 

 

 

X m (ξ)]= qn , n = 1, 2, 3,…

nm

 

IV

2

nm

′′

4

~nm

 

 

[Ay

 

X m

(ξ)+ 2λ

By

X m (ξ)

+λ

Cy

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.4.4)

 

 

Aymn = 1 Ym (η) Yn (η)dη ;

 

Bymn = 1 Ym′′(η) Yn (η)dη ;

где

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

~mn

1 IV

 

 

 

 

 

1 q(ξ,η)

 

 

 

 

 

Cy

 

= Ym (η) Yn (η)dη ;

qn =

 

 

 

 

Yn (η)dη .

 

 

 

 

 

q0

 

 

 

 

 

 

 

 

0

 

 

 

~mn

0

 

 

 

 

 

 

 

Проинтегрируем выражение

дважды по частям

 

 

Cy

 

 

 

~mn

 

1

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

IV

 

 

′′′

 

 

 

 

 

 

′′′

 

 

 

 

 

 

 

 

 

 

 

 

 

Cy = Ym (η) Yn (η)dη =Ym (η) Yn (η)

0 Ym (η) Yn (η)dη =

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

′′′

 

 

′′

 

 

 

1

 

1 ′′

(η)dη .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= [Ym (η) Yn (η)Ym (η) Yn

(η)]

0

Ym

(η) Yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

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

~mn

1

1

′′

 

 

 

IV

′′

mn

.

(6.4.5)

Cy

= Ym

(η) Yn (η)dη = Ym (η) Yn (η)dη = Cy

 

0

0

 

 

 

 

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

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

независимых дифференциальных уравнений для каждой функции

X m (ξ). Это возможно, когда система функций Ym (η)

решения

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

mn

mn

 

 

~mn

mn

= 0 при m n.

(6.4.6)

Ay

= By

= Cy

= Cy

Тогда получаем систему независимых обыкновенных

дифференциальных уравнений 4-го порядка

 

IV

(ξ)+ 2λ

2

 

′′

4

Cm X m (ξ)= Em ,

(6.4.7)

X m

 

Bm X m (ξ)+

λ

где

Bm =

Bymm

; Cm =

Cymm

; Em =

q

m

.

Aymm

Aymm

 

 

 

 

 

 

Aymm

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

коэффициент Bymn при m n для динамических балочных

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

Характеристическое уравнение дифференциального уравнения (6.4.7)

r4

+ 2λ4 B r 2

+ λ4C

m

= 0 .

(6.4.8)

m

m m

 

 

 

Как показали практические расчеты, это уравнение имеет два типа корней:

1)четыре комплексно-сопряженных корня (общий случай)

rm,14 = ±(αm ± iβm ),

(6.4.9)

 

где i - мнимая единица, а α и β

определяются формулами

αm = λ

Cm Bm ;

βm = λ

Cm + Bm .

(6.4.10)

 

2

 

2

 

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

Во втором случае решение получается при шарнирном опирании продольных сторон пластинки, и решение КанторовичаВласова совпадает с решением Леви.

Остановимся на общем случае характеристических корней. В этом случае общее решение дифференциального уравнения (6.4.7) имеет вид

X m (ξ)= 4

AmiΦim (αm ,βm ,ξ)+ X m0 (ξ),

(6.4.11)

 

i=1

 

 

 

 

где Φim (αm ,βm ,ξ) - независимые

частные решения однородного

дифференциального

уравнения

(6.4.7);

X m0 (ξ)

- частное

решение неоднородного

дифференциального

уравнения; Ami -

константы интегрирования, определяемые из условий опирания пластинки на поперечных краях пластинки ξ = 0, ξ = 1 (х = 0, х =

а).

6.5. Метод БубноваГалеркина в задачах изгиба пластинок.

В основу решения задач изгиба пластин методом БубноваГалеркина положено функциональное уравнение (6.2.7), полученное при доказательстве принципа Лагранжа.

Решение принимается в виде двойного ряда

 

q0a

4 ∞ ∞

 

w(ξ,η)=

 

∑∑Amn wmn (ξ,η),

(6.5.1)

D

 

 

 

m=1n=1

 

где ξ, η - безразмерные

координаты (см. формулы

(6.3.2));

а - характерный размер пластинки в плане; q0 - произвольное значение нагрузки; w(ξ,η) – функция, удовлетворяющая всем

граничным условиям опирания пластинки; Аm,n - неизвестные коэффициенты.

Функциональное уравнение (6.2.7) при переходе к безразмерным координатам имеет вид

~

4

w

q(ξ,ξ)a

4

wkl (ξ,η)dA = 0 ,

k, l = 1, 2, 3,… (6.5.2)

∫∫

 

 

 

 

 

 

 

 

 

D

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

4

.. =

 

4 ..

 

 

2λ

2

4 ..

 

4 ..

 

Здесь

 

 

 

 

 

 

+

 

 

+

 

- бигармонический

 

 

 

ξ 4

 

 

ξ 2η2

ξ 4

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

Подставляя в уравнение (6.5.2) решение (6.5.1) и интегрируя, получаем систему алгебраических уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bklmn Amn = qkl ,

k, l = 1, 2, 3,…

 

 

 

(6.5.3)

 

 

 

 

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

(ξ,η)

 

 

 

 

 

 

(ξ,η)

 

4 w (ξ,η)

 

mn

 

 

 

4 w

 

 

 

2

4 w

 

 

wkl (ξ,η)dξ dη ;

=

∫∫

 

 

mn

 

 

+ 2λ

 

mn

 

+

 

mn

 

 

Bkl

 

 

ξ

4

 

 

 

ξ

2

η

2

 

ξ

4

 

 

A

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

kl

∫∫

 

q(ξ,η)

w

kl

(ξ,η)dξ dη .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для прямоугольной пластинки функция w(ξ,η) обычно принимается в виде произведения функций независимых аргументов w(ξ,η) = Хm(ξ) Yn(η). Тогда коэффициенты системы алгебраических уравнений можно определять по формуле

 

mn

~mk

nl

mk

nl

+ λ

4

mk ~nl

;

(6.5.4)

 

Bkl

= Cx

Ay + 2λBx

By

 

Ax Cy

mk

1

 

 

 

mk

1

 

′′

 

 

где Ax

= X m (x) X k (x)dx ;

Bx

= X m (x) X k (x)dx ;

 

0

 

 

 

 

0

 

 

 

 

~mk

1

IV

 

1

′′

 

 

′′

mk

;

 

 

 

 

 

Cx

= X m (x) X k (x)dx

= X m (x) X k (x)dx = Cx

 

0

 

 

0

 

 

 

 

 

 

Коэффициенты с нижним индексом у получаются, соответствующей заменой функции Х(ξ) на функцию Y(η).

Решая систему алгебраических уравнений (6.5.3), определяем коэффициенты Аmn.

6.6. Пример расчета пластинки на изгиб.

Задача. Провести расчет прямоугольной пластинки методами РитцаТимошенко и Канторовича–Власова на равномерно

распределенную нагрузку q0

в первом приближении при

различных соотношениях сторон

λ =

a

, правый поперечный край

b

 

 

 

пластинки (х = 0) жестко защемлен, остальные шарнирно оперты (рис.6.4). Коэффициент Пуассона принять ν = 0,3.

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

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

Решение. В соответствии с граничными условиями опирания пластинки на левом поперечном шарнирно опертом крае (ξ = 0 )

имеем

w(0,η)= 0 ,

 

 

M x (0,η)= 0 2 w(0,η), на

правом

жестко

 

 

 

 

 

 

 

 

 

 

 

ξ 2

 

 

 

 

 

 

 

защемленном крае (ξ =1) w(1,η)= 0 ,

 

w(1,η) = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

а расчет методом РитцаТимошенко.

 

 

 

 

 

Граничные условия соответствуют условиям опирания балки-

полоски,

вырезанной из пластинки

в поперечном направлении:

X (0)= 0 ,

X

′′

 

 

 

 

,

 

 

 

Балочную функцию в

 

(0)= 0

X (1)= 0 , X (1)= 0 .

 

продольном направлении принимаем в виде:

 

 

 

 

 

 

 

 

3

 

 

 

4

 

2

 

 

′′

 

 

 

 

 

X (ξ)=ξ 3ξ + 2ξ

 

 

=ξ(1ξ) (1+ 2ξ);

 

X

 

 

 

 

 

 

 

 

 

(ξ)= −6ξ(3 4ξ);

 

 

 

 

2

 

 

 

3

= (1 ξ)(1 +ξ 8ξ

2

); X

′′′

 

 

 

 

 

X (ξ)=1 9ξ

+ 8ξ

(ξ)= −6(3 8ξ). (6.6.1)

 

 

 

 

Соответственно на продольных краях имеем граничные

условия опирания балки-полоски:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

′′

 

 

 

′′

 

 

 

 

 

 

 

 

 

 

 

Y (0)= 0 , Y (0)= 0 , Y (1)= 0 Y (1)= 0 ,

 

 

 

 

которым соответствует балочная функция:

 

 

 

 

 

 

 

 

 

 

 

3

 

 

4

=η(1 η)[1 +η(1 η)];

 

 

2

 

3

 

Y (η)

=η 2η

+η

 

Y (η)

=1 6η

+ 4η

;

 

 

 

 

 

′′

 

 

 

 

 

η)

;

′′′

 

 

 

 

 

 

(6.6.2)

Y (η)= −12η(1

Y (η)= −12(12η).

 

 

Легко убедиться, что все граничные условия опирания пластинки выполняются. В первом приближении решение

методом РитцаТимошенко ищем в виде

 

w(ξ,η)=

q0 a4

A X (ξ) Y (η).

(6.6.3)

D

 

 

 

 

 

Числовые индексы в первом приближении опускаем.

 

Коэффициент А определяется из уравнения

 

B A = q ;

A =

q1

,

(6.6.4)

 

1

 

B

 

 

 

 

 

где В – коэффициент, определяемый в соответствии с формулами

(6.3.9), (6.3.11)

 

при m = n = k = l =1:

 

 

 

 

 

 

B = C

x

A

y

+ 2λ2 B

B

y

+ λ4 A C

y

;

q = q

x

q

y

.

(6.6.5)

 

 

x

 

x

 

 

1

 

 

 

При q = q0 = const ;

 

qx = 1

X (ξ)dξ ;

 

qy = 1 Y (η)dη .

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

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

 

с

формулами

(6.5.4)

вычисляем

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

Ax = 1

X 2 (ξ)dξ = 1

ξ 2 (1 ξ)4 (1 + 2ξ)2 dξ =

2! 4!

+ 4

3!4!

+ 4

4!4!

=

19

;

7!

 

8!

 

9!

 

630

0

0

 

 

 

 

 

 

 

 

Bx = 1 X ′′(ξ)X (ξ)dξ = −61 ξ 2 (1 ξ)2 (3 + 2ξ 8ξ 2 )dξ =

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −6

 

 

2! 2!

+

3!2!

8

 

4!2!

= −

12

;

 

 

 

 

 

 

3

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

5!

 

6!

 

7!

35

1

2

1

2

 

 

2

 

 

 

 

 

9

 

24

 

16

 

 

36

 

 

Cx = X ′′

(ξ)dξ = 36ξ

(3

4ξ)

dξ = 36

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

=

 

 

;

 

 

 

3

4

 

5

5

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qx

= 1

X (ξ)dξ = 1

ξ(1ξ)2 (1+ 2ξ)dξ =

2!

+ 2

2!2!

=

 

3

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

4!

 

 

 

 

5!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ay = 1 Y 2 (η)dη = 1η2 (1 η)2 [1 + 2η(1 η)+η2 (1 η)2 ]dη =

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2! 4!

 

+ 2

3!3!

 

+ 4

4!4!

=

 

 

31

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7!

 

 

 

9!

 

 

630

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

By = 1 Y ′′(η)Y (η)d = −121 η2 (1 η)2 [1 +η(1 η)]dξ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2! 2!

 

 

 

 

 

3!3!

 

 

 

 

17

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −12

 

 

 

 

 

+

 

 

 

= −

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5!

 

 

7!

35

 

 

Cy

= 1 Y ′′2 (η)dη =1441

η2 (1 η)2 dη =144

 

2!2!

= 24

;

 

 

 

 

 

 

 

 

 

 

 

5!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qy

= 1 Y (η)dη = 1

η(1 η)[1 +η(1 η)]dη =

1

+

 

 

2!2!

=

1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5!

 

 

5

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

3!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1 =

3

1 = 0,03 ; B =

36

 

31

+ λ2

2

 

12

 

 

 

17

+ λ4

19

 

24

=

 

 

 

 

630

35

 

 

 

35

 

 

5

 

 

20

5

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

630

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0,3343 + 0,3331λ2 + 0,1448λ4 .

 

Значения коэффициентов

В, А

для различных отношений

сторон λ приведены в табл. 6.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 6.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ =

a

 

 

 

0,5

1,0

 

 

1,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,0

 

 

 

 

 

 

 

 

2,5

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

 

 

 

04266

0,8122

 

1,817

 

 

 

 

 

 

 

 

3,984

 

 

 

 

 

 

 

8,072

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

0,07032

0,03694

 

0,01651

 

 

 

 

0,007531

 

0,003716

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При вычисленных значениях коэффициента А, прогибы в пластинке определяются по формуле (6.6.4), а внутренние усилия с учетом перехода к безразмерным координатам по формулам:

2

 

 

′′

 

 

2

′′

 

 

 

 

 

 

 

 

M x (ξ,η)= −q0a A [X (ξ) Y (η)+

νλ X (ξ) Y (η)];

 

2

 

′′

 

 

2

′′

 

M y (ξ,η)= −q0a A [ν

 

 

 

 

X (ξ) Y (η)+ λ X (ξ) Y (η)];

 

Qx (ξ,η)= −q0a

A [X

′′′

2

′′

(6.6.6)

 

(ξ) Y (η)+ λ

 

 

X (ξ) Y

(η)];

Qy (ξ,η)= −q0a

 

′′

 

2

 

′′′

 

 

 

 

 

 

A [X (ξ) Y (η)+ λ

 

X (ξ) Y (η)].

 

б решение методом КанторовичаВласова.

Решение в первом приближении ищем в виде

w(ξ,η)=

q0a4

X (ξ) Y (η).

(6.6.7)

D

 

 

 

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