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

1781

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

а площадь сечения Fj , статический момент S j и момент инерции J j ребра относительно произвольной начальной поверхности при помощи дельта-функции x xj , учитывающей дискретный характер расположения ребер, находятся по формуле

 

 

 

 

 

1

 

 

Fj x xj ;

 

F

j

 

 

 

h

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

j

 

 

1

S j x xj ;

(2.5)

S

 

 

 

 

 

 

 

 

 

 

 

h

 

 

j

 

 

 

j

 

1

J j x xj .

 

J

 

 

 

 

 

 

 

 

 

 

 

h

 

 

j

 

Компоненты вектора перемещения точки срединной поверхности оболочки показаны на рис.2.3. При этом через U(x,s) обозначено продольное перемещение вдоль оси х, через V(x,s) – поперечное продольное перемещение вдоль оси s, через W(x,s) – нормальное перемещение вдоль оси n.

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

 

 

U (x,s) ;

 

 

 

1

x

 

 

 

 

 

 

 

 

2

V (x,s)

W (x,s)

;

 

 

 

s

 

R

 

 

 

V (x,s)

 

U (x,s)

;

(2.6)

 

 

x

 

s

 

 

æ1 2W (x,s) ;

x2

æ2 2W (x,s) ;

s

2W (x,s) .x s

21

Компоненты вектора перемещения согласно вариационному методу В.З. Власова представим в виде разложений

 

n

 

 

U (x,s) 1 Ui (s) i (x) ;

 

 

i 1

 

 

n

 

 

V (x,s) 2 Vj (s) j (x) ;

(2.7)

 

j 1

 

 

n

 

 

W (x,s) 3 Wd (s) fd (x) .

 

 

d 1

 

Число п1, п2,

п3, обобщенных Ui (s), Vj (s), Wd (s)

и единичных

i (x), j (x), fd (x)

функций перемещений произвольно.

 

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

i (x), j (x), fd (x) (i 1, ..., n1; j 1, ..., n2; d 1, ..., n3 )

Умножая первое из уравнений (2.2) на k (x), второе на g (x), а третье на fd (x) и интегрируя их по всей длине продольной полоски

выделенного участка оболочки l, получим интегральные условия равновесия элементарной полоски ds (см. рис.2.3)

 

l

N

l S

l

 

 

 

x1 k (x)dx s

k (x)dx q1 k (x)dx 0

 

 

0

 

 

0

0

 

l

S

l

N

l

(2.8)

x

g (x)dx

s2 g (x)dx q2 g (x)dx 0

0

 

 

0

 

0

 

l

 

2

M2 1

 

 

l

 

2

 

l

 

2

2

 

 

 

fl (x)dx 2

 

 

H fl (x)dx

 

 

M2 fl (x)dx

0

x

 

0

x s

0

 

s

 

 

 

 

l

N2

 

 

 

 

l

 

 

 

 

 

 

 

 

fl (x)dx qn fl (x)dx 0 .

 

 

 

 

R

 

 

 

 

0

 

 

 

 

 

0

 

 

 

22

Интегрируя по частям один раз члены уравнений (2.8), содержащие первую производную по координате x, и два раза члены, содержащие вторую производную по x, получим

 

 

 

l

S

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

k (x)dx

N1 k (x)dx

 

N1 k

 

 

 

 

 

s

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1 k (x)dx 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l N2

g (x)dx

l

 

 

 

 

 

 

 

 

 

 

 

l

 

(2.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

S g (x)dx

 

S g

 

0

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q2 g (x)dx 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

l

H

 

 

 

l

2M2

fl

(x)dx

 

M1 fl (x)dx 2

s

fl (x)dx

 

 

s

2

 

0

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

l

N2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

fl (x)dx

 

M1 fl

 

l0

 

Q1 fl

 

l0 qn fl (x)dx 0 .

 

 

 

 

 

 

 

 

R

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Выражения в прямых скобках означают работу нормальных усилий

N1 сдвигающих

S,

 

обобщенных

поперечных

 

 

 

когда

 

усилий Q1,

Q1 M1 2 H , и изгибающих моментов, М1, действующих по концам

x s

выделенной продольной полоски оболочки, на возможных перемещениях. Эти силовые факторы определяются граничными условиями на поперечных краях оболочки при x=0 и x=l и могут быть отнесены к свободным членам от нагрузки.

Отметим, что этими же выражениями в прямых скобках учитывается податливость ребер на сжатие и изгиб, при наличии последних. В этом случае, подобно составным стержням из работы [64], эти члены входят в выражения для подсчета коэффициентов приведенной ниже системы уравнений (2.10). Подробнее учет податливости разберем в последующих главах, в главе 3 – податливость на сжатие, в главе 4 – на изгиб.

23

Подставляя выражения для N1, N2 , S, M1, M 2 и Н из уравнений

(2.3) в (2.9) и учитывая при этом выражения для деформаций и перемещений (2.6-2.7), получим следующую систему разрешающих дифференциальных уравнений равновесия элементарной продольной полоски ds для каждого участка, на котором разбита пологая оболочка с отверстием:

n

 

 

n

 

 

 

 

 

 

n

 

 

 

1 akiUi (s)

1 bkiUi (s) 2 (ckj ckj )Vj (s)

 

i 1

 

 

i 1

 

 

 

 

 

 

j 1

 

 

 

n

 

 

 

 

n

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

(dkd d

kd )Wd (s) (ekd ekd )Wd (s) q1k 0 ;

 

d 1

 

 

 

 

d 1

 

 

 

 

 

 

n

 

 

 

 

n

 

 

 

 

n

 

1 (cgi cgi )Ui

(s) 2 fgjVj (s) 2 Vj (s)

 

i 1

 

 

 

j 1

 

 

 

 

j 1

 

n3

 

 

n3

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

(kgd

kgd

 

 

 

(s) q2g 0 ;

(2.10)

hgdWd (s)

kgd )Wd

d 1

 

 

d 1

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

1 (dli dli )Ui (s) 1 (eli eli )Ui (s)

 

 

i 1

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

n2

 

n2

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

hljVj (s) (klj klj klj )Vj (s)

 

 

j 1

 

 

j 1

 

 

 

 

 

 

n3

 

 

n3

 

 

 

 

 

 

 

 

 

 

 

lldWd (s)

(mld

m

ld

 

mˆld mld )Wd (s)

 

d 1

 

 

d 1

 

 

 

 

 

 

 

 

 

 

 

 

n3 (nld nld nˆld nld )Wd (s) q3l 0. d 1

(k 1, ..., n1; g 1, ..., n2; l 1, ..., n3)

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

24

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

a 1 l B (x) (x)dx , a a

 

ki

 

 

2

 

 

 

i

 

k

 

 

ki ik

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

bki B i (x) k (x)dx , bki bik

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

1 l

B (x) (x)dx ;

 

 

kj

 

 

 

 

 

2

 

 

j

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

ckj Cfd (x) k (x)dx ;

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

dkd Cfd (x) k (x)dx ;

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

kd (1 ) Cfd (x) k (x)dx ;

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

ekd

Cfd (x) k (x)dx ;

(2.11)

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

ekd

1

 

l

Bfd (x) k (x)dx ;

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

cgi B i (x) g (x)dx ;

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

c

 

 

 

1 l

B (x) (x)dx ;

 

 

 

 

gi

 

 

 

 

2

 

 

i

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fgj (B EFj ) j (x) g (x)dx ,

fgj f jg

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 l

 

 

 

 

 

 

 

 

 

 

ggj

2

 

 

 

B j (x) g (x)dx ,

ggj g jg

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

hgd (C ES j ) fd (x) g (x)dx ;

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 l

 

 

 

 

 

 

 

 

 

kgd

 

 

 

(B EFj ) fd (x) g (x)dx ;

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

25

l

kgd (1 ) Cfd (x) g (x)dx ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kgd Cfd

(x) g (x)dx ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dli (1 ) C i (x) fl (x)dx ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dli C i (x) fl (x)dx ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eli

 

 

1

l B i (x) fl (x)dx ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eli C i (x) fl (x)dx ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hlj (C ES j ) j (x) fl (x)dx ;

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

(B

EF

 

 

 

(x) f

(x)dx ;

 

(2.11)

 

 

 

lj

 

 

R

)

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

klj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 ) C j

(x) fl (x)dx ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

l

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

klj C j (x) fl

(x)dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lld (D EJ j ) fd (x) fl (x)dx ,

lld ldl ;

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mld 2(1 ) Dfd (x) fl (x)dx , mld mdl ;

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

 

l

(C ES ) f (x) f (x)dx ,

 

 

 

;

 

m

 

 

m

m

 

R

 

 

ld

 

 

 

 

 

 

 

j

 

d

 

l

 

 

 

 

 

 

ld

dl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mˆld Dfd (x) fl (x)dx ,

mˆld mˆ dl ;

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

26

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dfd

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mld

(x) fl

(x)dx , mld mdl ;

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

1

l

(B EF

) f

d

(x) f

(x)dx ,

 

n

n

;

 

 

ld

 

R2

0

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

ld

dl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ndl ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nld Dfd (x) fl (x)dx , nld

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nˆ

1

 

 

l Cf (x) f (x)dx , nˆ

 

 

 

n ;

 

 

 

 

 

 

 

 

ld

 

 

R

 

d

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

ld

 

dl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

1

 

 

l

Cf (x) f (x)dx , n

 

nˆ ;

 

 

 

 

R

 

 

 

 

 

ld

 

 

 

 

d

 

 

 

 

l

 

 

 

 

 

 

 

 

 

ld

 

dl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сkj cgi ,

 

 

ckj cgi

 

 

 

 

 

 

 

 

 

 

 

 

dkd dli ,

 

 

 

 

kd

 

 

 

 

 

 

 

k = i

 

 

 

 

 

 

 

 

d

dli

 

 

 

 

 

 

 

 

 

 

 

ekd eli ,

 

 

ekd eli

 

 

 

 

 

d = l

 

 

 

 

 

 

h

h ,

 

 

k

k

 

 

 

 

 

g = j

 

 

 

 

 

 

 

 

 

gd

 

 

 

lj

 

 

 

 

gd

 

 

 

 

 

lj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kgd klj ,

 

 

 

klj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kgd

 

 

 

 

 

 

 

Коэффициенты свободных членов от внешней

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l0 ;

 

 

 

 

 

 

q1k q1 k (x)dx

 

N1 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l0 ;

 

 

 

 

 

 

q2g q2 g (x)dx

S g

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

q3l qn fl (x)dx Q1 fl l0 M1 fl l0 .

0

(2.11)

нагрузки

(2.12)

Понимая квадратуры (2.11-2.12) в смысле интегралов Стильтьесса, легко учитывается дискретное расположение ребер. Для этого к интегралам при непрерывном распределении по всей длине выделенной

27

продольной полоски ds дифференциалов Bdx, Cdx, Ddx, qdx следует добавить величины, представляющие собой суммы произведений из сосредоточенных в определенных точках xj конечных факторов

( EFj , ES j , EJ j , сосредоточенных сил) и значений соответствующих подынтегральных функций в этих точках xj .

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

2.3.Решение разрешающих уравнений

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

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

 

 

(s) 0 .

(2.13)

Y (s) A Y

Согласно работе [68] будем искать решение системы в виде:

Y e s , Y

 

e s , ...,

Y

 

N

e s

(2.14)

1 1

2

2

 

N

 

 

 

Подставляя в систему (2.13) и сокращая на множитель e s , будем иметь для определения постоянных 1, 2 , ..., N систему уравнений

a11 1 a12 2 ... a1N N 0 ;

a2N 1

a22 2 ... a2N N 0 ;

(2.15)

…………………………………….

 

aN 1

aN

2 ... aNN N 0.

 

1

2

 

 

28

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

a11

а12 ...

а1N

 

 

 

 

а21

a22

а2N

=0.

(2.16)

......…………………….

 

 

аN1

аN2 aNN

 

 

 

 

 

 

 

Уравнение (2.16) есть алгебраическое уравнение степени N со старшим членом ( )N . Если это уравнение имеет N различных корней

1, 2 , ..., N ,

то, подставляя каждый из них j в коэффициенты системы (2.15),

будем иметь N однородных уравнений для соответствующих неизвестных 1, 2 , ..., N с определителем равным нулю, и сможем

получить для этих неизвестных решение, отличное от нулевого. Таким образом, по формулам (2.14) получим N линейно-независимых решений системы (2.13), и их линейная комбинация даст общий интеграл. Если же уравнение (2.16) имеет корень кратности k, где k 1, то этому корню должны соответствовать k линейно-независимых решений. Если удается получить линейно-независимых векторов 1, 2 , ..., N , где k, то все

решений соответствующих корню кратности k будут иметь вид (2.14), а оставшиеся k решений должны содержать еще множителем

полином от S.

Разберем последний случай более подробно. Пусть корню кратности k соответствует только одно решение вида (2.14). Тогда подобно работе [68], найдем матрицу (A k ) в степени k

 

 

 

a11

 

a12

a1N

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(A

)k =

 

a

21

a

a

.

 

 

(2.17)

k

 

 

 

 

22

2 N

 

 

 

 

 

 

 

……………………………………..

 

 

 

 

 

 

aN1

 

aN 2

aNN

 

 

 

29

Из матрицы (2.17), имеющей ранг (N–k), найдем k линейно-

независимых векторов. Таким образом, получим полную матрицу N

состоящую из N линейно-независимых

1

собственных векторов

1, 2 , ..., N , расположенных по столбцам, с определителем отличным

от нуля. Затем преобразованием подобия N 1

A N приведем матрицу А

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

имеющей

характеристическое число

k кратности

k. Составляем

матрицу

B (D k ) , имеющую

характеристическое

число =0,

кратности k. Определив наименьшее число l, такое что Bl 0 , найдем

матрицу

C B l 1 . Определяя ранг матрицы С, возьмем те векторы,

которые не удовлетворяют системе уравнений С 0 , и, совершая над

ними последовательное (l–1) раз преобразование В, построим серии новых ортов, состоящие из векторов .Если после этого в

подпространстве, определяемом написанной системой, останутся еще векторы, т.е., если k>l, то применяя к ним последовательное преобразование B, мы получим новые серии ортов и т.д. Построив,

таким образом, матрицу N2 k-го порядка, состоящую из k линейно-

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

1

 

 

J

( ),

J

 

(

 

),

..., J

 

(

 

)

 

,

(2.18)

N

A N

 

2

2

p

p

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

где матрицы J ( ) имеют вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

...0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J ( )

 

 

1

 

0

 

...0

 

 

 

.

 

 

 

 

 

 

(2.19)

 

 

 

 

 

 

 

0

1

 

 

...0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

0

 

...

 

 

 

 

 

 

 

 

 

 

 

Нижний значок указывает порядок матрицы, а аргумент дает число, стоящее на главной диагонали, т.е. для простого корня =1, для кратного >1. Для получения такой матрицы N , надо подматрицу, состоящую из k столбцов матрицы N1 , соответствующих корню

30

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