![](/user_photo/_userpic.png)
книги из ГПНТБ / Амиро И.Я. Ребристые цилиндрические оболочки
.pdfВ соответствии с принятыми выше предположениями вместо уравнения (11.87) будем рассматривать систему уравнений:
1 + |
7Сѵ2 - |
2Ч Х 2 + «а2 + Чс) С1 + Ус) - |
б2! %4 + |
|
||
+ |
2 а д 4 [%2ЬПз (х) + (1 + а Ѵ ) Fn, (х) — 2П%Пз(х)1 + |
|
||||
+ 2'/2 (1 — V2 + |
а У ) {УсЬПз(х) — 26сх$Пя(х) + |
\ Х ^ П] (х) + |
||||
|
+ 2асХ4 [Ln, (X) Fn, (X) - £;,3 (X)]} = |
0, |
(II. 89) |
|||
Здесь |
Dik — 0 |
(/ = л3 + 1, rt3 + |
2, .. .). |
|
(11.90) |
|
|
|
|
|
|
|
|
^л,(х) = У] bikLik, |
Soj (х) = 5] ЬЬІѵк, |
F n, (х) —Л bikFik- |
||||
|
/=I |
|
z=i |
|
*=i |
|
Численное решение системы (11.89), (11.90) в принципе может быть найдено всегда. Необходимо лишь выбрать п3так, чтобы корни системы уравнений (11.89), (11.90) были близки к корням исход ного уравнения (11.87).
Предположим, что жесткость ребер на растяжение (ус и бс) не оказывает существенного влияния на искомую величину п3. Тогда для выбора п3 можно использовать следующую систему неравенств:
|
2т1г (1 - v |
2 + a X K |
S |
b]kFlk « 1 , |
(11.91) |
|
|
|
l=n2-\r\ |
|
|
|
1 — V2 + а2Хр( |
|
|
||
|
1 — V2 + (а2 + щ ) %4 Хрі |
|
|||
|
1 - |
V2 + a Y |
■рі |
< 1. |
(И-92) |
|
|
||||
|
і+ 2 а дpl 1 — V 2 |
bpi |
|
^ |
|
|
+ ( а 2 + Т|( |
і,+і |
|
|
|
|
|
|
|
|
|
Здесь |
xt — /,-й корень |
уравнения |
(11.89), имеющего всего 8п3 + 4 |
||
корня; |
% — р-й корень |
/-го уравнения |
системы (11.90), |
каждое |
уравнение которой имеет восемь корней.
Условия (11.91) и (11.92) получены в результате сравнения урав нения (11.87) с уравнениями (11.89) и (11.90). При этом уравнения (П.87) и (11.89) сравнивались непосредственно, а для сравнения с уравнениями системы (II. 90) уравнение (11.87) было представлено в виде Ош+Тш = 0, где Тік— некоторая известная функция %. Параметр п3 предлагается определить так, чтобы условия (11.91), (11.92) удовлетворялись с заданной точностью.
Рассмотрим, каким образом изменится общее решение одно родной системы уравнений равновесия (11.86) и граничные усло вия (11.88) после замены уравнения (11.87) системой уравнений
60
(11.89), (11.90). Нетрудно заметить, |
что при Xl0s-^ X P1 выражения |
Lik, i/h, Kik, Mik, Fik стремятся к |
бесконечности. Считая, что и, |
V, w должны оставаться ограниченными при такой замене, пола-
гаем |
IW1 |
я. |
/~>S |
= |
С0 |
|
(%0s^ |
|
Т. |
/ТТ ОС\ |
результате |
||||||
при X0s — XPlС0 |
---- 57------- • |
Іогда из (11.86) в |
|||||||||||||||
предельного |
|
|
|
|
|
ZOlk |
получаем следующие формулы для |
||||||||||
перехода xJs^X |
|
||||||||||||||||
определения |
перемещений: |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
П |
а 8 П з + 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
“ |
= |
S |
S |
CSMjé.t.e*»*6 COS/Â0 + |
|
|
|
||||||
|
|
со |
|
|
/=0 і1=1 |
|
|
|
|
|
|
|
|
|
|||
|
+ |
/ 8л3-1-4 |
|
|
|
|
|
8 |
|
\ |
|
|
|
||||
|
£ |
|
2 С о ' А ° і Ы і е х ^ |
+ £ C lp iA°lPle%p n c o s |
Ik Q , |
|
|||||||||||
|
|
l= n ,+ l |
\ ( ,= I |
|
|
|
|
|
Pf= I |
J |
|
|
|
||||
|
|
|
|
|
|
|
n , 8П |
3+4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
v = |
Y, |
t |
c o'Bik.t,eXifi sin IkQ + |
|
|
|
|||||||
|
|
|
|
|
|
l—0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I8rts-b4 |
|
|
|
|
|
8 |
|
|
\ |
|
|
||
|
+ |
Ѵ |
( |
у |
С о'В°ІЫіех^1 |
+ |
£ |
C^B °lpie%Pil |
sin IkQ, |
(11.93) |
|||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l=n>+\ \ t, = 1 |
|
|
|
|
|
Т |
1 |
> |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
л, 8n3-f4 |
Co1C°ik,tie%i£ cos IkQ + |
|
|
|
||||||
|
|
|
|
w |
= |
Y i |
S |
|
|
|
|||||||
|
|
|
|
|
|
/=0 fI=l |
|
|
|
|
|
|
|
|
|
||
|
|
°® /8Л3+4 |
|
|
|
|
|
|
8 |
|
\ |
|
|
||||
|
+ |
£ |
( |
£ |
|
|
|
|
|
+ |
|
Z C tfiC l& M co slkQ . |
|
||||
Здесь |
/=л3+1\ |
(,=1 |
|
|
|
|
|
pl=1 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A°iP[ = 23 [x2, + |
a % k [x% - |
T4 |
v |
/2^2) j + |
44Xp< (Z*ft* + vX2;), |
||||||||||||
B%t = ~ d z {l + |
|
|
a*Ärt) /ÄXP/ + d4 [(V - |
(2 + |
V) X“ j lk, |
||||||||||||
|
|
|
|
C°iPl = |
d, [/V |
+ |
vXp^] xP; + |
d3A/b |
|
|
|
||||||
|
A» = |
( f i ' - |
/V )2, |
d. = |
- |
|
d5da (1 - |
V2 + |
cftc*), |
|
|||||||
|
|
dt = |
d8d1 (1 — V2 + a2xj ), |
dt |
— vd2, |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
* |
|
|
|
|
|
4 |
= VcX2 |
X2Pl + |
°&ik (x2P{ ~ |
j é |
v |
№ )] - |
6cX4P/ (VX2Z+ m |
, |
|||||||||
|
|
|
|
d2 = |
YcXo. ( ^ 2 + |
n |
2 ) — ÖCX2 А |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
СЛР, Ik' |
|
|
|
61
Переход от уравнения (11.87) к системе (11.89), (11.90) позво ляет упростить граничные условия задачи (11.88). Осуществив
предельный переход |
от |
|
к %ti и |
, |
из |
(11.88) |
получим |
такую |
||||||||||
систему |
граничных |
условий: |
|
|
|
|
|
|
|
|
|
|||||||
Eh |
- |
оо |
8лз+4 |
|
|
|
|
|
|
8л,+ 4 |
|
|
|
|
||||
£ |
2 |
Co'TUe^cos т |
+ ¥ с (Ѳ ) £ |
С ' о |
' Т ^ |
1 |
+ |
|||||||||||
(1 - |
Ѵ2)г |
|||||||||||||||||
/=о |
і,=1 |
|
|
|
|
|
|
t.=i |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
+ |
2 |
|
2 |
|
|
cos im |
= |
T ш, |
|
|
|
||||
|
|
|
|
|
|
P^=I |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
8л ,+ 4 |
|
|
|
|
|
|
|
|
|
|
|
|
2(1 |
|
|
( 2 |
2 |
|
|
sin /fee + |
|
|
|
||||||
|
|
|
|
|
|
|
і=і (,=і |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
оо |
8 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
+ |
2 |
|
2 |
|
e%p?c * lS U |
sin |
|
= |
s 10, |
|
|
|
|||
|
|
|
г=л,+і p;=i |
|
|
|
|
|
|
|
|
|
||||||
|
Eh |
■ со |
8яя+ 4 |
|
|
|
|
|
|
8na+ 4 |
|
|
|
|||||
|
|
2 ] |
2 |
C'0'Q\\lkex‘fi cos /fee - |
¥ c(Ѳ) £ |
CfrQi'e*'«6 + |
||||||||||||
( l - v > |
|
|||||||||||||||||
|
<=# (,=1 |
|
|
|
|
|
|
|
i,=l |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
- |
|
8 |
|
|
|
|
|
|
_ |
|
|
|
(H-94) |
|
|
|
+ |
2 |
|
2 ^ W X ^ c o s / f e e |
|
= |
Q 10’ |
|
|
|
|||||||
|
|
|
/ = л ,+ І p/ = l |
|
|
|
|
|
|
|
|
|
|
|||||
|
Eh |
' |
со |
8л,“j—1 |
|
|
|
|
|
|
8 n ,+ 4 |
|
|
|
||||
|
|
2; |
2 |
|
|
|
cosim- ¥ c (Ѳ) 2 |
G° G^ |
‘5 + |
|||||||||
1 — V2 |
|
|
|
|
||||||||||||||
|
|
i=o |
(,=1 |
8 |
|
|
|
|
|
|
<i=i |
|
|
|
||||
|
|
+ |
0 0 |
|
e V C ^ C ^ c o s /fe e |
|
|
|
|
|
|
|||||||
|
|
2 |
|
2 |
|
G 10' |
|
|
|
|||||||||
Значения 7 ^ ft, |
S ^ ft, |
Q^ifc, |
G\llk можно |
получить |
соответственно |
|||||||||||||
из T\'ik, |
S\jk, Qlли, G\,ik, |
заменив А%.3, |
B°ik,s, |
C°k.s |
на |
А?а,Р;, |
|
|||||||||||
с?,.рг. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Анализируя |
формулы |
(11.89), (11.90), (11.93) и (11.94), нетрудно |
убедиться, что решение сформулированной в начале параграфа краевой задачи сводится к определению корней полиномов 8 п3 + 4 и 8-й степени и к решению систем линейных алгебраических урав нений 8п3 •+■ 4 или 8-го порядка.
Следует указать, что изложенный общий метод определения на пряженно-деформированного состояния ребристых оболочек в ка ждом конкретном случае может быть видоизменен. Так, например,
62
всегда можно найти такое значение /2, при котором
а Ѵ 4 » 1. |
(11.95) |
Если /2 > п3, то уравнения (11.90) можно записать в таком виде:
( Â - f k 2)4 = 0. |
(ІІ.90а) |
Все корни этого уравнения определяются элементарно (%Р[= ± Щ.
Соответствующие им состояния оболочки можно трактовать как сумму состояний бесконечно длинных тонких полос \ нагруженных по краям периодическими нагрузками. Подчеркнутые слагаемые в соотношениях (11.93) следует в этом случае представить в таком виде:
в формуле для и
- |
\е~ Щ (Со2 - С1о' - |
IktCo) |
+ |
elokHC14 + Со3 + |
Щ С І 4)}, |
в формуле для V |
|
|
|
|
|
|
е -“* (2Cf + i t * |
Ci' - |
|
tt|C ? ) - |
|
|
- e " s (2CS* + i i L C ? |
+ |
i i L t t l C i ‘), |
(11.96) |
|
в формуле для w |
|
|
|
|
|
|
- l k \ ( C I 5 + l k ^ c l 6 ) + g U t { C 1 7 + |
|
|||
Вернемся теперь к уравнениям (11.89), (11.90) и неравенствам |
|||||
(11.91) , |
(11.92) Нетрудно заметить, |
что при /2< л 3 можно су |
щественно упростить уравнение (11.89). Действительно, в этом
случае, |
по крайней мере, |
при ус = |
бс = |
0 оно разлагается на мно- |
|||||||
жители. Выделив |
из |
него |
множитель |
]~[ (X2t — f k 2)2, можно |
тем |
||||||
самым |
понизить степень уравнения на 4 (л3— /2). Полученные та- |
||||||||||
|
|
|
|
|
|
|
|
|
Лз |
|
|
ким образом 4(л3 — /2) корня |
уравнения |
П |
— l2k2)2 = 0 дают |
||||||||
4(п3 — /2) частных |
решения вида |
(11.96) |
и |
позволяют соответст |
|||||||
вующим образом упростить формулы для вычисления и и ѵ. |
|
||||||||||
При |
a2k4(n3+ |
I)4 |
1 можно |
упростить |
неравенства (11.91) и |
||||||
(11.92) |
. Нетрудно заметить, |
что |
при |
аг < |
неравенства |
(11.92) |
|||||
1 Поскольку при |
выполнении условия |
(11.95) |
Ik |
всегда велико, то взаим |
ное влияние краев полос будет несущественным. В этом случае рассматрива емые состояния ребристой оболочки следует трактовать как сумму состояний бесконечно длинных тонких полуполос, нагруженных по краям периодически ми нагрузками.
63
удовлетворяются, если
(11.97)
Величина /( быстро убывает с возрастанием I и даже для оболочек, усиленных очень тонкими ребрами (b2[ k 1), становится пренебре
жимо малой по сравнению с единицей при / = |
> |
5. При I < |
/3 |
значение U может оказаться большим (Іх = |
1,577 |
при bik= |
1). |
Очевидно, что в этом случае неравенства (11.92) и (11.97) не выпол няются и формулы (11.93) и (11.96) не пригодны.
Если п3 < /3, а а2/г4п з> 1, то для вычисления больших корней уравнения (11.87) вместо (11.90) можно предложить такую систему уравнений Е
(р; = 5, 6, 7, 8), (/ = Яд -}- 1, П3-ф 2, . . .).
Уравнения (11.98) описывают растяжение и сдвиг неподкрепленной бесконечной полосы периодической нагрузкой, уравнения (11.99) — циклически симметричный изгиб бесконечной полосы, усиленной абсолютно жесткими ребрами. Для решения уравнения (11.99) предлагается использовать метод последовательных прибли
жений. Поскольку при малых I для реальных оболочек ЬЧк~ 1, то аР[ не будут зависеть от характеристик обшивки и ребер. При
нимается, что в нулевом приближении а ^ = 1. Последующие при ближения вычисляются по формуле
± (сг ± idi),
(11.100)
где т — номер приближения; і = У ■— 1 • Вычисления по формуле (11.100) выполнены для / = 2, 3, 4, 5.
Для этих значений I предложенный процесс приближений сходится достаточно быстро и при вычислении С[ и dt можно ограничиться первым приближением.
Результаты вычисления ct и Д (/ > 2) приведены в табл. 1. Сопоставляя значения сг и d2, найденные во втором приближении,1
1 Уравнения системы (11.99) выведены в предположении, что rjo ^ a 2.
64
с их значениями, полученными в первом приближении, нетрудно убедиться, что они отличаются несущественно. Очевидно, что при I > 2 первое и второе приближения будут отличаться еще меньше.
Необходимо отметить, что предложенный процесс последователь ных приближений приводит во втором приближении к противоре чию (корни уравнения с вещественными коэффициентами не явля ются комплексно сопряженными); однако поскольку первое и вто рое приближения различаются несущественно, то в дальнейшем всегда можно ограничиться первым приближением.
|
|
|
|
|
|
|
Т а б л и ц а 1 |
|
|
|
т-=0 |
т-=1 |
|
|
п 1=2 |
|
|
1 |
ч |
Ч |
ч |
ч |
|
|
|
d"l |
|
Ч |
d'l |
c"l |
|||||
2 |
1 |
0 |
1,01 |
0,16 |
1,00 |
1,03 |
0,17 |
0,16 |
3 |
1 |
0 |
1,00 |
0,07 |
— |
— |
— |
— |
4 |
1 |
0 |
1,00 |
0,04 |
— |
— |
— |
— |
5 |
1 |
0 |
1,00 |
0,02 |
— |
— |
— |
|
При 1 = 1 (/, = 1,577) не удается получить величины с; и d; по формуле (11.100), использовав небольшое число приближений. В связи с этим для вычисления корней уравнений (11.99) предла гается другой прием: корни этого уравнения при I = 1 отожде ствляются с минимальными корнями полинома, полученного из (11.99) путем пренебрежения в суммах по Іг слагаемыми с Іг > I' + + 2. При этом в полученных полиномах после приведения подобных членов сохраняются лишь свободные члены и слагаемые, содержа щие аРі во второй и в четвертой степени.
Для оценки погрешности, вносимой в расчет первым допуще нием, выполнены вычисления, в результате которых получено: при
I' = 0 , 1, 2, 3; сг и |
соответственно равны 0,67; |
0,63; 0,61; |
0,60 |
и 0,35; 0,31; 0,29; 0,29. Нетрудно заметить, что уже |
при I' = |
1 по |
|
лучены значения корней, близкие к предельным. |
|
|
Оценка второго допущения проведена путем сопоставления кор ней биквадратного уравнения и минимальных корней уравнения восьмой степени, полученного для I' = I. Последние определялись через корни биквадратного уравнения последовательными прибли
жениями, |
где было получено |
с), = |
0,66, |
d' = 0,34; |
с"=0,63, d" = |
|
=0,35; |
в |
результате точного решения |
уравнения |
восьмого по |
||
рядка |
получено сг = 0,67, |
dx = |
0,36. |
Сравнивая |
эти данные с |
данными, приведенными выше, нетрудно убедиться, что корни би квадратного уравнения близки к минимальным корням уравнения восьмой степени.
•В рассмотренном случае (l3 > п3, a2k4n43 1) для вычисления и и о можно использовать соотношения (11.93), в которых подчер
5—39 |
65 |
кнутые |
слагаемые |
заменены |
формулами (11.96). При вычислении |
|
w в соотношениях (11.93) подчеркнутые слагаемые (при /8 > I > п3) |
||||
следует |
выбирать |
в таком виде: |
|
|
п ІЬ —{ci+td^Lkl |
IQ —Uj— |
l7 (Cj+JdjjMfc |
lB (ci—idi)lkl |
|
C0 e |
-t C0 e |
+ CQe |
+ C0 e |
|
|
|
|
|
(II. 101) |
В том случае, если используются формулы (11.96), соответствующие изменения вносятся и в граничные условия, причем условия (11.94) по-прежнему можно трактовать как точные. При использовании формул (II.101) условия (II.94), полученные с помощью предель ного перехода (при / < /3), уже нельзя считать точными.
Следует отметить, что неравенство (II.95) для широкого класса оболочек, как правило, выполняется, начиная с / > 2. Поэтому, чтобы использовать полученные результаты при расчете этих обо
лочек, необходимо уточнить систему уравнений |
(11.98), (11.99) |
|
для 1 = 1 . Для рассматриваемых |
оболочек а2&4« |
1, в этом случае |
уравнение (II.90) не распадается |
на два уравнения четвертой сте |
пени. Все восемь корней этого уравнения предлагается вычислить из (II.99) с помощью способа, использованного выше при отыска нии корней уравнения (II.99) для 1 = 1 . Будем предполагать, что ус и 6Сне оказывают существенного влияния на величину больших
корней |
характеристического |
уравнения и |
положим |
ус = |
6С= 0. |
|||
Тогда из (П.89) |
получаем следующую систему уравнений *: |
|||||||
(1 - |
V2) ( 4 ) 4 + |
а2 [(ХІ5)2 - |
( I k ff {[(Xjs)2 - |
( I k ff |
+ X x(Xi,)} = 0, |
|||
где |
|
|
|
|
|
|
|
(11.102) |
|
|
|
|
|
|
|
|
|
X'(4s) |
______________ b% |
|
|
|
|
|||
1— v2 -f- (g2 + T|e ) (xjs)4 |
|
|
|
|
||||
|
|
2t1e (Xos)4 I1— v2 + |
a2 (Xo,)4] |
|
|
|
|
|
Приняв |
а2<C 11c f\c (X o s)4 > 1 , 1 |
6a2k4> |
1, b2k » 1 |
при |
L= 1, |
пред- |
||
ставим (11.102) в таком виде: |
|
|
|
|
|
|
||
(1 - |
V2) (Xis)4 + а2[(Xis)2 - |
Ä2]2 {[(Xis)2 - |
k2]2+хг(Xis)} = |
0. |
||||
В этом случае |
|
|
|
|
|
|
(ІІ.ЮЗ) |
|
|
|
|
|
|
|
|
||
|
*>(ХІ51 |
|
|
|
|
|
|
|
|
|
2 [ l - v 2 + |
I ч4. |
+ У |
|
|
|
|
|
|
g2 (Xi,)' |
£ 2[(xis)2- W |
|
|
1 Нетрудно заметить, что при a2Ä4/4^» 1 (11.102) можно заменить системой (11.98), (11.99).
66
Положив Хо5 = ^а Р,> переписываем (11.103) в виде
(S і)4с4 . + (а р( — i f |
|
|
|
|
+ ■ |
|
|
= о. |
(11.104) |
К +(S,)41 |
+ 1 |
(< |
~ 1У |
|
|
l,=2 |
|
Сравнивая (11.104) и (11.99), нетрудно заметить, что корни урав нения (11.104), в отличие от корней уравнения (11.99), зависят от параметров оболочки. Решение уравнения (11.104) получено в пер вом приближении: предполагалось, что в сумме по можно ограни читься первым слагаемым. Полученное таким образом уравнение 12-й степени решалось на ЭЦВМ точно. Из найденных в результате вычислений трех существенно различных корней для определения напряженного состояния оболочки предлагается использовать два меньших по модулю корня, предельные значения которых (при
Sl ->• 0) совпадают с найденными выше корнями уравнений (11.98) и (11.99). Результаты вычисления указанных корней для различ
ных значений |
приведены в табл. 2. |
|
|
||
|
|
|
|
|
Т а б л и ц а 2 |
|
S4 |
с'і |
d'l |
с"\ |
h |
|
|
||||
|
0,910 |
0,91 |
0,30 |
0,88 |
0,51 |
|
0,819 |
0,92 |
0,29 |
0,86 |
0,50 |
|
0,728 |
0,93 |
0,28 |
0,83 |
0,48 |
|
0,637 |
0,95 |
0,27 |
0,81 |
0,47 |
|
0,546 |
0,96 |
Qi,21 |
0,79 |
0,45 |
|
0,455 |
0,97 |
0,24 |
0,77 |
0,45 |
|
0,364 |
0,98 |
0,21 |
0,75 |
0,43 |
|
0,273 |
0,99 |
0,19 |
0,73 |
0,42 |
|
0,182 |
0,99 |
0,16 |
0,71 |
0,40 |
|
0,091 |
1,00 |
0,11 |
0,69 |
0,38 |
Для обозначения действительных |
и мнимых |
частей корней, |
|
приведенных |
в табл. 2, приняты |
такие обозначения: с| = |
|
Re(ap>)1>2j3j4, |
Im (ир^^.з^і И Re (а р,)5,б,7,8' |
^ (а р,)5,б,7,8' |
Для определения малых корней уравнения (11.87) возможен и другой путь. Так, приняв п3 = 5, нетрудно заметить, что из урав нения (11.89) необходимо определить не более 44 меньших по модулю корней исходного характеристического уравнения (11.87) (11 суще ственно различных корней). Для реальных оболочек неравенство
5* 67
(11.95), как правило, выполняется при /2 = |
/4 < 5. Отсюда сле |
дует, что 2 (5 — /4) существенно различных |
корней можно опре |
делить из уравнений (11.98), (11.99). Таким образом, для реальных оболочек из уравнения (11.89) необходимо определить 11 — 2 (5—■
— /4) минимальных существенно различных корней уравнения (11.87). Поскольку разыскиваются малые корни, то можно отказа ться от удовлетворения условия (11.91), а для определения корней воспользоваться методом последовательных приближений, опре деляя в каждом приближении корни полинома 8/га + 4 степени (11.89) для заданного п3. Этот способ вычисления корней проверен
на |
примере |
двух |
оболочек: |
ус |
= 6с = |
0, с = 0,1258, |
rj* = 60, |
||
V = 0,3; |
ус = |
6с = 0, с — 1, т]* |
= |
30, ѵ = |
0,3, где с = cPk*, г\’с = |
||||
= |
ті ft4. |
|
|
|
|
|
|
|
(п3 = 3), |
36 |
Последовательно |
определялись |
корни |
полиномов |
28 |
||||
(п3 = |
4), 44 (пд = |
5) степени. Корни полиномов 1 |
14 и |
18-й сте |
пени определялись на ЭЦВМ непосредственно, минимальные корни
полинома 22-й степени — последовательными приближениями: |
ис |
||||
ходное |
уравнение |
заменялось последовательно уравнениями |
16,18 |
||
и 20-й |
степени |
(табл. 4, ./V — степень |
уравнения). |
Результаты |
|
вычисления пяти |
минимальных корней |
для первой |
из рассмат |
риваемых оболочек и семи минимальных корней для второй
оболочки приведены в табл. 3 и 5 (аи , |
bt~— соответственно дей |
||||
ствительная часть корня и коэффициент при |
мнимой части |
||||
корня). |
|
|
|
|
|
|
|
|
|
Т а б л и ц а 3 |
|
п,=3 |
лз—4 |
т,=5 |
|||
% |
|
% |
ьи |
% |
ьи |
0,12449 |
0,12117 |
0,12438 |
0,12108 |
0,12434 |
0,12104 |
0,75110 |
0,30620 |
0,74800 |
0,30705 |
0,74678 |
0,30741 |
1,27724 |
1,05907 |
1,27819 |
1,06158 |
1,27871 |
1,06261 |
1,57326 |
0,46204 |
1,56839 |
0,46020 |
1,56745 |
0,45977 |
2,16233 |
0,80567 |
2,02868 |
0,26361 |
2,15567 |
0,85041 |
|
|
2,10768 |
1,26276 |
|
|
Следует отметить, что предложенные способы вычисления кор ней характеристического уравнения^ 1.87) позволяют определить все корни характеристического уравнения с необходимой точностью. Постановка полученных корней в формулы (11.93) и граничные условия (11.94) позволяет в принципе получить решение любой за дачи о напряженно-деформированном состоянии ребристой обо лочки при современном уровне возможностей ЭЦВМ. Точность,
1 Использовалась подстановка У=х2.
68
|
|
|
|
Т а б л и ц а |
4 |
|
N = 3 2 |
N = 3 3 |
N = 4 0 |
|
|||
% |
b h |
% |
b t l |
а и |
b h |
|
0,12434 |
0,12104 |
0,12434 |
0,12104 |
0,12434 |
0,12104 |
|
0,74677 |
0.30741 |
0,74678 |
0,30741 |
0,74678 |
0,30741 |
|
1,27853 |
1,06089 |
1,27872 |
1,06262 |
1,27871 |
1,06261 |
|
1,54889 |
0,45476 |
1,56743 |
0,45959 |
1,56745 |
0,45977 |
|
1,67010 |
0,73462 |
1,97446 |
0,64019 |
2,15567 |
0,85041 |
|
|
|
|
|
Т а б л и ц а |
5 |
|
п,=3 |
пѵ=4 |
Пі=Ъ |
|
|||
|
ьи |
% |
ьи |
% |
К |
|
0,23177 |
0,21741 |
0,23161 |
0,21731 |
0,23155 |
0,21727 |
|
0,89689 |
0,30999 |
0,89444 |
0,31034 |
0,89347 |
0,31049 |
|
0,88397 |
0,51123 |
0,88600 |
0,512/6 |
0,88684 |
0,51336 |
|
1,67882 |
0,43169 |
1,68237 |
0,43122 |
1,68355 |
0,43116 |
|
2,64772 |
0,46034 |
2,66519 |
0,46710 |
2,70534 |
0,45313 |
|
2,03681 |
0,40256 |
2,03437 |
0,41637 |
2,03698 |
0,42186 |
|
3,07957 |
0,33814 |
3,04408 |
0,41570 |
2,95640 |
0,41255 |
|
|
|
|
|
Т а б л и ц а 6 |
||
|
* Ѳ = о |
|
|
k0= n |
|
|
У ; |
и* |
|
|
|
• |
|
|
|
u * |
0 7 * — 1 |
г , |
|
|
0,000 |
0 000 |
0,170 |
0,000 |
—1,665 |
0,248 |
|
0,533 |
0,106 |
0,169 |
0,155 |
—1,664 |
0,249 |
|
0,10710 |
0,210 |
0,165 |
0,310 |
—1,658 |
0,252 |
|
0,160-10 |
0,314 |
0,161 |
0,464 |
—1,640 |
0,253 |
|
0,213-10 |
0,419 |
0,159 |
0,615 |
—1,590 |
0,250 |
|
0,267-10 |
0,535 |
0,167 |
0,762 |
—1,471 |
0,235 |
|
0,320-10 |
0,675 |
0,192 |
0,907 |
—1,231 |
0,207 |
|
0,373-ТО |
0,861 |
0,239 |
0,106-10 |
—0,828 |
0,166 |
|
0,427-10 |
0,112-10 |
0,304 |
0,124-10 |
—0,314 |
0,120 |
|
0,480-10 |
0,144-10 |
0,364 |
0,144-10 |
0 |
0,814-10_1 |
достигнутая при реализации предложенного метода, определяется точностью нахождения корней уравнения (11.87) из приближенных уравнений (11.89), (11.90) или (11.89), (11.98), (11.99); ее всегда мож но' повысить, уточнив приближенные уравнения.
Отличный от изложенного метод определения напряженнодеформированного состояния ребристых оболочек предложен в [72].
69