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

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

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

где

Хт — корни

уравнения

J l (X) =

0;

А %

В ? ,

— произвольные

постоянные.

 

 

 

 

 

 

 

 

 

 

 

Компоненты перемещений

vtz

на основе (2.73), (5,32) примут вид

 

 

B

P r -

S £

C{Jl]iXmc h 4 ^ = - J i M

+

 

 

 

i= 1

 

ш=1 i=l

 

 

V * i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ И

£ X '.- V * А Л

(

У

cos knz,

 

 

 

 

 

/1=1 t=l

 

 

 

 

 

 

 

 

 

uf = _ 2 1 м » A . , + j

£

 

К

s h - ^ J o ( M -

 

 

i=l

*

m=l i=l

 

V Hl

 

 

 

 

oo

2

 

 

 

 

 

 

 

 

 

 

~ ■S

S

 

( К х Л / ) sin £nz

(a„ =

 

,

(5.33)

 

n—l i=I

 

 

 

 

 

'

 

'

 

где

kc — корни

алгебраического

уравнения

(2.64).

 

 

 

Согласно (2.72) при переходе к изотропному телу функции ф{л (г, г)

иФг* (г, z) становятся линейно зависимыми и, следовательно, из

(5.33) не получаем вы раж ений, соответствующих (5.13).

В

случае

кратны х корней х, (хх =

х 2) функции Ф (Л (г, г) и Ф2Л (г, г)

такж е

становятся линейно зависимыми. Поэтому, оставляя, например,

Ф{л (г, z)

в прежнем виде (5.32) при

г =

I, функцию Фгу) (г, г) следует

выбирать

в

форме

 

 

 

 

оо

Ф ^ (г , z) — ВтРг2 +

 

 

 

оо

 

 

+

s

A & rh (У*шкпГ) cos knz +

£

C & s h - ^ / o M - (5.34)

 

 

п=1

 

m=l

К ^2

Если в соответствии с (2.70) корни х х и х 2 характеристического уравне­

ния (2.65) являю тся комплексно-сопряженными, то функции Ф/Л (г, г) могут быть выбраны в виде (5.32), однако произвольные постоянные

тож е

должны

быть комплексно-сопряженными, т. е.

Ап\ = А п\,

В ? =

с £ ,

= С‘л 2.

 

Если на основе функций (5.32) найти вы раж ения для

напряжений

по формулам типа (2.73), то краевые условия (5.5) для касательных

напряж ений

(с учетом того что J t (кт ) = 0 и sin knh =

0)

приводят к

зависимостям

 

 

 

 

 

 

С % - У ^ - sh - ^ 4 - = — С &

sh - ^ r -

с

Я.2

 

 

К «2

F x ,

Kxx

У*х

 

 

С44ЛОТ

(5.35)

 

 

 

 

 

 

 

 

y|l/>_ _

/I(/)

(1 + *1) v^«7/a (K x iy

 

 

 

 

л п.2-------/l/i.l

---------=----— ----- — ---------

 

 

 

 

 

 

(l + A2) / « 2/l(^«2ftn)

 

 

 

 

 

 

Ф$

 

 

 

 

 

 

«44 (И -

 

 

 

 

позволяющим

исключить из

рассмотрения

постоянные

А % и

 

151

Краевые условия (5.5) для нормальных напряж ений на основе р а з ­

лож ения

типа

(5.21),

(5.22) приводят

к системе алгебраических

у р ав ­

нений второго

порядка

относительно

постоянных

B p

(t =

1, 2)

 

 

 

 

 

 

 

2

 

 

 

оо

ф

 

 

 

 

 

 

 

 

 

 

-

 

( с „ - с ,о

S

fil' ’ +

т

2

-

г - ; '

ы

=

 

 

 

 

 

 

 

 

 

 

 

 

1=1

 

 

m=1

 

 

 

 

 

 

 

 

 

 

 

_

2 £

 

 

 

 

 

gjo +

2 £

(_

1}»

 

_

 

ц л

 

 

(5.36)

 

 

 

,t=l

 

 

*

 

 

 

 

п=1

 

 

л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и к бесконечной системе линейных алгебраических уравнений

 

 

Л!)

__ *44 (1 +

^l)

00

'

 

1

 

 

 

 

1

 

 

 

ш(/>

 

 

 

 

 

 

 

 

 

 

 

 

хя

 

 

Ли

--

 

т п

 

s

 

</S?(

 

 

Km

 

kl Ki +

 

.) +

 

Гя

 

 

 

 

 

 

 

m=l

. Йпх1 +

 

 

 

 

 

(5.37)

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

Ф<£

 

 

*44(1-M l)

 

 

«1

 

 

 

Xa

 

'

 

 

 

е й

- s

* ? (

 

 

 

 

 

 

 

m

*

 

-

 

 

 

 

 

 

 

 

 

fenX8 +

 

;) +

Sm

 

 

 

 

 

 

 

/i=I

. # » + £

 

 

 

 

 

 

где вместо

А р и

Cm!i

введены

новые

переменные

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.38)

 

 

 

 

уЯ

=

- т ё = - С Й Л

 

sh - i s

t

(>.»,)■

 

 

 

 

 

 

 

 

 

 

 

 

 

Л К «1

 

 

 

 

У

 

 

 

 

 

 

 

 

 

Вы раж ения

для

Т„, S m,

 

 

имеют вид

 

 

 

 

 

 

 

 

 

у

 

1

Г

h

 

k2

с1г — с1а .

сД4 (1 -|- fej)

/р (УмдАп)__

 

 

 

 

2 L i + L

 

 

 

 

 

 

 

 

h < y ^ k n)

 

 

 

 

 

 

 

 

 

 

 

*44 U + Л )

 

/р (У*^л) 1

 

 

 

 

 

 

 

(5.39)

 

 

 

 

 

 

 

 

ka Уха

/ 1 (У игАп) J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S - =

ж

г

Я

 

+ * •>

 

с4Л

 

 

^

dA

 

 

 

 

( - 1 ) " у У

= / й - < р й | -

Сп—

-У—

+

 

 

 

 

,

1

-

 

 

 

 

 

 

 

 

 

L

с41 (1 +

/г2) kn

У *аЛ (У « aU

J

 

 

 

 

 

 

 

 

 

2

(— l)m Ят Ф&

(^ш)»

 

 

 

 

 

 

 

 

 

 

 

-

i

s

 

ф . +

ч

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

+ < &

У щ cth

 

j 0 (к ,) -

г

s

ф й

 

 

 

,

 

 

 

 

 

 

 

 

 

 

У ха

 

 

 

л=.

 

 

* п « а + ^

где tpia, фгт. /Й .

 

 

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

Ф урье и Ф у р ье — Д ини

рядов

(5.18).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д л я бесконечной системы

(5.37)

выполняется асимптотический за ­

кон (5.27). Это устанавливается аналогично случаю

изотропного тел а

1191. Следовательно, для достаточно больших

полож ительны х

целы х

152

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чисел

N и М выполняю тся приближенные равенства (5.28). Это позво­

л я е т

(в предполож ении выполнения условий

парности касательных

н ап р яж ен и й по угловым линиям г — 1, z = ±

К) построить эффектив­

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

Замечание. Соответствующие результаты , относящиеся к краевы м задачам для полых изотропных и трансверсально изотропных цилиндров с осесимметричными выточками на боковых поверхностях,,

являю тся

более громоздкими

и

здесь не приводятся. Некоторые из=

них будут

приведены в п. 1.2

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

полых цилиндров с неплоскими

торцами.

1.2.Цилиндры с выточками на торцах. Рассмотрим полый кру­

говой цилиндр конечных размеров, боковые S 0, 5 t и торцевые S t поверхности которого в цилиндрической системе координат г, 0, z (г, z — безразмерные переменные, отнесенные к внешнему радиусу) опи­

сываю тся

 

уравнениями

 

 

 

 

 

 

 

 

S 0 '—' r

= a,

S

±

г = 1,

S t

z =

^ z h -f- £<в±/ (г).

(5.40)'

П араметры

е,

(о± ,

а

и

h

удовлетворяю т

следующим

неравенствам:

0 < е

1,

— 1

^

(о± ^

 

1, 0 <

а < 1,

h

= const >

0. Предполо­

жим, что требуется исследовать осесимметричное напряженно-дефор­ мированное состояние цилиндра при следую щ их граничных условиях:

{ОгтП? + <5гтП г)„± =

(/72 = Г, z),

Grm |s0 = Тт0»

(5.41)

&гт |s, — Тт [,

где t i t , n f — направляю щ ие косинусы единичных нормалей п* к тор­

цевым поверхностям

S t \

T t ,

Тт о, Тт \ — известные функции.

Если

решение задачи искать в виде рядов (5.4), то в /-м приближе­

нии можно получить следующие

краевые условия

на

координатных

поверхностях:

 

 

 

 

 

 

 

 

 

 

o g

|г=±Л =

Т ? »

-

£

[Q tu•о й -»

+ Q i l* o i - % = ±h =

f t и) (г),

 

 

 

S=1

 

 

 

 

 

 

 

а й

-

Г « л -

£

ю ? « а й -*

+

 

=

Ф?М «

 

 

 

 

S=1

 

 

 

 

 

(5.42)

 

 

 

 

 

 

 

 

 

 

 

 

ст/г |г=а =

т%

=

$

(г),

о(Л \г=а =

Т% =

ч4л (г),

 

4 ? U . =

т $

=

f i ] (2),

4/1 U i =

T /i =

ф£л (г),

где Q ?(s), Q t(s) — дифференциальные операторы, которые в общем слу­ чае имеют вид (3.38), а при s = 1, 2, 3 — вид (3.141). При записи вы­ раж ений для перемещений и напряж ений используем результаты [191.

В предположении, что изотропный цилиндр деформируется симмет­ рично относительно плоскости г = 0, вы ражения для перемещений,.

153-

удовлетворяю щ ие

уравнениям равновесия и

построенные аналогично

(5.13),

будут иметь

 

вид

 

 

 

 

 

 

 

 

 

 

 

 

_

 

п(/‘)

 

 

 

 

 

 

 

 

 

 

 

и ?

-

B f r

+

- 5 -

+ £

{A fzsh l {z +

C f ch V ) N { (V )

+

 

 

 

 

 

 

 

 

Г

1=1

 

 

 

 

 

 

 

 

 

 

 

+

£ ( |М Л ( M

 

+

r Q f (knr) — 4

* -~ v -

(knr )j cos Ляг,

 

 

 

 

 

 

oo

 

 

 

 

 

 

 

 

 

 

 

 

(5.43)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u f =

A f z -f-

(~3 ~^4v

sh A-,-2 — A f z ch h(z — C f sh %tzj X

 

 

 

 

X N 0 ( V ) -

2

[РФ (V ) +

r Q f {knr)} sin knz.

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

З д е с ь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P f

 

(knr) =

£ ^ / 0 ( V )

+

^ / C o

( V ) .

 

 

 

 

 

 

 

Qort ( V ) =

Gjploiknr) +

Я < Х

(V )>

 

 

 

 

 

 

knP\n M

=

 

d P ^ ( M

>

W

</)

 

) “

dQg *( M

 

 

 

 

 

 

 

' V

" '

( V

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

причем

A q \ B $ ,

D f , A f ,

C f E f ,

F f , G f,

H f — произвольные

по­

стоян ны е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующие

 

перемещениям

(5.43) вы раж ения д л я

компонен­

то в

напряж ений

с

использованием закона

Г ука

принимаю т

форму

 

 

 

1

_(Л

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£

( A f z sh \ z - f

 

 

 

 

 

 

2G

а "

 

 

1 — 2v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

OQ

 

 

 

 

 

 

 

 

 

 

-f- C f ch

X(z) —

N t (X,r) +

£

(i4i;,^ z sh

+

2 v A f c h ^ z

-f-

 

 

 

 

 

 

r

 

 

 

<=i

 

 

 

 

 

 

 

 

 

 

 

+

с!"»., ch м

 

jv„ ( V ) +

f

[

 

( V ) -

~

P f 1 M

-

 

 

 

 

 

 

 

 

 

 

 

 

/1=1 L

 

 

 

 

r

 

 

 

 

— (3 — 2v) Q(,y) (knr) +

knrQ p (knr) +

4

L~

v - Qly) ( V ) ]

cos

 

1

„(/)

 

4 /] +

vAp

 

f

 

2

(A fz sh %tz - f

o f

ch V ) - j -

 

( V )

+

2G

ff00 “

1 — 2v

 

•"

л2

 

 

t=l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 2v £ A f ch M 'f . t V ) + s f - Г W ( V ) +

 

 

 

 

 

 

 

 

f=i

 

 

 

 

 

« = iL

 

 

 

 

 

 

 

 

+

(1 -

2v) 0У> (fe„r) - 4

 

Q P (*„r)] cos

 

 

 

154

2vB(J )

 

 

 

 

 

 

 

(5.44)

 

— A \% z sh X(z — d / %

ch X£z] N 0 ( V ) — 2

[ M o ’1(knr) +

 

 

 

 

 

n=1

 

 

 

oo

+ 2vQo' ( M

- f ^„г<3|л ( V ) ] cos knz,

 

 

 

 

 

 

 

 

2G

аЙ =

2

[ 4 % г ch Xtz — (1 — 2v) A P sh Xtz +

CpXc sh Xtz) N t (X£r) —

 

t= I

 

 

 

 

 

 

 

OQ

 

 

 

 

 

 

-

2

[knp p ( M + knrQ P (knr) -

2 (1 - V) Q P ( M

l sin knz,

где

 

 

N 0 (X£r) — J i (^ui) У о ( \ r ) —

Ух (A.,) J о (X£r),

 

 

 

 

 

 

 

 

( V ) = J i (h ) Yг ( V ) - Уг (h ) Л M •

 

Здесь

Ут , Km, Ym — общ епринятые обозначения

цилиндриче­

ских функций (функций Бесселя); Хс — корни уравнения Мг (Ха) = 0.

Представим

функции, входящ ие в правые

части краевых условий

(5.42), рядами

Ф урье — Д ини, коэффициенты

которых определяются

по формулам

 

 

 

а

а

 

 

(5.46)

При удовлетворении граничным условиям (5.42) для касательных на­

пряж ений

при

г =

1, а

и

z = ±

h

с

учетом

равенств

(Х£а) =

= N t (Х() =

sin knh — 0

приходим

к

соотношениям

 

 

 

 

 

 

 

 

 

 

„(Л

 

 

СРХ£ =

А р (1 - 2v -

X(h cth Щ +

,

 

К Р Р ( М ) =

2 (1 —

v) Р р

(kna) — krlaQo] (kna) — фз«,

(5.47)

 

K P P

{kn) = 2 (1 -

v) P P

(kn) -

knQP (ft,,) - <pP,

 

позволяющим исключить из рассмотрения часть переменных.

 

Д л я перехода

от функциональных

уравнений

(5.42) к системе ал­

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

155

Соответствующий

ряд

для Qo1 (knr)

получается из

(5.48) заменой Р\Р

на Q\P с учетом,

что

knP \ ] (knr)

=

d P § {knr)!dr. У довлетворяя крае­

вым

условиям для Czz при z — ±

h и оУг при г =

1, а , получаем сис­

тему

трех уравнений

относительно

ЛоЛ, В о \ D o1

 

и бесконечную систему алгебраических уравнений д л я определения оставш ихся в результате зависимостей (5.47) неизвестны х постоянны х

где

i f f =

~ 2 G ( - 1 ) " k j j f f (kn),

= 2G ( - 1 ) " ak nQ\n ( M ) ,

 

 

A0 =

K t (kna) / , (k„) —

К г (kn) / , (kna),

(5.51)

 

 

A i =

K i (kna) I 0 (kn) +

K 0 (kn) I г {kna),

 

 

 

A2 =

K 0(kna) I г (kn) +

К г (kn) I 0(kna).

 

Члены

ф Ц

Фга,

Фз?

представляют определенные

комбинации

коэффициентов

разлож ения

функций,

задаю щ их граничные условия

156

(5.42):

 

 

 

 

 

 

 

 

 

ф Й

=

( - 1 ) ”

f(/')

 

^2(1

,

,_(/) / Af

I \

/ з « -

+

срз„

* ~ j

 

 

 

 

 

 

&пД0

 

 

 

 

 

 

_

J L

V

(/)

ljA/0 (k,g)

 

 

 

 

mw

 

 

 

 

 

 

 

h

2 j

*£ + *?

 

 

 

 

 

 

i=

 

 

ф р

 

 

 

 

 

ч>з«

 

ф2п

 

f2tl =

 

 

[ «

+

М Л .

 

т Ь “) ] ~

 

 

 

 

 

"'(~S7 +

 

 

 

 

 

 

(/)

M o (A,t-fl)

(5.52)

 

 

 

-

4

2

r f

<

+ »!

 

 

 

 

 

i=i

 

 

ф $

=

-

/ I X

(a,) -

ФЙ* cth k{h N 0 (M

+

,

V

.

n ,I .

Nl fa) Фз/t — aN« fa) ^0 ( М

Ф(2п

+

n=l

 

"

 

 

 

+

 

Найти асимптотические оценки для неизвестных величин бесконечной системы (5.50) значительно слож нее, чем в случае сплошных цилинд­ ров. Учитывая соответствующие результаты [19] для полых круговых цилиндров, торцы которых совпадают с координатными плоскостями, и тот факт, что в каж дом приближении решается краевая задача именно для кругового цилиндра с плоскими торцами, предполагается, что

свойства неизвестных

у п( , z(,P бесконечной системы (5.50) описываются

асимптотическими

равенствами

 

 

 

 

 

 

lim

ylC =

у о \

lim zj? =

zip.

 

(5.53)

 

 

rt-voo

 

 

 

 

гг—*-оо

 

 

 

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

для

больш их

значений N и М можно полагать

 

 

 

Уп1 =

у о \

Zn] = г(0п

( V n > N),

 

 

 

 

 

 

 

 

 

 

 

(5.54)

л п

_

1 ^

{Xl)

..</>

,

1

No (h) N0 (Х{а)

Jf)

(V i > M).

Xi

~

2

R,

У°

' f_

2

R(

 

 

 

С точки зрения построения численного алгоритма расчета осесимметрич­ ного напряженного состояния полых цилиндров с торцевыми выточ­ ками принципиальным отличием от соответствующего алгоритма рас­

чета

цилиндров с выточками

на

боковых поверхностях является то,

что

в правых частях уравнений

(5.42) (в отличие от уравнений (5.5))

ф игурирую т производные по

переменной z (а не по г) от напряжений

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

Рассмотрим, например, производную

 

 

N0 (М

дг

+

 

157

 

+

S

 

 

0 (^W) "f—

Qo iHmf) “t“ fotnfQi (^m^)j| Kn s *ri kmZ —

 

 

 

 

 

 

 

 

 

OO

 

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

— h £

 

Л4</ )

+

£

 

 

sin &mz.

 

 

 

(5.55)

 

 

 

 

 

 

 

 

i=i

 

 

/71=1

 

 

 

 

 

 

 

 

И спользуя асимптотические выражения

при

больш их

/,

п

для

ф унк­

ций

Бесселя

и

корней

уравнения

N x (Ха) =

0 [133],

получаем

 

 

М (Р

«

М*(/) -------- ехР

 

6е (h — ZU {уоЛ

[т»

 

 

+

 

,

3 (аа — 1)

 

 

, J .

(— !)f 4У)

Г

,

ч .

3 (а3 — 1)

rt

 

. Л \

+

— 8яш

' C0S

<r>J +

 

у=а---- [ *

W +

 

16я^—

C0S “ * W J) ’

 

 

 

 

 

 

8i =

- j 4

s - ’

а * (r) =

8t- (1 — г),

 

 

 

 

 

 

 

 

 

9 (1 — а) (I — г)2

 

 

 

, .

 

З г+ 1

.

м .

 

 

 

ъ

(г) =

— -— ё е т — ~

 

cos а{ W —

§ г “ sm а <

+

(5 -56>

 

+

~6l a r ri + 1 )~ C0S «*

+

-(l5 - - l 62

8

^ (1~ g- C0S «*• W +

 

 

 

 

 

 

 

 

+ 6tcos a { ( r ) ------J p

 

sin

a , (r),

 

 

 

 

 

 

 

 

 

 

 

,

 

f_IVяf/Л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lm «

I « w =

Y 'r

°

exp l— km (1 — r)]

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[*«(1

- r )

3r2 + 6r — 1

 

33r3 — 27r* + 3r — 9

 

 

 

 

 

 

 

 

8r

 

 

 

 

 

\28r2kn

 

] •

 

 

 

 

 

 

m “ U) =

£

m : u\

l* * (/) =

 

s

l;„(/) sin kmz.

 

 

 

 

 

 

 

 

 

 

i= 1

 

 

 

 

m=I

 

 

 

 

 

 

 

Значения

сумм

М *,(й и Г ' (Л

находят по формулам [126] и вычисление

производной

можно производить по формуле

 

 

 

 

 

 

 

 

 

 

 

d a ^

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

Зг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

£

( I #

-

L X ') sin /гт г

+

Г

,(й.

 

 

 

(5.57)

Аналогичным образом строятся

расчетные выражения и для

произ­

водных высших порядков. В результате проведенных исследований вместо бесконечных сумм (5.55) имеем конечные суммы (5.57). Это служ и т основанием для построения эффективного численного алго ­ ритма расчета напряженного состояния полых цилиндров с осесиммет­ ричными торцевыми выточками.

158

159

§

2. М етод возм ущ ения упругих свойств и МВФГ

в

неклассических краевы х зад ачах д л я неканонических областей

Краевые задачи для кусочно-однородных областей, близких к кано­ ническим, допускаю т эффективное приближенное решение (с требуе­ мой для приложений точности) с помощью развитых в гл. 2—4 вари­ антов М ВФГ в тех случаях, когда физико-механические свойства рассматриваемых тел позволяю т получить точное общее решение ос­ новных уравнений (равновесия, движ ения, теплопроводности и др.). В статической теории упругости к таковым, например, относятся изо­ тропные, трансверсально изотропные, частного типа неоднородные и некоторые другие тела, допускаю щ ие точное общее аналитическое ре­ шение уравнений равновесия. В более сложных случаях прямолиней­ ной и криволинейной анизотропии, физической нелинейности, неодно­ родности и других приближенное аналитическое решение краевых задач для неканонических областей можно получить на основе сов­ местного применения МВФГ и метода возмущения упругих свойств (если последние незначительно отклоняю тся от допускающих точное аналитическое решение уравнений равновесия в перемещениях в пря­ моугольных, круговых цилиндрических и сферических координатах). Иллюстрации такой возможности в общем и ряде частных случаев-

посвящен

настоящ ий параграф,

в основу которого положены

публи­

кации [65— 69, 77, 129].

 

 

2.1.

Общий подход. Предположим, что упругие свойства рассмат­

риваемого

многослойного тела

незначительно отклоняю тся

от тех

свойств тела, которые допускают точное общее решение уравнений рав­

новесия

в перемещ ениях. В связи с этим

представим

схематически за­

висимости

между

напряж ениями

ст,/,/ и

малыми деформациями

ецл в-

l-м

слое

в

виде

суммы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оцл

= o if,i‘

+ b K [,j.

 

 

 

 

 

(5.58)

Здесь

а,/,г — часть

уравнений

состояния,

допускающая точное общее

решение, например уравнений равновесия (2.13)

(для

изотропных

тел), или типа (2.62)

(для трансверсально

изотропных тел); к,-/./ —

возмущ ения,

характеризую щ ие те

отклонения

в уравнениях состоя­

ния,

которые

не

допускаю т

точное общее

аналитическое

решение-

уравнений

равновесия;

8 — безразмерный параметр

(0 ^

б ^

1).

 

Подстановка зависимостей (5.58) в уравнения равновесия в напря­

ж ениях

(2.12)

в отсутствие

объемных

сил

приводит к

уравнениям'

 

 

 

 

 

 

 

 

 

S NcfU,-, +

6FW =

0,

 

 

 

 

 

(5.59)

 

 

 

 

 

 

 

 

 

1=1

 

 

 

 

 

 

 

 

 

 

 

где

«/,/ — компоненты

 

перемещений /-го

слоя;

N (f — линейные диф­

ференциальные

операторы;

F,-,i — функции,

которые определяются

через

возмущения

х;/./

по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Гд (НоНяхи1)

д (Н1Н3,хХ!21)

д(Н хНгу,J3,/) 1

 

 

F l'1

=

l l j r j h

I

 

fai

^

 

да2

 

+

 

<>а3

 

J +

 

 

x,V

аЯ,

*13,/

дНл

*22,1

 

дНг

 

*зз,/

дНs

(1,

2,

3).

(5.60)

+

H JI2

даг

f

ЯхЯа

да3

Я ,Я а

дах

 

НХН3

дах

 

—---------

 

П р ед п о л о ж и м , что

требуется исследовать напряж енно-деф орм ирован ­

н о е состояние тел а

при заданны х на его неканонической поверхности

о

о

-S0 п ерем ещ ениях щ,о-, а на S n — уси ли ях тддг. С ледовательно, краевы е

у с л о в и я

на

гран и чн ы х

поверхностях

имеют

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

з

 

 

 

 

 

 

о

 

 

 

 

(5.61)

 

 

 

 

 

 

и1Л ISo =

и/,о,

Zj

a tf,Nn i,N |siw =

Т/.ЛГ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=l

 

 

 

 

 

 

 

 

 

 

 

 

 

:где

tii,N —• н ап равляю щ и е

косинусы

единичной

 

нормали

елг к

Sn-

В

сл у чае

идеального

 

контакта

на

поверхности

раздела

Si (/ =

• = 1, 2,

.... N

1)

условиям и

сопряж ения

будут

 

 

 

 

 

 

 

 

 

{Щ.1 —

w/,/+i)s, =

0,

Yj

(а и-1 — a i/./+i) п и Is, =

0.

 

(5.62)

 

 

 

 

 

 

 

 

 

1

 

 

 

f=.i

 

 

 

 

 

 

 

1

 

 

 

 

Д о п у с т и м ,

что граничны е

поверхности

S 0, S n ,

 

а так ж е

поверхности

р а зд е л а

S i

{I =

1, 2,

...,

N

1) являю тся иеоргогональны ми

и

опи­

сы ваю тся уравнен и ям и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a i =

 

Н~ еЛ

(а 2*

а з)

(<7 =

0, 1,

2,

 

. . . , N ),

 

 

(5.63)

где

 

f q ( а и

а 2) — непреры вны е

и

дифференцируемые

ф ункции,

гд —

м ал ы е

п арам етры

( |e e |< ^

1),

характеризую щ ие отклонение

поверх­

н о с те й

S Q от координатны х

поверхностей а х =

 

а°

=

const

(q =

0, 1,

2 , ..., N ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П ри

аналитическом

решении поставленной

 

задачи

возникаю т две

п р и н ц и п и альн ы е

м атематические трудности,

которые

вы нуж даю т р аз­

р аб аты в ать

 

приближ енны е

методы.

 

П ервая — слож ная

геометриче­

с к а я ф орм а

поверхностей

раздела

и

граничны х

поверхностей,

описы ­

ваем ы х

уравнениям и

(5.63). В таки х

случаях

 

на основе второго ва ­

р и ан та

М В Ф Г (гл.

3)

поставленная

задача

сводится

к

рекуррентной

последовательности

соответствую щ их

краевы х

 

задач для

составных

тел

 

с

координатны ми

поверхностями

раздела

<хг =

а ” (/ =

1, 2,

...,

N — 1) и граничны ми

 

поверхностями a L =

a<j,

 

= а%- Вторая —

слож ны е уравнения состояния, схематически

записанны е в виде (5.58),

в р езу л ьтате

чего

уравнения

равновесия (5.59)

не допускаю т точного

общ его

реш ен и я .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П о

этим

причинам

реш ение задачи

предлагается

искать

в

виде

двойны х рядов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U/.1 =

f

 

S

8kenu fjn\

о if,I

= У

У

6 V a f t7 \

 

(5.64)

 

 

 

 

 

 

k=0 /1=0

 

 

 

 

 

 

 

k—0 /1=0

 

 

 

 

 

 

где малый парам етр е выбран

по принципу (3.52), т. е. eo)e =

&q (—

1 <1

< ;

<о, < ; 1).

Д л я

 

удобства

интерпретации

разделим

ряды

(5.64) на

тр и

 

группы . П ервая

группа слагаем ы х (k

0, п =

0) с компонентами

u f j° \

crf/j

 

будет

характери зовать

 

реш ение

задачи

для

составного

тел а

с

услож ненны м и упругими

свойствами,

описываемыми

уравне­

ниям и

состояния

(5.58),

в

случае

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

(координатных) по­

в е р х н о с т е й

 

разд ела и граничны х поверхностей a L =

a° (q — 0,

1, 2,

160