24.07.2014 Views

Полная версия - Самарский государственный аэрокосмический ...

Полная версия - Самарский государственный аэрокосмический ...

Полная версия - Самарский государственный аэрокосмический ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Вестник Самарского государственного аэрокосмического университета, № 1, 2007<br />

Первые два частных решения<br />

( y , y , y , ) и ( , y , y , )<br />

11 21 31<br />

y41<br />

y найдем<br />

12 22 32<br />

y42<br />

путем интегрирования от начального значения<br />

ξ 0 с заданным шагом t при началь-<br />

0 =<br />

ных условиях по верхнему краю оболочки<br />

θ = α :<br />

y11 , 0<br />

= 1,<br />

y21,<br />

0<br />

= 0,<br />

y31,<br />

0<br />

= 0,<br />

y41,<br />

0<br />

= 1<br />

для первого решения и условиях<br />

y12 , 0<br />

= 0,<br />

y22,<br />

0<br />

= 1,<br />

y32,<br />

0<br />

= 1,<br />

y42,<br />

0<br />

= 0<br />

для второго решения.<br />

Полученные таким образом решения<br />

будут “возрастающими” по мере удаления от<br />

верхнего края оболочки вниз.<br />

Третье (<br />

13,<br />

y23,<br />

y33,<br />

y43)<br />

( , y , y , )<br />

14 24 34<br />

y44<br />

y и четвертое<br />

y решения целесообразно<br />

строить интегрированием от значения ξ = 0<br />

1<br />

с отрицательным шагом t ( t < 0 ) при начальных<br />

условиях<br />

y13 , 0<br />

= 1,<br />

y23,<br />

0<br />

= 0,<br />

y33,<br />

0<br />

= 0,<br />

y43,<br />

0<br />

= 1<br />

и, соответственно,<br />

y , y = 1,<br />

y = 1,<br />

0 .<br />

14 , 0<br />

= 0<br />

24,<br />

0 34,<br />

0<br />

y44,<br />

0<br />

=<br />

При этом получаются решения, “возрастающие”<br />

при движении от нижнего края<br />

оболочки вверх.<br />

2.3. Признаком линейной независимости<br />

решений ( ,...),...,( y ,...)<br />

y является неравенство<br />

нулю<br />

11 14<br />

определителя<br />

D( ξ)=<br />

y<br />

y<br />

y<br />

y<br />

11<br />

21<br />

31<br />

41<br />

y<br />

y<br />

y<br />

y<br />

12<br />

22<br />

32<br />

42<br />

y<br />

y<br />

y<br />

y<br />

13<br />

23<br />

33<br />

43<br />

y<br />

y<br />

y<br />

y<br />

14<br />

24<br />

34<br />

44<br />

В [3] показано, что<br />

ξ<br />

⎛<br />

⎞<br />

() ξ = D0 exp⎜∫( a11<br />

+ + a44<br />

) dξ⎟ ⎠<br />

. (22)<br />

D L ,<br />

⎝ 0<br />

где D<br />

0<br />

- значение определителя (22) при<br />

ξ = 0.<br />

В нашем случае согласно (3)<br />

a + a + a + a 0.<br />

11 22 33 44<br />

=<br />

Поэтому<br />

( ξ ) = D const<br />

D . (23)<br />

0 =<br />

Условие (23) можно использовать для<br />

контроля правильности вычисления значений<br />

функций y<br />

11,...,<br />

y44<br />

.<br />

2.4. Частное решение ( , y , y , y )<br />

y1q<br />

2q<br />

3q<br />

4q<br />

системы (1) , соответствующее заданной поверхностной<br />

нагрузке, целесообразно искать,<br />

применяя метод вариации постоянных [3]:<br />

y ∑ D y<br />

jq<br />

= 4 i=<br />

1<br />

i<br />

ji<br />

, j = 1 , 2,<br />

3,<br />

4 . (24)<br />

Здесь y<br />

ji - совокупность всех линейно<br />

независимых решений системы однородных<br />

уравнений (11), ( ξ)<br />

D - функции, определяемые<br />

из соотношений<br />

D<br />

i<br />

ξ<br />

∫<br />

1<br />

() ξ = D′<br />

()ξ ξ d<br />

для i= 1, 2 и<br />

D<br />

i<br />

ξ<br />

∫<br />

0<br />

() ξ = D′<br />

()ξ ξ d<br />

для i= 3, 4.<br />

i<br />

i<br />

i<br />

Сами производные ( ξ)<br />

′i<br />

(25)<br />

(26)<br />

D (i = 1,2,3,4)<br />

согласно методу вариации произвольных постоянных<br />

находятся из зависимостей<br />

4<br />

∑<br />

i=<br />

1<br />

или<br />

D ′ y =<br />

i<br />

ji<br />

f<br />

jq<br />

[ ][ Z] [ F]<br />

(27)<br />

Y = . (28)<br />

Здесь<br />

232

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!